Adiabatic potentials of linear three-site mixed-valence systems

We propose a method for calculating adiabatic potentials of linear or quasi-linear threesite systems of mixed valence. The method is based on transformation of the potential energy matrix to two effective interacting modes. The adiabatic potential is calculated and its extrema are found. The electronic wave functions at the minima of the adiabatic potential correspond to complex charge distributions in a cluster. The barrier to transition of the extra electron between the terminal atoms is lower than for the transition from the terminal atom to the central atom.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 25, No. 2, pp. 129–134, March–April, 1989.     

Dressed adiabatic and diabatic potentials to study conical intersections for F + H2

We follow a suggestion by Lipoff and Herschbach [Mol. Phys. 108, 1133 (2010)] and compare dressed and bare adiabatic potentials to get insight regarding the low-energy dynamics (e.g., cold reaction) taking place in molecular systems. In this particular case, we are interested to study the effect of conical intersections (ci) on the interacting atoms. For this purpose, we consider vibrational dressed adiabatic and vibrational dressed diabatic potentials in the entrance channel of reactive systems. According to our study, the most one should expect, in case of F + H(2), is a mild effect of the (1, 2) ci on its reactive/exchange process--an outcome also supported by experiment. This happens although the corresponding dressed and bare potential barriers (and the corresponding van der Waals potential wells) differ significantly from each other.     

Intracule densities in the strong-interaction limit of density functional theory

The correlation energy in density functional theory can be expressed exactly in terms of the change in the probability of finding two electrons at a given distance r(12) (intracule density) when the electron-electron interaction is multiplied by a real parameter lambda varying between 0 (Kohn-Sham system) and 1 (physical system). In this process, usually called adiabatic connection, the one-electron density is (ideally) kept fixed by a suitable local one-body potential. While an accurate intracule density of the physical system can only be obtained from expensive wavefunction-based calculations, being able to construct good models starting from Kohn-Sham ingredients would highly improve the accuracy of density functional calculations. To this purpose, we investigate the intracule density in the lambda --> infinity limit of the adiabatic connection. This strong-interaction limit of density functional theory turns out to be, like the opposite non-interacting Kohn-Sham limit, mathematically simple and can be entirely constructed from the knowledge of the one-electron density. We develop here the theoretical framework and, using accurate correlated one-electron densities, we calculate the intracule densities in the strong interaction limit for few atoms. Comparison of our results with the corresponding Kohn-Sham and physical quantities provides useful hints for building approximate intracule densities along the adiabatic connection of density functional theory.     

Range-dependent adiabatic connections

Recently, we have implemented a scheme for the calculation of the adiabatic connection linking the Kohn-Sham system to the physical, interacting system. This scheme uses a generalized Lieb functional, in which the electronic interaction strength is varied in a simple linear fashion, keeping the potential or the density fixed in the process. In the present work, we generalize this scheme further to accommodate arbitrary two-electron operators, allowing the calculation of adiabatic connections following alternative paths as outlined by Yang [J. Chem. Phys. 109, 10107 (1998)]. Specifically, we examine the error-function and Gaussian-attenuated error-function adiabatic connections. It is shown that while the error-function connection displays some promising features, making it amenable to the possible development of new exchange-correlation functionals by modeling the adiabatic connection integrand, the Gaussian-attenuated error-function connection is less promising. We explore the high-density and strong static correlation regimes for two-electron systems. Implications of this work for the utility of range-separated schemes are discussed.     

Generator coordinate method in time-dependent density-functional theory: memory made simple

