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:1912.11004 · Pure and Applied Analysis

Beyond Bogoliubov dynamics

Lea Boßmann, Sören Petrat, Peter Pickl, Avy Soffer

Published Dec 31, 2021• 27 citations
arXivHTMLPDF

Abstract

We consider a system of N interacting bosons in the mean-field scaling regime and construct corrections to the Bogoliubov dynamics that approximate the true N-body dynamics in norm to arbitrary precision. The N-independent corrections are given in terms of the solutions of the Bogoliubov and Hartree equations and satisfy a generalized form of Wick's theorem. We determine the n-point correlation functions of the excitations around the condensate, as well as the reduced densities of the N-body system, to arbitrary accuracy, given only the knowledge of the two-point functions of a quasi-free state and the solution of the Hartree equation. In this way, the complex problem of computing all n-point correlation functions for an interacting N-body system is essentially reduced to the problem of solving the Hartree equation and the PDEs for the Bogoliubov two-point functions.

DOI
10.2140/paa.2021.3.677

Related papers

Higher Order Corrections to the Mean-Field Description of the Dynamics of Interacting Bosons

Lea Boßmann, Nataša Pavlović, Peter Pickl +1 more

Abstract In this paper, we introduce a novel method for deriving higher order corrections to the mean-field description of the dynamics of interacting bosons. More precisely, we consider the dynamics of N $$d$$ d -dimensional bosons for large N. The bosons initially form a Bose–Einstein condensate and interact with each other via a pair potential of the form $$(N-1)^{-1}N^{d\beta }v(N^\beta \cdot )$$ (N-1)-1Ndβv(Nβ·) for $$\beta \in [0,\frac{1}{4d})$$ β∈[0,14d) . We derive a sequence of N-body functions which approximate the true many-body dynamics in $$L^2({\mathbb {R}}^{dN})$$ L2(RdN) -norm to arbitrary precision in powers of $$N^{-1}$$ N-1 . The approximating functions are constructed as Duhamel expansions of finite order in terms of the first quantised analogue of a Bogoliubov time evolution.

arXiv:1905.06164•Journal of Statistical Physics•Feb 21, 2020•24 citations

Numerical Bogoliubov approximation of bosonic many-body quantum dynamics

Yoann Le Hénaff, Christian Lubich, Peter Pickl

We propose an algorithm for the numerical solution of the time-dependent Schrödinger equation for N >> 1bosons interacting via a two-body interaction potential in the mean-field scaling regime. Unlike existing numerical algorithms for quantum dynamics, the proposed algorithm has a computational complexity that is independent of the number N of particles. It is based on an ansatz for the wave function inspired by the Bogoliubov theory, which is an asymptotic theory for large N. In a first part, we follow Bogoliubov theory in evolving the condensate function by the Hartree equation, but we use the Dirac--Frenkel time-dependent variational principle to derive nonlinear equations of motion for the excitation functions in the Bogoliubov ansatz. By carefully exploiting the bosonic symmetry, the so obtained approximate wave functions are shown to have derivatives bounded independently of N apart from a time-dependent global phase factor which rotates with a frequency proportional to N. In the second part, we numerically approximate the wave function (modulo the phase) by a composite algorithm, which employs an appropriate splitting of the equations of motion and uses a dynamical low-rank approximation of the fundamental two-particle excitation function by linear combinations of products of single-particle functions. The low-rank excitation function is evolved in time by a basis update and Galerkin (BUG) integrator that is robust to small singular values of the coefficient matrix and allows for rank adaptivity. Numerical results are provided to assess the performance of the proposed algorithm.

arXiv:2610.08956•Oct 6, 2026•0 citations

Derivation of mean-field equations and next-order corrections for bosons and fermions

David Mitrouskas

