No abstract

Determining the energy and properties of an N-electron molecule through a two-electron variational optimization has been a dream for more than half a century. While optimizations, using two-electron reduced density matrices constrained to represent N electrons, have recently been achieved, the computational costs are prohibitive. In this report an efficient algorithm with an order-of-magnitude reduction in floating-point operations and memory usage is presented. Because the optimization occurs on the space of two electrons, this method automatically treats strong, multireference correlation. Application is made to N2 and H6 where the method yields consistent accuracy at all geometries.

We present an implementation of the time-dependent configuration-interaction singles (TDCIS) method for treating atomic strong-field processes. In order to absorb the photoelectron wave packet when it reaches the end of the spatial grid, we add to the exact nonrelativistic many-electron Hamiltonian a radial complex absorbing potential (CAP). We determine the orbitals for the TDCIS calculation by diagonalizing the sum of the Fock operator and the CAP using a flexible pseudospectral grid for the radial degree of freedom and spherical harmonics for the angular degrees of freedom. The CAP is chosen such that the occupied orbitals in the Hartree-Fock ground state remain unaffected. Within TDCIS, the many-electron wave packet is expanded in terms of the Hartree-Fock ground state and its single excitations. The virtual orbitals satisfy nonstandard orthogonality relations, which must be taken into consideration in the calculation of the dipole and Coulomb matrix elements required for the TDCIS equations of motion. We employ a stable propagation scheme derived by second-order finite differencing of the TDCIS equations of motion in the interaction picture and subsequent transformation to the Schrödinger picture. Using the TDCIS wave packet, we calculate the expectation value of the dipole acceleration and the reduced density matrix of the residual ion. The technique implemented will allow one to study electronic channel-coupling effects in strong-field processes.

Molecular systems in chemistry often have wave functions with substantial contributions from two-or-more electronic configurations. Because traditional complete-active-space self-consistent-field (CASSCF) methods scale exponentially with the number N of active electrons, their applicability is limited to small active spaces. In this paper we develop an active-space variational two-electron reduced-density-matrix (2-RDM) method in which the expensive diagonalization is replaced by a variational 2-RDM calculation where the 2-RDM is constrained by approximate N-representability conditions. Optimization of the constrained 2-RDM is accomplished by large-scale semidefinite programming [Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. Because the computational cost of the active-space 2-RDM method scales polynomially as r(a)(6) where r(a) is the number of active orbitals, the method can be applied to treat active spaces that are too large for conventional CASSCF. The active-space 2-RDM method performs two steps: (i) variational calculation of the 2-RDM in the active space and (ii) optimization of the active orbitals by Jacobi rotations. For large basis sets this two-step 2-RDM method is more efficient than the one-step, low-rank variational 2-RDM method [Gidofalvi and Mazziotti, J. Chem. Phys. 127, 244105 (2007)]. Applications are made to HF, H(2)O, and N(2) as well as n-acene chains for n=2-8. When n>4, the acenes cannot be treated by conventional CASSCF methods; for example, when n=8, CASSCF requires optimization over approximately 1.47x10(17) configuration state functions. The natural occupation numbers of the n-acenes show the emergence of bi- and polyradical character with increasing chain length.

We present a constructive solution to the N-representability problem: a full characterization of the conditions for constraining the two-electron reduced density matrix to represent an N-electron density matrix. Previously known conditions, while rigorous, were incomplete. Here, we derive a hierarchy of constraints built upon (i) the bipolar theorem and (ii) tensor decompositions of model Hamiltonians. Existing conditions D, Q, G, T1, and T2, known classical conditions, and new conditions appear naturally. Subsets of the conditions are amenable to polynomial-time computations of strongly correlated systems.