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

Quarterly of Applied Mathematics

Evolution of the Boson gas at zero temperature: Mean-field limit and second-order correction

Manoussos G. Grillakis, Matei Machedon, Dionisios Margetis

Published Aug 24, 2016• 5 citations
Publisher

Abstract

A large system of N N integer-spin atoms, called Bosons, manifests one of the most coherent macroscopic quantum states known to date, the “Bose-Einstein condensate”, at extremely low temperatures. As N → ∞ N\to \infty , this system is usually described by a mean-field limit: a single-particle wave function, the condensate wave function , that satisfies a nonlinear Schrödinger-type equation. In this expository paper, we review kinetic aspects of the mean-field Boson evolution. Furthermore, we discuss recent advances in the rigorous study of second-order corrections to this mean-field limit. These corrections originate from the quantum-kinetic mechanism of pair excitation , which lies at the core of pioneering works in theoretical physics including ideas of Bogoliubov, Lee, Huang, Yang and Wu. In the course of our exposition, we revisit the formalism of Fock space, which is indispensable for the analysis of pair excitation.

DOI
10.1090/qam/1455

Related papers

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

Bose–Einstein Condensation beyond Mean Field: Many-Body Bound State of Periodic Microstructure

Dionisios Margetis

We study stationary quantum fluctuations around a mean field limit in trapped, dilute atomic gases of repulsively interacting bosons at zero temperature. Our goal is to describe quantum-mechanically the lowest macroscopic many-body bound state consistent with a microscopic Hamiltonian that accounts for inhomogeneous particle scattering processes. In the mean field limit, the wave function of the condensate (macroscopic quantum state) satisfies a defocusing cubic nonlinear Schrödinger-type equation, the Gross–Pitaevskii equation. We include consequences of pair excitation, i.e., the scattering of particles in pairs from the condensate to other states, proposed in [T. T. Wu, J. Math. Phys., 2 (1961), pp. 105–123]. Our derivations rely on an uncontrolled yet physically motivated assumption for the many-body wave function. By relaxing mathematical rigor, from a particle Hamiltonian with a spatially varying interaction strength we derive via heuristics an integro–partial differential equation for the pair collision kernel, K, under a stationary condensate wave function, $\Phi$. For a scattering length with periodic microstructure of subscale $\epsilon$, we formally describe via classical homogenization the lowest many-body bound state in terms of $\Phi$ and K up to second order in $\epsilon$. If the external potential is slowly varying, we solve the homogenized equations via boundary layer theory. As an application, we describe the partial depletion of the condensate.

Multiscale Modeling and Simulation•Jan 1, 2012•5 citations

Solvable model for pair excitation in trapped Boson gas at zero temperature

Dionisios Margetis

In Bose-Einstein condensation, a macroscopically large number of particles occupy the same single-particle quantum state. Our goal is to study time-dependent aspects of particle excitations to states other than the single-particle macroscopic state in trapped dilute atomic gases. We adopt the view that atoms are excited in pairs so that their scattering from the single-particle state to vector positions x and y at time t is described by the pair-excitation function, K0(x,y,t). We solve a nonlocal equation for K0 under a slowly varying external potential by assuming that the wave function of the macroscopic state satisfies a time-independent nonlinear Schroedinger equation. For zero initial excitation (K0=0 at t=0) and sufficiently large t, we evaluate asymptotically K0 in terms of the one-variable Lommel function for any distance |x-y|.

arXiv:0709.2645•Journal of Physics A Mathematical and Theoretical•Aug 22, 2008•4 citations

Bose-Einstein condensation at finite temperatures: Mean field laws with periodic microstructure

Dionisios Margetis

At finite temperatures below the phase transition point, the Bose-Einstein condensation, the macroscopic occupation of a single quantum state by particles of integer spin, is not complete. In the language of superfluid helium, this means that the superfluid coexists with the normal fluid. Our goal is to describe this coexistence in trapped, dilute atomic gases with repulsive interactions via mean field laws that account for a {\em spatially varying} particle interaction strength. By starting with the $N$-body Hamiltonian, $N\gg 1$, we formally derive a system of coupled, nonlinear evolution equations in $3+1$ dimensions for the following quantities: (i) the wave function of the macroscopically occupied state; and (ii) the single-particle wave functions of thermally excited states. For stationary (bound) states and a scattering length with {\em periodic microstructure} of subscale $ε$, we heuristically extract effective equations of motion via periodic homogenization up to second order in $ε$.

