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

Quenched dynamics of classical isolated systems: The spherical spin model with two-body random interactions or the Neumann integrable model

Leticia F. Cugliandolo, Gustavo Sergio Lozano, Nessi, Emilio Nicolás, Picco, Marcos Fernando, Alessandro Tartaglia

Published Jun 1, 2018• 13 citations
OpenAlex

Abstract

We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N - 1 integrals of motion, notably, their scaling with N, and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs-Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.

Related papers

Quenched dynamics of classical isolated systems: the spherical spin model with two-body random interactions or the Neumann integrable model

Leticia F. Cugliandolo, Gustavo Sergio Lozano, Nicolás Nessi +2 more

Abstract We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N − 1 integrals of motion, notably, their scaling with N , and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs–Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.

arXiv:1712.07688•Journal of Statistical Mechanics Theory and Experiment•Jun 1, 2018•13 citations

(Non-equilibrium) thermodynamics of integrable models: The Generalized Gibbs Ensemble description of the classical Neumann model

Damien Barbier, Leticia F. Cugliandolo, Gustavo Sergio Lozano +1 more

Abstract We study the motion of a classical particle subject to anisotropic harmonic forces and constrained to move on the sphere. In the integrable-systems literature this problem is known as the Neumann model. We choose the spring constants in a way that makes the connection with the so-called p = 2 spherical disordered system transparent. We tackle the problem in the N → ∞ limit by introducing a soft version in which the spherical constraint is imposed only on average over initial conditions. We show that the Generalized Gibbs Ensemble, constructed with N conserved charges in involution, captures the long-time averages of all relevant observables of the soft model after sudden changes in the parameters (quenches). We reveal the full dynamic phase diagram with four different phases in which the particles' position and momentum are both extended, only the position quasi-condenses or condenses, and both condense. The scaling properties of the fluctuations allow us to establish in which of these cases the strict and soft spherical constraints are equivalent. We thus confirm the validity of the GGE hypothesis for the Neumann model on a large portion of the dynamic phase diagram.

arXiv:2007.14364•Europhysics Letters (EPL)•Dec 1, 2020•7 citations

(Non equilibrium) Thermodynamics of Integrable models: The Generalized Gibbs Ensemble description of the classical Neumann Model

Damien Barbier, Leticia F. Cugliandolo, Gustavo Sergio Lozano +1 more

We study the motion of a classical particle subject to anisotropic harmonic forces and constrained to move on the SN-1 sphere. In the integrable-systems literature this problem is known as the Neumann model. We choose the spring constants in a way that makes the connection with the so-called p = 2 spherical disordered system transparent. We tackle the problem in the N → ∞ limit by introducing a soft version in which the spherical constraint is imposed only on average over initial conditions. We show that the Generalized Gibbs Ensemble, constructed with N conserved charges in involution, captures the long-time averages of all relevant observables of the soft model after sudden changes in the parameters (quenches). We reveal the full dynamic phase diagram with four different phases in which the particles' position and momentum are both extended, only the position quasi-condenses or condenses, and both condense. The scaling properties of the fluctuations allow us to establish in which of these cases the strict and soft spherical constraints are equivalent. We thus confirm the validity of the GGE hypothesis for the Neumann model on a large portion of the dynamic phase diagram.

•Oct 1, 2020•10 citations

Non equilibrium dynamics of isolated disordered systems: the classical Hamiltonian p-spin model

Leticia F. Cugliandolo, Gustavo Sergio Lozano, Nicolás Nessi

Abstract We study the dynamics of a classical disordered macroscopic model completely isolated from the environment reproducing, in a classical setting, the ‘quantum quench’ protocol. We show that, depending on the pre and post quench parameters, the system approaches equilibrium, succeeding to act as a bath on itself, or remains out of equilibrium, in two different ways. In the latter one, the system stays confined in a metastable state in which it undergoes stationary dynamics characterised by a single temperature. In the other, the system ages and its dynamics are characterised by two temperatures associated with observations made at short and long time differences (high and low frequencies). The parameter dependence of the asymptotic states is rationalised in terms of a dynamic phase diagram with one equilibrium and two out of equilibrium phases. Aspects of pre-thermalisation are observed and discussed. Similarities and differences with the dynamics of the dissipative model are also explained.

