Thank you for your reply. Without writing some code, I rather try another possible way you suggested, and I have some issues on format of read_data “body” section.
In the “in.file” contains the body style “body nparticle 2 2” and in the read_data, I write the “Bodies” like this:
I give “6 bodies” in the head lines and also for the “Bodies” section (24 lines total, 4 lines for each body), when the program read to the body section, it shows
“ERROR on proc 0: identifier in data file: 0 0 0 0 0 0 1.122640509 -0.779584938 -13.95639781 -0.155203758”
It can’t identify the 7th line. If I set “24 bodies”(but actually I have only 6 bodies) in the head lines, it would appear “Unexpected end of data file”.
When I run the example of body, it works. The body format of the data files are the same. So I can’t understand why it stopped.
(I deleted the space or Tab after the lines in the body section, but it doesn’t matter.) Maybe it just a simple problem, sorry about that.
I don’t understand why you are using atom_style body for ellipsoids.
You should just use atom_style ellipsoid. If you want to add
the extra particles at the ends of the ellipsoid, they can be
point particles. Atom_style ellipsoid allows that.
If I want to represent sphere1 as a point, it seems can’t use npt/asphere ensemble without setting shape to it.
Then I set small radius (shape) and small mass to sphere1. If I set R0=0 to the spring-like bond between the two sphere1(or the sphere1-ellipsoid bond), they would bounce fiercely, and the ellipsoids also rotate freely (either for the normal R0).
So I would like to use a 3-point segment as a rigid body which the central atom is assigned as a ellipsoid, then connect the bodies to form a chain. Perhaps the torque won’t take effect to the ellipsoid.
Now I use hybrid atom_style, is that possible to use ellipsoid & body style in the same time? (or there are something wrong when I use these styles, so the program can’t read)