3.2. Autocorrelation Functions and Spectral Density

The spectral density of the υ_{N-H} mode is given, within the linear response theory [27,28], by the Fourier transform of the autocorrelation function G(t) of transition moment operator of the fast mode :

(14)

Using the symmetry of the system, the ACF may be split into symmetric parts (g) and antisymmetric parts (u).

In the presence of Davydov coupling, the autocorrelation function (ACF) of the dipole moment operator of the fast mode may be written [25]:

(15)

Here, γ is the natural width of the excited state of the high frequency mode, the expression of which has been calculated by Rösh and Ratner [22]. [G (t)]_{g} is the reduced ACF of the g part of the system related to Hamiltonian described the indirect damping. It’s given by:

(16)

In these equation,

(17)

(18)

Note that by “reduced”, we mean that a partial trace has to be performed over the thermal bath because of the coupling between the symmetric Hamiltonian and the surrounding. In the last equation, [G^{+}(t)]_{u} and [G^{–}(t)]_{u} are the two (u) ACFs which are affected only by the Davydov coupling. They are given by [29]:

(19)

where, are the expansion coefficients of the eigenvectors on the basis of the eigenstates of the Hamiltonian of the quantum harmonic oscillator. One may observe that the angular frequency ω˚ of the high frequency mode must be decreased by a factor on D isotopic substitution of the proton involved in the H-bonds. Besides, according to Maréchal and Witkowski theory [21], the anharmonic coupling parameter α˚ must also be reduced by this factor upon this same substitution, whereas the frequency Ω of the H-bond bridge has no reason to be modified.

Following Equations (15) and (19), the spectral density of the dimer involving Davydov effect can be written as:

(20)

That may be written formally:

(21)

with respectively:

(22)

In view of the above equations, the components appearing here are given by the following equations [29]:

(23)

(24)

(25)

(26)

(27)

Recall that is the thermal average of the number occupation operator of the H-bond bridge vibration given by Equation (18), which is a function of the angular frequency Ω of the H-bond bridge and of the absolute temperature T. Besides, β is the effective anharmonic coupling parameter given by Equation (17), which is a function of the anharmonic coupling parameter α˚ between the slow and fast mode, and of the angular frequency Ω and the damping parameter γ of the H-bond bridge.

3.3. Situation with Fermi Resonances

The previous treatment was developed with the neglect of Fermi resonances. Now, suppose the situation where this later effect is taken into account. It is resulting from the interactions occurring between the first excited state of the high frequency mode and the first harmonic of some bending modes. As it was stated by Maréchal and Witkowski [21], if Fermi resonances are taken into account, one has to consider one fast mode, one slow mode and one bending mode for each hydrogen bond of the cyclic dimer. If we take into account Fermi resonances [29], they affect only the g state of the system. As a consequence, the autocorrelation functions [G^{–}(t)]_{u}_{ }are not modified. In the presence of Davydov coupling and Fermi resonances, the ACF can be written following [29]:

(28)

Here, {ω_{μ}}_{g} are the eigenvalues appearing in Equation (29) whereas, a_{{μ,0,m}g} are the expansion coefficients defined by Equation (30).

(29)

With

(30)

The g states involved in the above expansions are defined by:

(31)

(32)

(33)

Here, the kets

, and

are the g eigenstates of respectively the symmetrized high frequency quantum harmonic oscillator, the slow frequency quantum harmonic oscillator and the bending frequency quantum harmonic oscillator. The Fermi resonance mechanism, characterized by the coupling parameter f_{i}, is described by the following coupling operators ℏf which express the non-resonant exchanges between the state of the jth fast mode and second damped excited state of the jth bending mode. The introduction of the Fermi resonance coupling effects in the lineshape is presented by the complexes angular frequency gap Δ_{i}:

(34)

with

(35)

where is the frequency of the bending modes.. The imaginary part in this gap is related to the lifetime of the corresponding excited states. Recall that the Hamiltonian of the dimer involving Davydov coupling, Fermi resonances between the g excited state of the fast mode and the g first harmonics of the bending mode, together with the damping of theses excited states is:

(36)

with

(37)

and

(38)

In these two last equations, a^{+} and a are the boson operators obeying [a^{+}, a] = 1 and i^{2} = –1.

