Size distribution of particles in Saturn's rings from aggregation and fragmentation

Pattern formation by turbulent cascades

Fully developed turbulence is a universal and scale-invariant chaotic state characterized by an energy cascade from large to small scales at which the cascade is eventually arrested by dissipation1-6. Here we show how to harness these seemingly structureless turbulent cascades to generate patterns. Pattern formation entails a process of wavelength selection, which can usually be traced to the linear instability of a homogeneous state7. By contrast, the mechanism we propose here is fully nonlinear. It is triggered by the non-dissipative arrest of turbulent cascades: energy piles up at an intermediate scale, which is neither the system size nor the smallest scales at which energy is usually dissipated. Using a combination of theory and large-scale simulations, we show that the tunable wavelength of these cascade-induced patterns can be set by a non-dissipative transport coefficient called odd viscosity, ubiquitous in chiral fluids ranging from bioactive to quantum systems8-12. Odd viscosity, which acts as a scale-dependent Coriolis-like force, leads to a two-dimensionalization of the flow at small scales, in contrast with rotating fluids in which a two-dimensionalization occurs at large scales4. Apart from odd viscosity fluids, we discuss how cascade-induced patterns can arise in natural systems, including atmospheric flows13-19, stellar plasma such as the solar wind20-22, or the pulverization and coagulation of objects or droplets in which mass rather than energy cascades23-25.

Structured model conserving biomass for the size-spectrum evolution in aquatic ecosystems

Mathematical modelling of the evolution of the size-spectrum dynamics in aquatic ecosystems was discovered to be a powerful tool to have a deeper insight into impacts of human- and environmental driven changes on the marine ecosystem. In this article we propose to investigate such dynamics by formulating and investigating a suitable model. The underlying process for these dynamics is given by predation events, causing both growth and death of individuals, while keeping the total biomass within the ecosystem constant. The main governing equation investigated is deterministic and non-local of quadratic type, coming from binary interactions. Predation is assumed to strongly depend on the ratio between a predator and its prey, which is distributed around a preferred feeding preference value. Existence of solutions is shown in dependence of the choice of the feeding preference function as well as the choice of the search exponent, a constant influencing the average volume in water an individual has to search until it finds prey. The equation admits a trivial steady state representing a died out ecosystem, as well as-depending on the parameter-regime-steady states with gaps in the size spectrum, giving evidence to the well known cascade effect. The question of stability of these equilibria is considered, showing convergence to the trivial steady state in a certain range of parameters. These analytical observations are underlined by numerical simulations, with additionally exhibiting convergence to the non-trivial equilibrium for specific ranges of parameters.

Diffusion in multicomponent granular mixtures

We investigate diffusion in polydisperse granular media. We derive the mean-squared displacement of granular particles in a polydisperse granular gas in a homogeneous cooling state, containing an arbitrary amount of species of different sizes and masses. We investigate both models of constant and time-dependent restitution coefficients and obtain a universal law for the size dependence of the mean-squared displacement for steep size distributions.

Supercluster states and phase transitions in aggregation-fragmentation processes

We study the evolution of aggregates triggered by collisions with monomers that either lead to the attachment of monomers or the break-up of aggregates into constituting monomers. Depending on parameters quantifying addition and break-up rates, the system falls into a jammed or a steady state. Supercluster states (SCSs) are very peculiar nonextensive jammed states that also arise in some models. Fluctuations underlie the formation of the SCSs. Conventional tools, such as the van Kampen expansion, apply to small fluctuations. We go beyond the van Kampen expansion and determine a set of critical exponents quantifying SCSs. We observe continuous and discontinuous phase transitions between the states. Our theoretical predictions are in good agreement with numerical results.

Lattice models for ballistic aggregation: Cluster-shape-dependent exponents

