Vesa Hänninen

Calculation of the O−H Stretching Vibrational Overtone Spectrum of the Water Dimer

The O-H stretching vibrational overtone spectrum of the water dimer has been calculated with the dimer modeled as two individually vibrating monomer units. Vibrational term values and absorption intensities have been obtained variationally with a computed dipole moment surface and an internal coordinate Hamiltonian, which consists of exact kinetic energy operators within the Born-Oppenheimer approximation of the monomer units. Three-dimensional ab initio potential energy and dipole moment surfaces have been calculated using the internal coordinates of the monomer units using the coupled cluster method including single, double, and perturbative triple excitations [CCSD(T)] with the augmented correlation consistent valence triple zeta basis set (aug-cc-pVTZ). The augmented correlation consistent valence quadruple zeta basis set (aug-cc-pVQZ), counterpoise correction, basis set extrapolation to the complete basis set limit, relativistic corrections, and core and valence electron correlations effects have been included in one-dimensional potential energy surface cuts. The aim is both to investigate the level of ab initio and vibrational calculations necessary to produce accurate results when compared with experiment and to aid the detection of the water dimer under atmospheric conditions.

Calculation of spectroscopic parameters and vibrational overtones of methanol

A high-level quartic ab initio potential energy surface of methanol has been used to calculate spectroscopic constants of the 12CH3OH molecule. These include coefficients of quartic anharmonic resonance terms, Darling-Dennison constants, for stretching states. A model expressed in terms of dimensionless normal coordinates has been employed in the calculation of O—H and C—H stretching vibrational states in high-overtone regions. Both cubic Fermi and quartic Darling-Dennison anharmonic coupling terms have been included in the model in order to take into account strong resonances between different states. The nonlinear least-squares method has been used to optimize some of the model parameters employing experimental term values of 12CH3OH as data. Vibrational assignments are suggested for the first C—H stretching overtone region.

Anharmonic force field for methanol

An ab initio quartic anharmonic force field for methanol has been calculated at the equilibrium position using the CCSD(T) method for the structure and the harmonic potential energy surface, and the MP4(SDQ) method for the anharmonic part of the surface. A triple zeta basis set was employed with symmetrized curvilinear internal valence coordinates in all calculations. The internal coordinate force field constants have been transformed into force constants in the dimensionless normal coordinate representation for various isotopomers. Vibrational term values for CH3OH, CH3OD, CD3OH, and CD3OD have been obtained using second order perturbation theory. Particular care has been devoted to the inclusion of Fermi resonance interactions between different vibrational states. A good accuracy has been achieved in the calculation of the fundamentals for all the isotopomers, the mean absolute error being 5.8 cm−1.

The effect of large amplitude motions on the transition frequency redshift in hydrogen bonded complexes: A physical picture

We describe the vibrational transitions of the donor unit in water dimer with an approach that is based on a three-dimensional local mode model. We perform a perturbative treatment of the intermolecular vibrational modes to improve the transition wavenumber of the hydrogen bonded OH-stretching transition. The model accurately predicts the transition wavenumbers of the vibrations in water dimer compared to experimental values and provides a physical picture that explains the redshift of the hydrogen bonded OH-oscillator. We find that it is unnecessary to include all six intermolecular modes in the vibrational model and that their effect can, to a good approximation, be computed using a potential energy surface calculated at a lower level electronic structure method than that used for the unperturbed model.

Torsional motion and vibrational overtone spectroscopy of methanol

An internal coordinate Hamiltonian model has been constructed to model torsional motion in the OH stretching vibrational overtone region of methanol, CH3OH. The model includes harmonic couplings between OH and CH stretching vibrations and Fermi resonance interactions between OH stretches and COH bends and between CH stretches and CH2 bends. A symmetrized basis set has been used to form block diagonal Hamiltonian matrices with strong resonance couplings off-diagonal. Observed torsional levels of the excited vibrational states have been used as data in a least squares optimization of the model parameters, some of which have been estimated by ab initio calculations. The experimentally observed increase in the effective torsional barrier in moving to highly excited OH stretching states has been explained by the model.

The effect of large amplitude motions on the vibrational intensities in hydrogen bonded complexes.