4. Results and Discussions

In our theoretical approach, we have assumed that the 2-pyridone dimer crystal in alpha-phase is of approximately centrosymmetric type. In the case of the interprettation of the spectra of alpha-phase 2-pyridone crystal, this interpretation cannot be relied on since the N-H…O bonds are much shorter than the sulfur atoms containing hydrogen bonds. Most probably, the purely vibrational approach to solving the problem of the exciton interacttions in the hydrogen bond dimers should be abandoned. However, the contribution of electronic interactions to the vibrational exciton interactions between the hydrogen bonds in the dimers, expressed by the electronic coordinates, should be taken account.

From these considerations some consequences result for the electronic structure of 2-pyridone dimers in alpha-phase crystals. This is most probably the structure of a larger contribution of the zwitterion electronic structure. Such a continuous structure of the π-electron cloud in the pyridine rings, in cyclic N^{+}-H…O^{–} hydrogen bond dimers allows for an effective (“head-to-tail”) coupling of the hydrogen bonds in the cyclic dimers. The experimental and theoretical results allow us to state that, for alpha-phase 2-pyridone studied the υ_{N-H} absorption bands of their dimers observed in the IR spectra are similar in structure and are very broad. This means that the mechanism of formation of these bands should, first of all, involve the participation of the cyclic structure, responsible for the formation of the hydrogen bond. Insignificant changes in the shape of the dimers and their deuterated derivatives bands observed in alpha-phase 2-pyridone. We believe that the interaction of the two intermolecular bonds N-H…O in the cyclic complex plays a determining role in the formation of the broad absorption band.

The main object in this work is to reproduce theoretically the IR spectra of alpha-phase 2-pyridone at 300 K. In a first place we have reproduced the υ_{N-H(D)} IR spectra of alpha-phase 2-pyridone dimer crystal by comparing theoretical line shapes. The spectral densities are computed by Equation (19) after construct and diagonalize the full Hamiltonian in a truncated basis. The stability of the computed spectra with respect to the size of this basis set was carefully checked. The stability of the spectra was also checked versus the order of the Taylor expansion of the Morse potential. Thus, in the first place, we have considered the general physical situation in which multiple Fermi resonances are taken into account and for which the spectral density is given by Equations (19), (20) and (28). The procedure we have used is the fitting of the experimental line shapes by optimizing the values of the basic parameters.

Figures 2 and 3 present the theoretical spectra of alpha-phase 2-pyridone dimer crystal and its deuterated derivative where multiple Fermi resonances have been introduced in our theoretical approach. We have performed numerical experimentation by increasing progressively the number of Fermi resonances. The spectral densities, presented in this Figure, are computed in the presence of three Fermi resonances. Tables 1 and 2 show the parameters involved in the calculations.

The examination of the spectra obtained using our theoretical model shows that there is a good agreement between theory and experiment. It is important to emphasize that if Fermi resonances are acting, one observes the classical behavior of the Fermi resonance, which is well described in the literature: the displacement of the bands and the redistribution of the intensities in the bands. However, only the intensities at high frequencies are redistributed if one increases the number of Fermi resonances. In fact, a recent study on acetic acid dimers in the liquid phase [30] did not show any wavelength dependence of their dynamics.

Figures 4 and 5 show the theoretical line shapes of alpha-phase 2-pyridone dimer crystal and its deuterated derivative calculated by Equations (16), (19) and (20) when Fermi resonances are ignored. In both situations, the theoretical line shapes appear as red lines whereas the experimental ones are black lines.

Figure 3. Comparison between the experimental and theoretical spectra of deuterium-bonded in alpha-phase 2-pyridone at room temperature in the presence of three Fermi resonances.

Table 1. Parameters used for fitting the experimental spectra of alpha-phase 2-pyridone at room temperature in the presence of three Fermi resonances.

Table 2. Fermi coupling parameters (in cm^{–}^{1}) used for fitting experimental spectra of phase 2-pyridone at room temperature.

Figure 4. Comparison between the experimental and theoretical spectra of hydrogen-bonded in alpha-phase 2-pyridone at room temperature when the Fermi resonances are ignored.

