Hydrodynamics of a vibrated granular monolayer

Non-equilibrium coexistence between a fluid and a hotter or colder crystal of granular hard disks

Non-equilibrium phase coexistence is commonly observed in both biological and artificial systems, yet understanding it remains a significant challenge. Unlike equilibrium systems, where free energy provides a unifying framework, the absence of such a quantity in non-equilibrium settings complicates their theoretical understanding. Granular materials, driven out of equilibrium by energy dissipation during collisions, serve as an ideal platform to investigate these systems, offering insights into the parallels and distinctions between equilibrium and non-equilibrium phase behavior. For example, the coexisting dense phase is typically colder than the dilute phase, a result usually attributed to greater dissipation in denser regions. In this article, we demonstrate that this is not always the case. Using a simple numerical granular model, we show that a hot solid and a cold liquid can coexist in granular systems. This counterintuitive phenomenon arises because the collision frequency can be lower in the solid phase than in the liquid phase, consistent with equilibrium results for hard-disk systems. We further demonstrate that kinetic theory can be extended to accurately predict phase temperatures even at very high packing fractions, including within the solid phase. Our results highlight the importance of collisional dynamics and energy exchange in determining phase behavior in granular materials, offering new insights into non-equilibrium phase coexistence and the complex physics underlying granular systems.

Interplay between an Absorbing Phase Transition and Synchronization in a Driven Granular System

Absorbing phase transitions (APTs) are widespread in nonequilibrium systems, spanning condensed matter, epidemics, earthquakes, ecology, and chemical reactions. APTs feature an absorbing state in which the system becomes entrapped, along with a transition, either continuous or discontinuous, to an active state. Understanding which physical mechanisms determine the order of these transitions represents a challenging open problem in nonequilibrium statistical mechanics. Here, by numerical simulations and mean-field analysis, we show that a quasi-2D vibrofluidized granular system exhibits a novel form of APT. The absorbing phase is observed in the horizontal dynamics below a critical packing fraction, and can be continuous or discontinuous based on the emergent degree of synchronization in the vertical motion. Our results provide a direct representation of a feasible experimental scenario, showcasing a surprising interplay between dynamic phase transition and synchronization.

Ballistic propagation of density correlations and excess wall forces in quenched granular media

We investigate a granular gas in a shaken quasi-two-dimensional box in molecular dynamics computer simulations. After a sudden change (quench) of the shaking amplitude, transient density correlations are observed orders of magnitude beyond the steady-state correlation length scale. Propagation of the correlations is ballistic, in contrast to recently investigated quenches of Brownian particles that show diffusive propagation [Rohwer et al., Phys. Rev. Lett. 118, 015702 (2017)PRLTAO0031-900710.1103/PhysRevLett.118.015702, Rohwer et al., Phys. Rev. E 97, 032125 (2018)2470-004510.1103/PhysRevE.97.032125]. At sufficiently strong cooling of the fluid the effect is overlaid by clustering instability of the homogeneous cooling state with different scaling behavior. We are able to identify different quench regimes. In each regime correlations exhibit remarkably universal position dependence. In simulations performed with side walls we find confinement effects for temperature and pressure in steady-state simulations and an additional transient wall pressure contribution when changing the shaking amplitude. The transient contribution is ascribed to enhanced relaxation of the fluid in the presence of walls. From incompatible scaling behavior we conclude that the observed effects with and without side walls constitute distinct phenomena.

Understanding an instability in vibrated granular monolayers

We investigate the dynamics of an ensemble of smooth inelastic hard spheres confined between two horizontal plates separated by a distance smaller than twice the diameter of the particles, in such a way that the system is quasi-two-dimensional. The bottom wall is vibrating and, therefore, it injects energy into the system in the vertical direction and a stationary state is reached. It is found that if the size of the plates is small enough, the stationary state is homogeneous. Otherwise, a cluster of particles is developed. The instability is understood by using some effective hydrodynamic equations in the horizontal plane. Moreover, the theoretical prediction for the size of the system above which it becomes unstable agrees very well with molecular dynamics simulation results without any fitting parameter.

Nonequilibrium steady states, coexistence, and criticality in driven quasi-two-dimensional granular matter

