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 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, Luis Laín, Alicia Torre, Ofelia B. Oña, J. Dukelsky

Published Sep 25, 2023• 1 citation
Publisher

Abstract

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.

DOI
10.1063/5.0168585

Related papers

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

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

Adán Garros

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.

The Journal of Chemical Physics•Oct 1, 2024•1 citation

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

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

Direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions: The influence of three-index N-representability conditions

Diego R. Alcoba, Alicia Torre, Luis Laín +7 more

This work proposes the variational determination of two-electron reduced density matrices corresponding to the ground state of N-electron systems within the doubly occupied-configuration-interaction methodology. The P, Q, and G two-index N-representability conditions have been extended to the T1 and T2 (T2') three-index ones and the resulting optimization problem has been addressed using a standard semidefinite program. We report results obtained from the doubly occupied-configuration-interaction method, from the two-index constraint variational procedure and from the two- and three-index constraint variational treatment. The discussion of these results along with a study of the computational cost demanded shows the usefulness of our proposal.

The Journal of Chemical Physics•Jan 9, 2018•37 citations

Incorporating dynamic correlation into the variational determination method of the second‐order reduced density matrix in the doubly occupied configuration interaction space

Diego R. Alcoba, Alicia Torre, Luis Laín +3 more

Abstract The variational determination of the second‐order reduced density matrices arising from N‐electron doubly occupied configuration interaction wave functions has recently been proposed as a method to describe systems possessing static correlation in their ground states. In this work, we propose to combine this approach with the on‐top pair density functional theory in order to incorporate dynamic correlations and improve the description of systems possessing both types of electron correlation. The procedure requires ensuring that the second‐order reduced density matrix elements satisfy determined N‐representability variational conditions, as well as other constraints arising from the doubly‐occupied configuration interaction wave function features. An analysis of the results obtained through this method, in studies of strongly correlated many‐electron systems, shows the usefulness of our proposals.

International Journal of Quantum Chemistry•Apr 30, 2020•7 citations
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.
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 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.
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.
This work proposes the variational determination of two-electron reduced density matrices corresponding to the ground state of N-electron systems within the doubly occupied-configuration-interaction methodology. The P, Q, and G two-index N-representability conditions have been extended to the T1 and T2 (T2') three-index ones and the resulting optimization problem has been addressed using a standard semidefinite program. We report results obtained from the doubly occupied-configuration-interaction method, from the two-index constraint variational procedure and from the two- and three-index constraint variational treatment. The discussion of these results along with a study of the computational cost demanded shows the usefulness of our proposal.
Abstract The variational determination of the second‐order reduced density matrices arising from N‐electron doubly occupied configuration interaction wave functions has recently been proposed as a method to describe systems possessing static correlation in their ground states. In this work, we propose to combine this approach with the on‐top pair density functional theory in order to incorporate dynamic correlations and improve the description of systems possessing both types of electron correlation. The procedure requires ensuring that the second‐order reduced density matrix elements satisfy determined N‐representability variational conditions, as well as other constraints arising from the doubly‐occupied configuration interaction wave function features. An analysis of the results obtained through this method, in studies of strongly correlated many‐electron systems, shows the usefulness of our proposals.