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

Contributions to Plasma Physics

Nonequilibrium Green Functions Approach to Strongly Correlated Fermions in Lattice Systems

Niclas Schlünzen, M. Bönitz

Published Jan 1, 2016• 55 citations
Publisher

Abstract

Abstract Quantum dynamics in strongly correlated systems are of high current interest in many fields including dense plasmas, nuclear matter and condensed matter and ultracold atoms. An important model case are fermions in lattice systems that is well suited to analyze, in detail, a variety of electronic and magnetic properties of strongly correlated solids. Such systems have recently been reproduced with fermionic atoms in optical lattices which allow for a very accurate experimental analysis of the dynamics and of transport processes such as diffusion. The theoretical analysis of such systems far from equilibrium is very challenging since quantum and spin effects as well as correlations have to be treated non‐perturbatively. The only accurate method that has been successful so far are density matrix renormalization group (DMRG) simulations. However, these simulations are presently limited to one‐dimensional (1D) systems and short times. Extension of quantum dynamics simulations to two and three dimensions is commonly viewed as one of the major challenges in this field. Recently we have reported a breakthrough in this area [N. Schlünzen et al., Phys. Rev. B (2016)] where we were able to simulate the expansion dynamics of strongly correlated fermions in a Hubbard lattice following a quench of the confinement potential in 1D, 2D and 3D. The results not only exhibited excellent agreement with the experimental data but, in addition, revealed new features of the short‐time dynamics where correlations and entanglement are being build up. The method used in this work are nonequilibrium Green functions (NEGF) which are found to be very powerful in the treatment of fermionic lattice systems filling the gap presently left open by DMRG in 2D and 3D. In this paper we present a detailed introduction in the NEGF approach and its application to inhomogeneous Hubbard clusters. In detail we discuss the proper strong coupling approximation which is given by T ‐matrix selfenergies that sum up two‐particle scattering processes to infinite order. The efficient numerical implemen‐tation of the method is discussed in detail as it has allowed us to achieve dramatic performance gains. This has been the basis for the treatment of more than 100 particles over large time intervals. The numerical results presented in this paper concentrate on the diffusion in 1D to 3D lattices. We find that the expansion dynamics consist of three different phases that are linked with the build‐up of correlations. In the long time limit, a universal scaling with the particle number is revealed. By extrapolating the expansion velocities to the macroscopic limit, the obtained results show excellent agreement with recent experiments on ultracold fermions in optical lattices. Moreover we present results for the site‐resolved behavior of correlations and entanglement that can be directly compared with experiments using the recently developed atomic microscope technique. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

DOI
10.1002/ctpp.201610003

Related papers

Ultrafast dynamics of strongly correlated fermions—nonequilibrium Green functions and selfenergy approximations

Niclas Schlünzen, S. Hermanns, Miriam Scharnke +1 more

This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). We discuss applications of NEGF simulations to describe the femtosecond dynamics of various finite fermionic systems following an excitation out of equilibrium. This includes the expansion dynamics of ultracold atoms in optical lattices following a confinement quench and the excitation of strongly correlated electrons in a solid by the impact of a charged particle. NEGF, presently, are the only ab initio quantum approach that is able to study the dynamics of correlations for long times in two and three dimensions. However, until recently, NEGF simulations have mostly been performed with rather simple selfenergy approximations such as the second-order Born approximation (SOA). While they correctly capture the qualitative trends of the relaxation towards equilibrium, the reliability and accuracy of these NEGF simulations has remained open, for a long time. Here we report on recent tests of NEGF simulations for finite lattice systems against exact-diagonalization and density-matrix-renormalization-group benchmark data. The results confirm the high accuracy and predictive capability of NEGF simulations-provided selfenergies are used that go beyond the SOA and adequately include strong correlation and dynamical-screening effects. With an extended arsenal of selfenergies that can be used effectively, the NEGF approach has the potential of becoming a powerful simulation tool with broad areas of new applications including strongly correlated solids and ultracold atoms. The present review aims at making such applications possible. To this end we present a selfcontained introduction to the theory of NEGF and give an overview on recent numerical applications to compute the ultrafast relaxation dynamics of correlated fermions. In the second part we give a detailed introduction to selfenergies beyond the SOA. Important examples are the third-order approximation, the [Formula: see text] approximation, the T-matrix approximation and the fluctuating-exchange approximation. We give a comprehensive summary of the explicit selfenergy expressions for a variety of systems of practical relevance, starting from the most general expressions (general basis) and the Feynman diagrams, and including also the important cases of diagonal basis sets, the Hubbard model and the differences occuring for bosons and fermions. With these details, and information on the computational effort and scaling with the basis size and propagation duration, readers will be able to choose the proper basis set and straightforwardly implement and apply advanced selfenergy approximations to a broad class of systems.

