Electrostatically driven granular media: phase transitions and coarsening

Wave spectroscopy in a driven granular material

Driven granular media constitute model systems in out-of-equilibrium statistical physics. By assimilating the motions of granular particles to those of atoms, by analogy, one can obtain macroscopic equivalent of phase transitions. Here, we study fluid-like and crystal-like two-dimensional states in a driven granular material. In our experimental device, a tunable magnetic field induces and controls remote interactions between the granular particles. We use high-speed video recordings to analyse the velocity fluctuations of individual particles in stationary regime. Using statistical averaging, we find that the particles self-organize into collective excitations characterized by dispersion relations in the frequency-wavenumber space. These findings thus reveal that mechanical waves analogous to condensed matter phonons propagate in driven granular media. When the magnetic coupling is weak, the waves are longitudinal, as expected for a fluid-like phase. When the coupling is stronger, both longitudinal and transverse waves propagate, which is typically seen in solid-like phases. We model the dispersion relations using the spatial distribution of particles and their interaction potential. Finally, we infer the elastic parameters of the granular assembly from equivalent sound velocities, thus realizing the spectroscopy of a granular material.

Patterns in magnetic granular media at the crossover from two to three dimensions

We perform three-dimensional particle-based simulations of confined, vibrated, and magnetizable beads to study the effect of cell geometry on pattern selection. For quasi-two-dimensional systems, we reproduce previously observed macroscopic patterns such as hexagonal crystals and labyrinthine structures. For systems at the crossover from two to three dimensions, labyrinthine branches shorten and are replaced by triplets of beads forming upright triangles which self-organize into a herringbone pattern. This transition is associated with increases in both translational and orientational orders.

Plastic instabilities in charged granular systems: Competition between elasticity and electrostatics

Electrostatic theory preserves charges, but allows dipolar excitations. Elasticity theory preserves dipoles, but allows quadrupolar (Eshelby-like) plastic events. Charged amorphous granular systems are interesting in their own right; here we focus on their plastic instabilities and examine their mechanical response to external strain and to an external electric field, to expose the competition between elasticity and electrostatics. In this paper a generic model is offered, its mechanical instabilities are examined, and a theoretical analysis is presented. Plastic instabilities are discussed as saddle-node bifurcations that can be fully understood in terms of eigenvalues and eigenfunctions of the relevant Hessian matrix. This system exhibits moduli that describe how electric polarization and stress are influenced by strain and the electric field. Theoretical expression for these moduli are offered and compared to the measurements in numerical simulations.

Spreading of triboelectrically charged granular matter

We report on the spreading of triboelectrically charged glass particles on an oppositely charged surface of a plastic cylindrical container in the presence of a constant mechanical agitation. The particles spread via sticking, as a monolayer on the cylinder's surface. Continued agitation initiates a sequence of instabilities of this monolayer, which first forms periodic wavy-stripe-shaped transverse density modulation in the monolayer and then ejects narrow and long particle-jets from the tips of these stripes. These jets finally coalesce laterally to form a homogeneous spreading front that is layered along the spreading direction. These remarkable growth patterns are related to a time evolving frictional drag between the moving charged glass particles and the countercharges on the plastic container. The results provide insight into the multiscale time-dependent tribolelectric processes and motivates further investigation into the microscopic causes of these macroscopic dynamical instabilities and spatial structures.

Short-term forecasting of Taiwanese earthquakes using a universal model of fusion-fission processes

Predicting how large an earthquake can be, where and when it will strike remains an elusive goal in spite of the ever-increasing volume of data collected by earth scientists. In this paper, we introduce a universal model of fusion-fission processes that can be used to predict earthquakes starting from catalog data. We show how the equilibrium dynamics of this model very naturally explains the Gutenberg-Richter law. Using the high-resolution earthquake catalog of Taiwan between Jan 1994 and Feb 2009, we illustrate how out-of-equilibrium spatio-temporal signatures in the time interval between earthquakes and the integrated energy released by earthquakes can be used to reliably determine the times, magnitudes, and locations of large earthquakes, as well as the maximum numbers of large aftershocks that would follow.

Effect of grain size on the forces governing the dynamic behavior of electrostatically driven powder media

The aim of this paper is to identify the forces that govern the dynamic behavior of polydisperse granular systems under high electric fields in air media. In a previous paper, we reported the changes in the morphology and dynamics of grain structures under a frequency and amplitude variable voltage and summarized them in a diagram [L. M. Salvatierra et al., Chem. Phys. Lett. 481, 194 (2009)]. In the present article, we investigate the influence of the grain size distribution and the speed of voltage increase during tests. We explain the main transitions observed in the diagram in terms of the change in grain size and the frequency-dependent electrical response, related to electrophoretic and dielectrophoretic interactions.

