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.     

Theoretical calculation of positron annihilation spectrum using positron-electron correlation-polarization potential

The positron-electron correlation-polarization potential model is used to calculate annihilation spectra of carbon disulfide and benzene. We assume that the positron is captured in the vibrationally excited states of the target molecule through vibrational Feshbach resonances. Using the standard normal mode representation, we calculated the resonance energies and widths for each vibrational mode. The resonance widths were calculated with Fermi's Golden Rule approximation, where the time-dependent wave packet approach has been applied. We found that vibrational resonances of infrared-active modes play a dominant role in resonant annihilation; however, infrared-inactive modes also contribute to the annihilation spectrum through polarizability changes along normal mode coordinates.     

Ab initio calculation of the rotational spectrum of methane vibrational ground state

In a previous article we have introduced an alternative perturbation scheme to the traditional one starting from the harmonic oscillator, rigid rotator Hamiltonian, to find approximate solutions of the spectral problem for rotation-vibration molecular Hamiltonians. The convergence of our method for the methane vibrational ground state rotational energy levels was quicker than that of the traditional method, as expected, and our predictions were quantitative. In this second article, we study the convergence of the ab initio calculation of effective dipole moments for methane within the same theoretical frame. The first order of perturbation when applied to the electric dipole moment operator of a spherical top gives the expression used in previous spectroscopic studies. Higher orders of perturbation give corrections corresponding to higher centrifugal distortion contributions and are calculated accurately for the first time. Two potential energy surfaces of the literature have been used for solving the anharmonic vibrational problem by means of the vibrational mean field configuration interaction approach. Two corresponding dipole moment surfaces were calculated in this work at a high level of theory. The predicted intensities agree better with recent experimental values than their empirical fit. This suggests that our ab initio dipole moment surface and effective dipole moment operator are both highly accurate.     

Vibrational and correlation effects on properties of water

Large basis set, ab initio potential energy and property surfaces of water have been used with quantum Monte Carlo vibrational analysis in the evaluation of the molecule's rotational constants, zero-point energy, and dipole moment. While there are clearly differences in vibrational state parameters due to including correlation effects, the vibrational averaging effect on rotational constants is very nearly additive with the correlation effect. This has implications for evaluation and estimation of properties of molecules in specific vibrational states. Received: 9 February 1999 / Accepted: 11 February 1999 / Published online: 5 May 1999     

Automatic generation of force fields and property surfaces for use in variational vibrational calculations of anharmonic vibrational energies and zero-point vibrational averaged properties

An automatic and general procedure for the calculation of geometrical derivatives of the energy and general property surfaces for molecular systems is developed and implemented. General expressions for an n-mode representation are derived, where the n-mode representation includes only the couplings between n or less degrees of freedom. The general expressions are specialized to derivative force fields and property surfaces, and a scheme for calculation of the numerical derivatives is implemented. The implementation is interfaced to electronic structure programs and may be used for both ground and excited electronic states. The implementation is done in the context of a vibrational structure program and can be used in combination with vibrational self-consistent field (VSCF), vibrational configuration interaction (VCI), vibrational Moller-Plesset, and vibrational coupled cluster calculations of anharmonic wave functions and calculation of vibrational averaged properties at the VSCF and VCI levels. Sample calculations are presented for fundamental vibrational energies and vibrationally averaged dipole moments and frequency dependent polarizabilities and hyperpolarizabilities of water and formaldehyde.     

A virtual vibrational self-consistent-field method for efficient calculation of molecular vibrational partition functions and thermal effects on molecular properties

A new method is described for the calculation of molecular vibrational partition functions and thermal effects on molecular properties including an explicit account of anharmonicity. The approach is based on the vibrational self-consistent-field method. Partition functions and thermal averages of the energies calculated with the new method are generally in good agreement with the result of more accurate methods. At lower temperatures the method gives in addition good results for thermal averages of dipole moments and polarizabilities. The new method is much more efficient than explicit sum-over-states approaches previously used for calculation of thermal averages. Unlike the standard sum-over-states approach, the newly developed method is feasible for larger systems despite the formal exponential increase in the number of states with the size of the system. Thus, it is presently the only practical way for including an explicit treatment of anharmonicity in vibrational wave function based calculations of molecular vibrational partition functions and thermally averaged properties of larger molecules.     

The predicted spectrum of the hypermetallic molecule MgOMg

