Temperature
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 microcanonica
l 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$.
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 th
e 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.
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.
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 w
ith 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 $lambda$ of the SGR using the geometric observables. The geometric estimate always gives the correct order of magnitude of $lambda$ and is also quantitatively correct at small energy densities and, in the limit $Ntoinfty$, in the whole homogeneous phase.
