Annales Henri Poincaré

Mean Field Derivation of DNLS from the Bose–Hubbard Model

E. Picari, Antonio Ponno, Lorenzo Zanelli

Published Oct 1, 2021• 1 citation

Abstract

Abstract We prove that the flow of the discrete nonlinear Schrödinger equation (DNLS) is the mean field limit of the quantum dynamics of the Bose–Hubbard model for N interacting particles. In particular, we show that the Wick symbol of the annihilation operators evolved in the Heisenberg picture converges, as N becomes large, to the solution of the DNLS. A quantitative $$L^p$$ L p -estimate, for any $$p \ge 1$$ p ≥ 1 , is obtained with a linear dependence on time due to a Gaussian measure on initial data coherent states.

DOI: 10.1007/s00023-021-01112-6

Related papers

Mean-Field Dynamics of the Bose–Hubbard Model in High Dimension

Shahnaz Farhat, Denis Périce, Sören Petrat

arXiv:2501.05304Annales Henri PoincaréOct 10, 20260 citations

The mean-field limit of quantum Bose gases at positive temperature

Jürg Martin Fröhlich, Antti Knowles, Benjamin Schlein +1 more

arXiv:2001.01546Journal of the American Mathematical SocietyJun 16, 202115 citations

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

David Mitrouskas

Jan 1, 20173 citations

Mean-field approximation for a Bose–Hubbard dimer with complex interaction strength

Eva-Maria Graefe, C. Liverani

arXiv:1305.7160Journal of Physics A Mathematical and TheoreticalOct 24, 20136 citations

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

Elif Kuz

Jan 1, 20160 citations

Existence, dynamics and mobility of Quantum Compactons in an extended Bose-Hubbard model

Peter Jason, Magnus Johansson

Jan 1, 20130 citations

Annales Henri Poincaré

Mean Field Derivation of DNLS from the Bose–Hubbard Model

E. Picari, Antonio Ponno, Lorenzo Zanelli

Published Oct 1, 2021• 1 citation

Publisher

Abstract

L^p

L p -estimate, for any

p \ge 1

p ≥ 1 , is obtained with a linear dependence on time due to a Gaussian measure on initial data coherent states.

DOI: 10.1007/s00023-021-01112-6

Related papers

Mean-Field Dynamics of the Bose–Hubbard Model in High Dimension

Shahnaz Farhat, Denis Périce, Sören Petrat

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.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.

arXiv:2501.05304Annales Henri PoincaréOct 10, 20260 citations

The mean-field limit of quantum Bose gases at positive temperature

Jürg Martin Fröhlich, Antti Knowles, Benjamin Schlein +1 more

We prove that the grand canonical Gibbs state of an interacting quantum Bose gas converges to the Gibbs measure of a nonlinear Schrödinger equation in the mean-field limit, where the density of the gas becomes large and the interaction strength is proportional to the inverse density. Our results hold in dimensions d ⩽ 3 d \leqslant 3 . For d > 1 d > 1 the Gibbs measure is supported on distributions of negative regularity and we have to renormalize the interaction. More precisely, we prove the convergence of the relative partition function and of the reduced density matrices in the L r L^r -norm with optimal exponent r r . Moreover, we prove the convergence in the L ∞ L^\infty -norm of Wick-ordered reduced density matrices, which allows us to control correlations of Wick-ordered particle densities as well as the asymptotic distribution of the particle number. Our proof is based on a functional integral representation of the grand canonical Gibbs state, in which convergence to the mean-field limit follows formally from an infinite-dimensional stationary phase argument for ill-defined non-Gaussian measures. We make this argument rigorous by introducing a white-noise-type auxiliary field, through which the functional integral is expressed in terms of propagators of heat equations driven by time-dependent periodic random potentials and can, in turn, be expressed as a gas of interacting Brownian loops and paths. When the gas is confined by an external trapping potential, we control the decay of the reduced density matrices using excursion probabilities of Brownian bridges.We prove that the grand canonical Gibbs state of an interacting quantum Bose gas converges to the Gibbs measure of a nonlinear Schrödinger equation in the mean-field limit, where the density of the gas becomes large and the interaction strength is proportional to the inverse density. Our results hold in dimensions d ⩽ 3 d \leqslant 3 . For d > 1 d > 1 the Gibbs measure is supported on distributions of negative regularity and we have to renormalize the interaction. More precisely, we prove the convergence of the relative partition function and of the reduced density matrices in the L r L^r -norm with optimal exponent r r . Moreover, we prove the convergence in the L ∞ L^\infty -norm of Wick-ordered reduced density matrices, which allows us to control correlations of Wick-ordered particle densities as well as the asymptotic distribution of the particle number. Our proof is based on a functional integral representation of the grand canonical Gibbs state, in which convergence to the mean-field limit follows formally from an infinite-dimensional stationary phase argument for ill-defined non-Gaussian measures. We make this argument rigorous by introducing a white-noise-type auxiliary field, through which the functional integral is expressed in terms of propagators of heat equations driven by time-dependent periodic random potentials and can, in turn, be expressed as a gas of interacting Brownian loops and paths. When the gas is confined by an external trapping potential, we control the decay of the reduced density matrices using excursion probabilities of Brownian bridges.

