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

arXiv:2607.29396

Revisiting the Equation-of-Motion Method: A Universal Framework for Correlated Quantum Systems

Andrea Porro

nucl-thcond-mat.str-elphysics.chem-phquant-phSubmitted Jul 31, 2026• Announced Aug 3, 2026• 0 citations
arXivHTMLPDF

Abstract

A general implementation of the equation-of-motion (EOM) formalism for correlated many-body states is presented and applied to the description of collective excitations in atomic nuclei. While EOM approaches are traditionally formulated on top of independent-particle reference states, the present work extends the method to correlated reference states generated by modern many-body solvers. This formulation enables a consistent treatment of ground-state correlations and excited-state dynamics within a unified framework. Particular emphasis is placed on collective nuclear excitations employing chiral nuclear Hamiltonians in an ab initio context. The approach is motivated by the renewed interest in EOM techniques across several fields, including quantum chemistry and quantum computing, where they provide efficient and systematically improvable descriptions of excitation spectra. The present results demonstrate that the EOM framework offers a flexible and powerful tool for the microscopic description of nuclear spectroscopy beyond the traditional mean-field paradigm.

Comments
19 pages, 7 figures

Related papers

Ab initioexcited states from the in-medium similarity renormalization group

Nathan Parzuchowski, Titus Morris, S. K. Bogner

We present two new methods for performing ab initio calculations of excited states for closed-shell systems within the in-medium similarity renormalization group (IMSRG) framework. Both are based on combining the IMSRG with simple many-body methods commonly used to target excited states, such as the Tamm-Dancoff approximation (TDA) and equations-of-motion (EOM) techniques. In the first approach, a two-step sequential IMSRG transformation is used to drive the Hamiltonian to a form where a simple TDA calculation (i.e., diagonalization in the space of $1\mathrm{p}1\mathrm{h}$ excitations) becomes exact for a subset of eigenvalues. In the second approach, EOM techniques are applied to the IMSRG ground-state-decoupled Hamiltonian to access excited states. We perform proof-of-principle calculations for parabolic quantum dots in two dimensions and the closed-shell nuclei $^{16}\mathrm{O}$ and $^{22}\mathrm{O}$. We find that the TDA-IMSRG approach gives better accuracy than the EOM-IMSRG when calculations converge, but it is otherwise lacking the versatility and numerical stability of the latter. Our calculated spectra are in reasonable agreement with analogous EOM-coupled-cluster calculations. This work paves the way for more interesting applications of the EOM-IMSRG approach to calculations of consistently evolved observables such as electromagnetic strength functions and nuclear matrix elements, and extensions to nuclei within one or two nucleons of a closed shell by generalizing the EOM ladder operator to include particle-number nonconserving terms.

arXiv:1611.00661•Physical Review C•Apr 4, 2017•26 citations

Many-Body Perturbation Theory for Ab Initio Nuclear Structure

A. Tichai

The solution of the quantum many-body problem for medium-mass nuclei using realistic nuclear interactions poses a superbe challenge for nuclear structure research. Because an exact solution can only be provided for the lightest nuclei, one has to rely on approximate solutions when proceeding to heavier systems. Over the past years, tremendous progress has been made in the development and application of systematically improvable expansion methods and an accurate description of nuclear observables has become viable up to mass number $A \approx 100$. While closed-shell systems are consistently described via a plethora of different many-body methods, the extension to genuine open-shell systems still remains a major challenge and up to now there is no \textit{ab initio} many-body method which applies equally well to systems with even and odd mass numbers. The goal of this thesis is the development and implementation of innovative perturbative approaches with genuine open-shell capabilities. This requires the extension of well-known single-reference approaches to more general vacua. In this work we choose two complementary routes for the usage of generalized reference states. First, we derive a new \emph{ab initio} approach based on multi-configurational reference states that are conveniently derived from a prior no-core shell model calculation. Perturbative corrections are derived via second-order many-body perturbation theory, thus, merging configuration interaction and many-body perturbation theory. The generality of this ansatz enables for a treatment of medium-mass systems with arbitrary mass number, as well as the extension to low-lying excited states such that ground and excited states are treated on an equal footing. In a complementary approach, we use reference states that break a symmetry of the underlying Hamiltonian. In the simplest case this corresponds to the expansion around a particle-number-broken Hartree-Fock-Bogoliubov vacuum which is obtained from a mean-field calculation. Pairing correlations are already absorbed in the reference state. The mild scaling allows for investigations up to tin isotopic chains. All benchmarks were performed using state-of-the-art chiral two- plus three-nucleon interactions thus allowing for a universal description of nuclear observables in the medium-mass regime.

