Source author record

Lapo Casetti

Lapo Casetti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2023arXiv

Modified Thirring model beyond the excluded-volume approximation

Long-range interacting systems may exhibit ensemble inequivalence and can possibly attain equilibrium states under completely open conditions, for which energy, volume and number of particles simultaneously fluctuate. Here we consider a modified version of the Thirring model for self-gravitating systems with attractive and repulsive long-range interactions in which particles are treated as hard spheres in dimension d=1,2,3. Equilibrium states of the model are studied under completely open conditions, in the unconstrained ensemble, by means of both Monte Carlo simulations and analytical methods and are compared with the corresponding states at fixed number of particles, in the isothermal-isobaric ensemble. Our theoretical description is performed for an arbitrary local equation of state, which allows us to examine the system beyond the excluded-volume approximation. The simulations confirm the theoretical prediction of the possible occurrence of first-order phase transitions in the unconstrained ensemble. This work contributes to the understanding of long-range interacting systems exchanging heat, work and matter with the environment.

preprint2022arXiv

Symplectic coarse graining approach to the dynamics of spherical self-gravitating systems

We investigate the evolution of the phase-space distribution function around slightly perturbed stationary states and the process of violent relaxation in the context of the dissipationless collapse of an isolated spherical self-gravitating system. By means of the recently introduced symplectic coarse graining technique, we obtain an effective evolution equation that allows us to compute the scaling of the frequencies around a stationary state, as well as the damping times of Fourier modes of the distribution function, with the magnitude of the Fourier $k-$vectors themselves. We compare our analytical results with $N$-body simulations.

preprint2021arXiv

Violent relaxation in the Hamiltonian mean field model: II. Non-equilibrium phase diagrams

A classical long-range-interacting $N$-particle system relaxes to thermal equilibrium on time scales growing with $N$; in the limit $N\to \infty$ such a relaxation time diverges. However, a completely non-collisional relaxation process, known as violent relaxation, takes place on a much shorter time scale independent of $N$ and brings the system towards a non-thermal quasi-stationary state. A finite system will eventually reach thermal equilibrium, while an infinite system will remain trapped in the quasi-stationary state forever. For times smaller than the relaxation time the distribution function of the system obeys the collisionless Boltzmann equation, also known as the Vlasov equation. The Vlasov dynamics is invariant under time reversal so that it does not "naturally" describe a relaxational dynamics. However, as time grows the dynamics affects smaller and smaller scales in phase space, so that observables not depending upon small-scale details appear as relaxed after a short time. Herewith we present an approximation scheme able to describe violent relaxation in a one-dimensional toy-model, the Hamiltonian Mean Field (HMF). The approach described here generalizes the one proposed in G. Giachetti and L. Casetti, J. Stat. Mech.: Theory Exp. 2019, 043201 (2019), that was limited to "cold" initial conditions, to generic initial conditions, allowing us to to predict non-equilibrium phase diagrams that turn out to be in good agreement with those obtained from the numerical integration of the Vlasov equation.

preprint2020arXiv

Coarse-grained collisionless dynamics with long-range interactions

We present an effective evolution equation for a coarse-grained distribution function of a long-range-interacting system preserving the symplectic structure of the non-collisional Boltzmann, or Vlasov, equation. We first derive a general form of such an equation based on symmetry considerations only. Then, we explicitly derive the equation for one-dimensional systems, finding that it has the form predicted on general grounds. Finally, we use such an equation to predict the dependence of the damping times on the coarse-graining scale and numerically check it for some one-dimensional models, including the Hamiltonian Mean Field (HMF) model, a scalar field with quartic interaction, a 1-d self-gravitating system, and the Self-Gravitating Ring (SGR).

preprint2020arXiv

Discreteness effects, $N-$body chaos and the onset of radial-orbit instability

We study the stability of a family of spherical equilibrium models of self-gravitating systems, the so-called $γ-$models with Osipkov-Merritt velocity anisotropy, by means of $N-$body simulations. In particular, we analyze the effect of self-consistent $N-$body chaos on the onset of radial-orbit instability (ROI). We find that degree of chaoticity of the system associated to its largest Lyapunov exponent $Λ_{\rm max}$ has no appreciable relation with the stability of the model for fixed density profile and different values of radial velocity anisotropy. However, by studying the distribution of the Lyapunov exponents $λ_{\rm m}$ of the individual particles in the single-particle phase space, we find that more anisotropic systems have a larger fraction of orbits with larger $λ_{\rm m}$.

preprint2019arXiv

$N$-body chaos, phase-space transport and relaxation in numerical simulations

