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.
In these notes, based on a mini-course held at the summer school "Methods and Models of Kinetic Theory" that took place in Porto Ercole in June 2016, we review some of the recent developments in the derivation of effective evolution equations starting from many-body quantum mechanics. We discuss the derivation of the Hartree equation in the bosonic mean-field limit, of the Gross-Pitaevskii equation describing the dynamics of initially trapped Bose-Einstein condensates and of the Hartree-Fock equation for fermions in a joint mean-field and semiclassical limit.
Mean-field Hartree theory is a central tool for mapping interacting many-body dynamics into an effective nonlinear one-particle evolution. This approximation has been also employed when the Hamiltonian that governs the many-body dynamics is not Hermitian. Indeed, non-Hermitian Hamiltonians model particle gain/loss or the evolution of open quantum systems between consecutive quantum jumps. Furthermore, the validity of the Hartree approximation for generic non-Hermitian Hamiltonians lies at the basis of a quantum algorithm for nonlinear differential equations, where the nonlinear equation is mapped into a many-body Schr\"odinger equation, which may be non-Hermitian. In this application, the mean-field approximation can be seen as a linearization method, and at the end of the algorithm, the obtained linear system is solved using the Harrow-Hassidim-Lloyd (HHL) algorithm. The non-Hermitian mean-field Hartree approximation is far from being completely understood, and no rigorous results justify its applications. Instead, it is usually presented as a heuristic extension of the Hermitian case. In the first part of this thesis, we show that this approximation can fail. We consider the most general bosonic mean-field Hamiltonian with two-body interactions and we prove that, if the Hamiltonian is not Hermitian, the Hartree equation does not capture the correct second time derivative of the one-particle marginal state at
t=0 in the limit of infinitely many particles. We construct a series of examples that we can solve analytically or numerically, which support the failure of the Hartree approximation and give further insight into the new phenomena that arise in non-Hermitian physics. Furthermore, we consider the mean-field evolution of open quantum systems and we show how to use the mean-field approximation of open quantum systems as a linearization method of nonlinear differential equations. In the second part of this thesis, we study the mean-field dynamics of a Hermitian system of
N interacting bosons starting from an initially condensed state. For a broad class of mean-field Hamiltonians, including models with arbitrary bounded interactions and models with unbounded interaction potentials, we prove that the probability of having
n particles outside the condensate decays exponentially in
n for any finite evolution time. Our results strengthen previously known bounds that provide only polynomial control on the probability of having
n excitations. Finally, we consider the HHL algorithm. Several circuit-level variants have been proposed to simplify its implementation, including constructions in which the clock register is factorized via Hadamard gates to reduce circuit depth and entanglement. We provide a rigorous analysis of the convergence of the output towards the ideal solution, and we show that this commonly deployed variant can fundamentally alter the asymptotic convergence of the algorithm. This phenomenon is not a finite-size artifact and compromises the correctness guarantees of the factorized-clock implementation. We provide an analytic characterization of the underlying mechanism and propose a minimal and experimentally feasible modification that suffices to restore convergence while preserving the practical advantages of the simplified circuit. Our results clarify an overlooked correctness issue in a widely used HHL implementation pattern and delineate the conditions under which simplified realizations faithfully reproduce the intended algorithmic behavior.
We consider the time evolution of N bosons in the mean field regime for factorized initial data. In the limit of large N, the many body evolution can be approximated by the non-linear Hartree equation. In this paper we are interested in the fluctuations around the Hartree dynamics. We choose k self-adjoint one-particle operators O 1,…,O k on
L2(R3) , and we average their action over the N-particles. We show that, for every fixed
t∈R , expectations of products of functions of the averaged observables approach, as N→∞, expectations with respect to a complex Gaussian measure, whose covariance matrix can be expressed in terms of a Bogoliubov transformation describing the dynamics of quantum fluctuations around the mean field Hartree evolution. If the operators O 1,…,O k commute, the Gaussian measure is real and positive, and we recover a "classical" multivariate central limit theorem. All our results give explicit bounds on the rate of the convergence.
In these notes, we review some recent mathematical results concerning the derivation of effective evolution equations from many body quantum mechanics. In particular, we discuss the emergence of the Hartree equation in the so-called mean field regime (for example, for systems of gravitating bosons), and we show that the Gross-Pitaevskii equation approximates the dynamics of initially trapped Bose-Einstein condensates. We explain how effective evolution equations can be derived, on the one hand, by analyzing the so called BBGKY hierarchy, describing the time-evolution of reduced density matrices, and, on the other hand, by studying the dynamics of coherent initial states in a Fock-space representation of the many body system.