The spin-forbidden reaction CH(2∏) + N2→ HCN + N(4S) revisited I. Ab initio study of the potential energy surfaces

High-level ab initio electronic structure theories have been applied to investigate the detailed reaction mechanism of the spin-forbidden reaction CH(²∏) + N₂ → HCN + N(⁴S). The G2M(RCC) calculations provide accurate energies for the intermediates and transition states involved in the reaction, whereas the B3LYP/6-311G(d,p) method overestimates the stability of some intermediates by as much as about 10 kcal/mol. A few new structures have been found for both the doublet and quartet electronic states, which are mainly involved in the dative pathways. However, due to the higher energies of these structures, the dominant mechanism remains the one involving the C ₂ intersystem-crossing step. The C ₂ minima on the seam of crossing (MSX) structures and the spin-orbit coupling between the doublet and quartet electronic states are rather close to those found in previous studies. Vibrational frequencies orthogonal to the normal of the seam which have been applied in a separate publication to calculate the rate of the CH(²∏) + N₂ → HCN + N(⁴S) reaction with a newly proposed nonadiabatic transition-state theory for spin-forbidden reactions have been calculated at the MSX from first principles. Received: 23 June 1998 / Accepted: 21 September 1998 / Published online: 8 February 1999     

Diatomic potential well depths from analyses of high-resolution electron energy spectra for autoionizing collision complexes

High-resolution energy spectra of electrons released in Penning ionization collisions of metastable rare gas atoms Rg*(ns) (Rg = He, Ne, Ar, Kr, Xe) with several open-shell and closed-shell atoms are analyzed to determine the well depth of the potential energy curve which describes the respective autoionizing collision complex. We thereby elucidate trends in the chemical interaction of Rg* with closed-shell target atoms A(ns ²) and establish a basis for detailed comparison with the respective interactions involving open-shell, ground state alkali atoms A(ns). From electron energy spectra due to␣associative ionization (RgH⁺ formation) in Rg* + H(1s) collisions, we determine binding energies for the RgH⁺(¹Σ) ground state potential (Rg = Ne, Ar, Kr, Xe) with uncertainties around 0.03 eV. Received: 30 June 1998 / Accepted: 5 August 1998 / Published online: 28 October 1998     

Theoretical study of the valence π → π* excited states of polyacenes: anthracene and naphthacene

The valence π → π ^* excited states of anthracene and naphthacene are studied with multireference perturbation theory with complete active space self-consistent field reference functions. The predicted spectra provide a consistent assignment of all one- and two-photon spectra and T-T spectra of low-lying valence π → π ^* excited states of anthracene and naphthacene. The present theory predicts the valence π → π ^* excitation energies with an accuracy of 0.15 eV for anthracene and of 0.25 eV or better for naphthacene. The excited states of anthracene and naphthacene are compared with those of benzene and naphthalene studied previously. The present calculations predict that, going from anthracene to naphthacene, there is a symmetry reversal of the two lowest singlet state transitions, but not for the triplet, just as indicated by the experimental data. Some general trends of polyacene excited states are discussed based on the calculated results for benzene to naphthacene. Conclusive results obtained for anthracene and naphthacene can be used as a model for understanding the excited states of larger polyacenes. Received: 22 April 1998 / Accepted: 6 July 1998 / Published online: 28 September 1998     

Ab initio calculations find 2,2-disilylcyclopentane-1,3-diyl is a singlet diradical with a high barrier to ring closure

Ab initio calculations on the lowest singlet and triplet states of 2,2-disilylcyclopentane-1,3-diyl find that the singlet lies well below the triplet. The C ₂ singlet diradical is calculated to be a minimum on the potential energy surface with an enthalpic barrier to ring closure of ΔH ^‡ ₂₉₈ = 13.5 kcal/mol at the CASPT2/6-31G* level of theory. The energy of the 1,3-divinyl-substituted singlet diradical is calculated to be only 0.8 kcal/mol higher than that of 5,5-disilyl-1,3-divinylbicyclo[2.1.0]pentane at this level of theory, but the transition state for their equilibration is computed to be 12.8 kcal/mol above the diradical in energy. Received: 2 July 1998 / Accepted: 4 August 1998 / Published online: 16 November 1998     