Figure 5. Comparison between the experimental and theoretical spectra of deuterium-bonded in alpha-phase 2-pyridone at room temperature in the absence of Fermi resonances.

In order to improve the validity of the present model, let us comment the magnitude of the parameters that we have used. Both the anharmonic coupling parameter and the angular frequencies of the high frequency mode for alpha-phase 2-pyridone in all treated situations, their magnitudes is decreased by a factor which is different from. This result is in agreement with theory when passing from the H to the D species and when the slow mode is assumed to be of Morse type (anharmonic potential). Indeed, according to [30,31], the isotope ratio

for hydrogen bonds O-H…O depends on the bond strength, it is close to the harmonic value for weak bonds and decreases with the bond strength reaching minimum near for low-barrier bonds. In our cases, to obtain good agreement between the theoretical and experimental one, we have used:

(39)

It is important to note that the introduction modulation of the equilibrium positions of the fast modes (q_{e}) and the quadratic dependence of their frequencies on Q_{i} which can represented by two expansions to the second order:

(40)

Recall that in present work and following Maréchal and Witkowski [21], we have only used the first order dependence of the angular frequency of the fast modes on the coordinate of the slow modes (Q_{i}) (see Equation (40) and we have neglected the modulation of the equilibrium positions of the fast modes and the quadratic dependence of their frequencies on (Q_{i}) [32] as conformed by several experimental correlations [33,34]. We have shown in recent [32] work that the fine structure of the IR υ_{X-H} stretching band are connected with the new anharmonic coupling parameters (i.e. β, f˚ and g˚), whereas these parameters do not manifest markedly in effects of temperature and deuteration on the IR spectra. Besides, the account for all these parameters does not affect the similarity of the spectra in gas and condensed phases [35].

The isotope ratio χ of centers of gravity of band of light and deuterated alpha-phase 2-pyridone dimer is 1.306. The values obtained here are also in satisfactory agreement with the experimental data reported by Odinokov, et al. [36] for complexes of carboxylic acids with various bases and are comparable to those found in our recent work dealing with the H/D isotopic effects in H-bond spectra [33]. However, these ratios are different from those used by Blaise, et al. [37] since in their approaches dealing with theoretical interpretation of the IR line shapes of liquid and gaseous acetic acid [38] and gaseous propynoic and acrylic acid dimers [29], low and high-frequency hydrogen stretching vibrations in individual hydrogen bonds are assumed to be harmonic whereas in the present work we use a Morse potential in order to describe the anharmonicity of the H-bond Bridge. Recall that the removal of the harmonic approximation for the slow modes by introducing Morse potential in place of harmonic one has been done by Leviel and Maréchal [39] in a model similar to the present one, involving cyclic dimer, however without damping. They have shown that the value of the angular frequency of the slow mode which must be used to the experimental lineshape is more close to the experimental value when the anharmonicity of the bridge is introduced.

Now let us look to the strongest hydrogen bonds observed in alpha-phase 2-pyridone dimers. Generally this type of hydrogen bonds is observed in ion molecular complexes of the (AHA)^{–} and (BHB)^{+} types [40,41]; however, it is especially difficult to study the spectral properties of charged complexes under conditions of weak interactions with the medium. In the present approach, we take into account the natural width of the excited states of the fast mode due to the medium (direct relaxation) and the damping of the H-bond Bridge (indirect relaxation). The direct relaxation is included following quantum treatment of Rösch and Ratner [22] whereas the indirect one was taken into account via the approach of Louisell and Walker [42] dealing with the relaxation of driven damped quantum oscillators studied initially by Feynman and Vernon [43] and later by Louisell [44]. The values of the direct and indirect relaxation parameters reflecting the effect of the medium used presently for alpha-phase 2-pyridone are of the same magnitude as those used by Blaise, et al. [29,37,38] in their study dealing with acetic acid in the gas phase, whereas the indirect damping at 298 K for acetic acid in crystalline state is γ = 1 cm^{–1}. One may ask whether the indirect damping used for crystalline state is weaker than that used for the gaseous phase since the indirect relaxation is ought to be larger in the solid state. The reason is that in the gas phase the indirect damping is an effective one which is the result of the combination of the indirect damping and of the rotational structure [45]. In neutral systems, the strongest hydrogen bonds of the N-H…O type are formed in self-associates of alpha-phase 2-pyridone and, seemingly, phosphinic acids [46], and it is only due to the high thermal stability of these associates that one is able to observe their IR spectra in the gas phase, at which cyclic dimers are in equilibrium with monomeric molecules [47,48]. This made it possible to analyze the dependence of the bandshape of the υ_{N-H} band, which serves as an important criterion in the choice of theoretical models. The parameters of this broad band and its shape turned out to be almost independent of the bonding frequency.

In order to obtain a good agreement with the experimental line shapes, we have take into account some breaking of the IR selection rule for the centrosymmetric cyclic dimer, via a large amount (η = 0.5; 0.9) of forbidden Ag transition. Recall that in general way the quality of the fitting is weakly improved by taking small values for η which lying between 0 and 1. This assumption was initially introduced by Flakus [49]. This is a general trend which has been observed recently in the cases of the centrosymmetric cyclic dimers of gaseous acetic acid [37]. Note that, according to Flakus hypothesis, the lack of forbidden transition ought to be stronger in the solid state (value near η = 1) than in the gaseous one and we must keep in mind that the Flakus assumption has been seen by this author to be unavoidable in diverse crystalline H-bonded carboxylic acids such as for instance glutaric [50] and cinnamic [51] acids and particularly in centrosymmetric H-bonded dimers.

Figure 6 presents the theoretical spectra of alphaphase 2-pyridone dimer where multiple Fermi resonances have been introduced. The number of Fermi resonances n_{F} is increasing from zero (without Fermi resonances) to 3 (n_{F} = 0, 1, 2 and 3). The numerical results show that Fermi resonances appear to play an important role especially for the hydrogenated species whereas the Fermi resonances effect does not affect markedly the N-D derivative species. It may be explain by the fact that Fermi resonances mechanism may involve the N-H bond bending-in-plane vibrations in their first overtone states. So the coupling between the electronic systems with the electrons of the associated carboxyl groups implies a stabilization of the dimers. On the other hand, each deformation of the dimers provides a destruction of this

(a)(b)

Figure 6. Multiple Fermi resonance effect on the IR spectral densities of hydrogen (a) and deuterium (b) bonded in alpha-phase 2-pyridone. The number of Fermi resonances nF is increasing from zero (without Fermi resonances) to 3.

stabilization mechanism. This is the most probable source of the anharmonic coupling involving the proton stretching and the proton bending-in-plane vibrations in the dimer of carboxylic acids. We think that the sensitive improvement of the theoretical line shape of strongly bound dimers of alpha-phase 2-pyridone by introduction of Fermi resonances might be due to some whole effect in which the Fermi resonances assisted by the strong anharmonic coupling would be augmented by the combination of the Davydov coupling and of the quantum direct and indirect damping. It appears from the present theoretical study that the simple model of Davydov coupling taking into account especially Fermi Resonances and quantum direct and indirect damping, is able to reproduce the experimental line shapes. When one of the two identical fast modes is excited, then, because of the symmetry of the cyclic dimer, and of Davydov coupling V_{D} between the two degenerate fast mode excited states, an interaction (Davydov coupling) occurs leading to an exchange between the two identical excited parts of the dimer of alpha-phase 2-pyridone. It is interesting to note that the Davydov coupling parameters used in the present work are similar to those used by Maréchal and Witkowski in their pioneering work dealing with the adipic acid [21]. However, we shall introduce in the Davydov coupling some further flexibility by assuming it to be dependent on the H-bond Bridge coordinate.

The complexity of the task to the experimental line shape of the cyclic dimer of alpha-phase 2-pyridone and to reproduce the corresponding isotope effect is determined by several factors, acting simultaneously. To take them into account, one needs to know a great number of parameters. For this idea, we think that the small residual discrepancies between theory and experiment are related to the neglect of some of these parameters: 1) assumption of a linear dependence of the angular frequency of the fast mode on the coordinate of the H-bond bridge [32], 2) assumption of the independence of the equilibrium position of the fast mode on the coordinate of the H-bond bridge [32], 3) assumption of a constant term for the Davydov coupling parameter [52], 4) neglect of electrical anharmonicity, 5) neglect of eventual weak tunneling through the potential barrier separating the two H-bond bridge minima, and 6) neglect of eventual relaxation mechanisms of non-adiabatic nature.

