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

The Journal of Chemical Physics

Determination of reduced density matrices in the doubly occupied configuration interaction space: A Hellmann–Feynman theorem approach

Adán Garros

Published Oct 1, 2024• 1 citation
Publisher

Abstract

In this work, the Hellmann-Feynman theorem is extended within the doubly occupied configuration interaction space to enable practical calculations of reduced density matrices and expected values. This approach is straightforward, employing finite energy differences, yet remains reliable and accurate even with approximate energies from successive approximation methods. The method's validity is rigorously tested against the Richardson-Gaudin-Kitaev and reduced Bardeen-Cooper-Schrieffer models using approximate excitation energies procured from the Hermitian operator method within the same space, effectively proving the approach's reliability with median error rates for reduced density matrix calculations around 0.1%. These results highlight the procedure's potential as a practical tool for computing reduced density matrices and expected values, particularly valuable as an ad hoc method in scenarios where only system energies are easily available.

DOI
10.1063/5.0228431

Related papers

Determination of electronic excitation energies within the doubly occupied configuration interaction space by means of the Hermitian operator method

Adán Garros, Diego R. Alcoba, P. Capuzzi +4 more

In this work, we formulate the equations of motion corresponding to the Hermitian operator method in the framework of the doubly occupied configuration interaction space. The resulting algorithms turn out to be considerably simpler than the equations provided by that method in more conventional spaces, enabling the determination of excitation energies in N-electron systems under an affordable polynomial computational cost. The implementation of this technique only requires to know the elements of low-order reduced density matrices of an N-electron reference state, which can be obtained from any approximate method. We contrast our procedure against the reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev integrable models, pointing out the reliability of our proposal.

The Journal of Chemical Physics•Sep 25, 2023•1 citation

Variational determination of ground and excited-state two-electron reduced density matrices in the doubly occupied configuration space: A dispersion operator approach

Diego R. Alcoba, Ofelia B. Oña, Luis Laín +6 more

This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.

The Journal of Chemical Physics•Jun 10, 2021•12 citations

Benchmarking the Variational Reduced Density Matrix Theory in the Doubly Occupied Configuration Interaction Space with Integrable Pairing Models

Álvaro Rubio-García, Diego R. Alcoba, P. Capuzzi +1 more

The variational reduced density matrix theory has been recently applied with great success to models within the truncated doubly occupied configuration interaction space, which corresponds to the seniority zero subspace. Conservation of the seniority quantum number restricts the Hamiltonians to be based on the SU(2) algebra. Among them there is a whole family of exactly solvable Richardson-Gaudin pairing Hamiltonians. We benchmark the variational theory against two different exactly solvable models, the Richardson-Gaudin-Kitaev and the reduced BCS Hamiltonians. We obtain exact numerical results for the so-called [Formula: see text] N-representability conditions in both cases for systems that go from 10 to 100 particles. However, when random single-particle energies as appropriate for small superconducting grains are considered, the exactness is lost but still a high accuracy is obtained.

arXiv:1805.01322•Journal of Chemical Theory and Computation•Jun 15, 2018•15 citations

Influence of the N-representability conditions on the variational determination of two-electron reduced density matrices for ground and excited N-electron states in the doubly occupied configuration interaction space

Dario A. Corvalan, Ofelia B. Oña, Diego R. Alcoba +3 more

In this work, we analyze the effectiveness of different sets of well-known necessary N-representability conditions that are imposed on the variational determination of two-electron reduced density matrices corresponding to the ground and excited states of N-electron Hamiltonians within doubly occupied configuration interaction subspaces. The two-electron reduced density matrix elements are determined so that they are suitably related to those of higher-order reduced density matrices and other associated matrices, up to a given order, constituting positive semidefinite matrices. These conditions up to two-, three-, and four-order have been sequentially imposed on the variational calculation of the expectation values of the dispersion operator corresponding to several Hamiltonian models, obtaining the electronic spectra. These results have been contrasted with those arising from the full configuration interaction method, illuminating the sufficiency/insufficiency of the variational N-representability conditions applied.

The Journal of Chemical Physics•Nov 7, 2022•9 citations

Boson correlation energies and density matrices from reduced Hamiltonian interpolation

David A. Mazziotti, Dudley R. Herschbach

The ground-state energies of interacting bosons are computed beyond the mean-field approximation by a new method known as reduced Hamiltonian interpolation (RHI) [D. A. Mazziotti and D. R. Herschbach, Phys. Rev. Lett. 83, 5185 (1999)]. In the RHI, the N-particle Hamiltonian is represented through a sequence of p-particle expanded and reduced Hamiltonians that give upper and lower bounds on the true energy. Assimilating ideas from N representability and dimensional interpolation, the RHI technique interpolates over the number p of quasiparticles to determine the N-particle energy with close upper and lower bounds. With the Hellmann-Feynman theorem, we extend the RHI to compute the two-particle reduced density matrix (2RDM) as well as the energy. We examine both the computational advantages of RHI in comparison with traditional methods and the possibility of extending the RHI to treat fermion correlation. Applied to bosons with harmonic interactions as well as a two-level system, the RHI technique yields more than 99% of the correlation energy and an accurate correlation correction for the elements of the 2RDM.

Physical Review A•Sep 1, 2000•16 citations

Generalized spin in the variational determination of two-electron reduced density matrices within the doubly occupied configuration interaction framework

