Abstract
Following the previous theoretical developments to completely reproduce the IR spectra of weak hydrogen bond complexes within the framework of the linear response theory (LRT), the quantum theory of the high stretching mode spectral density (SD) of weak H-bonds is reconsidered. Within the LRT theory, the SD is the one sided Fourier transform of the autocorrelation function (ACF) of the high stretching mode dipole moment operator. In order to provide more accurate theoretical bandshapes, we have explored the equivalence between the SDs given in previous studies with respect to a new quantum one, and revealed that in place of the basic equations used in the precedent works for which the SD l(old)(omega) = 2 Re integral(infinity)(0) G(old) (t)e-(i omega t)dt where the ACF G(old)(t) (mu(0)mu(t)(+) = tr{rho{(0)}{mu(t)}(+)), one can use a new expression for the SD, given by 1(New)(omega)) = 2 omega Re integral(infinity)(0) G(New)(t)e-(I omega t)dt where G(New)(t) = <mu(0)mu(t)(+))> = 1/beta tr {rho(B) integral(beta)(0) {mu(0)} {mu(t + i lambda h)}(+) d lambda}. Here rho(B) is the Boltzmann density operator, mu(0) the dipole moment operator at initial time and ki(t) the dipole moment operator at time tin the Heisenberg picture, h is the Planck constant, beta is the inverse of the Boltzmann factor k(B)T where T is the absolute temperature and k(B) the Boltzmann constant. Using this formalism, we demonstrated that the new quantum approach gives the same final SD as used by previous models, and reduces to the Franck-Condon progression appearing in the Marechal and Witkowski's pioneering approach when the relaxation mechanisms are ignored. Results of this approach shed light on the equivalence between the quantum and classical IR SD approaches for weak H-bonds in absence of medium surroundings effect, which has been a subject of debate for decades. (C) 2018 Elsevier B.V. All rights reserved.