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.

Finite-size effects in addition and chipping processes.

We investigate analytically and numerically a system of clusters evolving via collisions with clusters of minimal mass (monomers). Each collision either leads to the addition of the monomer to the cluster or the chipping of a monomer from the cluster, and emerging behaviors depend on which of the two processes is more probable. If addition prevails, monomers disappear in a time that scales as lnN with the total mass N≫1, and the system reaches a jammed state. When chipping prevails, the system remains in a quasistationary state for a time that scales exponentially with N, but eventually, a giant fluctuation leads to the disappearance of monomers. In the marginal case, monomers disappear in a time that scales linearly with N, and the final supercluster state is a peculiar jammed state; i.e., it is not extensive.

Range-Controlled Random Walks.

We introduce range-controlled random walks with hopping rates depending on the range N, that is, the total number of previously distinct visited sites. We analyze a one-parameter class of models with a hopping rate N^{a} and determine the large time behavior of the average range, as well as its complete distribution in two limit cases. We find that the behavior drastically changes depending on whether the exponent a is smaller, equal, or larger than the critical value, a_{d}, depending only on the spatial dimension d. When a>a_{d}, the forager covers the infinite lattice in a finite time. The critical exponent is a_{1}=2 and a_{d}=1 when d≥2. We also consider the case of two foragers who compete for food, with hopping rates depending on the number of sites each visited before the other. Surprising behaviors occur in 1D where a single walker dominates and finds most of the sites when a>1, while for a<1, the walkers evenly explore the line. We compute the gain of efficiency in visiting sites by adding one walker.

Dynamical fluctuations in the Riesz gas.

We consider an infinite system of particles on a line performing identical Brownian motions and interacting through the |x-y|^{-s} Riesz potential, causing the overdamped motion of particles. We investigate fluctuations of the integrated current and the position of a tagged particle. We show that for 01, the interactions are effectively short-ranged, and the universal subdiffusive t^{1/4} growth emerges with only amplitude depending on the exponent s. We also show that the two-time correlations of the tagged-particle position have the same form as for fractional Brownian motion.

Competition in a system of Brownian particles: Encouraging achievers.

We introduce and analytically and numerically study a simple model of interagent competition, where underachievement is strongly discouraged. We consider N≫1 particles performing independent Brownian motions on the line. Two particles are selected at random and at random times, and the particle closest to the origin is reset to it. We show that, in the limit of N→∞, the dynamics of the coarse-grained particle density field can be described by a nonlocal hydrodynamic theory which was encountered in a study of the spatial extent of epidemics in a critical regime. The hydrodynamic theory predicts relaxation of the system toward a stationary density profile of the "swarm" of particles, which exhibits a power-law decay at large distances. An interesting feature of this relaxation is a nonstationary "halo" around the stationary solution, which continues to expand in a self-similar manner. The expansion is ultimately arrested by finite-N effects at a distance of order sqrt[N] from the origin, which gives an estimate of the average radius of the swarm. The hydrodynamic theory does not capture the behavior of the particle farthest from the origin-the current leader. We suggest a simple scenario for typical fluctuations of the leader's distance from the origin and show that the mean distance continues to grow indefinitely as sqrt[t]. Finally, we extend the inter-agent competition from n=2 to an arbitrary number n of competing Brownian particles (n≪N). Our analytical predictions are supported by Monte Carlo simulations.

Divergence and consensus in majority rule.

We investigate majority rule dynamics in a population with two classes of people, each with two opinion states ±1, and with tunable interactions between people in different classes. In an update, a randomly selected group adopts the majority opinion if all group members belong to the same class; if not, majority rule is applied with rate ε. Consensus is achieved in a time that scales logarithmically with population size if ε≥ε_{c}=1/9. For ε<ε_{c}, the population can get trapped in a polarized state, with one class preferring the +1 state and the other preferring -1. The time to escape this polarized state and reach consensus scales exponentially with population size.

Blast in a One-Dimensional Cold Gas: From Newtonian Dynamics to Hydrodynamics.

A gas composed of a large number of atoms evolving according to Newtonian dynamics is often described by continuum hydrodynamics. Proving this rigorously is an outstanding open problem, and precise numerical demonstrations of the equivalence of the hydrodynamic and microscopic descriptions are rare. We test this equivalence in the context of the evolution of a blast wave, a problem that is expected to be at the limit where hydrodynamics could work. We study a one-dimensional gas at rest with instantaneous localized release of energy for which the hydrodynamic Euler equations admit a self-similar scaling solution. Our microscopic model consists of hard point particles with alternating masses, which is a nonintegrable system with strong mixing dynamics. Our extensive microscopic simulations find a remarkable agreement with Euler hydrodynamics, with deviations in a small core region that are understood as arising due to heat conduction.

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.

Random walk through a fertile site.

