Explicitly correlated Gaussian calculations of the 2P(o) Rydberg spectrum of the lithium atom

Accurate quantum-mechanical nonrelativistic variational calculations are performed for the nine lowest members of the (2)P(o) Rydberg series (1s(2)np(1), n = 2, ..., 10) of the lithium atom. The effect of the finite nuclear mass is included in the calculations allowing for determining the isotopic shifts of the energy levels. The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian functions. The exponential parameters of the Gaussians are variationally optimized with the aid of the analytical energy gradient determined with respect to those parameters. The calculated state energies are compared with the available experimental data.     

1s22p3 and 1s22s23l, l = s,p,d, excited states of boron isoelectronic series from explicitly correlated wave functions

For some members of the boron isoelectronic series and starting from explicitly correlated wave functions, six low-lying excited states have been studied. Three of them arise from the 1s(2)2p(3) configuration, and the other three from the 1s(2)2s(2)3l, l = s,p,d, configurations. This work follows a previous one on both the 1s(2)2s(2)2p-(2)P ground state and the four excited states coming from the 1s(2)2s2p(2) configuration. Energies, one- and two-body densities in position space and some other two-body properties in position and momentum spaces have been obtained. A systematic analysis of the energetic ordering of the states as a function of the total orbital angular momentum and spin is performed in terms of the electron-nucleus and electron-electron potential energies and the role of the angular correlation is discussed. All calculations have been carried out by using the Monte Carlo algorithm.     

Relativistic corrections to the non-Born-Oppenheimer energies of the lowest singlet Rydberg states of 3He and 4He

In this work the authors present an approach to calculate the leading-order relativistic corrections for ground and excited states of helium isotopomers. In the calculations they used variational wave functions expanded in terms of explicitly correlated Gaussians obtained without assuming the Born-Oppenheimer approximation.     

Accurate one-dimensional potential energy curve of the linear (H2)2 cluster

We present a sub-0.3 K accuracy, ground-state one-dimensional potential energy curve of the metastable linear configuration of the (H(2))(2) cluster calculated exclusively with explicitly correlated Gaussian functions with shifted centers. The H(2) internuclear distance is kept at the isolated H(2) vibrational ground-state average value of 1.448 736 bohr and the intermonomer separation is varied between 2 and 100 bohrs. The analytical gradient of the energy with respect to the nonlinear parameters of the Gaussians (i.e., the exponents and the coordinates of the shifts) has been employed in the variational optimization of the wave function. Procedures for enlarging the basis set and for adjusting the centers of the Gaussians to the varying intermonomer separation have been developed and used in the calculations.     

Non-Born-Oppenheimer variational calculation of the ground-state vibrational spectrum of LiH+

Very accurate, rigorous, variational, non-Born-Oppenheimer (non-BO) calculations have been performed for the fully symmetric, bound states of the LiH(+) ion. These states correspond to the ground and excited vibrational states of LiH(+) in the ground (2)Sigma(+) electronic state. The non-BO wave functions of the states have been expanded in terms of spherical N-particle explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance and 5600 Gaussians were used for each state. The calculations that, to our knowledge, are the most accurate ever performed for a diatomic system with three electrons have yielded six bound states. Average interparticle distances and nucleus-nucleus correlation function plots are presented.     

Refinement of the experimental energy levels of higher 2D Rydberg states of the lithium atom with very accurate quantum mechanical calculations

Very accurate variational non-relativistic calculations are performed for four higher Rydberg (2)D states (1s(2)nd(1), n = 8,..., 11) of the lithium atom ((7)Li). The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian functions and finite nuclear mass is used. The exponential parameters of the Gaussians are optimized using the variational method with the aid of the analytical energy gradient determined with respect to those parameters. The results of the calculations allow for refining the experimental energy levels determined with respect to the (2)S 1s(2)2s(1) ground state.     

Explicitly correlated Gaussian calculations of the (2)D Rydberg states of the boron atom

Accurate non-relativistic variational calculations are performed for the seven lowest members of the (2)D Rydberg series (1s(2)2s2p(2), and 1s(2)2s(2)nd, n = 3, [ellipsis (horizontal)], 8) of the boron atom. The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian basis functions and the effect of the finite nuclear mass is directly included in the calculations allowing for determining the isotopic shifts of the energy levels. The Gaussian basis is optimized independently for each state with the aid of the analytic energy gradient with respect to the Gaussian parameters. The calculations represent the highest accuracy level currently achievable for the considered states. The computed energies are compared with the available experimental data.     