Antiphase synchronization of electrically shaken conducting beads

When a spherical conducting bead is placed in an electrode, it experiences an electric force. In a plane capacitor, it can undergo a periodic bouncing between the electrodes. Using a fast video camera, we measured the acceleration of the bead and the period of its motion as a function of the applied voltage. A mathematical model based on the hypothesis of electrostatic equilibrium is proposed to describe the dynamics of the system. We observe a stabilization of the trajectories: A bead bouncing between two electrodes tends to oscillate on a quasivertical trajectory, whatever its initial horizontal velocity. When two identical beads are placed together in a capacitor, they oscillate at the same frequency and an antiphase synchronization effect occurs. We propose a simple mechanism based on a Kuramoto-like model to explain it.

Hydrodynamics of a vibrated granular monolayer

We investigate the long-standing puzzle of phase separation in a granular monolayer vibrated from below. Although this system is three dimensional, an interesting dynamics occurs mostly in the horizontal plane, perpendicular to the direction of vibration. Experiments [Olafsen and Urbach, Phys. Rev. Lett. 81, 4369 (1998)] demonstrated that for a high amplitude of vibration the system is in the gaslike phase, but when the amplitude becomes smaller than a certain threshold, a phase separation occurs: A solidlike dense condensate of particles forms in the center of the system, surrounded by particles in the gaslike phase. We explain theoretically the experimentally observed coexistence of dilute and dense phases, employing Navier-Stokes granular hydrodynamics. We show that the phase separation is associated with a negative compressibility of granular gas.

Precision assembly of oppositely and like-charged nanoobjects mediated by charge-induced dipole interactions

The range of electrostatic interactions controls precisely the mutual orientations of assembling charged nanoobjects. For nonspherically symmetric particles, polarization effects and induced dipoles can dominate charge-charge interactions. These charge-induced dipole interactions mediate orientation-specific aggregation of both oppositely and like-charged particles.

Comment on "Long-lived giant number fluctuations in a swarming granular nematic"

Narayan et al. (Reports, 6 July 2007, p. 105) reported giant number fluctuations attributed to curvature-driven active currents specific for nonequilibrium nematic systems. We present data demonstrating that similar results can be found in systems of spherical particles due either to inelastic clustering or persistent density inhomogeneity, suggesting two alternative explanations for their results.

Sublimation of a vibrated granular monolayer: coexistence of gas and solid

The fluidization of a monolayer of glass beads in a horizontally and vertically vibrated annular container is studied. At peak forcing accelerations between 1.1 and 1.5 g, a solidlike and a gaslike domain coexist. The solid fraction decreases with increasing acceleration and shows hysteresis. The sharp boundaries between the two regions travel around the channel faster than the particles are transported. Complementary to our experimental studies, a molecular dynamics simulation is used to extract local granular temperature and number density. It is found that the number density in the solid phase is several times that in the gas, while the temperature is orders of magnitude lower.

Velocity distributions of granular gases with drag and with long-range interactions

We study velocity statistics of electrostatically driven granular gases. For two different experiments, (i) nonmagnetic particles in a viscous fluid and (ii) magnetic particles in air, the velocity distribution is non-Maxwellian, and its high-energy tail is exponential, P(upsilon) approximately exp(-/upsilon/). This behavior is consistent with the kinetic theory of driven dissipative particles. For particles immersed in a fluid, viscous damping is responsible for the exponential tail, while for magnetic particles, long-range interactions cause the exponential tail. We conclude that velocity statistics of dissipative gases are sensitive to the fluid environment and to the form of the particle interaction.

Anisotropic dynamics in a shaken granular dimer gas experiment

The dynamics, velocity fluctuations, and particle-plate interactions for a two-dimensional granular gas of shaken, nonspherical particles are studied experimentally. The experiment consists of a horizontal plate that is vertically oscillated to drive the dynamics of macroscopic dimers, spherical pairs that are loosely connected by a rod that couple the interaction each of the spheres has with the shaking plate. The extended nature of the particles results in more than one energy-momentum transfer between the plate and each dimer per shaking cycle. This complex interaction results in anisotropic behavior for the dimer that is a function of the shaking parameters.

Structure formation in electromagnetically driven granular media

We report structure formation in submonolayers of magnetic microparticles subjected to periodic electrostatic and magnetic excitations. Depending on the excitation parameters, we observe the formation of a rich variety of structures: clusters, rings, chains, and networks. The dynamics and shapes of the structures are strongly dependent on the amplitude and frequency of the external magnetic field. We find that for pure ac magnetic driving the low-frequency magnetic excitation favors compact clusters, whereas high frequency driving favors chains and netlike structures. An abrupt phase transition from chains to a network phase was observed for a high density of particles.

Far-from-equilibrium Ostwald ripening in electrostatically driven granular powders

