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 N2 diluted in liquid Ar and pure liquid N2, both at normal conditions. We have analyzed the solvent shifts influence of the linear and quadratic (in the vibrational coordinate) oscillator-bath 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 [17]. 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 non-linear () 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 N2 diluted in liquid Ar and pure N2 liquid, both at normal conditions, founding that the non-linear 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 N2 () [6, 7]; however, because the HCl anharmonicity constant is clearly much larger than the N2 one, we expect a major contribution of non-linear 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 non-linear 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 fifth-order Dunham anharmonic oscillator [6, 7]. Fourth- and fifth-order 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 diatomic-solvent interaction potential depending on the vibration. Like in the previous work [8] we use for such part of the HCl-rare 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 site-site HCl-Ar potential of the previous work [8] (obtained by fitting the site-site potential to the functional form proposed by Holmgren et al. [10]), while for HCl diluted in Kr and Xe we have employed the site-site HCl-rare atom potentials obtained by fitting the site-site 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 non-linear 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 oscillator-bath 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 HCl-rare atom interaction depending of the vibration, the fitting of the site-site 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 S-B interaction is zero. This condition always can be achieved by redefining the S-B 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 S-B 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 S-B 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 p-solvent 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 S-B interaction can be written as where with . Then, the bath average of the isotropic part of the S-B 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 second-order 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 non-linear solvent shift, also called anharmonic shift, deduced by Alessi et al. [7]. Using the relations ((13), (14)) we can write the linear and non-linear 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 HCl-solvent 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 site-site models which take into account the presence of the vibrational coordinate in the site-site distances, also it is assumed that the Lennard-Jones parameters and present a linear dependence with the vibrational coordinate [57]. 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 Lennard-Jones parameters of the diatomic-solvent 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 Leap-Frog-Verlet 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 106 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 Lennard-Jones potentials, in our case with the values of the parameters , for Ar-Ar, , for Kr-Kr, and , for Xe-Xe [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 site-site potential. Because actually there are more accurate potential models available for the HCl-Ar (Kr,Xe) interaction [10, 11]; the site-site potentials were fitted to these more precise potential models. For the HCl-Ar we adopt the result obtained by Medina et al. [16] which fitted the site-site model to the Holmgren et al. potential [10]; the values obtained for the site-site 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 HCl-Ar. In the HCl-Kr case we have fitted the Hutson and Howard potential [11] considering the interval [3.5,10.5]. For the HCl-Xe case we have fitted the Hutson potential [11] by considering the [3.7,11] interval. In Table 2 are shown the values of the site-site parameters obtained for HCl-Kr and HCl-Xe 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 liquid-vapor 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 vdW-711 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 liquid-vapor 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 liquid-vapor 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 abundance-weighted 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 HCl-Ar, HCl-Kr, and HCl-Xe, 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 HCl-Ar, and a 19% for the HCl-Kr. The non-linear 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 HCl-Ar where the theoretical experimental differences are less than 10% for all the temperatures. For the HCl-Kr these differences are less than 10% unless for its two maxima temperatures, while for the HCl-Xe 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 diatomic-atomic 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 non-linear 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 non-linear (20) contributions, and as can be appreciated, the Buckingham term is the principal contribution to the fundamental solvent shifts . However, the anharmonic non-linear terms for the HCl in rare gases are greater than the corresponding values for pure N2 liquid and N2 diluted in liquid Ar at normal conditions, where the non-linear 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 N2 molecule. So for the first ten overtones of N2 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 non-linear 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 non-linear 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 N2 liquid and N2 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 non-linear 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 non-linear 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 non-linear anharmonic term makes a modest contribution to the fundamental solvent shift; however, this contribution is more important than that observed in pure N2 liquid and N2 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.