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

arXiv:2511.09206 · Journal of Statistical Mechanics Theory and Experiment

Mean-field theory of the DNLS equation at positive and negative absolute temperatures

Michele Giusfredi, Stefano Iubini, Antonio Politi, Paolo Politi

Published May 1, 2026• 0 citations
arXivHTMLPDF

Abstract

Abstract The discrete nonlinear Schrödinger (DNLS) model, owing to the existence of two conserved quantities, exhibits an equilibrium transition from a homogeneous phase at positive absolute temperature to a localized phase at negative absolute temperature. Here, we provide a mean-field (MF) theory of DNLS through a suitable approximation of the grand canonical partition function, which makes it factorizable and can be used to describe the equilibrium state at positive temperatures and the metastable state at negative temperatures. Comparison of our MF results with numerically exact calculations shows that this approximation is good to excellent in the entire grand canonical phase diagram. Explicit approximate expressions for equilibrium observables are provided in the high-temperature limit. Our theory represents a clear advancement over the model that neglects the interaction between sites.

DOI
10.1088/1742-5468/ae625e

Related papers

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.

arXiv:2001.01546•Journal of the American Mathematical Society•Jun 16, 2021•15 citations

Mean Field Derivation of DNLS from the Bose–Hubbard Model

E. Picari, Antonio Ponno, Lorenzo Zanelli

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.

Annales Henri Poincaré•Oct 1, 2021•1 citation

Equilibrium temperatures of discrete nonlinear systems

Uri Levy, Yaron R. Silberberg

During the evolution of discrete nonlinear systems with dynamics dictated by the discrete nonlinear Schr\"odinger equation, two quantities are conserved: system energy (Hamiltonian) and system density (number of particles). It is then possible to analyze system evolution in relation to an energy-density phase diagram. Previous works have identified a ``thermalization zone'' on the phase diagram where regular statistical mechanics methods apply. Based on these statistical mechanics methods we have now assigned a specific equilibrium temperature to every point of the thermalization zone. Temperatures were derived in the grand canonical picture through an entropy-temperature relation, modified to suit the nonlinear lattice systems. Generally, everywhere in the thermalization zone of the phase diagram, temperatures along a fixed system-density line, grow monotonously from zero to infinity. Isotherms on the phase diagram are concave.

Physical review. B./Physical review. B•Aug 30, 2018•12 citations

Low temperature dynamics of the one-dimensional discrete nonlinear Schrödinger equation

Christian B. Mendl, Herbert Spohn

We study equilibrium time correlations for the discrete nonlinear Schrödinger equation on a one-dimensional lattice and unravel three dynamical regimes. There is a high temperature regime with density and energy as the only two conserved fields. Their correlations have zero velocity and spread diffusively. In the low temperature regime umklapp processes are rare with the consequence that phase differences appear as an additional (almost) conserved field. In an approximation where all umklapp is suppressed, while the equilibrium state remains untouched, one arrives at an anharmonic chain. Using the method of nonlinear fluctuating hydrodynamics we establish that the DNLS equilibrium time correlations have the same signature as a generic anharmonic chain, in particular KPZ broadening for the sound peaks and Lévy 5/3 broadening for the heat peak. In the, so far not sharply defined, ultra-low temperature regime the integrability of the dynamics becomes visible. As an illustration we simulate the completely integrable Ablowitz–Ladik model and confirm ballistic broadening of the time correlations.

arXiv:1505.04218•Journal of Statistical Mechanics Theory and Experiment•Aug 26, 2015•37 citations

Even-Odd Effect on the Thermal Equilibrium State of the Pairing Model. II: Mean Field Approximation and Renormalized Distribution

A. Kuriyama, João da Providência, Y. Tsue +1 more

We define a mean field approximation of the modified grand partition function of the pairing model, which has been defined on the modified grand canonical ensemble with even or odd particle number. Using numerical evaluations, it is shown that this approximation reproduces the even-odd effect exhibited in the case of the modified grand partition function at low temperature very well in a qualitative sense. The expectation values of one-body operators can be expressedsimply in terms of those obtainedwith the conventional mean field approximation when a Fermi-type distribution is replaced by a renormalized distribution. Although Wick’s theorm does not hold in this case, the dominant parts of the expectation values of two-body operators are also given using those obtained by the conventional mean fieldapproximation with the replacement mentionedabove.

Progress of Theoretical Physics•Sep 1, 2002•0 citations

Pairing Model and Mixed State Representation. II: Grand Partition Function and Its Mean Field Approximation

A. Kuriyama, João da Providência, Y. Tsue +1 more

