An ab initio multireference doubles excitation configuration interaction study of low-lying electronic states of Cd2 using Slater-type orbitals

Potential-energy curves for the ground state and lower excited states of the Cd₂ dimer have been calculated. They are obtained using a multireference doubles excitation configuration interaction procedure and employing Slater basis sets, previously optimized at the self-consistent-field level for excited states of the Cd atom. The spectroscopic constants and excitation energies for the bound states of Cd₂ have been compared with experimental as well as other theoretical results. The ground state of Cd₂ is essentially repulsive and presents a shallow van der Waals minimum. The computed adiabatic electronic transitions are in good agreement with the experimental ones. Received: 16 September 1999 / Accepted: 3 February 2000 / Published online: 2 May 2000     

CEPA calculations on open-shell molecules. III. Potential curves for the six lowest excited states of He2 in the vicinity of their equilibrium distances

CEPA-PNO and PNO-CI calculations have been performed for the potential energy curves of the He ₂ ⁺ ground state and the six lowest excited states of He₂ in the range of 1.4 a₀ ≤R ≤ 3.5 a₀. The calculated equilibrium distances as well as the spectroscopic constants are in very good agreement with molecular constants as derived experimentally from the rotation-vibration spectrum of He₂ by Ginter, except for thec ³∑ _g ⁺ state. This latter discrepancy is probably due to an “obligatory” hump in thec ³∑ _g ⁺ state occurring at 3.5 a₀ which cannot be properly treated in our calculation. The relative energetic positions of the six lowest states and their ionization energies are reproduced by our calculations with an accuracy of 0–400 cm⁻¹. Extrapolation of our results to infinite basis sets leads to estimates of the dissociation energies of He₂ excited states which cannot be measured spectroscopically because of the humps in all these states.     

Theoretical calculation of the low laying electronic states of the molecular ion RbH+

Using an ab initio method, the potential energy has been calculated for the 29 lowest molecular states of symmetries ²Σ⁺, ²Π, ²Δ for the molecular ion RbH⁺. The calculation is based on nonempirical pseudopotentials and parameterized ?‐dependent polarization potentials. Gaussian basis sets have been used for both atoms. The spectroscopic constants for 18 electronic sates have been calculated by fitting the calculated energy values to a polynomial in terms of the internuclear distance R. Through the canonical functions approach the eigenvalue E_v, the abscissas of the corresponding turning points (R_min and R_max) and the rotational constants B_v have been calculated up to 24 vibrational levels for the considered bound states. The comparison of the present results with those available in literature shows a very good agreement. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Quantum-mechanical treatment of molecules. I. Calculations of the potential energy curves and molecular constants of the ground and ionized states of N2

The self-consistent-field molecular-orbital method in LCAO (linear combination of atomic orbitals) approximation is applied to the ground and three ionized states of N₂ at a number of internuclear distances for the computation of the potential energy curves. In these calculations both the linear coefficients and the screening constants of the atomic orbitals have been optimized. The molecular constants ω_e, ω_ex_e, B_e, α_e, and R_e have also been calculated for the above states from the computed potential energy curves. The computed spectral results are compared with the experimental data as well as with the results reported by others from ab initio calculations.     

Theoretical calculation of the low laying electronic states of the molecular ion KH

Using an ab initio method the potential energy has been calculated for the 25 lowest molecular states of symmetries ²Σ⁺, ²Π, ²Δ for the molecular ion KH⁺. The calculation is based on nonempirical pseudopotentials and parameterized ℓ-dependent polarization potentials. Gaussian basis sets have been used for both atoms. The spectroscopic constants for 18 electronic sates have been calculated by fitting the calculated energy values to a polynomial in terms of the internuclear distance R. Through the canonical functions approach the eigenvalue E_v, the abscissas of the corresponding turning points (R_min and R_max) and the rotational constants B_v have been calculated up to 24 vibrational levels for the considered bound states. The comparison of the present results with those available in literature shows a very good agreement.     

Molecular orbital theory of the hydrogen bond. 25. Water–uracil complexes in excited n → π states

Hydrogen bonding of uracil with water in excited n → π* states has been investigated by means of ab initio SCF -CI calculations on uracil and water–uracil complexes. Two low-energy excited states arise from n → π* transitions in uracil. The first is due to excitation of the C₄? O group, while the second is associated with excitation of the C₂? O group. In the first n → π* state, hydrogen bonds at O₄ are broken, so that the open water–uracil dimer at O₄ dissociates. The “wobble” dimer, in which a water molecule is essentially free to move between its position in an open structure at N₃? H and a cyclic structure at N₃? H and O₄ in the ground state, collapses to a different “wobble” dimer at N₃? H and O₂ in the excited state. The third dimer, a “wobble” dimer at N₁? H and O₂, remains intact, but is destabilized relative to the ground state. Although hydrogen bonds at O₂ are broken in the second n → π* state, the three water–uracil dimers remain bound. The “wobble” dimer at N₁? H and O₂ changes to an excited open dimer at N₁? H. The “wobble” dimer at N₃? H and O₄ remains intact, and the open dimer at O₄ is further stabilized upon excitation. Dimer blue shifts of n → π* bands are nearly additive in 2:1 and 3:1 water:uracil structures. The fates of the three 2:1 water:uracil trimers and the 3:1 water:uracil tetramer in the first and second n → π* states are determined by the fates of the corresponding excited dimers in these states.     

Ab initio study of low-lying electronic states of the PF2 radical

The equilibrium geometries, excitation energies, force constants, and vibrational frequencies of the low-lying electronic states X²B₁, ²A₁, ²B₂, and ²A₂ of the PF₂ radical have been calculated at the MRSDCI level with a double zeta plus polarization basis set. Our calculated geometry, force constants, and vibrational frequencies for the X²B₁ state are in good agreement with experimental data. The electronic transition moments, oscillator strengths for the ²A₁ → X²B₁ and ²A₂ → X²B₁ transitions, and radiative lifetimes for the ²A₁ and ²A₂ states are calculated based on the MRSDCI wave functions. © 1994 by John Wiley & Sons, Inc.     

Ab initio calculations for the potential curves and spin–orbit coupling of Mg2

 The ground state and several low-lying excited states of the Mg₂ dimer have been studied by means of a combination of the complete-active-space multiconfiguration self-consistent-field (CASSCF)/CAS multireference second-order perturbation theory (CASPT2) method and coupled-cluster with single and double excitations and perturbative contribution of connected triple excitations [CCSD(T)] scheme. Reasonably good agreement with experiment has been obtained for the CCSD(T) ground-state potential curve but the dissociation energy of the only experimentally known A¹Σ _u ⁺ excited state of Mg₂ is somewhat overestimated at the CASSCF/CASPT2 level. The spectroscopic constants D _e, R _e and ω_e deduced from the calculated potential curves for other states are also reported. In addition, some spin–orbit matrix elements between the excited singlet and triplet states of Mg₂ have been evaluated as a function of internuclear separation. Received: 10 May 2001 / Accepted: 15 August 2001 / Published online: 30 October 2001     

Electrostatic and relaxation effects in ionization of molecular dimers and intermolecular complexes

The shifts in ionization energies which occur when a molecule is incorporated as an asymmetric dimer or in an intermolecular complex are analyzed theoretically. MO ? SCF calculations with 4–31G basis sets were performed on closed- and open-shell states of (HF)₂, H₂O·HF, and their valence–hole ions, as well as on the heterodimers incorporating the higher homologues CH₃F, CH₃OH, and (CH₃)₂O. The analysis concerns the influence of electrostatic, polarization, and charge transfer effects associated with complexation on the initial molecular state of each monomer system, as well as monomer–dimer differences in the electronic relaxation mechanism considered as a final state effect in the ionization process. The calculated ionization energy shifts which agree well with the experimental data available for (CH₃)₂O·HF, show that the shifts are dominated by electrostatic effects, but some effects arising from differences in molecular size and electric polarizability of the monomers can be discerned.     

Theoretical calculation of the electronic states with spin‐orbit effects of the molecule NaRb

The potential energy curves of the molecule NaRb have been calculated for the 60 low‐lying electronic states in the Ω‐representation. Using an ab‐initio method the calculation is based on nonempirical pseudo‐potential in the interval 3.0a_o ≤ R ≤ 44.0a_o of the internuclear distance. The spin‐orbit effects have been taken into account through a semiempirical spin‐orbit pseudo‐potential added to the electrostatic Hamiltonian with Gaussian basis sets for both atoms. The spectroscopic constants have been calculated for 42 states and the components of the spin‐orbit splitting have been identified for the states (1, 2, 5)³Π and (1, 2)³Δ. The comparison of the present results with those available in literature shows a good agreement, whereas the other results, to the best of our knowledge, are given here for the first time. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Ab Initio calculation of molecular, thermodynamic, and kinetic parameters of CF, CF2, CF3, and C2F5 radicals and CF4, CF3I, C2F4, and C2F6 molecules in gas

Enthalpies of formation of ground states of the gaseous particles CF, CF₂, C₂F₅, CF₄, CF₃I, C₂F₄, and C₂F₆ were calculated by ab initio method in the CCSD(T) approximation with extrapolation to the full basis and regard to the correlation energy. Their equilibrium geometrics, frequencies of normal vibrations, and other values were found by the B3LYP/aug-cc-pvdz method, from which thermodynamic functions within the range of 0–6000 K were calculated. Equilibrium constants were calculated from these functions, and then the information on the rate constants in the limit of high pressures was obtained.     

