Accurate ab initio pair potentials between helium and the heavier group 2 elements

Interactions between the heavier Group 2 metals (Ca, Sr, and Ba) and helium were studied using the well-tempered model core potential method. Accurate pair potentials, calculated at the coupled-cluster level of theory with very large basis sets, were used in bound state calculations. Three bound rovibrational states were found for each complex. The pair-potential parameters were used to predict how each of the metal atoms would be solvated by a helium nanodroplet. The Ca atom is not fully solvated by the droplet and the interaction between the helium and the metal decreases from Ca to Ba. This agrees with the experimental observation that the spectra of these atoms in a nanodroplet are intermediate between the spectra of the free atoms and the spectra in liquid helium.     

Interactions of transition metal atoms in high-spin states: Cr2, Sc-Cr, and Sc-Kr

The high-spin van der Waals states are examined for the following dimers: Cr(2) ((13)Sigma(g)(+)), Sc-Cr ((8)Sigma(+), (8)Pi, (8)Delta), and Sc-Kr ((2)Sigma(+), (2)Pi, (2)Delta). These three systems offer a wide range of van der Waals interactions: anomalously strong, intermediate, and typically weak. The single-reference [coupled cluster with single, double, and noniterative triple excitations, RCCSD(T)] method is used in the calculations for all three systems. In addition, a range of configuration-interaction based methods is applied in Cr(2) and Sc-Cr. The three dimers are shown to be bound by the dispersion interaction of varying strength. In a related effort, the dispersion energy and its exchange counterpart are calculated using the newly developed open-shell variant of the symmetry-adapted perturbation theory (SAPT). The restricted open-shell time-dependent Hartree-Fock linear response function is used in the calculations of the dispersion energy in Sc-Cr and Sc-Kr calculations, while the restricted open-shell time-dependent density functional linear response function is used for Cr(2). A hybrid method combining the repulsive restricted open-shell Hartree-Fock (or complete active space self-consistent field) interaction energy with the dispersion and exchange-dispersion terms is tested against the RCCSD(T) results for the three complexes. The Cr(2) ((13)Sigma(g)(+)) complex has the well depth of 807.8 cm(-1) at the equilibrium distance of 6.18a(0) and the dissociation energy of 776.8 cm(-1). The octet-state Sc-Cr is about four times more strongly bound with the order of well depths of (8)Delta>(8)Pi>(8)Sigma(+) and a considerable anisotropy. The enhanced bonding is attributed to the unusually strong dispersion interaction. Sc-Kr ((2)Sigma(+), (2)Pi, (2)Delta) is a typical van der Waals dimer with well depths in the range of 81 cm(-1) ((2)Delta), 84 cm(-1) ((2)Sigma(+)), and 86 cm(-1) ((2)Pi). The hybrid model based on SAPT leads to results which are in excellent qualitative agreement with RCCSD(T) for all three interactions.     

Infrared spectra of isotopic CO2-He complexes

Infrared spectra of three isotopic forms of the weakly bound CO(2)-He van der Waals complex have been studied in the region of the CO(2) nu(3) fundamental band around 2300 cm(-1), using a tunable diode laser to probe a pulsed supersonic expansion. The complex is a T-shaped near-oblate asymmetric rotor, and it is found that (18)O isotopic substitution is sufficient to interchange the a and b inertial axes. For the symmetric isotopes, such as the normal species and the (16)O(13)C(16)O and (18)O(13)C(18)O forms studied here, half of the normal rotational levels of the complex are missing due to the effects of (16)O (or (18)O) interchange symmetry. However, for asymmetrically substituted ones, such as (16)O(13)C(18)O, all rotational levels are present. Moreover, for the asymmetric isotope, both a- and b-type transitions were observed, so that the spectrum was much richer. The CO(2)-He system is of interest both as a benchmark for intermolecular potential energy surface calculations, and because CO(2) is a valuable probe molecule for helium cluster spectroscopy.     

Vibrational energies for the X1A1, A1B1, and B1A1 states of SiH2/SiD2 and related transition probabilities based on global potential energy surfaces