We study ballistic aggregation on a two-dimensional square lattice, where particles move ballistically in between momentum and mass conserving coalescing collisions. Three models are studied based on the shapes of the aggregates: In the first the aggregates remain point particles, in the second they retain the fractal shape at the time of collision, and in the third they assume a spherical shape. The exponents describing the power-law temporal decay of number of particles and energy as well as dependence of velocity correlations on mass are determined using large-scale Monte Carlo simulations. It is shown that the exponents are universal only for the point-particle model. In the other two cases, the exponents are dependent on the initial number density and correlations vanish at high number densities. The fractal dimension for the second model is close to 1.49.

Universal Framework for Record Ages under Restart

We propose a universal framework to compute record age statistics of a stochastic time series that undergoes random restarts. The proposed framework makes minimal assumptions on the underlying process and is furthermore suited to treat generic restart protocols going beyond the Markovian setting. After benchmarking the framework for classical random walks on the 1D lattice, we derive a universal criterion underpinning the impact of restart on the age of the nth record for generic time series with nearest-neighbor transitions. Crucially, the criterion contains a penalty of order n that puts strong constraints on restart expediting the creation of records, as compared to the simple first-passage completion. The applicability of our approach is further demonstrated on an aggregation-shattering process where we compute the typical growth rates of aggregate sizes. This unified framework paves the way to explore record statistics of time series under restart in a wide range of complex systems.

Anomalous aggregation regimes of temperature-dependent Smoluchowski equations

Temperature-dependent Smoluchowski equations describe the ballistic agglomeration. In contrast to the standard Smoluchowski equations for the evolution of cluster densities, with constant rate coefficients, the temperature-dependent equations describe both-the evolution of the densities as well as cluster temperatures, which determine the agglomeration rates. To solve these equations, we develop a Monte Carlo technique based on the low-rank approximation for the aggregation kernel. Using this highly effective approach, we perform a comprehensive study of the kinetic phase diagram of the system and reveal a few surprising regimes, including permanent temperature growth and "density separation" regime, with a large gap in the size distribution for middle-size clusters. We perform scaling analysis and classify the aggregation kernels for the temperature-dependent equations. Furthermore, we conjecture the lack of gelation in such systems. The results of the scaling theory agree well with the simulation data.

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.

Collision of nanoparticles of covalently bound atoms: Impact of stress-dependent adhesion

The impact of nanoparticles (NPs) composed of atoms with covalent bonding is investigated numerically and theoretically. We use recent models of covalent bonding of carbon atoms and elaborate a numerical model of amorphous carbon (a-C) NPs, which may be applied for modeling soot particles. We compute the elastic moduli of the a-C material which agree well with the available data. We reveal an interesting phenomenon-stress-dependent adhesion, which refers to stress-enhanced formation of covalent bonds between contacting surfaces. We observe that the effective adhesion coefficient linearly depends on the maximal stress between the surfaces and explain this dependence. We compute the normal restitution coefficient for colliding NPs and explore the dependence of the critical velocity, demarcating bouncing and aggregative collisions, on the NP radius. Using the obtained elastic and stress-dependent adhesive coefficients we develop a theory for the critical velocity. The predictions of the theory agree very well with the simulation results.

Nonextensive Supercluster States in Aggregation with Fragmentation

Systems evolving through aggregation and fragmentation may possess an intriguing supercluster state (SCS). Clusters constituting this state are mostly very large, so the SCS resembles a gelling state, but the formation of the SCS is controlled by fluctuations and in this aspect, it is similar to a critical state. The SCS is nonextensive, that is, the number of clusters varies sublinearly with the system size. In the parameter space, the SCS separates equilibrium and jamming (extensive) states. The conventional methods, such as, e.g., the van Kampen expansion, fail to describe the SCS. To characterize the SCS we propose a scaling approach with a set of critical exponents. Our theoretical findings are in good agreement with numerical results.

Curvature flows, scaling laws and the geometry of attrition under impacts