Degenerate symmetry-adapted perturbation theory of weak interactions between closed- and open-shell monomers: application to Rydberg states of helium hydride

Symmetry-adapted perturbation theory is extended to the (quasi) degenerate, open-shell case. The new formalism is tested in calculations of the interaction energies for a helium atom in the ground state interacting with an excited hydrogen atom. It is shown that the method gives satisfactory results if the coupling with higher Rydberg states of the dimer is small, as is the case for the A²Σ⁺,B²Π,E²Π,3²Π, and 1²Δ states of HeH. For the C²Σ⁺ state convergence of the method is very slow, but it can be improved by including the n=3 states in the model space. Received: 3 June 1998 / Accepted: 9 September 1998 / Published online: 7 December 1998     

A theoretical study on the first ionic states of vinyl fluoride, vinyl chloride, trifluoroethylene, and trichloroethylene with an analysis of the vibrational structures of the photoelectron spectra

Ab initio calculations have been performed to study the molecular structures and the vibrational levels of the first ionic states of vinyl fluoride, vinyl chloride, trifluoroethylene, and trichloroethylene. The equilibrium molecular structures and vibrational modes of these states are presented. The theoretical ionization intensity curves including the vibrational structures are also presented and compared with the photoelectron spectra. Received: 11 September 1998 / Accepted: 13 October 1998 / Published online: 1 February 1999     

Pseudopotential study of the ground and excited states of Yb2

The ground state of the van der Waals-type lanthanide dimer Yb₂ has been studied by means of relativistic energy-consistent ab initio pseudopotentials using three different core definitions. Electron correlation was treated by coupled-cluster theory, whereby core-valence correlation effects have been accounted for either explicitly by correlating the energetically highest coreorbitals or implicitly by means of an effective core-polarization potential. Results for the first and second atomic ionization potentials, the atomic dipole polarizability, and the spectroscopic constants of the molecular ground state are reported. Low-lying excited states have been investigated with spin-orbit configuration interaction calculations. It is also demonstrated for the whole lanthanide series that correlation effects due to the atomic-like, possibly open 4f-shell in lanthanides can be modeled effectively by adding a core-polarization potential to pseudopotentials attributing the 4f-shell to the core. Received: 3 April 1998 / Accepted: 27 July 1998 / Published online: 9 October 1998     

A theoretical study on the ionization of OCS with an analysis of vibrational structures of the photoelectron spectrum

 Ab initio calculations have been performed to study the molecular structures and vibrational levels of the four low-lying ionic states (1, 2²Π, and 1, 2²Σ⁺) of carbonyl sulfide. The global regions of the potential-energy surfaces have been obtained by multireference single and double excitation configuration interaction calculations. Vibrational calculations using explicit vibrational Hamiltonians have been used for vibrational analysis. The equilibrium molecular structures and a vibrational analysis of the four ionic states are presented. The theoretical ionization intensity curves including the vibrational structures of the ionic states are also presented and are compared with the photoelectron spectrum. Received: 20 January 2001 / Accepted: 22 August 2001 / Published online: 30 October 2001     

Dissociative recombination of HeH+. I. Rovibrational spectrum of HeH Rydberg states

The repulsive ground electronic state X²Σ⁺ of HeH is strongly coupled to the Rydberg states at small interatomic distances. Such large couplings also occur between some of the Rydberg states. HeH⁺ ions that capture an electron in a Rydberg state end up in separated He and H atoms by indirect predissociation. This paper presents a study of potential functions and pertinent matrix elements involving the lowest electronic states: the ²Σ⁺ states, X, A, C, and D, and the ²Π states B and E. Individual transition rates as well as total radiative and non-radiative lifetimes have been computed for the lowest vibrational and rotational levels. Received: 22 June 1998 / Accepted: 21 August 1998 / Published online: 12 October 1998     

Possible long-lived quartet resonance states of CO−