The present study of MgOMg is a continuation of our theoretical work on Group 2 M(2)O hypermetallic oxides. Previous ab initio calculations have shown that MgOMg has a linear (1)Σ(g)+ ground electronic state and a very low lying first excited triplet electronic state that is also linear; the triplet state has (3)Σ(u)+ symmetry. No gas phase spectrum of this molecule has been assigned, and here we simulate the infrared absorption spectrum for both states. We calculate the three-dimensional potential energy surface, and the electric dipole moment surfaces, of each of the two states using a multireference configuration interaction (MRCISD) approach based on full-valence complete active space self-consistent field (FV-CASSCF) wavefunctions with a cc-pCVQZ basis set. A variational MORBID calculation using our potential energy and dipole moment surfaces is performed to determine rovibrational term values and to simulate the infrared absorption spectrum of the two states. We also calculate the dipole polarizability of both states at their equilibrium geometry in order to assist in the interpretation of future beam deflection experiments. Finally, in order to assist in the analysis of the electronic spectrum, we calculate the vertical excitation energies, and electric dipole transition matrix elements, for six excited singlet states and five excited triplet states using the state-average full valence CASSCF-MRCISD/aug-cc-pCVQZ procedure.     

Thermally induced polarizabilities and dipole moments of small tin clusters

We study the influence of thermal excitation on the electric susceptibilities for Sn(6) and Sn(7) clusters by molecular beam electric deflection and Monte-Carlo simulations in conjunction with quantum-chemical calculations. At low temperatures (40 K), no field-induced broadening of the Sn(6) and Sn(7) cluster beams are observed, in agreement with vanishing permanent electric dipole moments due to their centro-symmetrical ground states. The electric polarizabilities of Sn(6) and Sn(7), as inferred from the field-induced molecular beam deflection, are in good agreement with the quantum-chemical predictions. At elevated temperatures of 50-100 K, increased polarizabilities of about 2-3 ?(3) are obtained. Also, we found indications of a field-induced beam broadening which points to the existence of permanent dipole moments of about 0.01-0.02 D per atom at higher temperatures. These results cannot be explained by thermal excitations within a harmonic oscillator model, which would yield a temperature-independent polarizability and fluxional, but not permanent, dipole moments. We analyze this behavior by Monte-Carlo simulations in order to compute average temperature-induced electric dipole moments. For that purpose, we developed a novel technique for predicting observables sampled on the quantum-chemical potential energy surface by an umbrella sampling correction of Monte-Carlo results obtained from simulations utilizing an empirical potential. The calculated, fluxional dipole moments are in tune with the observed beam broadenings. The cluster dynamics underlying the polarizability appear to be intermediate between rigid and floppy molecules which leads to the conclusion that the rotational, not the vibrational temperature seems to be the key parameter that determines the temperature dependence of the polarizability.     

Theoretical rate constants and kinetic isotope effects in the reaction of methane with H, D, T, and Mu atoms

The rate constants and kinetic isotope effects of the reaction of methane with four isotopes of hydrogen, protium (H), deuterium (D), tritium (T), and muonium (Mu), were studied using variational transition state theory with multidimensional tunneling on an analytical potential energy surface, PES-2002, previously constructed by our group. For the four isotopes, our kinetics results agree reasonably with available experimental measurements, improving previous theoretical results that used different potential energy surfaces and/or theoretical approaches. In the comparison of the reactivity between protium and muonium, which is the most severe test of the surface and theoretical method due to the large mass difference between the two isotopes, some sources of discrepancy between theory and experiment were analyzed. These were the zero-point energy, tunneling effect, and the role of the reactivity from methane excited vibrational states.     

A surface hopping method for chemical reaction dynamics in solution described by diabatic representation: an analysis of tunneling and thermal activation

We present a surface hopping method for chemical reaction in solution based on diabatic representation, where quantum mechanical time evolution of the vibrational state of the reacting nuclei as well as the reaction-related electronic state of the system are traced simultaneously together with the classical motion of the solvent. The method is effective in describing the system where decoherence between reactant and product states is rapid. The diabatic representation can also give a clear picture for the reaction mechanism, e.g., thermal activation mechanism and a tunneling one. An idea of molecular orbital theory has been applied to evaluate the solvent contribution to the electronic coupling which determines the rate of reactive transition between the reactant and product potential surfaces. We applied the method to a model system which can describe complex chemical reaction of the real system. Two numerical examples are presented in order to demonstrate the applicability of the present method, where the first example traces a chemical reaction proceeded by thermal activation mechanism and the second examines tunneling mechanism mimicking a proton transfer reaction.     