Impact induced attrition processes are, beyond being essential models of industrial ore processing, broadly regarded as the key to decipher the provenance of sedimentary particles. Here we establish the first link between microscopic, particle-based models and the mean field theory for these processes. Based on realistic computer simulations of particle-wall collision sequences we first identify the well-known damage and fragmentation energy phases, then we show that the former is split into the abrasion phase with infinite sample lifetime (analogous to Sternberg's Law) at finite asymptotic mass and the cleavage phase with finite sample lifetime, decreasing as a power law of the impact velocity (analogous to Basquin's Law). This splitting establishes the link between mean field models (curvature-driven partial differential equations) and particle-based models: only in the abrasion phase does shape evolution emerging in the latter reproduce with startling accuracy the spatio-temporal patterns (two geometric phases) predicted by the former.

Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. II. Stability analysis

Conditions for the stability under linear perturbations around the homogeneous cooling state are studied for dilute granular gases of inelastic and rough hard disks or spheres with constant coefficients of normal (α) and tangential (β) restitution. After a formally exact linear stability analysis of the Navier-Stokes-Fourier hydrodynamic equations in terms of the translational (d_{t}) and rotational (d_{r}) degrees of freedom, the transport coefficients derived in the companion paper [A. Megías and A. Santos, "Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. I. Transport coefficients" Phys. Rev. E 104, 034901 (2021)10.1103/PhysRevE.104.034901] are employed. Known results for hard spheres [Garzó, Santos, and Kremer, Phys. Rev. E 97, 052901 (2018)10.1103/PhysRevE.97.052901] are recovered by setting d_{t}=d_{r}=3, while novel results for hard disks (d_{t}=2, d_{r}=1) are obtained. In the latter case, a high-inelasticity peculiar region in the (α,β) parameter space is found, inside which the critical wave number associated with the longitudinal modes diverges. Comparison with event-driven molecular dynamics simulations for dilute systems of hard disks at α=0.2 shows that this theoretical region of absolute instability may be an artifact of the extrapolation to high inelasticity of the approximations made in the derivation of the transport coefficients, although it signals a shrinking of the conditions for stability. In the case of moderate inelasticity (α=0.7), however, a good agreement between the theoretical predictions and the simulation results is found.

Hopf bifurcation in addition-shattering kinetics

In aggregation-fragmentation processes, a steady state is usually reached. This indicates the existence of an attractive fixed point in the underlying infinite system of coupled ordinary differential equations. The next simplest possibility is an asymptotically periodic motion. Never-ending oscillations have not been rigorously established so far, although oscillations have been recently numerically detected in a few systems. For a class of addition-shattering processes, we provide convincing numerical evidence for never-ending oscillations in a certain region U of the parameter space. The processes which we investigate admit a fixed point that becomes unstable when parameters belong to U and never-ending oscillations effectively emerge through a Hopf bifurcation.

Packing fraction of clusters formed in free-falling granular streams based on flash x-ray radiography

We study the packing fraction of clusters in free-falling streams of spherical and irregularly shaped particles using flash x-ray radiography. The estimated packing fraction of clusters is low enough to correspond to coordination numbers less than 6. Such coordination numbers in numerical simulations correspond to aggregates that collide and grow without bouncing. Moreover, the streams of irregular particles evolved faster and formed clusters of larger sizes with lower packing fraction. This result on granular streams suggests that particle shape has a significant effect on the agglomeration process of granular materials.

Rheology of two-dimensional crushable granular materials

The rheology of two-dimensional crushable granular materials under shear is numerically studied using the discrete element method. We find that the mean fragment size changes as the shear strain increases while the shear stress is almost independent of this mean size. The fragment size distribution is found to follow a power law. In particular, the exponent in the intermediate fragment size regime becomes approximately – 11/6, which is almost independent of the shear rate.

Energy nonequipartition in a collisional model of a confined quasi-two-dimensional granular mixture

