Theoretical Simulation of the Infrared Absorption Spectrum of the Strong Hydrogen and Deuterium Bond in 2-Pyridone Dimer

This work presents a theoretical simulation of the infrared spectra of strong hydrogen bond in alpha-phase 2-pyridone dimers, as well as in their deuterium derivatives at room temperature. The theory takes into account an adiabatic anharmonic coupling between the high-frequency N-H(D) stretching and the low-frequency intermolecular N...O stretching modes by considering that the effective angular frequency of the fast mode N-H(D) is assumed to be strongly dependent on the slow mode stretching coordinate N...O, the intrinsic anharmonicity of the low-frequency N...O mode through a Morse potential, Davydov coupling triggered by resonance exchange between the excited states of the fast modes of the two hydrogen bonds involved in the cyclic dimer, multiple Fermi resonances between the N-H(D) stretching and the overtone of the N-H(D) bending vibrations and the direct and indirect damping of the fast stretching modes of the hydrogen bonds and of the bending modes. The IR spectral density is computed within the linear response theory by Fourier transform of the autocorrelation function of the transition dipole moment operator of the N-H(D) bond. The theoretical line shapes of the υN-H(D) band of alpha-phase 2-pyridone dimers are compared to the experimental ones. The effect of deuteration is successfully reproduced.


Introduction
The hydrogen bonding is responsible for the appearance of spectacular changes in IR spectra of associated molecules.This remark particularly concerns the υ X-H proton stretching vibration bands, which are the attribute of the X-H atomic groups in the X-H…Y bridges.The main effects depend on a considerable low-frequency υ X-H band shift, on the band integral intensity growth by even two orders of magnitude and on a noticeable increase of the band-width [1][2][3].Numerous theoretical studies on hydrogen-bonded systems such as carboxylic acid have to explain the changes in the infrared spectrum induced by the formation of the H-bond Bridge.
In this work we present a theoretical approach for vibrational couplings in moderately strong hydrogenbonded systems and use it for simulating experimental infrared spectra for strong hydrogen-bonded of alpha solid-state phase 2-pyridone crystal dimers.It is shown alpha phase of 2-pyridone forms approximately a centrosymmetric dimer.In fact, that from the results given in the original crystallographic work, it results that the dimers are not ideally centrosymmetric since one hydrogen bridge in a dimer is longer than the other.Without making significant errors, in the discussion of spectroscopic effects, one may assume that the dimers are approximately centrosymmetric.The alpha-phase crystals belong to the monoclinic system with Z = 8 [4].However, we think that for the IR spectroscopy reasons you may safely consider these dimers as centrosymmetric ones.Several studies show that 2-pyridone form centrosymmetric dimers [5,6].
The studies have proved that in the case of cyclic dimers in this phase of 2-pyridone a single hydrogen bond is considerably stronger than a single hydrogen bond in a chain of molecules in the beta-phase.This fact suggests that the N + -H…O -bonds exist only in the alpha-phase and that these bonds are obviously stronger than the N-H…O hydrogen bonds in the beta-phase, since ionic hydrogen bonds are generally known to be stronger than the neutral ones [7,8].2-pyridone is an important precursor of antibiotics, which are used as inhibitors of the DNA gyrase.Thus 2-pyridone [9] and its derivatives [10,11] are widely investigated and exploited.Wuest and coworkers have proved that dipyridone can be formed by linking the two 2-pyridone using the functional groups like the acetylene or amine [12,13].Asymmetric dipyridones are ready to form a dimer, whereas symmetric dipyridones are easy to self-assembly to form the linear and planar polymers in aprotic solvent [12].The formation of the polymers via hydrogen bonding is an important approach of preparing liquid crystal and other functional polymers.There have been only a few reports on the hydrogen-bonding complexes formed by dipyridones.In addition, 2-pyridone dimer is an important reference system because it represents the rare case of a closed complex with two linear H-bonds.Recently, Wójcik presented a theoretical simulation of the bandshape and fine structure of υ N-H stretching band for 2-pyridone-H and 2-pyridone-D [14], taking into account an adiabatic coupling between high-frequency N-H(D) stretching and low-frequency intermolecular N…O stretching modes, Davydov interactions and linear and quadratic distortions of potential energies for the low-frequency vibrations in the excited state of the N-H(D) stretching vibrations.The slow vibration modes were assumed as harmonic potentials.In other work, Wójcik [15] has presented a theoretical model for vibrational couplings in weak and moderately strong hydrogen-bonded systems and use it for modeling experimental infrared spectra for hydrogenbonded crystals and hexagonal ice.The model is based on vibronic-type couplings between high and low frequency modes in hydrogen bridges, Davydov interactions [16] and Fermi resonance [17,18].It allows calculation of energy and intensity distributions in the infrared spectra of hydrogen-bonded systems.The present theory is based on strong anharmonic coupling between the highfrequency hydrogen stretching vibration υ N-H and lowfrequency phonons, Davydov interactions and multiple Fermi resonances interaction between a fundamental vibration of a υ N-H and one overtone or a combination tone of intermolecular vibration.Besides, this theory incorporated the intrinsic anharmonicity of the slow frequency mode through a Morse potential [19,20] whereas the fast mode was considered as harmonic.Note that, the Morse potential is undoubtedly more realistic than the harmonic one to describe this slow frequency mode.The adiabatic approximation [21] have been performed for each separate part of the dimer together with a strong non-adiabatic correction via the resonant exchange between the excited states of the two fast mode moieties.Both quantum direct (relaxation of the high-frequency modes) and indirect (relaxation of the H-bond bridges) dampings of the systems [22,23] were taken into account.This theory allows calculation of intensity distributions in the infrared spectra of hydrogen-bonded systems.The main pur-pose is to reproduce the experimental υ N-H IR line shapes of hydrogen and deuterium bond in alpha-phase 2-pyridone crystal dimer at room temperature.We shall use infrared spectra of 2-pyridone in the alpha solid-state phase, measured by Flakus at the room temperature (Figures 1 and 2 from [24]).
The numerical results show that this theoretical approach allows fitting the experimental υ N-H infrared line shapes of cyclic alpha-phase 2-pyridone crystal dimer and its deuterium derivative by using a minimum number of parameters.With such tools, experimentalist should be able to compare experimental and theoretical data in an easy-to-use way.

