Stochastic thermodynamics of opinion dynamics models.

We show that models of opinion formation and dissemination in a community of individuals can be framed within stochastic thermodynamics from which we can build a nonequilibrium thermodynamics of opinion dynamics. This is accomplished by decomposing the original transition rate that defines an opinion model into two or more transition rates, each representing the contact with heat reservoirs at different temperatures, and postulating an energy function. As the temperatures are distinct, heat fluxes are present even at the stationary state and linked to the production of entropy, the fundamental quantity that characterizes nonequilibrium states. We apply the present framework to a generic-vote model including the majority-vote model in a square lattice and in a cubic lattice. The fluxes and the rate of entropy production are calculated by numerical simulation and by the use of a pair approximation.

Stochastic motion in phase space on a surface of constant energy.

We study closed systems of particles that are subject to stochastic forces in addition to the conservative forces. The stochastic equations of motion are set up in such a way that the energy is strictly conserved at all times. To ensure this conservation law, the evolution equation for the probability density is derived using an appropriate interpretation of the stochastic equation of motion that is not the Itô nor the Stratonovic interpretation. The trajectories in phase space are restricted to the surface of constant energy. Despite this restriction, the entropy is shown to increase with time, expressing irreversible behavior and relaxation to equilibrium. This main result of the present approach contrasts with that given by the Liouville equation, which also describes closed systems, but does not show irreversibility.

Analysis of earlier times and flux of entropy on the majority voter model with diffusion.

We study the properties of nonequilibrium systems modelled as spin models without defined Hamiltonian as the majority voter model. This model has transition probabilities that do not satisfy the condition of detailed balance. The lack of detailed balance leads to entropy production phenomena, which are a hallmark of the irreversibility. By considering that voters can diffuse on the lattice we analyze how the entropy production and how the critical properties are affected by this diffusion. We also explore two important aspects of the diffusion effects on the majority voter model by studying entropy production and entropy flux via time-dependent and steady-state simulations. This study is completed by calculating some critical exponents as function of the diffusion probability.

Stochastic thermodynamics and entropy production of chemical reaction systems.

We investigate the nonequilibrium stationary states of systems consisting of chemical reactions among molecules of several chemical species. To this end, we introduce and develop a stochastic formulation of nonequilibrium thermodynamics of chemical reaction systems based on a master equation defined on the space of microscopic chemical states and on appropriate definitions of entropy and entropy production. The system is in contact with a heat reservoir and is placed out of equilibrium by the contact with particle reservoirs. In our approach, the fluxes of various types, such as the heat and particle fluxes, play a fundamental role in characterizing the nonequilibrium chemical state. We show that the rate of entropy production in the stationary nonequilibrium state is a bilinear form in the affinities and the fluxes of reaction, which are expressed in terms of rate constants and transition rates, respectively. We also show how the description in terms of microscopic states can be reduced to a description in terms of the numbers of particles of each species, from which follows the chemical master equation. As an example, we calculate the rate of entropy production of the first and second Schlögl reaction models.

Stochastic approach to equilibrium and nonequilibrium thermodynamics.

We develop the stochastic approach to thermodynamics based on stochastic dynamics, which can be discrete (master equation) and continuous (Fokker-Planck equation), and on two assumptions concerning entropy. The first is the definition of entropy itself and the second the definition of entropy production rate, which is non-negative and vanishes in thermodynamic equilibrium. Based on these assumptions, we study interacting systems with many degrees of freedom in equilibrium or out of thermodynamic equilibrium and how the macroscopic laws are derived from the stochastic dynamics. These studies include the quasiequilibrium processes; the convexity of the equilibrium surface; the monotonic time behavior of thermodynamic potentials, including entropy; the bilinear form of the entropy production rate; the Onsager coefficients and reciprocal relations; and the nonequilibrium steady states of chemical reactions.

Stochastic dynamics of dengue epidemics.

We use a stochastic Markovian dynamics approach to describe the spreading of vector-transmitted diseases, such as dengue, and the threshold of the disease. The coexistence space is composed of two structures representing the human and mosquito populations. The human population follows a susceptible-infected-recovered (SIR) type dynamics and the mosquito population follows a susceptible-infected-susceptible (SIS) type dynamics. The human infection is caused by infected mosquitoes and vice versa, so that the SIS and SIR dynamics are interconnected. We develop a truncation scheme to solve the evolution equations from which we get the threshold of the disease and the reproductive ratio. The threshold of the disease is also obtained by performing numerical simulations. We found that for certain values of the infection rates the spreading of the disease is impossible, for any death rate of infected mosquitoes.

Entropy production in nonequilibrium systems at stationary states.