CASSCF and CASPT2 studies on the structures, transition energies, and dipole moments of ground and excited states for azulene

The electronic structure of azulene molecule has been studied. We have obtained the optimized structures of ground and singlet excited states by using the complete active space self-consistent-field (CASSCF) method, and calculated vertical and 0-0 transition energies between the ground and excited states with second-order M?ller-Plesset perturbation theory (CASPT2). The CASPT2 calculations indicate that the bond-equalized C(2v) structure is more stable than the bond-alternating C(s) structure in the ground state. For a physical understanding of electronic structure change from C(2v) to C(s), we have performed the CASSCF calculations of Duschinsky matrix describing mixing of the b(2) vibrational mode between the ground (1A(1)) and the first excited (1B(2)) states based on the Kekule-crossing model. The CASPT2 0-0 transition energies are in fairly good agreement with experimental results within 0.1-0.3 eV. The CASSCF oscillator strengths between the ground and excited states are calculated and compared with experimental data. Furthermore, we have calculated the CASPT2 dipole moments of ground and excited states, which show good agreement with experimental values.     

Time-resolved resonance Raman spectroscopic studies on the triplet excited state of thioxanthone

Thioxanthone has been investigated extensively owing to its unique photochemical and photophysical applications and its solvatochromic behavior. Here, we report the time-resolved resonance Raman studies on the structure of the lowest triplet excited state of thioxanthone in carbon tetrachloride. In addition, FT-IR and FT-Raman techniques have been used to study the vibrational structure in the ground state. To corroborate the experimental findings, density functional theory calculations have been carried out. Isotopic calculations and normal coordinate analysis have been used to help in assigning the observed bands to Raman vibrational modes. Structural information derived from this study is expected to help in better understanding the triplet state photochemistry of thioxanthone.     

Calculations of static dipole polarizabilities of alkali dimers: prospects for alignment of ultracold molecules

The rapid development of experimental techniques to produce ultracold alkali molecules opens the ways to manipulate them and to control their dynamics using external electric fields. A prerequisite quantity for such studies is the knowledge of their static dipole polarizability. In this paper, we computed the variations with internuclear distance and with vibrational index of the static dipole polarizability components of all homonuclear alkali dimers including Fr(2), and of all heteronuclear alkali dimers involving Li to Cs, in their electronic ground state and in their lowest triplet state. We use the same quantum chemistry approach as in our work on dipole moments [Aymar and Dulieu, J. Chem. Phys. 122, 204302 (2005)], based on pseudopotentials for atomic core representation, Gaussian basis sets, and effective potentials for core polarization. Polarizabilities are extracted from electronic energies using the finite-field method. For the heaviest species Rb(2), Cs(2), and Fr(2) and for all heteronuclear alkali dimers, such results are presented for the first time. The accuracy of our results on atomic and molecular static dipole polarizabilities is discussed by comparing our values with the few available experimental data and elaborate calculations. We found that for all alkali pairs, the parallel and perpendicular components of the ground state polarizabilities at the equilibrium distance R(e) scale as (R(e))(3), which can be related to a simple electrostatic model of an ellipsoidal charge distribution. Prospects for possible alignment and orientation effects with these molecules in forthcoming experiments are discussed.     

Vibrational contributions to indirect spin-spin coupling constants calculated via variational anharmonic approaches

Zero-point vibrational contributions to indirect spin-spin coupling constants for N2, CO, HF, H2O, C2H2, and CH4 are calculated via explicitly anharmonic approaches. Thermal averages of indirect spin-spin coupling constants are calculated for the same set of molecules and for C2X4, X = H, F, Cl. Potential energy surfaces have been calculated on a grid of points and analytic representations have been obtained by a linear least-squares fit in a direct product polynomial basis. Property surfaces have been represented by a fourth-order Taylor expansion around the equilibrium geometry. The electronic structure calculations employ density functional theory, and vibrational contributions to indirect spin-spin coupling constants are calculated employing vibrational self-consistent-field and vibrational configuration-interaction methods. The performance of vibrational perturbation theory and various approximate variational calculations are discussed. Thermal averages are computed by state-specific and virtual vibrational self-consistent-field methods.     

