Research Article  Open Access
A. Padilla, J. Pérez, "A Simulation Study of the Fundamental Vibrational Shifts of HCl Diluted in Ar, Kr, and Xe: Anharmonic Corrections Effects", International Journal of Spectroscopy, vol. 2012, Article ID 851923, 7 pages, 2012. https://doi.org/10.1155/2012/851923
A Simulation Study of the Fundamental Vibrational Shifts of HCl Diluted in Ar, Kr, and Xe: Anharmonic Corrections Effects
Abstract
We have calculated the vibrational solvent shifts of the fundamental bands of HCl diluted in Ar, Kr, and Xe solutions at different thermodynamic conditions by means of the molecular dynamics technique and a model for the isotropic part of the interaction depending on the vibration. The theoretical vibrational shifts, which were compared with the available experimental data, have been determined by considering both, the usual linear Buckingham terms and the nonlinear anharmonic corrections, and the latter omitted in a previous work for the HCl in Ar and Kr. We have found that the Buckingham contributions dominate the solvent shifts of the fundamental bands of HCl in Ar, Kr, and Xe, although the anharmonic shifts’ present significant greater values than those obtained previously for N_{2} diluted in liquid Ar and pure liquid N_{2}, both at normal conditions. We have analyzed the solvent shifts influence of the linear and quadratic (in the vibrational coordinate) oscillatorbath interaction terms and also the Dunham intramolecular potential effects on the anharmonic contributions.
1. Introduction
The study of the fundamental and overtones infrared and Raman vibrational solvent shifts gives relevant information about the molecular interaction depending of the intramolecular vibration [1–7]. In an earlier work, according to Buckingham, it was deduced that the diatomic solvent shift associated to the vibrational transition of a cubic anharmonic oscillator presents a linear dependence with the vibrational quantum number [1]. This result, which is valid in principle for the fundamental band or lower overtones of several diatomics, would break down for the fundamental bands of some molecules and high overtones of much ones. Recently, Alessi et al. [6, 7] has deduced the nonlinear () correction (anharmonic correction) to the linear Buckingham solvent shift for a diatomic molecule perturbed by an atomic or molecular solvent. Alessi et al. [6, 7] analyzed the solvent shift of N_{2} diluted in liquid Ar and pure N_{2} liquid, both at normal conditions, founding that the nonlinear solvent correction, also called anharmonic shift [6, 7], represents only the 1% of the fundamental shift of the isotropic Raman band. This result is due in part to the small anharmonicity parameter of N_{2} () [6, 7]; however, because the HCl anharmonicity constant is clearly much larger than the N_{2} one, we expect a major contribution of nonlinear solvent shift for the fundamental band of HCl diluted in Ar, Kr, and Xe.
We have made a previous study [8] of the vibrational solvent shifts of HCl diluted in Ar and Kr considering the usual linear Buckingham approach. In the present work we will study the solvent shifts of HCl also in Ar and Kr and additionally in liquid Xe by considering the nonlinear corrections to the solvent shift deduced by Alessi et al. [6, 7]. Like in the previous work [8] both the diatomic vibration and rotation will be treated with a quantum approach (small system S), while the translational degree of freedom, both of the diatomic and the solvent molecules, will be treated classically (the bath B). Unlike the previous study [8] where the vibration was modeled by means of a cubic anharmonic oscillator, in this work the vibration will be characterized by using a fifthorder Dunham anharmonic oscillator [6, 7]. Fourth and fifthorder terms of the intramolecular potential must be considered because they have relevant effects on the anharmonic solvent shift. On the other hand, classical evolution of the diatomic and solvent translational coordinates (the bath B) will be represented by using molecular dynamic simulation techniques for fixed values of temperature and density.
A key element in the solvent shift calculation is to know the functional form of the isotropic part of diatomicsolvent interaction potential depending on the vibration. Like in the previous work [8] we use for such part of the HClrare atom potential the model proposed by Marteau el al. [9], which depends only on a fitting parameter, that will be determined by comparing the theoretical and the experimental solvent shifts, at the different thermodynamic states in which the HCl absorption bands were measured.
For the solvent shifts of HCl diluted in dense Ar we have used the same sitesite HClAr potential of the previous work [8] (obtained by fitting the sitesite potential to the functional form proposed by Holmgren et al. [10]), while for HCl diluted in Kr and Xe we have employed the sitesite HClrare atom potentials obtained by fitting the sitesite functional form to the accurate potentials proposed by Hutson and Howard [11].
In the present study we have analyzed the relative contributions of the linear Buckingham shift and the nonlinear anharmonic term to the vibrational solvent shift of HCl in Ar, Kr, and Xe, for the different thermodynamic states considered in the experimental study [12]. Additionally, we have analyzed the contributions to the solvent shift of the linear and quadratic (in the diatomic vibrational coordinate) parts of the oscillatorbath interaction. The effects on the anharmonic solvent shift of the third, fourth, and fifth Dunham potential terms also are discussed.
The organization of the paper is as follows. In Section 2 we develop the basic theoretical aspects related with the vibrational solvent shifts. In Section 3 we introduce the model for the isotropic part of the HClrare atom interaction depending of the vibration, the fitting of the sitesite potentials, and the details of the molecular dynamics simulations. In Section 4 we compare the theoretical and the experimental solvent shift and discuss different aspects related with this magnitude. Finally, a brief summary and some conclusions are collected in Section 5.
2. Theoretical Background
In order to deduce the basic expressions of the solvent shift we make reference to a spectral theory previously developed [13], although the results obtained are more general than the proper theory. This theory [13] considers the evolution of a small system S, formed by the diatomic rotation and vibration, which interacts with a thermal bath B, constituted by the translational degree of freedom, both of the diatomic and solvent molecules. In this theory it is assumed that the bath average of the SB interaction is zero. This condition always can be achieved by redefining the SB interaction and the system S Hamiltonian as , and respectively. It is precisely the new Hamiltonian who defines the frequencies of the vibrorotational spectral lines, introducing the term the frequency shifts which are the object of our study.
The SB interaction Hamiltonian can be divided into the isotropic and the anisotropic parts: Because the bath is macroscopically isotropic, in lowest order approximation, when the anisotropic intermolecular forces are neglected in the bath average, the mean value of the anisotropic part of the interaction is zero [14]; therefore, in this approach, the frequency shift produced by the solvent is basically due to the isotropic part of the interaction and hence to the diatomic vibration. Although Alessi et al. [6, 7] included the anisotropic force effects on bath averages in their study of pure and diluted in liquid , we have neglected these effects for in rare gases in the present study; they will be considered in a future work.
The isotropic part of the SB interaction can be written as the sum of the binary isotropic potentials depending of the vibration : where is the vibrational coordinate and is the norm of the vector joining the center of mass of the diatomic and the psolvent molecule. The binary potential will be approximated by means of the first two terms of its power expansion in the vibrational coordinate: so the isotropic SB interaction can be written as where with . Then, the bath average of the isotropic part of the SB potential (4) is given by where is the bath average of the coupling coefficient (5). In [7] also the cubic and quartic terms of were considered, but it was found that their contributions to the solvent shift were small.
The vibrorotational Hamiltonian is modeled as a rigid rotor coupled to an anharmonic oscillator : The anharmonic oscillator is represented initially by the usual Morse model: where is the vibrational momentum, is the reduced mass of the diatomic, and and are the parameters of the Morse potential. The energy spectrum of the Morse Hamiltonian is [7] where the harmonic frequency is and the anharmonicity parameter is The expression (9) is usually used to fit spectroscopy data, by means of which obtain the values of the vibrational harmonic frequency and the anharmonicity .
When the Morse potential (8) is expanded in powers of the vibrational coordinate up to the fifth order we obtain the Dunham model of the anharmonic oscillator [7]: where In order to relating more easily the perturbation coefficient with the vibrational spectroscopy constants and , we rearrange the expressions ((10), (11)) as [7]
By making secondorder quantum mechanical perturbation theory calculations with the Hamiltonian ((6), (12)) we obtain the following expression for the solvent shift associated to the vibrational transition [7]: where is the linear solvent shift predicted by Buckingham’s theory [1], while is the nonlinear solvent shift, also called anharmonic shift, deduced by Alessi et al. [7]. Using the relations ((13), (14)) we can write the linear and nonlinear solvent shifts ((16), (17)) in terms of the spectroscopic constants and , that is, defining the coefficients: where is the equilibrium intermolecular distance of the diatomic; the solvent shifts can be expressed as where is the diatomic rotational constant. The first summand in (20) corresponds to the two terms of in (17), the second summand corresponds to the term in (17), and the last summand corresponds to the two terms of in (17). The first term of (20) is the result of the large cancelation effect of the two terms of (17) [7], being the unique summand that presents opposite sign to .
3. Interaction Potentials and Simulations
The HClsolvent binary interaction potential depending on the vibration is a key element in the calculation of the vibrational solvent shift of the diatomic molecule. Not many models for this part of the isotropic interaction depending of the vibration are available in the literature. For this kind of potential it is usual to consider sitesite models which take into account the presence of the vibrational coordinate in the sitesite distances, also it is assumed that the LennardJones parameters and present a linear dependence with the vibrational coordinate [5–7]. However, in this work, like in the previous one [8], we consider the model introduced by Marteau et al. [9] for the binary isotropic interaction potential depending on the interaction. By writing the vibrational components (3) as the binary isotropic components of Marteau et al. can be expressed: where and are the usual 6–12 LennardJones parameters of the diatomicsolvent interaction. and are given by where is the diatomic polarizability and , while and are given by where is a phenomenological parameter, called the repulsive factor, which will be obtained by comparing the theoretical solvent shift with the experimental data [13]. Unless , all the parameters involved in the isotropic components (22) are collected in Table 1.
For the Molecular dynamics simulations we consider a single HCl molecule and 250 solvent atoms collected in a cubic box with periodic boundary conditions. In each case the size of the box was chosen to give the experimental solvent density . The motion equations were solved by using a LeapFrogVerlet algorithm which includes a SHAKE procedure to fix the diatomic bond length to and a coupling to a thermal bath to keep the system temperature near to the experimental solvent value . The time step was chosen as , including each simulation run with a total of 10^{6} time steps. Previous to each simulation run, a long equilibration period was considered.
As it is usual, the Molecular dynamics simulations were carried out considering that the interaction between the solvent atoms can be described by means of the 6–12 LennardJones potentials, in our case with the values of the parameters , for ArAr, , for KrKr, and , for XeXe [16].
Also, as it is usual in the Molecular dynamics simulations, the binary interaction between the diatomic molecule and the solvent atom will be modeled by means of a sitesite potential. Because actually there are more accurate potential models available for the HClAr (Kr,Xe) interaction [10, 11]; the sitesite potentials were fitted to these more precise potential models. For the HClAr we adopt the result obtained by Medina et al. [16] which fitted the sitesite model to the Holmgren et al. potential [10]; the values obtained for the sitesite parameters are shown in Table 2. Logically, the results of the fitting depend of the interatomic interval considered for the potential, so as to have similar fitting conditions to the HClAr. In the HClKr case we have fitted the Hutson and Howard potential [11] considering the interval [3.5,10.5]. For the HClXe case we have fitted the Hutson potential [11] by considering the [3.7,11] interval. In Table 2 are shown the values of the sitesite parameters obtained for HClKr and HClXe with this fitting procedure.

