Curl condition for a four-state Born-Oppenheimer system employing the Mathieu equation

When a group of four states forms a subspace of the Hilbert space, i.e., appears to be strongly coupled with each other but very weakly interacts with all other states of the entire space, it is possible to express the nonadiabatic coupling (NAC) elements either in terms of s or in terms of electronic basis function angles, namely, mixing angles presumably representing the same sub-Hilbert space. We demonstrate that those explicit forms of the NAC terms satisfy the curl conditions--the necessary requirements to ensure the adiabatic-diabatic transformation in order to remove the NAC terms (could be often singular also at specific point(s) or along a seam in the configuration space) in the adiabatic representation of nuclear SE and to obtain the diabatic one with smooth functional form of coupling terms among the electronic states. In order to formulate extended Born-Oppenheimer (EBO) equations [J. Chem. Phys. 2006, 124, 074101] for a group of four states, we show that there should exist a coordinate independent ratio of the gradients for each pair of ADT/mixing angles leading to zero curls and, thereafter, provide a brief discussion on its analytical validity. As a numerical justification, we consider the first four eigenfunctions of the Mathieu equation to demonstrate the interesting features of nonadiabatic coupling (NAC) elements, namely, the validity of curl conditions and the nature of curl equations around CIs.     

On the construction of diabatic and adiabatic potential energy surfaces based on ab initio valence bond theory

A theoretical model is presented for deriving effective diabatic states based on ab initio valence bond self-consistent field (VBSCF) theory by reducing the multiconfigurational VB Hamiltonian into an effective two-state model. We describe two computational approaches for the optimization of the effective diabatic configurations, resulting in two ways of interpreting such effective diabatic states. In the variational diabatic configuration (VDC) method, the energies of the diabatic states are variationally minimized. In the consistent diabatic configuration (CDC) method, both the configuration coefficients and orbital coefficients are simultaneously optimized to minimize the adiabatic ground-state energy in VBSCF calculations. In addition, we describe a mixed molecular orbital and valence bond (MOVB) approach to construct the CDC diabatic and adiabatic states for a chemical reaction. Note that the VDC-MOVB method has been described previously. Employing the symmetric S(N)2 reaction between NH(3) and CH(3)NH(3)(+) as a test system, we found that the results from ab initio VBSCF and from ab initio MOVB calculations using the same basis set are in good agreement, suggesting that the computationally efficient MOVB method is a reasonable model for VB simulations of condensed phase reactions. The results indicate that CDC and VDC diabatic states converge, respectively, to covalent and ionic states as the molecular geometries are distorted from the minimum of the respective diabatic state along the reaction coordinate. Furthermore, the resonance energy that stabilizes the energy of crossing between the two diabatic states, resulting in the transition state of the adiabatic ground-state reaction, has a strong dependence on the overlap integral between the two diabatic states and is a function of both the exchange integral and the total diabatic ground-state energy.     