Using direct $N$-body simulations of self-gravitating systems we study the dependence of dynamical chaos on the system size $N$. We find that the $N$-body chaos quantified in terms of the largest Lyapunov exponent $Λ_{\rm max}$ decreases with $N$. The values of its inverse (the so-called Lyapunov time $t_λ$) are found to be smaller than the two-body collisional relaxation time but larger than the typical violent relaxation time, thus suggesting the existence of another collective time scale connected to many-body chaos.

preprint2016arXiv

Phase transitions in Thirring's model

In his pioneering work on negative specific heat, Walter Thirring in\-tro\-duced a model that is solvable in the microcanonical ensemble. Here, we give a complete description of the phase-diagram of this model in both the microcanonical and the canonical ensemble, highlighting the main features of ensemble inequivalence. In both ensembles, we find a line of first-order phase transitions which ends in a critical point. However, neither the line nor the point have the same location in the phase-diagram of the two ensembles. We also show that the microcanonical and canonical critical points can be analytically related to each other using a Landau expansion of entropy and free energy, respectively, in analogy with what has been done in [O. Cohen, D. Mukamel, J. Stat. Mech., P12017 (2012)]. Examples of systems with certain symmetries restricting the Landau expansion have been considered in this reference, while no such restrictions are present in Thirring's model. This leads to a phase diagram that can be seen as a prototype for what happens in systems of particles with kinematic degrees of freedom dominated by long-range interactions.

preprint2016arXiv

Surprises from quenches in long-range interacting systems: Temperature inversion and cooling

What happens when one of the parameters governing the dynamics of a long-range interacting system of particles in thermal equilibrium is abruptly changed (quenched) to a different value? While a short-range system, under the same conditions, will relax in time to a new thermal equilibrium with a uniform temperature across the system, a long-range system shows a fast relaxation to a nonequilibrium {\em quasistationary state} (QSS). The lifetime of such an off-equilibrium state diverges with the system size, and the temperature is non-uniform across the system. Quite surprisingly, the density profile in the QSS obtained after the quench is anticorrelated with the temperature profile in space, thus exhibiting the phenomenon of {\em temperature inversion}: denser regions are colder than sparser ones. We illustrate with extensive molecular dynamics simulations the ubiquity of this scenario in a prototypical long-range interacting system subject to a variety of quenching protocols, and in a model that mimics an experimental setup of atoms interacting with light in an optical cavity. We further demonstrate how a procedure of iterative quenching combined with filtering out the high-energy particles in the system may be employed to cool the system. Temperature inversion is observed in nature in some astrophysical settings; our results imply that such a phenomenon should be observable, and could even be exploitable to advantage, also in controlled laboratory experiments.

preprint2015arXiv

Temperature inversion in long-range interacting systems

Temperature inversions occur in nature, e.g., in the solar corona and in interstellar molecular clouds: somewhat counterintuitively, denser parts of the system are colder than dilute ones. We propose a simple and appealing way to spontaneously generate temperature inversions in systems with long-range interactions, by preparing them in inhomogeneous thermal equilibrium states and then applying an impulsive perturbation. In similar situations, short-range systems would typically relax to another thermal equilibrium, with uniform temperature profile. By contrast, in long-range systems, the interplay between wave-particle interaction and spatial inhomogeneity drives the system to nonequilibrium stationary states that generically exhibit temperature inversion. We demonstrate this mechanism in a simple mean-field model and in a two-dimensional self-gravitating system. Our work underlines the crucial role the range of interparticle interaction plays in determining the nature of steady states out of thermal equilibrium.

preprint2015arXiv

Thermodynamics of nonadditive systems

The usual formulation of thermodynamics is based on the additivity of macroscopic systems. However, there are numerous examples of macroscopic systems that are not additive, due to the long-range character of the interaction among the constituents. We present here an approach in which nonadditive systems can be described within a purely thermodynamics formalism. The basic concept is to consider a large ensemble of replicas of the system where the standard formulation of thermodynamics can be naturally applied and the properties of a single system can be consequently inferred. After presenting the approach, we show its implementation in systems where the interaction decays as $1/r^α$ in the interparticle distance $r$, with $α$ smaller than the embedding dimension $d$, and in the Thirring model for gravitational systems.

preprint2014arXiv

Critical energy density of O$(n)$ models in $d=3$