A collisional model of a confined quasi-two-dimensional granular mixture is considered to analyze homogeneous steady states. The model includes an effective mechanism to transfer the kinetic energy injected by vibration in the vertical direction to the horizontal degrees of freedom of grains. The set of Enskog kinetic equations for the velocity distribution functions of each component is derived first to analyze the homogeneous state. As in the one-component case, an exact scaling solution is found where the time dependence of the distribution functions occurs entirely through the granular temperature T. As expected, the kinetic partial temperatures T_{i} of each component are different and, hence, energy equipartition is broken down. In the steady state, explicit expressions for the temperature T and the ratio of partial kinetic temperatures T_{i}/T_{j} are obtained by considering Maxwellian distributions defined at the partial temperatures T_{i}. The (scaled) granular temperature and the temperature ratios are given in terms of the coefficients of restitution, the solid volume fraction, the (scaled) parameters of the collisional model, and the ratios of mass, concentration, and diameters. In the case of a binary mixture, the theoretical predictions are exhaustively compared with both direct simulation Monte Carlo and molecular dynamics simulations with a good agreement. The deviations are identified to be originated in the non-Gaussianity of the velocity distributions and on microsegregation patterns, which induce spatial correlations not captured in the Enskog theory.

Role of energy in ballistic agglomeration

We study a ballistic agglomeration process in the reaction-controlled limit. Cluster densities obey an infinite set of Smoluchowski rate equations, with rates dependent on the average particle energy. The latter is the same for all cluster species in the reaction-controlled limit and obeys an equation depending on densities. We express the average energy through the total cluster density that allows us to reduce the governing equations to the standard Smoluchowski equations. We derive basic asymptotic behaviors and verify them numerically. We also apply our formalism to the agglomeration of dark matter.

Plato's cube and the natural geometry of fragmentation

Plato envisioned Earth's building blocks as cubes, a shape rarely found in nature. The solar system is littered, however, with distorted polyhedra-shards of rock and ice produced by ubiquitous fragmentation. We apply the theory of convex mosaics to show that the average geometry of natural two-dimensional (2D) fragments, from mud cracks to Earth's tectonic plates, has two attractors: "Platonic" quadrangles and "Voronoi" hexagons. In three dimensions (3D), the Platonic attractor is dominant: Remarkably, the average shape of natural rock fragments is cuboid. When viewed through the lens of convex mosaics, natural fragments are indeed geometric shadows of Plato's forms. Simulations show that generic binary breakup drives all mosaics toward the Platonic attractor, explaining the ubiquity of cuboid averages. Deviations from binary fracture produce more exotic patterns that are genetically linked to the formative stress field. We compute the universal pattern generator establishing this link, for 2D and 3D fragmentation.

Size-polydisperse dust in molecular gas: Energy equipartition versus nonequipartition

We investigate numerically and analytically size-polydisperse granular mixtures immersed into a molecular gas. We show that the equipartition of granular temperatures of particles of different sizes is established; however, the granular temperatures significantly differ from the temperature of the molecular gas. This result is surprising since, generally, the energy equipartition is strongly violated in driven granular mixtures. Qualitatively, the obtained results do not depend on the collision model, being valid for a constant restitution coefficient ɛ, as well as for the ɛ for viscoelastic particles. Our findings may be important for astrophysical applications, such as protoplanetary disks, interstellar dust clouds, and comets.

Temperature distribution in driven granular mixtures does not depend on mechanism of energy dissipation

We study analytically and numerically the distribution of granular temperatures in granular mixtures for different dissipation mechanisms of inelastic inter-particle collisions. Both driven and force-free systems are analyzed. We demonstrate that the simplified model of a constant restitution coefficient fails to predict even qualitatively a granular temperature distribution in a homogeneous cooling state. At the same time we reveal for driven systems a stunning result - the distribution of temperatures in granular mixtures is universal. That is, it does not depend on a particular dissipation mechanism of inter-particles collisions, provided the size distributions of particles is steep enough. The results of the analytic theory are compared with simulation results obtained by the direct simulation Monte Carlo (DSMC). The agreement between the theory and simulations is perfect. The reported results may have important consequences for fundamental science as well as for numerous application, e.g. for the experimental modelling in a lab of natural processes.

Critical properties of cluster size distribution in an asymmetric diffusion-aggregation model

