Researcher profile

E. Ben-Naim

E. Ben-Naim contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
15topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2014arXiv

Fragmentation of Random Trees

We study fragmentation of a random recursive tree into a forest by repeated removal of nodes. The initial tree consists of N nodes and it is generated by sequential addition of nodes with each new node attaching to a randomly-selected existing node. As nodes are removed from the tree, one at a time, the tree dissolves into an ensemble of separate trees, namely, a forest. We study statistical properties of trees and nodes in this heterogeneous forest, and find that the fraction of remaining nodes m characterizes the system in the limit N --> infty. We obtain analytically the size density phi_s of trees of size s. The size density has power-law tail phi_s ~ s^(-alpha) with exponent alpha=1+1/m. Therefore, the tail becomes steeper as further nodes are removed, and the fragmentation process is unusual in that exponent alpha increases continuously with time. We also extend our analysis to the case where nodes are added as well as removed, and obtain the asymptotic size density for growing trees.

preprint2014arXiv

Persistence of Random Walk Records

We study records generated by Brownian particles in one dimension. Specifically, we investigate an ordinary random walk and define the record as the maximal position of the walk. We compare the record of an individual random walk with the mean record, obtained as an average over infinitely many realizations. We term the walk "superior" if the record is always above average, and conversely, the walk is said to be "inferior" if the record is always below average. We find that the fraction of superior walks, S, decays algebraically with time, S ~ t^(-beta), in the limit t --> infty, and that the persistence exponent is nontrivial, beta=0.382258.... The fraction of inferior walks, I, also decays as a power law, I ~ t^(-alpha), but the persistence exponent is smaller, alpha=0.241608.... Both exponents are roots of transcendental equations involving the parabolic cylinder function. To obtain these theoretical results, we analyze the joint density of superior walks with given record and position, while for inferior walks it suffices to study the density as function of position.

preprint2014arXiv

Slow Kinetics of Brownian Maxima

We study extreme-value statistics of Brownian trajectories in one dimension. We define the maximum as the largest position to date and compare maxima of two particles undergoing independent Brownian motion. We focus on the probability P(t) that the two maxima remain ordered up to time t, and find the algebraic decay P ~ t^(-beta) with exponent beta=1/4. When the two particles have diffusion constants D1 and D2, the exponent depends on the mobilities, beta=(1/pi)arctan[sqrt(D2/D1)]. We also use numerical simulations to investigate maxima of multiple particles in one dimension and the largest extension of particles in higher dimensions.

preprint2013arXiv

First Passage in Conical Geometry and Ordering of Brownian Particles

We survey recent results on first-passage processes in unbounded cones and their applications to ordering of particles undergoing Brownian motion in one dimension. We first discuss the survival probability S(t) that a diffusing particle, in arbitrary spatial dimension, remains inside a conical domain up to time t. In general, this quantity decays algebraically S ~ t^{-beta} in the long-time limit. The exponent beta depends on the opening angle of the cone and the spatial dimension, and it is root of a transcendental equation involving the associated Legendre functions. The exponent becomes a function of a single scaling variable in the limit of large spatial dimension. We then describe two first-passage problems involving the order of N independent Brownian particles in one dimension where survival probabilities decay algebraically as well. To analyze these problems, we identify the trajectories of the N particles with the trajectory of one particle in N dimensions, confined to within a certain boundary, and we use a circular cone with matching solid angle as a replacement for the confining boundary. For N=3, the confining boundary is a wedge and the approach is exact. In general, this "cone approximation" gives strict lower bounds as well as useful estimates or the first-passage exponents. Interestingly, the cone approximation becomes asymptotically exact when N-->infinity as it predicts the exact scaling function that governs the spectrum of first-passage exponents.

preprint2013arXiv

Recurrence Statistics of Great Earthquakes

We investigate the sequence of great earthquakes over the past century. To examine whether the earthquake record includes temporal clustering, we identify aftershocks and remove those from the record. We focus on the recurrence time, defined as the time between two consecutive earthquakes. We study the variance in the recurrence time and the maximal recurrence time. Using these quantities, we compare the earthquake record with sequences of random events, generated by numerical simulations, while systematically varying the minimal earthquake magnitude Mmin. Our analysis shows that the earthquake record is consistent with a random process for magnitude thresholds 7.0<=Mmin<=8.3, where the number of events is larger. Interestingly, the earthquake record deviates from a random process at magnitude threshold 8.4<=Mmin<= 8.5, where the number of events is smaller; however, this deviation is not strong enough to conclude that great earthquakes are clustered. Overall, the findings are robust both qualitatively and quantitatively as statistics of extreme values and moment analysis yield remarkably similar results.

preprint2013arXiv

Scaling Exponent for Incremental Records

We investigate records in a growing sequence of identical and independently distributed random variables. The record equals the largest value in the sequence, and our focus is on the increment, defined as the difference between two successive records. We investigate sequences in which all increments decrease monotonically, and find that the fraction I_N of sequences that exhibit this property decays algebraically with sequence length N, namely I_N ~ N^{-nu} as N --> infinity. We analyze the case where the random variables are drawn from a uniform distribution with compact support, and obtain the exponent nu = 0.317621... using analytic methods. We also study the record distribution and the increment distribution. Whereas the former is a narrow distribution with an exponential tail, the latter is broad and has a power-law tail characterized by the exponent nu. Empirical analysis of records in the sequence of waiting times between successive earthquakes is consistent with the theoretical results.