Tomás R. Ayala, Elías Ríos, Ofelia B. Oña +4 more

This work extends the variational determination of two-electron reduced density matrices to the generalized spin formulation within the doubly occupied configuration interaction treatment. We describe electronic states of N-electron systems by means of variational evaluation of the corresponding two-electron reduced density matrix elements, subject to determined N-representability conditions, arising from the restricted, unrestricted, and generalized spin approaches. The results are compared with those predicted by the counterpart treatments involving wave functions. We report potential energy curves for a series of hydrogen atom clusters with different geometries in the presence of an external uniform magnetic field. These results allow us to analyze and discuss crossings between potential energy curves involving transitions between electronic states of different spin multiplicities. The predictions provided by these methods have been contrasted with the results arising from the full configuration interaction procedure, confirming changes of ⟨Ŝ2⟩ and ⟨Ŝz⟩ expectation values.

The Journal of Chemical Physics•Sep 22, 2025•1 citation
In this work, the Hellmann-Feynman theorem is extended within the doubly occupied configuration interaction space to enable practical calculations of reduced density matrices and expected values. This approach is straightforward, employing finite energy differences, yet remains reliable and accurate even with approximate energies from successive approximation methods. The method's validity is rigorously tested against the Richardson-Gaudin-Kitaev and reduced Bardeen-Cooper-Schrieffer models using approximate excitation energies procured from the Hermitian operator method within the same space, effectively proving the approach's reliability with median error rates for reduced density matrix calculations around 0.1%. These results highlight the procedure's potential as a practical tool for computing reduced density matrices and expected values, particularly valuable as an ad hoc method in scenarios where only system energies are easily available.
In this work, we formulate the equations of motion corresponding to the Hermitian operator method in the framework of the doubly occupied configuration interaction space. The resulting algorithms turn out to be considerably simpler than the equations provided by that method in more conventional spaces, enabling the determination of excitation energies in N-electron systems under an affordable polynomial computational cost. The implementation of this technique only requires to know the elements of low-order reduced density matrices of an N-electron reference state, which can be obtained from any approximate method. We contrast our procedure against the reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev integrable models, pointing out the reliability of our proposal.
This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.
The variational reduced density matrix theory has been recently applied with great success to models within the truncated doubly occupied configuration interaction space, which corresponds to the seniority zero subspace. Conservation of the seniority quantum number restricts the Hamiltonians to be based on the SU(2) algebra. Among them there is a whole family of exactly solvable Richardson-Gaudin pairing Hamiltonians. We benchmark the variational theory against two different exactly solvable models, the Richardson-Gaudin-Kitaev and the reduced BCS Hamiltonians. We obtain exact numerical results for the so-called [Formula: see text] N-representability conditions in both cases for systems that go from 10 to 100 particles. However, when random single-particle energies as appropriate for small superconducting grains are considered, the exactness is lost but still a high accuracy is obtained.
In this work, we analyze the effectiveness of different sets of well-known necessary N-representability conditions that are imposed on the variational determination of two-electron reduced density matrices corresponding to the ground and excited states of N-electron Hamiltonians within doubly occupied configuration interaction subspaces. The two-electron reduced density matrix elements are determined so that they are suitably related to those of higher-order reduced density matrices and other associated matrices, up to a given order, constituting positive semidefinite matrices. These conditions up to two-, three-, and four-order have been sequentially imposed on the variational calculation of the expectation values of the dispersion operator corresponding to several Hamiltonian models, obtaining the electronic spectra. These results have been contrasted with those arising from the full configuration interaction method, illuminating the sufficiency/insufficiency of the variational N-representability conditions applied.
The ground-state energies of interacting bosons are computed beyond the mean-field approximation by a new method known as reduced Hamiltonian interpolation (RHI) [D. A. Mazziotti and D. R. Herschbach, Phys. Rev. Lett. 83, 5185 (1999)]. In the RHI, the N-particle Hamiltonian is represented through a sequence of p-particle expanded and reduced Hamiltonians that give upper and lower bounds on the true energy. Assimilating ideas from N representability and dimensional interpolation, the RHI technique interpolates over the number p of quasiparticles to determine the N-particle energy with close upper and lower bounds. With the Hellmann-Feynman theorem, we extend the RHI to compute the two-particle reduced density matrix (2RDM) as well as the energy. We examine both the computational advantages of RHI in comparison with traditional methods and the possibility of extending the RHI to treat fermion correlation. Applied to bosons with harmonic interactions as well as a two-level system, the RHI technique yields more than 99% of the correlation energy and an accurate correlation correction for the elements of the 2RDM.
This work extends the variational determination of two-electron reduced density matrices to the generalized spin formulation within the doubly occupied configuration interaction treatment. We describe electronic states of N-electron systems by means of variational evaluation of the corresponding two-electron reduced density matrix elements, subject to determined N-representability conditions, arising from the restricted, unrestricted, and generalized spin approaches. The results are compared with those predicted by the counterpart treatments involving wave functions. We report potential energy curves for a series of hydrogen atom clusters with different geometries in the presence of an external uniform magnetic field. These results allow us to analyze and discuss crossings between potential energy curves involving transitions between electronic states of different spin multiplicities. The predictions provided by these methods have been contrasted with the results arising from the full configuration interaction procedure, confirming changes of ⟨Ŝ2⟩ and ⟨Ŝz⟩ expectation values.