Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Chemical Physics Letters

Analytical nuclear derivatives for the parametric two-electron reduced density matrix method

Andrew J. Valentine, David A. Mazziotti

Published Jul 27, 2017• 5 citations
Publisher
DOI
10.1016/j.cplett.2017.07.064

Related papers

Methods for directly determining the two-electron reduced density matrix with applications to quantum phase transitions and chemical reactions

Christine A. Schwerdtfeger

•Jan 1, 2012•1 citation

The Treatment of Diradical Systems and Molecular Optimizations with Two-Electron Reduced Density Matrices

Andrew J. Valentine

The electronic Hamiltonian contains only pairwise interactions, allowing the energy of an electronic system to be expressed in terms of the two-electron reduced-density-matrix (2-RDM) in lieu of the many-electron wavefunction. The variable space for the exact N-electron wavefunction scales exponentially with the size of the system, while the 2-RDM is polynomial in scale. By using the 2-RDM as the primary variable in electronic structure calculations, it may be possible to obtain very accurate energies at a much more favorable scaling than wavefunction methods. In this thesis, we will use two existing 2-RDM methods to treat electronic systems. First, we will apply the active-space variational 2-RDM method, which directly minimizes the energy with respect to the 2-RDM, to a cadmium telluride polymer that was recently used to greatly enhance the conductivity of CdTe quantum dots. We find that this polymer is very highly correlated despite a deceptively simple structure. We will then turn to the parametric 2-RDM method (p2-RDM), which parameterizes the 2-RDM in terms of a truncated configuration interaction ansatz, but which includes additional flexibility in order to be size-extensive. We apply p2-RDM to the study of the olympicene molecule, which features both fully aromatic and diradical isomers. The parametric 2-RDM method predicts that all isomers are stable to dissociation, in contrast to coupled cluster methods which do not predict stable diradical states. We then present analytical nuclear gradients for p2-RDM, which greatly decrease the number of calculations required to perform geometry optimizations. We apply these gradients to the study of trans-polyacetylene, for which p2-RDM, unlike many wavefunction methods, is able to predict a bond length alternation (BLA) to within experimental values. Lastly, as a single-reference method, p2-RDM may encounter numerical difficulties when the reference wavefunction is of particularly poor quality. We propose a modification to the parameterization that may render the method more generally robust.

•Jan 1, 2017•0 citations

Density matrix methods in nuclear reactions

R. Parthasarathy

Institute of Mathematical Sciences eBooks•Jan 1, 1979•0 citations

Reduced density matrix approach to calculation of electronic structure

Bastiaan J. Braams, Shidong Jiang, Madhu V. Nayakkankuppam +2 more

APS•Nov 1, 1998•0 citations

Analytical gradients of variational reduced-density-matrix and wavefunction-based methods from an overlap-reweighted semidefinite program

Anthony W. Schlimgen, David A. Mazziotti

Analytical gradients of variational two-electron reduced-density matrix (2-RDM) methods are derived by transforming the atomic-orbital reduced-density matrices to remove the dependence of the N-representability conditions on the orbital-overlap matrix. The transformation, performed through a Cholesky decomposition of the geminal-overlap matrix, generates a Hellmann-Feynman-like expression for the gradient that only depends on the derivative of the transformed reduced Hamiltonian matrix. The formulation is applicable not only to the variational 2-RDM method but also to variational wavefunction methods like the full configuration interaction and complete active-space self-consistent-field. To illustrate, we apply the analytical gradients to perform geometry optimizations on several transition metal complexes, octahedral and trigonal prismatic CrF6 as well as the (ethylene-1,2-dithiolato)nickel, or Ni(edt)2, complex.

The Journal of Chemical Physics•Oct 28, 2018•7 citations

On the relationship between parametric two-electron reduced-density-matrix methods and the coupled electron pair approximation

A. Eugene DePrince, David A. Mazziotti