A relation between O$(n)$ models and Ising models has been recently conjectured [L. Casetti, C. Nardini, and R. Nerattini, Phys. Rev. Lett. 106, 057208 (2011)]. Such a relation, inspired by an energy landscape analysis, implies that the microcanonical density of states of an O$(n)$ spin model on a lattice can be effectively approximated in terms of the density of states of an Ising model defined on the same lattice and with the same interactions. Were this relation exact, it would imply that the critical energy densities of all the O$(n)$ models (i.e., the average values per spin of the O$(n)$ Hamiltonians at their respective critical temperatures) should be equal to that of the corresponding Ising model; it is therefore worth investigating how different the critical energies are and how this difference depends on $n$. We compare the critical energy densities of O$(n)$ models in three dimensions in some specific cases: the O$(1)$ or Ising model, the O$(2)$ or $XY$ model, the O$(3)$ or Heisenberg model, the O$(4)$ model and the O$(\infty)$ or spherical model, all defined on regular cubic lattices and with ferromagnetic nearest-neighbor interactions. The values of the critical energy density in the $n=2$, $n=3$, and $n=4$ cases are derived through a finite-size scaling analysis of data produced by means of Monte Carlo simulations on lattices with up to $128^3$ sites. For $n=2$ and $n=3$ the accuracy of previously known results has been improved. We also derive an interpolation formula showing that the difference between the critical energy densities of O$(n)$ models and that of the Ising model is smaller than $1\%$ if $n<8$ and never exceeds $3\%$ for any $n$.

preprint2014arXiv

Velocity filtration and temperature inversion in a system with long-range interactions

Temperature inversion due to velocity filtration, a mechanism originally proposed to explain the heating of the solar corona, is demonstrated to occur also in a simple paradigmatic model with long-range interactions, the Hamiltonian mean-field model. Using molecular dynamics simulations, we show that when the system settles into an inhomogeneous quasi-stationary state in which the velocity distribution has suprathermal tails, the temperature and density profiles are anticorrelated: denser parts of the system are colder than dilute ones. We argue that this may be a generic property of long-range interacting systems.

preprint2013arXiv

Density of states of the $XY$ model: an energy landscape approach

Among the stationary configurations of the Hamiltonian of a classical O$(n)$ lattice spin model, a class can be identified which is in one-to-one correspondence with all the the configurations of an Ising model defined on the same lattice and with the same interactions. Starting from this observation it has been recently proposed that the microcanonical density of states of an O$(n)$ model could be written in terms of the density of states of the corresponding Ising model. Later, it has been shown that a relation of this kind holds exactly for two solvable models, the mean-field and the one-dimensional $XY$ model, respectively. We apply the same strategy to derive explicit, albeit approximate, expressions for the density of states of the two-dimensional $XY$ model with nearest-neighbor interactions on a square lattice. The caloric curve and the specific heat as a function of the energy density are calculated and compared against simulation data, yielding a very good agreement over the entire energy density range. The concepts and methods involved in the approximations presented here are valid in principle for any O$(n)$ model.

preprint2012arXiv

Caloric curve of star clusters

Self-gravitating systems, like globular clusters or elliptical galaxies, are the prototypes of many-body systems with long-range interactions, and should be the natural arena where to test theoretical predictions on the statistical behaviour of long-range-interacting systems. Systems of classical self-gravitating particles can be studied with the standard tools of equilibrium statistical mechanics, provided the potential is regularized at small length scales and the system is confined in a box. The confinement condition looks rather unphysical in general, so that it is natural to ask whether what we learn with these studies is relevant to real self-gravitating systems. In order to provide a first answer to this question we consider a basic, simple, yet effective model of globular clusters, the King model. This model describes a self-consistently confined system, without the need of any external box, but the stationary state is a non-thermal one. In particular, we consider the King model with a short-distance cutoff on the interactions and we discuss how such a cutoff affects the caloric curve, i.e. the relation between temperature and energy. We find that the cutoff stabilizes a low-energy phase which is absent in the King model without cutoff; the caloric curve of the model with cutoff turns out to be very similar to that of previously studied confined and regularized models, but for the absence of a high-energy gas-like phase. We briefly discuss the possible phenomenological as well as theoretical implications of these results.

preprint2012arXiv

Density of states of continuous and discrete spin models: a case study

A relation between O(n) lattice spin models and Ising models defined on the same lattice was recently put forward [L. Casetti, C. Nardini, and R. Nerattini, Phys. Rev. Lett. 106, 057208 (2011)]. Such a relation, inspired by an energy landscape analysis, implies that the density of states of an O(n) spin model on a lattice can be effectively approximated, at least close to the phase transition, in terms of the density of states of an Ising model defined on the same lattice and with the same interactions. In the present paper we show that such a relation exactly holds, albeit in a slightly modified form, in the special cases of the mean-field XY model and of the one-dimensional XY model. We also discuss the possible consequences of this result for the general case.

preprint2012arXiv

Exploring the energy landscape of XY models

We investigate the energy landscape of two- and three-dimensional XY models with nearest-neighbor interactions by analytically constructing several classes of stationary points of the Hamiltonian. These classes are analyzed, in particular with respect to possible signatures of the thermodynamic phase transitions of the models. We find that, even after explicitly breaking the global O(2) symmetry of the XY spins, an exponentially large class of stationary points are singular and occur in continuous one-parameter families. This property may complicate the use of theoretical tools developed for the investigation of phase transitions based on stationary points of the energy landscape, and we discuss strategies to avoid these difficulties.