arXiv:1509.04364•Sep 15, 2015•1 citation

Quantitative Derivation of Effective Evolution Equations for the Dynamics of Bose-Einstein Condensates

Elif Kuz

This thesis proves certain results concerning an important question in non-equilibrium quantum statistical mechanics which is the derivation of effective evolution equations approximating the dynamics of a system of large number of bosons initially at equilibrium (ground state at very low temperatures). The dynamics of such systems are governed by the time-dependent linear many-body Schroedinger equation from which it is typically difficult to extract useful information due to the number of particles being large. We will study quantitatively (i.e. with explicit bounds on the error) how a suitable one particle non-linear Schroedinger equation arises in the mean field limit as number of particles N → ∞ and how the appropriate corrections to the mean field will provide better approximations of the exact dynamics. In the first part of this thesis we consider the evolution of N bosons, where N is large, with two-body interactions of the form N³ᵝv(Nᵝ⋅), 0≤β≤1. The parameter β measures the strength and the range of interactions. We compare the exact evolution with an approximation which considers the evolution of a mean field coupled with an appropriate description of pair excitations, see [18,19] by Grillakis-Machedon-Margetis. We extend the results for 0 ≤ β < 1/3 in [19, 20] to the case of β < 1/2 and obtain an error bound of the form p(t)/Nᵅ, where α>0 and p(t) is a polynomial, which implies a specific rate of convergence as N → ∞. In the second part, utilizing estimates of the type discussed in the first part, we compare the exact evolution with the mean field approximation in the sense of marginals. We prove that the exact evolution is close to the approximate in trace norm for times of the order o(1)√N compared to log(o(1)N) as obtained in Chen-Lee-Schlein [6] for the Hartree evolution. Estimates of similar type are obtained for stronger interactions as well.

•Jan 1, 2016•0 citations

Equilibrium and Non-Equilibrium Properties of Dilute Bose Gases

Oldenburg, Jakob Felix

Quantum systems consisting of many particles are both interesting from the physics and the mathematics perspective. One aspect that makes them fascinating from the physics point of view is that they can exhibit quantum effects on macroscopic scales such as Bose-Einstein condensation. This phenomenon, first predicted 100 years ago by Bose [19] and Einstein [44], occurs when a dilute gas of bosons is cooled down to very low temperatures. Below a certain critical temperature a condensate forms, meaning that a macroscopic fraction of the particles is in the same one-particle state, and thus, quantum phenomena become relevant on length scales much larger than the individual atom. The occurrence of Bose-Einstein condensation is also related to superfluidity as suggested in [88]. Bose-Einstein condensation has been achieved experimentally independently by the groups of Cornell and Wieman [3] and Ketterle [40] earning them a Nobel prize in 2001. While Bose and Einstein analyzed theoretically the non-interacting gas, an effective theory describing weakly interacting Bose gases was proposed by Bogoliubov in 1947 with the goal of explaining the occurrence of superfluidity in liquid helium [18]. Using assumptions such as that condensation indeed occurs, his theory yields predictions for the low-energy states of the system. Seen from the mathematics perspective, describing the Bose gas is an interesting challenge due to the strong correlations between particles introduced by the interaction making this case much harder to study than the non-interacting one. In particular, a proof that Bose-Einstein condensation occurs in the thermodynamic limit for continuous systems is still out of reach. However, there has been a lot of work regarding the rigorous analysis of bosonic systems in various relevant regimes. This thesis contains three new contributions establishing rigorous results on bosonic systems in different, physically relevant settings. It is based on the articles [37, 27, 31]. The first result concerns the time-evolution of initially trapped Bose gases in the ultra-dilute Gross-Pitaevskii regime. Under some assumptions on the initial data satisfied in particular by a class of states that are expected to approximate the ground state in a trap, we derive a quasi-free approximation that is close to the true solution in $L^2$-norm. In particular, we show that this approximate wave function is described in terms of a Bose-Einstein condensate in a one-particle state following the time-dependent Gross-Pitaevskii equation together with pair excitations with respect to this condensate. As a corollary, we establish a central limit theorem for bounded one-particle observables. The second result is a simple proof of Bose-Einstein condensation for interacting, ultra-dilute gases. More specifically, we consider scalings interpolating between the Gross-Pitaevskii regime and the thermodynamic limit. While condensation had been shown to hold in these regimes before, the advantage of the method used here is its simplicity. Lastly, the third result establishes an upper bound on the ground state energy of a dilute Bose gas in the thermodynamic limit. It captures rigorously the third order correction to the energy per unit volume for sufficiently low densities as predicted in the physics literature first by Wu and subsequently others. $\textit{Full-text embargoed until: 2026-01-01}$