Parametric two-electron reduced-density-matrix (p-2RDM) methods have enjoyed much success in recent years; the methods have been shown to exhibit accuracies greater than coupled cluster with single and double substitutions (CCSD) for both closed- and open-shell ground-state energies, properties, geometric parameters, and harmonic frequencies. The class of methods is herein discussed within the context of the coupled electron pair approximation (CEPA), and several CEPA-like topological factors are presented for use within the p-2RDM framework. The resulting p-2RDM/n methods can be viewed as a density-based generalization of CEPA/n family that are numerically very similar to traditional CEPA methodologies. We cite the important distinction that the obtained energies represent stationary points, facilitating the efficient evaluation of properties and geometric derivatives. The p-2RDM/n formalism is generalized for an equal treatment of exclusion-principle-violating (EPV) diagrams that occur in the occupied and virtual spaces. One of these general topological factors is shown to be identical to that proposed by Kollmar [C. Kollmar, J. Chem. Phys. 125, 084108 (2006)], derived in an effort to approximately enforce the D, Q, and G conditions for N-representability in his size-extensive density matrix functional.

arXiv:1009.0047•Aug 31, 2010•0 citations
The electronic Hamiltonian contains only pairwise interactions, allowing the energy of an electronic system to be expressed in terms of the two-electron reduced-density-matrix (2-RDM) in lieu of the many-electron wavefunction. The variable space for the exact N-electron wavefunction scales exponentially with the size of the system, while the 2-RDM is polynomial in scale. By using the 2-RDM as the primary variable in electronic structure calculations, it may be possible to obtain very accurate energies at a much more favorable scaling than wavefunction methods. In this thesis, we will use two existing 2-RDM methods to treat electronic systems. First, we will apply the active-space variational 2-RDM method, which directly minimizes the energy with respect to the 2-RDM, to a cadmium telluride polymer that was recently used to greatly enhance the conductivity of CdTe quantum dots. We find that this polymer is very highly correlated despite a deceptively simple structure. We will then turn to the parametric 2-RDM method (p2-RDM), which parameterizes the 2-RDM in terms of a truncated configuration interaction ansatz, but which includes additional flexibility in order to be size-extensive. We apply p2-RDM to the study of the olympicene molecule, which features both fully aromatic and diradical isomers. The parametric 2-RDM method predicts that all isomers are stable to dissociation, in contrast to coupled cluster methods which do not predict stable diradical states. We then present analytical nuclear gradients for p2-RDM, which greatly decrease the number of calculations required to perform geometry optimizations. We apply these gradients to the study of trans-polyacetylene, for which p2-RDM, unlike many wavefunction methods, is able to predict a bond length alternation (BLA) to within experimental values. Lastly, as a single-reference method, p2-RDM may encounter numerical difficulties when the reference wavefunction is of particularly poor quality. We propose a modification to the parameterization that may render the method more generally robust.
Analytical gradients of variational two-electron reduced-density matrix (2-RDM) methods are derived by transforming the atomic-orbital reduced-density matrices to remove the dependence of the N-representability conditions on the orbital-overlap matrix. The transformation, performed through a Cholesky decomposition of the geminal-overlap matrix, generates a Hellmann-Feynman-like expression for the gradient that only depends on the derivative of the transformed reduced Hamiltonian matrix. The formulation is applicable not only to the variational 2-RDM method but also to variational wavefunction methods like the full configuration interaction and complete active-space self-consistent-field. To illustrate, we apply the analytical gradients to perform geometry optimizations on several transition metal complexes, octahedral and trigonal prismatic CrF6 as well as the (ethylene-1,2-dithiolato)nickel, or Ni(edt)2, complex.
Parametric two-electron reduced-density-matrix (p-2RDM) methods have enjoyed much success in recent years; the methods have been shown to exhibit accuracies greater than coupled cluster with single and double substitutions (CCSD) for both closed- and open-shell ground-state energies, properties, geometric parameters, and harmonic frequencies. The class of methods is herein discussed within the context of the coupled electron pair approximation (CEPA), and several CEPA-like topological factors are presented for use within the p-2RDM framework. The resulting p-2RDM/n methods can be viewed as a density-based generalization of CEPA/n family that are numerically very similar to traditional CEPA methodologies. We cite the important distinction that the obtained energies represent stationary points, facilitating the efficient evaluation of properties and geometric derivatives. The p-2RDM/n formalism is generalized for an equal treatment of exclusion-principle-violating (EPV) diagrams that occur in the occupied and virtual spaces. One of these general topological factors is shown to be identical to that proposed by Kollmar [C. Kollmar, J. Chem. Phys. 125, 084108 (2006)], derived in an effort to approximately enforce the D, Q, and G conditions for N-representability in his size-extensive density matrix functional.