Cumulant reconstruction of the three-electron reduced density matrix in the anti-Hermitian contracted Schrödinger equation

DePrince III, A. Eugene; Mazziotti, David A.
September 2007
Journal of Chemical Physics;9/14/2007, Vol. 127 Issue 10, p104104
Academic Journal
Differing perspectives on the accuracy of three-electron reduced-density-matrix (3-RDM) reconstruction in nonminimal basis sets exist in the literature. This paper demonstrates the accuracy of cumulant-based reconstructions, developed by Valdemoro (V) [F. Colmenero et al., Phys. Rev. A 47, 971 (1993)], Nakatsuji and Yasuda (NY) [Phys. Rev. Lett. 76, 1039 (1996)], Mazziotti (M) [Phys. Rev. A 60, 3618 (1999)], and Valdemoro–Tel–Pérez–Romero (VTP) [Many-electron Densities and Density Matrices, edited by J. Cioslowski (Kluwer, Boston, 2000)]. Computationally, we extend previous investigations to study a variety of molecules, including LiH, HF, NH3, H2O, and N2, in Slater-type, double-zeta, and polarized double-zeta basis sets at both equilibrium and nonequilibrium geometries. The reconstructed 3-RDMs, compared with 3-RDMs from full configuration interaction, demonstrate in nonminimal basis sets the accuracy of the first-order expansion (V) as well as the important role of the second-order corrections (NY, M, and VTP). Calculations at nonequilibrium geometries further show that cumulant functionals can reconstruct the 3-RDM from a multireferenced 2-RDM with reasonable accuracy, which is relevant to recent multireferenced formulations of the anti-Hermitian contracted Schrödinger equation (ACSE) and canonical diagonalization. Theoretically, we perform a detailed perturbative analysis of the M functional to identify its second-order components. With these second-order components we connect the M, NY, and VTP reconstructions for the first time by deriving both the NY and VTP functionals from the M functional. Finally, these 3-RDM reconstructions are employed within the ACSE [D. Mazziotti, Phys. Rev. Lett. 97, 143002 (2006)] to compute ground-state energies which are compared with the energies from the contracted Schrödinger equation and several wave function methods.


Related Articles

  • Schrödinger equation in terms of linearly averaged position moments. Koga, Toshikatsu // Journal of Chemical Physics;10/15/1989, Vol. 91 Issue 8, p4774 

    Linearly averaged position moment (LAPM) [Rk11Rk22··· Rk3N3N] is defined as the linear (not quadratic) average of the position moment operator Rk11Rk22··· Rk3N3N over the N-electron wave function Ψ({Rj}), where {Rj} are 3N Cartesian coordinates of electrons and {kj} are...

  • Photoexcited conversion of gauche-1,3-butadiene to bicyclobutane via a conical intersection: Energies and reduced density matrices from the anti-Hermitian contracted Schrödinger equation. Snyder, James W.; Mazziotti, David A. // Journal of Chemical Physics;7/14/2011, Vol. 135 Issue 2, p024107 

    The photoexcited reaction pathway from gauche-1,3-butadiene to bicyclobutane via a conical intersection is analyzed through the direct calculation of two-electron reduced density matrices (2-RDMs) from solutions to the anti-Hermitian contracted Schrödinger equation (ACSE). The study utilizes...

  • Ultrashort metal–metal distances and extreme bond orders. Wagner, Frank R.; Noor, Awal; Kempe, Rhett // Nature Chemistry;Oct2009, Vol. 1 Issue 7, p529 

    Chemical bonding is at the very heart of chemistry. Although main-group-element E–E′ bond orders range up to triple bonds, higher formal bond orders are known between transition metals. Here we review recent developments related to the synthesis of formally quintuply bonded...

  • Variational method for solving the contracted Schro¨dinger equation through a projection of the N-particle power method onto the two-particle space. Mazziotti, David A. // Journal of Chemical Physics;1/22/2002, Vol. 116 Issue 4, p1239 

    The power method for solving N-particle eigenvalue equations is contracted onto the two-particle space to produce a reduced “variational” method for solving the contracted Schro¨dinger equation (CSE), also known as the density equation. In contrast to the methods which solve a...

  • Wigner’s (2 n + 1) rule for nonlinear Schrödinger equations. Ángyán, János G. // Journal of Mathematical Chemistry;Jul2009, Vol. 46 Issue 1, p1 

    Wigner’s (2 n + 1) rule is proven for general total energy functionals, constructed from the expectation value of a linear Hamiltonian and a nonlinear functional of the electron density or of the first order density matrix. Such functionals are common in independent particle models, like...

  • Vibrational predissociation of methylnitrite using phase-locked ultrashort laser pulses. Dateo, Christopher E.; Metiu, Horia // Journal of Chemical Physics;2/1/1993, Vol. 98 Issue 3, p1957 

    We solve numerically the time-dependent Schrödinger equation to study the behavior of a molecule interacting with two phase-locked ultrashort laser pulses. The two-dimensional model used in the calculations mimics the properties of the CH3ONO molecule. The two pulses are identical except for...

  • The intrapair electron correlation in natural orbital functional theory. Piris, M.; Matxain, J. M.; Lopez, X. // Journal of Chemical Physics;12/21/2013, Vol. 139 Issue 23, preceding p1 

    A previously proposed [M. Piris, X. Lopez, F. Ruipérez, J. M. Matxain, and J. M. Ugalde, J. Chem. Phys. 134, 164102 (2011)] formulation of the two-particle cumulant, based on an orbital-pairing scheme, is extended here for including more than two natural orbitals. This new approximation is...

  • On the h→0 limit of the Schrödinger equation. Ling, Song // Journal of Chemical Physics;5/15/1992, Vol. 96 Issue 10, p7869 

    Contrary to the previous belief, it is found that the time dependent Schrödinger equation is not always reduced to the Hamilton–Jacobi equation in the h→0 limit. An explicit example is given.

  • Nonlinear Schrödinger equations and the separation property. Goldin, Gerald A.; Svetlichny, George // Journal of Mathematical Physics;Jul94, Vol. 35 Issue 7, p3322 

    Hierarchies of nonlinear Schrödinger equations were investigated for multiparticle systems, satisfying the separation property, i.e., where product wave functions evolve by the separate evolution of each factor. Such a hierarchy defines a nonlinear derivation on tensor products of the...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics