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:2101.03878

Characteristics of Equilibrated Nonlinear Oscillator Systems

Uri Levy

Published Jan 11, 2021• 0 citations
arXivHTMLPDF

Abstract

During the evolution of coupled nonlinear oscillators on a lattice, with dynamics dictated by the discrete nonlinear Schrödinger equation (DNLSE systems), two quantities are conserved: system energy (Hamiltonian) and system density (number of particles). If the number of system oscillators is large enough, a significant portion of the array can be considered to be an "open system", in intimate energy and density contact with a "bath" - the rest of the array. Thus, as indicated in previous works, the grand canonical formulation can be exploited in order to determine equilibrium statistical properties of thermalized DNLSE systems. In this work, given the values of the two conserved quantities, we have calculated the necessary values of the two Lagrange parameters (typically designated $β,μ$) associated with the grand canonical partition function in two different ways. One is numerical and the other is analytic, based on a published approximate entropy expression. In addition we have accessed a purposely-derived approximate PDF expression of site-densities. Applying these mathematical tools we have generated maps of temperatures, chemical potentials, and field correlations for DNLSE systems over the entire thermalization zone of the DNLSE phase diagram, subjected to all system-nonlinearity levels. The end result is a rather complete picture, characterizing equilibrated large DNLSE systems.

DOI
10.48550/arxiv.2101.03878

Related papers

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

Theoretical and Numerical Studies of the Dynamics of Open Quantum Systems

Souvik Agasti

Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.

JYU dissertations•Jan 1, 2020•0 citations

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

Michele Giusfredi, Stefano Iubini, Antonio Politi +1 more

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.

arXiv:2511.09206•Journal of Statistical Mechanics Theory and Experiment•May 1, 2026•0 citations

Grand canonical finite-size numerical approaches: A route to measuring bulk properties in an applied field

Chisa Hotta, Naokazu Shibata

We exploit a prescription to observe directly the physical properties of the thermodynamic limit in a continuously applied field in one-dimensional quantum finite lattice systems. By systematically scaling down the energy of the Hamiltonian of the open system from center toward both ends, one could adopt the edge sites with a negligibly small energy scale as the grand canonical small particle bath, and equilibrium states with noninteger arbitrary conserved numbers, e.g., electron numbers or ${s}_{z}$, are realized in the main part of the system. This will enable the evaluation of response functions under a continuously varying external field in a small lattice without any fine-tuning or scaling of parameters while keeping the standard numerical accuracy. Demonstrations are given on quantum spin systems and on a Hubbard model by the density-matrix renormalization group.

arXiv:1307.3713•Physical Review B•Jul 24, 2012•42 citations

Localization and ‘classical entanglement’ in the discrete non-linear Schrödinger equation

Martina Giachello, Stefano Iubini, Roberto Livi +1 more

Abstract We perform a detailed numerical study of the very peculiar thermodynamic properties of the localized high-energy phase of the discrete non-linear Schrödinger equation (DNLSE). A numerical sampling of the microcanonical ensemble done by means of Hamiltonian dynamics reveals a new and subtle relation between the presence of the localized phase and a property of the system that we have called classical entanglement . Our main finding is that a quantity defined for our classical system in perfect analogy with the entanglement entropy of quantum ones, and that we have therefore called S ent , grows with the system size N in the localized phase as S ent ( N ) ∼ log ⁡ ( N ) , therefore revealing the presence of subtle non-local correlations between any finite portion of the system and the rest of it. This manifestation of classical entanglement beautifully captures the lack of system separability in the DNLSE localized phase, revealing how statistical correlations specific to the microcanonical ensemble and non-reproducible in the canonical one, may concur to determine a property totally analogous to the one produced by non-local quantum correlations.

arXiv:2503.14364•Journal of Statistical Mechanics Theory and Experiment•Nov 1, 2025•0 citations

Nonequilibrium discrete nonlinear Schrödinger equation

Stefano Iubini, Stefano Lepri A, Antonio Politi

