A. Prados

Thermodynamic description in a simple model for granular compaction

A simple model for the dynamics of a granular system under tapping is studied. The model can be considered as a particularization for short taps of a more general one-dimensional lattice model with facilitated dynamics. The steady state reached by the system is discussed and the results are shown to be consistent with the thermodynamic granular theory developed by Edwards and coworkers. In particular, the basic assumption of the theory, i.e., that the probability distribution depends only on the volume of the configuration, is verified.

STOCHASTIC RESONANCE IN A ONE-DIMENSIONAL ISING MODEL

Abstract The stochastic resonance phenomenon for a one-dimensional Ising model in an oscillating magnetic field is discussed. Within the linear field approximation, the amplitude of the induced magnetization presents a maximum as a function of the temperature. The behaviour of the phase shift is also studied.

A Dynamical Monte Carlo Algorithm for Master Equations with Time-Dependent Transition Rates

A Monte Carlo algorithm for simulating master equations with time-dependent transition rates is described. It is based on a waiting time image, and takes into account that the system can become frozen when the transition rates tend to zero fast enough in time. An analytical justification is provided. The algorithm reduces to the Bortz-Kalos-Lebowitz one when the transition rates are constant. Since the exact evaluation of waiting times is rather involved in general, a simple and efficient iterative method for accurately calculating them is introduced. As an example, the algorithm is applied to a one-dimensional Ising system with Glauber dynamics. It is shown that it reproduces the exact analytical results, being more efficient than the direct implementation of the Metropolis algorithm

Simple model with facilitated dynamics for granular compaction

A simple lattice model is used to study compaction in granular media. As in real experiments, we consider a series of taps separated by large enough waiting times. The relaxation of the density exhibits the characteristic inverse logarithmic law. Moreover, we have been able to identify analytically the relevant time scale, leading to a relaxation law independent of the specific values of the parameters. Also, an expression for the asymptotic density reached in the compaction process has been derived. The theoretical predictions agree fairly well with the results from the Monte Carlo simulation.

Hysteresis in vibrated granular media

Some general dynamical properties of models for compaction of granular media based on master equations are analyzed. In particular, a one-dimensional lattice model with short-ranged dynamical constraints is considered. The stationary state is consistent with Edwards theory of powders. The system is submitted to processes in which the tapping strength is monotonically increased and decreased. In such processes the behavior of the model resembles the reversible–irreversible branches which have been recently observed in experiments. This behavior is understood in terms of the general dynamical properties of the model, and related to the hysteresis cycles exhibited by structural glasses in thermal cycles. The existence of a “normal” solution, i.e., a special solution of the master equation which is monotonically approached by all the other solutions, plays a fundamental role in the understanding of the hysteresis effects.

Large fluctuations in driven dissipative media.

We analyze the fluctuations of the dissipated energy in a simple and general model where dissipation, diffusion, and driving are the key ingredients. The full dissipation distribution, which follows from hydrodynamic fluctuation theory, shows non-Gaussian tails and no negative branch, thus violating the fluctuation theorem as expected from the irreversibility of the dynamics. It exhibits simple scaling forms in the weak- and strong-dissipation limits, with large fluctuations favored in the former case but strongly suppressed in the latter. The typical path associated with a given dissipation fluctuation is also analyzed in detail. Our results, confirmed in extensive simulations, strongly support the validity of hydrodynamic fluctuation theory to describe fluctuating behavior in driven dissipative media.

Slow logarithmic relaxation in models with hierarchically constrained dynamics.

A general kind of model with hierarchically constrained dynamics is shown to exhibit logarithmic anomalous relaxation similar to a variety of complex strongly interacting materials. The logarithmic behavior describes most of the decay of the response function.

Typical and rare fluctuations in nonlinear driven diffusive systems with dissipation.

