Two-Dimension Study of Methanol Internal-Overall Rotation in Argon Matrix

In the approximation B3LYP/cc-pVTZ, the geometry of a methanol molecule surrounded by eight argon atoms has been optimized. By the independent rotation of the methyl and the hydroxyl groups at the fixed position of the C-O bond relative to the argon atoms there was obtained the two-dimensional grid of values of the internal-overall rotation energy. Despite the fact that, initially the energy was calculated for 65 points in the square 2 2 π π × , the presence of 2 / 3 π period for methyl group rotation has allowed to increase the number of points up to 195. The analytical approximation for internal rotation energy was found. Two dimensional Schrödinger equation for internal rotation – overall rotation of rotator with fixed axis was solved, energy levels, wave functions and transition probabilities were found. According to the results of these computations, degeneracy of the Е-type states is relieved with increase in splitting of the ground torsional state.


Introduction
It is known that low-temperature IR spectra of molecules in the matrix isolation characterized by narrow absorption bands due to the absence of intermolecular interactions and the rotational structure. However, matrix effects do exist [1][2][3] and they are often able to significantly complicate the interpretation of the spectra. These effects are mainly determined by stacking order of the matrix atoms around the trapped molecule. Therefore, modelling of the matrix structure and its influence on the vibrational spectra of the trapped molecules is an interesting and important task.
Theoretical study of the matrix influence on the trapped molecules began only two decades after the first experimental works. Effective interaction potentials were studied first using se mi -e mpirical a nd ab initio calculations [4,5]. Later the supermolecule approach was applied to examine the effects on the structure and spectra of the ammonia-hydrogen halide complexes of adding 3 Ne or 3 Ar atoms [6,7]. Significant effects were calculated, particularly for the HBr complex. Latajka [8] examined the e f fe ct s o f ad d i n g up to 4 N 2 mo lec u le s to t he ammonia-hydrogen chloride complex and found that 3 N 2 molecules gave similar results to a cavity model in a medium with a relative permittivity of 1.5. Matrix effects on proton transfer in hydrogen -bonded molecular complexes were studied by Barnes [9]. As demonstrated by the recent studies [10][11][12], the matrix affects to the values of the internal rotation barriers in a methanol molecule, that is reflected by changes in splitting of the torsional states if compare with gas phase. Previously [13], in FTIR spectra of methanol recorded at a temperature of 10 K in the argon matrix such an effect has been considered as one of the possible mechanisms responsible for the appearance of a multiplet structure of some absorption bands, in particular a doublet of the bands at 1033.25 and 1036.5 cm-1 caused by the valence vibration of С-О bond. As far as we know, at the moment there are no theoretical works that investigate the influence of the matrix environment on the torsional motion in methanol molecule. In [14] a geometric model of a methanol molecule in the environment of eight argon atoms was proposed and approximately, using a one-dimensional approach, the problem of finding the torsional energy levels was solved. According to the computational results obtained, splitting of the ground torsional state is somewhat increased. Besides, splitting of the Е-type degenerate states occurs due to the relieved degeneracy. However, the more accurate two-dimensional solution is possible and is presented in this paper.

Computation Method
In approximation B3LYP/cc-pVTZ a configuration of the complex including a molecule of methanol and eight argon atoms was optimized using the package GAMESS [15]. The gradient convergence tolerance parameter OPTTOL was set to 10 -5 . No imaginary frequencies were found for optimized configuration. Positions of Ar, C and O atoms were fixed and independent rotation of O-H and CH 3 groups with respect to the argon atoms was performed. Zero values for OH ϕ and 3 CH ϕ were taken for optimized configuration. As it was found that the internal forces in a methanol molecule were greater than the interacting forces with the matrix, on optimization for all the parameters characterizing a position of СН 3 ОН molecule in the matrix, the methyl group rotation relative to the argon atoms was accompanied by the hydroxyl group rotation, and vice versa. Because of this, in the case of a methyl group rotation by steps of 50º the hydroxyl group position was additionally fixed relative to the matrix, and vice versa -for rotation of a hydroxyl group the position of a methyl group was fixed. Besides, all other internal parameters of СН 3 ОН were optimized. Then for every value of ϕ OH = 0º, 50º, 100º, 150º, 200º, 250º, 300º, 350º relative argon atoms the value of 3 CH ϕ = 0º, 50º, 100º, 150º, 200º, 250º, 300º, 350º relative argon atoms too were taken.
As an approximation in the process of the subsequent computations, the geometric parameters of methyl and hydroxyl tops were considered to be constant. According to [14], the moments of inertia are as follows: Because of 120º period for rotation of a methyl group, cloning of the points was performed using the following relation for the internal rotation energy: , an analytical expression for the potential energy of internal rotation was derived in the following form: Substituting (2) into (1), we obtain the analytical expression of the internal rotation energy for methanol in the argon matrix by setting the coordinates s, t:    , but in expression (3) the coefficient of s is not a rational number and internal rotation of a methanol molecule in the argon matrix is aperiodic. At the same time, as seen from Fig. 1, the principal change in the internal energy of the molecule is caused by changing of an internal rotation angle for s.

