       making a fast projection of a crytsallographic unit cell

• To: mathgroup at smc.vnet.net
• Subject: [mg80779] making a fast projection of a crytsallographic unit cell
• From: Roger Bagula <rlbagula at sbcglobal.net>
• Date: Sat, 1 Sep 2007 00:32:47 -0400 (EDT)

```The problem is that I want to plot representations of crystallographic
cells fast in Mathematica.

The face centerd cubic unit cell can be represented in 3d as a set of
points.
Hermann-Mauguin space groupo notation: Fm3m
Sodium Chloride and many metal crystals use this crystal cell.
(* points for cubic close pack/ face centered cubic*)
p = {{1, 1, 1}, {1, 0, 0}, {2, 1, 0}, {0, 0, 1}, {2, 2, 1}, {1, 2,
2}, {0, 1, 2}, {0, 1, 0}, {1, 2, 0}, {0, 2, 1}, {2, 0, 1}, {
2, 1, 2}, {1, 0, 2}};
pieces =
Complement[
Flatten[Table[{i, j, k}, {i, 0, 2}, {j, 0, 2}, {k, 0, 2}],
2], p]
(* showing as points*)
Show[Graphics3D[Table[Point[pieces[[n]]], {
n, 1, Length[pieces]}]], Boxed -> False]

(* creating the equivalent of a graphic primitive like a Cuboid[] for a
sphere*)

Clear[sphere]
<< RealTime3D`
(* RealTime3D allows you to drag the graphic around to see from
different view points*)
sphere[n_]:=ParametricPlot3D[{Sin[t]*Sin[p]/2-pieces[[n]][[
1]],Cos[t]*Sin[p]/2-pieces[[n]][],Cos[p]/2-pieces[[n]][]},{t,-
Pi,Pi},{p,0,Pi},Axes\[Rule]False,Boxed\[Rule]False]

(* plotting them at the points them at the points*)
gr=Table[sphere[n],{n,1,Length[pieces]}];

Show[gr]

The problem isn't doing it, but doing it faster.
There are 230 Federov space groups in 3d.
Roger Bagula

```

• Prev by Date: Re: Re: problem with parsing InputForm scientific notation (eg. 3.`*^-7)
• Next by Date: Re: [Mathematica 6] Bug in Labeled, Center align option
• Previous by thread: Re: problem with parsing InputForm scientific notation (eg. 3.`*^-7)
• Next by thread: Re: making a fast projection of a crytsallographic unit cell