git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@10509 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp
2013-08-01 14:38:21 +00:00
parent bb72ec9d1f
commit 217a0a3b7f
2 changed files with 52 additions and 44 deletions

View File

@ -34,10 +34,12 @@ velocities, which may or may not be desirable. New atom IDs are
assigned to new atoms, as are molecule IDs. Bonds and other topology assigned to new atoms, as are molecule IDs. Bonds and other topology
interactions are created between pairs of new atoms as well as between interactions are created between pairs of new atoms as well as between
old and new atoms. This is done by using the image flag for each atom old and new atoms. This is done by using the image flag for each atom
to "unwrap" it out of the periodic box before replicating it. This to "unwrap" it out of the periodic box before replicating it.
means that molecular bonds you specify in the original data file that </P>
span the periodic box should be between two atoms with image flags <P>This means that any molecular bond you specify in the original data
that differ by 1. This will allow them to be unwrapped appropriately. file that crosses a periodic boundary should be between two atoms with
image flags that differ by 1. This will allow the bond to be
unwrapped appropriately.
</P> </P>
<P><B>Restrictions:</B> <P><B>Restrictions:</B>
</P> </P>
@ -47,25 +49,27 @@ that differ by 1. This will allow them to be unwrapped appropriately.
when replicating it in that dimension, as it may put atoms nearly on when replicating it in that dimension, as it may put atoms nearly on
top of each other. top of each other.
</P> </P>
<P>IMPORTANT NOTE: You cannot use the replicate command to replicate a <P>IMPORTANT NOTE: You cannot use the replicate command on a system which
molecule that is bonded to itself across a periodic boundary, so the has a molecule that spans the box and is bonded to itself across a
molecule is efffectively a loop. A simple example would be a linear periodic boundary, so that the molecule is efffectively a loop. A
polymer chain that spans the simulation box and bonds back to itself simple example would be a linear polymer chain that spans the
across the periodic boundary. More realistic examples would be a CNT simulation box and bonds back to itself across the periodic boundary.
(meant to be an infinitely long CNT) or graphene sheet. (Note that More realistic examples would be a CNT (meant to be an infinitely long
this only applies to such molecules that have permanent bonds as CNT) or a graphene sheet or a bulk periodic crystal where there are
specified in the data file. A CNT that is just atoms modeled with the explicit bonds specified between near neighbors. (Note that this only
<A HREF = "pair_airebo.html">AIREBO potential</A> has no such permanent bonds, so it applies to systems that have permanent bonds as specified in the data
is not a problem to replicate it.) The reason replicate does not work file. A CNT that is just atoms modeled with the <A HREF = "pair_airebo.html">AIREBO
with those systems is that the image flag settings described above potential</A> has no such permanent bonds, so it can be
cannot be made consistent. I.e. it is not possible to have every pair replicated.) The reason replication does not work with those systems
of bonded atoms have the same image flags, when they are unwrapped so is that the image flag settings described above cannot be made
as to be close to each other. The only way the replicate command consistent. I.e. it is not possible to define images flags so that
when every pair of bonded atoms is unwrapped (using the image flags),
they will be close to each other. The only way the replicate command
could work in this scenario is for it to break a bond, insert more could work in this scenario is for it to break a bond, insert more
atoms, and re-connect the loop for the larger simulation box. It is atoms, and re-connect the loop for the larger simulation box. But it
not clever enough to do this. So you will have to construct a larger is not clever enough to do this. So you will have to construct a
version of your molecule as a pre-processing step and input a new data larger version of your molecule as a pre-processing step and input a
file to LAMMPS. new data file to LAMMPS.
</P> </P>
<P>If the current simulation was read in from a restart file (before a <P>If the current simulation was read in from a restart file (before a
run is performed), there can have been no fix information stored in run is performed), there can have been no fix information stored in

View File

@ -31,10 +31,12 @@ velocities, which may or may not be desirable. New atom IDs are
assigned to new atoms, as are molecule IDs. Bonds and other topology assigned to new atoms, as are molecule IDs. Bonds and other topology
interactions are created between pairs of new atoms as well as between interactions are created between pairs of new atoms as well as between
old and new atoms. This is done by using the image flag for each atom old and new atoms. This is done by using the image flag for each atom
to "unwrap" it out of the periodic box before replicating it. This to "unwrap" it out of the periodic box before replicating it.
means that molecular bonds you specify in the original data file that
span the periodic box should be between two atoms with image flags This means that any molecular bond you specify in the original data
that differ by 1. This will allow them to be unwrapped appropriately. file that crosses a periodic boundary should be between two atoms with
image flags that differ by 1. This will allow the bond to be
unwrapped appropriately.
[Restrictions:] [Restrictions:]
@ -44,25 +46,27 @@ If a simulation is non-periodic in a dimension, care should be used
when replicating it in that dimension, as it may put atoms nearly on when replicating it in that dimension, as it may put atoms nearly on
top of each other. top of each other.
IMPORTANT NOTE: You cannot use the replicate command to replicate a IMPORTANT NOTE: You cannot use the replicate command on a system which
molecule that is bonded to itself across a periodic boundary, so the has a molecule that spans the box and is bonded to itself across a
molecule is efffectively a loop. A simple example would be a linear periodic boundary, so that the molecule is efffectively a loop. A
polymer chain that spans the simulation box and bonds back to itself simple example would be a linear polymer chain that spans the
across the periodic boundary. More realistic examples would be a CNT simulation box and bonds back to itself across the periodic boundary.
(meant to be an infinitely long CNT) or graphene sheet. (Note that More realistic examples would be a CNT (meant to be an infinitely long
this only applies to such molecules that have permanent bonds as CNT) or a graphene sheet or a bulk periodic crystal where there are
specified in the data file. A CNT that is just atoms modeled with the explicit bonds specified between near neighbors. (Note that this only
"AIREBO potential"_pair_airebo.html has no such permanent bonds, so it applies to systems that have permanent bonds as specified in the data
is not a problem to replicate it.) The reason replicate does not work file. A CNT that is just atoms modeled with the "AIREBO
with those systems is that the image flag settings described above potential"_pair_airebo.html has no such permanent bonds, so it can be
cannot be made consistent. I.e. it is not possible to have every pair replicated.) The reason replication does not work with those systems
of bonded atoms have the same image flags, when they are unwrapped so is that the image flag settings described above cannot be made
as to be close to each other. The only way the replicate command consistent. I.e. it is not possible to define images flags so that
when every pair of bonded atoms is unwrapped (using the image flags),
they will be close to each other. The only way the replicate command
could work in this scenario is for it to break a bond, insert more could work in this scenario is for it to break a bond, insert more
atoms, and re-connect the loop for the larger simulation box. It is atoms, and re-connect the loop for the larger simulation box. But it
not clever enough to do this. So you will have to construct a larger is not clever enough to do this. So you will have to construct a
version of your molecule as a pre-processing step and input a new data larger version of your molecule as a pre-processing step and input a
file to LAMMPS. new data file to LAMMPS.
If the current simulation was read in from a restart file (before a If the current simulation was read in from a restart file (before a
run is performed), there can have been no fix information stored in run is performed), there can have been no fix information stored in