arXiv:2001.01546Journal of the American Mathematical SocietyJun 16, 202115 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.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, 20173 citations

Mean-field approximation for a Bose–Hubbard dimer with complex interaction strength

Eva-Maria Graefe, C. Liverani

In the limit of large particle numbers and low densities systems of cold atoms can be effectively described as macroscopic single-particle systems in a mean-field approximation. In the case of a Bose–Hubbard system, modelling bosons on a discrete lattice with on-site interactions, this yields a discrete nonlinear Schrödinger equation of Gross–Pitaevskii type. It has been recently shown that the correspondence between the Gross–Pitaevskii equation and the Bose–Hubbard system breaks down for complex extensions. In particular, for a Bose–Hubbard dimer with complex on-site energy the mean-field approximation yields a generalized complex nonlinear Schrödinger equation. Conversely, a Gross–Pitaevskii equation with complex on-site energies arises as the mean-field approximation of many-particle Lindblad dynamics rather than a complex extension of the Bose–Hubbard system. Here we address the question of how the mean-field description is modified in the presence of a complex-valued particle interaction term for a Bose–Hubbard dimer. We derive the mean-field equations of motion leading to nonlinear dissipative Bloch dynamics, related to a nontrivial complex generalization of the nonlinear Schrödinger equation. The resulting dynamics are analysed in detail. It is shown that depending on the parameter values there can be up to six stationary states, and for small values of the interaction strength there are limit cycles. Furthermore, we show how a Gross–Pitaevskii equation with a complex interaction term can be derived as the mean-field approximation of a Bose–Hubbard dimer with an additional Lindblad term modelling two-particle losses.In the limit of large particle numbers and low densities systems of cold atoms can be effectively described as macroscopic single-particle systems in a mean-field approximation. In the case of a Bose–Hubbard system, modelling bosons on a discrete lattice with on-site interactions, this yields a discrete nonlinear Schrödinger equation of Gross–Pitaevskii type. It has been recently shown that the correspondence between the Gross–Pitaevskii equation and the Bose–Hubbard system breaks down for complex extensions. In particular, for a Bose–Hubbard dimer with complex on-site energy the mean-field approximation yields a generalized complex nonlinear Schrödinger equation. Conversely, a Gross–Pitaevskii equation with complex on-site energies arises as the mean-field approximation of many-particle Lindblad dynamics rather than a complex extension of the Bose–Hubbard system. Here we address the question of how the mean-field description is modified in the presence of a complex-valued particle interaction term for a Bose–Hubbard dimer. We derive the mean-field equations of motion leading to nonlinear dissipative Bloch dynamics, related to a nontrivial complex generalization of the nonlinear Schrödinger equation. The resulting dynamics are analysed in detail. It is shown that depending on the parameter values there can be up to six stationary states, and for small values of the interaction strength there are limit cycles. Furthermore, we show how a Gross–Pitaevskii equation with a complex interaction term can be derived as the mean-field approximation of a Bose–Hubbard dimer with an additional Lindblad term modelling two-particle losses.

