Galina M. Chaban

Calculation of Vibrational Transition Frequencies and Intensities in Water Dimer: Comparison of Different Vibrational Approaches

We have calculated frequencies and intensities of fundamental and overtone vibrational transitions in water and water dimer with use of different vibrational methods. We have compared results obtained with correlation-corrected vibrational self-consistent-field theory and vibrational second-order perturbation theory both using normal modes and finally with a harmonically coupled anharmonic oscillator local mode model including OH-stretching and HOH-bending local modes. The coupled cluster with singles, doubles, and perturbative triples ab initio method with augmented correlation-consistent triple-zeta Dunning and atomic natural orbital basis sets has been used to obtain the necessary potential energy and dipole moment surfaces. We identify the strengths and weaknesses of these different vibrational approaches and compare our results to the available experimental results.

Degenerate perturbation theory corrections for the vibrational self-consistent field approximation: Method and applications

A new algorithm for computing anharmonic vibrational states for polyatomic molecules is proposed. The algorithm starts with the vibrational self-consistent field (VSCF) method and uses degenerate perturbation theory to correct for effects of correlation between different vibrational modes. The algorithm is developed in a version that computes the anharmonic vibrational spectroscopy directly from potential energy surface points calculated by using ab initio codes. The method is applied to several molecules where near degeneracies occur for excited vibrational states, including HOOH, HSSH, and HOOOH. The method yields results in very good accordance with experiments and generally provides improvements over nondegenerate perturbation corrections for VSCF.

Combined ab initio and anharmonic vibrational spectroscopy calculations for rare gas containing fluorohydrides, HRgF

Abstract MP2 and CCSD(T) calculations are used to analyse the structures and vibrational spectra of HRgF molecules, where the rare gas atom is He, Ne, Ar, Kr, Xe or Rn. We extend the analysis of the vibrational spectra of these molecules to include anharmonic corrections for the most likely candidates for experimental detection, i.e., HArF, HKrF, HXeF, and their deuterated isotopomers. The anharmonic correlation-corrected vibrational self-consistent-field (CC-VSCF) calculations are used for this, and fundamental, overtone and combination frequencies and their absorption intensities are computed.

Vibrational spectroscopy and matrix-site geometries of HArF, HKrF, HXeCl, and HXeI in rare-gas solids

The vibrational spectroscopy and the matrix-site geometries of several novel rare-gas compounds in the matrix environment were computed theoretically, and compared with experiment. Ab initio calculations are used in the fitting of analytical potential surfaces for the HRgY molecules and for the interactions between HRgY and the matrix atoms Rg. With these potentials, matrix-site geometries for the molecule in the solid are computed. Finally, the vibrational spectroscopy of HRgY in the Rg matrix is computed using the vibrational self-consistent field (VSCF) method. The VSCF includes anharmonic effects, that are essential in this case. The version of VSCF used here includes coupling between HRgY and the vibrations of the solid atoms. The vibrations of 72 matrix atoms are treated. The main results are: (1) The matrix shifts are considerably greater than typically found for neutral, strongly bond molecules, but are much smaller than discrepancies between theory and experiment. This can be attributed to the in...

A natural orbital diagnostic for multiconfigurational character in correlated wave functions

The natural orbitals and their corresponding occupation numbers are constructed for several interesting problems to demonstrate that the existence of negative natural orbital occupation numbers for single reference correlation methods provides a simple diagnostic for the need for a multiconfigurational description of the wave function.

Theoretical study of decomposition pathways for HArF and HKrF

To provide theoretical insights into the stability and dynamics of the new rare gas compounds HArF and HKrF, reaction paths for decomposition processes HRgF ! Rg þ HF and HRgF ! H þ Rg þ F (Rg ¼ Ar, Kr) are calculated using ab initio electronic structure methods. The bending channels, HRgF ! Rg þ HF, are described by single-configurational MP2 and CCSD(T) electronic structure methods, while the linear decomposition paths, HRgF ! H þ Rg þ F, require the use of multi-configurational wave functions that include dynamic correlation and are size extensive. HArF and HKrF molecules are found to be energetically stable with respect to atomic dissociation products (H + Rg + F) and separated by substantial energy barriers from Rg + HF products, which ensure their kinetic stability. The results are compatible with experimental data on these systems. 2002 Elsevier Science B.V. All rights reserved.

The Transition from Hydrogen Bonding to Ionization in (HCI)n(NH3)n and (HCI)n(H2O)n Clusters: Consequences for Anharmonic Vibrational Spectroscopy

