next up previous contents
Next: Constraint Methods Up: Mechanics Previous: Constant-pressure integration with a   Contents


Constant pressure in 3 axes (NPTxyz)

There is one additional extended system integrator which is somewhat simpler than the NPTf method described above. In this case, the three axes have independent barostats which each attempt to preserve the target pressure along the box walls perpendicular to that particular axis. The lengths of the box axes are allowed to fluctuate independently, but the angle between the box axes does not change. The equations of motion are identical to those described above, but only the diagonal elements of $ \overleftrightarrow{\eta}$ are computed. The off-diagonal elements are set to zero (even when the pressure tensor has non-zero off-diagonal elements).

It should be noted that the NPTxyz integrator is not known to preserve any Hamiltonian of interest to the chemical physics community. The integrator is extremely useful, however, in generating initial conditions for other integration methods. It is suitable for use with liquid simulations, or in cases where there is orientational anisotropy in the system (i.e. in lipid bilayer simulations).



Copyright © 2006 - University of Notre Dame

Updated on January 16, 2006