We present a stochastic approach to nonequilibrium thermodynamics based on the expression of the entropy production rate advanced by Schnakenberg for systems described by a master equation. From the microscopic Schnakenberg expression we get the macroscopic bilinear form for the entropy production rate in terms of fluxes and forces. This is performed by placing the system in contact with two reservoirs with distinct sets of thermodynamic fields and by assuming an appropriate form for the transition rate. The approach is applied to an interacting lattice gas model in contact with two heat and particle reservoirs. On a square lattice, a continuous symmetry breaking phase transition takes place such that at the nonequilibrium ordered phase a heat flow sets in even when the temperatures of the reservoirs are the same. The entropy production rate is found to have a singularity at the critical point of the linear-logarithm type.

Probability currents and entropy production in nonequilibrium lattice systems.

The structure of probability currents is studied for the dynamical network after consecutive contraction on two-state, nonequilibrium lattice systems. This procedure allows us to investigate the transition rates between configurations on small clusters and highlights some relevant effects of lattice symmetries on the elementary transitions that are responsible for entropy production. A method is suggested to estimate the entropy production for different levels of approximations (cluster sizes) as demonstrated in the two-dimensional contact process with mutation.

Entropy production in irreversible systems described by a Fokker-Planck equation.

We analyze the irreversibility and the entropy production in nonequilibrium interacting particle systems described by a Fokker-Planck equation by the use of a suitable master equation representation. The irreversible character is provided either by nonconservative forces or by the contact with heat baths at distinct temperatures. The expression for the entropy production is deduced from a general definition, which is related to the probability of a trajectory in phase space and its time reversal, that makes no reference a priori to the dissipated power. Our formalism is applied to calculate the heat conductance in a simple system consisting of two Brownian particles each one in contact to a heat reservoir. We show also the connection between the definition of entropy production rate and the Jarzynski equality.

Critical behavior of the susceptible-infected-recovered model on a square lattice.

By means of numerical simulations and epidemic analysis, the transition point of the stochastic asynchronous susceptible-infected-recovered model on a square lattice is found to be c0=0.1765005(10), where c is the probability a chosen infected site spontaneously recovers rather than tries to infect one neighbor. This point corresponds to an infection/recovery rate of λ(c)=(1-c0)/c0=4.665 71(3) and a net transmissibility of (1-c0)/(1+3c0)=0.538 410(2), which falls between the rigorous bounds of the site and bond thresholds. The critical behavior of the model is consistent with the two-dimensional percolation universality class, but local growth probabilities differ from those of dynamic percolation cluster growth, as is demonstrated explicitly.

Role of noise in population dynamics cycles.

Noise is an intrinsic feature of population dynamics and plays a crucial role in oscillations called phase-forgetting quasicycles by converting damped into sustained oscillations. This function of noise becomes evident when considering Langevin equations whose deterministic part yields only damped oscillations. We formulate here a consistent and systematic approach to population dynamics, leading to a Fokker-Planck equation and the associate Langevin equations in accordance with this conceptual framework, founded on stochastic lattice-gas models that describe spatially structured predator-prey systems. Langevin equations in the population densities and predator-prey pair density are derived in two stages. First, a birth-and-death stochastic process in the space of prey and predator numbers and predator-prey pair number is obtained by a contraction method that reduces the degrees of freedom. Second, a van Kampen expansion in the inverse of system size is then performed to get the Fokker-Planck equation. We also study the time correlation function, the asymptotic behavior of which is used to characterize the transition from the cyclic coexistence of species to the ordinary coexistence.

Entropy production in the majority-vote model.

We analyzed the entropy production in the majority-vote model by using a mean-field approximation and Monte Carlo simulations. The dynamical rules of the model do not obey detailed balance so that entropy is continuously being produced. This nonequilibrium stochastic model is known to have a critical behavior belonging to the universality class of the equilibrium Ising model. We show that the entropy production exhibits a singularity at the critical point whose exponent is estimated numerically.

Stationary distribution of finite-size systems with absorbing states.

We introduce a procedure that allows us to obtain nontrivial stationary distributions of finite-size models with absorbing states. Two models are studied: the contact process and the sandpile model with height restriction. To avoid the permanence of the system in the absorbing state we create a small perturbation that drives the system out of the absorbing state. In the former model a particle is created, in the latter an active site is created. The stationary distributions around the critical point are analyzed by the use of finite-size scaling.

Sandpiles with height restrictions.

We study stochastic sandpile models with a height restriction in one and two dimensions. A site can topple if it has a height of two, as in Manna's model, but, in contrast to previously studied sandpiles, here the height (or number of particles per site), cannot exceed two. This yields a considerable simplification over the unrestricted case, in which the number of states per site is unbounded. Two toppling rules are considered: in one, the particles are redistributed independently, while the other involves some cooperativity. We study the fixed-energy system (no input or loss of particles) using cluster approximations and extensive simulations, and find that it exhibits a continuous phase transition to an absorbing state at a critical value zeta(c) of the particle density. The critical exponents agree with those of the unrestricted Manna sandpile.