We have developed a model to calculate accurately the intensity of the hydrogen bonded XH-stretching vibrational transition in hydrogen bonded complexes. In the Local Mode Perturbation Theory (LMPT) model, the unperturbed system is described by a local mode (LM) model, which is perturbed by the intermolecular modes of the hydrogen bonded system that couple with the intramolecular vibrations of the donor unit through the potential energy surface. We have applied the model to three complexes containing water as the donor unit and different acceptor units, providing a series of increasing complex binding energy: H2O⋯N2, H2O⋯H2O, and H2O⋯NH3. Results obtained by the LMPT model are presented and compared with calculated results obtained by other vibrational models and with previous results from gas-phase and helium-droplet experiments. We find that the LMPT model reduces the oscillator strengths of the fundamental hydrogen bonded OH-stretching transition relative to the simpler LM model.

Ab Initio Structural and Vibrational Investigation of Sulfuric Acid Monohydrate

We employ ab initio methods to find stable geometries and to calculate potential energy surfaces and vibrational wavenumbers for sulfuric acid monohydrate. Geometry optimizations are carried out with the explicitly correlated coupled-cluster approach that includes single, double, and perturbative triple excitations (CCSD(T)-F12a) with a valence double-ζ basis set (VDZ-F12). Four different stable geometries are found, and the two lowest are within 0.41 kJ mol(-1) (or 34 cm(-1)) of each other. Vibrational harmonic wavenumbers are calculated at both the density-fitted local spin component scaled second-order Møller-Plesset perturbation theory (DF-SCS-LMP2) with the aug-cc-pV(T+d)Z basis set and the CCSD-F12/VDZ-F12 level. Water O-H stretching vibrations and two highly anharmonic large-amplitude motions connecting the three lowest potential energy minima are considered by limiting the dimensionality of the corresponding potential energy surfaces to small two- or three-dimensional subspaces that contain only strongly coupled vibrational degrees of freedom. In these anharmonic domains, the vibrational problem is solved variationally using potential energy surfaces calculated at the CCSD(T)-F12a/VDZ-F12 level.

Simulation of inversion motion and N–H stretching overtone spectra of aniline

A curvilinear internal coordinate Hamiltonian is used to simulate the N-H stretching overtone spectra and the associated inversion splittings in aniline. A simple local mode type model is applied to the N-H stretching and H-N-H bending modes. Geometric algebra is employed to derive the kinetic energy operator for the large amplitude inversion motion. Electronic structure calculations at the Moller-Plesset second order perturbation theory and correlation consistent aug-cc-pVTZ basis set level are used to obtain model parameters, some of which have been optimized with the least-squares method using experimental vibrational term values as data. The observed N-H stretching overtone vibrational levels and the inversional tunneling splittings are well reproduced with our approach.

The (n00),n=3, 4, and 6, local mode states of H3SiD: Fourier transform infrared and laser photoacoustic spectra and ab initio calculations of spectroscopic parameters

The rotational structure of the local mode Si–H stretching vibrational bands (n00 A1/E), n=3, 4, and 6, of H328SiD have been studied by high-resolution Fourier transform infrared and by photoacoustic laser spectroscopy. The recorded bands have been rotationally analyzed with a Hamiltonian model which makes use of simple arithmetic relations between some of the rovibrational parameters. While the (300 A1/E) states were found to be unperturbed, severe perturbations by unknown dark states affect the (400 A1/E) and (600 A1/E) states for J values exceeding 8. Ab initio calculations have been performed to form the quadratic and the cubic potential energy surfaces which have been used to calculate spectroscopic parameters for the Si–H stretching fundamentals. These results, together with the local mode relations, have been successfully used to model the vibrational dependence of effective rovibrational parameters in the excited local mode states.

Dispersion Interactions in Small Zinc, Cadmium, and Mercury Clusters

Dispersion interactions in small clusters of group XII (Zn, Cd, Hg) metal atoms are studied at the CCSD(T) level with triple-ζ basis sets. A pair potential model together with a least-squares fit to the interaction potential energy surface is used to calculate interatomic dispersion coefficients, which are found to be in good agreement with atomistic calculations. The angular dependence of the dispersion interaction, extracted by explicitly accounting for the leading order electrostatic and induction terms, is determined with the aid of a coupled spherical harmonic expansion. The distance dependence of the orientation averaged dispersion energy is modeled by an integral average of our potential model, which is able to robustly account for the average dispersion energy to a remarkable degree.