Theoretical spectroscopy of the N2HAr+ complex

The six-dimensional potential energy surface of the electronic ground state of N2HAr+ is determined by ab initio computations at the CCSD(T) level of theory. The potential energy surface is used to derive a set of spectroscopic data for N2HAr+ and N2DAr+ using second order perturbation theory. Full six-dimensional (6-D) rotation-vibration computations are also carried out using an analytical representation of the surface for J=0 and 1, in order to deduce the rovibrational spectra of N2HAr+ and its deuterated isotopomer. Our variationally determined anharmonic wavenumbers differ by less than 15 cm(-1) from the most accurate experimental values. Strong anharmonic resonances are found between the rovibrational levels of both cations even at low energies.     

Ab initio potential energy surface and bound states for the Kr-OCS complex

The first ab initio potential energy surface of the Kr-OCS complex is developed using the coupled-cluster singles and doubles with noniterative inclusion of connected triples [CCSD(T)]. The mixed basis sets, aug-cc-pVTZ for the O, C, and S atom, and aug-cc-pVQZ-PP for the Kr atom, with an additional (3s3p2d1f) set of midbond functions are used. A potential model is represented by an analytical function whose parameters are fitted numerically to the single point energies computed at 228 configurations. The potential has a T-shaped global minimum and a local linear minimum. The global minimum occurs at R = 7.146 a(0), θ = 105.0° with energy of -270.73 cm(-1). Bound state energies up to J = 9 are calculated for three isotopomers (82)Kr-OCS, (84)Kr-OCS, and (86)Kr-OCS. Analysis of the vibrational wavefunctions and energies suggests the complex can exist in two isomeric forms: T-shaped and quasi-linear. The calculated transition frequencies and spectroscopic constants of the three isotopomers are in good agreement with the experimental values.     

Ab initio potential energy surface and bound states of the Xe-CO complex

The first two-dimensional potential energy surface for the Xe-CO van der Waals interaction is calculated by the single and double excitation coupled-cluster theory with noniterative treatment of triple excitations. Mixed basis sets, aug-cc-pVQZ for the C and O atoms, and aug-cc-pVQZ-PP for the Xe atom, with an additional (3s3p2d2f1g) set of midbond functions, are used. Our potential energy surface has a single, nearly T-shaped minimum of -131.87 cm(-1) at R(e)=7.80a(0) and theta(e)=102.5 degrees. Based on the potential, the bound state energies are calculated for seven isotopomers of the Xe-(12)C(16)O complex, seven isotopomers of the Xe-(13)C(16)O complex, and three isotopomers of the Xe-(13)C(18)O complex. Compared with available experimental data, the predicted transition frequencies and spectroscopic constants are in good agreement with the experimental results.     

Rovibrational spectra of ammonia. I. Unprecedented accuracy of a potential energy surface used with nonadiabatic corrections

In this work, we build upon our previous work on the theoretical spectroscopy of ammonia, NH(3). Compared to our 2008 study, we include more physics in our rovibrational calculations and more experimental data in the refinement procedure, and these enable us to produce a potential energy surface (PES) of unprecedented accuracy. We call this the HSL-2 PES. The additional physics we include is a second-order correction for the breakdown of the Born-Oppenheimer approximation, and we find it to be critical for improved results. By including experimental data for higher rotational levels in the refinement procedure, we were able to greatly reduce our systematic errors for the rotational dependence of our predictions. These additions together lead to a significantly improved total angular momentum (J) dependence in our computed rovibrational energies. The root-mean-square error between our predictions using the HSL-2 PES and the reliable energy levels from the HITRAN database for J = 0-6 and J = 7∕8 for (14)NH(3) is only 0.015 cm(-1) and 0.020∕0.023 cm(-1), respectively. The root-mean-square errors for the characteristic inversion splittings are approximately 1∕3 smaller than those for energy levels. The root-mean-square error for the 6002 J = 0-8 transition energies is 0.020 cm(-1). Overall, for J = 0-8, the spectroscopic data computed with HSL-2 is roughly an order of magnitude more accurate relative to our previous best ammonia PES (denoted HSL-1). These impressive numbers are eclipsed only by the root-mean-square error between our predictions for purely rotational transition energies of (15)NH(3) and the highly accurate Cologne database (CDMS): 0.00034 cm(-1) (10 MHz), in other words, 2 orders of magnitude smaller. In addition, we identify a deficiency in the (15)NH(3) energy levels determined from a model of the experimental data.     

