Phonostat: Thermostatting phonons in molecular dynamics simulations
Author(s)Raghunathan, Rajamani; Greaney, P. Alex; Grossman, Jeffrey C.
MetadataShow full item record
Thermostat algorithms in a molecular dynamics simulation maintain an average temperature of a system by regulating the atomic velocities rather than the internal degrees of freedom. Herein, we present a “phonostat” algorithm that can regulate the total energy in a given internal degree of freedom. In this algorithm, the modal energies are computed at each time step using a mode-tracking scheme and then the system is driven by an external driving force of desired frequency and amplitude. The rate and amount of energy exchange between the phonostat and the system is controlled by two distinct damping parameters. Two different schemes for controlling the external driving force amplitude are also presented. In order to test our algorithm, the method is applied initially to a simple anharmonic oscillator for which the role of various phonostat parameters can be carefully tested. We then apply the phonostat to a more realistic (10,0) carbon nanotube system and show how such an approach can be used to regulate energy of highly anharmonic modes.
DepartmentMassachusetts Institute of Technology. Department of Materials Science and Engineering
The Journal of Chemical Physics
American Institute of Physics (AIP)
Raghunathan, Rajamani, P. Alex Greaney, and Jeffrey C. Grossman. “Phonostat: Thermostatting phonons in molecular dynamics simulations.” The Journal of Chemical Physics 134, no. 21 (2011): 214117. © 2011 American Institute of Physics
Final published version