•Dec 9, 2025•0 citations
A large system of N N integer-spin atoms, called Bosons, manifests one of the most coherent macroscopic quantum states known to date, the “Bose-Einstein condensate”, at extremely low temperatures. As N → ∞ N\to \infty , this system is usually described by a mean-field limit: a single-particle wave function, the condensate wave function , that satisfies a nonlinear Schrödinger-type equation. In this expository paper, we review kinetic aspects of the mean-field Boson evolution. Furthermore, we discuss recent advances in the rigorous study of second-order corrections to this mean-field limit. These corrections originate from the quantum-kinetic mechanism of pair excitation , which lies at the core of pioneering works in theoretical physics including ideas of Bogoliubov, Lee, Huang, Yang and Wu. In the course of our exposition, we revisit the formalism of Fock space, which is indispensable for the analysis of pair excitation.
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 study stationary quantum fluctuations around a mean field limit in trapped, dilute atomic gases of repulsively interacting bosons at zero temperature. Our goal is to describe quantum-mechanically the lowest macroscopic many-body bound state consistent with a microscopic Hamiltonian that accounts for inhomogeneous particle scattering processes. In the mean field limit, the wave function of the condensate (macroscopic quantum state) satisfies a defocusing cubic nonlinear Schrödinger-type equation, the Gross–Pitaevskii equation. We include consequences of pair excitation, i.e., the scattering of particles in pairs from the condensate to other states, proposed in [T. T. Wu, J. Math. Phys., 2 (1961), pp. 105–123]. Our derivations rely on an uncontrolled yet physically motivated assumption for the many-body wave function. By relaxing mathematical rigor, from a particle Hamiltonian with a spatially varying interaction strength we derive via heuristics an integro–partial differential equation for the pair collision kernel, K, under a stationary condensate wave function, Φ\PhiΦ. For a scattering length with periodic microstructure of subscale ϵ\epsilonϵ, we formally describe via classical homogenization the lowest many-body bound state in terms of Φ\PhiΦ and K up to second order in ϵ\epsilonϵ. If the external potential is slowly varying, we solve the homogenized equations via boundary layer theory. As an application, we describe the partial depletion of the condensate.
In Bose-Einstein condensation, a macroscopically large number of particles occupy the same single-particle quantum state. Our goal is to study time-dependent aspects of particle excitations to states other than the single-particle macroscopic state in trapped dilute atomic gases. We adopt the view that atoms are excited in pairs so that their scattering from the single-particle state to vector positions x and y at time t is described by the pair-excitation function, K0(x,y,t). We solve a nonlocal equation for K0 under a slowly varying external potential by assuming that the wave function of the macroscopic state satisfies a time-independent nonlinear Schroedinger equation. For zero initial excitation (K0=0 at t=0) and sufficiently large t, we evaluate asymptotically K0 in terms of the one-variable Lommel function for any distance |x-y|.
At finite temperatures below the phase transition point, the Bose-Einstein condensation, the macroscopic occupation of a single quantum state by particles of integer spin, is not complete. In the language of superfluid helium, this means that the superfluid coexists with the normal fluid. Our goal is to describe this coexistence in trapped, dilute atomic gases with repulsive interactions via mean field laws that account for a {\em spatially varying} particle interaction strength. By starting with the NNN-body Hamiltonian, N≫1N\gg 1N≫1, we formally derive a system of coupled, nonlinear evolution equations in 3+13+13+1 dimensions for the following quantities: (i) the wave function of the macroscopically occupied state; and (ii) the single-particle wave functions of thermally excited states. For stationary (bound) states and a scattering length with {\em periodic microstructure} of subscale εεε, we heuristically extract effective equations of motion via periodic homogenization up to second order in εεε.
This thesis proves certain results concerning an important question in non-equilibrium quantum statistical mechanics which is the derivation of effective evolution equations approximating the dynamics of a system of large number of bosons initially at equilibrium (ground state at very low temperatures). The dynamics of such systems are governed by the time-dependent linear many-body Schroedinger equation from which it is typically difficult to extract useful information due to the number of particles being large. We will study quantitatively (i.e. with explicit bounds on the error) how a suitable one particle non-linear Schroedinger equation arises in the mean field limit as number of particles N → ∞ and how the appropriate corrections to the mean field will provide better approximations of the exact dynamics. In the first part of this thesis we consider the evolution of N bosons, where N is large, with two-body interactions of the form N³ᵝv(Nᵝ⋅), 0≤β≤1. The parameter β measures the strength and the range of interactions. We compare the exact evolution with an approximation which considers the evolution of a mean field coupled with an appropriate description of pair excitations, see [18,19] by Grillakis-Machedon-Margetis. We extend the results for 0 ≤ β < 1/3 in [19, 20] to the case of β < 1/2 and obtain an error bound of the form p(t)/Nᵅ, where α>0 and p(t) is a polynomial, which implies a specific rate of convergence as N → ∞. In the second part, utilizing estimates of the type discussed in the first part, we compare the exact evolution with the mean field approximation in the sense of marginals. We prove that the exact evolution is close to the approximate in trace norm for times of the order o(1)√N compared to log(o(1)N) as obtained in Chen-Lee-Schlein [6] for the Hartree evolution. Estimates of similar type are obtained for stronger interactions as well.
Quantum systems consisting of many particles are both interesting from the physics and the mathematics perspective. One aspect that makes them fascinating from the physics point of view is that they can exhibit quantum effects on macroscopic scales such as Bose-Einstein condensation. This phenomenon, first predicted 100 years ago by Bose [19] and Einstein [44], occurs when a dilute gas of bosons is cooled down to very low temperatures. Below a certain critical temperature a condensate forms, meaning that a macroscopic fraction of the particles is in the same one-particle state, and thus, quantum phenomena become relevant on length scales much larger than the individual atom. The occurrence of Bose-Einstein condensation is also related to superfluidity as suggested in [88]. Bose-Einstein condensation has been achieved experimentally independently by the groups of Cornell and Wieman [3] and Ketterle [40] earning them a Nobel prize in 2001. While Bose and Einstein analyzed theoretically the non-interacting gas, an effective theory describing weakly interacting Bose gases was proposed by Bogoliubov in 1947 with the goal of explaining the occurrence of superfluidity in liquid helium [18]. Using assumptions such as that condensation indeed occurs, his theory yields predictions for the low-energy states of the system. Seen from the mathematics perspective, describing the Bose gas is an interesting challenge due to the strong correlations between particles introduced by the interaction making this case much harder to study than the non-interacting one. In particular, a proof that Bose-Einstein condensation occurs in the thermodynamic limit for continuous systems is still out of reach. However, there has been a lot of work regarding the rigorous analysis of bosonic systems in various relevant regimes. This thesis contains three new contributions establishing rigorous results on bosonic systems in different, physically relevant settings. It is based on the articles [37, 27, 31]. The first result concerns the time-evolution of initially trapped Bose gases in the ultra-dilute Gross-Pitaevskii regime. Under some assumptions on the initial data satisfied in particular by a class of states that are expected to approximate the ground state in a trap, we derive a quasi-free approximation that is close to the true solution in L2L^2L2-norm. In particular, we show that this approximate wave function is described in terms of a Bose-Einstein condensate in a one-particle state following the time-dependent Gross-Pitaevskii equation together with pair excitations with respect to this condensate. As a corollary, we establish a central limit theorem for bounded one-particle observables. The second result is a simple proof of Bose-Einstein condensation for interacting, ultra-dilute gases. More specifically, we consider scalings interpolating between the Gross-Pitaevskii regime and the thermodynamic limit. While condensation had been shown to hold in these regimes before, the advantage of the method used here is its simplicity. Lastly, the third result establishes an upper bound on the ground state energy of a dilute Bose gas in the thermodynamic limit. It captures rigorously the third order correction to the energy per unit volume for sufficiently low densities as predicted in the physics literature first by Wu and subsequently others. Full-text embargoed until: 2026-01-01\textit{Full-text embargoed until: 2026-01-01}Full-text embargoed until: 2026-01-01