Solvation of copper ions by acetone. structures and sequential binding energies of Cu+(acetone)x,x = 1-4 from collision-induced dissociation and theoretical studies

Collision-induced dissociation of Cu+(acetone)(x), x = 1-4, with Xe is studied as a function of kinetic energy using guided ion beam mass spectrometry. In all cases, the primary and lowest energy dissociation channel observed is endothermic loss of one acetone molecule. The primary cross section thresholds are interpreted to yield 0 and 298 K bond energies after accounting for the effects of multiple ion-neutral collisions, internal energy of the complexes, and dissociation lifetimes. Density functional calculations at the B3LYP/6-31G* level of theory are used to determine the structures of these complexes and provide molecular constants necessary for the thermodynamic analysis of the experimental data. Theoretical bond dissociation energies are determined from single point calculations at the B3LYP/6-311+G(2d,2p) and MP2(full)/6-311+G(2d,2p) levels, using the B3LYP/6-31G* optimized geometries. The experimental bond energies determined here are in good agreement with previous experimental measurements made in a high-pressure mass spectrometer for the sum of the first and second bond energy (i.e., Cu+(acetone)2 --> Cu+ + 2 acetone) when these results are properly anchored. The agreement between theory and experiment is reasonable in all cases, but varies both with the size of the cluster and the level of theory employed. B3LYP does an excellent job for the x = 1 and 3 clusters, but is systematically low for the x = 2 and 4 clusters such that the overall trends in sequential binding energies are not parallel. In contrast, all MP2 values are somewhat low, but the overall trends parallel the measured values for all clusters. The trends in the measured Cu+(acetone), binding energies are explained in terms of 4s-3d sigma hybridization effects and ligand-ligand repulsion in the clusters.     

Dual-level direct dynamics studies for the reactions of dimethyl ether with hydrogen atom and methyl radical

The dual-level direct dynamics approach is employed to study the dynamics of the CH(3)OCH(3) + H (R1) and CH(3)OCH(3) + CH(3) (R2) reactions. Low-level calculations of the potential energy surface are carried out at the MP2/6-311+G(d,p) level of theory. High-level energetic information is obtained at the QCISD(T) level of theory with the 6-311+G(3df,3pd) basis set. The dynamics calculations are performed using variational transition state theory (VTST) with the interpolated single-point energies (ISPE) method, and small-curvature tunneling (SCT) is included. It is shown that the reaction of CH(3)OCH(3) with H (R1) may proceed much easier and with a lower barrier height than the reaction with CH(3) radical (R2). The calculated rate constants and activation energies are in good agreement with the experimental values. The calculated rate constants are fitted to k(R1) = 1.16 x 10(-19) T(3) exp(-1922/T) and k(R2) = 1.66 x 10(-28) T(5) exp(-3086/T) cm(3) mol(-1) s(-1) over a temperature range 207-2100 K. Furthermore, a small variational effect and large tunneling effect in the lower temperature range are found for the two reactions.     

Density functional theory study of the Jahn-Teller effect and spin-orbit coupling for copper and gold trimers

The Born-Oppenheimer potential energy hypersurfaces of copper and gold trimers were calculated using density functional theory with an analytic potential. The calculated Jahn-Teller distortion energies, pseudorotation barriers, dissociation, and isomerization energies for the two trimers are discussed. Global minima from the surfaces were optimized using the density functional theory method as well as the coupled cluster-singles-doubles-with-triples energies technique. The agreement of the optimized structures with the analytic potential was very good. The Mulliken population analysis compared favorably with the experimental electron spin resonance results. Spin-orbit coupling was subsequently included and the effect was significant for gold, but negligible for copper. The spin-orbit effect suppressed the Jahn-Teller distortion of the gold trimer, and the potential surface with the spin-orbit effect included was also obtained. The spin-orbit splitting for the D(3h) geometry of the gold trimer was in excellent agreement with the most recent infrared spectroscopic results.     

