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.
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.
We establish a theoretical framework for exploring the quantum dynamics of finite ultracold bosonic ensembles based on the Born-Bogoliubov-Green-Kirkwood-Yvon (BBGKY) hierarchy of 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 symmetry conserving closure approximation is developed on the basis of a recursively formulated cluster expansion for these finite systems. 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. We apply this theoretical framework to both a tunneling and an interaction-quench scenario. Due to our efficient formulation of the theory, we can reach truncation orders as large as twelve and thereby systematically study the impact of the truncation order on the results. While the short-time dynamics is found to be excellently described with controllable accuracy, significant deviations occur on a longer time-scale in sufficiently far off-equilibrium situations. Theses deviations are accompanied by exponential-like instabilities leading to unphysical results. The phenomenology of these instabilities is investigated in detail and we show that the minimal-invasive correction algorithm of the equation of motion can indeed stabilize the BBGKY hierarchy truncated at the second order.
We establish a theoretical framework for exploring 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 symmetry conserving closure approximation is developed on the basis of a recursively formulated cluster expansion for these finite systems. In order to enforce necessary representability conditions, two minimally 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. We apply this theoretical framework to both a tunneling and an interaction-quench scenario. Due to our efficient formulation of the theory, we can reach truncation orders as large as twelve and thereby systematically study the impact of the truncation order on the results. While the short-time dynamics is found to be excellently described with controllable accuracy, significant deviations occur on a longer timescale in sufficiently far off-equilibrium situations. Theses deviations are accompanied by exponential-like instabilities leading to unphysical results. The phenomenology of these instabilities is investigated in detail and we show that the minimally invasive correction algorithm of the equation of motion can indeed stabilize the BBGKY hierarchy truncated at the second order.
The nonequilibrium description of quantum systems requires, for more than two or three particles, the use of a reduced description to be numerically tractable. Two possible approaches are based on either reduced density matrices or nonequilibrium Green functions (NEGF). Both concepts are formulated in terms of hierarchies of coupled equations—the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy for the reduced density operators and the Martin-Schwinger-hierarchy (MS) for the Green functions, respectively. In both cases, similar approximations are introduced to decouple the hierarchy, yet still many questions regarding the correspondence of both approaches remain open. Here we analyze this correspondence by studying the generalized Kadanoff–Baym ansatz (GKBA) that reduces the NEGF to a single-time theory. Starting from the BBGKY-hierarchy we present the approximations that are necessary to recover the GKBA result both, with Hartree-Fock propagators (HF-GKBA) and propagators in second Born approximation. To test the quality of the HF-GKBA, we study the dynamics of a 4-electron Hubbard nanocluster starting from a strong nonequilibrium initial state and compare to exact results and the Wang-Cassing approximation to the BBGKY hierarchy presented recently by Akbari et al. [1].
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 time evolution of -spin reduced density operators is studied for a class of Heisenberg-type quantum spin models with long-range interactions. In the framework of the quantum Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) hierarchy, we introduce an unconventional representation, different from the usual cluster expansion, which casts the hierarchy into the form of a second-order recursion. This structure suggests a scaling of the expansion coefficients and the corresponding time scales in powers of N 1/2 with the system size N , implying a separation of time scales in the large-system limit. For special parameter values and initial conditions, we can show analytically that closing the BBGKY hierarchy by neglecting -spin correlations never leads to equilibration, but gives rise to quasi-periodic time evolution with at most independent frequencies. Moreover, for the same special parameter values and in the large- N limit, we solve the complete recursion relation (the full BBGKY hierarchy), observing a superexponential decay to equilibrium in rescaled time τ = tN − 1/2 .
An accurate description of the nonequilibrium dynamics of systems with coupled spin and bosonic degrees of freedom remains theoretically challenging, especially for large system sizes and in higher than one dimension. Phase-space methods such as the truncated Wigner approximation (TWA) have the advantage of being easily scalable and applicable to arbitrary dimensions. In this work we adapt the TWA to generic spin-boson models by making use of recently developed algorithms for discrete phase spaces [J. Schachenmayer, A. Pikovski, and A. M. Rey, Phys. Rev. X 5, 011022 (2015)]. Furthermore we go beyond the standard TWA approximation by applying a scheme based on the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy of equations to our coupled spin-boson model. This allows us, in principle, to study how systematically adding higher-order corrections improves the convergence of the method. To test various levels of approximation we study an exactly solvable spin-boson model, which is particularly relevant for trapped-ion arrays. Using TWA and its BBGKY extension we accurately reproduce the time evolution of a number of one- and two-point correlation functions in several dimensions and for an arbitrary number of bosonic modes.