Nonequilibrium steady states of vibrated inelastic frictionless spheres are investigated in quasi-two-dimensional confinement via molecular dynamics simulations. The phase diagram in the density-amplitude plane exhibits a fluidlike disordered and an ordered phase with threefold symmetry, as well as phase coexistence between the two. A dynamical mechanism exists that brings about metastable traveling clusters and at the same time stable clusters with anisotropic shapes at low vibration amplitude. Moreover, there is a square bilayer state which is connected to the fluid by BKTHNY-type two-step melting with an intermediate tetratic phase. The critical behavior of the two continuous transitions is studied in detail. For the fluid-tetratic transition, critical exponents of γ[over ̃]=1.73, η_{4}≈1/4, and z=2.05 are obtained.

Effective two-dimensional model for granular matter with phase separation

Granular systems confined in vertically vibrated shallow horizontal boxes (quasi-two-dimensional geometry) present a liquid-to-solid phase transition when the frequency of the periodic forcing is increased. An effective model, where grains move and collide in two-dimensions is presented, which reproduces the aforementioned phase transition. The key element is that besides the two-dimensional degrees of freedom, each grain has an additional variable ɛ that accounts for the kinetic energy stored in the vertical motion in the real quasi-two-dimensional motion. This energy grows monotonically during free flight, mimicking the energy gained by collisions with the vibrating walls and, at collisions, this energy is instantaneously transferred to the horizontal degrees of freedom. As a result, the average values of ɛ and the kinetic temperature are decreasing functions of the local density, giving rise to an effective pressure that can present van der Waals loops. A kinetic theory approach predicts the conditions that must satisfy the energy growth function to obtain the phase separation, which are verified with molecular dynamics simulations. Notably, the effective equation of state and the critical points computed considering the velocity-time-of-flight correlations differ only slightly from those obtained by simple kinetic theory calculations that neglect those correlations.

Boltzmann kinetic equation for a strongly confined gas of hard spheres

A Boltzmann-like kinetic equation for a quasi-two-dimensional gas of hard spheres is derived. The system is confined between two parallel hard plates separated a distance between one and two particle diameters. An entropy Lyapunov function for the equation is identified. In addition to the usual Boltzmann expression, it contains a contribution associated to the confinement of the particles. The steady properties of the system agree with equilibrium statistical mechanics results. Equations describing the energy transfer between the degrees of freedom parallel and perpendicular to the confining plates are obtained for some simple initial configurations. The theoretical predictions are compared with molecular dynamics simulation data and a good agreement is found.

Noise induces rare events in granular media

The granular Leidenfrost effect [B. Meerson, et al., Phys. Rev. Lett. 91, 024301 (2003)PRLTAO0031-900710.1103/PhysRevLett.91.024301; P. Eshuis et al., Phys. Rev. Lett. 95, 258001 (2005)PRLTAO0031-900710.1103/PhysRevLett.95.258001] is the levitation of a mass of granular matter when a wall below the grains is vibrated, giving rise to a hot granular gas below the cluster. We find by simulation that for a range of parameters the system is bistable: the levitated cluster can occasionally break and give rise to two clusters and a hot granular gas above and below. We use techniques from the theory of rare events to compute the mean transition time for breaking to occur. This requires the introduction of a two-component reaction coordinate.

Stability analysis of the homogeneous hydrodynamics of a model for a confined granular gas

The linear hydrodynamic stability of a model for confined quasi-two-dimensional granular gases is analyzed. The system exhibits homogeneous hydrodynamics, i.e., there are macroscopic evolution equations for homogeneous states. The stability analysis is carried out around all these states and not only the homogeneous steady state reached eventually by the system. It is shown that in some cases the linear analysis is not enough to reach a definite conclusion on the stability, and molecular dynamics simulation results are presented to elucidate these cases. The analysis shows the relevance of nonlinear hydrodynamic contributions to describe the behavior of spontaneous fluctuations occurring in the system, that lead even to the transitory formation of clusters of particles. The conclusion is that the system is always stable. The relevance of the results for describing the instabilities of confined granular gases observed experimentally is discussed.