5. Conclusions

The presented theoretical approach is dealing within the strong anharmonic coupling theory according to which the high frequency mode and the H-bond bridge are anharmonically coupled through a linear dependence of the frequency of the fast mode on the elongation of the Hbond bridge and takes into account Davydov coupling, Fermi resonances, anharmonicity of the H-bond Bridge and direct and indirect quantum relaxations. The present approach contains as special instances the majority precedent theoretical approach [53] dealing with the subject.

This approach, which is applied to reproduce the υ_{N-H(D)} IR spectra of centrosymmetric strongly bound dimers of alpha-phase 2-pyridone and their deuterated derivatives in the solid-state, is comprised in the linear response theory calculates the line shape by aid of the Fourier transform of the autocorrelation function of the dipole moment operator. The model has been applied to alphaphase 2-pyridone dimer crystal. It has been found that it is possible to correctly the experimental line shape of the hydrogenated compound and to predict satisfactorily the deuterium effect by using a set of spectral parameters. It appears from the present theoretical study that the simple model of Davydov coupling taking into account quantum indirect damping and the anharmonicity of the H-bond Bridge, is able to reproduce the experimental line shapes especially when the Fermi resonances were not ignored. In the end, we can say that the Fermi resonances appear play an important role to reproduce the experimental line shapes of alpha-phase 2-pyridone dimer crystal.