We consider fluctuations of the dissipated energy in nonlinear driven diffusive systems subject to bulk dissipation and boundary driving. With this aim, we extend the recently introduced macroscopic fluctuation theory to nonlinear driven dissipative media, starting from the fluctuating hydrodynamic equations describing the system mesoscopic evolution. Interestingly, the action associated with a path in mesoscopic phase space, from which large-deviation functions for macroscopic observables can be derived, has the same simple form as in nondissipative systems. This is a consequence of the quasielasticity of microscopic dynamics, required in order to have a nontrivial competition between diffusion and dissipation at the mesoscale. Euler-Lagrange equations for the optimal density and current fields that sustain an arbitrary dissipation fluctuation are also derived. A perturbative solution thereof shows that the probability distribution of small fluctuations is always Gaussian, as expected from the central limit theorem. On the other hand, strong separation from the Gaussian behavior is observed for large fluctuations, with a distribution which shows no negative branch, thus violating the Gallavotti-Cohen fluctuation theorem, as expected from the irreversibility of the dynamics. The dissipation large-deviation function exhibits simple and general scaling forms for weakly and strongly dissipative systems, with large fluctuations favored in the former case but heavily suppressed in the latter. We apply our results to a general class of diffusive lattice models for which dissipation, nonlinear diffusion, and driving are the key ingredients. The theoretical predictions are compared to extensive numerical simulations of the microscopic models, and excellent agreement is found. Interestingly, the large-deviation function is in some cases nonconvex beyond some dissipation. These results show that a suitable generalization of macroscopic fluctuation theory is capable of describing in detail the fluctuating behavior of nonlinear driven dissipative media.

Nonlinear driven diffusive systems with dissipation: fluctuating hydrodynamics.

We consider a general class of nonlinear diffusive models with bulk dissipation and boundary driving and derive its hydrodynamic description in the large size limit. Both the average macroscopic behavior and the fluctuating properties of the hydrodynamic fields are obtained from the microscopic dynamics. This analysis yields a fluctuating balance equation for the local energy density at the mesoscopic level, characterized by two terms: (i) a diffusive term, with a current that fluctuates around its average behavior given by nonlinear Fouriers law, and (ii) a dissipation term which is a general function of the local energy density. The quasielasticity of microscopic dynamics, required in order to have a nontrivial competition between diffusion and dissipation in the macroscopic limit, implies a noiseless dissipation term in the balance equation, so dissipation fluctuations are enslaved to those of the density field. The microscopic complexity is thus condensed in just three transport coefficients-the diffusivity, the mobility, and a new dissipation coefficient-which are explicitly calculated within a local equilibrium approximation. Interestingly, the diffusivity and mobility coefficients obey an Einstein relation despite the fully nonequilibrium character of the problem. The general theory here presented is applied to a particular albeit broad family of systems, the simplest nonlinear dissipative variant of the so-called KMP model for heat transport. The theoretical predictions are compared to extensive numerical simulations, and an excellent agreement is found.

Lattice models for granular-like velocity fields: finite-size effects

Long-range spatial correlations in the velocity and energy fields of a granular fluid are discussed in the framework of a 1d lattice model. The dynamics of the velocity field occurs through nearest-neighbour inelastic collisions that conserve momentum but dissipate energy. A set of equations for the fluctuating hydrodynamics of the velocity and energy mesoscopic fields give a first approximation for (i) the velocity structure factor and (ii) the finite-size correction to the Haff law, both in the homogeneous cooling regime. At a more refined level, we have derived the equations for the two-site velocity correlations and the total energy fluctuations. First, we seek a perturbative solution thereof, in powers of the inverse of system size. On the one hand, when scaled with the granular temperature, the velocity correlations tend to a stationary value in the long time limit. On the other hand, the scaled standard deviation of the total energy diverges, that is, the system shows multiscaling. Second, we find an exact solution for the velocity correlations in terms of the spectrum of eigenvalues of a certain matrix. The results of numerical simulations of the microscopic model confirm our theoretical results, including the above described multiscaling phenomenon.