Two states of CO⁻ which are expected to live longer than the well-known ²Π shape resonance are examined and the results obtained are compared with the analogous states of the isoelectronic N₂ ⁻ system. The ⁴Σ⁻ and ⁴Π states are found to be stable with respect to dissociation. The latter state is also stable with respect to single-electron loss to its parent neutral state. The lifetime of this state is supposed to be long, but not necessarily long enough to allow detection in a mass spectrometer. Received: 28 May 1998 / Accepted: 23 July 1998 / Published online: 9 October 1998     

Theoretical treatment of predissociation of the CO (3sσ) B and (3pσ) C1Σ+ Rydberg states based on a rigorous adiabatic representation

Ab initio multireference configuration interaction calculations for adiabatic potential curves, nonadiabatic couplings 〈φ _i (R,r)|d/dR|φ _j (R,r)〉 and 〈φ _i (R,r)|d²/dR ²|φ _j (R,r)〉, and nuclear kinetic energy corrections 〈dφ _i (R,r)/dR|dφ _i (R,r)/dR〉 for the (3sσ) B and (3pσ) C¹Σ⁺ Rydberg states of the CO molecule have been carried out. The energy positions and predissociation linewidths for the observed vibrational levels of these two states have been determined in a rigorous adiabatic representation by the complex scaling method employing a basis of complex scaled harmonic vibrational functions in conjunction with the Gauss-Hermite quadrature method to evaluate the complex Hamiltonian matrix elements. The present treatment correctly reproduces the observed trends in energies and line broadening for vibrational levels of the B¹Σ⁺ state and represents an improvement over the previous treatment in literature. The errors in the determined spacings of the v = 0–4 vibrational levels of the C¹Σ⁺ state are less than 2% compared with measured data. The predissociation linewidths for the v=3,4 levels of the C¹Σ⁺ state are found to be 4.9 and 8.9 cm⁻¹, respectively, in good agreement with the observed values. Received: 23 March 1998 / Accepted: 27 July 1998 / Published online: 9 October 1998     

Micellar properties of alkyldimethylphenylammonium bromides in water

The synthesis and methods applied for the purification of dodecyl-, tetradecyl-, and hexadecyldimethylphenylammonium bromides are described. The results of surface tension measurements of aqueous solutions of these surfactants show that slight amounts of strongly surface-active nonionic impurities are persistent in the crystalline materials presumably due to their low thermal stability. The Critical micelle concentration (cmc) and the degree of ionization (β) of the micelles of the salts studied in aqueous solutions were determined at 25 °C from specific conductivity versus molality plots. The temperature dependence of the cmc and of β of the tetradecyl homologue was measured in the range 4–34 °C. A minimum cmc amounting to 1.20 mmol/kg was determined at about 14 °C. The values of β were found to grow linearly with temperature. From these results, the standard Gibbs energy, the enthalpy and the entropy of the process of micellization were obtained by application of the pseudo-phase-separation model. Enthalpy and entropy show a compensation effect in their contribution to the Gibbs energy. At low temperatures the process of micellization is driven mainly by the entropic term, whereas with increasing temperature the enthalpic term becomes predominant. At the temperature of the minimum cmc, the value of the enthalpy is far from being zero because of the important contribution of the (β/T)R ln X _cmc term. Received: 27 July 1998 Accepted in revised form: 15 December 1998     

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     

A theoretical study of the second band of the photoelectron spectrum of benzene with an analysis of the vibrational structure

Ab initio calculations have been performed to study the second band of the photoelectron spectrum of benzene with analysis of the vibrational structure. The ²E_{2 g} and ²A_{2 u} states contribute to this band. In this study each contribution to the band is discussed. We propose that the onset of the band should be assigned to the 0-0 transition of the ²A_{2 u} state. Received: 20 August 1998 / Accepted: 15 January 1999 / Published online: 7 June 1999     

Generalization of coupled-cluster response theory to multireference expansion spaces: application of the coupled-cluster singles and doubles effective Hamiltonian