We report an experimental study of cluster size distributions in electrostatically driven granular submonolayers. The cluster size distribution in this far-from-equilibrium process exhibits dynamic scaling behavior characteristic of the (nearly equilibrium) Ostwald ripening, controlled by the attachment and detachment of the "gas" particles. The scaled size distribution, however, is different from the classical Wagner distribution obtained in the limit of a vanishingly small area fraction of the clusters. A much better agreement is found with the theory of Phys. Rev. E 65, 046117 (2002)] which accounts for the cluster merger.

Self-assembly and vortices formed by microparticles in weak electrolytes

We carried out experimental studies of the self-assembly of metallic micron-size particles in poorly conducting liquid subject to a constant electric field. Depending on the experimental conditions, the particles self-assemble into long chains directed along the electric field lines and form vortices and other structures. The vortices perform Brownian-type random motion due to self-induced chaotic hydrodynamic flows. We measured the diffusivity constant of the vortices and the conductivity and mechanical stiffness of the chains.

Pressure measurement in two-dimensional horizontal granular gases

A two-dimensional granular gas is produced by vibrating vertically a partial layer of beads on a horizontal plate. Measurements of the force applied by the granular gas to the sidewalls of the container, or granular pressure, are used to study the effect of the shaking strength, density, bead-plate restitution coefficient, and particle size on the steady properties of the gas.

Theory of pattern formation of metallic microparticles in poorly conducting liquids

We develop a continuum theory of self-assembly and pattern formation in metallic microparticles immersed in a poorly conducting liquid in a dc electric field. The theory is formulated in terms of conservation laws for the densities of immobile particles (precipitate) and bouncing particles (gas) coupled to the Navier-Stokes equation for the liquid. This theory successfully reproduced the correct topology of the phase diagram and primary patterns observed in the experiment [Phys. Rev. Lett. 90, 114301 (2003)]]: static crystals, honeycombs, dynamic pulsating rings, and rotating multipetal vortices.

Exchange-driven growth

We study a class of growth processes in which clusters evolve via exchange of particles. We show that depending on the rate of exchange there are three possibilities: (I) Growth-clusters grow indefinitely, (II) gelation-all mass is transformed into an infinite gel in a finite time, and (III) instant gelation. In regimes I and II, the cluster size distribution attains a self-similar form. The large size tail of the scaling distribution is Phi(x) approximately exp(-x(2-nu)), where nu is a homogeneity degree of the rate of exchange. At the borderline case nu=2, the distribution exhibits a generic algebraic tail, Phi(x) approximately x(-5). In regime III, the gel nucleates immediately and consumes the entire system. For finite systems, the gelation time vanishes logarithmically, T approximately [lnN](-(nu-2)), in the large system size limit N--> infinity. The theory is applied to coarsening in the infinite range Ising-Kawasaki model and in electrostatically driven granular layers.

Dynamic self-assembly and patterns in electrostatically driven granular media

We show that granular media, consisting of metallic microparticles immersed in a poorly conducting liquid in a strong dc electric field, self-assemble into a rich variety of novel phases. These phases include static precipitates: honeycombs and Wigner crystals; and novel dynamic condensates: toroidal vortices and pulsating rings. The observed structures are explained by the interplay between charged granular gas and electrohydrodynamic convective flows in the liquid.

Model of coarsening and vortex formation in vibrated granular rods

Neicu et al. observed experimentally spontaneous formation of the long-range orientational order and large-scale vortices in a system of vibrated macroscopic rods. We propose a phenomenological theory of this phenomenon based on a coupled system of equations for local rods density and tilt. The density evolution is described by the modified Cahn-Hilliard equation, while the tilt is described by the Ginzburg-Landau type equation. Our analysis shows that, in accordance with the Cahn-Hilliard dynamics, islands of the ordered phase appear spontaneously and grow due to coarsening. The generic vortex solutions of the Ginzburg-Landau equation for the tilt correspond to the vortical motion of the rods around the cores which are located near the centers of the islands.

Clustering transitions in vibrofluidized magnetized granular materials

We study the effects of long-range interactions on the phases observed in cohesive granular materials. At high vibration amplitudes, a gas of magnetized particles is observed with velocity distributions similar to nonmagnetized particles. Below a transition temperature compact clusters are observed to form and coexist with single particles. The cluster growth rate is consistent with a classical nucleation process. However, the temperature of the particles in the clusters is significantly lower than the surrounding gas, indicating a breakdown of equipartition. If the system is quenched to low temperatures, a metastable network of connected chains self-assemble due to the anisotropic nature of magnetic interactions between particles.

Coarsening of granular clusters: Two types of scaling behaviors