Transition probabilities were evaluated for the X(1)A(1)-A(1)B(1) and A(1)B(1)-B(1)A(1) systems of SiH(2) and SiD(2) to analyze the X-->A-->B photoexcitation. The Franck-Condon factors (FCFs) and Einstein's B coefficients were computed by quantum vibrational calculations using the three-dimensional potential energy surfaces (PESs) of the SiH(2)(X(1)A(1),A(1)B(1),B(1)A(1)) electronic states and the electronic transition moments for the X-A, X-B, and A-B system. The global PESs were determined by the multireference configuration interaction calculations with the Davidson correction and the interpolant moving least-squares method combined with the Shepard interpolation. The obtained FCFs for the X-A and A-B systems exhibit that the bending mode is strongly enhanced in the excitation since the equilibrium bond angle greatly varies with the three states; the barrier to linearity is evaluated to be 21,900 cm(-1) for the X state, 6400 cm(-1) for the A state, and 230-240 cm(-1) for the B state. The theoretical lifetimes for the pure bending levels of the A and B states were calculated from the fluorescence decay rates for the A-X, B-A, and B-X emissions.     

Excited state electronic structures and dynamics of NOCl: a new potential function set, absorption spectrum, and photodissociation mechanism

A set of analytical potential energy surfaces (PESs) for six singlet excited states of NOCl are constructed based on multireference configuration interaction calculations. The total absorption cross section at the energy range of 2-7 eV is calculated by quantum dynamics calculations with the present PESs and transition dipole moments. The calculated absorption spectrum agrees well with the experiment. It is also found that the A band with the absorption maximum at 6.3 eV is attributed to the transition to the 4 1A' state, though the excitations to the 3 1A' and 3 1A" states contribute to the spectrum at the energy range between 4 and 5 eV. The spin-forbidden transitions are concluded to be negligibly weak. The mechanism of photodissociation reaction at the energy region corresponding to the A band is examined. The nonadiabatic transition rates from the 4 1A' state to lower singlet and triplet states are estimated by Fermi's golden rule, and the transitions to the 1 1A' and 3 1A' states induced by vibronic coupling are found to be the predominant dissociation pathways. The experimentally observed energy dependence of the recoil anisotropy of the fragments is discussed based on the calculated nonadiabatic transition rates.     

Symmetry-switching molecular Fe(O2)n(+) clusters

Experimental and theoretical studies based on mass spectrometry, collision-induced dissociation, and ab initio calculations are performed on the formation and stability of FeO(n)(+) clusters, as well as on their structural, electronic, and magnetic properties. In the mass spectra, clusters with an even number of oxygen atoms show increased stability, most prominently for FeO(10)(+). The extra stability of this cluster is confirmed by measurements of fragmentation cross sections through crossed molecular beam experiments. In addition, the calculations indicate a structural phase transition at this size, and most importantly, the FeO(n)(+) clusters show unique magnetic features, exhibiting isoenergetic low-spin (LS) and high-spin (HS) ground states. In the LS state, the magnetic moments of the O atoms adopt an antiferromagnetic alignment with respect to the magnetic moment of Fe(+), whereas in the HS state, the alignment is ferromagnetic. FeO(10)(+) is the largest thermodynamically stable complex, with the highest magnetic moment among the FeO(n)(+) clusters (13 μ(B) in HS).     

Scaled Born cross sections for excitations of H2 by electron impact

