-
Catalytic Coagulation
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We introduce an autocatalytic aggregation model in which the rate at which two clusters merge to form a cluster is controlled by the presence of a third "catalytic" cluster whose mass must equal to the mass of one of the reaction partners. The catalyst is unaffected by the joining event and is available to either participate in or catalyze subsequent reactions. This model is meant to mimic the sel…
▽ More
We introduce an autocatalytic aggregation model in which the rate at which two clusters merge to form a cluster is controlled by the presence of a third "catalytic" cluster whose mass must equal to the mass of one of the reaction partners. The catalyst is unaffected by the joining event and is available to either participate in or catalyze subsequent reactions. This model is meant to mimic the self-replicating reactions that occur in models for the origin of life. We solve the kinetics of this catalytic coagulation model for the case of mass-independent rates and show that the total cluster density decays as $t^{-1/3}$, while the density of clusters of any fixed mass decays as $t^{-2/3}$. These behaviors contrast with the corresponding $t^{-1}$ and $t^{-2}$ scalings for classic aggregation. We extend our model to mass-dependent reaction rates, to situations where only "magic" mass clusters can catalyze reactions, and to include steady monomer input.
△ Less
Submitted 25 June, 2024; v1 submitted 25 April, 2024;
originally announced April 2024.
-
Simple evolving random graphs
Authors:
P. L. Krapivsky
Abstract:
We study the evolution of graphs densifying by adding edges: Two vertices are chosen randomly, and an edge is (i) established if each vertex belongs to a tree; (ii) established with probability $p$ if only one vertex belongs to a tree; (iii) an attempt fails if both vertices belong to unicyclic components. Emerging simple random graphs contain only trees and unicycles. In the thermodynamic limit o…
▽ More
We study the evolution of graphs densifying by adding edges: Two vertices are chosen randomly, and an edge is (i) established if each vertex belongs to a tree; (ii) established with probability $p$ if only one vertex belongs to a tree; (iii) an attempt fails if both vertices belong to unicyclic components. Emerging simple random graphs contain only trees and unicycles. In the thermodynamic limit of an infinite number of vertices, the fraction of vertices in unicycles undergoes a phase transition resembling a percolation transition in classical random graphs. In classical random graphs, a complex giant component born at the transition eventually engulfs all finite components and densifies forever. The evolution of simple random graphs freezes when trees disappear. We quantify simple random graphs in the supercritical phase and the properties of the frozen state.
△ Less
Submitted 25 January, 2024; v1 submitted 5 December, 2023;
originally announced December 2023.
-
Epidemic Forecast Follies
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We introduce a simple multiplicative model to describe the temporal behavior and the ultimate outcome of an epidemic. Our model accounts, in a minimalist way, for the competing influences of imposing public-health restrictions when the epidemic is severe, and relaxing restrictions when the epidemic is waning. Our primary results are that different instances of an epidemic with identical starting p…
▽ More
We introduce a simple multiplicative model to describe the temporal behavior and the ultimate outcome of an epidemic. Our model accounts, in a minimalist way, for the competing influences of imposing public-health restrictions when the epidemic is severe, and relaxing restrictions when the epidemic is waning. Our primary results are that different instances of an epidemic with identical starting points have disparate outcomes and each epidemic temporal history is strongly fluctuating.
△ Less
Submitted 23 September, 2023;
originally announced September 2023.
-
Random Maps with Sociological Flavor
Authors:
P. L. Krapivsky
Abstract:
A map of a set to itself admits a representation by a graph with vertices being the elements of the set and an edge between every vertex and its image. Communities defined as the maximal connected components are uni-cyclic. The distributions of the sizes of communities and lengths of cycles for unconstrained random maps is a classical subject. We call experts the images and followers the remaining…
▽ More
A map of a set to itself admits a representation by a graph with vertices being the elements of the set and an edge between every vertex and its image. Communities defined as the maximal connected components are uni-cyclic. The distributions of the sizes of communities and lengths of cycles for unconstrained random maps is a classical subject. We call experts the images and followers the remaining vertexes, and we further define prophets, egocentrics, and introverts. We introduce and analyze classes of random maps with sociological flavor.
△ Less
Submitted 25 January, 2024; v1 submitted 15 September, 2023;
originally announced September 2023.
-
Steady oscillations in aggregation-fragmentation processes
Authors:
N. V. Brilliantov,
W. Otieno,
S. A. Matveev,
A. P. Smirnov,
E. E. Tyrtyshnikov,
P. L. Krapivsky
Abstract:
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)…
▽ More
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 simulations results when $ θ\equiv ν-μ<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 two 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.
△ Less
Submitted 17 July, 2023;
originally announced July 2023.
-
The Magic of Networks Grown by Redirection
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We highlight intriguing features of complex networks that are grown by \emph{redirection}. In this mechanism, a target node is chosen uniformly at random from the pre-existing network nodes and the new node attaches either to this initial target or to a neighbor of this target. This exceedingly simple algorithm generates preferential attachment networks in an algorithmic time that is linear in the…
▽ More
We highlight intriguing features of complex networks that are grown by \emph{redirection}. In this mechanism, a target node is chosen uniformly at random from the pre-existing network nodes and the new node attaches either to this initial target or to a neighbor of this target. This exceedingly simple algorithm generates preferential attachment networks in an algorithmic time that is linear in the number of network nodes $N$. Even though preferential attachment ostensibly requires \emph{global knowledge} of the network, redirection requires only \emph{local knowledge}. We also show that changing just a \emph{single} attachment rate in linear preferential attachment leads to a non-universal degree distribution. Finally, we present unexpected consequences of redirection in networks with undirected links, where highly modular and non-sparse networks arise.
△ Less
Submitted 16 August, 2023; v1 submitted 17 May, 2023;
originally announced May 2023.
-
Blast waves in the zero temperature hard sphere gas: double scaling structure
Authors:
Sahil Kumar Singh,
Subhadip Chakraborti,
Abhishek Dhar,
P. L. Krapivsky
Abstract:
We study the blast generated by sudden localized release of energy in a cold gas. Specifically, we consider one-dimensional hard-rod gas and two-dimensional hard disc gas. For this problem, the Taylor-von Neumann-Sedov (TvNS) solution of Euler equations has a self-similar form. The shock wave remains infinitely strong for the zero-temperature gas, so the solution applies indefinitely. The TvNS sol…
▽ More
We study the blast generated by sudden localized release of energy in a cold gas. Specifically, we consider one-dimensional hard-rod gas and two-dimensional hard disc gas. For this problem, the Taylor-von Neumann-Sedov (TvNS) solution of Euler equations has a self-similar form. The shock wave remains infinitely strong for the zero-temperature gas, so the solution applies indefinitely. The TvNS solution ignores dissipation, however. We show that this is erroneous in the core region which, in two dimensions, expands as $t^{2/5}$ while the shock wave propagates as $t^{1/2}$. A new self-similar solution depending on the scaling variable $r/t^{2/5}$ describes the core, while the TvNS solution describes the bulk. We demonstrate this from a numerical solution of the Navier-Stokes (NS) equations and from molecular dynamics simulations for a gas of hard discs in two dimensions and hard rods in one dimension. In both cases, the shock front position predicted by NS equations and by the TvNS solution agrees with that predicted by molecular dynamics simulations. However, the NS equations fail to describe the near-core form of the scaling functions.
△ Less
Submitted 12 July, 2023; v1 submitted 18 November, 2022;
originally announced November 2022.
-
Birds on a Wire
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We investigate the occupancy statistics of birds on a wire and on higher-dimensional substrates. In one dimension, birds land one by one on a wire and rest where they land. Whenever a newly arriving bird lands within a fixed distance of already resting birds, these resting birds immediately fly away. We determine the steady-state occupancy of the wire, the distribution of gaps between neighboring…
▽ More
We investigate the occupancy statistics of birds on a wire and on higher-dimensional substrates. In one dimension, birds land one by one on a wire and rest where they land. Whenever a newly arriving bird lands within a fixed distance of already resting birds, these resting birds immediately fly away. We determine the steady-state occupancy of the wire, the distribution of gaps between neighboring birds, and other basic statistical features of this process. We discuss conjectures for corresponding observables in higher dimensions.
△ Less
Submitted 26 July, 2022; v1 submitted 2 May, 2022;
originally announced May 2022.
-
A splash in a one-dimensional cold gas
Authors:
Subhadip Chakraborti,
Abhishek Dhar,
P. L. Krapivsky
Abstract:
We consider a set of hard point particles distributed uniformly with a specified density on the positive half-line and all initially at rest. The particle masses alternate between two values, $m$ and $M$. The particles interact via collisions that conserve energy and momentum. We study the cascade of activity that results when the left-most particle is given a positive velocity. At long times we f…
▽ More
We consider a set of hard point particles distributed uniformly with a specified density on the positive half-line and all initially at rest. The particle masses alternate between two values, $m$ and $M$. The particles interact via collisions that conserve energy and momentum. We study the cascade of activity that results when the left-most particle is given a positive velocity. At long times we find that this leads to two fascinating features in the observed dynamics. First, in the bulk of the gas, a shock front develops separating the cold gas from a thermalized region. The shock-front travels sub-ballistically, with the bulk described by self-similar solutions of Euler hydrodynamics. Second, there is a splash region formed by the recoiled particles which move ballistically with negative velocities. The splash region is completely non-hydrodynamic and we propose two conjectures for the long time particle dynamics in this region. We provide a detailed analytic understanding of these coexisting regimes. These are supported by the results of molecular dynamics simulations.
△ Less
Submitted 13 July, 2022; v1 submitted 21 March, 2022;
originally announced March 2022.
-
Diffusion in a fluid flow generated by a source at the apex of a wedge
Authors:
P. L. Krapivsky
Abstract:
We consider a particle diffusing inside a wedge with absorbing boundaries and driven by a radial flow of incompressible fluid generated by a source at the apex. The survival probability decays as (time)^{-b} with exponent depending on the opening angle of the wedge and the Reynolds number associated with the hydrodynamic flow. The computation of the decay exponent reduces to finding the ground sta…
▽ More
We consider a particle diffusing inside a wedge with absorbing boundaries and driven by a radial flow of incompressible fluid generated by a source at the apex. The survival probability decays as (time)^{-b} with exponent depending on the opening angle of the wedge and the Reynolds number associated with the hydrodynamic flow. The computation of the decay exponent reduces to finding the ground state energy of the quantum particle in an infinitely deep potential well with shape determined by the radial flow velocity.
△ Less
Submitted 1 April, 2022; v1 submitted 30 December, 2021;
originally announced December 2021.
-
Aggregation in non-uniform systems with advection and localized source
Authors:
Rishat R. Zagidullin,
Alexander P. Smirnov,
Sergey A. Matveev,
Nikolai V. Brilliantov,
P. L. Krapivsky
Abstract:
We explore analytically and numerically agglomeration driven by advection and localized source. The system is inhomogeneous in one dimension, viz. along the direction of advection. We analyze a simplified model with mass-independent advection velocity, diffusion coefficient, and reaction rates. We also examine a model with mass-dependent coefficients describing aggregation with sedimentation. For…
▽ More
We explore analytically and numerically agglomeration driven by advection and localized source. The system is inhomogeneous in one dimension, viz. along the direction of advection. We analyze a simplified model with mass-independent advection velocity, diffusion coefficient, and reaction rates. We also examine a model with mass-dependent coefficients describing aggregation with sedimentation. For the simplified model, we obtain an exact solution for the stationary spatially dependent agglomerate densities. In the model describing aggregation with sedimentation, we report a new conservation law and develop a scaling theory for the densities. For numerical efficiency we exploit the low-rank approximation technique; this dramatically increases the computational speed and allows simulations of large systems. The numerical results are in excellent agreement with the predictions of our theory.
△ Less
Submitted 28 December, 2021; v1 submitted 25 December, 2021;
originally announced December 2021.
-
Non-extensive super-cluster states in aggregation with fragmentation
Authors:
Nikolai V. Brilliantov,
Wendy Otieno,
P. L. Krapivsky
Abstract:
Systems evolving through aggregation and fragmentation may possess an intriguing super-cluster state (SCS). Clusters constituting this state are mostly large, so the SCS resembles a gelling state. The formation of the SCS is controlled by fluctuations and in this aspect, it is similar to a critical state. The SCS is non-extensive as the number of clusters varies sub-linearly with the system size.…
▽ More
Systems evolving through aggregation and fragmentation may possess an intriguing super-cluster state (SCS). Clusters constituting this state are mostly large, so the SCS resembles a gelling state. The formation of the SCS is controlled by fluctuations and in this aspect, it is similar to a critical state. The SCS is non-extensive as the number of clusters varies sub-linearly with the system size. In the parameter space, the SCS separates equilibrium and jamming (extensive) states. The conventional methods such as 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.
△ Less
Submitted 1 November, 2021;
originally announced November 2021.
-
Divergence and Consensus in Majority Rule
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We investigate majority rule dynamics in a population with two classes of people, each with two opinion states $\pm 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 probability $ε$. Consensus is achieved in a time that scales…
▽ More
We investigate majority rule dynamics in a population with two classes of people, each with two opinion states $\pm 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 probability $ε$. Consensus is achieved in a time that scales logarithmically with population size if $ε\geq ε_c=\frac{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.
△ Less
Submitted 6 March, 2021;
originally announced March 2021.
-
Blast in a One-Dimensional Cold Gas: From Newtonian Dynamics to Hydrodynamics
Authors:
Subhadip Chakraborti,
Santhosh Ganapa,
P. L. Krapivsky,
Abhishek Dhar
Abstract:
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 expecte…
▽ More
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.
△ Less
Submitted 21 June, 2021; v1 submitted 16 February, 2021;
originally announced February 2021.
-
The Taylor-von Neumann-Sedov blast-wave solution: comparisons with microscopic simulations of a one-dimensional gas
Authors:
Santhosh Ganapa,
Subhadip Chakraborti,
P. L. Krapivsky,
Abhishek Dhar
Abstract:
We study the response of an infinite system of point particles on the line initially at rest on the instantaneous release of energy in a localized region. We make a detailed comparison of the hydrodynamic variables predicted by Euler equations for non-dissipative ideal compressible gas and the results of direct microscopic simulations. At long times the profiles of the three conserved variables ev…
▽ More
We study the response of an infinite system of point particles on the line initially at rest on the instantaneous release of energy in a localized region. We make a detailed comparison of the hydrodynamic variables predicted by Euler equations for non-dissipative ideal compressible gas and the results of direct microscopic simulations. At long times the profiles of the three conserved variables evolve to self-similar scaling forms, with a scaling exponent as predicted by the Taylor-von Neumann-Sedov (TvNS) blast-wave solution. The scaling functions obtained from the microscopic dynamics show a remarkable agreement with the TvNS predictions, except at the blast core, where the TvNS solution predicts a diverging temperature which is not observed in simulations. We show that the effect of heat conduction becomes important and present results from a numerical solution of the full Navier-Stokes-Fourier equations. A different scaling form is observed in the blast core and this is carefully analyzed. Our microscopic model is the one-dimensional alternate mass hard-particle gas which has the ideal gas equation of state but is non-integrable and known to display fast equilibration.
△ Less
Submitted 17 January, 2021; v1 submitted 29 October, 2020;
originally announced October 2020.
-
The role of energy in ballistic agglomeration
Authors:
N. V. Brilliantov,
A. I. Osinsky,
P. L. Krapivsky
Abstract:
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 allo…
▽ More
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.
△ Less
Submitted 2 October, 2020;
originally announced October 2020.
-
Infinite System of Random Walkers: Winners and Losers
Authors:
P. L. Krapivsky
Abstract:
We study an infinite system of particles initially occupying a half-line $y\leq 0$ and undergoing random walks on the entire line. The right-most particle is called a leader. Surprisingly, every particle except the original leader may never achieve the leadership throughout the evolution. For the equidistant initial configuration, the $k^{\text{th}}$ particle attains the leadership with probabilit…
▽ More
We study an infinite system of particles initially occupying a half-line $y\leq 0$ and undergoing random walks on the entire line. The right-most particle is called a leader. Surprisingly, every particle except the original leader may never achieve the leadership throughout the evolution. For the equidistant initial configuration, the $k^{\text{th}}$ particle attains the leadership with probability $e^{-2} k^{-1} (\ln k)^{-1/2}$ when $k\gg 1$. This provides a quantitative measure of the correlation between earlier misfortune (represented by $k$) and eternal failure. We also show that the winner defined as the first walker overtaking the initial leader has label $k\gg 1$ with probability decaying as $\exp\!\left[-\tfrac{1}{2}(\ln k)^2\right]$.
△ Less
Submitted 28 September, 2020;
originally announced September 2020.
-
An infection process near criticality: Influence of the initial condition
Authors:
P. L. Krapivsky
Abstract:
We investigate how the initial number of infected individuals affects the behavior of the critical susceptible-infected-recovered process. We analyze the outbreak size distribution, duration of the outbreaks, and the role of fluctuations.
We investigate how the initial number of infected individuals affects the behavior of the critical susceptible-infected-recovered process. We analyze the outbreak size distribution, duration of the outbreaks, and the role of fluctuations.
△ Less
Submitted 3 December, 2020; v1 submitted 18 September, 2020;
originally announced September 2020.
-
Epidemics with containment measures
Authors:
Ginestra Bianconi,
Pavel L. Krapivsky
Abstract:
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, mimic…
▽ More
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.
△ Less
Submitted 28 October, 2020; v1 submitted 8 April, 2020;
originally announced April 2020.
-
Where Should You Park Your Car? The $\frac{1}{2}$ Rule
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We investigate parking in a one-dimensional lot, where cars enter at a rate $λ$ and each attempts to park close to a target at the origin. Parked cars also depart at rate 1. An entering driver cannot see beyond the parked cars for more desirable open spots. We analyze a class of strategies in which a driver ignores open spots beyond $τL$, where $τ$ is a risk threshold and $L$ is the location of th…
▽ More
We investigate parking in a one-dimensional lot, where cars enter at a rate $λ$ and each attempts to park close to a target at the origin. Parked cars also depart at rate 1. An entering driver cannot see beyond the parked cars for more desirable open spots. We analyze a class of strategies in which a driver ignores open spots beyond $τL$, where $τ$ is a risk threshold and $L$ is the location of the most distant parked car, and attempts to park at the first available spot encountered closer than $τL$. When all drivers use this strategy, the probability to park at the best available spot is maximal when $τ=\frac{1}{2}$, and parking at the best available spot occurs with probability $\frac{1}{4}$.
△ Less
Submitted 9 June, 2020; v1 submitted 23 March, 2020;
originally announced March 2020.
-
Simple Parking Strategies
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We investigate simple strategies that embody the decisions that one faces when trying to park near a popular destination. Should one park far from the target (destination), where finding a spot is easy, but then be faced with a long walk, or should one attempt to look for a desirable spot close to the target, where spots may be hard to find? We study an idealized parking process on a one-dimension…
▽ More
We investigate simple strategies that embody the decisions that one faces when trying to park near a popular destination. Should one park far from the target (destination), where finding a spot is easy, but then be faced with a long walk, or should one attempt to look for a desirable spot close to the target, where spots may be hard to find? We study an idealized parking process on a one-dimensional geometry where the desired target is located at $x=0$, cars enter the system from the right at a rate $λ$ and each car leaves at a unit rate. We analyze three parking strategies---meek, prudent, and optimistic---and determine which is optimal.
△ Less
Submitted 23 February, 2020; v1 submitted 13 April, 2019;
originally announced April 2019.
-
Optimally Frugal Foraging
Authors:
O. Benichou,
U. Bhat,
P. L. Krapivsky,
S. Redner
Abstract:
We introduce the \emph{frugal foraging} model in which a forager performs a discrete-time random walk on a lattice, where each site initially contains $\mathcal{S}$ food units. The forager metabolizes one unit of food at each step and starves to death when it last ate $\mathcal{S}$ steps in the past. Whenever the forager decides to eat, it consumes all food at its current site and this site remain…
▽ More
We introduce the \emph{frugal foraging} model in which a forager performs a discrete-time random walk on a lattice, where each site initially contains $\mathcal{S}$ food units. The forager metabolizes one unit of food at each step and starves to death when it last ate $\mathcal{S}$ steps in the past. Whenever the forager decides to eat, it consumes all food at its current site and this site remains empty (no food replenishment). The crucial property of the forager is that it is \emph{frugal} and eats only when encountering food within at most $k$ steps of starvation. We compute the average lifetime analytically as a function of frugality threshold and show that there exists an optimal strategy, namely, a frugality threshold $k^*$ that maximizes the forager lifetime.
△ Less
Submitted 9 November, 2017;
originally announced November 2017.
-
Trapping and Escape in a Turbid Medium
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We investigate the absorption of diffusing molecules in a fluid-filled spherical beaker that contains many small reactive traps. The molecules are absorbed either by hitting a trap or by escaping via the beaker walls. In the physical situation where the number $N$ of traps is large and their radii $a$ are small compared to the beaker radius $R$, the fraction of molecules $E$ that escape to the bea…
▽ More
We investigate the absorption of diffusing molecules in a fluid-filled spherical beaker that contains many small reactive traps. The molecules are absorbed either by hitting a trap or by escaping via the beaker walls. In the physical situation where the number $N$ of traps is large and their radii $a$ are small compared to the beaker radius $R$, the fraction of molecules $E$ that escape to the beaker wall and the complementary fraction $T$ that eventually are absorbed by the traps depend only on the dimensionless parameter combination $λ= Na/R$. We compute $E$ and $T$ as a function of $λ$ for a spherical beaker and for beakers of other three-dimensional shapes. The asymptotic behavior is found to be universal: $1- E\sim λ$ for $λ\to 0$ and $E\simλ^{-1/2}$ for $λ\to\infty$.
△ Less
Submitted 13 October, 2017;
originally announced October 2017.
-
Stability of a Giant Connected Component in a Complex Network
Authors:
Maksim Kitsak,
Alexander A. Ganin,
Daniel A. Eisenberg,
Pavel L. Krapivsky,
Dmitri Krioukov,
David L. Alderson,
Igor Linkov
Abstract:
We analyze the stability of the network's giant connected component under impact of adverse events, which we model through the link percolation. Specifically, we quantify the extent to which the largest connected component of a network consists of the same nodes, regardless of the specific set of deactivated links. Our results are intuitive in the case of single-layered systems: the presence of la…
▽ More
We analyze the stability of the network's giant connected component under impact of adverse events, which we model through the link percolation. Specifically, we quantify the extent to which the largest connected component of a network consists of the same nodes, regardless of the specific set of deactivated links. Our results are intuitive in the case of single-layered systems: the presence of large degree nodes in a single-layered network ensures both its robustness and stability. In contrast, we find that interdependent networks that are robust to adverse events have unstable connected components. Our results bring novel insights to the design of resilient network topologies and the reinforcement of existing networked systems.
△ Less
Submitted 23 January, 2018; v1 submitted 28 September, 2017;
originally announced September 2017.
-
Oscillations in aggregation-shattering processes
Authors:
S. A. Matveev,
P. L. Krapivsky,
A. P. Smirnov,
E. E. Tyrtyshnikov,
N. V. Brilliantov
Abstract:
We observe never-ending oscillations in systems undergoing aggregation and collision-controlled 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 syst…
▽ More
We observe never-ending oscillations in systems undergoing aggregation and collision-controlled 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 to 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 the shattering rate 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.
△ Less
Submitted 4 August, 2017;
originally announced August 2017.
-
Emergent Network Modularity
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We introduce a network growth model based on complete redirection: a new node randomly selects an existing target node, but attaches to a random neighbor of this target. For undirected networks, this simple growth rule generates unusual, highly modular networks. Individual network realizations typically contain multiple macrohubs---nodes whose degree scales linearly with the number of nodes $N$. T…
▽ More
We introduce a network growth model based on complete redirection: a new node randomly selects an existing target node, but attaches to a random neighbor of this target. For undirected networks, this simple growth rule generates unusual, highly modular networks. Individual network realizations typically contain multiple macrohubs---nodes whose degree scales linearly with the number of nodes $N$. The size of the network "nucleus"---the set of nodes of degree greater than one---grows sublinearly with $N$ and thus constitutes a vanishingly small fraction of the network. The network therefore consists almost entirely of leaves (nodes of degree one) as $N\to\infty$.
△ Less
Submitted 20 June, 2017; v1 submitted 5 June, 2017;
originally announced June 2017.
-
Densification and Structural Transitions in Networks that Grow by Node Copying
Authors:
U. Bhat,
P. L. Krapivsky,
R. Lambiotte,
S. Redner
Abstract:
We introduce a growing network model---the copying model---in which a new node attaches to a randomly selected target node and, in addition, independently to each of the neighbors of the target with copying probability $p$. When $p<\frac{1}{2}$, this algorithm generates sparse networks, in which the average node degree is finite. A power-law degree distribution also arises, with a non-universal ex…
▽ More
We introduce a growing network model---the copying model---in which a new node attaches to a randomly selected target node and, in addition, independently to each of the neighbors of the target with copying probability $p$. When $p<\frac{1}{2}$, this algorithm generates sparse networks, in which the average node degree is finite. A power-law degree distribution also arises, with a non-universal exponent whose value is determined by a transcendental equation in $p$. In the sparse regime, the network is "normal", e.g., the relative fluctuations in the number of links are asymptotically negligible. For $p\geq \frac{1}{2}$, the emergent networks are dense (the average degree increases with the number of nodes $N$) and they exhibit intriguing structural behaviors. In particular, the $N$-dependence of the number of $m$-cliques (complete subgraphs of $m$ nodes) undergoes $m-1$ transitions from normal to progressively more anomalous behavior at a $m$-dependent critical values of $p$. Different realizations of the network, which start from the same initial state, exhibit macroscopic fluctuations in the thermodynamic limit---absence of self averaging. When linking to second neighbors of the target node can occur, the number of links asymptotically grows as $N^2$ as $N\to\infty$, so that the network is effectively complete as $N\to \infty$.
△ Less
Submitted 5 October, 2016;
originally announced October 2016.
-
Structural Transitions in Dense Networks
Authors:
R. Lambiotte,
P. L. Krapivsky,
U. Bhat,
S. Redner
Abstract:
We introduce an evolving network model in which a new node attaches to a randomly selected target node and also to each of its neighbors with probability $p$. The resulting network is sparse for $p<\frac{1}{2}$ and dense (average degree increasing with number of nodes $N$) for $p\geq \frac{1}{2}$. In the dense regime, individual networks realizations built by this copying mechanism are disparate a…
▽ More
We introduce an evolving network model in which a new node attaches to a randomly selected target node and also to each of its neighbors with probability $p$. The resulting network is sparse for $p<\frac{1}{2}$ and dense (average degree increasing with number of nodes $N$) for $p\geq \frac{1}{2}$. In the dense regime, individual networks realizations built by this copying mechanism are disparate and not self-averaging. Further, there is an infinite sequence of structural anomalies at $p=\frac{2}{3}$, $\frac{3}{4}$, $\frac{4}{5}$, etc., where the dependences on $N$ of the number of triangles (3-cliques), 4-cliques, undergo phase transitions. When linking to second neighbors of the target can occur, the probability that the resulting graph is complete---where all nodes are connected---is non-zero as $N\to\infty$.
△ Less
Submitted 13 July, 2016;
originally announced July 2016.
-
Emergence of Clustering in an Acquaintance Model without Homophily
Authors:
Uttam Bhat,
P. L. Krapivsky,
S. Redner
Abstract:
We introduce an agent-based acquaintance model in which social links are created by processes in which there is no explicit homophily. In spite of the homogeneous nature of the social interactions, highly-clustered social networks can arise. The crucial feature of our model is that of variable transitive interactions. Namely, when an agent introduces two unconnected friends, the rate at which a co…
▽ More
We introduce an agent-based acquaintance model in which social links are created by processes in which there is no explicit homophily. In spite of the homogeneous nature of the social interactions, highly-clustered social networks can arise. The crucial feature of our model is that of variable transitive interactions. Namely, when an agent introduces two unconnected friends, the rate at which a connection actually occurs between them depends on the number of their mutual acquaintances. As this transitive interaction rate is varied, the social network undergoes a dramatic clustering transition. Close to the transition, the network consists of a collection of well-defined communities. As a function of time, the network can also undergo an \emph{incomplete} gelation transition, in which the gel, or giant cluster, does not constitute the entire network, even at infinite time. Some of the clustering properties of our model also arise, but in a more gradual manner, in Facebook networks. Finally, we discuss a more realistic variant of our original model in which there is a soft cutoff in the rate of transitive interactions. With this variant, one can construct network realizations that quantitatively match Facebook networks.
△ Less
Submitted 12 November, 2014; v1 submitted 27 August, 2014;
originally announced August 2014.
-
Large Deviations in Single File Diffusion
Authors:
P. L. Krapivsky,
Kirone Mallick,
Tridib Sadhu
Abstract:
We apply macroscopic fluctuation theory to study the diffusion of a tracer in a one-dimensional interacting particle system with excluded mutual passage, known as single-file diffusion. In the case of Brownian point particles with hard-core repulsion, we derive the cumulant generating function of the tracer position and its large deviation function. In the general case of arbitrary inter-particle…
▽ More
We apply macroscopic fluctuation theory to study the diffusion of a tracer in a one-dimensional interacting particle system with excluded mutual passage, known as single-file diffusion. In the case of Brownian point particles with hard-core repulsion, we derive the cumulant generating function of the tracer position and its large deviation function. In the general case of arbitrary inter-particle interactions, we express the variance of the tracer position in terms of the collective transport properties, viz. the diffusion coefficient and the mobility. Our analysis applies both for fluctuating (annealed) and fixed (quenched) initial configurations.
△ Less
Submitted 28 July, 2014; v1 submitted 5 May, 2014;
originally announced May 2014.
-
Highly Dispersed Networks Generated by Enhanced Redirection
Authors:
Alan Gabel,
P. L. Krapivsky,
S. Redner
Abstract:
We analyze growing networks that are built by enhanced redirection. Nodes are sequentially added and each incoming node attaches to a randomly chosen 'target' node with probability 1-r, or to the parent of the target node with probability r. When the redirection probability r is an increasing function of the degree of the parent node, with r-->1 as the parent degree diverges, networks grown via th…
▽ More
We analyze growing networks that are built by enhanced redirection. Nodes are sequentially added and each incoming node attaches to a randomly chosen 'target' node with probability 1-r, or to the parent of the target node with probability r. When the redirection probability r is an increasing function of the degree of the parent node, with r-->1 as the parent degree diverges, networks grown via this enhanced redirection mechanism exhibit unusual properties, including: (i) multiple macrohubs---nodes with degrees proportional to the number of network nodes N; (ii) non-extensivity of the degree distribution in which the number of nodes of degree k, N_k, scales as N^{nu-1}/k^{nu}, with 1<nu<2; (iii) lack of self-averaging, with large fluctuations between individual network realizations. These features are robust and continue to hold when the incoming node has out-degree greater than 1 so that networks contain closed loops. The latter networks are strongly clustered; for the specific case of the double attachment, the average local clustering coefficient is <C_i>=4(ln2)-2=0.77258...
△ Less
Submitted 13 March, 2014; v1 submitted 30 December, 2013;
originally announced December 2013.
-
Choice-Driven Phase Transition in Complex Networks
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We investigate choice-driven network growth. In this model, nodes are added one by one according to the following procedure: for each addition event a set of target nodes is selected, each according to linear preferential attachment, and a new node attaches to the target with the highest degree. Depending on precise details of the attachment rule, the resulting networks has three possible outcomes…
▽ More
We investigate choice-driven network growth. In this model, nodes are added one by one according to the following procedure: for each addition event a set of target nodes is selected, each according to linear preferential attachment, and a new node attaches to the target with the highest degree. Depending on precise details of the attachment rule, the resulting networks has three possible outcomes: (i) a non-universal power-law degree distribution; (ii) a single macroscopic hub (a node whose degree is of the order of N, the number of network nodes), while the remainder of the nodes comprises a non-universal power-law degree distribution; (iii) a degree distribution that decays as (k ln k)^{-2} at the transition between cases (i) and (ii). These properties are robust when attachment occurs to the highest-degree node from at least two targets. When attachment is made to a target whose degree is not the highest, the degree distribution has the ultra-narrow double-exponential form exp(-const. x e^k), from which the largest degree grows only as ln(ln N).
△ Less
Submitted 24 March, 2014; v1 submitted 30 December, 2013;
originally announced December 2013.
-
Highly Dispersed Networks
Authors:
Alan Gabel,
P. L. Krapivsky,
S. Redner
Abstract:
We introduce a new class of networks that grow by enhanced redirection. Nodes are introduced sequentially, and each either attaches to a randomly chosen target node with probability 1-r or to the ancestor of the target with probability r, where r an increasing function of the degree of the ancestor. This mechanism leads to highly-dispersed networks with unusual properties: (i) existence of multipl…
▽ More
We introduce a new class of networks that grow by enhanced redirection. Nodes are introduced sequentially, and each either attaches to a randomly chosen target node with probability 1-r or to the ancestor of the target with probability r, where r an increasing function of the degree of the ancestor. This mechanism leads to highly-dispersed networks with unusual properties: (i) existence of multiple macrohubs---nodes whose degree is a finite fraction of the total number of network nodes N, (ii) lack of self averaging, and (iii) anomalous scaling, in which N_k, the number of nodes of degree k scales as N_k N^{nu-1}/k^{nu}, with 1<nu<2.
△ Less
Submitted 14 July, 2013;
originally announced July 2013.
-
Distinct Degrees and Their Distribution in Complex Networks
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We investigate a variety of statistical properties associated with the number of distinct degrees that exist in a typical network for various classes of networks. For a single realization of a network with N nodes that is drawn from an ensemble in which the number of nodes of degree k has an algebraic tail, N_k ~ N/k^nu for k>>1, the number of distinct degrees grows as N^{1/nu}. Such an algebraic…
▽ More
We investigate a variety of statistical properties associated with the number of distinct degrees that exist in a typical network for various classes of networks. For a single realization of a network with N nodes that is drawn from an ensemble in which the number of nodes of degree k has an algebraic tail, N_k ~ N/k^nu for k>>1, the number of distinct degrees grows as N^{1/nu}. Such an algebraic growth is also observed in scientific citation data. We also determine the N dependence of statistical quantities associated with the sparse, large-k range of the degree distribution, such as the location of the first hole (where N_k=0), the last doublet (two consecutive occupied degrees), triplet, dimer (N_k=2), trimer, etc.
△ Less
Submitted 23 June, 2013; v1 submitted 6 April, 2013;
originally announced April 2013.
-
Molecular Spiders on the Plane
Authors:
Tibor Antal,
P. L. Krapivsky
Abstract:
Synthetic bio-molecular spiders with "legs" made of single-stranded segments of DNA can move on a surface covered by single-stranded segments of DNA called substrates when the substrate DNA is complementary to the leg DNA. If the motion of a spider does not affect the substrates, the spider behaves asymptotically as a random walk. We study the diffusion coefficient and the number of visited sites…
▽ More
Synthetic bio-molecular spiders with "legs" made of single-stranded segments of DNA can move on a surface covered by single-stranded segments of DNA called substrates when the substrate DNA is complementary to the leg DNA. If the motion of a spider does not affect the substrates, the spider behaves asymptotically as a random walk. We study the diffusion coefficient and the number of visited sites for spiders moving on the square lattice with a substrate in each lattice site. The spider's legs hop to nearest-neighbor sites with the constraint that the distance between any two legs cannot exceed a maximal span. We establish analytic results for bipedal spiders, and investigate multileg spiders numerically. In experimental realizations legs usually convert substrates into products (visited sites). The binding of legs to products is weaker, so the hopping rate from the substrates is smaller. This makes the problem non-Markovian and we investigate it numerically. We demonstrate the emergence of a counter-intuitive behavior - the more spiders are slowed down on unvisited sites, the more motile they become.
△ Less
Submitted 20 July, 2012; v1 submitted 28 April, 2012;
originally announced April 2012.
-
Popularity-Driven Networking
Authors:
E. Ben-Naim,
P. L. Krapivsky
Abstract:
We investigate the growth of connectivity in a network. In our model, starting with a set of disjoint nodes, links are added sequentially. Each link connects two nodes, and the connection rate governing this random process is proportional to the degrees of the two nodes. Interestingly, this network exhibits two abrupt transitions, both occurring at finite times. The first is a percolation transiti…
▽ More
We investigate the growth of connectivity in a network. In our model, starting with a set of disjoint nodes, links are added sequentially. Each link connects two nodes, and the connection rate governing this random process is proportional to the degrees of the two nodes. Interestingly, this network exhibits two abrupt transitions, both occurring at finite times. The first is a percolation transition in which a giant component, containing a finite fraction of all nodes, is born. The second is a condensation transition in which the entire system condenses into a single, fully connected, component. We derive the size distribution of connected components as well as the degree distribution, which is purely exponential throughout the evolution. Furthermore, we present a criterion for the emergence of sudden condensation for general homogeneous connection rates.
△ Less
Submitted 30 November, 2011;
originally announced December 2011.
-
Exact solution of a two-type branching process: Models of tumor progression
Authors:
Tibor Antal,
P. L. Krapivsky
Abstract:
An explicit solution for a general two-type birth-death branching process with one way mutation is presented. This continuous time process mimics the evolution of resistance to treatment, or the onset of an extra driver mutation during tumor progression. We obtain the exact generating function of the process at arbitrary times, and derive various large time scaling limits. In the simultaneous smal…
▽ More
An explicit solution for a general two-type birth-death branching process with one way mutation is presented. This continuous time process mimics the evolution of resistance to treatment, or the onset of an extra driver mutation during tumor progression. We obtain the exact generating function of the process at arbitrary times, and derive various large time scaling limits. In the simultaneous small mutation rate and large time scaling limit, the distribution of the mutant cells develops some atypical properties, including a power law tail and diverging average.
△ Less
Submitted 5 May, 2011;
originally announced May 2011.
-
Reinforcement-Driven Spread of Innovations and Fads
Authors:
P. L. Krapivsky,
S. Redner,
D. Volovik
Abstract:
We propose kinetic models for the spread of permanent innovations and transient fads by the mechanism of social reinforcement. Each individual can be in one of M+1 states of awareness 0,1,2,...,M, with state M corresponding to adopting an innovation. An individual with awareness k<M increases to k+1 by interacting with an adopter. Starting with a single adopter, the time for an initially unaware p…
▽ More
We propose kinetic models for the spread of permanent innovations and transient fads by the mechanism of social reinforcement. Each individual can be in one of M+1 states of awareness 0,1,2,...,M, with state M corresponding to adopting an innovation. An individual with awareness k<M increases to k+1 by interacting with an adopter. Starting with a single adopter, the time for an initially unaware population of size N to adopt a permanent innovation grows as ln(N) for M=1, and as N^{1-1/M} for M>1. The fraction of the population that remains clueless about a transient fad after it has come and gone changes discontinuously as a function of the fad abandonment rate lambda for M>1. The fad dies out completely in a time that varies non-monotonically with lambda.
△ Less
Submitted 29 September, 2011; v1 submitted 20 April, 2011;
originally announced April 2011.
-
Wealth Distributions in Asset Exchange Models
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
How do individuals accumulate wealth as they interact economically? We outline the consequences of a simple microscopic model in which repeated pairwise exchanges of assets between individuals build the wealth distribution of a population. This distribution is determined for generic exchange rules --- transactions that involve a fixed amount or a fixed fraction of individual wealth, as well as ran…
▽ More
How do individuals accumulate wealth as they interact economically? We outline the consequences of a simple microscopic model in which repeated pairwise exchanges of assets between individuals build the wealth distribution of a population. This distribution is determined for generic exchange rules --- transactions that involve a fixed amount or a fixed fraction of individual wealth, as well as random or greedy exchanges. In greedy multiplicative exchange, a continuously evolving power law wealth distribution arises, a feature that qualitatively mimics empirical observations.
△ Less
Submitted 29 August, 2010; v1 submitted 23 June, 2010;
originally announced June 2010.
-
Random Ancestor Trees
Authors:
E. Ben-Naim,
P. L. Krapivsky
Abstract:
We investigate a network growth model in which the genealogy controls the evolution. In this model, a new node selects a random target node and links either to this target node, or to its parent, or to its grandparent, etc; all nodes from the target node to its most ancient ancestor are equiprobable destinations. The emerging random ancestor tree is very shallow: the fraction g_n of nodes at dist…
▽ More
We investigate a network growth model in which the genealogy controls the evolution. In this model, a new node selects a random target node and links either to this target node, or to its parent, or to its grandparent, etc; all nodes from the target node to its most ancient ancestor are equiprobable destinations. The emerging random ancestor tree is very shallow: the fraction g_n of nodes at distance n from the root decreases super-exponentially with n, g_n=e^{-1}/(n-1)!. We find that a macroscopic hub at the root coexists with highly connected nodes at higher generations. The maximal degree of a node at the nth generation grows algebraically as N^{1/beta_n} where N is the system size. We obtain the series of nontrivial exponents which are roots of transcendental equations: beta_1= 1.351746, beta_2=1.682201, etc. As a consequence, the fraction p_k of nodes with degree k has algebraic tail, p_k ~ k^{-gamma}, with gamma=beta_1+1=2.351746.
△ Less
Submitted 10 April, 2010;
originally announced April 2010.
-
Shepherd Model for Knot-Limited Polymer Ejection from a Capsid
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner
Abstract:
We construct a tractable model to describe the rate at which a knotted polymer is ejected from a spherical capsid via a small pore. Knots are too large to fit through the pore and must reptate to the end of the polymer for ejection to occur. The reptation of knots is described by symmetric exclusion on the line, with the internal capsid pressure represented by an additional biased particle that…
▽ More
We construct a tractable model to describe the rate at which a knotted polymer is ejected from a spherical capsid via a small pore. Knots are too large to fit through the pore and must reptate to the end of the polymer for ejection to occur. The reptation of knots is described by symmetric exclusion on the line, with the internal capsid pressure represented by an additional biased particle that drives knots to the end of the chain. We compute the exact ejection speed for a finite number of knots L and find that it scales as 1/L. We also construct a continuum theory for many knots that matches the exact discrete theory for large L.
△ Less
Submitted 11 August, 2009; v1 submitted 6 May, 2009;
originally announced May 2009.
-
Exciting Hard Spheres
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner
Abstract:
We investigate the collision cascade that is generated by a single moving incident particle on a static hard-sphere gas. We argue that the number of moving particles at time t grows as t^{xi} and the number collisions up to time t grows as t^{eta}, with xi=2d/(d+2) and eta=2(d+1)/(d+2) and d the spatial dimension. These growth laws are the same as those from a hydrodynamic theory for the shock w…
▽ More
We investigate the collision cascade that is generated by a single moving incident particle on a static hard-sphere gas. We argue that the number of moving particles at time t grows as t^{xi} and the number collisions up to time t grows as t^{eta}, with xi=2d/(d+2) and eta=2(d+1)/(d+2) and d the spatial dimension. These growth laws are the same as those from a hydrodynamic theory for the shock wave emanating from an explosion. Our predictions are verified by molecular dynamics simulations in d=1 and 2. For a particle incident on a static gas in a half-space, the resulting backsplatter ultimately contains almost all the initial energy.
△ Less
Submitted 25 May, 2008;
originally announced May 2008.
-
Transition from small to large world in growing networks
Authors:
S. N. Dorogovtsev,
P. L. Krapivsky,
J. F. F. Mendes
Abstract:
We examine the global organization of growing networks in which a new vertex is attached to already existing ones with a probability depending on their age. We find that the network is infinite- or finite-dimensional depending on whether the attachment probability decays slower or faster than $(age)^{-1}$. The network becomes one-dimensional when the attachment probability decays faster than…
▽ More
We examine the global organization of growing networks in which a new vertex is attached to already existing ones with a probability depending on their age. We find that the network is infinite- or finite-dimensional depending on whether the attachment probability decays slower or faster than $(age)^{-1}$. The network becomes one-dimensional when the attachment probability decays faster than $(age)^{-2}$. We describe structural characteristics of these phases and transitions between them.
△ Less
Submitted 23 November, 2007; v1 submitted 19 September, 2007;
originally announced September 2007.
-
Dynamics of Microtubule Instabilities
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner
Abstract:
We investigate the dynamics of an idealized model of microtubule growth that evolves by: (i) attachment of guanosine triphosphate (GTP) at rate lambda, (ii) conversion of GTP to guanosine diphosphate (GDP) at rate 1, and (iii) detachment of GDP at rate mu. As a function of these rates, a microtubule can grow steadily or its length can fluctuate wildly. For mu=0, we find the exact tubule and GTP ca…
▽ More
We investigate the dynamics of an idealized model of microtubule growth that evolves by: (i) attachment of guanosine triphosphate (GTP) at rate lambda, (ii) conversion of GTP to guanosine diphosphate (GDP) at rate 1, and (iii) detachment of GDP at rate mu. As a function of these rates, a microtubule can grow steadily or its length can fluctuate wildly. For mu=0, we find the exact tubule and GTP cap length distributions, and power-law length distributions of GTP and GDP islands. For mu=infinity, we argue that the time between catastrophes, where the microtubule shrinks to zero length, scales as exp(lambda). We also find the phase boundary between a growing and shrinking microtubule.
△ Less
Submitted 23 November, 2010; v1 submitted 14 March, 2007;
originally announced March 2007.
-
Dynamics of Microtubule Growth and Catastrophe
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner,
M. Mailman,
B. Chakraborty
Abstract:
We investigate a simple model of microtubule dynamics in which a microtubule evolves by: (i) attachment of guanosine triphosphate (GTP) to its end at rate lambda, (ii) GTP converting irreversibly to guanosine diphosphate (GDP) at rate 1, and (iii) detachment of GDP from the end of a microtubule at rate mu. As a function of these elemental rates, the microtubule can grow steadily or its length ca…
▽ More
We investigate a simple model of microtubule dynamics in which a microtubule evolves by: (i) attachment of guanosine triphosphate (GTP) to its end at rate lambda, (ii) GTP converting irreversibly to guanosine diphosphate (GDP) at rate 1, and (iii) detachment of GDP from the end of a microtubule at rate mu. As a function of these elemental rates, the microtubule can grow steadily or its length can fluctuate wildly. A master equation approach is developed to characterize these intriguing features. For mu=0, we find exact expressions for tubule and GTP cap length distributions, as well as a power-law length distributions of GTP and GDP islands. For mu=oo, we find the average time between catastrophes, where the microtubule shrinks to zero length, and extend this approach to also determine the size distribution of avalanches (sequence of consecutive GDP detachment events). We obtain the phase diagram for general rates and verify our predictions by numerical simulations.
△ Less
Submitted 22 April, 2008; v1 submitted 1 March, 2007;
originally announced March 2007.
-
Probing Non-Integer Dimensions
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We show that two-dimensional convection-diffusion problems with a radial sink or source at the origin may be recast as a pure diffusion problem in a fictitious space in which the spatial dimension is continuously-tunable with the Peclet number. This formulation allows us to probe various diffusion-controlled processes in non-integer dimensions.
We show that two-dimensional convection-diffusion problems with a radial sink or source at the origin may be recast as a pure diffusion problem in a fictitious space in which the spatial dimension is continuously-tunable with the Peclet number. This formulation allows us to probe various diffusion-controlled processes in non-integer dimensions.
△ Less
Submitted 18 September, 2006;
originally announced September 2006.
-
Social Balance on Networks: The Dynamics of Friendship and Enmity
Authors:
T. Antal,
P. L. Krapivsky,
S. Redner
Abstract:
How do social networks evolve when both friendly and unfriendly relations exist? Here we propose a simple dynamics for social networks in which the sense of a relationship can change so as to eliminate imbalanced triads--relationship triangles that contains 1 or 3 unfriendly links. In this dynamics, a friendly link changes to unfriendly or vice versa in an imbalanced triad to make the triad bala…
▽ More
How do social networks evolve when both friendly and unfriendly relations exist? Here we propose a simple dynamics for social networks in which the sense of a relationship can change so as to eliminate imbalanced triads--relationship triangles that contains 1 or 3 unfriendly links. In this dynamics, a friendly link changes to unfriendly or vice versa in an imbalanced triad to make the triad balanced. Such networks undergo a dynamic phase transition from a steady state to "utopia"--all friendly links--as the amount of network friendliness is changed. Basic features of the long-time dynamics and the phase transition are discussed.
△ Less
Submitted 21 May, 2006;
originally announced May 2006.
-
Random Walk with Shrinking Steps
Authors:
P. L. Krapivsky,
S. Redner
Abstract:
We outline basic properties of a symmetric random walk in one dimension, in which the length of the nth step equals lambda^n, with lambda<1. As the number of steps N-->oo, the probability that the endpoint is at x, P_{lambda}(x;N), approaches a limiting distribution P_{lambda}(x) that has many beautiful features. For lambda<1/2, the support of P_{lambda}(x) is a Cantor set. For 1/2<=lambda<1, th…
▽ More
We outline basic properties of a symmetric random walk in one dimension, in which the length of the nth step equals lambda^n, with lambda<1. As the number of steps N-->oo, the probability that the endpoint is at x, P_{lambda}(x;N), approaches a limiting distribution P_{lambda}(x) that has many beautiful features. For lambda<1/2, the support of P_{lambda}(x) is a Cantor set. For 1/2<=lambda<1, there is a countably infinite set of lambda values for which P_{lambda}(x) is singular, while P_{lambda}(x) is smooth for almost all other lambda values. In the most interesting case of lambda=(sqrt{5}-1)/2=g, P_g(x) is riddled with singularities and is strikingly self-similar. The self-similarity is exploited to derive a simple form for the probability measure M(a,b)= int_a^b P_g(x) dx.
△ Less
Submitted 18 October, 2003; v1 submitted 9 April, 2003;
originally announced April 2003.