We report on an experimental study of small cluster dynamics during the coarsening process in driven granular submonolayers of 120-microm bronze particles. The techniques of electrostatic and vertical mechanical vibration were employed to excite the granular gas. We measure the scaling exponent for the evaporation of small clusters during coarsening. It was found that the surface area of small clusters S vs time t behaves as S to (t(0)-t)(2/3) for lower frequencies and S to (t(0)-t) for higher frequencies. We argue that the change in the scaling exponent is related to the transition from three-dimensional (3D) to 2D character of motion in the granular gas.

Velocity fluctuations in electrostatically driven granular media

We study experimentally the particle velocity fluctuations in an electrostatically driven dilute granular gas. The velocity distributions have strong deviations from a Maxwellian form over a wide range of parameters. We have found that the tails of the distribution functions are consistent with a stretched exponential law with typical exponents of the order 3/2. Molecular dynamic simulations shows qualitative agreement with experimental data. Our results suggest that this non-Gaussian behavior is typical of most inelastic gases with both short- and long-range interactions.

Symmetry breaking and coarsening of clusters in a prototypical driven granular gas

Granular hydrodynamics predicts symmetry-breaking instability in a two-dimensional ensemble of nearly elastically colliding smooth hard disks driven, at zero gravity, by a rapidly vibrating sidewall. Supercritical and subcritical symmetry-breaking bifurcations of the stripe state are identified, and the supercritical bifurcation curve is computed. The cluster dynamics proceed as a coarsening process mediated by the gas phase. Well above the bifurcation point the final steady state, selected by coarsening, represents a single strongly localized densely packed "droplet."

Phase separation and coarsening in electrostatically driven granular media

A continuum model for the phase separation and coarsening in electrostatically driven granular media is formulated in terms of a Ginzburg-Landau equation subject to conservation of the total number of grains. In the regime of well-developed clusters, the continuum model is used to derive "sharp-interface" equations that govern the dynamics of the interphase boundary. The model captures the essential physics of this system.

Phase ordering with a global conservation law: Ostwald ripening and coalescence

Globally conserved phase ordering dynamics is investigated in systems with short range correlations at t=0. A Ginzburg-Landau equation with a global conservation law is employed as the phase field model. The conditions are found under which the sharp-interface limit of this equation is reducible to the area-preserving motion by curvature. Numerical simulations show that, for both critical and off-critical quench, the equal-time pair correlation function exhibits dynamic scaling, and the characteristic coarsening length obeys l(t) approximately t(1/2). For the critical quench, our results are in excellent agreement with earlier results. For off-critical quench (Ostwald ripening) we investigate the dynamics of the size distribution function of the minority phase domains. The simulations show that, at large times, this distribution function has a self-similar form with growth exponent 1/2. The scaled distribution, however, strongly differs from the classical Wagner distribution. We attribute this difference to coalescence of domains. A theory of Ostwald ripening is developed that takes into account binary coalescence events. The theoretical scaled distribution function agrees well with that obtained in the simulations.

Normal scaling in globally conserved interface-controlled coarsening of fractal clusters

We find that globally conserved interface-controlled coarsening of diffusion-limited aggregates exhibits dynamic scale invariance (DSI) and normal scaling. This is demonstrated by a numerical solution of the Ginzburg-Landau equation with a global conservation law. The general sharp-interface limit of this equation is introduced and reduced to volume preserving motion by mean curvature. A simple example of globally conserved interface-controlled coarsening system: the sublimation/deposition dynamics of a solid and its vapor in a small closed vessel, is presented in detail. The results of the numerical simulations show that the scaled form of the correlation function has a power-law tail accommodating the fractal initial condition. The coarsening length exhibits normal dynamic scaling. A decrease of the cluster radius with time, predicted by DSI, is observed. The difference between global and local conservation is discussed.

Dynamics of electrostatically driven granular media: effects of humidity

We performed experimental studies of the effect of humidity on the dynamics of electrostatically driven granular materials. Both conducting and dielectric particles undergo a phase transition from an immobile state (granular solid) to a fluidized state (granular gas) with increasing applied field. Spontaneous precipitation of solid clusters from the gas phase occurs as the external driving is decreased. The clustering dynamics in conducting particles is primarily controlled by screening of the electric field but is aided by cohesion due to humidity. It is shown that humidity effects dominate the clustering process with dielectric particles.

Noise to order

Patterns in natural systems abound, from the stripes on a zebra to ripples in a riverbed. In many of these systems, the appearance of an ordered state is not unexpected as the outcome of an underlying ordered process. Thus crystal growth, honeycomb manufacture and floret evolution generate regular and predictable patterns. Intrinsically noisy and disordered processes such as thermal fluctuations or mechanically randomized scattering generate surprisingly similar patterns. Here we discuss some of the underlying mechanisms believed to be at the heart of these similarities.