No Rotation of Methanol Molecule as a Whole with Respect to ArgonMatrix
Experimental spectra for methanol in the argon matrix demonstrate the absence of the rotational structure of vibrational absorption bands as actually there is no rotation of the molecule as a whole, yet exhibiting the bands due to torsional motion in СН 3 ОН. Because of this, it is desirable to consider the case of pure internal rotation in this molecule. According to the computational results, a minimum on the potential surface (2), on the surface ( , ) U s t it is also attained at the point with the coordinates s = t = 0. When there is no rotation of the molecule as a whole, an internal rotation is realized along the axis s. We can derive the form of a curve for the internal rotation potential energy from U(s,t) setting t equal to zero. The curve for U(s,0), its periodic part and perturbation, is given in Fig. 2. A Schrödinger equation for internal rotation takes the form In (4) the function U(s,0) has more than 30 terms. As a perturbation potential we will take all terms in U(s,0) apart from two periodic ones, which are represented by the function U S (s): Equation (4) assumes the form In this case finding of eigenvalues and eigenfunctions in the zeroth approximation is completely equivalent to [14]. Then we calculate the Hamiltonian matrix with the elements of the following form By diagonalization of the Hamiltonian matrix, we have found the energies and wave functions for torsional states of a methanol molecule in the argon matrix. The energy values are listed in Tab. 1; their positions with respect to the potential curve are given in Fig. 3 together with some wave functions.
It is obvious that a change in the molecular dipole moment is due to rotation of a hydroxyl group. As the dipole moment projection onto the С-О axis is invariable with rotation of this group, changing of its components occurs in the plane perpendicular to the rotation axis only. We designate this

Rotation of Methanol Molecule in Argon Matrix with Respect to С-О Axis
If we assume that there is no internal rotation in methanol, but there is rotation of the molecule as whole with respect to the axis coincident with С-О bond, then such a rotation is no longer free in the surroundings of argon atoms. A function for the molecular rotation potential energy may be derived from U(s,t) provided s is zero. The curve for U(0,t) is shown in Fig. 5.
To find the energy levels, we solve a Schrödinger equation of the following form: Potential energy is then of the following form: Solution of (10-12) is similar to that of equations (6,7).
Eigenvalues of the energy for pure rotational levels are listed in Tab. 3. Positions of energy levels and some squared wave functions are demonstrated in Fig. 6 As seen, positions of pure rotational levels are significantly differing from those in the case of free rotation of the fixed-axis rotator for which the following is valid: E k = Dk 2 . (13) Using (13) and considering that D = 4.209 cm -1 , one can easily find the energies of rotational levels for the fixed-axis rotator.

The Case of Motion of Two Types
The results presented in Sections 3 and 4 show that the potential barriers of internal rotation are considerably higher than those associated with rotation of the methanol molecule as a whole relative to С-О bond, in the surroundings of argon atoms. In principle, there are no factors capable to strictly forbid both types of motion simultaneously. Because of this, we must consider this case as well. As the potential energy is a periodic function for both the coordinates OH ϕ and , a problem of two-dimensional large-amplitude motion is readily solved with the use of precisely these coordinates. A Schrödinger equation in this case takes the form: where Let the potential energy be given in the more general form then (1): , , Substituting (15) and (16) into (14), we obtain: 3 3 , (18) Instead of (17), we have: 3 3 , Then we construct the finite matrix with the dimensions 2 2 (2 1) (2 1) c c + × + ; c ∈  . This means that n and m are varying within the limits from -c to c per unity. From (19) we derive: 3 3 , Now we take (20) Fig. 7.
Matrix elements of the dipole moment components can be computed with the use of . The following expression is used: The squared matrix elements of the dipole moment operator are given in Tab. 5. Fig. 8 shows the computed IR absorption spectrum.

