Ballistic aggregation in systems of inelastic particles: Cluster growth, structure, and aging

Hydrodynamic effects in kinetics of phase separation in binary fluids: Critical versus off-critical compositions

Via hydrodynamics-preserving molecular dynamics simulations we study growth phenomena in a phase-separating symmetric binary mixture model. We quench high-temperature homogeneous configurations to state points inside the miscibility gap, for various mixture compositions. For compositions at the symmetric or critical value we capture the rapid linear viscous hydrodynamic growth due to advective transport of material through tubelike interconnected domains. For state points very close to any of the branches of the coexistence curve, the growth in the system, following nucleation of disconnected droplets of the minority species, occurs via a coalescence mechanism. Using state-of-the-art techniques, we have identified that these droplets, between collisions, exhibit diffusive motion. The value of the exponent for the power-law growth, related to this diffusive coalescence mechanism, has been estimated. While the exponent nicely agrees with that for the growth via the well-known Lifshitz-Slyozov particle diffusion mechanism, the amplitude is stronger. For the intermediate compositions we observe initial rapid growth that matches the expectations for viscous or inertial hydrodynamic pictures. However, at later times these types of growth cross over to the exponent that is decided by the diffusive coalescence mechanism.

Active particles in explicit solvent: Dynamics of clustering for alignment interaction

We study the dynamics of clustering in systems containing active particles that are immersed in an explicit solvent. For this, we have adopted a hybrid simulation method, consisting of molecular dynamics and multiparticle collision dynamics. In our model, the overlap-avoiding passive interaction of an active particle with another active particle or a solvent particle has been taken care of via variants of the Lennard-Jones potential. Dynamic interactions among the active particles have been incorporated via a Vicsek-like alignment rule in self-propulsion that facilitates clustering. We quantify the effects of activity and importance of hydrodynamics on the dynamics of clustering via variations of relevant system parameters. Results are obtained for low overall density of active particles, for which the state point is close to the vapor branch of the coexistence curve, and thus the morphology consists of disconnected clusters. In such a situation, the mechanism of growth switches among particle diffusion, diffusive coalescence, and ballistic aggregation, depending upon the presence or absence of active and hydrodynamic interactions providing different kinds of mobilities to the clusters. Corresponding growth laws have been quantified and discussed in the context of appropriate theoretical pictures. Our results suggest that multiparticle collision dynamics is an effective method for the investigation of hydrodynamic phenomena in phase-separating active matter systems.

How do clusters in phase-separating active matter systems grow? A study for Vicsek activity in systems undergoing vapor-solid transition

Via molecular dynamics simulations, we have studied the kinetics of vapor-"solid" phase transition in an active matter model in which self-propulsion is introduced via the well-known Vicsek rule. The overall density of the particles is chosen in such a way that the evolution morphology consists of disconnected clusters that are defined as regions of high density of particles. Our focus has been on understanding the influence of the above-mentioned self-propulsion on structure and growth of these clusters by comparing the results with those for the passive limit of the model that also exhibits vapor-"solid" transition. While in the passive case growth occurs due to a standard diffusive mechanism, the Vicsek activity leads to very rapid growth, via a process that is practically equivalent to the ballistic aggregation mechanism. The emerging growth law in the latter case has been accurately estimated and explained by invoking information on velocity and structural aspects of the clusters into a relevant theory. Some of these results are also discussed with reference to a model for active Brownian particles.

Kinetics of domain growth and aging in a two-dimensional off-lattice system

We have used molecular dynamics simulations for a comprehensive study of phase separation in a two-dimensional single-component off-lattice model where particles interact through the Lennard-Jones potential. Via state-of-the-art methods we have analyzed simulation data on structure, growth, and aging for nonequilibrium evolutions in the model. These data were obtained following quenches of well-equilibrated homogeneous configurations, with density close to the critical value, to various temperatures inside the miscibility gap, having vapor-"liquid" as well as vapor-"solid" coexistence. For the vapor-liquid phase separation we observe that ℓ, the average domain length, grows with time (t) as t^{1/2}, a behavior that has connection with hydrodynamics. At low-enough temperature, a sharp crossover of this time dependence to a much slower, temperature-dependent, growth is identified within the timescale of our simulations, implying "solid"-like final state of the high-density phase. This crossover is, interestingly, accompanied by strong differences in domain morphology and other structural aspects between the two situations. For aging, we have presented results for the order-parameter autocorrelation function. This quantity exhibits data collapse with respect to ℓ/ℓ_{w}, ℓ, and ℓ_{w} being the average domain lengths at times t and t_{w} (≤t), respectively, the latter being the age of a system. Corresponding scaling function follows a power-law decay: ∼(ℓ/ℓ_{w})^{-λ} for t≫t_{w}. The decay exponent λ, for the vapor-liquid case, is accurately estimated via the application of an advanced finite-size scaling method. The obtained value is observed to satisfy a bound.

Aging exponents for nonequilibrium dynamics following quenches from critical points

Via Monte Carlo simulations we study nonequilibrium dynamics in the nearest-neighbor Ising model, following quenches to points inside the ordered region of the phase diagram. With the broad objective of quantifying the nonequilibrium universality classes corresponding to spatially correlated and uncorrelated initial configurations, in this paper we present results for the decay of the order-parameter autocorrelation function for quenches from the critical point. This autocorrelation is an important probe for the aging dynamics in far-from-equilibrium systems and typically exhibits power-law scaling. From the state-of-the-art analysis of the simulation results, we quantify the corresponding exponents (λ) for both conserved and nonconserved (order-parameter) dynamics of the model in space dimension d=3. Via structural analysis we demonstrate that the exponents satisfy a bound. We also revisit the d=2 case to obtain more accurate results. It appears that irrespective of the dimension, λ is approximately the same for both conserved and nonconserved dynamics.