Darwin and mass-velocity relativistic corrections in the non-Born-Oppenheimer calculations of pure vibrational states of H2

The Darwin and mass-velocity relativistic corrections have been calculated for all pure vibrational states of the H2 using the perturbation theory and very accurate variational wave functions obtained without assuming the Born-Oppenheimer (BO) approximation. Expansions in terms of explicitly correlated Gaussians with premultipliers in the form of even powers of the internuclear distance were used for the wave functions. With the inclusion of the two relativistic corrections to the non-BO energies the transition energies for the highest states agree more with the experimental results.     

The 1 (2)A(1), 1 (2)B(2), and 1 (2)A(2) states of the SO(2) (+) ion studied using multiconfiguration second-order perturbation theory

The 1 (2)A(1), 1 (2)B(2), and 1 (2)A(2) electronic states of the SO(2) (+) ion have been studied using multiconfiguration second-order perturbation theory (CASPT2) and two contracted atomic natural orbital basis sets, S[6s4p3d1f]/O[5s3p2d1f] (ANO-L) and S[4s3p2d]/O[3s2p1d] (ANO-S), and the three states were considered to correspond to the observed X, B, and A states, respectively, in the previous experimental and theoretical studies. Based on the CASPT2/ANO-L adiabatic excitation energy calculations, the X, A, and B states of SO(2) (+) are assigned to 1 (2)A(1), 1 (2)B(2), and 1 (2)A(2), respectively, and our assignments of the A and B states are contrary to the previous assignments (A to (2)A(2) and B to (2)B(2)). The CASPT2/ANO-L energetic calculations also indicate that the 1 (2)A(1), 1 (2)B(2), and 1 (2)A(2) states are, respectively, the ground, first excited, and second excited states at the ground-state (1 (2)A(1)) geometry of the ion and at the geometry of the ground-state SO(2) molecule. Based on the CASPT2/ANO-L results for the geometries, we realize that the experimental geometries (determined by assuming the bond lengths to be the same as the neutral ground state of SO(2)) were not accurate. The CASPT2/ANO-S calculations for the potential energy curves as functions of the OSO angle confirm that the 1 (2)B(2) and 1 (2)A(2) states are the results of the Renner-Teller effect in the degenerate (2)Pi(g) state at the linear geometry, and it is clearly shown that the 1 (2)B(2) curve, as the lower component of the Renner splitting, lies below the 1 (2)A(2) curve. The UB3LYP/cc-pVTZ adiabatic excitation energy calculations support the assignments (A to (2)B(2) and B to (2)A(2)) based on the CASPT2/ANO-L calculations.     

The enthalpies of formation of AsX(n) molecules, where X=H, F or Cl, and n=1, 2 or 3, by RCCSD(T) and UCCSD(T)-F12x calculations

RCCSD(T) and UCCSD(T)-F12x calculations were performed on AsX(n) molecules, where X = H, F or Cl, and n = 1, 2 or 3, and related species, in order to evaluate their enthalpies of formation (ΔH(f)(?)). The recommended ΔH(f)(?) values obtained from the present investigation are AsH, 57.7(2); AsF, -7.9(3); AsCl, 27.2(4); AsH(2), 39.8(4); AsF(2), -96.6(9); AsCl(2), -17.8(10); AsH(3), 17.1(4); AsF(3)-196.0(5) and AsCl(3), -59.1(27) kcal mole(-1). These values are anchored only on one thermodynamic quantity, namely, ΔH(f)(?)(As) (= 70.3 kcal mole(-1)). In the calculations, the fully-relativistic small-core effective core potential (ECP10MDF) was used for As. Contributions from outer core correlation of As 3d(10) electrons were computed explicitly in both RCCSD(T) and UCCSD(T)-F12 calculations with additional tight basis functions designed for As 3d(10) electrons. Basis sets of up to augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality were used in RCCSD(T) calculations and computed relative electronic energies were extrapolated to the complete basis set (CBS) limit. For the simplified, explicitly correlated UCCSD(T)-F12x calculations, basis sets of up to quadruple-zeta (QZ) quality were employed. Based on the RCCSD(T)/CBS benchmark values, the reliability of available theoretical and experimental values have been assessed.     

Theoretical investigation of ground and excited states of the methylene amidogene radical (H(2)CN)