The generator coordinate (GC) method is a variational approach to the quantum many-body problem in which interacting many-body wave functions are constructed as superpositions of (generally nonorthogonal) eigenstates of auxiliary Hamiltonians containing a deformation parameter. This paper presents a time-dependent extension of the GC method as a new approach to improve existing approximations of the exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT). The time-dependent GC method is shown to be a conceptually and computationally simple tool to build memory effects into any existing adiabatic XC potential. As an illustration, the method is applied to driven parametric oscillations of two interacting electrons in a harmonic potential (Hooke's atom). It is demonstrated that a proper choice of time-dependent generator coordinates in conjunction with the adiabatic local-density approximation reproduces the exact linear and nonlinear two-electron dynamics quite accurately, including features associated with double excitations that cannot be captured by TDDFT in the adiabatic approximation.     

Angle dependent total cross sections and the optical theorem

Cross sections are either represented by generalized asymptotical partial wave expansions or obtained as a spherical average of an appropriate differential cross section. In these cases it is usually assumed that the total scattering cross section, as a property of a scattering object, does not depend on the incident angles. This viewpoint is supported by common knowledge in connection with low energy scattering. However this unconscious belief is not always correct. In the present paper we will show that a non-spherical scatterer may exhibit strong dependence on the incident direction. To do this we will represent the scattering data of the most general potential, separable in ellipsoidal coordinates, in perturbed ellipsoidal (Lamé) wave functions. These functions arise when variables in the Schr?dinger equation are separated in an ellipsoidal coordinate system. The Lamé wave functions are analogous to spherical- and Bessel functions in the case of spherical symmetry. We will expand the total scattering cross section and derive the optical theorem explicitly demonstrating the incident angle dependence for such a class of potentials. As an illustration we will present and display some calculations of the total cross section versus incident direction. Unexpected behavior will be discussed and explained. We also use results from classical acoustic scattering by a triaxial ellipsoid. The general character of the ellipsoidal coordinate system is emphasized.     

Comparison of the localization of an electron as determined by the two-particle distribution function and by the single-particle sharing index

A comparison of the measure of the delocalization of a particle based on the two-particle distribution function and that based on the single-particle density matrix is made using a simple set of wave functions which span states ranging from single determinant ground and doubly excited states through states mimicking correlated states and which include the singly excited state for electrons and for bosons replacing electrons in H2. The comparison further includes an analysis of the application of the measures to a classical ideal gas and a compressible fluid. It is found that the values of the integrated atom-atom measures agree for a range of wave functions involving combinations of the two single determinant (and equivalent Bose) wave functions but disagree for a different range of these wave functions and for the singly excited wave functions. Aside from the single determinant (and equivalent Bose) wave functions, the two sets of point-point measures that underlie the integrated measures all differ. For the sets of wave functions considered, the values of the measures are identical for electrons and bosons. When applied to a closed classical ideal gas and to a closed compressible fluid, the delocalization measure based on the two-particle distribution has a residual long range term, whereas the sharing index in the classical limit gives a completely localized particle. In general, the two measures describe different aspects of the behavior of the particles. The measures based on the two-particle distribution function give only two-particle properties and the single-particle density, and the sharing quantities give only single-particle properties. The latter includes, however, the quantitative measures of the delocalization of a single particle, the point-point sharing index and the sharing amplitude.     

Hamiltonian matrix and reduced density matrix construction with nonlinear wave functions

An efficient procedure to compute Hamiltonian matrix elements and reduced one- and two-particle density matrices for electronic wave functions using a new graphical-based nonlinear expansion form is presented. This method is based on spin eigenfunctions using the graphical unitary group approach (GUGA), and the wave function is expanded in a basis of product functions (each of which is equivalent to some linear combination of all of the configuration state functions), allowing application to closed- and open-shell systems and to ground and excited electronic states. In general, the effort required to construct an individual Hamiltonian matrix element between two product basis functions H(MN) = M|H|N scales as theta (beta n4) for a wave function expanded in n molecular orbitals. The prefactor beta itself scales between N0 and N2, for N electrons, depending on the complexity of the underlying Shavitt graph. Timings with our initial implementation of this method are very promising. Wave function expansions that are orders of magnitude larger than can be treated with traditional CI methods require only modest effort with our new method.     

Method of Calculating Two-Center Two-Particle Matrix Elements from a Screened Coulomb Potential in an Exponential AO Basis

An analytical method is suggested for calculating matrix elements from an exponentially screened two-particle Coulomb potential in a basis set of exponential functions specified on different centers in space. The two-particle potentials of this type may be used to approximate the short-acting parts of atom–atomic interactions in molecules in quantum-chemical structural calculations.     

Constrained adiabatic trajectory method: a global integrator for explicitly time-dependent Hamiltonians

The constrained adiabatic trajectory method (CATM) is reexamined as an integrator for the Schro?dinger equation. An initial discussion places the CATM in the context of the different integrators used in the literature for time-independent or explicitly time-dependent Hamiltonians. The emphasis is put on adiabatic processes and within this adiabatic framework the interdependence between the CATM, the wave operator, the Floquet, and the (t, t') theories is presented in detail. Two points are then more particularly analyzed and illustrated by a numerical calculation describing the H(2)(+) ion submitted to a laser pulse. The first point is the ability of the CATM to dilate the Hamiltonian spectrum and thus to make the perturbative treatment of the equations defining the wave function possible, possibly by using a Krylov subspace approach as a complement. The second point is the ability of the CATM to handle extremely complex time-dependencies, such as those which appear when interaction representations are used to integrate the system.     

Effects of higher order Jahn-Teller coupling on the nuclear dynamics

In this paper effects of higher order Jahn-Teller coupling terms on the nonadiabatic dynamics are studied. Of particular interest is the case when the potential energy surfaces of the degenerate state show pronounced anharmonicity. In order to demonstrate the effects a two-dimensional E multiply sign in circle e Jahn-Teller model system is treated which is based on the e(') stretching vibration of the photoactive (2)E(') state of NO(3) as a realistic example. The sixth order E multiply sign in circle e Jahn-Teller Hamiltonian is derived in the diabatic representation which is valid for any system with a C(3) rotation axis. This diabatization scheme is compared to lower-order Jahn-Teller Hamiltonians and to symmetry adapted as well as ad hoc approximations. Lower-order potentials result in pronounced quantitative and qualitative differences in the dynamics, including differences in the evolution of mean values, the autocorrelation functions (and thus the corresponding spectra), and the electronic population evolution. In the particular example treated, the results of fourth and fifth order potentials are very similar to the sixth order reference system. In contrast, the approximate sixth order Hamiltonians, though the corresponding adiabatic surfaces seem to be nearly identical, results in pronounced differences. The possible consequences for the dynamics of realistic systems with higher dimensionality are briefly discussed.     

Adiabatic passage by light-induced potentials in polyatomic molecules

In this paper we study the first application of adiabatic passage by light-induced potentials in polyatomic molecules. We analyze the effects of increasing the dimensionality of the system on the adiabatic requirements of the method and the role of intramolecular coupling among the vibrational modes. By using a model of two-dimensional displaced harmonic oscillators with or without rotation of the normal mode axis of the excited states (Duschinsky effect) we find that (1) it is possible to selectively transfer the vibrational population by adiabatic elongation of the bonds, (2) the adiabatic demands depend mainly on the energy barrier between the ground and excited electronic configurations, and not on the dimension of the system, (3) in the presence of intramolecular couplings the selective transfer can be achieved but at the cost of increasing the duration and/or the intensity of the pulses, which are needed to overcome small avoided crossings, and (4) the problem of selectivity becomes more important as the vibrational energy of the initial wave function increases.     

An ab initio study of the fcc and hcp structures of helium

The hexagonal close packed (hcp) and face centered cubic (fcc) structures of helium are studied by using a new ab initio computational model for large complexes comprising small subsystems. The new model is formulated within the framework of the energy incremental scheme. In the calculation of intra- and intersystem energies, model systems are introduced. To each subsystem associated is a set of partner subsystems defined by a vicinity criterion. In the independent calculations of intra- and intersystem energies, the calculations are performed on model subsystems defined by the subsystems considered and their partner subsystems. A small and a large basis set are associated with each subsystem. For partner subsystems in a model system, the small basis set is adopted. By introducing a particular decomposition scheme, the intermolecular potential is written as a sum of effective one-body potentials. The binding energy per atom in an infinite crystal of atoms is the negative value of this one-body potential. The one-body potentials for hcp and fcc structures are calculated for the following nearest neighbor distances (d0): 4.6, 5.1, 5.4, 5.435, 5.5, 5.61, and 6.1 a.u. The equilibrium distance is 5.44 a.u. for both structures. The equilibrium dimer distance is 5.61 a.u. For the larger distances, i.e., d0 > 5.4 a.u., the difference of the effective one-body potentials for the two structures is less than 0.2 microE(h). However, the hcp structure has the lowest effective one-body potential for all the distances considered. For the smallest distance the difference in the effective one-body potential is 3.9 microE(h). Hence, for solid helium, i.e., helium under high pressure, the hcp structure is the preferred one. The error in the calculated effective one-body potential for the distance d0 = 5.61 a.u. is of the order of 1 microE(h) (approximately 0.5%).     

The properties of off-shell coulomb radial wavefunctions

Approximations to exact wave functions for the scattering of few-particle systems often involve components corresponding to the interaction of two of the particles “off the energy-shell”. Several examples arising in the collision of ions and photons with atoms are given. An expansion in partial waves leads to an off-shell radial wave function. The defining differential equation is solved here numerically with particular emphasis on the behaviour arising from two-body potentials of long-range Coulomb form. The transition to shell of the radial wave functions, Jost functions and solutions andT-matrix elements is discussed for both short-range and Coulomb potentials. It is shown that the approximation of a Coulomb potential by a shorter-range form involves little error when sufficiently far off the energy shell.     

Nonlocal Wigner-like correlation energy density functional: parametrization and tests on two-electron systems

Reparametrization of Wigner's correlation energy density functional yields a very close fit to the correlation energies of the helium isoelectronic sequence. However, a quite different reparametrization is required to obtain an equally close fit to the isoelectronic sequence of Hooke's atom. In an attempt to avoid having to reparametrize the functional for different choices of the one-body potential, we propose a parametrization that depends on global characteristics of the ground-state electron density as quantified by scale-invariant combinations of expectation values of local one-body operators. This should be viewed as an alternative to the density-gradient paradigm, allowing one to introduce the nonlocal dependence of the density functional on the density in a possibly more effective way. Encouraging results are obtained for two-electron systems with one-body potentials of the form r(zeta) with zeta=-12,+12,1, which span the range between the Coulomb potential (zeta=-1) and the Hooke potential (zeta=2).     

On Theoretical Study of a Molecular Dimer System

In this paper, the vibronic structure of a dimer system is studied both theoretically and numerically. To construct adiabatic potential surfaces and electronic and vibrational wave functions for a dimer system, the adiabatic approximation is applied to two identical molecules, each of which has two electronic states with one vibrational mode. In this scheme, the excitonic splitting results not only from the electronic coupling of two molecules, but also from the vibronic coupling in each molecule. By using the resulting wavefunctions and the corresponding energies, the absorption and fluorescence spectra are studied. The effect of temperature on these spectra is also studied.     

Effective local potentials for excited states

The constrained variational Hartree-Fock method for excited states of the same symmetry as the ground state [Chem. Phys. Lett. 287, 189 (1998)] is combined with the effective local potential (ELP) method [J. Chem. Phys. 125, 081104 (2006)] to generate Kohn-Sham-type exact-exchange potentials for singly excited states of many-electron systems. Illustrative examples include the three lowest (2)S states of the Li and Na atoms and the three lowest (3)S states of He and Be. For the systems studied, excited-state ELPs differ from the corresponding ground-state potentials in two respects: They are less negative and have small additional "bumps" in the outer electron region. The technique is general and can be used to approximate excited-state exchange-correlation potentials for other orbital-dependent functionals.     

Local effective potential theory: nonuniqueness of potential and wave function

In local effective potential energy theories such as the Hohenberg-Kohn-Sham density functional theory (HKS-DFT) and quantal density functional theory (Q-DFT), electronic systems in their ground or excited states are mapped to model systems of noninteracting fermions with equivalent density. From these models, the equivalent total energy and ionization potential are also obtained. This paper concerns (i) the nonuniqueness of the local effective potential energy function of the model system in the mapping from a nondegenerate ground state, (ii) the nonuniqueness of the local effective potential energy function in the mapping from a nondegenerate excited state, and (iii) in the mapping to a model system in an excited state, the nonuniqueness of the model system wave function. According to nondegenerate ground state HKS-DFT, there exists only one local effective potential energy function, obtained as the functional derivative of the unique ground state energy functional, that can generate the ground state density. Since the theorems of ground state HKS-DFT cannot be generalized to nondegenerate excited states, there could exist different local potential energy functions that generate the excited state density. The constrained-search version of HKS-DFT selects one of these functions as the functional derivative of a bidensity energy functional. In this paper, the authors show via Q-DFT that there exist an infinite number of local potential energy functions that can generate both the nondegenerate ground and excited state densities of an interacting system. This is accomplished by constructing model systems in configurations different from those of the interacting system. Further, they prove that the difference between the various potential energy functions lies solely in their correlation-kinetic contributions. The component of these functions due to the Pauli exclusion principle and Coulomb repulsion remains the same. The existence of the different potential energy functions as viewed from the perspective of Q-DFT reaffirms that there can be no equivalent to the ground state HKS-DFT theorems for excited states. Additionally, the lack of such theorems for excited states is attributable to correlation-kinetic effects. Finally, they show that in the mapping to a model system in an excited state, there is a nonuniqueness of the model system wave function. Different wave functions lead to the same density, each thereby satisfying the sole requirement of reproducing the interacting system density. Examples of the nonuniqueness of the potential energy functions for the mapping from both ground and excited states and the nonuniqueness of the wave function are provided for the exactly solvable Hooke's atom. The work of others is also discussed.