Systematic kinetic study of H2 release from the dimer of lithium amidoborane (LiNH2BH3)2

The four-step dehydrogenation of lithium amidoborane dimer (LiNH₂BH₃)₂ has been systematically simulated for the first time, and the respective rate constants have been calculated. Density functional theory has been used to optimize the molecular structure and ab initio direct kinetic theory has been applied to identify dehydrogenation mechanisms. The transition states were confirmed by intrinsic reaction coordinate calculations to insure the validity of our simulation and the barrier associated with each reaction was calculated. The Arrhenius equations of the four-step reactions (two pathways in all) were then obtained. The result indicated the dissociation maybe dimer way different from the traditional views. Our study has indicated a lower activation energy for dehydrogenation of the dimer compared to that of the monomer. The simulation is consistent with experimental observation because each step of the process requires increasingly higher energy. The study provides useful information on the properties and dehydrogenation mechanisms of metal-amidoborane compounds.     

Removal of dependencies from nearly complete basis sets. Calculations on the helium dimer

A method is devised for dealing with almost linearly dependent basis sets that contain large sets of bond functions. Using the largest of such basis sets, LARSAT, the second-order Møller-Plesset polarization dispersion energy of the helium dimer is calculated to be - 17.08 K at R = 5.6 bohrs. MR-SDCI calculations, employing a set of 37 reference configurations, were performed for the helium dimer with several basis sets at 4.0 and 5.6 bohrs. Size-extensivity corrections were included to take into account the R dependency of the size-extensivity error in MR-SDCI calculations. The He₂ interaction energies computed with basis LARSAT are - 10.92 K at 5.6 bohrs and 295.1 K at 4.0 bohrs. The 37-MR-SDCI calculations with basis LARSAT almost reproduce the He₂ full configuration interaction (CI) interaction energies computed with the same basis, at notably smaller cost. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 805–815, 1997     

Ab initio valence bond calculations. VII. HF,HF+, and H2F+

Ab initio valence bond calculations for the ground and excited states of HF and HF⁺ are presented. Total energies, equilibrium geometries, dissociation energies, dipole moments, and spectroscopic constants for HF and HF⁺ have been calculated. The photoelectron spectrum of HF has been examined and interpreted by means of the valence bond formalism. The ground state of the protonated species H₂F⁺ has been investigated.     

Triplet excited states of nitronaphthalenes. IV. 1,2- and 1,8-dinitronaphthalenes

Nanosecond flash photolysis of 1,2- and 1,8-dinitronaphthalenes (1,2-DNO₂N; 1,8-DNO₂N) in nonpolar and polar solvents shows transient species with absorption maxima and lifetimes dependent on solvent polarity. In deaerated n-hexane the absorption maxima and lifetimes (1/K) are 490 nm and 1.0 μsec for 1,2-DNO₂N and 550 nm and 2.5 μsec for 1,8-DNO₂N. In deaerated ethanol the corresponding values are 550 nm and 4.3 μsec for 1,2-DNO₂N and 590 nm and 5.3 μsec for 1,8-DNO₂N. The transient absorptions are attributed to the lowest triplet excited states T₁ of the 1,2- and 1,8-DNO₂N. The observed red shifts in the absorption maxima of the T₁ states are indicative of the extent to which electronic charge is transferred intramolecularly during the T₁ → T_n transitions. Furthermore, the increased lifetime of the T₁ states with increasing solvent polarity indicates the intramolecular charge transfer character of the T₁ states. Changes of dipole moments accompanying the T₁ → T_n transitions as well as rate constants for electron or proton transfer and hydrogen abstraction reactions involving the T₁ states of 1,2- and 1,8-DNO₂N and tributyl tin hydride (Bu₃SnH) as the hydrogen donor were determined together with the activation energy of the hydrogen abstraction reaction for the case of 1,2-DNO₂N. The spectroscopic and kinetic data obtained in this work demonstrate that the triplet states of 1,2- and 1,8-DNO₂N behave like n → π* states in nonpolar media while in polar solvents the n → π* character of these states is reduced with a simultaneous increase in their intramolecular charge transfer character.     

Semiempirical local spin: Theory and implementation of the ZILSH method for predicting Heisenberg exchange constants of polynuclear transition metal complexes

