TY - JOUR
T1 - Can classical equations simulate quantum‐mechanical behavior? a molecular dynamics investigation of a diatomic molecule with a morse potential
AU - Schlick, Tamar
AU - Peskin, Charles S.
PY - 1989/12
Y1 - 1989/12
N2 - In a recent paper, we presented a new computational method for molecular dynamics which uses the Backward‐Euler scheme to solve the classical Langevin dynamics equations. Parameters for the simulation include a target temperature T, a time step Δt, and a cutoff frequency ωc. We showed for a harmonic oscillator system that the cutoff frequency can be set as ωc = kT/h in order to mimic quantum‐mechanical behavior. We now continue this investigation for a nonlinear case: a diatomic molecule governed by a Morse bond potential. Since approximate quantum‐mechanical energy levels are explicitly known for this model, a comparison of energies can be made with molecular dynamics results. By performing dynamics runs for a wide range of temperatures and calculating mean energies, we find a very good agreement between these energies and quantum mechanical predictions. Vibrational excitation begins at temperatures around 800 K, and for higher temperatures both energy curves (molecular dynamics and quantum mechanics) approach the classical prediction of 7/2kT energy per molecule. Future investigations will focus on more general nonlinear potential functions employed in force fields of nucleic acids and proteins.
AB - In a recent paper, we presented a new computational method for molecular dynamics which uses the Backward‐Euler scheme to solve the classical Langevin dynamics equations. Parameters for the simulation include a target temperature T, a time step Δt, and a cutoff frequency ωc. We showed for a harmonic oscillator system that the cutoff frequency can be set as ωc = kT/h in order to mimic quantum‐mechanical behavior. We now continue this investigation for a nonlinear case: a diatomic molecule governed by a Morse bond potential. Since approximate quantum‐mechanical energy levels are explicitly known for this model, a comparison of energies can be made with molecular dynamics results. By performing dynamics runs for a wide range of temperatures and calculating mean energies, we find a very good agreement between these energies and quantum mechanical predictions. Vibrational excitation begins at temperatures around 800 K, and for higher temperatures both energy curves (molecular dynamics and quantum mechanics) approach the classical prediction of 7/2kT energy per molecule. Future investigations will focus on more general nonlinear potential functions employed in force fields of nucleic acids and proteins.
UR - http://www.scopus.com/inward/record.url?scp=84990622745&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84990622745&partnerID=8YFLogxK
U2 - 10.1002/cpa.3160420807
DO - 10.1002/cpa.3160420807
M3 - Article
AN - SCOPUS:84990622745
SN - 0010-3640
VL - 42
SP - 1141
EP - 1163
JO - Communications on Pure and Applied Mathematics
JF - Communications on Pure and Applied Mathematics
IS - 8
ER -