-
Inference of compressed Potts graphical models
Authors:
Francesca Rizzato,
Alice Coucke,
Eleonora de Leonardis,
J. P. Barton,
Jérôme Tubiana,
Remi Monasson,
Simona Cocco
Abstract:
We consider the problem of inferring a graphical Potts model on a population of variables, with a non-uniform number of Potts colors (symbols) across variables. This inverse Potts problem generally involves the inference of a large number of parameters, often larger than the number of available data, and, hence, requires the introduction of regularization. We study here a double regularization sch…
▽ More
We consider the problem of inferring a graphical Potts model on a population of variables, with a non-uniform number of Potts colors (symbols) across variables. This inverse Potts problem generally involves the inference of a large number of parameters, often larger than the number of available data, and, hence, requires the introduction of regularization. We study here a double regularization scheme, in which the number of colors available to each variable is reduced, and interaction networks are made sparse. To achieve this color compression scheme, only Potts states with large empirical frequency (exceeding some threshold) are explicitly modeled on each site, while the others are grouped into a single state. We benchmark the performances of this mixed regularization approach, with two inference algorithms, the Adaptive Cluster Expansion (ACE) and the PseudoLikelihood Maximization (PLM) on synthetic data obtained by sampling disordered Potts models on an Erdos-Renyi random graphs. We show in particular that color compression does not affect the quality of reconstruction of the parameters corresponding to high-frequency symbols, while drastically reducing the number of the other parameters and thus the computational time. Our procedure is also applied to multi-sequence alignments of protein families, with similar results.
△ Less
Submitted 3 January, 2020; v1 submitted 30 July, 2019;
originally announced July 2019.
-
Nonequilibrium Statistical Mechanics of Systems with Long-Range Interactions: Ubiquity of Core-Halo Distributions
Authors:
Yan Levin,
Renato Pakter,
Felipe B. Rizzato,
Tarcísio N. Teles,
Fernanda P. da C. Benetti
Abstract:
Systems with long-range (LR) forces, for which the interaction potential decays with the interparticle distance with an exponent smaller than the dimensionality of the embedding space, remain an outstanding challenge to statistical physics. The internal energy of such systems lacks extensivity and additivity. Although the extensivity can be restored by scaling the interaction potential with the nu…
▽ More
Systems with long-range (LR) forces, for which the interaction potential decays with the interparticle distance with an exponent smaller than the dimensionality of the embedding space, remain an outstanding challenge to statistical physics. The internal energy of such systems lacks extensivity and additivity. Although the extensivity can be restored by scaling the interaction potential with the number of particles, the non-additivity still remains. Lack of additivity leads to inequivalence of statistical ensembles. Before relaxing to thermodynamic equilibrium, isolated systems with LR forces become trapped in out-of-equilibrium quasi-stationary state (qSS), the lifetime of which diverges with the number of particles. Therefore, in thermodynamic limit LR systems will not relax to equilibrium. The qSSs are attained through the process of collisionless relaxation. Density oscillations lead to particle-wave interactions and excitation of parametric resonances. The resonant particles escape from the main cluster to form a tenuous halo. Simultaneously, this cools down the core of the distribution and dampens out the oscillations. When all the oscillations die out the ergodicity is broken and a qSS is born. In this report, we will review a theory which allows us to quantitatively predict the particle distribution in the qSS. The theory is applied to various LR interacting systems, ranging from plasmas to self-gravitating clusters and kinetic spin models.
△ Less
Submitted 3 October, 2013;
originally announced October 2013.
-
Statistical Mechanics of Unbound Two Dimensional Self-Gravitating Systems
Authors:
Tarcísio N. Teles,
Yan Levin,
Renato Pakter,
Felipe B. Rizzato
Abstract:
We study, using both theory and molecular dynamics simulations, the relaxation dynamics of a microcanonical two dimensional self-gravitating system. After a sufficiently large time, a gravitational cluster of N particles relaxes to the Maxwell-Boltzmann distribution. The time to reach the thermodynamic equilibrium, however, scales with the number of particles. In the thermodynamic limit,…
▽ More
We study, using both theory and molecular dynamics simulations, the relaxation dynamics of a microcanonical two dimensional self-gravitating system. After a sufficiently large time, a gravitational cluster of N particles relaxes to the Maxwell-Boltzmann distribution. The time to reach the thermodynamic equilibrium, however, scales with the number of particles. In the thermodynamic limit, $N\to\infty$ at fixed total mass, equilibrium state is never reached and the system becomes trapped in a non-ergodic stationary state. An analytical theory is presented which allows us to quantitatively described this final stationary state, without any adjustable parameters.
△ Less
Submitted 19 April, 2010; v1 submitted 1 April, 2010;
originally announced April 2010.
-
Thermodynamics of driven collisionless systems
Authors:
Felipe B. Rizzato,
Renato Pakter,
Yan Levin
Abstract:
A statistical theory is presented which allows to calculate the stationary state achieved by a driven system after a process of collisionless relaxation. The theory is applied to study an electron beam driven by an external electric field. The Vlasov equation with appropriate boundary conditions is solved analytically and compared with the molecular dynamics simulation. A perfect agreement is fo…
▽ More
A statistical theory is presented which allows to calculate the stationary state achieved by a driven system after a process of collisionless relaxation. The theory is applied to study an electron beam driven by an external electric field. The Vlasov equation with appropriate boundary conditions is solved analytically and compared with the molecular dynamics simulation. A perfect agreement is found between the theory and the simulations. The full current-voltage phase diagram is constructed.
△ Less
Submitted 1 April, 2009;
originally announced April 2009.
-
Collisionless relaxation in gravitational systems: From violent relaxation to gravothermal collapse
Authors:
Yan Levin,
Renato Pakter,
Felipe B. Rizzato
Abstract:
Theory and simulations are used to study collisionless relaxation of a gravitational $N$-body system. It is shown that when the initial one particle distribution function satisfies the virial condition -- potential energy is minus twice the kinetic energy -- the system quickly relaxes to a metastable state described {\it quantitatively} by the Lynden-Bell distribution with a cutoff. If the initi…
▽ More
Theory and simulations are used to study collisionless relaxation of a gravitational $N$-body system. It is shown that when the initial one particle distribution function satisfies the virial condition -- potential energy is minus twice the kinetic energy -- the system quickly relaxes to a metastable state described {\it quantitatively} by the Lynden-Bell distribution with a cutoff. If the initial distribution function does not meet the virial requirement, the system undergoes violent oscillations, resulting in a partial evaporation of mass. The leftover particles phase separate into a core-halo structure. The theory presented allows us to quantitatively predict the amount and the distribution of mass left in the central core, without any adjustable parameters. On a longer time scale $τ_G \sim N$ collisionless relaxation leads to a gravothermal collapse.
△ Less
Submitted 22 August, 2008;
originally announced August 2008.
-
Superconducting pipes and levitating magnets
Authors:
Yan Levin,
Felipe B. Rizzato
Abstract:
Motivated by a beautiful demonstration of the Faraday's and Lenz's law in which a small neodymium magnet falls slowly through a conducting non-ferromagnetic tube, we consider the dynamics of a magnet falling through a superconducting pipe. Unlike the case of normal conducting pipes, in which the magnet quickly reaches the terminal velocity, inside a superconducting tube the magnet falls freely.…
▽ More
Motivated by a beautiful demonstration of the Faraday's and Lenz's law in which a small neodymium magnet falls slowly through a conducting non-ferromagnetic tube, we consider the dynamics of a magnet falling through a superconducting pipe. Unlike the case of normal conducting pipes, in which the magnet quickly reaches the terminal velocity, inside a superconducting tube the magnet falls freely. On the other hand, to enter the pipe the magnet must overcome a large electromagnetic energy barrier. For sufficiently strong magnets, the barrier is so large that the magnet will not be able to penetrate it and will be suspended over the front edge. We calculate the work that must done to force the magnet to enter a superconducting tube. The calculations show that superconducting pipes are very efficient at screening magnetic fields. For example, the magnetic field of a dipole at the center of a short pipe of radius $a$ and length $L \approx a$ decays, in the axial direction, with a characteristic length $ξ\approx 0.26 a$. The efficient screening of the magnetic field might be useful for shielding highly sensitive superconducting quantum interference devices, SQUIDs. Finally, the motion of the magnet through a superconducting pipe is compared and contrasted to the flow of ions through a trans-membrane channel.
△ Less
Submitted 15 September, 2006;
originally announced September 2006.