•Jan 1, 2017•1 citation

Equations-of-Motion Method and the Extended Shell Model

DAVID J. ROWE

This paper presents the equations-of-motion method as a useful and flexible tool in the study of nuclear spectroscopy. It is partly a review, but also it introduces a new and much more powerful equations-of-motion technique which supercedes the older linearization methods. The older methods worked with operator equations. To obtain closed expressions they had to be linearized in a rather arbitrary manner. The present approach works with the ground-state expectation of operator equations and thereby avoids all problems of linearization. Thus, like the Green's function method, the equations-of-motion method becomes potentially exact. It has many advantages over Green's function methods, however, among which are its greater compactness, simplicity, and the physical insight it yields.The method is first applied to rederive the random phase approximation (RPA) and the quasi-particle RPA (QRPA) and to show precisely what terms they neglect. It is demonstrated that some of these terms have coherent phases. A higher RPA and QRPA are then derived to include these terms. The corrections have some interesting effects: notably, there is a reduction of the effective interaction strength and a stabilization of the nucleus against sudden phase transitions. The equations-of-motion method is also used to generalize, in a very simple and natural way, the Hartree-Fock (HF) and Hartree-Bogolyubov (HB) concepts of independent particles and quasi-particles to nonsimple ground states.The equations-of-motion method is presented as a simple extension of the shell model to the treatment of excitationg of a correlated ground state. By concentrating on the quantities of direct physical interest, the complexity of workins with correlated wavefunctions is avoided.

Reviews of Modern Physics•Jan 1, 1968•1,127 citations

Unified theory of ground state properties of nuclei and nuclear dynamics

VICTOR A. KHODEL, Vasily R. Shaginyan

A new microscopic approach to the theory of dynamic correlation, ground state and collective excitations of many‐body systems is suggested. These dynamic and ground state properties are known as being interconnected and should be considered on the same base. A useful tool for this goal is a response function formalism. An exact functional equation for effective interaction and linear response function is obtained. The reduction method of the functional equation to a closed set of equations is given. In the simplest case the functional equation may be transformed into nonlinear integrodifferential equation. Few examples of many‐body systems are considered on the base of this equation. The developed method is applied to gas‐liquid phase transition, for evaluation of transport coefficients and collective excitation width.

AIP conference proceedings•Jan 1, 1992•0 citations

Toward an accurate strongly coupled many-body theory within the equation-of-motion framework

Елена Литвинова, Peter Schuck

Nonperturbative aspects of the quantum many-body problem are revisited, discussed and advanced in the equation of motion framework. We compare the approach to the two-fermion response function truncated on the two-body level by the cluster expansion of the dynamical interaction kernel to the approach known as time blocking approximation. Such a comparison leads to an extended many-body theory with nonperturbative treatment of high-order configurations. The present implementation of the advanced theory introduces a new class of solutions for the response functions, which include explicitly beyond-mean-field correlations between up to six fermions. The novel approach, which includes configurations with two quasiparticles coupled to two phonons ($2\mathrm{q}\ensuremath{\bigotimes}2\mathrm{phonon}$), is discussed in detail for the particle-hole nuclear response and applied to medium-mass nuclei. The proposed developments are implemented numerically on the basis of the relativistic effective meson-nucleon Lagrangian and compared to the models confined by two-fermion and four-fermion configurations, which are considered as state-of-the-art for the response theory in nuclear structure calculations. The results obtained for the dipole response of $^{42,48}\mathrm{Ca}$ and $^{68}\mathrm{Ni}$ nuclei in comparison to available experimental data show that the higher configurations are necessary for a successful description of both gross and fine details of the spectra in both high-energy and low-energy sectors.

arXiv:1908.05384•Physical Review C•Dec 24, 2019•45 citations

Equation-of-Motion Block-Correlated Coupled Cluster Method for Excited Electronic States of Strongly Correlated Systems

Haodong Zhang, Jingxiang Zou, Xiaochuan Ren +1 more

An equation-of-motion block-correlated coupled cluster method based on the generalized valence bond wave function (EOM-GVB-BCCC) is proposed to describe low-lying excited states for strongly correlated systems. The EOM-GVB-BCCC2b method with up to two-pair correlation has been implemented and tested for a few strongly correlated systems. For a water hexamer with stretched O-H bonds, which is beyond the capability of the CASSCF method, EOM-GVB-BCCC2b provides very close results as the density matrix renormalization group (DMRG). For four conjugated diradical species with triplet ground states, we found that their vertical S-T gaps from EOM-GVB-BCCC2b are also quite consistent with the DMRG results. This new method is expected to be a promising theoretical tool for describing the low-lying excited states of strongly correlated systems with large active spaces.

The Journal of Physical Chemistry Letters•Jul 21, 2023•10 citations
A general implementation of the equation-of-motion (EOM) formalism for correlated many-body states is presented and applied to the description of collective excitations in atomic nuclei. While EOM approaches are traditionally formulated on top of independent-particle reference states, the present work extends the method to correlated reference states generated by modern many-body solvers. This formulation enables a consistent treatment of ground-state correlations and excited-state dynamics within a unified framework. Particular emphasis is placed on collective nuclear excitations employing chiral nuclear Hamiltonians in an ab initio context. The approach is motivated by the renewed interest in EOM techniques across several fields, including quantum chemistry and quantum computing, where they provide efficient and systematically improvable descriptions of excitation spectra. The present results demonstrate that the EOM framework offers a flexible and powerful tool for the microscopic description of nuclear spectroscopy beyond the traditional mean-field paradigm.
We present two new methods for performing ab initio calculations of excited states for closed-shell systems within the in-medium similarity renormalization group (IMSRG) framework. Both are based on combining the IMSRG with simple many-body methods commonly used to target excited states, such as the Tamm-Dancoff approximation (TDA) and equations-of-motion (EOM) techniques. In the first approach, a two-step sequential IMSRG transformation is used to drive the Hamiltonian to a form where a simple TDA calculation (i.e., diagonalization in the space of 1p1h1\mathrm{p}1\mathrm{h}1p1h excitations) becomes exact for a subset of eigenvalues. In the second approach, EOM techniques are applied to the IMSRG ground-state-decoupled Hamiltonian to access excited states. We perform proof-of-principle calculations for parabolic quantum dots in two dimensions and the closed-shell nuclei 16O^{16}\mathrm{O}16O and 22O^{22}\mathrm{O}22O. We find that the TDA-IMSRG approach gives better accuracy than the EOM-IMSRG when calculations converge, but it is otherwise lacking the versatility and numerical stability of the latter. Our calculated spectra are in reasonable agreement with analogous EOM-coupled-cluster calculations. This work paves the way for more interesting applications of the EOM-IMSRG approach to calculations of consistently evolved observables such as electromagnetic strength functions and nuclear matrix elements, and extensions to nuclei within one or two nucleons of a closed shell by generalizing the EOM ladder operator to include particle-number nonconserving terms.
The solution of the quantum many-body problem for medium-mass nuclei using realistic nuclear interactions poses a superbe challenge for nuclear structure research. Because an exact solution can only be provided for the lightest nuclei, one has to rely on approximate solutions when proceeding to heavier systems. Over the past years, tremendous progress has been made in the development and application of systematically improvable expansion methods and an accurate description of nuclear observables has become viable up to mass number A≈100A \approx 100A≈100. While closed-shell systems are consistently described via a plethora of different many-body methods, the extension to genuine open-shell systems still remains a major challenge and up to now there is no \textit{ab initio} many-body method which applies equally well to systems with even and odd mass numbers. The goal of this thesis is the development and implementation of innovative perturbative approaches with genuine open-shell capabilities. This requires the extension of well-known single-reference approaches to more general vacua. In this work we choose two complementary routes for the usage of generalized reference states. First, we derive a new \emph{ab initio} approach based on multi-configurational reference states that are conveniently derived from a prior no-core shell model calculation. Perturbative corrections are derived via second-order many-body perturbation theory, thus, merging configuration interaction and many-body perturbation theory. The generality of this ansatz enables for a treatment of medium-mass systems with arbitrary mass number, as well as the extension to low-lying excited states such that ground and excited states are treated on an equal footing. In a complementary approach, we use reference states that break a symmetry of the underlying Hamiltonian. In the simplest case this corresponds to the expansion around a particle-number-broken Hartree-Fock-Bogoliubov vacuum which is obtained from a mean-field calculation. Pairing correlations are already absorbed in the reference state. The mild scaling allows for investigations up to tin isotopic chains. All benchmarks were performed using state-of-the-art chiral two- plus three-nucleon interactions thus allowing for a universal description of nuclear observables in the medium-mass regime.
This paper presents the equations-of-motion method as a useful and flexible tool in the study of nuclear spectroscopy. It is partly a review, but also it introduces a new and much more powerful equations-of-motion technique which supercedes the older linearization methods. The older methods worked with operator equations. To obtain closed expressions they had to be linearized in a rather arbitrary manner. The present approach works with the ground-state expectation of operator equations and thereby avoids all problems of linearization. Thus, like the Green's function method, the equations-of-motion method becomes potentially exact. It has many advantages over Green's function methods, however, among which are its greater compactness, simplicity, and the physical insight it yields.The method is first applied to rederive the random phase approximation (RPA) and the quasi-particle RPA (QRPA) and to show precisely what terms they neglect. It is demonstrated that some of these terms have coherent phases. A higher RPA and QRPA are then derived to include these terms. The corrections have some interesting effects: notably, there is a reduction of the effective interaction strength and a stabilization of the nucleus against sudden phase transitions. The equations-of-motion method is also used to generalize, in a very simple and natural way, the Hartree-Fock (HF) and Hartree-Bogolyubov (HB) concepts of independent particles and quasi-particles to nonsimple ground states.The equations-of-motion method is presented as a simple extension of the shell model to the treatment of excitationg of a correlated ground state. By concentrating on the quantities of direct physical interest, the complexity of workins with correlated wavefunctions is avoided.
A new microscopic approach to the theory of dynamic correlation, ground state and collective excitations of many‐body systems is suggested. These dynamic and ground state properties are known as being interconnected and should be considered on the same base. A useful tool for this goal is a response function formalism. An exact functional equation for effective interaction and linear response function is obtained. The reduction method of the functional equation to a closed set of equations is given. In the simplest case the functional equation may be transformed into nonlinear integrodifferential equation. Few examples of many‐body systems are considered on the base of this equation. The developed method is applied to gas‐liquid phase transition, for evaluation of transport coefficients and collective excitation width.
Nonperturbative aspects of the quantum many-body problem are revisited, discussed and advanced in the equation of motion framework. We compare the approach to the two-fermion response function truncated on the two-body level by the cluster expansion of the dynamical interaction kernel to the approach known as time blocking approximation. Such a comparison leads to an extended many-body theory with nonperturbative treatment of high-order configurations. The present implementation of the advanced theory introduces a new class of solutions for the response functions, which include explicitly beyond-mean-field correlations between up to six fermions. The novel approach, which includes configurations with two quasiparticles coupled to two phonons ($2\mathrm{q}\ensuremath{\bigotimes}2\mathrm{phonon}$), is discussed in detail for the particle-hole nuclear response and applied to medium-mass nuclei. The proposed developments are implemented numerically on the basis of the relativistic effective meson-nucleon Lagrangian and compared to the models confined by two-fermion and four-fermion configurations, which are considered as state-of-the-art for the response theory in nuclear structure calculations. The results obtained for the dipole response of 42,48Ca^{42,48}\mathrm{Ca}42,48Ca and 68Ni^{68}\mathrm{Ni}68Ni nuclei in comparison to available experimental data show that the higher configurations are necessary for a successful description of both gross and fine details of the spectra in both high-energy and low-energy sectors.
An equation-of-motion block-correlated coupled cluster method based on the generalized valence bond wave function (EOM-GVB-BCCC) is proposed to describe low-lying excited states for strongly correlated systems. The EOM-GVB-BCCC2b method with up to two-pair correlation has been implemented and tested for a few strongly correlated systems. For a water hexamer with stretched O-H bonds, which is beyond the capability of the CASSCF method, EOM-GVB-BCCC2b provides very close results as the density matrix renormalization group (DMRG). For four conjugated diradical species with triplet ground states, we found that their vertical S-T gaps from EOM-GVB-BCCC2b are also quite consistent with the DMRG results. This new method is expected to be a promising theoretical tool for describing the low-lying excited states of strongly correlated systems with large active spaces.