Relativistic calculations of ground and excited states of LiYb molecule for ultracold photoassociation spectroscopy studies

We report a series of quantum-chemical calculations for the ground and some of the low-lying excited states of an isolated LiYb molecule by the spin-orbit multistate complete active space second-order perturbation theory (SO-MS-CASPT2). Potential energy curves, spectroscopic constants, and transition dipole moments (TDMs) at both spin-free and spin-orbit levels are obtained. Large spin-orbit effects especially in the TDMs of the molecular states dissociating to Yb((3)P(0,1,2)) excited states are found. To ensure the reliability of our calculations, we test five types of incremental basis sets and study their effect on the equilibrium distance and dissociation energy of the ground state. We also compare CASPT2 and CCSD(T) results for the ground state spectroscopic constants at the spin-free relativistic level. The discrepancies between the CASPT2 and CCSD(T) results are only 0.01 ? in equilibrium bond distance (R(e)) and 200 cm(-1) in dissociation energy (D(e)). Our CASPT2 calculation in the supermolecular state (R=100 a.u.) with the largest basis set reproduces experimental atomic excitation energies within 3% error. Transition dipole moments of the super molecular state (R=100 a.u.) dissociating to Li((2)P) excited states are quite close to experimental atomic TDMs as compared to the Yb((3)P) and Yb((1)P) excited states. The information obtained from this work would be useful for ultracold photoassociation experiments on LiYb.     

A new "spectroscopic" potential energy surface for formaldehyde in its ground electronic state

We report a new "spectroscopic" potential energy surface (PES) of formaldehyde (H(2)(12)C(16)O) in its ground electronic state, obtained by refining an ab initio PES in a least-squares fitting to the experimental spectroscopic data for formaldehyde currently available in the literature. The ab initio PES was computed using the CCSD(T)/aug-cc-pVQZ method at 30 840 geometries that cover the energy range up to 44 000 cm(-1) above equilibrium. Ro-vibrational energies of formaldehyde were determined variationally for this ab initio PES by means of the program TROVE [Theoretical ROtation-Vibration Energies; S. N. Yurchenko, W. Thiel, and P. Jensen, J. Mol. Spectrosc. 245, 126 (2007)]. The parameter values in the analytical representation of the PES were optimized in fittings to 319 ro-vibrational energies with J = 0, 1, 2, and 5. The initial parameter values in the fittings were those of the ab initio PES, the ro-vibrational eigenfunctions obtained from this PES served as a basis set during the fitting process, and constraints were imposed to ensure that the refined PES does not deviate unphysically from the ab initio one in regions of configuration space not sampled by the experimental data. The resulting refined PES, referred to as H(2)CO-2011, reproduces the available experimental J ≤ 5 data with a root-mean-square error of 0.04 cm(-1).     

Accurate ab initio ro-vibronic spectroscopy of the X̃2Π CCN radical using explicitly correlated methods

Explicitly correlated CCSD(T)-F12b calculations have been carried out with systematic sequences of correlation consistent basis sets to determine accurate near-equilibrium potential energy surfaces for the X(2)Π and a(4)Σ(-) electronic states of the CCN radical. After including contributions due to core correlation, scalar relativity, and higher order electron correlation effects, the latter utilizing large-scale multireference configuration interaction calculations, the resulting surfaces were employed in variational calculations of the ro-vibronic spectra. These calculations also included the use of accurate spin-orbit and dipole moment matrix elements. The resulting ro-vibronic transition energies, including the Renner-Teller sub-bands involving the bending mode, agree with the available experimental data to within 3 cm(-1) in all cases. Full sets of spectroscopic constants are reported using the usual second-order perturbation theory expressions. Integrated absorption intensities are given for a number of selected vibronic band origins. A computational procedure similar to that used in the determination of the potential energy functions was also utilized to predict the formation enthalpy of CCN, ΔH(f)(0K) = 161.7 ± 0.5 kcal/mol.     

Dissociation energies of six NO2 isotopologues by laser induced fluorescence spectroscopy and zero point energy of some triatomic molecules