This article describes the scaling of plane-wave Born cross sections for the excitation of the H(2) molecule to four low-lying electronic states (B (1)Sigma(u) (+), C (1)Pi(u), B(') (1)Sigma(u) (+), and D (1)Pi(u)) by electron impact. The same BE and BEf scaling methods used on atoms were found to be equally effective for H(2) in converting Born cross sections into cross sections in good agreement with available experiments. These scaling methods are applicable only to dipole-allowed excitations. The possibility of using these scaling methods, as was done in atoms, to estimate the contribution of inner-shell excitations to the total ionization cross section via the excitation-autoionization mechanism is discussed, though this type of indirect ionization in molecules is not as common as in atoms.     

Ab initio investigation of the electronic structure and bonding of the HC(N2)x(+) and HC(CO)x(+) cations, x = 1, 2

Employing the coupled-cluster approach and correlation consistent basis sets of triple and quadruple cardinality, we have investigated the electronic structure and bonding of the HC(N2)x(+) and HC(CO)x(+), x = 1, 2, molecular cations. We report geometries, binding energies and potential energy profiles. The ground states of HC(N2)+, HC(CO)+ and HC(N2)2(+), HC(CO)2(+) are of 3sigma- and 1A1 symmetries, respectively. All four charged species are well bound with binding energies ranging from 81 [HC(N2)+ (X3sigma-) --> CH+(a3pi) + N2(X1sigma(g)+)] to 178 [HC(CO)2(+)(X1A1) --> CH+(X1sigma+) + 2CO(X1sigma+)] kcal/mol. It is our belief that the X1A1 states of HC(N2)2(+) and HC(CO)2(+) are isolable in the solid state if combined with appropriate counteranions.     

Ar(n)HF van der Waals clusters revisited: II. Energetics and HF vibrational frequency shifts from diffusion Monte Carlo calculations on additive and nonadditive potential-energy surfaces for n=1-12

The ground-state energies and HF vibrational frequency shifts of Ar(n)HF clusters have been calculated on the nonadditive potential-energy surfaces (PESs) for n=2-7 and on the pairwise-additive PESs for the clusters with n=1-12, using the diffusion Monte Carlo (DMC) method. For n>3, the calculations have been performed for the lowest-energy isomer and several higher-lying isomers which are the closest in energy. They provide information about the isomer dependence of the HF redshift, and enable direct comparison with the experimental data recently obtained in helium nanodroplets. The agreement between theory and experiment is excellent, in particular, for the nonadditive DMC redshifts. The relative, incremental redshifts are reproduced accurately even at the lower level of theory, i.e., the DMC and quantum five-dimensional (rigid Ar(n)) calculations on the pairwise-additive PESs. The nonadditive interactions make a significant contribution to the frequency shift, on the order of 10%-12%, and have to be included in the PESs in order for the theory to yield accurate magnitude of the HF redshift. The energy gaps between the DMC ground states of the cluster isomers are very different from the energy separation of their respective minima on the PES, due to the considerable variations in the intermolecular zero-point energy of different Ar(n)HF isomers.     

Three-state trajectory surface hopping studies of the photodissociation dynamics of formaldehyde on ab initio potential energy surfaces

Full-dimensional, three-state, surface hopping calculations of the photodissociation dynamics of formaldehyde are reported on ab initio potential energy surfaces (PESs) for electronic states S(1), T(1), and S(0). This is the first such study initiated on S(1) with ab initio-calculated spin-orbit couplings among the three states. We employ previous PESs for S(0) and T(1), and a new PES for S(1), which we describe here, as well as new spin-orbit couplings. The time-dependent electronic state populations and the branching ratio of radical products produced from S(0) and T(1) states and that of total radical products and molecular products at three total energies are calculated. Details of the surface hopping dynamics are described, and a novel pathway for isomerization on T(1) via S(0) is reported. Final translational energy distributions of H + HCO products from S(0) and T(1) are also reported as well as the translational energy distribution and final rovibrational distributions of H(2) products from the molecular channel. The present results are compared to previous trajectory calculations initiated from the global minimum of S(0). The roaming pathway leading to low rotational distribution of CO and high vibrational population of H(2) is observed in the present calculations.     

Ab initio study of the low lying electronic states of ZnF and ZnF-

Highly correlated ab initio calculations have been performed for an accurate determination of the electronic structure and of the spectroscopy of the low lying electronic states of the ZnF system. Using effective core pseudopotentials and aug-cc-pVQZ basis sets for both atoms, the potential curves, the dipole moment functions, and the transition dipole moments between relevant electronic states have been calculated at the multireference-configuration-interaction level. The spectroscopic constants calculated for the X(2)Sigma(+) ground state are in good agreement with the most recent theoretical and experimental values. It is shown that, besides the X(2)Sigma(+) ground state, the B(2)Sigma(+), the C(2)Pi, and the D(2)Sigma(+) states are bound. The A(2)Pi state, which has been mentioned in previous works, is not bound but its potential presents a shoulder in the Franck-Condon region of the X(2)Sigma(+) ground state. All of the low lying quartet states are found to be repulsive. The absorption transitions from the v=0 level of the X(2)Sigma(+) ground state toward the three bound states have been evaluated and the spectra are presented. The potential energy of the ZnF(-) molecular anion has been determined in the vicinity of its equilibrium geometry and the electronic affinity of ZnF (EA=1.843 eV with the zero energy point correction) has been calculated in agreement with the photoelectron spectroscopy experiments.     

Rotational excitations of N2O in small helium clusters and the role of Bose permutation symmetry

We present a detailed study of the energetics, structures, and Bose properties of small clusters of (4)He containing a single nitrous oxide (N(2)O) molecule, from N=1 (4)He up to sizes corresponding to completion of the first solvation shell around N(2)O (N=16 (4)He). Ground state properties are calculated using the importance-sampled rigid-body diffusion Monte Carlo method, rotational excited state calculations are made with the projection operator imaginary time spectral evolution method, and Bose permutation exchange and associated superfluid properties are calculated with the finite temperature path integral method. For N< or =5 the helium atoms are seen to form an equatorial ring around the molecular axis, at N=6 helium density starts to occupy the second (local) minimum of the N(2)O-He interaction at the oxygen side of the molecule, and N=9 is the critical size at which there is onset of helium solvation all along the molecular axis. For N> or =8 six (4)He atoms are distributed in a symmetric, quasirigid ring around N(2)O. Path integral calculations show essentially complete superfluid response to rotation about the molecular axis for N> or =5, and a rise of the perpendicular superfluid response from zero to appreciable values for N> or =8. Rotational excited states are computed for three values of the total angular momentum, J=1-3, and the energy levels fitted to obtain effective spectroscopic constants that show excellent agreement with the experimentally observed N dependence of the effective rotational constant B(eff). The non-monotonic behavior of the rotational constant is seen to be due to the onset of long (4)He permutation exchanges and associated perpendicular superfluid response of the clusters for N> or =8. We provide a detailed analysis of the role of the helium solvation structure and superfluid properties in determining the effective rotational constants.     

Path integral Monte Carlo study of 4He clusters doped with alkali and alkali-earth ions

Path integral Monte Carlo calculations of (4)He nanodroplets doped with alkali (Na(+), K(+) and Cs(+)) and alkali-earth (Be(+) and Mg(+)) ions are presented. We study the system at T = 1 K and between 14 and 128 (4)He atoms. For all studied systems, we find that the ion is well localized at the center of the droplet with the formation of a "snowball" of well-defined shells of localized (4)He atoms forming solid-like order in at least the first surrounding shell. The number of surrounding helium shells (two or three) and the number of atoms per shell and the degree of localization of the helium atoms are sensitive to the type of ion. The number of (4)He atoms in the first shell varies from 12 for Na(+) to 18 for Mg(+) and depends weakly on the size of the droplet. The study of the density profile and of the angular correlations shows that the local solid-like order is more pronounced for the alkali ions with Na(+) giving a very stable icosahedral order extending up to three shells.     

Microsolvation of LiH+ in helium clusters: many-body effects and additivity models for the interaction forces

The ab initio calculation of the interaction forces between the LiH+ molecular ion, at its equilibrium geometry, and several He atoms is carried out in order to isolate and assess the importance of many-body contributions in the search for realistic energy and geometry data. The full potential energy surface (PES) with a single helium partner is obtained first by using an aug-cc-pVQZ basis set for He and higher quality ones for Li and H. The calculations were performed at the CAS-SCF plus MRCI level for the lowest potential energy surface over a total of 480 grid points of the two intermolecular Jacobi coordinates, whereas the excited state surface has also been examined in order to exclude the presence of any significant nonadiabatic interaction between the two PESs. A numerical fit of the lower surface is presented and the general physical changes of the ionic interaction when going from the lower to the upper of the two potentials are described and discussed. The fairly limited importance of many-body effects for such systems is seen from further ab initio calculations including several He atoms: our results suggest that, at least in the present case, no strong charge migration occurs after He attachment, and therefore, one could realistically model larger clusters by implementing a sum-of-potentials approach via the presently computed PES.     

Avoided crossings, conical intersections, and low-lying excited states with a single reference method: The restricted active space spin-flip configuration interaction approach

The restricted active space spin-flip CI (RASCI-SF) performance is tested in the electronic structure computation of the ground and the lowest electronically excited states in the presence of near-degeneracies. The feasibility of the method is demonstrated by analyzing the avoided crossing between the ionic and neutral singlet states of LiF along the molecular dissociation. The two potential energy surfaces (PESs) are explored by means of the energies of computed adiabatic and approximated diabatic states, dipole moments, and natural orbital electronic occupancies of both states. The RASCI-SF methodology is also used to study the ground and first excited singlet surface crossing involved in the double bond isomerization of ethylene, as a model case. The two-dimensional PESs of the ground (S(0)) and excited (S(1)) states are calculated for the complete configuration space of torsion and pyramidalization molecular distortions. The parameters that define the state energetics in the vicinity of the S(0)/S(1) conical intersection region are compared to complete active space self-consistent field (CASSCF) results. These examples show that it is possible to describe strongly correlated electronic states using a single reference methodology without the need to expand the wavefunction to high levels of collective excitations. Finally, RASCI is also examined in the electronic structure characterization of the ground and 2(1)A(g) (-), 1(1)B(u) (+), 1(1)B(u) (-), and 1(3)B(u) (-) states of all-trans polyenes with two to seven double bonds and beyond. Transition energies are compared to configuration interaction singles, time-dependent density functional theory (TDDFT), CASSCF, and its second-order perturbation correction calculations, and to experimental data. The capability of RASCI-SF to describe the nature and properties of each electronic state is discussed in detail. This example is also used to expose the properties of different truncations of the RASCI wavefunction and to show the possibility to use an excitation operator with any number of α-to-β electronic promotions.     

Bosonic helium droplets with cationic impurities: onset of electrostriction and snowball effects from quantum calculations

Variational Monte Carlo and diffusion Monte Carlo calculations have been carried out for cations such as Li(+), Na(+), and K(+) as dopants of small helium clusters over a range of cluster sizes up to about 12 solvent atoms. The interaction has been modeled through a sum-of-potential picture that disregards higher order effects beyond atom-atom and atom-ion contributions. The latter were obtained from highly correlated ab initio calculations over a broad range of interatomic distances. This study focuses on two of the most striking features of the microsolvation in a quantum solvent of a cationic dopant: electrostriction and snowball effects. They are discussed here in detail and in relation with the nanoscopic properties of the interaction forces at play within a fully quantum picture of the cluster features.     

化学反应的高精度从头算势能面

势能面是化学反应动力学研究的基础。近年来随着理论方法的发展与计算技术的进步,不但含三、四个原子反应体系的电子基态势能面的构建精度进一步提高,一些反应体系的多电子态耦合势能面的构建和含六个原子以上反应体系的高维从头算势能面的构建也取得了重要进展。本文结合若干典型体系势能面的构建工作,主要介绍了高精度电子基态势能面,包括Renner-Teller、旋轨耦合等非绝热效应的耦合势能面以及高维势能面方面的研究进展。     

RgBF2(+) complexes (Rg = Ar, Kr, and Xe): the cations with large stabilities

Rare gas containing cations with general formula [Rg, B, 2F](+) have been investigated theoretically by second-order Mo?ller-Plesset perturbation, coupled cluster, and complete active space self-consistent field levels of theory with correlation-consistent basis sets. Totally two types of minima, i.e., boron centered C(2) (v) symmetried RgBF(2) (+) (Rg = Ar, Kr, and Xe) which can be viewed as loss of F(-) from FRgBF(2) and linear FRgBF(+) (Rg = Kr and Xe) are obtained at the CCSD(T)∕aug-cc-pVTZ∕SDD and CASSCF(10,8)∕aug-cc-pVTZ∕SDD levels, respectively. It is shown that the RgBF(2) (+) are global minima followed by FRgBF(+) at 170.9 and 142.2 kcal∕mol on the singlet potential-energy surfaces of [Rg, B, 2F](+) (Rg = Kr and Xe) at the CASPT2(10,8) ∕aug-cc-pVTZ∕SDD∕∕CASSCF(10,8)∕aug-cc-pVTZ∕SDD, respectively. The interconversion barrier heights between RgBF(2) (+) and FRgBF(+) (Rg = Kr and Xe) are at least 39 kcal∕mol. In addition, no dissociation transition state associated with RgBF(2) (+) and FRgBF(+) can be found. This suggests that RgBF(2) (+) (Rg = Ar, Kr, and Xe) can exist as both thermodynamically and kinetically stable species, while linear FRgBF(+) (Rg = Kr and Xe) can exist as metastable species compared with the lowest dissociation limit energies just like isoelectronic linear FRgBO and FRgBN(-). From natural bond orbital and atoms-in-molecules calculations, it is found that the positive charge is mainly located on Rg and boron atoms for both types of minima, the Rg-B bonds of ArBF(2) (+), KrBF(2) (+), and XeBF(2) (+) are mostly electrostatic, thus can be viewed as ion-induced dipole interaction; while that of linear FKrBF(+) and FXeBF(+) are covalent in nature. The previous experimental observation of ArBF(2) (+) by Pepi et al. [J. Phys. Chem. B. 110, 4492 (2006)] should correspond to C(2) (v) minimum. The presently predicted spectroscopies of KrBF(2) (+), XeBF(2) (+), FKrBF(+), and FXeBF(+) should be helpful for their experimental identification in the future.