This thesis is about the derivation of effective mean field equations and their next-order corrections starting from nonrelativistic many-body quantum theory. Mean field equations provide an approximate ansatz for the description of interacting many-particle systems. In this ansatz the interaction between the particles is replaced by a self-consistent external potential leading to an effective one-body description of the many-particle system. Next-order corrections provide an approximation which goes one step further and tries to capture also subleading effects that are not resolved by the mean field ansatz. We present mathematical proofs for the validity of such effective theories for different models that are motivated, e.g., from the theory of ultracold atoms (the bosonic Hartree equation and the corresponding Bogoliubov theory) and from plasma physics (the motion of a tracer particle through a degenerate and dense electron gas). Starting from a many-body Schrödinger equation, our goal is to show that the solutions converge in a particular limit to the solutions to an effective mean field equation and its next-order corrections. After a short introduction and a summary in Chapter one, we present the main part of this work in three self-contained chapters. In Chapter two we analyze the dynamics of a large number N of nonrelativistic bosons in the weak coupling limit, i.e., for a coupling constant g_N=1/N. It is well known that in the limit of infinite particle number, the Hartree equation emerges as an effective one-particle theory of the Bose gas. This is closely related to the remarkable physical phenomenon of Bose-Einstein condensation at low temperature, namely that the majority of particles in a Bose gas occupies the same copy of a single one-particle quantum state. Our emphasis lies in the description of the few particles that fluctuate around the Bose-Einstein condensate. We show convergence of the fully interacting dynamics to an auxiliary time evolution in the norm of the N-particle space. This result allows us to prove several other assertions. Among other things, it is used to derive the Hartree equation with optimal speed of convergence 1/N for initial states that are close to ground states of interacting systems and also to prove convergence of the N-particle solution towards a time evolution obtained from the Bogoliubov Hamiltonian on Fock space. Chapter three is about the low energy properties of the weakly interacting homogeneous Bose gas. Here, we derive a novel estimate for low energy eigenfunctions stating that the probability for finding $l$ particles out of their total number N not in the condensate is exponentially small in the number $l$. Using this bound, we then prove that the ground state wave function of the microscopic model satisfies certain quasifree type properties. The exponential decay is moreover used to provide an alternative proof for the validity of Bogoliubov's approximation for the low-lying energy eigenvalues. Bogoliubov theory states that the excitation energies of the Bose gas are given by excitations of free quasiparticles obeying an effective energy-momentum dispersion relation which is linear for small momentum. The linearity of the effective dispersion relation is an essential ingredient for the explanation of the superfluid character of the Bose gas. In Chapter four we study the time evolution of a single particle coupled through a pair potential to a dense and homogeneous ideal Fermi gas in two spatial dimensions. This type of model is well known in plasma physics where it is used to describe the energy loss of ions moving through a dense and degenerate electron gas. We analyze the model for a coupling parameter g=1 and prove closeness of the time evolution to an effective dynamics for large densities of the gas and for long time scales of the order of some power of the density. The effective dynamics is generated by the free Hamiltonian with a large but constant energy shift. To leading order, this energy shift is given by the spatially homogeneous mean field potential produced by the gas particles, whereas at next-to-leading order, one has to consider an additional correction to the mean field energy which is due to so-called recollision processes.

•Jan 1, 2017•3 citations

Quantum fluctuations of many-body dynamics around the Gross–Pitaevskii equation

Cristina Caraci, Jakob Oldenburg, Benjamin Schlein

We consider the evolution of a gas of N bosons in the three-dimensional Gross–Pitaevskii regime (in which particles are initially trapped in a volume of order one and interact through a repulsive potential with scattering length of order 1/N ). We construct a quasi-free approximation of the many-body dynamics, whose distance to the solution of the Schrödinger equation converges to zero, as N \to \infty , in the L^{2} (\mathbb{R}^{3N}) -norm. To achieve this goal, we let the Bose–Einstein condensate evolve according to a time-dependent Gross–Pitaevskii equation. After factoring out the microscopic correlation structure, the evolution of the orthogonal excitations of the condensate is governed instead by a Bogoliubov dynamics, with a time-dependent generator quadratic in creation and annihilation operators. As an application, we show a central limit theorem for fluctuations of bounded observables around their expectation with respect to the Gross–Pitaevskii dynamics.

arXiv:2308.11687•Annales de l Institut Henri Poincaré C Analyse Non Linéaire•Jun 27, 2024•7 citations

Bogoliubov corrections and trace norm convergence for the Hartree dynamics

David Mitrouskas, Sören Petrat, Peter Pickl

We consider the dynamics of a large number [Formula: see text] of nonrelativistic bosons in the mean field limit for a class of interaction potentials that includes Coulomb interaction. In order to describe the fluctuations around the mean field Hartree state, we introduce an auxiliary Hamiltonian on the [Formula: see text]-particle space that is similar to the one obtained from Bogoliubov theory. We show convergence of the auxiliary time evolution to the fully interacting dynamics in the norm of the [Formula: see text]-particle space. This result allows us to prove several other results: convergence of reduced density matrices in trace norm with optimal rate, convergence in energy trace norm, and convergence to a time evolution obtained from the Bogoliubov Hamiltonian on Fock space with expected optimal rate. We thus extend and quantify several previous results, e.g., by providing the physically important convergence rates, including time-dependent external fields and singular interactions, and allowing for more general initial states, e.g., those that are expected to be ground states of interacting systems.

arXiv:1609.06264•Reviews in Mathematical Physics•Feb 8, 2019•41 citations

A Central Limit Theorem in Many-Body Quantum Dynamics

Gérard Ben Arous, Kay Kirkpatrick, Benjamin Schlein

We study the many body quantum evolution of bosonic systems in the mean field limit. The dynamics is known to be well approximated by the Hartree equation. So far, the available results have the form of a law of large numbers. In this paper we go one step further and we show that the fluctuations around the Hartree evolution satisfy a central limit theorem. Interestingly, the variance of the limiting Gaussian distribution is determined by a time-dependent Bogoliubov transformation describing the dynamics of initial coherent states in a Fock space representation of the system.

arXiv:1111.6999•Nov 29, 2011•3 citations
We consider a system of N interacting bosons in the mean-field scaling regime and construct corrections to the Bogoliubov dynamics that approximate the true N-body dynamics in norm to arbitrary precision. The N-independent corrections are given in terms of the solutions of the Bogoliubov and Hartree equations and satisfy a generalized form of Wick's theorem. We determine the n-point correlation functions of the excitations around the condensate, as well as the reduced densities of the N-body system, to arbitrary accuracy, given only the knowledge of the two-point functions of a quasi-free state and the solution of the Hartree equation. In this way, the complex problem of computing all n-point correlation functions for an interacting N-body system is essentially reduced to the problem of solving the Hartree equation and the PDEs for the Bogoliubov two-point functions.
Abstract In this paper, we introduce a novel method for deriving higher order corrections to the mean-field description of the dynamics of interacting bosons. More precisely, we consider the dynamics of N ddd d -dimensional bosons for large N. The bosons initially form a Bose–Einstein condensate and interact with each other via a pair potential of the form (N−1)−1Ndβv(Nβ⋅)(N-1)^{-1}N^{d\beta }v(N^\beta \cdot )(N−1)−1Ndβv(Nβ⋅) (N-1)-1Ndβv(Nβ·) for β∈[0,14d)\beta \in [0,\frac{1}{4d})β∈[0,4d1​) β∈[0,14d) . We derive a sequence of N-body functions which approximate the true many-body dynamics in L2(RdN)L^2({\mathbb {R}}^{dN})L2(RdN) L2(RdN) -norm to arbitrary precision in powers of N−1N^{-1}N−1 N-1 . The approximating functions are constructed as Duhamel expansions of finite order in terms of the first quantised analogue of a Bogoliubov time evolution.
We propose an algorithm for the numerical solution of the time-dependent Schrödinger equation for N >> 1bosons interacting via a two-body interaction potential in the mean-field scaling regime. Unlike existing numerical algorithms for quantum dynamics, the proposed algorithm has a computational complexity that is independent of the number N of particles. It is based on an ansatz for the wave function inspired by the Bogoliubov theory, which is an asymptotic theory for large N. In a first part, we follow Bogoliubov theory in evolving the condensate function by the Hartree equation, but we use the Dirac--Frenkel time-dependent variational principle to derive nonlinear equations of motion for the excitation functions in the Bogoliubov ansatz. By carefully exploiting the bosonic symmetry, the so obtained approximate wave functions are shown to have derivatives bounded independently of N apart from a time-dependent global phase factor which rotates with a frequency proportional to N. In the second part, we numerically approximate the wave function (modulo the phase) by a composite algorithm, which employs an appropriate splitting of the equations of motion and uses a dynamical low-rank approximation of the fundamental two-particle excitation function by linear combinations of products of single-particle functions. The low-rank excitation function is evolved in time by a basis update and Galerkin (BUG) integrator that is robust to small singular values of the coefficient matrix and allows for rank adaptivity. Numerical results are provided to assess the performance of the proposed algorithm.
This thesis is about the derivation of effective mean field equations and their next-order corrections starting from nonrelativistic many-body quantum theory. Mean field equations provide an approximate ansatz for the description of interacting many-particle systems. In this ansatz the interaction between the particles is replaced by a self-consistent external potential leading to an effective one-body description of the many-particle system. Next-order corrections provide an approximation which goes one step further and tries to capture also subleading effects that are not resolved by the mean field ansatz. We present mathematical proofs for the validity of such effective theories for different models that are motivated, e.g., from the theory of ultracold atoms (the bosonic Hartree equation and the corresponding Bogoliubov theory) and from plasma physics (the motion of a tracer particle through a degenerate and dense electron gas). Starting from a many-body Schrödinger equation, our goal is to show that the solutions converge in a particular limit to the solutions to an effective mean field equation and its next-order corrections. After a short introduction and a summary in Chapter one, we present the main part of this work in three self-contained chapters. In Chapter two we analyze the dynamics of a large number N of nonrelativistic bosons in the weak coupling limit, i.e., for a coupling constant g_N=1/N. It is well known that in the limit of infinite particle number, the Hartree equation emerges as an effective one-particle theory of the Bose gas. This is closely related to the remarkable physical phenomenon of Bose-Einstein condensation at low temperature, namely that the majority of particles in a Bose gas occupies the same copy of a single one-particle quantum state. Our emphasis lies in the description of the few particles that fluctuate around the Bose-Einstein condensate. We show convergence of the fully interacting dynamics to an auxiliary time evolution in the norm of the N-particle space. This result allows us to prove several other assertions. Among other things, it is used to derive the Hartree equation with optimal speed of convergence 1/N for initial states that are close to ground states of interacting systems and also to prove convergence of the N-particle solution towards a time evolution obtained from the Bogoliubov Hamiltonian on Fock space. Chapter three is about the low energy properties of the weakly interacting homogeneous Bose gas. Here, we derive a novel estimate for low energy eigenfunctions stating that the probability for finding lll particles out of their total number N not in the condensate is exponentially small in the number lll. Using this bound, we then prove that the ground state wave function of the microscopic model satisfies certain quasifree type properties. The exponential decay is moreover used to provide an alternative proof for the validity of Bogoliubov's approximation for the low-lying energy eigenvalues. Bogoliubov theory states that the excitation energies of the Bose gas are given by excitations of free quasiparticles obeying an effective energy-momentum dispersion relation which is linear for small momentum. The linearity of the effective dispersion relation is an essential ingredient for the explanation of the superfluid character of the Bose gas. In Chapter four we study the time evolution of a single particle coupled through a pair potential to a dense and homogeneous ideal Fermi gas in two spatial dimensions. This type of model is well known in plasma physics where it is used to describe the energy loss of ions moving through a dense and degenerate electron gas. We analyze the model for a coupling parameter g=1 and prove closeness of the time evolution to an effective dynamics for large densities of the gas and for long time scales of the order of some power of the density. The effective dynamics is generated by the free Hamiltonian with a large but constant energy shift. To leading order, this energy shift is given by the spatially homogeneous mean field potential produced by the gas particles, whereas at next-to-leading order, one has to consider an additional correction to the mean field energy which is due to so-called recollision processes.
We consider the evolution of a gas of N bosons in the three-dimensional Gross–Pitaevskii regime (in which particles are initially trapped in a volume of order one and interact through a repulsive potential with scattering length of order 1/N ). We construct a quasi-free approximation of the many-body dynamics, whose distance to the solution of the Schrödinger equation converges to zero, as N \to \infty , in the L^{2} (\mathbb{R}^{3N}) -norm. To achieve this goal, we let the Bose–Einstein condensate evolve according to a time-dependent Gross–Pitaevskii equation. After factoring out the microscopic correlation structure, the evolution of the orthogonal excitations of the condensate is governed instead by a Bogoliubov dynamics, with a time-dependent generator quadratic in creation and annihilation operators. As an application, we show a central limit theorem for fluctuations of bounded observables around their expectation with respect to the Gross–Pitaevskii dynamics.
We consider the dynamics of a large number [Formula: see text] of nonrelativistic bosons in the mean field limit for a class of interaction potentials that includes Coulomb interaction. In order to describe the fluctuations around the mean field Hartree state, we introduce an auxiliary Hamiltonian on the [Formula: see text]-particle space that is similar to the one obtained from Bogoliubov theory. We show convergence of the auxiliary time evolution to the fully interacting dynamics in the norm of the [Formula: see text]-particle space. This result allows us to prove several other results: convergence of reduced density matrices in trace norm with optimal rate, convergence in energy trace norm, and convergence to a time evolution obtained from the Bogoliubov Hamiltonian on Fock space with expected optimal rate. We thus extend and quantify several previous results, e.g., by providing the physically important convergence rates, including time-dependent external fields and singular interactions, and allowing for more general initial states, e.g., those that are expected to be ground states of interacting systems.
We study the many body quantum evolution of bosonic systems in the mean field limit. The dynamics is known to be well approximated by the Hartree equation. So far, the available results have the form of a law of large numbers. In this paper we go one step further and we show that the fluctuations around the Hartree evolution satisfy a central limit theorem. Interestingly, the variance of the limiting Gaussian distribution is determined by a time-dependent Bogoliubov transformation describing the dynamics of initial coherent states in a Fock space representation of the system.