Employing separate cluster ansatz in time-independent and time-dependent wave-operators, coupled-cluster (CC) response theory is generalized to multireference (MR) expansion spaces. For state energies, this corresponds to the MR secular problem with an arbitrary similarity-transformed effective Hamiltonian, H˜=Ω⁻¹ HΩ. The effective Hamiltonian can be generated via size-extensive CC methods. Thus the states in MR linear response theory (MRLRT) maintain the usual CC core-extensive properties. We have used the Gelfand unitary group basis of the spin-adapted configurations to construct the matrix of H˜ in the MR excitation space. As a preliminary application, the CC singles and doubles effective Hamiltonian is applied to excitation and photoionization energies of the CH⁺ and N₂ molecules, and is compared with experimental results and results from other numerical procedures including conventional CC linear response theory (CC-LRT), MR and full configuration interaction (MRCI and FCI) methods. The numerical results indicate that MRLRT reproduces valence and external excited states quantitatively, combining the best features of CC-LRT and MRCI. Received: 2 July 1998 / Accepted: 28 August 1998 / Published online: 11 November 1998     

Analytic UHF-CCSD(T) second derivatives: implementation and application to the calculation of the vibration-rotation interaction constants of NCO and NCS

An implementation of analytic open-shell UHF-CCSD(T) second derivatives is presented. To demonstrate applicability and test the accuracy of the UHF-CCSD(T) approach for the determination of spectroscopical parameters, vibration-rotation interaction constants were calculated for the ground (1²Π) and first electronically excited (1²Σ) states of the NCO and NCS radicals. In addition, harmonic vibrational frequencies for both states, the Renner-Teller parameter for the ground state, as well as the 1²Π→1²Σ excitation energy are reported. While the computed values are in good agreement with reliable experimental information for NCO, most of the data presented for NCS are predictions of quantities not well known from experiment. Received: 4 May 1998 / Accepted: 27 July 1998 / Published online: 12 October 1998     

Perspective on “An extended Hückel theory. I. Hydrocarbons”

 The size-consistent self-consistent matrix dressing method has been applied on an open-shell single-configuration reference state. Once the reference state is converged, several low-lying roots can be obtained for the dressed configuration interaction (CI) matrices of appropriate symmetry. The CI matrices were built with a complete-active-space singles and doubles CI method in order to deal properly with multiconfiguration excited states. The vertical ionization and ionization–excitation transitions are obtained from the difference to the closed shell ground-state energy of the neutral molecule. The method has been applied to NH⁺ ₃ and N⁺ ₂ using atomic natural orbital basis sets and state-average adapted molecular orbitals. Two ²A₁ states, very similar and showing great mixing of the (2a _l ⁻¹) and (3a _l ⁻²5a _l ¹) determinants, can be assigned to the broad asymmetric band at 27.6 ± 2 eV in the photoelectron spectrum of NH₃. The possible contribution of a ²Π _g (3σ _g ⁻²1π _g ¹) state to the A shake-up peak of N₂ at 24.6 eV is also discussed. Other states, doublets and quadruplets, are reported for both systems up to 30 eV for NH₃ and 37 eV for N₂. Received: 16 September 1999 / Accepted: 3 February 2000 / Published online: 2 May 2000     

Ion channels of biological membranes: prediction of single channel conductance

The Poisson-Boltzmann equation was solved numerically for models of the pore regions of the Shaker K⁺ channel and of two glycoporins (LamB and ScrY) to yield electrostatic potential profiles along the pore axes. From these potential profiles, single-channel current-voltage (I–V) relations were calculated. The importance of a proper treatment of the ionisation state of two rings of aspartate sidechains at the mouth of the K⁺ channel pore emerged from such calculations. The calculated most likely state, in which only two of the eight aspartate sidechains were deprotonated, yielded better agreement with experimental conductance data. An approximate calculation of single-channel conductances based simply on pore geometry yielded very similar conductance values for the two glycoporins. This differed from an␣experimentally determined conductance ratio of ScrY:LamB=10:1. Preliminary electrostatics calculations appeared to reproduce the observed difference in conductance between the two glycoporins, confirming that single-channel conductance is determined by electrostatic as well as geometric considerations. Received: 25 May 1998 / Accepted: 4 August 1998 / Published online: 2 November 1998