In order to compare the approach of the mixed state representation of the pairing model with the exact evaluation of the grand partition function, we introduce the fermion space of a grand canonical ensemble. For the definition of the grand partition function, we adopt a quasi-spin space and evaluate its multiplicity. then we show how to derive the mixed state representation, which yields results that are approximately the same as those of the thermal Hartree-Fock-Bogoliubov theory for the thermal equilibrium state, from the grand partition function under the mean field approximation. We compare numerical results of the mean field approximation with those obtained from the exact evaluation of the grand partition function. We find that the mean field approximation gives good results. However, we find a distinctive difference between the two sets of results in the condensed phase.

Progress of Theoretical Physics•Jan 1, 2002•1 citation
Abstract The discrete nonlinear Schrödinger (DNLS) model, owing to the existence of two conserved quantities, exhibits an equilibrium transition from a homogeneous phase at positive absolute temperature to a localized phase at negative absolute temperature. Here, we provide a mean-field (MF) theory of DNLS through a suitable approximation of the grand canonical partition function, which makes it factorizable and can be used to describe the equilibrium state at positive temperatures and the metastable state at negative temperatures. Comparison of our MF results with numerically exact calculations shows that this approximation is good to excellent in the entire grand canonical phase diagram. Explicit approximate expressions for equilibrium observables are provided in the high-temperature limit. Our theory represents a clear advancement over the model that neglects the interaction between sites.
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.
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 LpL^pLp L p -estimate, for any p≥1p \ge 1p≥1 p ≥ 1 , is obtained with a linear dependence on time due to a Gaussian measure on initial data coherent states.
During the evolution of discrete nonlinear systems with dynamics dictated by the discrete nonlinear Schr\"odinger equation, two quantities are conserved: system energy (Hamiltonian) and system density (number of particles). It is then possible to analyze system evolution in relation to an energy-density phase diagram. Previous works have identified a ``thermalization zone'' on the phase diagram where regular statistical mechanics methods apply. Based on these statistical mechanics methods we have now assigned a specific equilibrium temperature to every point of the thermalization zone. Temperatures were derived in the grand canonical picture through an entropy-temperature relation, modified to suit the nonlinear lattice systems. Generally, everywhere in the thermalization zone of the phase diagram, temperatures along a fixed system-density line, grow monotonously from zero to infinity. Isotherms on the phase diagram are concave.
We study equilibrium time correlations for the discrete nonlinear Schrödinger equation on a one-dimensional lattice and unravel three dynamical regimes. There is a high temperature regime with density and energy as the only two conserved fields. Their correlations have zero velocity and spread diffusively. In the low temperature regime umklapp processes are rare with the consequence that phase differences appear as an additional (almost) conserved field. In an approximation where all umklapp is suppressed, while the equilibrium state remains untouched, one arrives at an anharmonic chain. Using the method of nonlinear fluctuating hydrodynamics we establish that the DNLS equilibrium time correlations have the same signature as a generic anharmonic chain, in particular KPZ broadening for the sound peaks and Lévy 5/3 broadening for the heat peak. In the, so far not sharply defined, ultra-low temperature regime the integrability of the dynamics becomes visible. As an illustration we simulate the completely integrable Ablowitz–Ladik model and confirm ballistic broadening of the time correlations.
We define a mean field approximation of the modified grand partition function of the pairing model, which has been defined on the modified grand canonical ensemble with even or odd particle number. Using numerical evaluations, it is shown that this approximation reproduces the even-odd effect exhibited in the case of the modified grand partition function at low temperature very well in a qualitative sense. The expectation values of one-body operators can be expressedsimply in terms of those obtainedwith the conventional mean field approximation when a Fermi-type distribution is replaced by a renormalized distribution. Although Wick’s theorm does not hold in this case, the dominant parts of the expectation values of two-body operators are also given using those obtained by the conventional mean fieldapproximation with the replacement mentionedabove.
In order to compare the approach of the mixed state representation of the pairing model with the exact evaluation of the grand partition function, we introduce the fermion space of a grand canonical ensemble. For the definition of the grand partition function, we adopt a quasi-spin space and evaluate its multiplicity. then we show how to derive the mixed state representation, which yields results that are approximately the same as those of the thermal Hartree-Fock-Bogoliubov theory for the thermal equilibrium state, from the grand partition function under the mean field approximation. We compare numerical results of the mean field approximation with those obtained from the exact evaluation of the grand partition function. We find that the mean field approximation gives good results. However, we find a distinctive difference between the two sets of results in the condensed phase.