We study the dynamics of random walks hopping on homogeneous hypercubic lattices and multiplying at a fertile site. In one and two dimensions, the total number N(t) of walkers grows exponentially at a Malthusian rate depending on the dimensionality and the multiplication rate µ at the fertile site. When d>d_{c}=2, the number of walkers may remain finite forever for any µ; it surely remains finite when µ≤µ_{d}. We determine µ_{d} and show that 〈N(t)〉 grows exponentially if µ>µ_{d}. The distribution of the total number of walkers remains broad when d≤2, and also when d>2 and µ>µ_{d}. We compute 〈N^{m}〉 explicitly for small m, and show how to determine higher moments. In the critical regime, 〈N〉 grows as sqrt[t] for d=3, t/lnt for d=4, and t for d>4. Higher moments grow anomalously, 〈N^{m}〉∼〈N〉^{2m-1}, in the critical regime; the growth is normal, 〈N^{m}〉∼〈N〉^{m}, in the exponential phase. The distribution of the number of walkers in the critical regime is asymptotically stationary and universal, viz., it is independent of the spatial dimension. Interactions between walkers may drastically change the behavior. For random walks with exclusion, if d>2, there is again a critical multiplication rate, above which 〈N(t)〉 grows linearly (not exponentially) in time; when d≤d_{c}=2, the leading behavior is independent on µ and 〈N(t)〉 exhibits a sublinear growth.

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.

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.

Epidemics with containment measures.

We propose a tractable epidemic model that includes containment measures. In the absence of containment measures, the epidemics spread exponentially fast whenever the infectivity rate is positive λ>0. The containment measures are modeled by considering a time-dependent modulation of the bare infectivity λ leading to effective infectivity that decays in time for each infected individual, mimicking, for instance, the combined effect of the asymptomatic onset of the disease, testing policies, and quarantine. We consider a wide range of temporal kernels for effective infectivity, and we investigate the effect of the considered containment measures. We find that not all kernels are able to push the epidemic dynamics below the epidemic threshold with some containment measures only able to reduce the rate of the exponential growth of newly infected individuals. We also propose a pandemic model caused by a growing number of separated foci.

Large deviations in one-dimensional random sequential adsorption.

In random sequential adsorption (RSA), objects are deposited randomly, irreversibly, and sequentially; attempts leading to an overlap with previously deposited objects are discarded. The process continues until the system reaches a jammed state when no further additions are possible. We analyze a class of lattice RSA models in which landing on an empty site in a segment is allowed when at least b neighboring sites on the left and the right are unoccupied. For the minimal model (b=1), we compute the full counting statistics of the occupation number. We reduce the determination of the full counting statistics to a Riccati equation that appears analytically solvable only when b=1. We develop a perturbation procedure which, in principle, allows one to determine cumulants consecutively, and we compute the variance of the occupation number for all b.

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.

Exclusion in junction geometries.

We investigate the dynamics of the asymmetric exclusion process at a junction. When two input roads are initially fully occupied and a single output road is initially empty, the ensuing rarefaction wave has a rich spatial structure. The density profile also changes dramatically as the initial densities are varied. Related phenomenology arises when one road feeds into two. Finally, we determine the phase diagram of the open system, where particles are fed into two roads at rate α for each road, the two roads merge into one, and particles are extracted from the single output road at rate ß.

Occupation-time statistics of a gas of interacting diffusing particles.

The time that a diffusing particle spends in a certain region of space is known as the occupation time, or the residence time. Recently, the joint occupation-time statistics of an ensemble of noninteracting particles was addressed using the single-particle statistics. Here we employ the macroscopic fluctuation theory (MFT) to study the occupation-time statistics of many interacting particles. We find that interactions can significantly change the statistics and, in some models, even cause a singularity of the large-deviation function describing these statistics. This singularity can be interpreted as a dynamical phase transition. We also point out to a close relation between the MFT description of the occupation-time statistics of noninteracting particles and the level 2 large deviation formalism which describes the occupation-time statistics of a single particle.

Kinetics of deposition in the diffusion-controlled limit.

The adsorption of particles diffusing in a half space bounded by the substrate and irreversibly sticking to the substrate upon contacts is investigated. We show that when absorbing particles are planar disks diffusing in the three-dimensional half space, the coverage approaches its saturated "jamming" value as t^{-1} in the large time limit (generally as t^{-1/(d-1)} when the substrate is d dimensional and d>1, and as e^{-t/ln(t)} when d=1). We also analyze the asymptotic behavior when particles are spherical and when particles are planar aligned squares.

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.

Extinction and survival in two-species annihilation.

We study diffusion-controlled two-species annihilation with a finite number of particles. In this stochastic process, particles move diffusively, and when two particles of opposite type come into contact, the two annihilate. We focus on the behavior in three spatial dimensions and for initial conditions where particles are confined to a compact domain. Generally, one species outnumbers the other, and we find that the difference between the number of majority and minority species, which is a conserved quantity, controls the behavior. When the number difference exceeds a critical value, the minority becomes extinct and a finite number of majority particles survive, while below this critical difference, a finite number of particles of both species survive. The critical difference Δ_{c} grows algebraically with the total initial number of particles N, and when N≫1, the critical difference scales as Δ_{c}∼N^{1/3}. Furthermore, when the initial concentrations of the two species are equal, the average number of surviving majority and minority particles, M_{+} and M_{-}, exhibit two distinct scaling behaviors, M_{+}∼N^{1/2} and M_{-}∼N^{1/6}. In contrast, when the initial populations are equal, these two quantities are comparable M_{+}∼M_{-}∼N^{1/3}.

Optimally frugal foraging.

We introduce the frugal foraging model in which a forager performs a discrete-time random walk on a lattice in which each site initially contains S food units. The forager metabolizes one unit of food at each step and starves to death when it last ate S steps in the past. Whenever the forager eats, it consumes all food at its current site and this site remains empty forever (no food replenishment). The crucial property of the forager is that it is frugal and eats only when encountering food within at most k steps of starvation. We compute the average lifetime analytically as a function of the frugality threshold and show that there exists an optimal strategy, namely, an optimal frugality threshold k^{*} that maximizes the forager lifetime.