arXiv:1902.07038•Journal of Physics Condensed Matter•Jun 27, 2019•55 citations

Toward a Nonequilibrium Green Functions Approach to Diffusion in Strongly Coupled Finite Quantum Systems

M. Bönitz, Niclas Schlünzen, S. Hermanns

Abstract Transport properties of strongly correlated quantum systems are of central interest in condensed matter, ultracold atoms and in dense plasmas. There, the proper treatment of strong correlations poses a great challenge to theory. Here, we apply a Nonequilibrium Green Functions approach using a lattice model as a basic system. This allow us to treat a finite spatially inhomogeneous system with an arbitrary nonequilibrium initial state. Placing all particles initially to one side of the system allows for a nonequilibrium study of diffusion. Strong correlation effects are incorporated via T‐matrix selfenergies. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

arXiv:1409.2252•Contributions to Plasma Physics•Nov 27, 2014•7 citations

Nonequilibrium dynamics in the one-dimensional Fermi-Hubbard model: Comparison of the nonequilibrium Green-functions approach and the density matrix renormalization group method

Niclas Schlünzen, Jan‐Philip Joost, Fabian Heidrich‐Meisner +1 more

Experiments with ultracold fermionic atoms in optical lattices have attracted much interest in condensed matter physics as they serve as a model for electrons in solids. While experiments now allow the observation of nonequilibrium transport processes with single-site resolution, the theoretical description, especially for strongly correlated fermions, remains very challenging. Standard methods, such as the density matrix renormalization group (DMRG) and nonequilibrium Green functions (NEGF), have been primarily applied to one-dimensional systems and weak to moderate coupling, respectively. Here, the authors perform a detailed comparison of DMRG and NEGF and demonstrate that both can be benchmarked against each other. They observe complementary applicability ranges suggesting that a combination of both will allow to substantially expand the range and duration of first-principles quantum dynamics simulations for fermionic lattice systems.

arXiv:1702.01956•Physical review. B./Physical review. B•Apr 25, 2017•80 citations

Dynamics of strongly correlated fermions:Ab initioresults for two and three dimensions

Niclas Schlünzen, S. Hermanns, M. Bönitz +1 more

Quantum transport of strongly correlated fermions is of central interest in condensed matter physics. While the stationary expansion dynamics have recently been measured with cold atoms in 2D optical lattices, ab initio simulations have been limited to 1D setups so far. Here, we present the first precise fermionic quantum dynamics simulations for 2D and 3D. The simulations are based on nonequilibrium Green functions and incorporate strong correlations via $T$-matrix self-energies. The simulations predict the short-time dynamics, and we discover a universal scaling of the expansion velocity with the particle number. Our predictions can be verified experimentally using the recently developed fermionic atom microscopes.

arXiv:1508.02947•Physical review. B./Physical review. B•Jan 11, 2016•63 citations

Dynamics of Hubbard Nano‐Clusters Following Strong Excitation

M. Bönitz, S. Hermanns, Karsten Balzer

Abstract The Hubbard model is a prototype for strongly correlated electrons in condensed matter, for molecules and fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the excitation and relaxation dynamics following a perturbation of the system are only poorly explored. Here, we present results for the dynamics of electrons following nonlinear strong excitation that are based on a nonequilibrium Green functions approach. We focus on small systems—“Hubbard nano‐clusters”—that contain just a few particles where, in addition to the correlation effects, finite size effects and spatial inhomegeneity can be studied systematically. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

arXiv:1309.4574•Contributions to Plasma Physics•Dec 1, 2013•22 citations

Hubbard nanoclusters far from equilibrium

S. Hermanns, Niclas Schlünzen, M. Bönitz