The local spin formalism ( 3 ) for computing expectation values 〈S_A · S_B〉 that appear in the Heisenberg spin model has been extended to semiempirical single determinant wave functions. An alternative derivation of expectation values in restricted and unrestricted cases is given that takes advantage of the zero differential overlap (ZDO) approximation. A formal connection between single determinant wave functions (which are not in general spin eigenfunctions) and the Heisenberg spin model was established by demonstrating that energies of single determinants that are eigenfunctions of the local spin operators with eigenvalues corresponding to high‐spin radical centers are given by the same Heisenberg coupling constants {J_AB} that describe the true spin states of the system. Unrestricted single determinant wave functions for transition metal complexes are good approximations of local spin eigenfunctions when the metal d orbitals are local in character and all unpaired electrons on each metal have the same spin (although spins on different metals might be reversed). Good approximations of the coupling constants can then be extracted from local spin expectation values 〈S_A · S_B〉 energies of the single determinant wave functions. Once the coupling constants are obtained, diagonalization of the Heisenberg spin Hamiltonian provides predictions of the energies and compositions of the spin states. A computational method is presented for obtaining coupling constants and spin‐state energies in this way for polynuclear transition metal complexes using the intermediate neglect of differential overlap Hamiltonian parameterized for optical spectroscopy (INDO/S) in the ZINDO program. This method is referred to as ZILSH, derived from ZINDO, Davidson's local spin formalism, and the Heisenberg spin model. Coupling constants and spin ground states obtained for 10 iron complexes containing from 2 to 6 metals are found to agree well with experimental results in most cases. In the case of the complex [Fe₆O₃(OAc)₉(OEt)₂(bpy)₂]⁺, a priori predictions of the coupling constants yield a ground‐state spin of zero, in agreement with variable‐temperature magnetization data, and corroborate spin alignments proposed earlier on the basis of structural considerations. This demonstrates the potential of the ZILSH method to aid in understanding magnetic interactions in polynuclear transition metal complexes. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Theoretical electronic structure of the molecule ScI

Theoretical investigation of the 18 lowest electronic states of the molecule ScI in the representation ^2S+1Λ^(±) has been performed via CASSCF and MRCI (single and double excitation with Davidson correction) calculations. To the best of our knowledge these calculated electronic states are the first ones from ab initio methods. Thirteen electronic states between 4,500 cm^?1 and 21,000 cm^?1 have been studied for the first time and have not yet been observed experimentally. The harmonic frequency ω_e, the internuclear distance R_e, the electronic transition energy with respect to the ground state T_e, and the rotational constant B_e have been calculated for the considered electronic states. By using the canonical functions approach the eigenvalues E_υ and the rotational constants B_υ have also been calculated for the six lowest‐lying electronic states. The comparison of these results with the theoretical and the experimental data available in the literature shows a good agreement. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Ab initio studies of the electronic structure and density of states of metallic beryllium

Hartree–Fock–Roothaan studies are reported for low-lying electronic states of metallic beryllium as modeled by a moiety of 135 beryllium atoms. The system corresponds to 16 coordination shells of a central Be with internuclear separations derived from the lattice constants of the bulk metal. The calculations become tractable by use of the full D_3h symmetry of the system at both the integrals and self-consistent-field stages and by employing ab initio effective potentials for the 1s electrons of each beryllium atom. Ionization potentials, binding energies, orbital energies, electric field gradients, nuclear-electrostatic potentials, diamagnetic shielding constants, second moments, and Mulliken populations are calculated for selected electronic states. The calculated ionization potential for the lowest state agrees to within 10% of the experimental bulk work function. A density-of-states analysis for that state is reported and compared with band structure calculations.     

Ground states of molecules. 53.MNDO calculations for molecules containing chlorine

MNDO calculations of heats of formation, dipole moments, ionization potentials, and structures are reported for a wide range of compounds containing chlorine in its characteristic valence state (Cl^I) and one or more of the elements H, B, Be, C, N, O, and F. The calculated errors in the heats of formation and the dipole moments are not significantly greater than those previously reported for compounds containing no chlorine. First vertical ionization potentials were on average 0.95 eV too high. The ordering of higher cationic states was found to be correct, even for species such as Cl₂O, Cl₂, and HOCl, where ab initio–Koopmans' theorem calculations predict the incorrect ordering. The calculated energies and geometries of compounds such as CIF₃ are qualitatively incorrect, probably because of the lack of 3d atomic orbitals in the orbital basis set.     

Equilibrium structure and spectroscopic constants of maleic anhydride

The quadratic, cubic, and semi-diagonal quartic force fields of maleic anhydride have been calculated at the MP2 level of theory employing the cc-pVTZ basis set. The spectroscopic constants derived from the force field are in excellent agreement with the corresponding experimental values. The semi-experimental equilibrium structure has been derived from experimental ground state rotational constants and rovibrational corrections calculated from the cubic force field. This semi-experimental equilibrium structure is in excellent agreement with the ab initio structures computed at the CCSD(T) level of theory and it is closer to the ab initio structure than the purely experimental (or empirical) structures r ₀, r _m⁽¹⁾, and r _m⁽²⁾ obtained by microwave spectroscopy as well as the equilibrium structure derived from gas-phase electron diffraction data.