preprint2013arXiv

Statistics of Superior Records

We study statistics of records in a sequence of random variables. These identical and independently distributed variables are drawn from the parent distribution rho. The running record equals the maximum of all elements in the sequence up to a given point. We define a superior sequence as one where all running records are above the average record, expected for the parent distribution rho. We find that the fraction of superior sequences S_N decays algebraically with sequence length N, S_N ~ N^{-beta} in the limit N-->infty. Interestingly, the decay exponent beta is nontrivial, being the root of an integral equation. For example, when rho is a uniform distribution with compact support, we find beta=0.450265. In general, the tail of the parent distribution governs the exponent beta. We also consider the dual problem of inferior sequences, where all records are below average, and find that the fraction of inferior sequences I_N decays algebraically, albeit with a different decay exponent, I_N ~ N^{-alpha}. We use the above statistical measures to analyze earthquake data.

preprint2012arXiv

Discrete Analog of the Burgers Equation

We propose the set of coupled ordinary differential equations dn_j/dt=(n_{j-1})^2-(n_j)^2 as a discrete analog of the classic Burgers equation. We focus on traveling waves and triangular waves, and find that these special solutions of the discrete system capture major features of their continuous counterpart. In particular, the propagation velocity of a traveling wave and the shape of a triangular wave match the continuous behavior. However, there are some subtle differences. For traveling waves, the propagating front can be extremely sharp as it exhibits double exponential decay. For triangular waves, there is an unexpected logarithmic shift in the location of the front. We establish these results using asymptotic analysis, heuristic arguments, and direct numerical integration.

preprint2012arXiv

Randomness in Competitions

We study the effects of randomness on competitions based on an elementary random process in which there is a finite probability that a weaker team upsets a stronger team. We apply this model to sports leagues and sports tournaments, and compare the theoretical results with empirical data. Our model shows that single-elimination tournaments are efficient but unfair: the number of games is proportional to the number of teams N, but the probability that the weakest team wins decays only algebraically with N. In contrast, leagues, where every team plays every other team, are fair but inefficient: the top $\sqrt{N}$ of teams remain in contention for the championship, while the probability that the weakest team becomes champion is exponentially small. We also propose a gradual elimination schedule that consists of a preliminary round and a championship round. Initially, teams play a small number of preliminary games, and subsequently, a few teams qualify for the championship round. This algorithm is fair and efficient: the best team wins with a high probability and the number of games scales as $N^{9/5}$, whereas traditional leagues require N^3 games to fairly determine a champion.

preprint2012arXiv

Scaling Behavior of Threshold Epidemics

We study the classic Susceptible-Infected-Recovered (SIR) model for the spread of an infectious disease. In this stochastic process, there are two competing mechanism: infection and recovery. Susceptible individuals may contract the disease from infected individuals, while infected ones recover from the disease at a constant rate and are never infected again. Our focus is the behavior at the epidemic threshold where the rates of the infection and recovery processes balance. In the infinite population limit, we establish analytically scaling rules for the time-dependent distribution functions that characterize the sizes of the infected and the recovered sub-populations. Using heuristic arguments, we also obtain scaling laws for the size and duration of the epidemic outbreaks as a function of the total population. We perform numerical simulations to verify the scaling predictions and discuss the consequences of these scaling laws for near-threshold epidemic outbreaks.

preprint2011arXiv

Dynamics of Random Graphs with Bounded Degrees

We investigate the dynamic formation of regular random graphs. In our model, we pick a pair of nodes at random and connect them with a link if both of their degrees are smaller than d. Starting with a set of isolated nodes, we repeat this linking step until a regular random graph, where all nodes have degree d, forms. We view this process as a multivariate aggregation process, and formally solve the evolution equations using the Hamilton-Jacoby formalism. We calculate the nontrivial percolation thresholds for the emergence of the giant component when d>=3. Also, we estimate the number of steps until the giant component spans the entire system and the total number of steps until the regular random graph forms. These quantities are non self-averaging, namely, they fluctuate from realization to realization even in the thermodynamic limit.

preprint2011arXiv

Kinetics of Ring Formation

We study reversible polymerization of rings. In this stochastic process, two monomers bond and as a consequence, two disjoint rings may merge into a compound ring, or, a single ring may split into two fragment rings. This aggregation-fragmentation process exhibits a percolation transition with a finite-ring phase in which all rings have microscopic length and a giant-ring phase where macroscopic rings account for a finite fraction of the entire mass. Interestingly, while the total mass of the giant rings is a deterministic quantity, their total number and their sizes are stochastic quantities. The size distribution of the macroscopic rings is universal, although the span of this distribution increases with time. Moreover, the average number of giant rings scales logarithmically with system size. We introduce a card-shuffling algorithm for efficient simulation of the ring formation process, and present numerical verification of the theoretical predictions.