We consider a stochastic dynamics for a system of diffusing hard-core particles on a periodic chain with asymmetric diffusion rules. A given cluster of particles can diffuse to the right as a whole but the particle located on the left boundary of the cluster is allowed to break off from it and diffuse to the left. Clusters of particles can eventually merge with other clusters. These rules allow for the creation of clusters of different sizes. We discuss the size distribution of the clusters in the long time or steady state limit, as a function of the particle concentration c. We consider the general time dependent master equation based on Smoluchowski's theory for local cluster merging or fragmentation and diffusion processes, and study the solutions using the generating function in the large size limit. We found that there exists a critical density c^{*}=sqrt[2]-1 for which the cluster distribution decays like a power law with exponent 5/2.

Jamming and tiling in fragmentation of rectangles

We investigate a stochastic process where a rectangle breaks into smaller rectangles through a series of horizontal and vertical fragmentation events. We focus on the case where both the vertical size and the horizontal size of a rectangle are discrete variables. Because of this constraint, the system reaches a jammed state where all rectangles are sticks, that is, rectangles with minimal width. Sticks are frozen as they cannot break any further. The average number of sticks in the jammed state, S, grows as S≃A/sqrt[2πlnA] with rectangle area A in the large-area limit, and remarkably, this behavior is independent of the aspect ratio. The distribution of stick length has a power-law tail, and further, its moments are characterized by a nonlinear spectrum of scaling exponents. We also study an asymmetric breakage process where vertical and horizontal fragmentation events are realized with different probabilities. In this case, there is a phase transition between a weakly asymmetric phase where the length distribution is independent of system size and a strongly asymmetric phase where this distribution depends on system size.

Dynamical Collective Memory in Fluidized Granular Materials

Recent experiments with rotational diffusion of a probe in a vibrated granular media revealed a rich scenario, ranging from a dilute gas to a dense liquid with cage effects and an unexpected superdiffusive behavior at large times. Here we set up a simulation that reproduces quantitatively the experimental observations and allows us to investigate the properties of the host granular medium, a task not feasible in the experiment. We discover a persistent collective rotational mode which emerges at a high density and a low granular temperature: a macroscopic fraction of the medium slowly rotates, randomly switching direction after very long times. Such a rotational mode of the host medium is the origin of the probe's superdiffusion. Collective motion is accompanied by a kind of dynamical heterogeneity at intermediate times (in the cage stage) followed by a strong reduction of fluctuations at late times, when superdiffusion sets in.

Electrification in granular gases leads to constrained fractal growth

The empirical observation of aggregation of dielectric particles under the influence of electrostatic forces lies at the origin of the theory of electricity. The growth of clusters formed of small grains underpins a range of phenomena from the early stages of planetesimal formation to aerosols. However, the collective effects of Coulomb forces on the nonequilibrium dynamics and aggregation process in a granular gas - a model representative of the above physical processes - have so far evaded theoretical scrutiny. Here, we establish a hydrodynamic description of aggregating granular gases that exchange charges upon collisions and interact via the long-ranged Coulomb forces. We analytically derive the governing equations for the evolution of granular temperature, charge variance, and number density for homogeneous and quasi-monodisperse aggregation. We find that, once the aggregates are formed, the granular temperature of the cluster population, the charge variance of the cluster population and the number density of the cluster population evolve in such a way that their non-dimensional combination obeys a physical constraint of nearly constant dimensionless ratio of characteristic electrostatic to kinetic energy. This constraint on the collective evolution of charged clusters is confirmed both by  our theory and our detailed molecular dynamics simulations. The inhomogeneous aggregation of monomers and clusters in their mutual electrostatic field proceeds in a fractal manner. Our theoretical framework is extendable to more precise charge exchange mechanisms, a current focus of extensive experimentation. Furthermore, it illustrates the collective role of long-ranged interactions in dissipative gases and can lead to novel designing principles in particulate systems.

Steady oscillations in aggregation-fragmentation processes