In the Molecular dynamics simulations it is necessary to specify the values of the solvent density and temperature. However, our experimental information is referred to the pressure and temperature [12], so the density must be obtained by an indirect method. For the liquidvapor coexistence states the densities were determined by means of a cubic splines interpolation of the experimental data ( and ) reported by Bertsev and Zelikina [18]. While for the supercritical states the densities were calculated by using the vdW711 equation of state [19]. The experimental values of the pressure and temperature and the values obtained for the density are shown in Tables 3, 4, and 5 for Ar, Kr, and Xe, respectively. In all the cases the superscripts in the temperatures denote supercritical states, while the remaining ones are in the liquidvapor coexistence curve.



4. Vibrational Solvent Shifts
The experimental spectral profile of the fundamental absorption bands of HCl diluted in Ar, Kr, and Xe along the liquidvapor coexistence curve and in supercritical states has been reported by Pérez et al. [12]. From these absorption line shapes the band origin frequencies have been determined by fitting the vibrorotational spectral line frequencies of the rotational resonances [12]. From the band origin frequencies we have obtained the vibrational solvent shifts by using the isotope abundanceweighted HCl molecular constants [20] cm^{−1} and ; also we have considered the weighted value of the rotational constant cm^{−1} [20] which has been used in the theoretical calculations of the solvent shift ((19), (20)). The values of the experimental solvent shifts obtained in this way are collected in Tables 3, 4, and 5 for HClAr, HClKr, and HClXe, respectively.
We have calculated the theoretical solvent shifts of HCl by using ((15), (19), and (20)) calculating the bath averages of the coupling coefficients (5) by means of the molecular dynamics simulations and employing the values of the repulsive factor (24) which give the best fit of the theoretical to the experimental solvent shifts . The repulsive factors obtained in the fitting procedure are for the Ar, for the Kr, and for the Xe. These values are greater than those obtained with the Buckingham cubic model of the previous paper [8], a 18% major for the HClAr, and a 19% for the HClKr. The nonlinear Dunham model, ((15), (19), and (20)), produces red solvent shifts greater in magnitude than those of the linear Buckingham model of the previous paper [8], so it is necessary to increase the repulsive factor , in order to produce a decrease of the magnitude of the vibrational shifts and obtain theoretical solvent shifts near to the experimental values.
The theoretical solvent shifts of the HCl in Ar, Kr, and Xe solutions are collected in Tables 3, 4, and 5, respectively. Also, in Figure 1 we have plotted the theoretical and the experimental solvent shifts of HCl in Ar, Kr, and Xe for different values of the temperature. As it can be appreciated in general terms there is a good agreement between the theoretical and the experimental solvent shifts; the better accord is found for HClAr where the theoretical experimental differences are less than 10% for all the temperatures. For the HClKr these differences are less than 10% unless for its two maxima temperatures, while for the HClXe the differences are less than 10% unless for its three maxima temperatures. In all the cases the vibrational shift is a red shift which increases its magnitude following the sequence Ar, Kr, and Xe, that is, the sequence in which increasing the diatomicatomic interaction (see Table 1). On the other hand, we can observe that for the three solvents the vibrational shifts decrease their magnitude with increasing temperature and decreasing density, or in other words, with increasing mean interatomic distance.
In all the cases studied in this work the nonlinear anharmonic shift (20) contributes an 8% to the solvent shift (15), so for the HCl in rare gases the linear Buckingham approximation (19) dominates with a 92% of the solvent shift of the fundamental absorption bands. In Figure 2 we have plotted the theoretical solvent shifts together with the linear (19) and nonlinear (20) contributions, and as can be appreciated, the Buckingham term is the principal contribution to the fundamental solvent shifts . However, the anharmonic nonlinear terms for the HCl in rare gases are greater than the corresponding values for pure N_{2} liquid and N_{2} diluted in liquid Ar at normal conditions, where the nonlinear term only represents the 1% of the vibrational shifts [7]. These differences are due in part to the anharmonicity parameter [7] which is smaller for the N_{2} molecule. So for the first ten overtones of N_{2} the solvent shifts present a relatively linear behavior [7], such as predicting the Buckingham theory [1]. However, for the first ten overtones of the HCl in rare gases the solvent shifts present a nonlinear behavior such as is what shown in Figure 3, where we have plotted the vibrational shifts of the HCl in Ar together with the Buckingham approach . In this figure we can observe as the solvent shift diverge of the Buckingham linear theory conforms the quantum number is increased, so we can verify that for high overtones the anharmonic terms make very relevant contributions to the solvent shifts. In this way, we must take into account the nonlinear solvent shifts in the study of the vibrational spectra, particularly for the vibrational overtones, as is the case of the first overtone absorption bands of in fluid , for which there are experimental data at different thermodynamics conditions [21].
In contrast to pure N_{2} liquid and N_{2} in liquid Ar where the term is negative and the term is positive [7], for the HCl in rare gases both and terms ((19), (20)) are negative, so in the Buckingham approach the terms represent the 71% of the linear contribution , while the terms represent the 29% of the of the Buckingham solvent shift. For the nonlinear anharmonic shift the terms represent the 33% and the terms represent the 67%. In this last case the positive contributions of represent the 2% of , the contribution of represents the 36% of , and finally the terms of represent 67% of .
5. Summary and Conclusions
We have calculated the vibrational solvent shifts of the fundamental absorption bands of the HCl diluted in Ar, Kr, and Xe by means of the molecular dynamics simulation and the model proposed by Marteau et al. [9] for the isotropic part of the interaction depending of the vibration. The Marteau potential depends of a phenomenological parameter, the repulsive factor , which has been obtained by comparing the theoretical and the experimental solvent shifts. In this work, in contrast with the previous paper [8], we have taken into account the nonlinear anharmonic contribution to the solvent shift, obtaining a reasonable agreement between the theoretical and the experimental solvent shifts.
We have found that for HCl in rare gases the nonlinear anharmonic term makes a modest contribution to the fundamental solvent shift; however, this contribution is more important than that observed in pure N_{2} liquid and N_{2} diluted in Ar liquid [7]. On the other hand, for medium and high overtones the anharmonic shifts make an important contribution to the vibrational solvent shifts of HCl in rare gases.
Finally, we have quantified the effects on the vibrational solvent shifts of the linear and quadratic terms of the isotropic part of interaction depending on the vibration, and of the Dunham terms on the anharmonic solvent shifts.
Acknowledgments
The authors thank Pedro N. Leal for careful reading the paper. We also acknowledge the financial support from the Consejería de Educación Cultura y Deporte del Gobierno de Canarias.
References
 A. D. Buckingham, “Solvent effects in vibrational spectroscopy,” Transactions of the Faraday Society, vol. 56, pp. 753–760, 1960. View at: Publisher Site  Google Scholar
 D. W. Oxtoby, “Dephasing of molecular vibrations in liquids,” Advances in Chemical Physics, vol. 40, 1979. View at: Publisher Site  Google Scholar
 D. W. Oxtoby, “Vibrational relaxation in liquids,” Annual Review of Physical Chemistry, vol. 32, pp. 77–101, 1981. View at: Publisher Site  Google Scholar
 J. P. J. Michels, M. I. M. Scheerboom, and J. A. Schouten, “Computer simulations of the linewidth of the Raman Qbranch in fluid nitrogen,” Journal of Chemical Physics, vol. 103, no. 19, pp. 8338–8345, 1995. View at: Google Scholar
 K. F. Everitt and J. L. Skinner, “Isotropic Raman line shapes of N_{2} and O_{2} along their liquidgas coexistence lines,” Journal of Chemical Physics, vol. 115, no. 18, pp. 8531–8539, 2001. View at: Publisher Site  Google Scholar
 N. Alessi, I. S. Tolokh, S. Goldman, and C. G. Gray, “Simulation study of solvent shifts of vibrational overtone spectra,” Molecular Physics, vol. 102, no. 1920, pp. 2037–2047, 2004. View at: Publisher Site  Google Scholar
 N. Alessi, I. S. Tolokh, S. Goldman, and C. G. Gray, “Simulation study of N_{2} overtone solvent shifts using improved potentials,” Molecular Physics, vol. 103, no. 17, pp. 2381–2396, 2005. View at: Publisher Site  Google Scholar
 A. Padilla and J. Pérez, “The fundamental vibrational shift of HCl diluted in dense Ar and Kr,” Molecular Physics, vol. 106, no. 5, pp. 687–691, 2008. View at: Publisher Site  Google Scholar
 P. Marteau, C. Boulet, and D. Robert, “Finite duration of collisions and vibrational dephasing effects on the Ar broadened HF infrared line shapes: asymmetric profiles,” Journal of Chemical Physics, vol. 80, no. 8, pp. 3632–3639, 1984. View at: Google Scholar
 S. L. Holmgren, M. Waldman, and W. Klemperer, “Internal dynamics of van der Waals complexes. II. Determination of a potential energy surface for ArHCl,” Journal of Chemical Physics, vol. 69, no. 4, pp. 1661–1669, 1978. View at: Google Scholar
 J. M. Hutson and B. J. Howard, “Anisotropic intermolecular forces,” Molecular Physics, vol. 45, no. 4, pp. 769–790, 1982. View at: Publisher Site  Google Scholar
 J. Pérez, A. Padilla, W. A. Herrebout, B. J. Van Der Veken, A. C. Hernández, and M. O. Bulanin, “Experimental analysis and modified rotor description of the infrared fundamental band of HCl in Ar, Kr, and Xe solutions,” Journal of Chemical Physics, vol. 122, no. 19, Article ID 194507, 9 pages, 2005. View at: Publisher Site  Google Scholar
 A. Padilla, J. Pérez, and A. Calvo Hernández, “Vibrorotational Raman and infrared spectra of polar diatomic molecules in inert solutions. I. Spectral theory,” Journal of Chemical Physics, vol. 111, no. 24, pp. 11015–11025, 1999. View at: Google Scholar
 C. G. Gray and K. E. Gubbins, Theory of Molecular Fluids, vol. I, Claredon Press, Oxford, UK, 1984.
 L. Bonamy and P. N. M. Hoang, “Far infrared absorption of diatomic polar molecules in simple liquids and statistical properties of the interactions. I. Spectral theory,” Journal of Chemical Physics, vol. 67, pp. 4423–4430, 1977. View at: Google Scholar
 A. Medina, J. M. M. Roco, A. C. Hernández et al., “Vibrationrotation spectra of HCI in raregas liquid mixtures: molecular dynamics simulations of Qbranch absorption,” Journal of Chemical Physics, vol. 116, no. 12, pp. 5058–5065, 2002. View at: Publisher Site  Google Scholar
 M. O. Bulanin and K. G. Tokhadze, Optics and Spectroscopy, vol. 61, p. 198, 1986.
 V. V. Bertsev and G. Ya. Zelikina, “Electronic absorption spectra of some cryosystems,” in Molecular Cryospectroscopy, R. J. H. Clark and R. E. Hester, Eds., John Wiley & Sons, Chichester, UK, 1995. View at: Google Scholar
 D. P. Tassios, Applied Chemical Engineering Thermodynamics, Springer, Berlin, Germany, 1993.
 G. Guelachvili, P. Niay, and P. Bernage, “Infrared bands of HCl and DCl by Fourier transform spectroscopy. Dunham coefficients for HCl, DCl, and TCl,” Journal of Molecular Spectroscopy, vol. 85, no. 2, pp. 271–281, 1981. View at: Google Scholar
 M. O. Bulanin, K. Kerl, A. Padilla, J. Pérez, and A. Calvo Hernández, “First vibrational overtone bandshape of HCl in fluid SF_{6}: an experimental and theoretical study,” Physical Chemistry Chemical Physics, vol. 2, pp. 5375–5382, 2000. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 A. Padilla and J. Pérez. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.