1
$\begingroup$

After I have done with some Monte Carlo simulation of Gay-Berne model, I want to visualize how the N molecules look like in the domain. Here the N molecules all have the same ellipsoidal shape. The centre positions $\mathbf{r}_i$ and the orientations $\mathbf{u}_i$ are given. Is there any software that can plot the configuration of these N molecules based on the given information? Either 3D or 2D is fine.

Thanks!

Improve this question
$\endgroup$

1 Answer 1

Reset to default
0
$\begingroup$

Almost any mathematical package will suffice: Matlab, Mathematica, Sage, ... Here is an image made in Mathematica, with the ellipses drawn by the Circle[] function, and then oriented with Rotate[]:
          ellipses
One can also draw 3D ellipsoids in these packages. See the earlier MO question, "Maximual tetrahedra inscribed in an ellipsoid" for an example (using Rotate[Scale[Sphere[]]]).

Improve this answer
$\endgroup$
5
  • $\begingroup$ I found that pdeellip() in Matlab has similar functionality. Thanks a lot! $\endgroup$
    – Chong Luo
    Commented Oct 25, 2011 at 11:50
  • $\begingroup$ Yes. Here is the man page on that Matlab function, in case anyone wants to use it: mathworks.com/help/toolbox/pde/ug/pdeellip.html $\endgroup$
    – Joseph O'Rourke
    Commented Oct 25, 2011 at 12:11
  • $\begingroup$ For N=1089 particles, calling pdeellip() makes Matlab run out of memory. Strange. I'll try Mathematica and see if it's better. $\endgroup$
    – Chong Luo
    Commented Oct 26, 2011 at 17:08
  • $\begingroup$ $10^3$ particles should be no challenge; $10^6$ might be a challenge. I don't use Matlab as regularly as Mathematica, so I cannot help you. But you could post a question on a Matlab forum, and get quick expert advice. $\endgroup$
    – Joseph O'Rourke
    Commented Oct 26, 2011 at 18:17
  • $\begingroup$ I found Mathematica can handle $10^3$ particles quite easily, and produce nice graphs. I'll just use Mathematica. Thanks a lot! $\endgroup$
    – Chong Luo
    Commented Oct 29, 2011 at 10:33

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .