Syntax:
fix ID group-ID rigid keyword values
single values = none molecule values = none group values = list of group IDs
Examples:
fix 1 clump rigid single fix 1 polychains rigid molecule fix 2 fluid rigid group clump1 clump2 clump3
Description:
Treat one or more sets of atoms as an independent rigid body. This means that each timestep the total force and torque on each rigid body is computed and the coordinates and velocities of the atoms in each body are updated so that they move as a rigid body. This can be useful for freezing one or more portions of a large biomolecule, or for simulating a system of colloidal particles.
This fix updates the positions and velocities of the rigid atoms with a constant-energy time integration, so you should not update the same atoms via other fixes (e.g. nve, nvt, npt, temp/rescale, langevin).
Each body must have two or more atoms. Which atoms are in which bodies can be defined via several options.
For option single the entire group of atoms is treated as one rigid body.
For option molecule, each set of atoms in the group with a different molecule ID is treated as a rigid body.
For option group, each of the listed groups is treated as a separate rigid body. Note that only atoms that are also in the fix group are included in each rigid body.
For computational efficiency, you should also turn off pairwise and bond interactions within each rigid body, as they no longer contribute to the motion. The neigh_modify exclude and delete_bonds commands are used to do this.
For computational efficiency, you should define one fix rigid which includes all the desired rigid bodies. LAMMPS will allow multiple rigid fixes to be defined, but it is more expensive.
The degrees-of-freedom removed by rigid bodies are accounted for in temperature and pressure computations. Similary, the rigid body contribution to the pressure virial is also accounted for. The latter is only correct if forces within the bodies have been turned off, and there is only a single fix rigid defined. For each linear rigid body of three or more atoms, one degree-of-freedom must be added using a temp_modify command (i.e. for a simulation of 10 such rigid bodies, use "temp_modify mine extra -7", where 3 is the default setting and -10 degrees-of-freedom are subtracted).
Note that this fix uses constant-energy integration, so you need to impose additional constraints to control the temperature of an ensemble of rigid bodies. You can use fix langevin for this purpose to treat the system as effectively immersed in an implicit solvent, i.e. a Brownian dynamics model. Or you can thermostat the atoms of an explicit solvent directly.
Restrictions:
This fix performs an MPI_Allreduce each timestep that is proportional in length to the number of rigid bodies. Hence it will not scale well in parallel if large numbers of rigid bodies are simulated.
If the atoms in a single rigid body initially straddle a periodic boundary, the input data file must define the image flags for each atom correctly, so that LAMMPS can "unwrap" the atoms into a valid rigid body.
You should not use this fix if you just want to hold group of atoms stationary. A better way to do this is to not include those atoms in your time integration fix. E.g. use "fix 1 mobile nve" instead of "fix 1 all nve", where "mobile" is the group of atoms that you want to move.
Related commands:
delete_bonds, neigh_modify exclude
Default: none