How and why could an interacting system of many particles be described as if all particles were independent and identically distributed? This question is at least as old as statistical mechanics itself. Its quantum version has been rejuvenated by the birth of cold atoms physics. In particular, the experimental creation of Bose–Einstein condensates leads to the following variant: why and how can a large assembly of very cold interacting bosons (quantum particles deprived of the Pauli exclusion principle) all populate the same quantum state? In this text I review the various mathematical techniques allowing to prove that the lowest energy state of a bosonic system forms, in a reasonable macroscopic limit of large particle number, a Bose–Einstein condensate. This means that indeed in the relevant limit all particles approximately behave as if independent and identically distributed, according to a law determined by minimizing a non-linear Schrödinger energy functional. This is a particular instance of the justification of the mean-field approximation in statistical mechanics, starting from the basic many-body Schrödinger Hamiltonian.
How and why could an interacting system of many particles be described as if all particles were independent and identically distributed? This question is at least as old as statistical mechanics itself. Its quantum version has been rejuvenated by the birth of cold atoms physics. In particular, the experimental creation of Bose–Einstein condensates leads to the following variant: why and how can a large assembly of very cold interacting bosons (quantum particles deprived of the Pauli exclusion principle) all populate the same quantum state? In this text I review the various mathematical techniques allowing to prove that the lowest energy state of a bosonic system forms, in a reasonable macroscopic limit of large particle number, a Bose–Einstein condensate. This means that indeed in the relevant limit all particles approximately behave as if independent and identically distributed, according to a law determined by minimizing a non-linear Schrödinger energy functional. This is a particular instance of the justification of the mean-field approximation in statistical mechanics, starting from the basic many-body Schrödinger Hamiltonian.
These notes deal with the mean-field approximation for equilibrium states of N-body systems in classical and quantum statistical mechanics. A general strategy for the justification of effective models based on statistical independence assumptions is presented in details. The main tools are structure theorems {à} la de Finetti, describing the large N limits of admissible states for these systems. These rely on the symmetry under exchange of particles, due to their indiscernability. Emphasis is put on quantum aspects, in particular the mean-field approximation for the ground states of large bosonic systems, in relation with the Bose-Einstein condensation phenomenon. Topics covered in details include: the structure of reduced density matrices for large bosonic systems, Fock-space localization methods, derivation of effective energy functionals of Hartree or non-linear Schr{ö}dinger type, starting from the many-body Schr{ö}dinger Hamiltonian.
These lecture notes treat the mean-field approximation for equilibrium states of N body systems in classical and quantum statistical mechanics. A general strategy to justify effective models based on assumptions of statistical independence of the particles is in presented in detail. The main tools are a structure theorems of de Finetti that describe large N limits of states accessible to the systems in question, exploiting the indistinguishablity of particles. The focus is on quantum aspects, particularly the mean-field approximation for the ground state of a large system of bosons, in connection with Bose-Einstein condensation: structure of reduced density matrices of a large bosonic system, localization methods in Fock space, derivation of Hartree and non-linear Schrödinger effective energy functionals.
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 L2-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
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.
This thesis is devoted to the rigorous study of dilute Bose gases with two-body and three-body interactions at zero temperature in two and three dimensions. Bose gases are many-body quantum systems composed of bosons, which notably do not obey the Pauli exclusion principle. As a consequence, at extremely low temperatures, a macroscopic fraction of the particles can occupy the same one-body quantum state. This phenomenon, known as Bose–Einstein condensation, was predicted in the 1920s by Bose and Einstein for noninteracting Bose gases, and first realised experimentally by the groups of Cornell and Wieman, and Ketterle, earning them the Nobel Prize for Physics in 2001. However, proving the occurrence of Bose–Einstein condensation in general interacting bosonic systems remains a major open problem and constitutes a central theme in this thesis. Another important property of bosonic matter is that it is generally more unstable than fermionic matter. In particular, for systems interacting via the Coulomb force, the Pauli exclusion principle is essential for explaining the stability of the second kind, which states that the energy and volume scale linearly with the number of particles. Nonetheless, bosonic systems with attractive but less singular interactions can also be stable, and the rigorous proof of this fact forms a second focus of this thesis. In three dimensions, we consider predominantly repulsive interactions. Firstly, we prove a leading order upper bound on the thermodynamic ground state energy of Bose gases with hard-core three-body interactions, which matches an existing lower bound. In addition, we obtain a simpler proof for integrable potentials. Secondly, we derive the leading order of the thermodynamic ground state energy of Bose gases with combined two-body and three-body interactions. More specifically, we show that, at the leading order, the ground state energy is entirely characterised by the scattering length of the two-body potential and the scattering hypervolume of the three-body potential. Thirdly, we study the time evolution of Bose–Einstein condensates with three-body interactions in the so-called Gross–Pitaevskii regime. We show that Bose–Einstein condensation is preserved under many-body evolution and that the condensate wavefunction evolves according to the quintic Gross–Pitaevskii equation. In particular, we show that the effective coupling constant is universal in the sense that it depends only on the scattering hypervolume of the interaction potential. As a by-product of our proof, we derive a stability of the second kind result for Bose gases with attractive two-body interactions and repulsive three-body interactions. In two dimensions, we consider trapped Bose gases with mean-field-type interactions, including the purely attractive case. We prove the stability of the second kind of the system and the convergence of the ground state energy per particle to that of a non-linear Schrödinger energy functional. As a consequence of the stability of the second kind, we also obtain Bose–Einstein condensation for the many-body ground states.
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.