preprint2012arXiv

Geometry of the energy landscape of the self-gravitating ring

We study the global geometry of the energy landscape of a simple model of a self-gravitating system, the self-gravitating ring (SGR). This is done by endowing the configuration space with a metric such that the dynamical trajectories are identified with geodesics. The average curvature and curvature fluctuations of the energy landscape are computed by means of Monte Carlo simulations and, when possible, of a mean-field method, showing that these global geometric quantities provide a clear geometric characterization of the collapse phase transition occurring in the SGR as the transition from a flat landscape at high energies to a landscape with mainly positive but fluctuating curvature in the collapsed phase. Moreover, curvature fluctuations show a maximum in correspondence with the energy of a possible further transition, occurring at lower energies than the collapse one, whose existence had been previously conjectured on the basis of a local analysis of the energy landscape and whose effect on the usual thermodynamic quantities, if any, is extremely weak. We also estimate the largest Lyapunov exponent $λ$ of the SGR using the geometric observables. The geometric estimate always gives the correct order of magnitude of $λ$ and is also quantitatively correct at small energy densities and, in the limit $N\to\infty$, in the whole homogeneous phase.

preprint2010arXiv

On a microcanonical relation between continuous and discrete spin models

A relation between a class of stationary points of the energy landscape of continuous spin models on a lattice and the configurations of a Ising model defined on the same lattice suggests an approximate expression for the microcanonical density of states. Based on this approximation we conjecture that if a O(n) model with ferromagnetic interactions on a lattice has a phase transition, its critical energy density is equal to that of the n = 1 case, i.e., a system of Ising spins with the same interactions. The conjecture holds true in the case of long-range interactions. For nearest-neighbor interactions, numerical results are consistent with the conjecture for n=2 and n=3 in three dimensions. For n=2 in two dimensions (XY model) the conjecture yields a prediction for the critical energy of the Berezinskij-Kosterlitz-Thouless transition, which would be equal to that of the two-dimensional Ising model. We discuss available numerical data in this respect.

preprint2010arXiv

Solvable model of a self-gravitating system

We introduce and discuss an effective model of a self-gravitating system whose equilibrium thermodynamics can be solved in both the microcanonical and the canonical ensemble, up to a maximization with respect to a single variable. Such a model can be derived from a model of self-gravitating particles confined on a ring, referred to as the self-gravitating ring (SGR) model, allowing a quantitative comparison between the thermodynamics of the two models. Despite the rather crude approximations involved in its derivation, the effective model compares quite well with the SGR model. Moreover, we discuss the relation between the effective model presented here and another model introduced by Thirring forty years ago. The two models are very similar and can be considered as examples of a class of minimal models of self-gravitating systems.

preprint2009arXiv

Energy landscape and phase transitions in the self-gravitating ring model

We apply a recently proposed criterion for the existence of phase transitions, which is based on the properties of the saddles of the energy landscape, to a simplified model of a system with gravitational interactions, referred to as the self-gravitating ring model. We show analytically that the criterion correctly singles out the phase transition between a homogeneous and a clustered phase and also suggests the presence of another phase transition, not previously known. On the basis of the properties of the energy landscape we conjecture on the nature of the latter transition.

preprint2004arXiv

Weak and strong chaos in FPU models and beyond

We briefly review some of the most relevant results that our group obtained in the past, while investigating the dynamics of the Fermi-Pasta-Ulam (FPU) models. A first result is the numerical evidence of the existence of two different kinds of transitions in the dynamics of the FPU models: i) a Stochasticity Threshold (ST), characterized by a value of the energy per degree of freedom below which the overwhelming majority of the phase space trajectories are regular (vanishing Lyapunov exponents). It tends to vanish as the number N of degrees of freedom is increased. ii) a Strong Stochasticity Threshold (SST), characterized by a value of the energy per degree of freedom at which a crossover appears between two different power laws of the energy dependence of the largest Lyapunov exponent, which phenomenologically corresponds to the transition between weakly and strongly chaotic regimes. It is stable with N. A second result is the development of a Riemannian geometric theory to explain the origin of Hamiltonian chaos. The starting of this theory has been motivated by the inadequacy of the approach based on homoclinic intersections to explain the origin of chaos in systems of arbitrarily large N, or arbitrarily far from quasi-integrability, or displaying a transition between weak and strong chaos. Finally, a third result stems from the search for the transition between weak and strong chaos in systems other than FPU. Actually, we found that a very sharp SST appears as the dynamical counterpart of a thermodynamic phase transition, which in turn has led, in the light of the Riemannian theory of chaos, to the development of a topological theory of phase transitions.