Experimental Spectra
It is of interested to note that, 2-pyridone compound which has been used in this investigation was a commercial substance (Sigma-Aldrich).The single crystals of the alpha phase formed rectangular plates.For the purpose of the experiment they were exposed by using on a tin diaphragm with a 1.5 mm hole diameter.It was used without further purification.Crystals of 2-pyridone suitable for spectral studies were obtained by crystallization from  melt between two closely spaced CaF2 windows.The solid-state spectra were measured by a transmission method at room temperature with the help of the FT-IR Nicolet Magna 560 spectrometer using a non polarized beam.The spectra were measured at 2 cm -1 resolution.Measurements of spectra were completed in a similar way for the crystals of deuterium derivative of 2-pyridone, which was synthesized by evaporation of the solution in D 2 O under reduced pressure.

Theory
The present theory was developed in the framework of the adiabatic approximation [25].The adiabatic approximation of the N-H stretching vibrations leads to the description of each moiety by effective Hamiltonians of the H-bond Bridge: for a single H-bond bridge, this effective Hamiltonian is either that of a harmonic oscillator, if the fast mode is in its ground state, or that of a driven harmonic oscillator if the fast mode is excited.When one of the two identical fast modes is excited, then, because of the symmetry of the cyclic dimer and because of coupling V 0 between the two degenerate fast-mode excited states, an interaction occurs (Davydov coupling) leading to an exchange between the two identical excited parts of the dimer, as considered by Maréchal and Witkowski in their pioneering work [21].Of course, this interaction between degenerate excited states is of non-adiabatic nature, although the adiabatic approximation was performed to separate the high-and low-frequency motions.Figure 1 present the geometry of the 2-pyridone dimer optimized at the HF/6-311 ++ G(d,p) level [26], which is formed by two hydrogen bonds (lengths of N6-H19… O20 and N21-H15…O22 are 2.89559 Å and 2.89558 Å, respectively).
Let q 1 and q 2 , be the coordinates of the high-frequency N-H stretching vibrations, in the first and second hydrogen bond, and Q 1 and Q 2 coordinates of the two low-frequency intermolecular N…O stretching modes.The two moieties of the system are exchanged by the symmetry C2 parity operator.
For the theory dealing with this system, the basic physical parameters are 1) vibration angular frequency ω˚ of the two degenerates fast modes moieties when the Hbond bridge is at the equilibrium, 2) vibration angular frequency Ω of the two degenerates H-bond bridge moieties, 3) dimensionless anharmonic coupling parameter α˚ between the high frequency mode of one moiety and the H-bond bridge coordinate of the same moiety, 4) Davydov coupling V parameter between the degenerate first excited state of the two moieties high frequency mode, 5) direct and indirect damping parameters γ and γ, 6) coupling parameters fi involved in the Fermi resonances coupling between the first harmonics of some bending modes of angular frequency and the first excited state of the g symmetrized high frequency mode, 7) relaxations parameter of the first harmonic bending modes, and 8) absolute temperature T.