arXiv:1305.7160Journal of Physics A Mathematical and TheoreticalOct 24, 20136 citations

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.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, 20160 citations

Existence, dynamics and mobility of Quantum Compactons in an extended Bose-Hubbard model

Peter Jason, Magnus Johansson

Lattice Compactons, discrete breathers with compact support, were found for a discrete nonlinear Schrodinger (DNLS) equation extended with nearest neighbour intersite nonlinearities [1], a model originally studied with waveguide arrays in mind. These compactons were shown to exhibit very good mobility if the parameters are tuned close to the compactons stability boundary. The DNLS can also be used to model the behaviour of Bose-Einstein condensates in optical lattices, and the remarkable control over the experiments in this field of research has made it possible to study the quantum mechanics of strongly correlated atoms.We will define the concept of a Quantum Lattice Compacton [2] and discuss the existence and dynamics, with special emphasis on mobility [3], of these in an extended Bose-Hubbard model corresponding to above-mentioned extended DNLS equation in the quantum mechanical limit. The compactons is given 'a kick' by means of a phase-gradient and it is shown that the size of this phase is crucial for the mobility of the compactons. For small phase-gradients, corresponding to a slow coherent motion in the classical model, the time-scales of the quantum tunnelings become of the same order as the time-scale of the translational motion and the classical mobility is destroyed by quantum fluctuations. For large phase-gradients, corresponding to rapid classical motion, the classical and quantum time-scales separate so that a mobile, localized coherent quantum state can be translated many sites in the lattice already for small particle numbers of the order of 10 [3].Acknowledgements: This project has been financed by the Swedish Research Council.References[1] M. Oster, M. Johansson, and A. Eriksson 2003 Phys. Rev. E 67 056606[2] P. Jason and M. Johansson 2012 Phys. Rev. A 85 011603(R)[3] P. Jason and M. Johansson 2013 Phys. Rev. A 88 033605Lattice Compactons, discrete breathers with compact support, were found for a discrete nonlinear Schrodinger (DNLS) equation extended with nearest neighbour intersite nonlinearities [1], a model originally studied with waveguide arrays in mind. These compactons were shown to exhibit very good mobility if the parameters are tuned close to the compactons stability boundary. The DNLS can also be used to model the behaviour of Bose-Einstein condensates in optical lattices, and the remarkable control over the experiments in this field of research has made it possible to study the quantum mechanics of strongly correlated atoms.We will define the concept of a Quantum Lattice Compacton [2] and discuss the existence and dynamics, with special emphasis on mobility [3], of these in an extended Bose-Hubbard model corresponding to above-mentioned extended DNLS equation in the quantum mechanical limit. The compactons is given 'a kick' by means of a phase-gradient and it is shown that the size of this phase is crucial for the mobility of the compactons. For small phase-gradients, corresponding to a slow coherent motion in the classical model, the time-scales of the quantum tunnelings become of the same order as the time-scale of the translational motion and the classical mobility is destroyed by quantum fluctuations. For large phase-gradients, corresponding to rapid classical motion, the classical and quantum time-scales separate so that a mobile, localized coherent quantum state can be translated many sites in the lattice already for small particle numbers of the order of 10 [3].Acknowledgements: This project has been financed by the Swedish Research Council.References[1] M. Oster, M. Johansson, and A. Eriksson 2003 Phys. Rev. E 67 056606[2] P. Jason and M. Johansson 2012 Phys. Rev. A 85 011603(R)[3] P. Jason and M. Johansson 2013 Phys. Rev. A 88 033605

Jan 1, 20130 citations

L^p

L p -estimate, for any

p \ge 1

p ≥ 1 , is obtained with a linear dependence on time due to a Gaussian measure on initial data coherent states.

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.