Initial correlation dependence of aging in phase separating solid binary mixtures and ordering ferromagnets

Following quenches of initial configurations having long range spatial correlations, prepared at the demixing critical point, to points inside the miscibility gap, we study aging phenomena in solid binary mixtures. Results on the decay of the two-time order-parameter autocorrelation functions, obtained from Monte Carlo simulations of the two-dimensional Ising model, with Kawasaki exchange kinetics, are analyzed via state-of-the art methods. The outcome is compared with that obtained for the ordering in uniaxial ferromagnets. For the latter, we have performed Monte Carlo simulations of the same model using the Glauber mechanism. For both types of systems we provide comparative discussion of our results with reference to those concerning quenches with configurations having no spatial correlation. We also discuss the role of structure on the decay of these correlations.

Aging phenomena during phase separation in fluids: decay of autocorrelation for vapor-liquid transitions

We performed molecular dynamics simulations to study relaxation phenomena during vapor-liquid transitions in a single component Lennard-Jones system. Results from two different overall densities are presented: one in the neighborhood of the vapor branch of the coexistence curve and the other being close to the critical density. The nonequilibrium morphologies, growth mechanisms and growth laws in the two cases are vastly different. In the low density case growth occurs via diffusive coalescence of droplets in a disconnected morphology. On the other hand, the elongated structure in the higher density case grows via advective transport of particles inside the tube-like liquid domains. The objective in this work has been to identify how the decay of the order-parameter autocorrelation, an important quantity to understand aging dynamics, differs in the two cases. In the case of the disconnected morphology, we observe a very robust power-law decay, as a function of the ratio of the characteristic lengths at the observation time and at the age of the system, whereas the results for the percolating structure appear rather complex. To quantify the decay in the latter case, unlike the standard method followed in a previous study, here we have performed a finite-size scaling analysis. The outcome of this analysis shows the presence of a strong preasymptotic correction, while revealing that in this case also, albeit in the asymptotic limit, the decay follows a power-law. Even though the corresponding exponents in the two cases differ drastically, this study, combined with a few recent ones, suggests that power-law behavior of this correlation function is rather universal in coarsening dynamics.

Dimension dependence of clustering dynamics in models of ballistic aggregation and freely cooling granular gas

Via event-driven molecular dynamics simulations we study kinetics of clustering in assemblies of inelastic particles in various space dimensions. We consider two models, viz., the ballistic aggregation model (BAM) and the freely cooling granular gas model (GGM), for each of which we quantify the time dependence of kinetic energy and average mass of clusters (that form due to inelastic collisions). These quantities, for both the models, exhibit power-law behavior, at least in the long time limit. For the BAM, corresponding exponents exhibit strong dimension dependence and follow a hyperscaling relation. In addition, in the high packing fraction limit the behavior of these quantities become consistent with a scaling theory that predicts an inverse relation between energy and mass. On the other hand, in the case of the GGM we do not find any evidence for such a picture. In this case, even though the energy decay, irrespective of packing fraction, matches quantitatively with that for the high packing fraction picture of the BAM, it is inversely proportional to the growth of mass only in one dimension, and the growth appears to be rather insensitive to the choice of the dimension, unlike the BAM.

Increasing temperature of cooling granular gases

The kinetic energy of a force-free granular gas decays monotonously due to inelastic collisions of the particles. For a homogeneous granular gas of identical particles, the corresponding decay of granular temperature is quantified by Haff’s law. Here, we report that for a granular gas of aggregating particles, the granular temperature does not necessarily decay but may even increase. Surprisingly, the increase of temperature is accompanied by the continuous loss of total gas energy. This stunning effect arises from a subtle interplay between decaying kinetic energy and gradual reduction of the number of degrees of freedom associated with the particles’ dynamics. We derive a set of kinetic equations of Smoluchowski type for the concentrations of aggregates of different sizes and their energies. We find scaling solutions to these equations and a condition for the aggregation mechanism predicting growth of temperature. Numerical direct simulation Monte Carlo results confirm the theoretical predictions.

Granular gases—dilute systems composed of dissipatively colliding particles—exhibit anomalous dynamics and numerous surprising phenomena. Here, Brilliantov et al. show that the aggregation mechanism can induce increase of the gas temperature despite the fact that the total kinetic energy decreases.

Early-stage aggregation in three-dimensional charged granular gas

Neutral grains made of the same dielectric material can attain considerable charges due to collisions and generate long-range interactions. We perform molecular dynamic simulations in three dimensions for a dilute, freely cooling granular gas of viscoelastic particles that exchange charges during collisions. As compared to the case of clustering of viscoelastic particles solely due to dissipation, we find that the electrostatic interactions due to collisional charging alter the characteristic size, morphology, and growth rate of the clusters. The average cluster size grows with time as a power law, whose exponent is relatively larger in the charged gas than the neutral case. The growth of the average cluster size is found to be independent of the ratio of characteristic Coulomb to kinetic energy, or equivalently, of the typical Bjerrum length. However, this ratio alters the crossover time of the growth. Both simulations and mean-field calculations based on Smoluchowski's equation suggest that a suppression of particle diffusion due to the electrostatic interactions helps in the aggregation process.