We report surprising steady oscillations in aggregation-fragmentation processes. Oscillating solutions are observed for the class of aggregation kernels K_{i,j}=i^{ν}j^{μ}+j^{ν}i^{μ} homogeneous in masses i and j of merging clusters and fragmentation kernels, F_{ij}=λK_{ij}, with parameter λ quantifying the intensity of the disruptive impacts. We assume a complete decomposition (shattering) of colliding partners into monomers. We show that an assumption of a steady-state distribution of cluster sizes, compatible with governing equations, yields a power law with an exponential cutoff. This prediction agrees with simulation results when θ≡ν-μ<1. For θ=ν-μ>1, however, the densities exhibit an oscillatory behavior. While these oscillations decay for not very small λ, they become steady if θ is close to 2 and λ is very small. Simulation results lead to a conjecture that for θ<1 the system has a stable fixed point, corresponding to the steady-state density distribution, while for any θ>1 there exists a critical value λ_{c}, such that for λ<λ_{c}, the system has an attracting limit cycle. This is rather striking for a closed system of Smoluchowski-like equations, lacking any sinks and sources of mass.

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.

Stationary mass distribution and nonlocality in models of coalescence and shattering

We study the asymptotic properties of the steady state mass distribution for a class of collision kernels in an aggregation-shattering model in the limit of small shattering probabilities. It is shown that the exponents characterizing the large and small mass asymptotic behavior of the mass distribution depend on whether the collision kernel is local (the aggregation mass flux is essentially generated by collisions between particles of similar masses) or nonlocal (collision between particles of widely different masses give the main contribution to the mass flux). We show that the nonlocal regime is further divided into two subregimes corresponding to weak and strong nonlocality. We also observe that at the boundaries between the local and nonlocal regimes, the mass distribution acquires logarithmic corrections to scaling and calculate these corrections. Exact solutions for special kernels and numerical simulations are used to validate some nonrigorous steps used in the analysis. Our results show that for local kernels, the scaling solutions carry a constant flux of mass due to aggregation, whereas for the nonlocal case there is a correction to the constant flux exponent. Our results suggest that for general scale-invariant kernels, the universality classes of mass distributions are labeled by two parameters: the homogeneity degree of the kernel and one further number measuring the degree of the nonlocality of the kernel.

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.

Aggregation-fragmentation and individual dynamics of active clusters

A remarkable feature of active matter is the propensity to self-organize. One striking instance of this ability to generate spatial structures is the cluster phase, where clusters broadly distributed in size constantly move and evolve through particle exchange, breaking or merging. Here we propose an exhaustive description of the cluster dynamics in apolar active matter. Exploiting large statistics gathered on thousands of Janus colloids, we measure the aggregation and fragmentation rates and rationalize the resulting cluster size distribution and fluctuations. We also show that the motion of individual clusters is entirely consistent with a model positing random orientation of colloids. Our findings establish a simple, generic model of cluster phase, and pave the way for a thorough understanding of clustering in active matter.

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.

Phase transitions in systems with aggregation and shattering

We consider a system of clusters made of elementary building blocks, monomers, and evolving via collisions between diffusing monomers and immobile composite clusters. In our model, the cluster-monomer collision can lead to the attachment of the monomer to the cluster (addition process) or to the total breakup of the cluster (shattering process). A phase transition, separating qualitatively different behaviors, occurs when the probability of shattering events exceeds a certain threshold. The novel feature of the phase transition is the dramatic dependence on the initial conditions.

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

We study far-from-equilibrium dynamics in models of freely cooling granular gas and ballistically aggregating compact clusters. For both the cases, from event-driven molecular dynamics simulations, we have presented detailed results on structure and dynamics in space dimensions d=1 and 2. Via appropriate analyses it has been confirmed that the ballistic aggregation mechanism applies in d=1 granular gases as well. Aging phenomena for this mechanism, in both the dimensions, have been studied via the two-time density autocorrelation function. This quantity is demonstrated to exhibit scaling property similar to that in the standard phase transition kinetics. The corresponding functional forms have been quantified and the outcomes have been discussed in connection with the structural properties. Our results on aging establish a more complete equivalence between the granular gas and the ballistic aggregation models in d=1.

