We consider the time evolution of quantum states by many-body Schrödinger dynamics and study the rate of convergence of their reduced density matrices in the mean field limit. If the prepared state at initial time is of coherent or factorized type and the number of particles
n is large enough then it is known that
1/n is the correct rate of convergence at any time. We show in the simple case of bounded pair potentials that the previous rate of convergence holds in more general situations with possibly correlated prepared states. In particular, it turns out that the coherent structure at initial time is unessential and the important fact is rather the speed of convergence of all reduced density matrices of the prepared states. We illustrate our result with several numerical simulations and examples of multi-partite entangled quantum states borrowed from quantum information.
Abstract We introduce the number-conserving quantum phase space description as a versatile tool to address fundamental aspects of quantum many-body systems. Using phase space methods we prove two alternative versions of the quantum de Finetti theorem for finite-dimensional bosonic quantum systems, which states that a reduced density matrix of a many-body quantum state can be approximated by a convex combination of product states where the error is proportional to the inverse particle number. This theorem provides a formal justification for the mean-field description of many-body quantum systems, as it shows that quantum correlations can be neglected for the calculation of few-body observables when the particle number is large. Furthermore we discuss methods to derive the exact evolution equations for quantum phase space distribution functions as well as upper and lower bounds for the ground state energy. As an important example, we consider the Bose–Hubbard model and show that the mean-field dynamics is given by a classical phase space flow equivalent to the discrete Gross–Pitaevskii equation.
Mean-field Hartree theory is a central tool for reducing many-body dynamics to an effective one-particle evolution. This approximation has been employed when the Hamiltonian governing 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. In this work, we show that this approximation can fail.
First, we consider the most general bosonic mean-field Hamiltonian with two-body interactions and we prove that, for non-Hermitian Hamiltonian, the Hartree equation does not capture the correct second time derivative of the one-particle marginal state at
t=0 in the thermodynamic limit. Then, we analytically solve a model of
N bosonic qubits with a purely anti-Hermitian Hamiltonian, compute the limit for
N→∞ of the one-particle marginal state and show that it does not agree with the solution of the non-Hermitian Hartree equation. We further show that there exists an initial condition such that the exact one-particle marginal state undergoes a finite-time transition to a mixed state, a phenomenon that is completely absent in the case of Hermitian Hamiltonians. Finally, we simulate a non-Hermitian Hamiltonian with real spectrum and show that, also in this case, there are initial conditions for which the Hartree equation does not reproduce the correct one-particle marginal state in the limit
N→∞.
Our findings challenge the validity of the mean-field Hartree approximation for non-Hermitian Hamiltonians, and call for additional conditions for the validity of the mean-field regime to model the dynamics of particle gain and loss and the open-system dynamics in bosonic many-body systems.
Abstract This paper is a first step toward rigorously establishing the validity of mean-field approximations for the dynamics of quantum systems in a large dimensional limit. Our result is a derivation of the large d mean-field limit of the Bose–Hubbard model on a square lattice. The Bose–Hubbard model effectively describes bosons on a lattice with on-site interactions and nearest-neighbour hopping, serving as a foundational framework for understanding strong particle interactions and the superfluid to Mott insulator phase transition. We prove a trace norm estimate between the one-lattice-site reduced density of the Schrödinger dynamics and the mean-field dynamics in the limit of large dimension. Here, the mean-field approximation is in the hopping amplitude and not in the interaction, leading to a very rich and non-trivial mean-field equation. This mean-field equation does not only describe the condensate, as is the case when the mean-field description comes from a large particle number limit averaging out the interaction, but it allows for a phase transition since it contains the full non-trivial interaction. Our work is a rigorous justification of a simple case of the highly successful dynamical mean-field theory (DMFT) for bosons, which somewhat surprisingly yields many qualitatively correct results in three dimensions.
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 a mean-field model to describe the dynamics of [Formula: see text] bosons of species one and [Formula: see text] bosons of species two in the limit as [Formula: see text] and [Formula: see text] go to infinity. We embed this model into Fock space and use it to describe the time evolution of coherent states which represent two-component condensates. Following this approach, we obtain a microscopic quantum description for the dynamics of such systems, determined by the Schrödinger equation. Associated to the solution to the Schrödinger equation, we have a reduced density operator for one particle in the first component of the condensate and one particle in the second component. In this paper, we estimate the difference between this operator and the projection onto the tensor product of two functions that are solutions of a system of equations of Hartree type. Our results show that this difference goes to zero as [Formula: see text] and [Formula: see text] go to infinity.