6. Acknowledgements

The authors thank Professor Henryk T Flakus from Institute of Chemistry, University of Silesia, 9 Szkolna Street, Pl-40 006 Katowice, Poland for measuring the experimental spectra.

REFERENCES

- G. C. Pimentel and A. L. McClellan, “The Hydrogen Bond,” W. H. Freeman, San Francisco, 1960.
- P. Schuster, G. Zundel and C. Sandorfy, “The Hydrogen Bond: Recent Developments in the Theory and Experiment, Parts I, II and III,” Elsevier, Amsterdam, 1976.
- Y. Marechal, “The Hydrogen Bond and the Water Molecule,” Elsevier, Amsterdam, Oxford, 2006.
- H. D. Arman, P. Poplaukhin and E. R. T. Tiekinkc, “2-Pyridone: Monoclinic Polymorph,” Acta Crystallographica Section E: Structure Reports Online, Vol. 65, No. 12, 2009, p. 3187. doi:10.1107/S1600536809049496
- M. T. Boisdon, S. Castillo, J. F. Brazier, J. Favrot and C. J. Marsden, “2(1H)-Pyridinone(2-Pyridone): Self-Association and Association with Water Spectral and Structural Characteristics: Infrared Study and Ab Initio Calculations,” Spectrochimica Acta Part A, Vol. 59, No. 14, 2003, pp. 3363-3377. doi:10.1016/S1386-1425(03)00162-8
- M. T. Boisdon, S. Castillo, J. F. Brazier and J. Favrot, “Vibrational Study of 2(1H)-Pyridinone(2-pyridone) in H
_{2}O and 1-D-2(1H)-Pyridinone(2-pyridone ND) in D_{2}O,” Spectrochimica Acta Part A, Vol. 55, No. 7-8, 1999, pp. 1379-1388. doi:10.1016/S1386-1425(98)00301-1 - H. Ratajczak and W. J. Orville-Thomas, “Molecular Interactions,” John Wiley & Sons Inc., New York, 1980.
- G. R. Desiraju and T. Steiner, “The Weak Hydrogen Bond in Structural Chemistry and Biology,” Oxford University Press, Oxford, 1999.
- M. Gallant, M. T. P. Viet and J. D. Wuest, “HydrogenBonded Dimers. Direct Study of the Interconversion of Pyridone Dimers and Hydroxypyridine Monomers by Low-Temperature Nuclear Magnetic Resonance Spectroscopy,” Journal of the American Chemical Society, Vol. 113, No. 2, 1991, pp. 721-723. doi:10.1021/ja00002a077
- C. Fuke, T. Arao, Y. Morinaga, H. Takaesu, K. Ameno and T. Miyazaki, “Analysis of Paraquat, Diquat and Two Diquat Metabolites in Biological Materials by High-Performance Liquid Chromatography,” Legal Medicine, Vol. 4, No. 3, 2002, pp. 156-163. doi:10.1016/S1344-6223(02)00011-1
- M. El-Kemary, J. A. Organero and A. Douhal, “Assessment of Solvent Effect on the Relaxation Dynamics of Milrinone, ” Journal of Photochemistry and Photobiology A: Chemistry, Vol. 187, No. 2-3, 2007, pp. 339-347. doi:10.1016/j.jphotochem.2006.10.032
- Y. Ducharme and J. D. Wuest, “Use of Hydrogen Bonds to Control Molecular Aggregation. Extensive, Self-Complementary Arrays of Donors and Acceptors,” Journal of Organic Chemistry, Vol. 53, No. 24, 1988, pp. 5787-5789. doi:10.1021/jo00259a037
- M. Gallant, M. T. P. Viet and J. D. Wuest, “Use of Hydrogen Bonds to Control Molecular Aggregation. Association of Dipyridones Joined by Flexible Spacers,” Journal of Organic Chemistry, Vol. 56, No. 7, 1991 pp. 2284- 2286. doi:10.1021/jo00007a007
- M. J. Wójcik, W. Tatara, M. Boczar, A. Apola and S. Ikeda, “Spectroscopic and Theoretical Study of Vibrational Spectra of Hydrogen-Bonded 2-Pyridone,” Journal of Molecular Structure, Vol. 596, No. 1-3, 2001, pp. 207- 214. doi:10.1016/S0022-2860(01)00719-0
- M. J. Wójcik, “Theoretical Modeling of Vibrational Spectra and Multidimensional Proton Tunneling in Hydrogen-Bonded Systems,” Journal of Molecular Liquids, Vol. 141, No. 1-2, 2008, pp. 39-46. doi:10.1016/j.molliq.2008.03.002
- M. J. Wójcik, “Theory of the Infrared Spectra of the Hydrogen Bond in Molecular Crystals,” International Journal of Quantum Chemistry, Vol. 10, No. 5, 1976, pp. 747- 760. doi:10.1002/qua.560100506
- A. Witkowski and M. J. Wójcik, “Infrared Spectra of Hydrogen Bond a General Theoretical Model,” Chemical Physics, Vol. 1, No. 1, 1973, pp. 9-16. doi:10.1016/0301-0104(73)87017-X
- M. J. Wójcik, “Fermi Resonance in Dimers: A Model Study,” Molecular Physics, Vol. 36, No. 6, 1978, pp. 1757- 1767. doi:10.1080/00268977800102741
- P. M. Morse, “Diatomic Molecules According to the Wave Mechanics. II. Vibrational Levels,” Physics Reviews, Vol. 34, No. 1, 1929, pp. 57-64. doi:10.1103/PhysRev.34.57
- N. Issaoui, N. Rekik, B. Oujia and M. J. Wójcik, “Theoretical Infrared Line Shapes of H-Bonds within the Strong Anharmonic Coupling Theory and Fermi Resonances Effects,” International Journal of Quantum Chemistry, Vol. 110, No. 14, 2010, pp. 2583-2602. doi:10.1002/qua.22395
- Y. Maréchal and A. Witkowski, “Infrared Spectra of HBonded Systems,” Journal of Chemical Physics, Vol. 48, No. 8, 1968, pp. 3697-3704. doi:10.1063/1.1669673
- N. Rösh and M. Ratner, “Model for the Effects of a Condensed Phase on the Infrared Spectra of HydrogenBonded Systems,” Journal of Chemical Physics, Vol. 61, No. 8, 1974, pp. 3344-3351. doi:10.1063/1.1682497
- B. Boulil, O. H. Rousseau and P. Blaise, “Infrared Spectra of Hydrogen Bonded Species in Solution,” Chemical Physics, Vol. 126, No. 2-3, 1988, pp. 263-290. doi:10.1016/0301-0104(88)85038-9
- H. T. Flakus, A. Tyl and A. Maślankiewicz, “ElectronInduced Phase Transition in Hydrogen-Bonded SolidState 2-Pyridone,” Journal of Physical Chemistry A, Vol. 115, No. 6, 2011, pp. 1027-1039. doi:10.1021/jp108717v
- P. Blaise, M. J. Wójcik and O. H. Rousseau, “Theoretical Interpretation of the Line shape of the Gaseous Acetic Acid Dimer,” Journal of Chemical Physics, Vol. 122, No. 6, 2005, pp. 64306-64317. doi:10.1063/1.1847491
- M. J. Frisch, “Gaussian 03, Revision B.03,” Gaussian, Inc., Pittsburgh, 2003.
- R. Kubo, “Statistical-Mechanical Theory of Irreversible Processes. I. General Theory and Simple Applications to Magnetic and Conduction Problems,” Journal of Physics Society Japan, Vol. 12, No. 6, 1957, pp. 570-586. doi:10.1143/JPSJ.12.570
- R. Kubo, “Lectures in Theoretical Physics I,” In: W. E. Brittin and L. G. DUnham, Eds., Interscience, 1st Annual Summer Institute for Theoretical Physics, Colorado, 16 June-22 August 1958, pp. 120-203.
- M. El-A. Benmalti, D. Chamma, P. Blaise and O. H. Rousseau, “Theoretical Interpretation of the Infrared Lineshape of Gaseous Propynoic and Acrylic Acid Dimers,” Journal of Molecular Structure, Vol. 785, No. 1-3, 2006, pp. 27-31. doi:10.1016/j.molstruc.2005.09.036
- K. Heyne, N. Huse, J. Dreyer, E. T. J. Nibbering, T. Elsaesser and S. Mukamel, “Coherent Low-Frequency Motions of Hydrogen Bonded Acetic Acid Dimers in the Liquid Phase,” Journal of Chemical Physics, Vol. 121, No. 2, 2004, pp. 902-913. doi:10.1063/1.1762873
- Z. Mielke and L. Sobczyk, “Vibrational Isotope Effects in Hydrogen Bonds,” In: A. Kohen and H.-H. Limbach, Eds., Isotope Effects in Chemistry and Biology, CRC Press, Taylor and Francis Group, West Palm Beach, London, 2006, pp. 281-304.
- N. Issaoui, N. Rekik, B. Oujia and M. J. Wójcik, “Anharmonic Effects on Theoretical IR Line Shapes of Medium Strong H(D)Bonds,” International Journal of Quantum Chemistry, Vol. 109, No. 3, 2009, pp. 483-499. doi:10.1002/qua.21839
- A. Novak, “Hydrogen Bonding in Solids Correlation of Spectroscopic and Crystallographic Data,” Structure and bonding, Vol. 18, 1974, pp. 177-216.
- I. Olovsson and P. G. Jönsson, “Hydrogen Bond,” In: P. Scheuster, G. Zundel and C. Sandorfy, Eds., Elsevier, Amsterdam, 1976, pp. 393-456.
- N. Rekik, B. Oujia and M. J. Wójcik, “Theoretical Infrared Spectral Density of H-bonds in Liquid and Gas Phases: Anharmonicities and Dampings Effect,” Chemical Physics, Vol. 352, No. 1-3, 2008, pp. 65-76. doi:10.1016/j.chemphys.2008.05.009
- S. E. Odinokov and A. V. Iogansen, “Torsional γ(OH) Vibrations, Fermi Resonance [2γ(OH) ν(OH)] and Isotopic Effects in i.r. Spectra of H-Complexes of Carboxylic Acids with Strong Bases,” Spectrochimica Acta Part A: Molecular Spectroscopy, Vol. 28, No. 12, 1972, pp. 2343-2350. doi:10.1016/0584-8539(72)80214-9
- P. Blaise, M. J. Wójcik and O. H. Rousseau, “Theoretical Interpretation of the Line Shape of the Gaseous Acetic Acid Cyclic Dimer,” Journal of Chemical Physics, Vol. 122, No. 6, 2005, p. 064306. doi:10.1063/1.1847491
- M. El-A. Benmalti, P. Blaise, H. T. Flakus and O. H. Rousseau, “Theoretical Interpretation of the Infrared Line Shape of Liquid and Gaseous Acetic Acid,” Chemical Physics, Vol. 320, No. 2-3, 2006, pp. 267-274. doi:10.1016/j.chemphys.2005.07.032
- J. L. Leviel and Y. Marechal, “Infrared Spectra of HBonded Systems: Anharmonicity oh the H-Bond Vibrations in Cyclic Dimers,” The Journal of Chemical Physics, Vol. 54 No. 3, 1971, pp. 1104-1107. doi:10.1063/1.1674943
- D. Hadzi, “Infrared Spectra of Strongly Hydrogen-Bonded Systems,” Pure and Applied Chemistry, Vol. 11, No. 3, 1965, pp. 435-454. doi:10.1351/pac196511030435
- J. Emsley, “Very Strong Hydrogen Bonding,” Chemical Society Reviews, Vol. 9, No. 1, 1980, pp. 91-124. doi:10.1039/cs9800900091
- W. H. Louisell and L. R. Walker, “Density-Operator Theory of Harmonic Oscillator Relaxation,” Physical Review, Vol. 137, No. 1B, 1965, pp. B204-B211. doi:10.1103/PhysRev.137.B204
- R. Feynmann and F. Vernon, “The Theory of a General Quantum System Interacting with a Linear Dissipative System,” Annals of Physics, Vol. 24, 1963, pp.118-173. doi:10.1016/0003-4916(63)90068-X
- W. H. Louisell, “Quantum Statistical Properties of Radiations,” Wiley, New York, 1973.
- K. Belhayara, D. Chamma, A. Velcescu and O. H. Rousseau, “On the Similarity of the IR Lineshapes of Weak H-Bonds in the Gas and Liquid Phases: Quantum Combined Effects of Strong an Harmonic Coupling, Multiple Fermi Resonances, Weak Dampings and Rotational Structure,” Journal of Molecular Structure, Vol. 833, No. 1-3, 2007, pp. 65-73. doi:10.1016/j.molstruc.2006.09.001
- L. C. Thomas, R. A. Chittenden and H. E. Hartley, “Hydrogen Bonding in Alkyl Phosphonic and Alkyl Phosphonothionic Acids,” Nature, Vol. 192, 1961, pp. 1283-1284. doi:10.1038/1921283a0
- K. G. Tokhadze, G. S. Denisov, M. Wierzejewska and M. Drozd, “First Example of the ABC ν(OH) Absorption Structure for Both Gaseous and Crystalline Phase: Infrared Studies of Dimethylphosphinic Acid,” Journal of Molecular Structure, Vol. 404, No. 1-2, 1997, pp. 55-62. doi:10.1016/S0022-2860(96)09360-X
- R. E Asfin, G. S. Denisov and K. G. Tokhadze, “The Infrared Spectra and Enthalpies of Strongly Bound Dimers of Phosphinic Acids in the Gas Phase. (CH
_{2}Cl)_{2}POOH and (C_{6}H_{5})_{2}POOH,” Journal of Molecular Structure, Vol. 608, No. 2-3, 2002, pp. 161-168. doi:10.1016/S0022-2860(01)00925-5 - H. T. Flakus, “On the Vibrational Transition Selection Rules for the Centrosymmetric Hydrogen-Bonded Dimeric Systems,” Journal of Molecular Structure: Theochem, Vol. 187, 1989, pp. 35-53.
- H. T. Flakus and A. Miros, “Infrared Spectra of the Hydrogen Bonded Glutaric Acid Crystals: Polarization and Temperature Effects,” Journal of Molecular Structure, Vol. 484, No. 1-3, 1999, pp. 103-115.
- H. T. Flakus and M. Jabonska, “Study of Hydrogen Bond Polarized IR Spectra of Cinnamic Acid Crystals,” Journal of Molecular Structure, Vol. 707, No. 1-3, 2004, pp. 97- 108.
- D. Chamma and O. H. Rousseau, “Infrared Spectra of weak H-Bonds: Beyond an Adiabatic Description of Fermi Resonances,” Chemical Physics, Vol. 248, No. 1, 1999, pp. 91-104.
- N. Issaoui, H. Ghalla and B. Oujia, “Theoretical Study of Hydrogen and Deuterium Bond in Glutaric Acid Crystal Dimer,” International Journal of Quantum Chemistry, Vol. 112, No. 4, 2012, pp. 1006-1015.

NOTES

^{*}Corresponding author.