We have measured the rotationless photodissociation threshold of six isotopologues of NO2 containing 14N, 15N, 16O, and 18O isotopes using laser induced fluorescence detection and jet cooled NO2 (to avoid rotational congestion). For each isotopologue, the spectrum is very dense below the dissociation energy while fluorescence disappears abruptly above it. The six dissociation energies ranged from 25 128.56 cm(-1) for 14N16O2 to 25 171.80 cm(-1) for 15N18O2. The zero point energy for the NO2 isotopologues was determined from experimental vibrational energies, application of the Dunham expansion, and from canonical perturbation theory using several potential energy surfaces. Using the experimentally determined dissociation energies and the calculated zero point energies of the parent NO2 isotopologue and of the NO product(s) we determined that there is a common De = 26 051.17+/-0.70 cm(-1) using the Born-Oppenheimer approximation. The canonical perturbation theory was then used to calculate the zero point energy of all stable isotopologues of SO2, CO2, and O3, which are compared with previous determinations.     

Theoretical study of the inner‐sphere energy barrier of the transition‐metal complex M(H2O) in electron‐transfer process

Two theoretical models, a reorganization model and an activation model, are presented for accurately determining the energy barrier of the type M(H₂O) of the transition‐metal complexes in the electron‐transfer process. Ab initio calculations are carried out at UMP2/6‐311G level for several redox pairs M(H₂O) (M=V, Cr, Mn, Fe, and Co) to calculate their inner‐sphere reorganization energies and activation energies according to the models presented in this article. The values of theoretical inner‐sphere reorganization energies and activational energies are comparable with the experimental results obtained from the vibration spectroscopic data. The theoretical reorganization energy of the every redox pair is four times as much as its activation energy, which agrees with Marcus' electron‐transfer theory. The fact proved that the theoretical models presented in this article are scientific and available for studying the electron‐transfer process of the transition‐metal complex. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 119–126, 1999     

Three-dimensional Potential Energy Surface and Bound States of the Ar2-Ne Complex

The first three-dimensional interaction potential energy surface (PES) of the Ar₂-Ne complex is developed using the single and double excitation coupled cluster theory with noniterative treatment of triple excitations CCSD(T). The aug-cc-pVQZ basis sets are employed for all atoms, including an additional (3s3p2d2f1g) set of midpoint bond functions. The calculated single point energies are fitted to an analytic two-dimensional potential model at each of seven fixed r_Ar2 values. The seven model potentials are then used to construct the three-dimensional PES by interpolating along (r—r_e) using a sixth-order polynomial. The PES is used in the following rovibrational energy levels calculations. The comparisons of theoretical transition frequencies and spectroscopic constants with the experimental results are given.     

A fully ab initio potential curve of near-spectroscopic quality for OH- ion: importance of connected quadruple excitations and scalar relativistic effects

A benchmark study has been carried out on the ground-state potential curve of the hydroxyl anion, OH-, including detailed calibration of both the l-particle and n-particle basis sets. The CCSD(T) basis set limit overestimates omega(e) by about 10 cm(-1), which is only remedied by inclusion of connected quadruple excitations in the coupled cluster expansion--or, equivalently, the inclusion of the 2pi orbitals in the active space of a multireference calculation. Upon inclusion of scalar relativistic effects (-3 cm(-1) on omega(e)), a potential curve of spectroscopic quality (sub-cm(-1) accuracy) is obtained. Our best computed EA(OH), 1.828 eV, agrees to three decimal places with the best available experimental value. Our best computed dissociation energies, D0(OH-) = 4.7796 eV and D0(OH) = 4.4124 eV, suggest that the experimental D0(OH) = 4.392 eV may possibly be about 0.02 eV too low.     

A simple variational-numerical approach to the ro-vibrational spectrum of diatomic molecules. An application to CH+

We describe a simple way of obtaining numerically the manifold of energies for ro-vibrational transitions for a centrifugally distorted oscillator, starting from the potential energy of the non-rotating oscillator calculated by an accurate ab initio method. It is shown that the energies so obtained compare well with those obtained variationally. The species of astrophysical interest methylidyne ion, CH(+), has been selected as an example that allow us to show the computational efficiency of the method with respect to the variational one. It is applied for the determination of ro-vibrational levels up J=6, and the spectroscopic parameters corresponding to the ground electronic state X(1)Sigma(+). From the potential energy surface computed at the MRCI/cc-pV5Z level, the fundamental frequency, B(0) and D(0) are determined to be 2724.8, 13.85688 and -1.53322x10(-3)cm(-1), respectively. We provide also an estimation of anharmonic constants.     