A gradient extremal walking algorithm

Gradient extremals define stream beds connecting stationary points on molecular potential energy surfaces. Using this concept we have developed an algorithm to determine transition states. We initiate walks at equilibrium geometries and follow the gradient extremals until a stationary point is reached. As an illustration we have investigated the mechanism for exchange of protons on carbon in methylenimine (H₂C=NH) using a multi-reference self-consistent-field wave function.     

On the performance of molecular polarization methods. I. Water and carbon tetrachloride close to a point charge

The three main methods to implement molecular polarization (point dipoles, fluctuating charges, and shell model) are tested against high level ab initio calculations for a molecule (water, carbon tetrachloride) close to a point charge (at the distance of a lithium or magnesium ion). The goal is to check whether an approximation (linear polarization) strictly valid at large intermolecular distances is sufficiently accurate for liquid state molecular dynamics simulations, where strong polarization effects are to be expected at short separations. The monitored observable is the molecular dipole moment as a function of the charge-molecule distance for selected molecular orientations. Analytic formulas are derived for the components of the molecular polarization tensor, facilitating the optimization of the performance for each polarization method as a function of its underlying parameters. Overall, the methods studied provide a remarkably good representation of the induced dipole, with no divergences appearing even at the shortest distances. For water close to a monovalent point charge the point dipole model, implemented with one or three dipoles, accurately reproduces the water dipole moment at all distances. Deficiencies appear as the molecular polarizability and/or charge increase: basically, the ab initio induced moments grow faster at intermediate distances than the linear increase characteristic of the phenomenological polarization methods, suggesting that nonlinear effects (hyperpolarizability) cannot be neglected in these cases. Regarding the capabilities of each method, the point dipole method is the one that performs best overall, with the shell model achieving acceptable results in most instances. The fluctuating charge method shows some noticeable limitations for implementations of comparable complexity (in terms of the number of sites required).     

Ab initio studies on hydrogen-transfer tunneling for Cl + HCl abstraction hydrogen reaction

This article presents a treatment scheme of the tunneling of hydrogen between two molecular centers (Cl…Cl). The purpose is to calculate the tunneling probabilities of hydrogen atom transfer from the initial (the proceeding complex) to the final-state energy minima (the succeeding complex) in two anharmonic vibrational states (0 → 0 and 1 → 1) in terms of the time-dependent perturbation theory expression and to see whether spectroscopic signatures of tunneling persist in the form of splittings of the vibrational modes. The analysis uses the realistic potential energy function calculated at the HF/6−31 + G** self-consistent-field basis-set level for the interaction between transferred hydrogen and its molecular skeleton (Cl…H…Cl). This potential energy surface is calibrated by comparing its properties with those from sf-POLCI and the LEPS potential-energy surfaces. The anharmonic vibrational state is characterized by the corrected vibrational energy levels and a set of linear combination coefficients obtained via perturbation theory. The tunneling probabilities for two transitions (0 → 0 and 1 → 1) were calculated and compared with those from Gamow's equation. Applicability of the time-dependent perturbation theory expression and Gamow's equation to the [Cl BOND H…Cl] system is discussed. The vibrational splitting energies are obtained, and a spectroscopic signature caused by tunneling is expected and should be observable. © 1996 John Wiley & Sons, Inc.     

Molecular dynamic structure and dimerization of polyacetylene

The vibrational self-consistent-field approximation is used to calculate excited vibrational energy levels of the water molecule in hyperspherical coordinates. The calculations are made for a global realistic Sorbie–Murrell-type potential surface for which exactum quantum variational results are known for comparison. The coupled SCF equations are solved using the discrete variable representation (DVR ) method, which allows computation of the coupled multidimensional integrals in a very simple and efficient way. The results are in good agreement with exactum quantum calculations and are more accurate than SCF energy eigenvalues obtained using normal mode coordinates.     

Communication: Multi-state analysis of the OCS ultraviolet absorption including vibrational structure

The first absorption band of OCS (carbonyl sulfide) is analyzed using potential energy surfaces and transition dipole moment functions of the lowest four singlet and the lowest four triplet states. Excitation of the 2 (1)A' state is predominant except at very low photon energies. It is shown that the vibrational structures in the center of the band are due to excitation of the 2 (3)A' triplet state, whereas the structures at very low energies are caused by bending excitation in the potential wells of states 2 (1)A' and 1 (1)A'.