The excited states and the absorption spectrum of the methylene amidogene radical are studied by high-level ab initio calculations. The multireference configuration interaction method was used in combination with different basis sets and basis set extrapolation to compute equilibrium geometries, harmonic frequencies, and excitation energies of the four lowest doublet electronic states of the title species. Potential curves and transition dipole moment functions were determined along the normal mode coordinates of the electronic ground state. These functions were employed to determine vibronic absorption spectra. The intensities of dipole forbidden but vibronically allowed transitions were calculated by explicitly evaluating integrals over the vibrational wave functions and the transition dipole functions of the involved electronic states. By this method the oscillator strengths of the dipole allowed (2)A(1)<--(2)B(2) and the dipole forbidden (2)B(1)<--(2)B(2) bands were computed. It turns out that the dipole forbidden transition is two orders of magnitude weaker than the dipole allowed one. The 0-0 excitation energies are found to be 30 256 cm(-1) for the (2)B(1) state and 34,646 cm(-1) for the (2)A(1) state. From the combined results of the excitation energies and oscillator strengths it is concluded that the experimentally observed peaks must be due to the (2)A(1) state, in contradiction to earlier assignments.     

Electron-impact excitation of the transition 1s 22s 22p 5(2 P 3 2/0 −2 P 1 2/0 ) of Mg IV

Effective collision strengths for electron-impact excitation of the ground state fine structure transition² P ₃ ^2/0 ?² P ₁ ^2/0 have been calculated by using theR-matrix method. Twelve lowest target states, represented by configuration interaction wavefunctions are included in the scattering calculations.M1 andE2 transition probabilities are also calculated by employing the Breit-Pauli Hamiltonian.     

Computational study of the rovibrational spectrum of (OCS)2

In this paper, we report a new intermolecular potential energy surface and rovibrational transition frequencies and line strengths computed for the OCS dimer. The potential is made by fitting energies obtained from explicitly correlated coupled-cluster calculations and fit using an interpolating moving least squares method. The rovibrational Schroedinger equation is solved with a symmetry-adapted Lanczos algorithm and an uncoupled product basis set. All four intermolecular coordinates are included in the calculation. On the potential energy surface we find, previously unknown, cross-shaped isomers and also polar and non-polar isomers. The associated wavefunctions and energy levels are presented. To identify polar and cross states we use both calculations of line strengths and vibrational parent analysis. Calculated rotational constants differ from their experimental counterparts by less than 0.001 cm(-1).     

HFS constants of Be I 1s22s2p 3p0, B I 1s22s2p24P and B I 1s22s2p22D obtained from the non-closed shell many-electron theory for excited states

Using the non-closed shell many-electron theory for excited states of Sinano?lu we compute the hyperfine structure constants of Be I 1s²2s2p ³P⁰, B I 1s²2s2p²⁴P and B I 1s²2s2p²²D. The quadrupole moment of ⁹Be is also newly determined to be, Q = 0.05494 barns.     

Accurate non-Born-Oppenheimer calculations of the complete pure vibrational spectrum of D2 with including relativistic corrections

In this work we report very accurate variational calculations of the complete pure vibrational spectrum of the D(2) molecule performed within the framework where the Born-Oppenheimer (BO) approximation is not assumed. After the elimination of the center-of-mass motion, D(2) becomes a three-particle problem in this framework. As the considered states correspond to the zero total angular momentum, their wave functions are expanded in terms of all-particle, one-center, spherically symmetric explicitly correlated Gaussian functions multiplied by even non-negative powers of the internuclear distance. The nonrelativistic energies of the states obtained in the non-BO calculations are corrected for the relativistic effects of the order of α(2) (where α = 1/c is the fine structure constant) calculated as expectation values of the operators representing these effects.     

Calculations of the ground states of BeH and BeH+ without the Born-Oppenheimer approximation

Non-Born-Oppenheimer variational calculations employing explicitly correlated Gaussian basis functions have been performed for the ground states of the beryllium monohydride molecule (BeH) and its ion (BeH+), as well as for the beryllium atom (Be) and its ion (Be+). An approach based on the analytical energy gradient calculated with respect to the Gaussian exponential parameters was employed. The calculated energies were used to determine the ionization potential of BeH and the dissociation energies of BeH and BeH+. Also, the generated wave functions were used to compute various expectation values, such as the average interparticle distances and the nucleus-nucleus correlation functions.     

Accurate theoretical study of PS(q) (q = 0,+1,-1) in the gas phase

