Molecular dynamics simulation of planar elongational flow at constant pressure and constant temperature

Molecular dynamics simulations of liquid systems under planar elongational flow have mainly been performed in the NVT ensemble. However, in most material processing techniques and common experimental settings, at least one surface of the fluid is kept in contact with the atmosphere, thus maintaining the sample in the NpT ensemble. For this reason, an implementation of the Nosé-Hoover integral-feedback mechanism for constant pressure is presented, implemented via the SLLOD algorithm for elongational flow. The authors test their procedure for an atomic liquid and compare the viscosity obtained with that in the NVT ensemble. The scheme is easy to implement, self-starting and reliable, and can be a useful tool for the simulation of more complex liquid systems, such as polymer melts and solutions.