Vibrational frequencies in Car-Parrinello molecular dynamics
Sheau Wei Ong1, Eng Soon Tok, Hway Chuan Kang
1Department of Chemistry, National University of Singapore, Singapore.
Abstract:
Car-Parrinello molecular dynamics (CPMD) are widely used to investigate the dynamical properties of molecular systems. An important issue in such applications is the dependence of dynamical quantities such as molecular vibrational frequencies upon the fictitious orbital mass μ. Although it is known that the correct Born-Oppenheimer dynamics are recovered at zero μ, it is not clear how these dynamical quantities are to be rigorously extracted from CPMD calculations. Our work addresses this issue for vibrational frequencies. We show that when the system is sufficiently close to the ground state the calculated ionic vibrational frequencies are ω(M) = ω(0M)[1 -C(μ/M)] for small μ/M, where ω(0M) is the Born-Oppenheimer ionic frequency, M the ionic mass, and C a constant that depends upon the ion-orbital coupling force constants. Our analysis also provides a quantitative understanding of the orbital oscillation amplitudes, leading to a relationship between the adiabaticity of a system and the ion-orbital coupling constants. In particular, we show that there is a significant systematic dependence of calculated vibrational frequencies upon how close the CPMD trajectory is to the Born-Oppenheimer surface. We verify our analytical results with numerical simulations for N(2), Sn(2), and H/Si(100)-(2×1).
More Related Videos
Related Concept Videos
IR Spectroscopy: Hooke's Law Approximation of Molecular Vibration
According to Hooke's law, the vibrational frequency is directly proportional to the...
IR Spectroscopy: Molecular Vibration Overview
Stretching vibrations are vibrational motions that occur along the bond line, changing the bond length or distance between two bonded atoms. They are further distinguished as symmetric or asymmetric. In symmetric stretching, the...
¹H NMR: Interpreting Distorted and Overlapping Signals
As Δν decreases and the signals move closer, the doublets appear increasingly distorted. The intensities of the inner lines increase at the cost of those of the outer lines as the signals are slanted or...
Spin–Spin Coupling Constant: Overview
Qualitatively, any spin plus-half nucleus polarizes the spins of its electrons to the minus-half state. Consequently, the paired electron in the hydrogen–carbon bond must have a...
¹H NMR of Conformationally Flexible Molecules: Temporal Resolution
Distribution of Molecular Speeds