Highly correlated ab initio methods were used in order to generate the potential energy curves and spin-orbit couplings of electronic ground and excited states of PS and PS(+). We also computed those of the bound parts of the electronic states of the PS(-) anion. We used standard coupled cluster CCSD(T) level with augmented correlation-consistent basis sets, internally contacted multi-reference configuration interaction, and the newly developed CCSD(T)-F12 methods in connection with the explicitly correlated basis sets. Core-valence correction and scalar relativistic effects were examined. Our data consist of a set of spectroscopic parameters (equilibrium geometries, harmonic vibrational frequencies, rotational constants, spin-orbit, and spin-spin constants), adiabatic ionization energies, and electron affinities. For the low laying electronic states, our calculations are consistent with previous works whereas the high excited states present rather different shapes. Based on these new computations, the earlier ultraviolet bands of PS and PS(+) were reassigned. For PS(-) and in addition to the already known anionic three bound electronic states (i.e., X(3)Σ(-), (1)Δ, and 1(1)Σ(+)), our calculations show that the (1)Σ(-), (3)Σ(+), and the (3)Δ states are energetically below their quartet parent neutral state (a(4)Π). The depletion of the J = 3 component of PS(-)((3)Δ) will mainly occur via weak interactions with the electron continuum wave.     

Comparison of density functionals for energy and structural differences between the high- [5T2g: (t2g)4(eg)2] and low- [1A1g: (t2g)6(eg)0] spin states of the hexaquoferrous cation [Fe(H2O)6]2+

A comparison of density functionals is made for the calculation of energy and geometry differences for the high- [(5)T(2g): (t(2g))(4)(e(g))(2)] and low- [(1)A(1g): (t(2g))(6)(e(g))(0)] spin states of the hexaquoferrous cation [Fe(H(2)O)(6)](2+). Since very little experimental results are available (except for crystal structures involving the cation in its high-spin state), the primary comparison is with our own complete active-space self-consistent field (CASSCF), second-order perturbation theory-corrected complete active-space self-consistent field (CASPT2), and spectroscopy-oriented configuration interaction (SORCI) calculations. We find that generalized gradient approximations (GGAs) and the B3LYP hybrid functional provide geometries in good agreement with experiment and with our CASSCF calculations provided sufficiently extended basis sets are used (i.e., polarization functions on the iron and polarization and diffuse functions on the water molecules). In contrast, CASPT2 calculations of the low-spin-high-spin energy difference DeltaE(LH)=E(LS)-E(HS) appear to be significantly overestimated due to basis set limitations in the sense that the energy difference of the atomic asymptotes ((5)D-->(1)I excitation of Fe(2+)) are overestimated by about 3000 cm(-1). An empirical shift of the molecular DeltaE(LH) based upon atomic calculations provides a best estimate of 12 000-13 000 cm(-1). Our unshifted SORCI result is 13 300 cm(-1), consistent with previous comparisons between SORCI and experimental excitation energies which suggest that no such empirical shift is needed in conjunction with this method. In contrast, after estimation of incomplete basis set effects, GGAs with one exception underestimate this value by 3000-4000 cm(-1) while the B3LYP functional underestimates it by only about 1000 cm(-1). The exception is the GGA functional RPBE which appears to perform as well as or better than the B3LYP functional for the properties studied here. In order to obtain a best estimate of the molecular DeltaE(LH) within the context of density functional theory (DFT) calculations we have also performed atomic excitation energy calculations using the multiplet sum method. These atomic DFT calculations suggest that no empirical correction is needed for the DFT calculations.     

Statistical thermodynamics of 1-butanol, 2-methyl-1-propanol, and butanal

The purpose of the present investigation is to calculate partition functions and thermodynamic quantities, viz., entropy, enthalpy, heat capacity, and Gibbs free energies, for 1-butanol, 2-methyl-1-propanol, and butanal in the vapor phase. We employed the multi-structural (MS) anharmonicity method and electronic structure calculations including both explicitly correlated coupled cluster theory and density functional theory. The calculations are performed using all structures for each molecule and employing both the local harmonic approximation (MS-LH) and the inclusion of torsional anharmonicity (MS-T). The results obtained from the MS-T calculations are in excellent agreement with experimental data taken from the Thermodynamics Research Center data series and the CRC Handbook of Chemistry and Physics, where available. They are also compared with Benson's empirical group additivity values, where available; in most cases, the present results are more accurate than the group additivity values. In other cases, where experimental data (but not group additivity values) are available, we also obtain good agreement with experiment. This validates the accuracy of the electronic structure calculations when combined with the MS-T method for estimating the thermodynamic properties of systems with multiple torsions, and it increases our confidence in the predictions made with this method for molecules and temperatures where experimental or empirical data are not available.