preprint2011arXiv

Popularity-Driven Networking

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.

preprint2010arXiv

First Passage Properties of the Polya Urn Process

We study first passage statistics of the Polya urn model. In this random process, the urn contains two types of balls. In each step, one ball is drawn randomly from the urn, and subsequently placed back into the urn together with an additional ball of the same type. We derive the probability G_n that the two types of balls are equal in number, for the first time, when there is a total of 2n balls. This first passage probability decays algebraically, G_n ~ n^{-2}, when n is large. We also derive the probability that a tie ever happens. This probability is between zero and one, so that a tie may occur in some realizations but not in others. The likelihood of a tie is appreciable only if the initial difference in the number balls is of the order of the square-root of the total number of balls.

preprint2010arXiv

First-Passage Exponents of Multiple Random Walks

We investigate first-passage statistics of an ensemble of N noninteracting random walks on a line. Starting from a configuration in which all particles are located in the positive half-line, we study S_n(t), the probability that the nth rightmost particle remains in the positive half-line up to time t. This quantity decays algebraically, S_n (t) ~ t^{-beta_n}, in the long-time limit. Interestingly, there is a family of nontrivial first-passage exponents, beta_1<beta_2<...<beta_{N-1}; the only exception is the two-particle case where beta_1=1/3. In the N-->infinity limit, however, the exponents attain a scaling form, beta_n(N)--> beta(z) with z=(n-N/2)/sqrt{N}. We also demonstrate that the smallest exponent decays exponentially with N. We deduce these results from first-passage kinetics of a random walk in an N-dimensional cone and confirm them using numerical simulations. Additionally, we investigate the family of exponents that characterizes leadership statistics of multiple random walks and find that in this case, the cone provides an excellent approximation.

preprint2010arXiv

Granular gases under extreme driving

We study inelastic gases in two dimensions using event-driven molecular dynamics simulations. Our focus is the nature of the stationary state attained by rare injection of large amounts of energy to balance the dissipation due to collisions. We find that under such extreme driving, with the injection rate much smaller than the collision rate, the velocity distribution has a power-law high energy tail. The numerically measured exponent characterizing this tail is in excellent agreement with predictions of kinetic theory over a wide range of system parameters. We conclude that driving by rare but powerful energy injection leads to a well-mixed gas and constitutes an alternative mechanism for agitating granular matter. In this distinct nonequilibrium steady-state, energy cascades from large to small scales. Our simulations also show that when the injection rate is comparable with the collision rate, the velocity distribution has a stretched exponential tail.

preprint2010arXiv

Kinetics of First Passage in a Cone

We study statistics of first passage inside a cone in arbitrary spatial dimension. The probability that a diffusing particle avoids the cone boundary decays algebraically with time. The decay exponent depends on two variables: the opening angle of the cone and the spatial dimension. In four dimensions, we find an explicit expression for the exponent, and in general, we obtain it as a root of a transcendental equation involving associated Legendre functions. At large dimensions, the decay exponent depends on a single scaling variable, while roots of the parabolic cylinder function specify the scaling function. Consequently, the exponent is of order one only if the cone surface is very close to a plane. We also perform asymptotic analysis for extremely thin and extremely wide cones.

preprint2010arXiv

On the Mixing of Diffusing Particles

We study how the order of N independent random walks in one dimension evolves with time. Our focus is statistical properties of the inversion number m, defined as the number of pairs that are out of sort with respect to the initial configuration. In the steady-state, the distribution of the inversion number is Gaussian with the average <m>~N^2/4 and the standard deviation sigma N^{3/2}/6. The survival probability, S_m(t), which measures the likelihood that the inversion number remains below m until time t, decays algebraically in the long-time limit, S_m t^{-beta_m}. Interestingly, there is a spectrum of N(N-1)/2 distinct exponents beta_m(N). We also find that the kinetics of first-passage in a circular cone provides a good approximation for these exponents. When N is large, the first-passage exponents are a universal function of a single scaling variable, beta_m(N)--> beta(z) with z=(m-<m>)/sigma. In the cone approximation, the scaling function is a root of a transcendental equation involving the parabolic cylinder equation, D_{2 beta}(-z)=0, and surprisingly, numerical simulations show this prediction to be exact.

preprint2007arXiv

Singular Energy Distributions in Granular Media

We study the kinetic theory of driven granular gases, taking into account both translational and rotational degrees of freedom. We obtain the high-energy tail of the stationary bivariate energy distribution, depending on the total energy E and the ratio x=sqrt{E_w/E} of rotational energy E_w to total energy. Extremely energetic particles have a unique and well-defined distribution f(x) which has several remarkable features: x is not uniformly distributed as in molecular gases; f(x) is not smooth but has multiple singularities. The latter behavior is sensitive to material properties such as the collision parameters, the moment of inertia and the collision rate. Interestingly, there are preferred ratios of rotational-to-total energy. In general, f(x) is strongly correlated with energy and the deviations from a uniform distribution grow with energy. We also solve for the energy distribution of freely cooling Maxwell Molecules and find qualitatively similar behavior.