TY - JOUR
T1 - Extended ensemble molecular dynamics method for constant strain rate uniaxial deformation of polymer systems
AU - Yang, Liu
AU - Srolovitz, David J.
AU - Yee, Albert F.
PY - 1997/9/15
Y1 - 1997/9/15
N2 - We describe a novel molecular dynamics (MD) method to simulate the uniaxial deformation of an amorphous polymer. This method is based on a rigorously defined statistical mechanics ensemble appropriate for describing an isothermal, displacement controlled, uniaxial stress mechanical test. The total number of particles is fixed and the normal stresses in the direction normal to the applied strain are constant, i.e., an NTLxσyyσzz ensemble. By using the Lagrangian of the extended system (i.e., including additional variables corresponding to the temperature and cross-sectional area fluctuations), we derive a set of equations of motion for the atomic coordinates and the additional variables appropriate to this ensemble. In order to avoid the short MD time step appropriate for the stiff covalent bonds along the polymer chains, we introduce bond length constraints. This is achieved using a variation of the commonly used SHAKE [J. P. Ryckaert, G. Ciccotti, and H. J. C. Berendsen, J. Comp. Phys. 23, 327 (1977)] algorithm. A numerical method for integrating the equations of motion with constraints via a modification of the velocity Verlet [W. C. Swope, H. C. Andersen, P. H. Berens, and K. R. Wilson, J. Chem. Phys. 76, 637 (1982)] algorithm is presented. We apply this new algorithm to the constant strain rate deformation of an amorphous polyethylene in a model containing several distinct polymer chains. To our knowledge, this is the first time that bond length constraints were applied to a macromolecular system together with an extended ensemble in which the simulation cell shape is allowed to fluctuate.
AB - We describe a novel molecular dynamics (MD) method to simulate the uniaxial deformation of an amorphous polymer. This method is based on a rigorously defined statistical mechanics ensemble appropriate for describing an isothermal, displacement controlled, uniaxial stress mechanical test. The total number of particles is fixed and the normal stresses in the direction normal to the applied strain are constant, i.e., an NTLxσyyσzz ensemble. By using the Lagrangian of the extended system (i.e., including additional variables corresponding to the temperature and cross-sectional area fluctuations), we derive a set of equations of motion for the atomic coordinates and the additional variables appropriate to this ensemble. In order to avoid the short MD time step appropriate for the stiff covalent bonds along the polymer chains, we introduce bond length constraints. This is achieved using a variation of the commonly used SHAKE [J. P. Ryckaert, G. Ciccotti, and H. J. C. Berendsen, J. Comp. Phys. 23, 327 (1977)] algorithm. A numerical method for integrating the equations of motion with constraints via a modification of the velocity Verlet [W. C. Swope, H. C. Andersen, P. H. Berens, and K. R. Wilson, J. Chem. Phys. 76, 637 (1982)] algorithm is presented. We apply this new algorithm to the constant strain rate deformation of an amorphous polyethylene in a model containing several distinct polymer chains. To our knowledge, this is the first time that bond length constraints were applied to a macromolecular system together with an extended ensemble in which the simulation cell shape is allowed to fluctuate.
UR - http://www.scopus.com/inward/record.url?scp=0031571941&partnerID=8YFLogxK
UR - https://www.scopus.com/record/pubmetrics.uri?eid=2-s2.0-0031571941&origin=recordpage
U2 - 10.1063/1.474781
DO - 10.1063/1.474781
M3 - RGC 21 - Publication in refereed journal
VL - 107
SP - 4396
EP - 4407
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
SN - 0021-9606
IS - 11
ER -