On the (Emultiply sign in circlee)-Jahn-Teller conical intersections in the 3p(E') and 3d(E") Rydberg electronic states of triatomic hydrogen

The static and dynamic aspects of the Jahn-Teller (JT) interactions in the 3p(E') and 3d(E") Rydberg electronic states of H3 are analyzed theoretically. The static aspects are discussed based on recent ab initio quantum chemistry results, and the dynamic aspects are examined in terms of the vibronic spectra and nonradiative decay behavior of these states. The adiabatic potential-energy surfaces of these degenerate electronic states are derived from extensive ab initio calculations. The calculated adiabatic potential-energy surfaces are diabatized following our earlier study on this system in its 2p(E') ground electronic state. The nuclear dynamics on the resulting conically intersecting manifold of electronic states is studied by a time-dependent wave-packet approach. Calculations are performed both for the uncoupled and coupled state situations in order to understand the importance of nonadiabatic interactions due to the JT conical intersections in these excited Rydberg electronic states.     

Systematic corrections to the Born-Oppenheimer approximation

A method for providing systematic diabatic corrections to the Born-Oppenheimer approximation is presented. We begin with an adiabatic expansion of the exact vibronic wavefunctions and, via the molecular Hamiltonian, develop expressions for the diabatic terms in the Schrödinger equation. We then derive recursion relations which allow one to introduce the diabatic interactions to any desired degree of approximation. As an example, the first approximation (beyond the Born-Oppenheimer approximation) is discussed explicitly. In passing, we assess some of the common misconceptions associated with the Born-Oppenheimer approximation.     

Adiabatic versus diabatic descriptions of the lowest Rydberg and valence 1Σ+ states of HCl

In this contribution we first report new ab initio self-consistent field configuration interaction calculations of the first excited adiabatic potential of (1)Σ(+) symmetry, the 2(1)Σ(+) or B(1)Σ(+) state, which presents two minima and can thus be seen as made up of the Rydberg E(1)Σ(+) and the valence V(1)Σ(+) states. Based on the computed 2(1)Σ(+) potential, we devised a theoretical procedure to compute the vibronic structure in order to try to explain the energy levels observed in the region above 76 254.4 cm(-1) which display an irregular vibrational structure, indicative of spectral perturbations. We try to find out which representation of the electronic states, the diabatic or the adiabatic one, is best suited to replicate the lowest observed vibronic levels of the E and V states. To this end, we deduce, from the 2(1)Σ(+) potential and its complementary adiabatic potential, two diabatic potentials. We then carry out a coupled equation treatment based on these diabatic potentials. The results of this treatment indicate that, in the present case, the adiabatic representation is better than the diabatic one to describe the observed vibronic levels. This is due, as expected, to the existence of a strong electrostatic interaction between the two diabatic potentials.     

A rigorous approach to the formulation of extended Born‐Oppenheimer equation for a three‐state system

If a coupled three‐state electronic manifold forms a sub‐Hilbert space, it is possible to express the non‐adiabatic coupling (NAC) elements in terms of adiabatic–diabatic transformation (ADT) angles. Consequently, we demonstrate: (a) Those explicit forms of the NAC terms satisfy the Curl conditions with non‐zero Divergences; (b) The formulation of extended Born‐Oppenheimer (EBO) equation for any three‐state BO system is possible only when there exists coordinate independent ratio of the gradients for each pair of ADT angles leading to zero Curls at and around the conical intersection(s). With these analytic advancements, we formulate a rigorous EBO equation and explore its validity as well as necessity with respect to the approximate one (Sarkar and Adhikari, J Chem Phys 2006, 124, 074101) by performing numerical calculations on two different models constructed with different chosen forms of the NAC elements. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Ab initio treatment of the chemical reaction precursor complex Br(2P)-HCN. 1. Adiabatic and diabatic potential surfaces

The three adiabatic potential surfaces of the Br(2P)-HCN complex that correlate to the 2P ground state of the Br atom were calculated ab initio. With the aid of a geometry-dependent diabatic mixing angle, also calculated ab initio, these adiabatic potential surfaces were transformed into a set of four diabatic potential surfaces required to define the full 3 x 3 matrix of diabatic potentials. Each of these diabatic potential surfaces was expanded in terms of the appropriate spherical harmonics in the atom-linear molecule Jacobi angle theta. The dependence of the expansion coefficients on the distance R between Br and the HCN center of mass and on the CH bond length was fit to an analytic form. For HCN in its equilibrium geometry, the global minimum with De = 800.4 cm(-1) and Re = 6.908a0 corresponds to a linear Br-NCH geometry, with an electronic ground state of Sigma symmetry. A local minimum with De = 415.1 cm-1, Re = 8.730a0, and a twofold degenerate Pi ground state is found for the linear Br-HCN geometry. The binding energy, De, depends strongly on the CH bond length for the Br-HCN complex and much less strongly for the Br-NCH complex, with a longer CH bond giving stronger binding for both complexes. Spin-orbit coupling was included and diabatic states were constructed that correlate to the ground 2P3/2 and excited 2P1/2 spin-orbit states of the Br atom. For the ground spin-orbit state with electronic angular momentum j = (3/2) the minimum in the potential for projection quantum number omega = +/-(3/2) coincides with the local minimum for linear Br-HCN of the spin-free case. The minimum in the potential for projection quantum number omega = +/-(1/2) occurs for linear Br-NCH but is considerably less deep than the global minimum of the spin-free case. According to the lowest spin-orbit coupling included adiabatic potential the two linear isomers, Br-NCH and Br-HCN, are about equally stable. In the subsequent paper, we use these potentials in calculations of the rovibronic states of the Br-HCN complex.     

Extended Born-Oppenheimer equation for a three-state system

We present explicit forms of nonadiabatic coupling (NAC) elements of nuclear Schrodinger equation (SE) for a coupled three-state electronic manifold in terms of mixing angles of real electronic basis functions. If the adiabatic-diabatic transformation (ADT) angles are the mixing angles of electronic bases, ADT matrix transforms away the NAC terms and brings diabatic form of SE. ADT and NAC matrices are shown to satisfy a curl condition with nonzero divergence. We have demonstrated that the formulation of extended Born-Oppenheimer (EBO) equation from any three-state BO system is possible only when there exists a coordinate-independent ratio of the gradients for each pair of mixing angles. On the contrary, since such relations among the mixing angles lead to zero curl, we explore its validity analytically around conical intersection(s) and support numerically considering two nuclear-coordinate-dependent three surface BO models. Numerical calculations are performed by using newly derived diabatic and EBO equations and expected transition probabilities are obtained.     

Direct calculation of coupled diabatic potential-energy surfaces for ammonia and mapping of a four-dimensional conical intersection seam

We used multiconfiguration quasidegenerate perturbation theory and the fourfold-way direct diabatization scheme to calculate ab initio potential-energy surfaces at 3600 nuclear geometries of NH3. The calculations yield the adiabatic and diabatic potential-energy surfaces for the ground and first electronically excited singlet states and also the diabatic coupling surfaces. The diabatic surfaces and coupling were fitted analytically to functional forms to obtain a permutationally invariant 2 x 2 diabatic potential-energy matrix. An analytic representation of the adiabatic potential-energy surfaces is then obtained by diagonalizing the diabatic potential-energy matrix. The analytic representation of the surfaces gives an analytic representation of the four-dimensional conical intersection seam which is discussed in detail.     

Ab initio treatment of the chemical reaction precursor complex Cl(2P)-HF. 1. Three-dimensional diabatic potential energy surfaces

The three adiabatic potential surfaces of the Cl(2P)-HF complex that correlate with the 2P ground state of the Cl atom were calculated with the ab initio RCCSD(T) method (partially spin-restricted coupled cluster theory including single and double excitations and perturbative correction for the triples). With the aid of a geometry-dependent diabatic mixing angle, calculated by the complete active space self-consistent field (CASSCF) and multireference configuration-interaction (MRCI) methods, these adiabatic potential surfaces were converted to a set of four distinct diabatic potential surfaces required to define the full 3 x 3 matrix of diabatic potentials. Each of these diabatic potential surfaces was expanded in terms of the appropriate spherical harmonics in the angle theta between the HF bond axis r and the Cl-HF intermolecular axis R. The dependence of the expansion coefficients on the Cl-HF distance R and the HF bond length r(HF) was fit to an analytic form. The strongest binding occurs for the hydrogen-bonded linear Cl-HF geometry, with D(e) = 676.5 cm(-1) and R(e) = 6.217 a0 when r(HF) = r(e) = 1.7328 a0. This binding energy D(e) depends strongly on r(HF), with larger r(HF) causing stronger binding. An important contribution to the binding energy is provided by the interaction between the quadrupole moment of the Cl(2P) atom and the dipole of HF. In agreement with this electrostatic picture, the ground state of linear Cl-HF is a 2-fold degenerate electronic Pi state. For the linear Cl-FH geometry the states are in opposite order, i.e., the Sigma state is lower in energy than the Pi state. The following paper in this issue describes full three-dimensional computations of the bound states of the Cl-HF complex, based on the ab initio diabatic potentials of this paper.     

Structural and spectroscopic study of the LiRb molecule beyond the Born-Oppenheimer approximation

Adiabatic and diabatic potential energy curves and the permanent and transition dipole moments of the low-lying electronic states of the LiRb molecule dissociating into Rb(5s, 5p, 4d, 6s, 6p, 5d, 7s, 6d) + Li(2s, 2p) have been investigated. The molecular calculations are performed with an ab initio approach based on nonempirical pseudopotentials for Rb(+) and Li(+) cores, parametrized l-dependent core polarization potentials and full configuration interaction calculations. The derived spectroscopic constants (R(e), D(e), T(e), ω(e), ω(e)x(e), and B(e)) of the ground state and lower excited states are in good agreement with the available theoretical works. However, the 8-10(1)Σ(+), 8-10(3)Σ(+), 6(1,3)Π, and 3(1,3)Δ excited states are studied for the first time. In addition, to the potential energy, accurate permanent and transition dipole moments have been determined for a wide interval of internuclear distances. The permanent dipole moment of LiRb has revealed ionic characters both relating to electron transfer and yielding Li(-)Rb(+) and Li(+)Rb(-) arrangements. The diabatic potential energy for the (1,3)Σ(+), (1,3)Π, and (1,3)Δ symmetries has been performed for this molecule for the first time. The diabatization method is based on variational effective Hamiltonian theory and effective metric, where the adiabatic and diabatic states are connected by an appropriate unitary transformation.     

A time-dependent wave packet study of the vibronic and spin-orbit interactions in the dynamics of Cl((2)P)+H(2)-->HCl(X (1)Sigma(g) (+))+H((2)S) reaction

We investigate the vibronic and spin-orbit (SO) coupling effects in the state-selected dynamics of the title reaction with the aid of a time-dependent wave packet approach. The ab initio potential energy surfaces of Capecchi and Werner [Science 296, 715 (2002)] have been employed for this purpose. Collinear approach of the Cl((2)P) atom to the H(2) molecule splits the degeneracy of the (2)P state and gives rise to (2)Sigma and (2)Pi electronic states. These two surfaces form a conical intersection at this geometry. These states transform as 1 (2)A('), 1 (2)A("), and 2 (2)A('), respectively, at the nonlinear configurations of the nuclei. In addition, the SO interaction due to Cl atom further splits these states into (2)Sigma(1/2), (2)Pi(3/2), and (2)Pi(1/2) components at the linear geometry. The ground-state reagent Cl((2)P(3/2))+H(2) correlates with (2)Sigma(1/2) and (2)Pi(3/2), where as the SO excited reagent Cl(*)((2)P(1/2))+H(2) correlates with (2)Pi(1/2) at the linear geometry. In order to elucidate the impact of the vibronic and SO coupling effects on the initial state-selected reactivity of these electronic states we carry out quantum scattering calculations based on a flux operator formalism and a time-dependent wave packet approach. In this work, total reaction probabilities and the time dependence of electronic population of the system by initiating the reaction on each of the above electronic states are presented. The role of conical intersection alone on the reaction dynamics is investigated with a coupled two-state model and for the total angular momentum J=0 (neglecting the electronic orbital angular momentum) both in a diabatic as well as in the adiabatic electronic representation. The SO interaction is then included and the dynamics is studied with a coupled three-state model comprising six diabatic surfaces for the total angular momentum J=0.5 neglecting the Coriolis Coupling terms of the Hamiltonian. Companion calculations are carried out for the uncoupled adiabatic and diabatic surfaces in order to explicitly reveal the impact of two different surface coupling mechanisms in the dynamics of this prototypical reaction.     

Interpolation of diabatic potential energy surfaces

A method is presented for constructing diabatic potential energy matrices from ab initio quantum chemistry data. The method is similar to that reported previously for single adiabatic potential energy surfaces, but correctly accounts for the nuclear permutation symmetry of diabatic potential energy matrices and other complications that arise from the derivative coupling of electronic states. The method is tested by comparison with an analytic model for the two lowest energy states of H(3).     

A diabatic representation including both valence nonadiabatic interactions and spin-orbit effects for reaction dynamics

A diabatic representation is convenient in the study of electronically nonadiabatic chemical reactions because the diabatic energies and couplings are smooth functions of the nuclear coordinates and the couplings are scalar quantities. A method called the fourfold way was devised in our group to generate diabatic representations for spin-free electronic states. One drawback of diabatic states computed from the spin-free Hamiltonian, called a valence diabatic representation, for systems in which spin-orbit coupling cannot be ignored is that the couplings between the states are not zero in asymptotic regions, leading to difficulties in the calculation of reaction probabilities and other properties by semiclassical dynamics methods. Here we report an extension of the fourfold way to construct diabatic representations suitable for spin-coupled systems. In this article we formulate the method for the case of even-electron systems that yield pairs of fragments with doublet spin multiplicity. For this type of system, we introduce the further simplification of calculating the triplet diabatic energies in terms of the singlet diabatic energies via Slater's rules and assuming constant ratios of Coulomb to exchange integrals. Furthermore, the valence diabatic couplings in the triplet manifold are taken equal to the singlet ones. An important feature of the method is the introduction of scaling functions, as they allow one to deal with multibond reactions without having to include high-energy diabatic states. The global transformation matrix to the new diabatic representation, called the spin-valence diabatic representation, is constructed as the product of channel-specific transformation matrices, each one taken as the product of an asymptotic transformation matrix and a scaling function that depends on ratios of the spin-orbit splitting and the valence splittings. Thus the underlying basis functions are recoupled into suitable diabatic basis functions in a manner that provides a multibond generalization of the switch between Hund's cases in diatomic spectroscopy. The spin-orbit matrix elements in this representation are taken equal to their atomic values times a scaling function that depends on the internuclear distances. The spin-valence diabatic potential energy matrix is suitable for semiclassical dynamics simulations. Diagonalization of this matrix produces the spin-coupled adiabatic energies. For the sake of illustration, diabatic potential energy matrices are constructed along bond-fission coordinates for the HBr and the BrCH(2)Cl molecules. Comparison of the spin-coupled adiabatic energies obtained from the spin-valence diabatics with those obtained by ab initio calculations with geometry-dependent spin-orbit matrix elements shows that the new method is sufficiently accurate for practical purposes. The method formulated here should be most useful for systems with a large number of atoms, especially heavy atoms, and/or a large number of spin-coupled electronic states.     

Ab initio adiabatic and diabatic energies and dipole moments of the CaH+ molecular ion

The diabatic and adiabatic potential-energy curves and permanent and transition dipole moments of the highly excited states of the CaH(+) molecular ion have been computed as a function of the internuclear distance R for a large and dense grid varying from 2.5 to 240 au. The adiabatic results are determined by an ab initio approach involving a nonempirical pseudopotential for the Ca core, operatorial core-valence correlation, and full valence configuration interaction. The molecule is thus treated as a two-electron system. The diabatic potential energy curves have been calculated using an effective metric combined to the effective Hamiltonian theory. The diabatic potential-energy curves and their permanent dipole moments for the (1)∑(+) symmetry are examined and corroborate the high imprint of the ionic state in the adiabatic representation. Taking the benefit of the diabatization approach, correction of hydrogen electron affinity was taken into account leading to improved results for the adiabatic potentials but also the permanent and transition electric dipole moments.     

Jahn-Teller effect in van der Waals complexes; Ar-C6H6 + and Ar-C6D6 +

The two asymptotically degenerate potential energy surfaces of argon interacting with the X (2)E(1g) ground state benzene(+) cation were calculated ab initio from the interaction energy of the neutral Ar-benzene complex given by Koch et al. [J. Chem. Phys. 111, 198 (1999)] and the difference of the geometry-dependent ionization energies of the complex and the benzene monomer computed by the outer valence Green's function method. Coinciding minima in the two potential surfaces of the ionic complex occur for Ar on the C(6v) symmetry axis of benzene(+) (the z axis) at z(e)=3.506 A. The binding energy D(e) of 520 cm(-1) is only 34% larger than the value for the neutral Ar-benzene complex. The higher one of the two surfaces is similar in shape to the neutral Ar-benzene potential, the lower potential is much flatter in the (x,y) bend direction. Nonadiabatic (Jahn-Teller) coupling was taken into account by transformation of the two adiabatic potentials to a two-by-two matrix of diabatic potentials. This transformation is based on the assumption that the adiabatic states of the Ar-benzene(+) complex geometrically follow the Ar atom. Ab initio calculations of the nonadiabatic coupling matrix element between the adiabatic states with the two-state-averaged CAS-SCF(5,6) method confirmed the validity of this assumption. The bound vibronic states of both Ar-C(6)H(6) (+) and Ar-C(6)D(6) (+) were computed with this two-state diabatic model in a basis of three-dimensional harmonic oscillator functions for the van der Waals modes. The binding energy D(0)=480 cm(-1) of the perdeuterated complex agrees well with the experimental upper bound of 485 cm(-1). The ground and excited vibronic levels and wave functions were used, with a simple model dipole function, to generate a theoretical far-infrared spectrum. Strong absorption lines were found at 10.1 cm(-1) (bend) and 47.9 cm(-1) (stretch) that agree well with measurements. The unusually low bend frequency is related to the flatness of the lower adiabatic potential in the (x,y) direction. The van der Waals bend mode of e(1) symmetry is quadratically Jahn-Teller active and shows a large splitting, with vibronic levels of A(1), E(2), and A(2) symmetry at 1.3, 10.1, and 50.2 cm(-1). The level at 1.3 cm(-1) leads to a strong absorption line as well, which could not be measured because it is too close to the monomer line. The level at 50.2 cm(-1) gives rise to weaker absorption. Several other weak lines in the frequency range of 10 to 60 cm(-1) were found.     

Allene and pentatetraene cations as models for intramolecular charge transfer: vibronic coupling Hamiltonian and conical intersections

We consider the vibronic coupling effects involving cationic states with degenerate components that can be represented as charge localized at either end of the short cumulene molecules allene and pentatetraene. Our aim is to simulate dynamically the charge transfer process when one component is artificially depopulated. We model the Jahn-Teller vibronic interaction within these states as well as their pseudo-Jahn-Teller coupling with some neighboring states. For the manifold of these states, we have calculated cross sections of the ab initio adiabatic potential energy surfaces along all nuclear degrees of freedom, including points at large distances from the equilibrium to increase the physical significance of our model. Ab initio calculations for the cationic states of allene and pentatetraene were based on the fourth-order M?ller-Plesset method and the outer valence Green's function method. In some cases we had to go beyond this method and use the more involved third-order algebraic diagrammatic construction method to include intersections with satellite states. The parameters for a five-state, all-mode diabatic vibronic coupling model Hamiltonian were least-square fitted to these potentials. The coupling parameters for the diabatic model Hamiltonian are such that, in comparison to allene, an enhanced preference for indirect charge transfer is predicted for pentatetraene.     

Ab initio Ehrenfest dynamics

We present an ab initio direct Ehrenfest dynamics scheme using a three time-step integrator. The three different time steps are implemented with nuclear velocity Verlet, nuclear-position-coupled midpoint Fock integrator, and time-dependent Hartree-Fock with a modified midpoint and unitary transformation algorithm. The computational cost of the ab initio direct Ehrenfest dynamics presented here is found to be only a factor of 2-4 larger than that of Born-Oppenheimer (BO) dynamics. As an example, we compute the vibration of the NaCl molecule and the intramolecular torsional motion of H2C=NH2+ by Ehrenfest dynamics compared with BO dynamics. For the vibration of NaCl with an initial kinetic energy of 1.16 eV, Ehrenfest dynamics converges to BO dynamics with the same vibrational frequency. The intramolecular rotation of H2C=NH2+ produces significant electronic excitation in the Ehrenfest trajectory. The amount of nonadiabaticity, suggested by the amplitude of the coherent progression of the excited and ground electronic states, is observed to be directly related to the strength of the electron-nuclear coupling. Such nonadiabaticity is seen to have a significant effect on the dynamics compared with the adiabatic approximation.     

On the peculiarities of the diabatic framework: new insight

In this article we consider the electronic diabatic presentation of a two-state system with the aim of earning insight regarding the distribution of conical intersections in a given region. In this process we revealed explicit relationship between the diabatic potentials and the locations of conical intersections. The study is accompanied with numerical examples as worked out for a model and ab initio potential energy surfaces of the Na+H2 system.