We study nonequilibrium steady states of the one-dimensional discrete nonlinear Schrödinger equation. This system can be regarded as a minimal model for the stationary transport of bosonic particles such as photons in layered media or cold atoms in deep optical traps. Due to the presence of two conserved quantities, namely, energy and norm (or number of particles), the model displays coupled transport in the sense of linear irreversible thermodynamics. Monte Carlo thermostats are implemented to impose a given temperature and chemical potential at the chain ends. As a result, we find that the Onsager coefficients are finite in the thermodynamic limit, i.e., transport is normal. Depending on the position in the parameter space, the "Seebeck coefficient" may be either positive or negative. For large differences between the thermostat parameters, density and temperature profiles may display an unusual nonmonotonic shape. This is due to the strong dependence of the Onsager coefficients on the state variables.

arXiv:1204.2470•Physical Review E•Jul 6, 2012•62 citations
During the evolution of coupled nonlinear oscillators on a lattice, with dynamics dictated by the discrete nonlinear Schrödinger equation (DNLSE systems), two quantities are conserved: system energy (Hamiltonian) and system density (number of particles). If the number of system oscillators is large enough, a significant portion of the array can be considered to be an "open system", in intimate energy and density contact with a "bath" - the rest of the array. Thus, as indicated in previous works, the grand canonical formulation can be exploited in order to determine equilibrium statistical properties of thermalized DNLSE systems. In this work, given the values of the two conserved quantities, we have calculated the necessary values of the two Lagrange parameters (typically designated β,μβ,μβ,μ) associated with the grand canonical partition function in two different ways. One is numerical and the other is analytic, based on a published approximate entropy expression. In addition we have accessed a purposely-derived approximate PDF expression of site-densities. Applying these mathematical tools we have generated maps of temperatures, chemical potentials, and field correlations for DNLSE systems over the entire thermalization zone of the DNLSE phase diagram, subjected to all system-nonlinearity levels. The end result is a rather complete picture, characterizing equilibrated large DNLSE systems.
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.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.
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 exploit a prescription to observe directly the physical properties of the thermodynamic limit in a continuously applied field in one-dimensional quantum finite lattice systems. By systematically scaling down the energy of the Hamiltonian of the open system from center toward both ends, one could adopt the edge sites with a negligibly small energy scale as the grand canonical small particle bath, and equilibrium states with noninteger arbitrary conserved numbers, e.g., electron numbers or sz{s}_{z}sz​, are realized in the main part of the system. This will enable the evaluation of response functions under a continuously varying external field in a small lattice without any fine-tuning or scaling of parameters while keeping the standard numerical accuracy. Demonstrations are given on quantum spin systems and on a Hubbard model by the density-matrix renormalization group.
Abstract We perform a detailed numerical study of the very peculiar thermodynamic properties of the localized high-energy phase of the discrete non-linear Schrödinger equation (DNLSE). A numerical sampling of the microcanonical ensemble done by means of Hamiltonian dynamics reveals a new and subtle relation between the presence of the localized phase and a property of the system that we have called classical entanglement . Our main finding is that a quantity defined for our classical system in perfect analogy with the entanglement entropy of quantum ones, and that we have therefore called S ent , grows with the system size N in the localized phase as S ent ( N ) ∼ log ⁡ ( N ) , therefore revealing the presence of subtle non-local correlations between any finite portion of the system and the rest of it. This manifestation of classical entanglement beautifully captures the lack of system separability in the DNLSE localized phase, revealing how statistical correlations specific to the microcanonical ensemble and non-reproducible in the canonical one, may concur to determine a property totally analogous to the one produced by non-local quantum correlations.
We study nonequilibrium steady states of the one-dimensional discrete nonlinear Schrödinger equation. This system can be regarded as a minimal model for the stationary transport of bosonic particles such as photons in layered media or cold atoms in deep optical traps. Due to the presence of two conserved quantities, namely, energy and norm (or number of particles), the model displays coupled transport in the sense of linear irreversible thermodynamics. Monte Carlo thermostats are implemented to impose a given temperature and chemical potential at the chain ends. As a result, we find that the Onsager coefficients are finite in the thermodynamic limit, i.e., transport is normal. Depending on the position in the parameter space, the "Seebeck coefficient" may be either positive or negative. For large differences between the thermostat parameters, density and temperature profiles may display an unusual nonmonotonic shape. This is due to the strong dependence of the Onsager coefficients on the state variables.