Anharmonic vibrational frequencies and intensities are calculated for 1:1 and 2:2 (HCl)n(NH3)n and (HCl)n(H2O)n complexes, employing the correlation-corrected vibrational self-consistent field method with ab initio potential surfaces at the MP2/TZP computational level. In this method, the anharmonic coupling between all vibrational modes is included, which is found to be important for the systems studied. For the 4:4 (HCl)n(H2O)n complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Just as the (HCl)n(NH3)n structure switches from hydrogen-bonded to ionic for n = 2, the (HCl)n(H2O)n switches to ionic structure for n = 4. For (HCl)2(H2O)2, the lowest energy structure corresponds to the hydrogen-bonded form. However, configurations of the ionic form are separated from this minimum by a barrier of less than an O−H stretching quantum. This suggests the possibility of experiments on ionization dynamics using infrared excitation of the hydrogen-bonded form....

Anharmonic vibrational spectroscopy of the glycine–water complex: Calculations for ab initio, empirical, and hybrid quantum mechanics/molecular mechanics potentials

Effects of intermolecular hydrogen bonding between glycine and one water molecule on the vibrational spectrum are investigated, using ab initio (at the level of second order Moller–Plesset perturbation theory), empirical (OPLS-AA), and mixed ab initio/empirical quantum mechanics/molecular mechanics (QM/MM) potentials. Vibrational spectroscopy is calculated using the correlation corrected vibrational self-consistent field method that accounts for anharmonicities and couplings between different vibrational normal modes. The intermolecular hydrogen bonding interactions are found to be very strong and to affect vibrational frequencies and infrared intensities of both the glycine and the water molecule to a very large extent. The predicted ab initio anharmonic spectra can be used to identify amino acids in complexes with water in experimental studies. The OPLS-AA potential is found to describe hydrogen bonding between glycine and water incorrectly, and to predict erroneous vibrational spectra. Hybrid (QM/MM) t...

Spectroscopically-tested, improved, semi-empirical potentials for biological molecules: Calculations for glycine, alanine and proline

A modification of the semi-empirical PM3 electronic structure method is proposed. It employs a coordinate scaling procedure, such that the harmonic frequencies from the modified PM3 potentials for lower-energy conformers of glycine (conformer I), alanine (conformers I and II) and proline (conformer II), fit more closely with ab initio (MP2/DZP) harmonic frequencies. The anharmonic frequencies are then calculated using the modified PM3 surfaces with the Vibrational Self-Consistent Field (VSCF) and Correlation-Corrected VSCF (CC-VSCF) methods. The computed anharmonic frequencies are in very good accord with spectroscopic experiments for the three amino acids. The results are much superior to those obtained from standard (unscaled) PM3 potentials, indicating that the modified PM3 potentials may be used as high quality potentials for biological molecules, at least in the configuration ranges pertinent to vibrational spectroscopy. The scaling parameters computed for the lowest energy conformers listed above were tested for transferability: they were used in computing the anharmonic spectra of two other conformers (glycine II and proline I). The good agreement of the resulting frequencies with observed frequencies, indicates the transferability of the scaling parameters. It is concluded from this study that the improved PM3 potentials offer accurate and computationally efficient force fields for vibrational spectroscopy calculations of biological molecules. Possible additional applications of the new potentials are discussed.

Vibrational spectroscopy of (SO42−)∙(H2O)n clusters, n=1–5: Harmonic and anharmonic calculations and experiment

The vibrational spectroscopy of (SO4(2-)).(H2O)n is studied by theoretical calculations for n=1-5, and the results are compared with experiments for n=3-5. The calculations use both ab initio MP2 and DFT/B3LYP potential energy surfaces. Both harmonic and anharmonic calculations are reported, the latter with the CC-VSCF method. The main findings are the following: (1) With one exception (H2O bending mode), the anharmonicity of the observed transitions, all in the experimental window of 540-1850 cm(-1), is negligible. The computed anharmonic coupling suggests that intramolecular vibrational redistribution does not play any role for the observed linewidths. (2) Comparison with experiment at the harmonic level of computed fundamental frequencies indicates that MP2 is significantly more accurate than DFT/B3LYP for these systems. (3) Strong anharmonic effects are, however, calculated for numerous transitions of these systems, which are outside the present observation window. These include fundamentals as well as combination modes. (4) Combination modes for the n=1 and n=2 clusters are computed. Several relatively strong combination transitions are predicted. These show strong anharmonic effects. (5) An interesting effect of the zero point energy (ZPE) on structure is found for (SO4(2-)).(H2O)(5): The global minimum of the potential energy corresponds to a C(s) structure, but with incorporation of ZPE the lowest energy structure is C2v, in accordance with experiment. (6) No stable structures were found for (OH-).(HSO4-).(H2O)n, for n<or=5.