Full Hamiltonian of the System
The evaluation of the spectral density of the hydrogen bond system requires the knowledge of the full Hamiltonian of the hydrogen bond system.For this purpose, it is important to describe the basic various vibration modes in the dimer.In our present work, for each part of the dimer, we have described the slow-frequency mode by a Morse potential, which can be written as follows: D e is the dissociation energy of the Hydrogen Bond Bridge and e β is given by: The fast frequency-mode is considered to be harmonic.It is important to note that in the majority of recent works, the slow-frequency mode was assumed as harmonic [25].
Within the strong anharmonic coupling theory and the anharmonic approximation for the H-bond bridge, the corresponding Hamiltonians of the slow and high-frequency modes of the two moieties of the dimer are, respectively given, using dimensionless operators, by: In these equations, P i are the dimensionless conjugate momenta of the H-bond bridges coordinates Q i of the two moieties, whereas p i and q i are the dimensionless coordinates and the conjugate momenta of the two degenerate high frequency modes of the two moieties.Ω is the angular frequency of the H-bond Bridge, whereas ω(Q i ) is that of the high frequency mode which is supposed to depend on the coordinate of the H bond bridge.
Using a Taylor development of the Morse potential, the Hamiltonian of the slow frequency modes given in Equation ( 3) can be rewritten as the sum of the Hamiltonian of a quantum harmonic oscillator and an anharmonic potential V: where [H Slow ] i and V are respectively given by: Expansion to first-order of the angular frequency of the fast mode with respect to the coordinate of the H-bond bridge leads to write: where is the angular frequency of the two degenerate fast modes when the corresponding H-bond Bridge coordinates are at equilibrium, whereas is a dimensionless parameter which will appear to be an anharmonic coupling parameter.

ω α
In presence of damping, the thermal bath may be figured, by an infinite set of harmonic oscillators, and its coupling with the H-bond bridge are described by terms which are linear in the position coordinates of the bridge and of the bath oscillators: Here, are the dimensionless position coordinate operators of the oscillators of the bath, are the corresponding conjugate moments, obeying the usual quantum commutation rules, are the corresponding angular frequencies and r r q  r p  r ω g are the coupling between the H-bond bridges and the oscillators of the bath.Within the adiabatic approximation, the Hamiltonian of each moiety of the dimer takes the form of sum of effective Hamiltonians which are depending on the degree of excitation of the fast mode according to: 1, 2; 0,1 An excitation, of the fast mode of one moiety of the dimer is resonant with the excitation of the other moiety.Thus, a strong non-adiabatic correction [21] i.e. a Davydov coupling V 12 , occurs between the two resonant states after excitation of one of the two fast modes, so that the full Hamiltonian of the two moieties is given by the equations: with,

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 : 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]: 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: In these equation, <n> and β are respectively the thermal average of the occupation number of the quantum harmonic oscillator describing the H-bond bridge and the effective dimensionless anharmonic coupling parameter related to α˚.These terms are given by: 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]: where, , n C    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 1 2 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: That may be written formally: with respectively: In view of the above equations, the components appearing here are given by the following equations [29]: ,   , 2 , Recall that n 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.

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]: Here, {ω μ } g are the eigenvalues appearing in Equation ( 29) whereas, a {μ,0,m}g are the expansion coefficients defined by Equation (30).
The g states involved in the above expansions are defined by: Here, the kets and   g l 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   1 j of the jth fast mode and second damped excited state   2 j 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 : with 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: with and In these two last equations, a + and a are the boson operators obeying [a + , a] = 1 and i 2 = -1.

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.Table 1.Parameters used for fitting the experimental spectra of alpha-phase 2-pyridone at room temperature in the presence of three Fermi resonances.
Species ω˚ (cm -1 ) Ω (cm -1 ) α˚V˚ γ˚ (cm -1 ) γ (cm  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 2 .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 2 for weak bonds and decreases with the bond strength reaching minimum near 0.9

 
for low-barrier bonds.In our cases, to obtain good agreement between the theoretical and experimental one, we have used: 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: ; 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 deal-ing 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 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.

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.

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

Figure 3 .
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.

Figure 4 .Figure 5 .
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 6 .
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.