Aggregation and fragmentation in liquids with dispersed nanoparticles

Nanoparticle-induced aggregation and fragmentation phenomena in liquid media are investigated by applying a model of preferential attachment of dispersing molecules to randomly chosen nanoparticles and larger particles, each containing a single nanoparticle.

Oscillations in Aggregation-Shattering Processes

We observe never-ending oscillations in systems undergoing collision-controlled aggregation and shattering. Specifically, we investigate aggregation-shattering processes with aggregation kernels K_{i,j}=(i/j)^{a}+(j/i)^{a} and shattering kernels F_{i,j}=λK_{i,j}, where i and j are cluster sizes, and parameter λ quantifies the strength of shattering. When 0≤a<1/2, there are no oscillations, and the system monotonically approaches a steady state for all values of λ; in this region, we obtain an analytical solution for the stationary cluster size distribution. Numerical solutions of the rate equations show that oscillations emerge in the 1/2<a≤1 range. When λ is sufficiently large, oscillations decay and eventually disappear, while for λ<λ_{c}(a), oscillations apparently persist forever. Thus, never-ending oscillations can arise in closed aggregation-shattering processes without sinks and sources of particles.

Kinetics of Vapor-Solid Phase Transitions: Structure, Growth, and Mechanism

The kinetics of the separation between low and high density phases in a single component Lennard-Jones model is studied via molecular dynamics simulations, at very low temperatures, in the space dimension d=2. For densities close to the vapor branch of the coexistence curve, disconnected nanoscale clusters of the high density phase exhibit essentially ballistic motion. Starting from nearly circular shapes, at the time of nucleation, these clusters grow via sticky collisions, gaining filamentlike nonequilibrium structure at a later time, with a very low fractal dimensionality. The origin of the latter is shown to lie in the low mobility of the constituent particles, in the corresponding cluster reference frame, due to the (quasi-long-range) crystalline order. Standard self-similarity in the domain pattern, typically observed in the kinetics of phase transitions, is found to be absent. This invalidates the common method, that provides a growth law comparable to that in solid mixtures, of quantifying growth. An appropriate alternative approach, involving the fractality, quantifies the growth of the characteristic "length" to be a power law with time, the exponent being strongly temperature dependent. The observed growth law is in agreement with the outcome of a nonequilibrium kinetic theory.

Genuine non-self-averaging and ultraslow convergence in gelation

In irreversible aggregation processes droplets or polymers of microscopic size successively coalesce until a large cluster of macroscopic scale forms. This gelation transition is widely believed to be self-averaging, meaning that the order parameter (the relative size of the largest connected cluster) attains well-defined values upon ensemble averaging with no sample-to-sample fluctuations in the thermodynamic limit. Here, we report on anomalous gelation transition types. Depending on the growth rate of the largest clusters, the gelation transition can show very diverse patterns as a function of the control parameter, which includes multiple stochastic discontinuous transitions, genuine non-self-averaging and ultraslow convergence of the transition point. Our framework may be helpful in understanding and controlling gelation.

Coefficient of restitution mediated by a diamagnetic interaction

We study how a magnetic bead bounces onto a horizontal diamagnetic conducting plane. The bead, falling down by gravity from a certain height, produces an Eddy current that creates a repelling force. For low velocities the bead is trapped by the surface, for intermediate ones it escapes. In such a case the induced current changes its sign, and so does the force. The balance between diamagnetic and viscoelastic interactions determines the bouncing dynamics. We find experimentally the restitution coefficient as a function of the impact speed of the bead and develop, taking into account simple energetic considerations, a model able to reproduce our findings.

Kinetic theory for dilute cohesive granular gases with a square well potential

We develop the kinetic theory of dilute cohesive granular gases in which the attractive part is described by a square well potential. We derive the hydrodynamic equations from the kinetic theory with the microscopic expressions for the dissipation rate and the transport coefficients. We check the validity of our theory by performing the direct simulation Monte Carlo.