The Hubbard model is a prototype for strongly correlated many-particle systems, including electrons in condensed matter and molecules, as well as for fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the nonequilibrium dynamics following a strong nonperturbative excitation only recently came into the focus of experiments and theory. It is of particular interest how the dynamics depend on the coupling strength and on the particle number and whether there exist universal features in the time evolution. Here, we present results for the dynamics of finite Hubbard clusters based on a self-consistent nonequilibrium Green functions (NEGF) approach invoking the generalized Kadanoff-Baym ansatz (GKBA). We discuss the conserving properties of the GKBA with Hartree-Fock propagators in detail and present a generalized form of the energy conservation criterion of Baym and Kadanoff for NEGF. Furthermore, we demonstrate that the HF-GKBA cures some artifacts of prior two-time NEGF simulations. Besides, this approach substantially speeds up the numerical calculations and thus presents the capability to study comparatively large systems and to extend the analysis to long times allowing for an accurate computation of the excitation spectrum via time propagation. Our data obtained within the second Born approximation compare favorably with exact diagonalization results (available for up to 13 particles) and are expected to have predictive capability for substantially larger systems in the weak-coupling limit.

arXiv:1402.7300•Physical Review B•Sep 5, 2014•80 citations
Abstract Quantum dynamics in strongly correlated systems are of high current interest in many fields including dense plasmas, nuclear matter and condensed matter and ultracold atoms. An important model case are fermions in lattice systems that is well suited to analyze, in detail, a variety of electronic and magnetic properties of strongly correlated solids. Such systems have recently been reproduced with fermionic atoms in optical lattices which allow for a very accurate experimental analysis of the dynamics and of transport processes such as diffusion. The theoretical analysis of such systems far from equilibrium is very challenging since quantum and spin effects as well as correlations have to be treated non‐perturbatively. The only accurate method that has been successful so far are density matrix renormalization group (DMRG) simulations. However, these simulations are presently limited to one‐dimensional (1D) systems and short times. Extension of quantum dynamics simulations to two and three dimensions is commonly viewed as one of the major challenges in this field. Recently we have reported a breakthrough in this area [N. Schlünzen et al., Phys. Rev. B (2016)] where we were able to simulate the expansion dynamics of strongly correlated fermions in a Hubbard lattice following a quench of the confinement potential in 1D, 2D and 3D. The results not only exhibited excellent agreement with the experimental data but, in addition, revealed new features of the short‐time dynamics where correlations and entanglement are being build up. The method used in this work are nonequilibrium Green functions (NEGF) which are found to be very powerful in the treatment of fermionic lattice systems filling the gap presently left open by DMRG in 2D and 3D. In this paper we present a detailed introduction in the NEGF approach and its application to inhomogeneous Hubbard clusters. In detail we discuss the proper strong coupling approximation which is given by T ‐matrix selfenergies that sum up two‐particle scattering processes to infinite order. The efficient numerical implemen‐tation of the method is discussed in detail as it has allowed us to achieve dramatic performance gains. This has been the basis for the treatment of more than 100 particles over large time intervals. The numerical results presented in this paper concentrate on the diffusion in 1D to 3D lattices. We find that the expansion dynamics consist of three different phases that are linked with the build‐up of correlations. In the long time limit, a universal scaling with the particle number is revealed. By extrapolating the expansion velocities to the macroscopic limit, the obtained results show excellent agreement with recent experiments on ultracold fermions in optical lattices. Moreover we present results for the site‐resolved behavior of correlations and entanglement that can be directly compared with experiments using the recently developed atomic microscope technique. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). We discuss applications of NEGF simulations to describe the femtosecond dynamics of various finite fermionic systems following an excitation out of equilibrium. This includes the expansion dynamics of ultracold atoms in optical lattices following a confinement quench and the excitation of strongly correlated electrons in a solid by the impact of a charged particle. NEGF, presently, are the only ab initio quantum approach that is able to study the dynamics of correlations for long times in two and three dimensions. However, until recently, NEGF simulations have mostly been performed with rather simple selfenergy approximations such as the second-order Born approximation (SOA). While they correctly capture the qualitative trends of the relaxation towards equilibrium, the reliability and accuracy of these NEGF simulations has remained open, for a long time. Here we report on recent tests of NEGF simulations for finite lattice systems against exact-diagonalization and density-matrix-renormalization-group benchmark data. The results confirm the high accuracy and predictive capability of NEGF simulations-provided selfenergies are used that go beyond the SOA and adequately include strong correlation and dynamical-screening effects. With an extended arsenal of selfenergies that can be used effectively, the NEGF approach has the potential of becoming a powerful simulation tool with broad areas of new applications including strongly correlated solids and ultracold atoms. The present review aims at making such applications possible. To this end we present a selfcontained introduction to the theory of NEGF and give an overview on recent numerical applications to compute the ultrafast relaxation dynamics of correlated fermions. In the second part we give a detailed introduction to selfenergies beyond the SOA. Important examples are the third-order approximation, the [Formula: see text] approximation, the T-matrix approximation and the fluctuating-exchange approximation. We give a comprehensive summary of the explicit selfenergy expressions for a variety of systems of practical relevance, starting from the most general expressions (general basis) and the Feynman diagrams, and including also the important cases of diagonal basis sets, the Hubbard model and the differences occuring for bosons and fermions. With these details, and information on the computational effort and scaling with the basis size and propagation duration, readers will be able to choose the proper basis set and straightforwardly implement and apply advanced selfenergy approximations to a broad class of systems.
Abstract Transport properties of strongly correlated quantum systems are of central interest in condensed matter, ultracold atoms and in dense plasmas. There, the proper treatment of strong correlations poses a great challenge to theory. Here, we apply a Nonequilibrium Green Functions approach using a lattice model as a basic system. This allow us to treat a finite spatially inhomogeneous system with an arbitrary nonequilibrium initial state. Placing all particles initially to one side of the system allows for a nonequilibrium study of diffusion. Strong correlation effects are incorporated via T‐matrix selfenergies. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
Experiments with ultracold fermionic atoms in optical lattices have attracted much interest in condensed matter physics as they serve as a model for electrons in solids. While experiments now allow the observation of nonequilibrium transport processes with single-site resolution, the theoretical description, especially for strongly correlated fermions, remains very challenging. Standard methods, such as the density matrix renormalization group (DMRG) and nonequilibrium Green functions (NEGF), have been primarily applied to one-dimensional systems and weak to moderate coupling, respectively. Here, the authors perform a detailed comparison of DMRG and NEGF and demonstrate that both can be benchmarked against each other. They observe complementary applicability ranges suggesting that a combination of both will allow to substantially expand the range and duration of first-principles quantum dynamics simulations for fermionic lattice systems.
Quantum transport of strongly correlated fermions is of central interest in condensed matter physics. While the stationary expansion dynamics have recently been measured with cold atoms in 2D optical lattices, ab initio simulations have been limited to 1D setups so far. Here, we present the first precise fermionic quantum dynamics simulations for 2D and 3D. The simulations are based on nonequilibrium Green functions and incorporate strong correlations via TTT-matrix self-energies. The simulations predict the short-time dynamics, and we discover a universal scaling of the expansion velocity with the particle number. Our predictions can be verified experimentally using the recently developed fermionic atom microscopes.
Abstract The Hubbard model is a prototype for strongly correlated electrons in condensed matter, for molecules and fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the excitation and relaxation dynamics following a perturbation of the system are only poorly explored. Here, we present results for the dynamics of electrons following nonlinear strong excitation that are based on a nonequilibrium Green functions approach. We focus on small systems—“Hubbard nano‐clusters”—that contain just a few particles where, in addition to the correlation effects, finite size effects and spatial inhomegeneity can be studied systematically. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
The Hubbard model is a prototype for strongly correlated many-particle systems, including electrons in condensed matter and molecules, as well as for fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the nonequilibrium dynamics following a strong nonperturbative excitation only recently came into the focus of experiments and theory. It is of particular interest how the dynamics depend on the coupling strength and on the particle number and whether there exist universal features in the time evolution. Here, we present results for the dynamics of finite Hubbard clusters based on a self-consistent nonequilibrium Green functions (NEGF) approach invoking the generalized Kadanoff-Baym ansatz (GKBA). We discuss the conserving properties of the GKBA with Hartree-Fock propagators in detail and present a generalized form of the energy conservation criterion of Baym and Kadanoff for NEGF. Furthermore, we demonstrate that the HF-GKBA cures some artifacts of prior two-time NEGF simulations. Besides, this approach substantially speeds up the numerical calculations and thus presents the capability to study comparatively large systems and to extend the analysis to long times allowing for an accurate computation of the excitation spectrum via time propagation. Our data obtained within the second Born approximation compare favorably with exact diagonalization results (available for up to 13 particles) and are expected to have predictive capability for substantially larger systems in the weak-coupling limit.