Discussion
Despite the fact that in the case of torsional-rotational motion there is a greater number of the computed vibrational frequencies over the spectral region 0-500 cm -1 compared to the case of pure torsional motion, with due regard for the dipole moment matrix elements, the spectra shown in Figs. 4 and 8 are found very similar. Inclusion of the rotation leads to a blue shift of the absorption bands doublet in the region 0-50 cm -1 and to a considerable increase of its relative intensity. A similar blue shift is also observed for the doublet of absorption bands in the interval 200-250 cm -1 . The presence of a high-intensity absorption band in the spectral region 275-300 cm -1 seems to be the principal feature distinguishing the torsional spectrum as compared to the torsional-rotational IR spectrum. With inclusion of a rotation motion, in this interval a low-intensity doublet of absorption bands is observed. Though the proposed model including the argon matrix effect on a torsional-rotational spectrum for a methanol molecule is far from being perfect, it is interesting to correlate the computational results and experimental data. Unfortunately, in [13] IR spectra for a methanol molecule in the argon matrix have been recorded over the region 400-4000 cm -1 . In [17] the spectrum given for methanol in the argon matrix exhibits an absorption band at 271.5 cm -1 and reveals a growth of the intensity in the region of 213 cm -1 , where recording of this spectrum comes abruptly to an end. As a result, computations of normal vibrations were performed considering the only band at 271.5 cm -1 . In [18] an IR spectrum for a methanol molecule in the argon matrix is given in the region 40-4000 cm -1 . Taking into account limitations of the model and the assumptions made in the process of computations, a similarity of the computational results and experimental spectra is remarkable. Some of the absorption bands in the interval 40-120 cm -1 are assigned by the authors as phonon bands. In the spectral region 40-400 cm -1 , apart from these bands, two more absorption bands are observed at 223 and 272 cm -1 , the first of them having a higher intensity. Besides, a half-width of the band at 223 cm -1 is no less than 10 cm -1 , and this is not typical for spectra for the isolated matrix which are recorded at 10 К. The band profile seems to be indicative of a doublet of the overlapping absorption bands. As an intensity of the absorption band at 272 cm -1 is rather high the computed torsional spectrum for a methanol molecule seems to be closer to the experimental one. Thus, the band at 272 cm -1 may be correlated to the computed absorption band at 288.06 cm -1 , and a wide high-intensity band at 223 cm -1 -to a doublet of the computed bands at 204.82 and 217.76 cm -1 . It is possible to assume that a good agreement between the theoretical and experimental spectra is due to the fact that in all the cases we deal with practically undisturbed torsional motion in a molecule of methanol. But, according to [14,[19][20][21][22], a pure torsional spectrum of a methanol molecule recorded at a temperature below 10 К and constructed on the basis of the selection rules ( A A ⇔ , E E ⇔ ) in the spectral interval 0-400 cm -1 should be represented by two absorption bands with the frequencies 294.5 and 353.2 cm -1 . In the case of an insignificant population of the first degenerate torsional state of the E-type there is a possibility for observation of a low-intensity absorption band with the frequency 199.8 cm -1 . It is clear that the matrix effect leads both to the frequency shift and to a significant changing of the transition probability, and this is supported by the presence of two high-intensity absorption bands in an IR spectrum of methanol in the argon matrix over the region 200-300 cm -1 .

Conclusions
In the case when a molecule of methanol is surrounded by Wave number  cm  1  Relative intensities eight argon atoms the authors have considered internal rotation and rotation of the molecule as a whole with respect to the С-О bond, and also their joint motion. Splitting of the deepest double-degenerate torsional state was found to be 29.7 cm -1 . Rotation of the methanol molecule as a whole with respect to С-О bond takes place in a small decelerating potential, height of the barrier being no more than 60 cm -1 .
Because of the elongated form, rotation of a methanol molecule with respect to the axes perpendicular to С-О bond is unlikely. The computations carried out point to the fact that a methanol molecule may be considered as a probe in analysis of the matrix properties. By simulation of the positions of argon atoms surrounding a molecule of methanol we can vary a two-dimensional potential surface and hence a theoretical IR spectrum of the molecule in the matrix. Correlation with the experiment is still required to reveal which of the types of the argon atoms surroundings is most probable.