Joint photoelectron and theoretical study of (Rh(m)Co(n))⁻ (m=1-5, n=1-2) cluster anions and their neutral counterparts

Anion photoelectron spectroscopic experiments and calculations based on density functional theory have been used to investigate and uniquely identify the structural, electronic, and magnetic properties of both neutral and anionic (Rh(m)Co(n)) and (Rh(m)Co(n))(-) (m=1-5, n=1-2) clusters, respectively. Negative ion photoelectron spectra are presented for electron binding energies up to 3.493 eV. The calculated electron affinities and vertical detachment energies are in good agreement with the measured values. Computational results for geometric structures and magnetic moments of both cluster anions and their neutrals are presented.     

Energetics and mechanisms of C-H bond activation by a doubly charged metal ion: guided ion beam and theoretical studies of Ta2+ + CH4

A guided ion beam tandem mass spectrometer is used to study the kinetic-energy dependence of doubly charged atomic tantalum cations (Ta(2+)) reacting with CH4 and CD4. As for the analogous singly charged system, the dehydrogenation reaction to form TaCH2(2+) + H2 is exothermic. The charge-transfer reaction to form Ta(+) + CH4(+) and the charge-separation reaction to form TaH(+) + CH3(+) are also observed at low energies in exothermic processes, as is a secondary reaction of TaCH2(2+) to form TaCH3(+) + CH3(+). At higher energies, other doubly charged products, TaC(2+) and TaCH3(2+), are observed, although no formation of TaH(2+) was observed. Modeling of the endothermic cross sections provides 0 K bond dissociation energies (in electronvolts) of D0(Ta(2+)-C) = 5.42 +/- 0.19 and D0(Ta(2+)-CH3) = 3.40 +/- 0.16. These experimental bond energies are in poor agreement with density functional calculations at the B3LYP/HW+/6-311++G(3df,3p) level of theory. However, the Ta(2+)-C bond energy is in good agreement with calculations at the QCISD(T) level of theory, and the Ta(2+)-CH3 bond energy is in good agreement with density functional calculations at the BHLYP level of theory. Theoretical calculations reveal the geometric and electronic structures of all product ions and are used to map the potential energy surface, which describes the mechanism of the reaction and key intermediates. Both experimental and theoretical results suggest that TaH(+), TaCH2(2+), and TaCH3(2+) are formed through a H-Ta(2+)-CH3 intermediate.     

DFT and TD-DFT calculations on the electronic structures and spectroscopic properties of cyclometalated platinum(II) complexes

The electronic structures and spectroscopic properties of the three tridentate cyclometalated Pt(II) complexes Pt(N/\N/\C)C(triple bond)CPh (N/\N/\CH = 6-phenyl-2,2'-bipyridine) (1), Pt(N/\N/\S)C(triple bond)CPh (N/\N/\SH = 6-thienyl-2,2'-bipyridine) (2), and Pt(N/\N/\O)C(triple bond)CPh (N/\N/\OH = 6-furyl-2,2'-bipyridine) (3) were investigated theoretically using the density functional theory (DFT) method. The geometric structures of the complexes in the ground and excited states were explored at the B3LYP and UB3LYP levels, respectively. The absorption and emission spectra of the complexes in CH2Cl2 and CH3CN solutions were calculated by time-dependent density functional theory (TD-DFT) with the PCM solvent model. The calculated energies of the lowest singlet state and lowest triplet state in the three complexes are in good agreement with the results of experimental absorption and luminescence studies. All of the lowest-lying transitions were categorized as LLCT combined with MLCT transitions. The 623-nm emission of 1 from the 3A' --> 1A' transition was assigned as 3LLCT and 3MLCT transitions, whereas the 657- and 681-nm emissions of 2 and 3, respectively, were attributed to 3ILCT perturbed by 3MLCT transitions. NLO response calculations revealed that the nonzero values of the static first hyperpolarizability (beta0) for 1-3 are greatly enhanced through the introduction of the metal Pt(II) into the cyclometalated ligands, an effect that is determined by MLCT and LLCT transitions.