Correlated real-time dynamics in large, spatially inhomogeneous quantum systems remain difficult to access with nonequilibrium many-body methods. Two-time nonequilibrium Green functions (NEGF) retain dynamical correlations but their computational runtime grows cubically with the number of time steps $N_\mathrm{t}$ . This scaling bottleneck could recently be overcome by introducing the G1--G2 scheme that is linear in $N_\mathrm{t}$ , but requires propagation of a two-particle correlation function and may suffer from numerical instabilities. This has restricted simulations to small systems with $N_\mathrm{b} \sim 10^2$ basis states. Here we introduce a quantum-fluctuation formulation of nonequilibrium Green functions, denoted $δ$ NEGF, that represents dynamical two-particle correlations through fluctuations of field-operator products, $δ\hat G$ . This guarantees stable dynamics by preserving the positivity of the reduced density matrices, avoids the explicit storage of the two-particle Green function, and reduces the propagation to a finite ensemble of Hartree-Fock-like trajectories. Combined with a stochastic low-rank decomposition of the correlation functions, the method retains time-linear scaling while extending dynamical $GW$ and particle-particle and particle-hole $T$ -matrix simulations to basis sizes of order $N_\mathrm{b}\sim 10^4$ . We benchmark $δ$ NEGF against exact and HF-GKBA results for lattice systems, finding stable correlated dynamics also at strong coupling. We further demonstrate large-scale simulations of diffusion in two-dimensional Hubbard lattices and ultrafast relaxation in graphene nanoribbon heterostructures with long-range Coulomb interactions. These results establish $δ$ NEGF as a scalable route to dynamical self-energy simulations of large, spatially inhomogeneous correlated quantum systems beyond the reach of existing NEGF implementations.Correlated real-time dynamics in large, spatially inhomogeneous quantum systems remain difficult to access with nonequilibrium many-body methods. Two-time nonequilibrium Green functions (NEGF) retain dynamical correlations but their computational runtime grows cubically with the number of time steps $N_\mathrm{t}$. This scaling bottleneck could recently be overcome by introducing the G1--G2 scheme that is linear in $N_\mathrm{t}$, but requires propagation of a two-particle correlation function and may suffer from numerical instabilities. This has restricted simulations to small systems with $N_\mathrm{b} \sim 10^2$ basis states. Here we introduce a quantum-fluctuation formulation of nonequilibrium Green functions, denoted $δ$NEGF, that represents dynamical two-particle correlations through fluctuations of field-operator products, $δ\hat G$. This guarantees stable dynamics by preserving the positivity of the reduced density matrices, avoids the explicit storage of the two-particle Green function, and reduces the propagation to a finite ensemble of Hartree-Fock-like trajectories. Combined with a stochastic low-rank decomposition of the correlation functions, the method retains time-linear scaling while extending dynamical $GW$ and particle-particle and particle-hole $T$-matrix simulations to basis sizes of order $N_\mathrm{b}\sim 10^4$. We benchmark $δ$NEGF against exact and HF-GKBA results for lattice systems, finding stable correlated dynamics also at strong coupling. We further demonstrate large-scale simulations of diffusion in two-dimensional Hubbard lattices and ultrafast relaxation in graphene nanoribbon heterostructures with long-range Coulomb interactions. These results establish $δ$NEGF as a scalable route to dynamical self-energy simulations of large, spatially inhomogeneous correlated quantum systems beyond the reach of existing NEGF implementations.

arXiv:2606.10773Jun 9, 20260 citations

Time‐Resolved Optical Response of the Dicke's Model via the Nonequilibrium Green's Function Approach

Megha Gopalakrishna, Y. Pavlyukh, Claudio Verdozzi

Due to their conceptual appeal and computational convenience, two‐level systems (TLSs) and generalizations are used to investigate nonlinear behavior in quantum optics and assess the applicability of theoretical methods. Herein, the focus is on second‐harmonic generation (SHG) and, as system of interest, on the Dicke model, which consists of several TLSs inside an optical cavity. The main aspect addressed is the scope of nonequilibrium Green's function (NEGF) to describe the effect of inhomogeneities and electron–electron (e–e) interactions on the SHG signal. For benchmarking purposes, exact diagonalization (ED) results are also presented and discussed. SHG spectra obtained with NEGF and ED are found to be in very good mutual agreement in most situations. Furthermore, inhomogeneity in the TLS and e–e interactions reduce the SHG signal, and the reduction is stronger with inhomogeneity than with interactions. This trend is consistently noted across different (small to large) system sizes. Finally, a modified NEGF approach is proposed to account for cavity leakage, where the quantum photon fields are coupled to a bath of classical oscillators. As expected, within this mixed quantum‐classical scheme, a decrease in the intensity of the fluorescent spectra takes place depending on the entity of leakage.Due to their conceptual appeal and computational convenience, two‐level systems (TLSs) and generalizations are used to investigate nonlinear behavior in quantum optics and assess the applicability of theoretical methods. Herein, the focus is on second‐harmonic generation (SHG) and, as system of interest, on the Dicke model, which consists of several TLSs inside an optical cavity. The main aspect addressed is the scope of nonequilibrium Green's function (NEGF) to describe the effect of inhomogeneities and electron–electron (e–e) interactions on the SHG signal. For benchmarking purposes, exact diagonalization (ED) results are also presented and discussed. SHG spectra obtained with NEGF and ED are found to be in very good mutual agreement in most situations. Furthermore, inhomogeneity in the TLS and e–e interactions reduce the SHG signal, and the reduction is stronger with inhomogeneity than with interactions. This trend is consistently noted across different (small to large) system sizes. Finally, a modified NEGF approach is proposed to account for cavity leakage, where the quantum photon fields are coupled to a bath of classical oscillators. As expected, within this mixed quantum‐classical scheme, a decrease in the intensity of the fluorescent spectra takes place depending on the entity of leakage.

arXiv:2312.13874physica status solidi (b)May 23, 20242 citations

Nonequilibrium Green’s functions formalism

WORLD SCIENTIFIC eBooksMay 17, 20170 citations

Nonequilibrium Green's functions formalism

WORLD SCIENTIFIC eBooksJun 1, 20100 citations

N_\mathrm{t}

. This scaling bottleneck could recently be overcome by introducing the G1--G2 scheme that is linear in

N_\mathrm{t}

, but requires propagation of a two-particle correlation function and may suffer from numerical instabilities. This has restricted simulations to small systems with

N_\mathrm{b} \sim 10^2

basis states. Here we introduce a quantum-fluctuation formulation of nonequilibrium Green functions, denoted

δ

NEGF, that represents dynamical two-particle correlations through fluctuations of field-operator products,

δ\hat G

. This guarantees stable dynamics by preserving the positivity of the reduced density matrices, avoids the explicit storage of the two-particle Green function, and reduces the propagation to a finite ensemble of Hartree-Fock-like trajectories. Combined with a stochastic low-rank decomposition of the correlation functions, the method retains time-linear scaling while extending dynamical

GW

and particle-particle and particle-hole

T

-matrix simulations to basis sizes of order

N_\mathrm{b}\sim 10^4

. We benchmark

δ

NEGF against exact and HF-GKBA results for lattice systems, finding stable correlated dynamics also at strong coupling. We further demonstrate large-scale simulations of diffusion in two-dimensional Hubbard lattices and ultrafast relaxation in graphene nanoribbon heterostructures with long-range Coulomb interactions. These results establish

δ

NEGF as a scalable route to dynamical self-energy simulations of large, spatially inhomogeneous correlated quantum systems beyond the reach of existing NEGF implementations.