arXiv:1703.10080•Journal of Statistical Mechanics Theory and Experiment•Aug 1, 2017•18 citations

Langevin dynamics towards Generalized Gibbs Ensembles

Rafael González-Albaladejo, Leticia F. Cugliandolo, Gustavo Sergio Lozano +1 more

In recent years, statistical physics tools have been used to study the dynamics of macroscopic quantum and classical integrable models. These studies have shown that, in the thermodynamic limit, long term time-averages of local observables equal statistical averages calculated with a Generalized Gibbs Ensemble (GGE). In the latter, the weight in phase space involves all constants of motion multiplied by Lagrange multipliers which are fixed by their initial values. We argue that the GGE can be interpreted as a proper Boltzmann measure of a different (non-standard) mechanical model. This measure can therefore be sampled with stochastic methods implemented, for example, as a Langevin process. We tailor the deterministic force in the Langevin equations so that at long times the stochastic trajectories sample the desired measure. We use the integrable Neumann Model -- the model of a classical particle constrained to move on a sphere under anisotropic harmonic potentials -- as a benchmark. We test the method in a setup where the Lagrange multipliers are fixed by the GGE governing the asymptotic dynamics of the isolated Neumann model following a quench. This problem has different phases depending on the initial conditions and quench protocol and, in all cases, with the stochastic evolution of an adequately tuned Langevin process we reproduce the asymptotic values of relevant observables.

arXiv:2610.37129•Oct 27, 2026•0 citations

Quantum quench dynamics in the transverse field Ising model of fully connected spins

Ingo Homrighausen

The transverse field Ising model of N fully connected spins provides a mathematical accessible, yet non-trivial, system to investigate the non-thermal dynamics after a sudden quantum quench. The focus in this thesis is three-fold. First, in a 1/N expansion the dynamics of the mean magnetization within its variance is obtained analytically in the large N limit. The variance constitutes a leading order correction to the mean field limit at infinite N. While mean field theories are one of the most accessible approximations to deal with the quantum complexity, its validity in time out of equilibrium has to be checked on a case by case basis. By studying the dynamics of the variance, we see that the mean field approximation can break down at surprisingly early times scaling with the square root of N. The underlying mechanism in the fully connected Ising model is identified as a dephasing effect between nearby effective orbits. The dynamics of the magnetization variance discriminates four qualitatively different regions in the dynamical phase diagram that are indistinguishable by just looking at the mean magnetization. Second, bipartite entanglement between disjoint groups of the spins in the fully connected Ising model is studied. In the large N limit the entanglement Hamiltonian is shown to be a harmonic oscillator. A quantitative relation between the angular frequency of the oscillator and a squeezing measure of the spin state is derived. Explicit expressions for the Renyi entanglement entropies are given. In contrast to the variance, which vanishes in the thermodynamic limit, the entanglement entropy saturate to an N-independent positive value. The dynamics of the entanglement Hamiltonian and the entanglement entropies is discussed. Linear growth as a function of time is linked to the unstable fixed point and the homoclinic orbit of a pitchfork bifurcation, while logarithmic increase is explained due to a dephasing mechanism of nearby effective orbits. A generalization to other models that permit a semiclassical effective description (on a possibly higher dimensional phase space) is given. As a byproduct the connection between entangled Gaussian states and the symplectic eigenvalues of its reduced covariance is derived. An upper bound on the von Neumann entanglement entropy of any bipartite pure state is given in terms of the symplectic properties of its reduced covariance. Third, the existence of dynamical phase transitions (DPTs) in the fully connected Ising model is argued by two different approaches. DPTs are defined as non-analyticities in the time-dependence of the return probability rate function in the thermodynamic limit. Like entanglement entropy, the behavior of the return rate function is an observable-independent quantity used to characterize non-equilibrium states. A numerically based connection between the DPT and the analytic property of the infinite time averaged order parameter as a function of quench strength is discussed. The concept of Feynman history states is used to devise the notion of history entanglement. An information theoretic interpretation of the history entanglement entropy as the precision of a quantum clock follows from Holevo’s bound. The dependence of this entanglement entropy on the observation time constitutes an operator-independent characterization of the quantum dynamics in closed systems. For observation times comparable to the Heisenberg time, the history entanglement becomes sensitive to the statistics of the energy gaps, and thus distinguishes integrable and non-integrable dynamics. This is demonstrated for two examples, a random matrix ensemble, and a one-dimensional hardcore boson model with integrability breaking next nearest neighbor interaction.

