Source author record

Sergio Davis

Sergio Davis 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

18works
9topics
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

18 published item(s)

preprint2022arXiv

The $q$-canonical ensemble as a consequence of Bayesian superstatistics

Superstatistics is a generalization of equilibrium statistical mechanics that describes systems in nonequilibrium steady states. Among the possible superstatistical distributions, the $q$-canonical ensemble (also known as Tsallis' statistics, and in plasma physics as Kappa distributions) is probably the most widely used, however the current explanations of its origin are not completely consistent. In this work it is shown that, under a Bayesian interpretation of superstatistics, the origin of the $q$-canonical ensemble can be explained as the superstatistical distribution with maximum Shannon-Jaynes entropy under noninformative constraints. The $q$-canonical distributions are singled out by the mathematical structure of superstatistics itself, and thus no assumptions about the physics of the systems of interest, or regarding their complexity or range of interactions, are needed. These results support the thesis that the success of the $q$-canonical ensemble is information-theoretical in nature, and explainable in terms of the original maximum entropy principle by Jaynes.

preprint2021arXiv

A portable and flexible implementation of the Wang--Landau algorithm in order to determine the Density of States

In this work we develop an implementation of the Wang--Landau algorithm [Phys. Rev. Lett. \textbf{86}, 2050-2053 (2001)]. This algorithm allows us to find the density of states (DOS), a function that, for a given system, describes the proportion of states that have a certain energy. The implementation uses the Python language for the algorithm itself, and it can take advantage of any library, such as the powerful LAMMPS library, for the computation of energy. Therefore, the resulting implementation is simple and flexible without sacrificing efficiency. This implementation also considers recent developments in the parallelization of the code for faster computation. We establish the soundness and effectiveness of our implementation by studying well-known systems such as the Ising model, the Lennard--Jones and EAM solids. We have found that our implementation can find the DOS with very good precision in a reasonable amount of time. Therefore, we are equipped with a very powerful and flexible implementation that can be easily used in order to study more realistic models of matter.

preprint2021arXiv

Fluctuating temperature outside superstatistics: thermodynamics of small systems

The existence of fluctuations of temperature has been a somewhat controversial topic in thermodynamics but nowadays it is recognized that they must be taken into account in small, finite systems. Although for nonequilibrium steady states superstatistics is becoming the \textit{de facto} framework for expressing such temperature fluctuations, some recent results put into question the idea of temperature as a phase space observable. In this work we present and explore the statistics that describes a part of an isolated system, small enough to have well-defined uncertainties in energy and temperature, but lacking a superstatistical description. These results motivate the use of the so-called fundamental temperature as an observable and may be relevant for the statistical description of small systems in physical chemistry.

preprint2020arXiv

Solving equations of motion by using Monte Carlo Metropolis: Novel method via Random Paths and Maximum Caliber Principle

A permanent challenge in physics and other disciplines is to solve partial differential equations, thereby a beneficial investigation is to continue searching for new procedures to do it. In this Letter, a novel Monte-Carlo Metropolis framework is presented for solving the equations of motion in Lagrangian systems. The implementation lies in sampling the paths space with a probability functional obtained by using the maximum caliber principle. The methodology was applied to the free particle and the harmonic oscillator problems, where the numerically-averaged path obtained from the Monte-Carlo simulation converges to the analytical solution from classical mechanics, in an analogous way with a canonical system where energy is minimized by sampling the state space and computing the average state for each system. Thus, we expect that this procedure can be general enough to solve other differential equations in physics and to be a useful tool to calculate the time-dependent properties of dynamical systems in order to understand the non-equilibrium behavior of statistical mechanical systems.

preprint2020arXiv

Superheated solid state induced by a single collision event

High-energy incident particles from both pulsed and continuous radiation sources can induce significant damage to the structure of a material by creating vacancy-interstitial pairs and other more complex defects, and this leads typically to localized melting. In this work, we present evidence via molecular dynamics simulations of a superheated solid state in BCC tungsten induced by single PKA events of $\sim$ 1.5 keV of energy. Despite the striking difference between this highly inhomogeneous energy injection and homogeneous melting, the lifetime of the obtained superheated solid state, reaching up to 200 ps, is described through a waiting time distribution in agreement with previous studies on superheating in the Z-method methodology.

preprint2019arXiv

A general statistical model for waiting times until collapse of a system

The distribution of waiting times until the occurrence of a critical event is a crucial statistical problem across several disciplines in Science. In this work we present a statistical model in which a relevant quantity X accumulates until overcoming a threshold X*, which defines the collapse. The obtained waiting time distribution is a mixture of gamma distributions, which in turn can be approximated as an effective gamma distribution.

preprint2019arXiv

On the possible distributions of temperature in nonequilibrium steady states

Superstatistics is a framework in nonequilibrium statistical mechanics that successfully describes a wide variety of complex systems, including hydrodynamic turbulence, weakly-collisional plasmas, cosmic rays, power grid fluctuations, among several others. In this work we analyze the class of nonequilibrium steady-state systems consisting of a subsystem and its environment, and where the subsystem is described by the superstatistical framework. In this case we provide an answer to the mechanism by which a broad distribution of temperature arises, namely due to correlation between subsystem and environment. We prove that there is a unique microscopic definition $\mathcal{B}$ of inverse temperature compatible with superstatistics, in the sense that all moments of $\mathcal{B}$ and $β=1/(k_B T)$ coincide. The function $\mathcal{B}$ however, cannot depend on the degrees of freedom of the system itself, only on the environment, in full agreement with our previous impossibility theorem [Physica A \textbf{505}, 864-870 (2018)]. The present results also constrain the possible joint ensembles of system and environment compatible with superstatistics.

preprint2016arXiv

Applications of the divergence theorem in Bayesian inference and MaxEnt

Given a probability density $P({\bf x}|{\boldsymbol λ})$, where $\bf x$ represents continuous degrees of freedom and $λ$ a set of parameters, it is possible to construct a general identity relating expectations of observable quantities, which is a generalization of the equipartition theorem in Thermodynamics. In this work we explore some of the consequences of this relation, both in the context of sampling distributions and Bayesian posteriors, and how it can be used to extract some information without the need for explicit calculation of the partition function (or the Bayesian evidence, in the case of posterior expectations). Together with the general family of fluctuation theorems it constitutes a powerful tool for Bayesian/MaxEnt problems.

preprint2016arXiv

Continuity equation for probability as a requirement of inference over paths

In this work we present the fundamental ideas of inference over paths, and show how this formalism implies the continuity equation, which is central for the derivation of the main partial differential equations that constitute non-equilibrium statistical mechanics. Equations such as the Liouville equation, Fokker-Planck equation, among others can be recovered as particular cases of the continuity equation, under different probability fluxes. We derive the continuity equation in its most general form through what we call the \bf time-slicing equation, which lays down the procedure to go from the representation in terms of a path probability functional $ρ[X()]$ to a time-dependent probability density $ρ(x; t)$. The original probability functional $ρ[X()]$ can in principle be constructed from different methods of inference; in this work we sketch an application using the maximum path entropy or maximum caliber principle.

preprint2016arXiv

Exponential Family Models from Bayes' Theorem under Expectation Constraints

It is shown that a consistent application of Bayesian updating from a prior probability density to a posterior using evidence in the form of expectation constraints leads to exactly the same results as the application of the maximum entropy principle, namely a posterior belonging to the exponential family. The Bayesian updating procedure presented in this work is not expressed as a variational principle, and does not involve the concept of entropy. Therefore it conceptually constitutes a complete alternative to entropic methods of inference.

preprint2016arXiv

Liouville's Theorem from the Principle of Maximum Caliber in Phase Space

One of the cornerstones in non--equilibrium statistical mechanics (NESM) is Liouville's theorem, a differential equation for the phase space probability $ρ(q,p; t)$. This is usually derived considering the flow in or out of a given surface for a physical system (composed of atoms), via more or less heuristic arguments. In this work, we derive the Liouville equation as the partial differential equation governing the dynamics of the time-dependent probability $ρ(q, p; t)$ of finding a "particle" with Lagrangian $L(q, \dot{q}; t)$ in a specific point $(q, p)$ in phase space at time $t$, with $p=\partial L/\partial \dot{q}$. This derivation depends only on considerations of inference over a space of continuous paths. Because of its generality, our result is valid not only for "physical" systems but for any model depending on constrained information about position and velocity, such as time series.

preprint2016arXiv

Plausible inference from the idea of estimation

The probability axioms by R. T. Cox can be regarded as the modern foundations of Bayesian inference, the idea of assigning degrees of belief to logical propositions in a manner consistent with Boolean logic. In this work it is shown that one can start from an alternative point of view, postulating the existence of an operation for \textit{estimating magnitudes given a state of knowledge}. The properties of this operation can be established through basic consistency requirements. Application of this estimation operation to the truth values of propositions in an integer representation directly yields non-negative, normalized quantities following the sum and product rules of probability.

preprint2015arXiv

Hamiltonian formalism and path entropy maximization

Maximization of the path information entropy is a clear prescription for constructing models in non-equilibrium statistical mechanics. Here it is shown that, following this prescription under the assumption of arbitrary instantaneous constraints on position and velocity, a Lagrangian emerges which determines the most probable trajectory. Deviations from the probability maximum can be consistently described as slices in time by a Hamiltonian, according to a nonlinear Langevin equation and its associated Fokker-Planck equation. The connections unveiled between the maximization of path entropy and the Langevin/Fokker-Planck equations imply that missing information about the phase space coordinate never decreases in time, a purely information-theoretical version of the Second Law of Thermodynamics. All of these results are independent of any physical assumptions, and thus valid for any generalized coordinate as a function of time, or any other parameter. This reinforces the view that the Second Law is a fundamental property of plausible inference.

preprint2014arXiv

A Bayesian interpretation of abrupt phase transitions

The formalism used in describing the thermodynamics of abrupt (or first-order) phase transitions is reviewed as an application of maximum entropy inference. In this treatment, we show that the concepts of transition temperature, latent heat and entropy difference between phases will inevitably have an equivalent in any problem of inferring the result of a yes/no question, given information in the form of expectation values.

preprint2014arXiv

A maximum entropy model for opinions in social groups

We study how the opinions of a group of individuals determine their spatial distribution and connectivity, through an agent-based model. The interaction between agents is described by a Potts-like Hamiltonian in which agents are allowed to move freely without an underlying lattice (the average network topology connecting them is determined from the parameters). This kind of model was derived using maximum entropy statistical inference under fixed expectation values of certain probabilities that (we propose) are relevant to social organization. Control parameters emerge as Lagrange multipliers of the maximum entropy problem, and they can be associated with the level of consequence between the personal beliefs and external opinions, and the tendency to socialize with peers of similar or opposing views. These parameters define a phase diagram for the social system, which we studied using Monte Carlo Metropolis simulations. Our model presents both first and second-order phase transitions, depending on the ratio between the internal consequence and the interaction with others. We have found a critical value for the level of internal consequence, below which the personal beliefs of the agents seem to be irrelevant.

preprint2014arXiv

Fitting of interatomic potentials without forces: a parallel particle swarm optimization algorithm

We present a methodology for fitting interatomic potentials to ab initio data, using the particle swarm optimization (PSO) algorithm, needing only a set of positions and energies as input. The prediction error of energies associated with the fitted parameters can be close to 1 meV/atom or lower, for reference energies having a standard deviation of about 0.5 eV/atom. We tested our method by fitting a Sutton-Chen potential for copper from \emph{ab initio} data, which is able to recover structural and dynamical properties, and obtain a better agreement of the predicted melting point versus the experimental value, as compared to the prediction of the standard Sutton-Chen parameters.

preprint2014arXiv

Statistical distribution of thermal vacancies close to the melting point

A detailed description of the statistical distribution of thermal vacancies near the melting point is presented, using copper as an example. As the temperature is increased, the average number of thermal vacancies generated by atoms migrating to neighboring sites also increase, according to Arrhenius' law. We present for the first time a model for the distribution of thermal vacancies, which according to our results follow a Gamma distribution. All the simulations are carried out by classical molecular dynamics and the recognition of vacancies is achieved via a recently developed algorithm. Our results could be useful in the further development of a theory explaining the mechanism of homogeneous melting, which seems to be mediated (at least in part) by the accumulation of thermal vacancies near the melting point.

preprint2013arXiv

Newtonian Dynamics from the principle of Maximum Caliber

The foundations of Statistical Mechanics can be recovered almost in their entirety from the Principle of Maximum Entropy. In this work we show that its non-equilibrium generalization, the Principle of Maximum Caliber (Jaynes, 1980), when applied to the unknown trajectory followed by a particle, leads to Newton's second law under two quite intuitive assumptions (the expected square displacement in one step and the spatial probability distribution of the particle are known at all times). Our derivation explicitly highlights the role of mass as an emergent measure of the fluctuations in velocity (inertia) and the origin of potential energy as a manifestation of spatial correlations. According to our findings, the application of Newton's equations is not limited to mechanical systems, and therefore could be used in modelling ecological, financial and biological systems, among others.