diff --git a/doc/Manual.html b/doc/Manual.html index 4e9bbd96d1..ba5f912d7d 100644 --- a/doc/Manual.html +++ b/doc/Manual.html @@ -179,7 +179,7 @@ it gives quick access to documentation for all LAMMPS commands.
6.13 NEMD simulations
- 6.14 Extended spherical and aspherical particles + 6.14 Finite-size spherical and aspherical particles
6.15 Output from LAMMPS (thermo, dumps, computes, fixes, variables)
diff --git a/doc/Manual.txt b/doc/Manual.txt index d5a57cdca8..08e17dbb2d 100644 --- a/doc/Manual.txt +++ b/doc/Manual.txt @@ -126,7 +126,7 @@ it gives quick access to documentation for all LAMMPS commands. 6.11 "Visualizing LAMMPS snapshots"_howto_11 :b 6.12 "Triclinic (non-orthogonal) simulation boxes"_howto_12 :b 6.13 "NEMD simulations"_howto_13 :b - 6.14 "Extended spherical and aspherical particles"_howto_14 :b + 6.14 "Finite-size spherical and aspherical particles"_howto_14 :b 6.15 "Output from LAMMPS (thermo, dumps, computes, fixes, variables)"_howto_15 :b 6.16 "Thermostatting, barostatting, and compute temperature"_howto_16 :b 6.17 "Walls"_howto_17 :b diff --git a/doc/Section_example.html b/doc/Section_example.html index 94e2d95d64..dfc19906dc 100644 --- a/doc/Section_example.html +++ b/doc/Section_example.html @@ -31,15 +31,21 @@ Site.

These are the sample problems in the examples sub-directories:

+ + + + + + @@ -51,11 +57,13 @@ Site. - + - + +
body body particles, 2d system
colloid big colloid particles in a small particle solvent, 2d system
comb models using the COMB potential
crack crack propagation in a 2d solid
dipole point dipolar particles, 2d system
dreiding methanol via Dreiding FF
eim NaCl using the EIM potential
ellipse ellipsoidal particles in spherical solvent, 2d system
flow Couette and Poiseuille flow in a 2d channel
friction frictional contact of spherical asperities between 2d surfaces
gpu use of the GPU package for GPU acceleration
hugoniostat Hugoniostat shock dynamics
indent spherical indenter into a 2d solid
kim use of potentials in Knowledge Base for Interatomic Models (KIM)
line line segment particles in 2d rigid bodies
meam MEAM test for SiC and shear (same as shear examples)
melt rapid melt of 3d LJ system
micelle self-assembly of small lipid-like molecules into 2d bilayers
peptide dynamics of a small solvated peptide chain (5-mer)
peri Peridynamic model of cylinder impacted by indenter
pour pouring of granular particles into a 3d box, then chute flow
prd parallel replica dynamics of a vacancy diffusion in bulk Si
prd parallel replica dynamics of vacancy diffusion in bulk Si
reax RDX and TATB models using the ReaxFF
rigid rigid bodies modeled as independent or coupled
shear sideways shear applied to 2d solid, with and without a void
srd stochastic rotation dynamics (SRD) particles as solvent +
srd stochastic rotation dynamics (SRD) particles as solvent
tad temperature-accelerated dynamics of vacancy diffusion in bulk Si
tri triangular particles in rigid bodies

Here is how you might run and visualize one of the sample problems: diff --git a/doc/Section_example.txt b/doc/Section_example.txt index 947c7957c7..21b5aa8363 100644 --- a/doc/Section_example.txt +++ b/doc/Section_example.txt @@ -27,15 +27,21 @@ Site"_lws. These are the sample problems in the examples sub-directories: +body: body particles, 2d system colloid: big colloid particles in a small particle solvent, 2d system comb: models using the COMB potential crack: crack propagation in a 2d solid dipole: point dipolar particles, 2d system +dreiding: methanol via Dreiding FF eim: NaCl using the EIM potential ellipse: ellipsoidal particles in spherical solvent, 2d system flow: Couette and Poiseuille flow in a 2d channel friction: frictional contact of spherical asperities between 2d surfaces +gpu: use of the GPU package for GPU acceleration +hugoniostat: Hugoniostat shock dynamics indent: spherical indenter into a 2d solid +kim: use of potentials in Knowledge Base for Interatomic Models (KIM) +line: line segment particles in 2d rigid bodies meam: MEAM test for SiC and shear (same as shear examples) melt: rapid melt of 3d LJ system micelle: self-assembly of small lipid-like molecules into 2d bilayers @@ -47,11 +53,13 @@ obstacle: flow around two voids in a 2d channel peptide: dynamics of a small solvated peptide chain (5-mer) peri: Peridynamic model of cylinder impacted by indenter pour: pouring of granular particles into a 3d box, then chute flow -prd: parallel replica dynamics of a vacancy diffusion in bulk Si +prd: parallel replica dynamics of vacancy diffusion in bulk Si reax: RDX and TATB models using the ReaxFF rigid: rigid bodies modeled as independent or coupled shear: sideways shear applied to 2d solid, with and without a void -srd: stochastic rotation dynamics (SRD) particles as solvent :tb(s=:) +srd: stochastic rotation dynamics (SRD) particles as solvent +tad: temperature-accelerated dynamics of vacancy diffusion in bulk Si +tri: triangular particles in rigid bodies :tb(s=:) Here is how you might run and visualize one of the sample problems: diff --git a/doc/Section_howto.html b/doc/Section_howto.html index 04c9dc6812..385c605a0e 100644 --- a/doc/Section_howto.html +++ b/doc/Section_howto.html @@ -26,7 +26,7 @@ 6.11 Visualizing LAMMPS snapshots
6.12 Triclinic (non-orthogonal) simulation boxes
6.13 NEMD simulations
-6.14 Extended spherical and aspherical particles
+6.14 Finite-size spherical and aspherical particles
6.15 Output from LAMMPS (thermo, dumps, computes, fixes, variables)
6.16 Thermostatting, barostatting and computing temperature
6.17 Walls
@@ -163,7 +163,7 @@ so that any forces induced by other fixes will be zeroed out.

Many of the example input scripts included in the LAMMPS distribution are for 2d models.

-

IMPORTANT NOTE: Some models in LAMMPS treat particles as extended +

IMPORTANT NOTE: Some models in LAMMPS treat particles as finite-size spheres, as opposed to point particles. In 2d, the particles will still be spheres, not disks, meaning their moment of inertia will be the same as in 3d. @@ -1012,7 +1012,7 @@ profile consistent with the applied shear strain rate.


-

6.14 Extended spherical and aspherical particles +

6.14 Finite-size spherical and aspherical particles

Typical MD models treat atoms or particles as point masses. Sometimes it is desirable to have a model with finite-size particles such as @@ -1028,8 +1028,12 @@ particles. The following aspects are discussed in turn:

  • pair potentials
  • time integration
  • computes, thermodynamics, and dump output -
  • rigid bodies composed of extended particles +
  • rigid bodies composed of finite-size particles +

    Example input scripts for these kinds of models are in the body, +colloid, dipole, ellipse, line, peri, pour, and tri directories of the +examples directory in the LAMMPS distribution. +

    Atom styles

    There are several atom styles that allow for @@ -1042,13 +1046,14 @@ particles store an angular velocity (omega) and can be acted upon by torque. The "set" command can be used to modify the diameter and mass of individual particles, after then are created.

    -

    The dipole style does not actually define extended particles, but is -often used in conjunction with spherical particles, via a command like +

    The dipole style does not actually define finite-size particles, but +is often used in conjunction with spherical particles, via a command +like

    atom_style hybrid sphere dipole 
     

    This is because when dipoles interact with each other, they induce -torques, and a particle must be extended (i.e. have a moment of +torques, and a particle must be finite-size (i.e. have a moment of inertia) in order to respond and rotate. See the atom_style dipole command for details. The "set" command can be used to modify the orientation and length of the dipole moment of @@ -1094,30 +1099,29 @@ diameter is set to 0.0, it will be a point particle. In the line or tri style, if the lineflag or triflag is specified as 0, then it will be a point particle.

    -

    Many of the pair styles used to compute pairwise interactions between -extended particles typically compute the correct interaction in these -simplified (cheaper) cases. e.g. the interaction between a point -particle and an extended particle or between two point particles. If -necessary, pair_style hybrid can be used to insure -the correct interactions are computed for the appropriate style of -interactions. Likewise, using groups to partition particles -(ellipsoids versus spheres versus point particles) will allow you to -use the appropriate time integrators and temperature computations for -each class of particles. See the doc pages for various commands for -details. +

    Some of the pair styles used to compute pairwise interactions between +finite-size particles also compute the correct interaction with point +particles as well, e.g. the interaction between a point particle and a +finite-size particle or between two point particles. If necessary, +pair_style hybrid can be used to insure the correct +interactions are computed for the appropriate style of interactions. +Likewise, using groups to partition particles (ellipsoids versus +spheres versus point particles) will allow you to use the appropriate +time integrators and temperature computations for each class of +particles. See the doc pages for various commands for details.

    Also note that for 2d simulations, atom styles sphere and ellipsoid still use 3d particles, rather than as circular disks or -ellipses. This means they have the same moment of inertia as a 3d -extended object. When temperature is computed, the correct degrees of -freedom are used for rotation in a 2d versus 3d system. +ellipses. This means they have the same moment of inertia as the 3d +object. When temperature is computed, the correct degrees of freedom +are used for rotation in a 2d versus 3d system.

    Pair potentials
    -

    When a system with extended particles is defined, the particles will -only rotate and experience torque if the force field computes such -interactions. These are the various pair styles -that generate torque: +

    When a system with finite-size particles is defined, the particles +will only rotate and experience torque if the force field computes +such interactions. These are the various pair +styles that generate torque:

    -

    Note that for mixtures of point and extended particles, these +

    Note that for mixtures of point and finite-size particles, these integration fixes can only be used with groups which -contain extended particles. +contain finite-size particles.

    Computes, thermodynamics, and dump output
    @@ -1188,15 +1192,15 @@ rotational energy of spherical or ellipsoidal particles:

    These include rotational degrees of freedom in their computation. If you wish the thermodynamic output of temperature or pressure to use -one of these computes (e.g. for a system entirely composed of extended -particles), then the compute can be defined and the +one of these computes (e.g. for a system entirely composed of +finite-size particles), then the compute can be defined and the thermo_modify command used. Note that by default thermodynamic quantities will be calculated with a temperature that only includes translational degrees of freedom. See the thermo_style command for details.

    -

    These commands can be used to output various attributes -of extended particles: +

    These commands can be used to output various attributes of finite-size +particles: