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

Physical Review C

New truncation scheme for a time-dependent density-matrix approach applied to the ground state of O16

M. Tohyama

Published Jan 8, 2015• 11 citations
Publisher

Abstract

The ground state of $^{16}\mathrm{O}$ is calculated by using a time-dependent density-matrix approach derived from a new truncation scheme of the Bogoliubov--Born--Green--Kirkwood--Yvon hierarchy for reduced density matrices, where a three-body density matrix is approximated by an antisymmetrized product of two-body density matrices. The new scheme is compared with a simpler truncation scheme previously used for the calculation of the ground state of $^{16}\mathrm{O}$ where the three-body density matrix is neglected and only two-particle--two-hole elements of the two-body density matrix are considered. It is shown that the results obtained from the two truncation schemes agree well with the exact solution.

DOI
10.1103/physrevc.91.017301

Related papers

Applications of Time-Dependent Density-Matrix Approach

M. Tohyama

The equations of motion for reduced density matrices form a coupled chain known as the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. To close the coupled chain at the two-body level, approximations for a three-body density matrix with one-body and two-body density matrices are needed. The time-dependent density-matrix theory (TDDM) assumes that the three-body density matrix is given by the antisymmetrized products of the one-body and two-body density matrices. In this review the truncation schemes of the BBGKY hierarchy beyond TDDM are discussed and a formulation for the study of excited states which is derived from the time-dependent approach is explained. The truncation schemes and the formulation for excited states are applied to the Lipkin model and the Hubbard model to corroborate their validity. Two realistic applications of the TDDM approaches are also presented. One is the dipole and quadrupole excitations of $^{40}$Ca and $^{48}$Ca and the other the fusion reactions of $^{16}$O + $^{16}$O.

Frontiers in Physics•Mar 25, 2020•18 citations

Challenges in truncating the hierarchy of time-dependent reduced density matrices equations

A. Akbari, M. J. Hashemi, Ángel Rubio +2 more

In this work, we analyze the Born, Bogoliubov, Green, Kirkwood, and Yvon (BBGKY) hierarchy of equations for describing the full time evolution of a many-body fermionic system in terms of its reduced density matrices (at all orders). We provide an exhaustive study of the challenges and open problems linked to the truncation of such a hierarchy of equations to make them practically applicable. We restrict our analysis to the coupled evolution of the one- and two-body reduced density matrices, where higher-order correlation effects are embodied into the approximation used to close the equations. We prove that within this approach, the number of electrons and total energy are conserved, regardless of the employed approximation. Further, we demonstrate that although most of the truncation schemes available in the literature give acceptable ground-state energy, when applied to describe driven electron dynamics, they exhibit undesirable and unphysical behavior, e.g., violation and even divergence in local electronic density, both in weakly and strongly correlated regimes. We illustrate and analyze these problems within the few-site Hubbard model. The model can be solved exactly and provides a unique reference for our detailed study of electron dynamics for different values of interaction, different initial conditions, and the large set of approximations considered here. Moreover, we study the role of compatibility between two hierarchical equations and positive semidefiniteness of reduced density matrices in the instability of electron dynamics. We show that even if the used approximation holds the compatibility, electron dynamics can still diverge when positive definitiveness is violated. We propose some partial solutions of such problems and point out the main paths for future work in order to make this approach applicable for the description of the correlated electron dynamics in complex systems.

arXiv:1204.4395•Physical Review B•Jun 12, 2012•48 citations

Progress in many-body theory with the equation of motion method: Time-dependent density matrix meets self-consistent RPA and applications to solvable models

Peter Schuck, M. Tohyama

The Bogoliubov-Born-Green-Kirkwood-Yvon or time-dependent density matrix (TDDM) hierarchy of equations for higher density matrices is truncated at the three-body level in approximating the three-body correlation function by a quadratic form of two-body ones, closing the equations in this way. The procedure is discussed in detail and it is shown in nontrivial model cases that the approximate inclusion of three-body correlation functions is very important to obtain precise results. A small amplitude approximation of this time-dependent nonlinear equation for the two-body correlation function is performed (STDDM*-b) and it is shown that the one-body sector of this generalized nonlinear second random phase approximation (RPA) equation is equivalent to the self-consistent RPA (SCRPA) approach which had been derived previously by different techniques. It is discussed in which way SCRPA also contains the three-body correlations. TDDM and SCRPA are tested versus exactly solvable model cases.

arXiv:1603.08276•Physical review. B./Physical review. B•Apr 14, 2016•33 citations

Extended time-dependent mean-field approximation

D. Portes, Takeshi Kodama, A. F. R. de Toledo Piza

The time-dependent mean-field approximation for two dynamically coupled subsystems is extended to include correlation effects between the subsystems, allowing for decorrelation processes to develop in the reduced density matrices. The extended scheme is formulated in terms of the truncation to M terms of the Schmidt decomposition of the full density matrix. This M natural orbitals truncation scheme is compared to the exact numerical solution for a system of two coupled anharmonic oscillators in a factorized initial state. It is found that the approximation M=3 gives a good approximation to the exact results over several characteristic times of the system. \textcopyright{} 1996 The American Physical Society.

Physical Review A•Sep 1, 1996•1 citation

The BBGKY hierarchy for ultracold bosonic systems: I. Theoretical framework

Sven Krönke, Peter Schmelcher

We establish a theoretical framework for the exploration of the quantum dynamics of finite ultracold bosonic ensembles based on the Born-Bogoliubov-Green-Kirkwood-Yvon (BBGKY) hierarchy for equations of motion for few-particle reduced density matrices (RDMs). The theory applies to zero as well as low temperatures and is formulated in a highly efficient way by utilizing dynamically optimized single-particle basis states and representing the RDMs in terms of permanents with respect to those. An energy, RDM compatibility and, if applicable, symmetry conserving closure approximation is developed on the basis of a cluster expansion for these finite systems, using suitably normalized RDMs and symmetrization operators. The underlying cluster expansion is formulated recursively by means of super-operators, which allows for going to large truncation orders. In order to enforce necessary representability conditions, two novel, minimal-invasive and energy-conserving correction algorithms are proposed, involving the dynamical purification of the solution of the truncated BBGKY hierarchy and the correction of the equations of motion themselves, respectively. For gaining conceptual insights, the impact of two-particle correlations on the dynamical quantum depletion is studied analytically. The theoretical framework provided here is extensively applied and evaluated in our immediate follow-up work [1], involving truncation orders as large as twelve.

•Dec 3, 2017•0 citations

Operator-based truncation scheme based on the many-body fermion density matrix

Siew Ann Cheong, Christopher L. Henley

In an earlier work [S. A. Cheong and C. L. Henley, preceding paper], we derived an exact formula for the many-body density matrix ${\ensuremath{\rho}}_{B}$ of a block of B sites cut out from an infinite chain of noninteracting spinless fermions, and found that the many-particle eigenvalues and eigenstates of ${\ensuremath{\rho}}_{B}$ can all be constructed out of the one-particle eigenvalues and one-particle eigenstates, respectively. In this paper we improved upon this understanding, and developed a statistical-mechanical analogy between the density-matrix eigenstates and the many-body states of a system of noninteracting fermions. Each density-matrix eigenstate corresponds to a particular set of occupation of single-particle pseudoenergy levels, and the density-matrix eigenstate with the largest weight, having the structure of a Fermi sea ground state, unambiguously defines a pseudo-Fermi level. Based on this analogy, we outlined the main ideas behind an operator-based truncation of the density-matrix eigenstates, where single-particle pseudoenergy levels far away from the pseudo-Fermi level are removed as degrees of freedom. We report numerical evidence for scaling behaviors in the single-particle pseudoenergy spectrum for different block sizes B and different filling fractions $\mathrm{n\ifmmode \bar{}\else \={}\fi{}}.$ With the aid of these scaling relations, which tell us that the block size B plays the role of an inverse temperature in the statistical-mechanical description of the density-matrix eigenstates and eigenvalues, we looked into the performance of our operator-based truncation scheme in minimizing the discarded density-matrix weight and the error in calculating the dispersion relation for elementary excitations. This performance was compared against that of the traditional density-matrix-based truncation scheme, as well as against an operator-based plane-wave truncation scheme, and found to be very satisfactory.

