1
|
Dyachenko RR, Matveev SA, Krapivsky PL. Finite-size effects in addition and chipping processes. Phys Rev E 2023; 108:044119. [PMID: 37978711 DOI: 10.1103/physreve.108.044119] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Accepted: 09/22/2023] [Indexed: 11/19/2023]
Abstract
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.
Collapse
Affiliation(s)
- R R Dyachenko
- Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, Moscow, 119991, Russia
- Marchuk Institute of Numerical Mathematics RAS, Moscow, 119333, Russia
| | - S A Matveev
- Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, Moscow, 119991, Russia
- Marchuk Institute of Numerical Mathematics RAS, Moscow, 119333, Russia
| | - P L Krapivsky
- Department of Physics, Boston University, Boston, Massachusetts 02215, USA
- Santa Fe Institute, Santa Fe, New Mexico 87501, USA
| |
Collapse
|
2
|
Esenturk E, Connaughton C. Role of zero clusters in exchange-driven growth with and without input. Phys Rev E 2020; 101:052134. [PMID: 32575316 DOI: 10.1103/physreve.101.052134] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2019] [Accepted: 04/16/2020] [Indexed: 11/07/2022]
Abstract
The exchange-driven growth model describes the mean-field kinetics of a population of composite particles (clusters) subject to pairwise exchange interactions. Exchange in this context means that upon interaction of two clusters, one loses a constituent unit (monomer) and the other gains this unit. Two variants of the exchange-driven growth model appear in applications. They differ in whether clusters of zero size are considered active or passive. In the active case, clusters of size zero can acquire a monomer from clusters of positive size. In the passive case they cannot, meaning that clusters reaching size zero are effectively removed from the system. We show that the large-time behavior is very different for the two variants of the model. We first consider an isolated system. In the passive case, the cluster size distribution tends towards a self-similar evolution and the typical cluster size grows as a power of time. In the active case, we identify a broad class of kernels for which the the cluster size distribution tends to a nontrivial time-independent equilibrium in which the typical cluster size is finite. We next consider a nonisolated system in which monomers are input at a constant rate. In the passive case, the cluster size distribution again attains a self-similar profile in which the typical cluster size grows as a power of time. In the active case, a surprising new behavior is found: the cluster size distribution asymptotes to the same equilibrium profile found in the isolated case but with an amplitude that increases linearly with time.
Collapse
Affiliation(s)
- Emre Esenturk
- Mathematics Institute, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom.,Warwick Manufacturing Group, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Colm Connaughton
- Mathematics Institute, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom.,Centre for Complexity Science, University of Warwick, Coventry CV4 7AL, United Kingdom.,London Mathematical Laboratory, 8 Margravine Gardens, London W6 8RH, United Kingdom
| |
Collapse
|
3
|
Grosskinsky S, Jatuviriyapornchai W. Derivation of mean-field equations for stochastic particle systems. Stoch Process Their Appl 2019. [DOI: 10.1016/j.spa.2018.05.006] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
4
|
Matsoukas T. Abrupt percolation in small equilibrated networks. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:052105. [PMID: 26066117 DOI: 10.1103/physreve.91.052105] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2014] [Indexed: 06/04/2023]
Abstract
Networks can exhibit an abrupt transition in the form of a spontaneous self-organization of a sizable fraction of the population into a giant component of connected members. This behavior has been demonstrated in random graphs under suppressive rules that passively or actively attempt to delay the formation of the giant cluster. We show that suppressive rules are not a necessary condition for a sharp transition at the percolation threshold. Rather, a finite system with aggressive tendency to form a giant cluster may exhibit an instability at the percolation threshold that is relieved through an abrupt and discontinuous transition to the stable branch. We develop the theory for a class of equilibrated networks that produce this behavior and find that the discontinuous jump is especially pronounced in small networks but disappears when the size of the system is infinite.
Collapse
Affiliation(s)
- Themis Matsoukas
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
5
|
Cremer P, Löwen H. Scaling of cluster growth for coagulating active particles. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:022307. [PMID: 25353472 DOI: 10.1103/physreve.89.022307] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2013] [Indexed: 06/04/2023]
Abstract
Cluster growth in a coagulating system of active particles (such as microswimmers in a solvent) is studied by theory and simulation. In contrast to passive systems, the net velocity of a cluster can have various scalings dependent on the propulsion mechanism and alignment of individual particles. Additionally, the persistence length of the cluster trajectory typically increases with size. As a consequence, a growing cluster collects neighboring particles in a very efficient way and thus amplifies its growth further. This results in unusual large growth exponents for the scaling of the cluster size with time and, for certain conditions, even leads to "explosive" cluster growth where the cluster becomes macroscopic in a finite amount of time.
Collapse
Affiliation(s)
- Peet Cremer
- Institut für Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universität Düsseldorf, D-40225 Düsseldorf, Germany
| | - Hartmut Löwen
- Institut für Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universität Düsseldorf, D-40225 Düsseldorf, Germany
| |
Collapse
|
6
|
Waclaw B, Evans MR. Explosive condensation in a mass transport model. PHYSICAL REVIEW LETTERS 2012; 108:070601. [PMID: 22401186 DOI: 10.1103/physrevlett.108.070601] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2011] [Indexed: 05/31/2023]
Abstract
We study a far-from-equilibrium system of interacting particles, hopping between sites of a 1D lattice with a rate which increases with the number of particles at interacting sites. We find that clusters of particles, which initially spontaneously form in the system, begin to move at increasing speed as they gain particles. Ultimately, they produce a moving condensate which comprises a finite fraction of the mass in the system. We show that, in contrast with previously studied models of condensation, the relaxation time to steady state decreases as an inverse power of lnL with system size L and that condensation is instantaneous for L→∞.
Collapse
Affiliation(s)
- Bartłomiej Waclaw
- SUPA, School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom
| | | |
Collapse
|
7
|
Ball RC, Connaughton C, Stein THM, Zaboronski O. Instantaneous gelation in Smoluchowski's coagulation equation revisited. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:011111. [PMID: 21867117 DOI: 10.1103/physreve.84.011111] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2010] [Revised: 04/26/2011] [Indexed: 05/31/2023]
Abstract
We study the solutions of the Smoluchowski coagulation equation with a regularization term which removes clusters from the system when their mass exceeds a specified cutoff size, M. We focus primarily on collision kernels which would exhibit an instantaneous gelation transition in the absence of any regularization. Numerical simulations demonstrate that for such kernels with monodisperse initial data, the regularized gelation time decreases as M increases, consistent with the expectation that the gelation time is zero in the unregularized system. This decrease appears to be a logarithmically slow function of M, indicating that instantaneously gelling kernels may still be justifiable as physical models despite the fact that they are highly singular in the absence of a cutoff. We also study the case when a source of monomers is introduced in the regularized system. In this case a stationary state is reached. We present a complete analytic description of this regularized stationary state for the model kernel, K(m(1),m(2)) = max{m(1),m(2)}(ν), which gels instantaneously when M → ∞ if ν>1. The stationary cluster size distribution decays as a stretched exponential for small cluster sizes and crosses over to a power law decay with exponent ν for large cluster sizes. The total particle density in the stationary state slowly vanishes as [(ν-1)log M](-1/2) when M → ∞. The approach to the stationary state is nontrivial: Oscillations about the stationary state emerge from the interplay between the monomer injection and the cutoff, M, which decay very slowly when M is large. A quantitative analysis of these oscillations is provided for the addition model which describes the situation in which clusters can only grow by absorbing monomers.
Collapse
Affiliation(s)
- Robin C Ball
- Centre for Complexity Science, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom.
| | | | | | | |
Collapse
|
8
|
Aarão Reis FDA, Stinchcombe RB. Unusual features of coarsening when detachment rates decrease with cluster mass. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:041411. [PMID: 18517621 DOI: 10.1103/physreve.77.041411] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2007] [Indexed: 05/26/2023]
Abstract
We study conserved one-dimensional models of particle diffusion, attachment, and detachment from clusters, where the detachment rates decrease with increasing cluster size as gamma(m) approximately m(-k), k>0 . Heuristic scaling arguments based on random walk properties show that the typical cluster size scales as (t/ln t)z , with z=1/(k+2) . The coarsening of neighboring clusters is characterized by initial symmetric flux of particles between them followed by an effectively asymmetric flux due to the unbalanced detachment rates, which leads to the above logarithmic corrections. Small clusters have densities of order t(-mz)(1), with z(1)=k/(k+2) . Thus for k<1 , the small clusters (mass of order unity) are statistically dominant and the average cluster size does not scale as the size of typically large clusters does. We also solve the master equation of the model under an independent interval approximation, which yields cluster distributions and exponent relations and gives the correct dominant coarsening exponent after suitable changes to incorporate effects of correlations. The coarsening of typical large clusters is described by the distribution Pt(m) approximately 1/ty f(m/tz) , with y=2z . All results are confirmed by simulation, which also illustrates the unusual features of cluster size distributions, with a power-law decay for small masses and a negatively skewed peak in the scaling region. The detachment rates considered here can apply in the presence of strong attractive interactions, and recent applications suggest that even more rapid rate decays are also physically realistic.
Collapse
Affiliation(s)
- F D A Aarão Reis
- Instituto de Física, Universidade Federal Fluminense, Avenida Litorânea s/n, 24210-340 Niterói RJ, Brazil.
| | | |
Collapse
|
9
|
Wang H, Lin Z, Ke J. Competition between the catalyzed birth and death in the exchange-driven growth. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 75:046108. [PMID: 17500962 DOI: 10.1103/physreve.75.046108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2006] [Indexed: 05/15/2023]
Abstract
We propose a three-species ( A , B , and C ) exchange-driven aggregate growth model with competition between catalyzed birth and catalyzed death. In the system, exchange-driven aggregation occurs between any two aggregates of the same species with the size-dependent rate kernel Kn(k,j)=Knkj (n=1,2,3) , and, meanwhile, monomer birth and death of species A occur under the catalysis of species B and C with the catalyzed birth and catalyzed death rate kernels I(k,j)=Ikjv and J(k,j)=Jkjv , respectively. The kinetic behavior is investigated by means of the mean-field rate equation approach. The form of the aggregate size distribution ak(t) of species A is found to depend crucially on the competition between species- B -catalyzed birth of species A and species- C -catalyzed death of species A , as well as the exchange-driven growth. The results show that (i) when exchange-driven aggregation dominates the process, ak(t) satisfies the conventional scaling form; (ii) when catalyzed birth dominates the process, ak(t) takes the conventional or generalized scaling form; and (iii) when catalyzed death dominates the process, the aggregate size distribution of species A evolves only according to some modified scaling forms.
Collapse
Affiliation(s)
- Haifeng Wang
- Department of Physics, Wenzhou University, Wenzhou 325027, China and Department of Physics, Qufu Normal University, Qufu 273165, China
| | | | | |
Collapse
|
10
|
Ke J, Chen X, Lin Z, Zheng Y, Lu W. Kinetics of migration-driven aggregation processes on scale-free networks. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:056102. [PMID: 17279963 DOI: 10.1103/physreve.74.056102] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/09/2006] [Indexed: 05/13/2023]
Abstract
We propose a solvable model for the migration-driven aggregate growth on completely connected scale-free networks. A reversible migration system is considered with the produce rate kernel K(k;l|i;j) approximately k(u)i(upsilon)(lj)(nu) or the generalized kernel K(k;l|i;j) approximately (k(upsilon)i(omega)+k(omega)i(upsilon)(lj)(nu), at which an i-mer aggregate locating on the node with j links gains one monomer from a k-mer aggregate locating on the node with l links. It is found that the evolution behavior of the system depends crucially on the details of the rate kernel. In some cases, the aggregate size distribution approaches a scaling form and the typical size S(t,l) of the aggregates locating on the nodes with l links grows infinitely with time; while in other cases, a gelation transition may emerge in the system at a finite critical time. We also introduce a simplified model, in which the aggregates independently gain or lose one monomer at the rate I(1)(k;l)=I(2)(k;l) proportional to k(omega)l(nu), and find the similar results. Most intriguingly, these models exhibit that the evolution behavior of the total distribution of the aggregates with the same size is drastically different from that for the corresponding system in normal space. We test our analytical results with the population data of all counties in the U.S. during the past century and find good agreement between the theoretical predictions and the realistic data.
Collapse
Affiliation(s)
- Jianhong Ke
- School of Physics and Electronic Information, Wenzhou University, Wenzhou 325027, China.
| | | | | | | | | |
Collapse
|
11
|
Lin Z, Ke J, Ye G. Mutually catalyzed birth of population and assets in exchange-driven growth. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:046113. [PMID: 17155139 DOI: 10.1103/physreve.74.046113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2005] [Revised: 07/28/2006] [Indexed: 05/12/2023]
Abstract
We propose an exchange-driven aggregation growth model of population and assets with mutually catalyzed birth to study the interaction between the population and assets in their exchange-driven processes. In this model, monomer (or equivalently, individual) exchange occurs between any pair of aggregates of the same species (population or assets). The rate kernels of the exchanges of population and assets are K(k,l) = Kkl and L(k,l) = Lkl , respectively, at which one monomer migrates from an aggregate of size k to another of size l. Meanwhile, an aggregate of one species can yield a new monomer by the catalysis of an arbitrary aggregate of the other species. The rate kernel of asset-catalyzed population birth is I(k,l) = Iklmu [and that of population-catalyzed asset birth is J(k,l) = Jklnu], at which an aggregate of size k gains a monomer birth when it meets a catalyst aggregate of size l . The kinetic behaviors of the population and asset aggregates are solved based on the rate equations. The evolution of the aggregate size distributions of population and assets is found to fall into one of three categories for different parameters mu and nu: (i) population (asset) aggregates evolve according to the conventional scaling form in the case of mu < or = 0 (nu < or = 0), (ii) population (asset) aggregates evolve according to a modified scaling form in the case of nu = 0 and mu > 0 (mu = 0 and nu > 0 ), and (iii) both population and asset aggregates undergo gelation transitions at a finite time in the case of mu = nu > 0.
Collapse
Affiliation(s)
- Zhenquan Lin
- Department of Physics, Wenzhou University, Wenzhou 325035, China.
| | | | | |
Collapse
|
12
|
Ke J, Lin Z, Zheng Y, Chen X, Lu W. Migration-driven aggregate growth on scale-free networks. PHYSICAL REVIEW LETTERS 2006; 97:028301. [PMID: 16907482 DOI: 10.1103/physrevlett.97.028301] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/09/2006] [Indexed: 05/11/2023]
Abstract
We study the kinetics of migration-driven aggregate growth on completely connected scale-free networks. A reversible migration system is considered with the size-dependent rate kernel K(k; l/i;j) approximately k(u)i(v)(lj)(v), at which an i-mer aggregate located on the node with j links gains one monomer from a k-mer aggregate on the node with l links. The results show that the evolution behavior of the aggregate size distribution is drastically different from that for the corresponding same system in normal space. This model can be used to mimic some phenomena such as the distribution of city populations. Moreover, we verify our analytic results in good agreement with the data of the population distributions of all U.S. counties.
Collapse
Affiliation(s)
- Jianhong Ke
- School of Physics and Electronic Information, Wenzhou University, Wenzhou 325027, China.
| | | | | | | | | |
Collapse
|
13
|
Ben-Naim E, Machta B, Machta J. Power-law velocity distributions in granular gases. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 72:021302. [PMID: 16196551 DOI: 10.1103/physreve.72.021302] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2005] [Indexed: 05/04/2023]
Abstract
The kinetic theory of granular gases is studied for spatially homogeneous systems. At large velocities, the equation governing the velocity distribution becomes linear, and it admits stationary solutions with a power-law tail, f (v) approximately v(-sigma) . This behavior holds in arbitrary dimension for arbitrary collision rates including both hard spheres and Maxwell molecules. Numerical simulations show that driven steady states with the same power-law tail can be realized by injecting energy into the system at very high energies. In one dimension, we also obtain self-similar time-dependent solutions where the velocities collapse to zero. At small velocities there is a steady state and a power-law tail but at large velocities, the behavior is time dependent with a stretched exponential decay.
Collapse
Affiliation(s)
- E Ben-Naim
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| | | | | |
Collapse
|
14
|
Aarão Reis FDA, Stinchcombe RB. Nonuniversal coarsening and universal distributions in far-from-equilibrium systems. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 71:026110. [PMID: 15783381 DOI: 10.1103/physreve.71.026110] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2004] [Indexed: 05/24/2023]
Abstract
Anomalous coarsening in far-from-equilibrium one-dimensional systems is investigated by applying simulation and analytic techniques to minimal hard-core particle (exclusion) models. They contain mechanisms of aggregated particle diffusion, with rates epsilon<<1 , particle deposition into cluster gaps, but suppressed for the smallest gaps, and breakup of clusters that are adjacent to large gaps. Cluster breakup rates vary with the cluster length x as k x(alpha) . The domain growth law x approximately (epsilont)(z) , with z=1/ (2+alpha) for alpha>0 , is explained by a simple scaling picture involving the time for two particles to coalesce and a new particle to be deposited. The density of double vacancies, at which deposition and cluster breakup are allowed, scales as 1/ [t (epsilont)(z) ] . Numerical simulations for several values of alpha and epsilon confirm these results. A fuller approach is presented which employs a mapping of cluster configurations to a column picture and an approximate factorization of the cluster configuration probability within the resulting master equation. The equation for a one-variable scaling function explains the above average cluster length scaling. The probability distributions of cluster lengths x scale as P (x) = [1/ (epsilont)(z) ] g (y) , with y identical with x/ (epsilont)(z) , which is confirmed by simulation. However, those distributions show a universal tail with the form g (y) approximately exp (- y(3/2) ) , which is explained by the connection of the vacancy dynamics with the problem of particle trapping in an infinite sea of traps. The high correlation of surviving particle displacement in the latter problem explains the failure of the independent cluster approximation to represent those rare events.
Collapse
Affiliation(s)
- F D A Aarão Reis
- Instituto de Física, Universidade Federal Fluminense, Avenida Litorânea s/n, 24210-340 Niterói RJ, Brazil.
| | | |
Collapse
|
15
|
Sapozhnikov MV, Peleg A, Meerson B, Aranson IS, Kohlstedt KL. Far-from-equilibrium Ostwald ripening in electrostatically driven granular powders. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 71:011307. [PMID: 15697595 DOI: 10.1103/physreve.71.011307] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/11/2004] [Indexed: 05/24/2023]
Abstract
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.
Collapse
Affiliation(s)
- M V Sapozhnikov
- Argonne National Laboratory, 9700 S. Cass Avenue, Argonne, Illinois 60439, USA
| | | | | | | | | |
Collapse
|
16
|
Ben-Naim E, Krapivsky PL. Size of outbreaks near the epidemic threshold. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 69:050901. [PMID: 15244801 DOI: 10.1103/physreve.69.050901] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2004] [Indexed: 05/24/2023]
Abstract
The spread of infectious diseases near the epidemic threshold is investigated. Scaling laws for the size and the duration of outbreaks originating from a single infected individual in a large susceptible population are obtained. The maximal size of an outbreak n(*) scales as N(2/3) with N the population size. This scaling law implies that the average outbreak size [n]scales as N(1/3). Moreover, the maximal and the average duration of an outbreak grow as t(*) approximately N(1/3) and [t] approximately ln N, respectively.
Collapse
Affiliation(s)
- E Ben-Naim
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| | | |
Collapse
|