We propose a numerical
technique to accurately simulate the vibrations
of organic molecules in the gas phase, when pairs of atoms (or, in
general, groups of degrees of freedom) are artificially decoupled,
so that their motion is instantaneously decorrelated. The numerical
technique we have developed is a symplectic integration algorithm
that never requires computation of the force but requires estimates
of the Hessian matrix. The theory we present to support our technique
postulates a pair-decoupling Hamiltonian function, which parametrically
depends on a decoupling coefficient α ∈ [0, 1]. The closer
α is to 0, the more decoupled the selected atoms. We test the
correctness of our numerical method on small molecular systems, and
we apply it to study the vibrational spectroscopic features of salicylic
acid at the Density Functional Theory
ab initio
level
on a fitted potential. Our pair-decoupled simulations of salicylic
acid show that decoupling hydrogen-bonded atoms do not significantly
influence the frequencies of stretching modes, but enhance enormously
the out-of-plane wagging and twisting motions of the hydroxyl and
carboxyl groups to the point that the carboxyl and hydroxyl groups
may overcome high potential energy barriers and change the salicylic
acid conformation after a short simulation time. In addition, we found
that the acidity of salicylic acid is more influenced by the dynamical
couplings of the proton of the carboxylic group with the carbon ring
than with the hydroxyl group.