arXiv:cond-mat/0307172•Physical Review B•Feb 27, 2004•28 citations
The ground state of 16O^{16}\mathrm{O}16O is calculated by using a time-dependent density-matrix approach derived from a new truncation scheme of the Bogoliubov--Born--Green--Kirkwood--Yvon hierarchy for reduced density matrices, where a three-body density matrix is approximated by an antisymmetrized product of two-body density matrices. The new scheme is compared with a simpler truncation scheme previously used for the calculation of the ground state of 16O^{16}\mathrm{O}16O where the three-body density matrix is neglected and only two-particle--two-hole elements of the two-body density matrix are considered. It is shown that the results obtained from the two truncation schemes agree well with the exact solution.
The equations of motion for reduced density matrices form a coupled chain known as the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. To close the coupled chain at the two-body level, approximations for a three-body density matrix with one-body and two-body density matrices are needed. The time-dependent density-matrix theory (TDDM) assumes that the three-body density matrix is given by the antisymmetrized products of the one-body and two-body density matrices. In this review the truncation schemes of the BBGKY hierarchy beyond TDDM are discussed and a formulation for the study of excited states which is derived from the time-dependent approach is explained. The truncation schemes and the formulation for excited states are applied to the Lipkin model and the Hubbard model to corroborate their validity. Two realistic applications of the TDDM approaches are also presented. One is the dipole and quadrupole excitations of 40^{40}40Ca and 48^{48}48Ca and the other the fusion reactions of 16^{16}16O + 16^{16}16O.
In this work, we analyze the Born, Bogoliubov, Green, Kirkwood, and Yvon (BBGKY) hierarchy of equations for describing the full time evolution of a many-body fermionic system in terms of its reduced density matrices (at all orders). We provide an exhaustive study of the challenges and open problems linked to the truncation of such a hierarchy of equations to make them practically applicable. We restrict our analysis to the coupled evolution of the one- and two-body reduced density matrices, where higher-order correlation effects are embodied into the approximation used to close the equations. We prove that within this approach, the number of electrons and total energy are conserved, regardless of the employed approximation. Further, we demonstrate that although most of the truncation schemes available in the literature give acceptable ground-state energy, when applied to describe driven electron dynamics, they exhibit undesirable and unphysical behavior, e.g., violation and even divergence in local electronic density, both in weakly and strongly correlated regimes. We illustrate and analyze these problems within the few-site Hubbard model. The model can be solved exactly and provides a unique reference for our detailed study of electron dynamics for different values of interaction, different initial conditions, and the large set of approximations considered here. Moreover, we study the role of compatibility between two hierarchical equations and positive semidefiniteness of reduced density matrices in the instability of electron dynamics. We show that even if the used approximation holds the compatibility, electron dynamics can still diverge when positive definitiveness is violated. We propose some partial solutions of such problems and point out the main paths for future work in order to make this approach applicable for the description of the correlated electron dynamics in complex systems.
The Bogoliubov-Born-Green-Kirkwood-Yvon or time-dependent density matrix (TDDM) hierarchy of equations for higher density matrices is truncated at the three-body level in approximating the three-body correlation function by a quadratic form of two-body ones, closing the equations in this way. The procedure is discussed in detail and it is shown in nontrivial model cases that the approximate inclusion of three-body correlation functions is very important to obtain precise results. A small amplitude approximation of this time-dependent nonlinear equation for the two-body correlation function is performed (STDDM*-b) and it is shown that the one-body sector of this generalized nonlinear second random phase approximation (RPA) equation is equivalent to the self-consistent RPA (SCRPA) approach which had been derived previously by different techniques. It is discussed in which way SCRPA also contains the three-body correlations. TDDM and SCRPA are tested versus exactly solvable model cases.
The time-dependent mean-field approximation for two dynamically coupled subsystems is extended to include correlation effects between the subsystems, allowing for decorrelation processes to develop in the reduced density matrices. The extended scheme is formulated in terms of the truncation to M terms of the Schmidt decomposition of the full density matrix. This M natural orbitals truncation scheme is compared to the exact numerical solution for a system of two coupled anharmonic oscillators in a factorized initial state. It is found that the approximation M=3 gives a good approximation to the exact results over several characteristic times of the system. \textcopyright{} 1996 The American Physical Society.
We establish a theoretical framework for the exploration of the quantum dynamics of finite ultracold bosonic ensembles based on the Born-Bogoliubov-Green-Kirkwood-Yvon (BBGKY) hierarchy for equations of motion for few-particle reduced density matrices (RDMs). The theory applies to zero as well as low temperatures and is formulated in a highly efficient way by utilizing dynamically optimized single-particle basis states and representing the RDMs in terms of permanents with respect to those. An energy, RDM compatibility and, if applicable, symmetry conserving closure approximation is developed on the basis of a cluster expansion for these finite systems, using suitably normalized RDMs and symmetrization operators. The underlying cluster expansion is formulated recursively by means of super-operators, which allows for going to large truncation orders. In order to enforce necessary representability conditions, two novel, minimal-invasive and energy-conserving correction algorithms are proposed, involving the dynamical purification of the solution of the truncated BBGKY hierarchy and the correction of the equations of motion themselves, respectively. For gaining conceptual insights, the impact of two-particle correlations on the dynamical quantum depletion is studied analytically. The theoretical framework provided here is extensively applied and evaluated in our immediate follow-up work [1], involving truncation orders as large as twelve.
In an earlier work [S. A. Cheong and C. L. Henley, preceding paper], we derived an exact formula for the many-body density matrix ${\ensuremath{\rho}}_{B}$ of a block of B sites cut out from an infinite chain of noninteracting spinless fermions, and found that the many-particle eigenvalues and eigenstates of ${\ensuremath{\rho}}_{B}$ can all be constructed out of the one-particle eigenvalues and one-particle eigenstates, respectively. In this paper we improved upon this understanding, and developed a statistical-mechanical analogy between the density-matrix eigenstates and the many-body states of a system of noninteracting fermions. Each density-matrix eigenstate corresponds to a particular set of occupation of single-particle pseudoenergy levels, and the density-matrix eigenstate with the largest weight, having the structure of a Fermi sea ground state, unambiguously defines a pseudo-Fermi level. Based on this analogy, we outlined the main ideas behind an operator-based truncation of the density-matrix eigenstates, where single-particle pseudoenergy levels far away from the pseudo-Fermi level are removed as degrees of freedom. We report numerical evidence for scaling behaviors in the single-particle pseudoenergy spectrum for different block sizes B and different filling fractions $\mathrm{n\ifmmode \bar{}\else \={}\fi{}}.$ With the aid of these scaling relations, which tell us that the block size B plays the role of an inverse temperature in the statistical-mechanical description of the density-matrix eigenstates and eigenvalues, we looked into the performance of our operator-based truncation scheme in minimizing the discarded density-matrix weight and the error in calculating the dispersion relation for elementary excitations. This performance was compared against that of the traditional density-matrix-based truncation scheme, as well as against an operator-based plane-wave truncation scheme, and found to be very satisfactory.