•Jan 1, 2021•0 citations
We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N - 1 integrals of motion, notably, their scaling with N, and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs-Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.
Abstract We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N − 1 integrals of motion, notably, their scaling with N , and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs–Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.
Abstract We study the motion of a classical particle subject to anisotropic harmonic forces and constrained to move on the sphere. In the integrable-systems literature this problem is known as the Neumann model. We choose the spring constants in a way that makes the connection with the so-called p = 2 spherical disordered system transparent. We tackle the problem in the N → ∞ limit by introducing a soft version in which the spherical constraint is imposed only on average over initial conditions. We show that the Generalized Gibbs Ensemble, constructed with N conserved charges in involution, captures the long-time averages of all relevant observables of the soft model after sudden changes in the parameters (quenches). We reveal the full dynamic phase diagram with four different phases in which the particles' position and momentum are both extended, only the position quasi-condenses or condenses, and both condense. The scaling properties of the fluctuations allow us to establish in which of these cases the strict and soft spherical constraints are equivalent. We thus confirm the validity of the GGE hypothesis for the Neumann model on a large portion of the dynamic phase diagram.
We study the motion of a classical particle subject to anisotropic harmonic forces and constrained to move on the SN-1 sphere. In the integrable-systems literature this problem is known as the Neumann model. We choose the spring constants in a way that makes the connection with the so-called p = 2 spherical disordered system transparent. We tackle the problem in the N → ∞ limit by introducing a soft version in which the spherical constraint is imposed only on average over initial conditions. We show that the Generalized Gibbs Ensemble, constructed with N conserved charges in involution, captures the long-time averages of all relevant observables of the soft model after sudden changes in the parameters (quenches). We reveal the full dynamic phase diagram with four different phases in which the particles' position and momentum are both extended, only the position quasi-condenses or condenses, and both condense. The scaling properties of the fluctuations allow us to establish in which of these cases the strict and soft spherical constraints are equivalent. We thus confirm the validity of the GGE hypothesis for the Neumann model on a large portion of the dynamic phase diagram.
Abstract We study the dynamics of a classical disordered macroscopic model completely isolated from the environment reproducing, in a classical setting, the ‘quantum quench’ protocol. We show that, depending on the pre and post quench parameters, the system approaches equilibrium, succeeding to act as a bath on itself, or remains out of equilibrium, in two different ways. In the latter one, the system stays confined in a metastable state in which it undergoes stationary dynamics characterised by a single temperature. In the other, the system ages and its dynamics are characterised by two temperatures associated with observations made at short and long time differences (high and low frequencies). The parameter dependence of the asymptotic states is rationalised in terms of a dynamic phase diagram with one equilibrium and two out of equilibrium phases. Aspects of pre-thermalisation are observed and discussed. Similarities and differences with the dynamics of the dissipative model are also explained.
In recent years, statistical physics tools have been used to study the dynamics of macroscopic quantum and classical integrable models. These studies have shown that, in the thermodynamic limit, long term time-averages of local observables equal statistical averages calculated with a Generalized Gibbs Ensemble (GGE). In the latter, the weight in phase space involves all constants of motion multiplied by Lagrange multipliers which are fixed by their initial values. We argue that the GGE can be interpreted as a proper Boltzmann measure of a different (non-standard) mechanical model. This measure can therefore be sampled with stochastic methods implemented, for example, as a Langevin process. We tailor the deterministic force in the Langevin equations so that at long times the stochastic trajectories sample the desired measure. We use the integrable Neumann Model -- the model of a classical particle constrained to move on a sphere under anisotropic harmonic potentials -- as a benchmark. We test the method in a setup where the Lagrange multipliers are fixed by the GGE governing the asymptotic dynamics of the isolated Neumann model following a quench. This problem has different phases depending on the initial conditions and quench protocol and, in all cases, with the stochastic evolution of an adequately tuned Langevin process we reproduce the asymptotic values of relevant observables.
The transverse field Ising model of N fully connected spins provides a mathematical accessible, yet non-trivial, system to investigate the non-thermal dynamics after a sudden quantum quench. The focus in this thesis is three-fold. First, in a 1/N expansion the dynamics of the mean magnetization within its variance is obtained analytically in the large N limit. The variance constitutes a leading order correction to the mean field limit at infinite N. While mean field theories are one of the most accessible approximations to deal with the quantum complexity, its validity in time out of equilibrium has to be checked on a case by case basis. By studying the dynamics of the variance, we see that the mean field approximation can break down at surprisingly early times scaling with the square root of N. The underlying mechanism in the fully connected Ising model is identified as a dephasing effect between nearby effective orbits. The dynamics of the magnetization variance discriminates four qualitatively different regions in the dynamical phase diagram that are indistinguishable by just looking at the mean magnetization. Second, bipartite entanglement between disjoint groups of the spins in the fully connected Ising model is studied. In the large N limit the entanglement Hamiltonian is shown to be a harmonic oscillator. A quantitative relation between the angular frequency of the oscillator and a squeezing measure of the spin state is derived. Explicit expressions for the Renyi entanglement entropies are given. In contrast to the variance, which vanishes in the thermodynamic limit, the entanglement entropy saturate to an N-independent positive value. The dynamics of the entanglement Hamiltonian and the entanglement entropies is discussed. Linear growth as a function of time is linked to the unstable fixed point and the homoclinic orbit of a pitchfork bifurcation, while logarithmic increase is explained due to a dephasing mechanism of nearby effective orbits. A generalization to other models that permit a semiclassical effective description (on a possibly higher dimensional phase space) is given. As a byproduct the connection between entangled Gaussian states and the symplectic eigenvalues of its reduced covariance is derived. An upper bound on the von Neumann entanglement entropy of any bipartite pure state is given in terms of the symplectic properties of its reduced covariance. Third, the existence of dynamical phase transitions (DPTs) in the fully connected Ising model is argued by two different approaches. DPTs are defined as non-analyticities in the time-dependence of the return probability rate function in the thermodynamic limit. Like entanglement entropy, the behavior of the return rate function is an observable-independent quantity used to characterize non-equilibrium states. A numerically based connection between the DPT and the analytic property of the infinite time averaged order parameter as a function of quench strength is discussed. The concept of Feynman history states is used to devise the notion of history entanglement. An information theoretic interpretation of the history entanglement entropy as the precision of a quantum clock follows from Holevo’s bound. The dependence of this entanglement entropy on the observation time constitutes an operator-independent characterization of the quantum dynamics in closed systems. For observation times comparable to the Heisenberg time, the history entanglement becomes sensitive to the statistics of the energy gaps, and thus distinguishes integrable and non-integrable dynamics. This is demonstrated for two examples, a random matrix ensemble, and a one-dimensional hardcore boson model with integrability breaking next nearest neighbor interaction.