Source author record

Nadav M. Shnerb

Nadav M. Shnerb appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

26 published item(s)

preprint2022arXiv

The time-identity tradeoff

Distinguishability plays a major role in quantum and statistical physics. When particles are identical their wave function must be either symmetric or antisymmetric under permutations and the number of microscopic states, which determines entropy, is counted up to permutations. When the particles are distinguishable, wavefunctions have no symmetry and each permutation is a different microstate. This binary and discontinuous classification raises a few questions: one may wonder what happens if particles are almost identical, or when the property that distinguishes between them is irrelevant to the physical interactions in a given system. Here I sketch a general answer to these questions. For any pair of non-identical particles there is a timescale, $τ_d$, required for a measurement to resolve the differences between them. Below $τ_d$, particles seem identical, above it - different, and the uncertainty principle provides a lower bound for $τ_d$. Thermal systems admit a conjugate temperature scale, $T_d$. Above this temperature the system appears to equilibrate before it resolves the differences between particles, below this temperature the system identifies these differences before equilibration. As the physical differences between particles decline towards zero, $τ_d \to \infty$ and $T_d \to 0$.

preprint2020arXiv

Population dynamics in stochastic environments

Populations are made up of an integer number of individuals and are subject to stochastic birth-death processes whose rates may vary in time. Useful quantities, like the chance of ultimate fixation, satisfy an appropriate difference (master) equation, but closed-form solutions of these equations are rare. Analytical insights in fields like population genetics, ecology and evolution rely, almost exclusively, on an uncontrolled application of the diffusion approximation (DA) which assumes the smoothness of the relevant quantities over the set of integers. Here we combine asymptotic matching techniques with a first-order (controlling-factor) WKB method to obtain a theory whose range of applicability is much wider. This allows us to rederive DA from a more general theory, to identify its limitations, and to suggest alternative analytical solutions and scalable numerical techniques when it fails. We carry out our analysis for the calculation of the fixation probability in a fluctuating environment, highlighting the difference between (on average) deleterious and beneficial mutant invasion and the intricate distinction between weak and strong selection.

preprint2020arXiv

Stochasticity-induced stabilization in ecology and evolution: a new synthesis

The ability of random environmental variation to stabilize competitor coexistence was pointed out long ago and, in recent years, has received considerable attention. Analyses have focused on variations in the log-abundances of species, with mean logarithmic growth rates when rare, $\mathbb{E}[r]$, used as metrics for persistence. However, invasion probabilities and the times to extinction are not single-valued functions of $\mathbb{E}[r]$ and, in some cases, decrease as $\mathbb{E}[r]$ increases. Here, we present a synthesis of stochasticity-induced stabilization (SIS) phenomena based on the ratio between the expected arithmetic growth $μ$ and its variance $g$. When the diffusion approximation holds, explicit formulas for invasion probabilities and persistence times are single valued, monotonic functions of $μ/g$. The storage effect in the lottery model, together with other well-known examples drawn from population genetics, microbiology and ecology (including discrete and continuous dynamics, with overlapping and non-overlapping generations), are placed together, reviewed, and explained within this new, transparent theoretical framework. We also clarify the relationships between life-history strategies and SIS, and study the dynamics of extinction when SIS fails.

preprint2020arXiv

Taming the diffusion approximation through a controlling-factor WKB method

The diffusion approximation (DA) is widely used in the analysis of stochastic population dynamics, from population genetics to ecology and evolution. DA is an uncontrolled approximation that assumes the smoothness of the calculated quantity over the relevant state space and fails when this property is not satisfied. This failure becomes severe in situations where the direction of selection switches sign. Here we employ the WKB (large-deviations) method, which requires only the logarithm of a given quantity to be smooth over its state space. Combining the WKB scheme with asymptotic matching techniques, we show how to derive the diffusion approximation in a controlled manner and how to produce better approximations, applicable for much wider regimes of parameters. We also introduce a scalable (independent of population size) WKB-based numerical technique. The method is applied to a central problem in population genetics and evolution, finding the chance of ultimate fixation in a zero-sum, two-types competition.

preprint2016arXiv

Communities as cliques

High-diversity assemblages are very common in nature, and yet the factors allowing for the maintenance of biodiversity remain obscure. The competitive exclusion principle and May's complexity-diversity puzzle both suggest that a community can support only a small number of species, turning the spotlight at the dynamics of local patches or islands, where stable and uninvadable (SU) subsets of species play a crucial role. Here we map the community SUs question to the geometric problem of finding maximal cliques of the corresponding graph. We solve for the number of SUs as a function of the species richness in the regional pool, $N$, showing that this growth is subexponential, contrary to long-standing wisdom. We show that symmetric systems relax rapidly to an SU, where the system stays until a regime shift takes place. In asymmetric systems the relaxation time grows much faster with $N$, suggesting an excitable dynamics under noise.

preprint2016arXiv

Spatial neutral dynamics

Neutral models, in which individual agents with equal fitness undergo a birth-death-mutation process, are very popular in population genetics and community ecology. Usually these models are applied to populations and communities with spatial structure, but the analytic results presented so far are limited to well-mixed or mainland-island scenarios. Here we present a new technique, based on interface dynamics analysis, and apply it to the neutral dynamics in one, two and three spatial dimensions. New results are derived for the correlation length and for the main characteristics of the community, like total biodiversity and the species abundance distribution above the correlation length. Our results are supported by extensive numerical simulations, and provide qualitative and quantitative insights that allow for a rigorous comparison between model predictions and empirical data.

preprint2016arXiv

Stability of two-species communities: drift, environmental stochasticity, storage effect and selection

The dynamics of two competing species in a finite size community is one of the most studied problems in population genetics and community ecology. Stochastic fluctuations lead, inevitably, to the extinction of one of the species, but the relevant timescale depends on the underlying dynamics. The persistence time of the community has been calculated for neutral models, where the only drive of the system is drift (demographic stochasticity) and for models with strong selection. Following recent analyses that stress the importance of environmental stochasticity in empirical systems, we present here a general theory of persistence time of two-species community where drift, environmental variations and time independent selective advantage are all taken into account.

preprint2016arXiv

The effect of environmental stochasticity on species richness in neutral communities

Environmental stochasticity is known to be a destabilizing factor, increasing abundance fluctuations and extinction rates of populations. However, the stability of a community may benefit from the differential response of species to environmental variations due to the storage effect. This paper provides a systematic and comprehensive discussion of these two contradicting tendencies, using the metacommunity version of the recently proposed time-average neutral model of biodiversity which incorporates environmental stochasticity and demographic noise and allows for extinction and speciation. We show that the incorporation of demographic noise into the model is essential to its applicability, yielding realistic behavior of the system when fitness variations are relatively weak. The dependence of species richness on the strength of environmental stochasticity changes sign when the correlation time of the environmental variations increases. This transition marks the point at which the storage effect no longer succeeds in stabilizing the community.

preprint2015arXiv

Age-abundance relationships for neutral communities

Neutral models for the dynamics of a system of competing species are used, nowadays, to describe a wide variety of empirical communities. These models are used in many situations, ranging from population genetics and ecological biodiversity to macroevolution and cancer tumors. One of the main issues discussed within this framework is the relationships between the abundance of a species and its age. Here we provide a comprehensive analysis of the age-abundance relationships for fixed-size and growing communities. Explicit formulas for the average and the most likely age of a species with abundance $n$ are given, together with the full probability distribution function. We further discuss the universality of these results and their applicability to the tropical forest community.

preprint2015arXiv

Fronts and fluctuations at a critical surface

The properties of a front between two different phases in the presence of a smoothly inhomogeneous external field that takes its critical value at the crossing point is analyzed. Two generic scenarios are studied. In the first, the system admits a bistable solution and the external field governs the rate in which one phase invades the other. The second mechanism corresponds to a second order transition that, in the case of reactive systems, takes the form of a transcritical bifurcation at the crossing point. We solve for the front shape and its response to external white noise, showing that static properties and also some of the dynamics features cannot distinguish between the two scenarios. The only reliable indicator turns out to be the fluctuation statistics. These take a Gaussian form in the bifurcation case and a double-peak shape in a bistable system. The results of a recent analysis of the morphogenesis process in Drosophila embryos are reanalyzed and we show, in contrast with the interpretation suggested by Krotov et. al., that the plausible underlying dynamics is bistable and not bifurcational.

preprint2015arXiv

Neutral dynamics with environmental noise: age-size statistics and species lifetimes

Neutral dynamics, where taxa are assumed to be demographically equivalent and their abundance is governed solely by the stochasticity of the underlying birth-death process, has proved itself as an important minimal model that accounts for many empirical datasets in genetics and ecology. However, the restriction of the model to demographic [${\cal{O}} ({\sqrt N})$] noise yields relatively slow dynamics that appears to be in conflict with both short-term and long-term characteristics of the observed systems. Here we analyze two of these problems - age size relationships and species extinction time - in the framework of a neutral theory with both demographic and environmental stochasticity. It turns out that environmentally induced variations of the demographic rates control the long-term dynamics and modify dramatically the predictions of the neutral theory with demographic noise only, yielding much better agreement with empirical data. We consider two prototypes of "zero mean" environmental noise, one which is balanced with regard to the arithmetic abundance, another balanced in the logarithmic (fitness) space, study their species lifetime statistics and discuss their relevance to realistic models of community dynamics.

preprint2015arXiv

Predicting catastrophic shifts

Catastrophic transitions, where a system shifts abruptly between alternate steady states, are a generic feature of many nonlinear systems. Recently these regime shift were suggested as the mechanism underlies many ecological catastrophes, such as desertification and coral reef collapses, which are considered as a prominent threat to sustainability and to the well-being of millions. Still, the methods proposed so far for the prediction of an imminent transition are quite ineffective, and some empirical and theoretical studies suggest that actual transitions may occur smoothly, without an abrupt shift. Here we present a new diagnostic tool, based on monitoring the dynamics of clusters through time. Our technique discriminates between systems with local positive feedback, where the transition is abrupt, and systems with negative density dependence, where the transition is smooth. Analyzing the spatial dynamics of these two generic scenarios, we show that changes in the critical cluster size provide a reliable early warning indicator for both transitions. Our method may allow for the prediction, and thus hopefully the prevention of such transitions, avoiding their destructive outcomes.

preprint2014arXiv

A generalized model of island biodiversity

The dynamics of a local community of competing species with weak immigration from a static regional pool is studied. Implementing the generalized competitive Lotka-Volterra model with demographic noise, a rich dynamics structure with four qualitatively distinct phases is unfolded. When the overall interspecies competition is weak, the island species are a sample of the mainland species. For higher values of the competition parameter the system still admit an equilibrium community, but now some of the mainland species are absent on the island. Further increase in competition leads to an intermittent "chaotic" phase, where the dynamics is controlled by invadable combinations of species and the turnover rate is governed by the migration. Finally, the strong competition phase is glassy, dominated by uninvadable state and noise-induced transitions. Our model contains, as a spatial case, the celebrated neutral island theories of Wilson-MacArthur and Hubbell. Moreover, we show that slight deviations from perfect neutrality may lead to each of the phases, as the Hubbell point appears to be quadracritical.

preprint2014arXiv

A generalized neutral theory explains static and dynamic properties of biotic communities

Understanding the forces shaping ecological communities is crucially important to basic science and conservation. In recent years, considerable progress was made in explaining communities using simple and general models, with neutral theory as a prominent example. However, while successful in explaining static patterns such as species abundance distributions, the neutral theory was criticized for making unrealistic predictions of fundamental dynamic patterns. Here we incorporate environmental stochasticity into the neutral framework, and show that the resulting generalized neutral theory is capable of predicting realistic patterns of both population and community dynamics. Applying the theory to real data (the tropical forest of Barro-Colorado Island), we find that it better fits the observed distribution of short-term fluctuations, the temporal scaling of such fluctuations, and the decay of compositional similarity with time, than the original theory, while retaining its power to explain static patterns of species abundance. Importantly, although the proposed theory is neutral (all species are functionally equivalent) and stochastic, it is a niche-based theory in the sense that species differ in their demographic responses to environmental variation. Our results show that this integration of niche forces and stochasticity within a minimalistic neutral framework is highly successful in explaining fundamental static and dynamic characteristics of ecological communities.

preprint2014arXiv

Spatial patterns in the tropical forest reveal connections between negative feedback, aggregation and abundance

The spatial arrangement of trees in a tropical forest reflects the interplay between aggregating processes, like dispersal limitation, and negative feedback that induces effective repulsion among individuals. Monitoring the variance-mean ratio for conspecific individuals along length-scales, we show that the effect of negative feedback is dominant at short scales, while aggregation characterizes the large-scale patterns. A comparison of different species indicates, surprisingly, that both aggregation and negative feedback scales are related to the overall abundance of the species. This suggests a bottom-up control mechanism, in which the negative feedback dictates the dispersal kernel and the overall abundance.

preprint2014arXiv

The glocal forest

Ecological spatial patterns reflect the underlying processes that shape the structure of species and communities. Mechanisms like inter and intra species competition, dispersal and host-pathogen interactions are believed to act over a wide range of scales, and the inference of the process from the pattern is, despite its popularity, a challenging task. Here we call attention to a quite unexpected phenomenon in the extensively studied tropical forest at the Barro-Colorado Island (BCI): the spatial deployment of (almost) all tree species is statistically equivalent, once distances are normalized by , the typical distance between neighboring conspecific trees. Correlation function, cluster statistics and nearest-neighbor distance distribution become species-independent after this rescaling. Global observables (species frequencies) and local spatial structure appear to be interrelated. This "glocality" suggests a radical interpretation of recent experiments that show a correlation between species' abundance and the negative feedback among conspecifics. For the forest to be glocal, the negative feedback must govern spatial patterns over all scales.

preprint2013arXiv

Evolutionarily Stable Density-Dependent Dispersal

An ab-initio numerical study of the density-dependent, evolutionary stable dispersal strategy is presented. The simulations are based on a simple discretei generation island model with four processes: reproduction, dispersal, competition and local catastrophe. We do not impose any a priori constraints on the dispersal schedule, allowing the entire schedule to evolve. We find that the system converges at long times to a unique nontrivial dispersal schedule such that the dispersal probability is a monotonically increasing function of the density. We have explored the dependence of the selected dispersal strategy on the various system parameters: mean number of offspring, site carrying capacity, dispersal cost and system size. A few general scaling laws are seen to emerge from the data.

preprint2013arXiv

Stochastic Desertification

The process of desertification is usually modeled as a first order transition, where a change of an external parameter (e.g. precipitation) leads to a catastrophic bifurcation followed by an ecological regime shift. However, vegetation elements like shrubs and trees undergo a stochastic birth-death process with an absorbing state; such a process supports a second order continuous transition with no hysteresis. We present a numerical study of a minimal model that supports bistability and catastrophic shift on spatial domain with demographic noise and an absorbing state. When the external parameter varies adiabatically the transition is continuous and the front velocity renormalizes to zero at the extinction transition. Below the transition one may identify three modes of desertification: accumulation of local catastrophes, desert invasion and global collapse. A catastrophic regime shift occurs as a dynamical hysteresis, when the pace of environmental variations is too fast. We present some empirical evidence, suggesting that the mid-holocene desertification of the Sahara was, indeed, continuous.

preprint2012arXiv

A scaling theory for the quasi-deterministic limit

Deterministic rate equations are widely used in the study of stochastic, interacting particles systems. This approach assumes that the inherent noise, associated with the discreteness of the elementary constituents, may be neglected when the number of particles $N$ is large. Accordingly, it fails close to the extinction transition, when the amplitude of stochastic fluctuations is comparable with the size of the population. Here we present a general scaling theory of the transition regime for spatially extended systems. Two fundamental models for out-of-equilibrium phase transitions are considered: the Susceptible-Infected-Susceptible (SIS) that belongs to the directed percolation equivalence class, and the Susceptible-Infected-Recovered (SIR) model belonging to the dynamic percolation class. Implementing the Ginzburg criteria we show that the width of the fluctuation-dominated region scales like $N^{-κ}$, where $N$ is the number of individuals per site and $κ= 2/(d_u-d)$, $d_u$ is the upper critical dimension. Other exponents that control the approach to the deterministic limit are shown to depend on $κ$. The theory is supported by the results of extensive numerical simulations for systems of various dimensionalities.

preprint2012arXiv

Neutral selection

Hubbell's neutral theory of biodiversity has successfully explained the observed composition of many ecological communities but it relies on strict demographic equivalence among species and provides no room for evolutionary processes like selection, adaptation and speciation. Here we show how to embed the neutral theory within the Darwinian framework. In a fitness landscape with a quadratic maximum, typical of quantitative traits, selection restricts the extant species to have traits close to optimal, so that the fitness differences between surviving species are small. For sufficiently small mutation steps, the community structure fits perfectly to the Fisher log-series species abundance distribution. The theory is relatively insensitive to moderate amounts of environmental noise, wherein the location of the fitness maximum changes by amounts of order the width of the noise-free distribution. Adding very large environmental noise to the model qualitatively changes the abundance distributions, converting the exponential fall-off of large species to a power-law decay, typical of a neutral model with environmental noise.

preprint2011arXiv

Noise Effects in Nonlinear Biochemical Signaling

It has been generally recognized that stochasticity can play an important role in the information processing accomplished by reaction networks in biological cells. Most treatments of that stochasticity employ Gaussian noise even though it is a priori obvious that this approximation can violate physical constraints, such as the positivity of chemical concentrations. Here, we show that even when such nonphysical fluctuations are rare, an exact solution of the Gaussian model shows that the model can yield unphysical results. This is done in the context of a simple incoherent-feedforward model which exhibits perfect adaptation in the deterministic limit. We show how one can use the natural separation of time scales in this model to yield an approximate model, that is analytically solvable, including its dynamical response to an environmental change. Alternatively, one can employ a cutoff procedure to regularize the Gaussian result.

preprint2010arXiv

Contact process with a defect: universal oasis, nonuniversal scaling

The extinction transition in the presence of a localized quenched defect is studied numerically. When the bulk is at criticality, the correlation length diverges and even an infinite system cannot "decouple" from the defect. The results presented here suggest that, in 1+1 dimensions, the critical exponent $δ$ that controls the asymptotic power-law decay depends on the strength of the local perturbation. On the other hand, the exponent was found to be independent of the local arrangement of the defect. In higher dimensions the defect seems to induce a transient behavior that decays algebraically in time.

preprint2010arXiv

The Birth-Death-Mutation process: a new paradigm for fat tailed distributions

Fat tailed statistics and power-laws are ubiquitous in many complex systems. Usually the appearance of of a few anomalously successful individuals (bio-species, investors, websites) is interpreted as reflecting some inherent "quality" (fitness, talent, giftedness) as in Darwin's theory of natural selection. Here we adopt the opposite, "neutral", outlook, suggesting that the main factor explaining success is merely luck. The statistics emerging from the neutral birth-death-mutation (BDM) process is shown to fit marvelously many empirical distributions. While previous neutral theories have focused on the power-law tail, our theory economically and accurately explains the entire distribution. We thus suggest the BDM distribution as a standard neutral model: effects of fitness and selection are to be identified by substantial deviations from it.

preprint2009arXiv

Globally coupled chaotic maps and demographic stochasticity

The affect of demographic stochasticity of a system of globally coupled chaotic maps is considered. A two-step model is studied, where the intra-patch chaotic dynamics is followed by a migration step that coupled all patches; the equilibrium number of agents on each site, $N$, controls the strength of the discreteness-induced fluctuations. For small $N$ (large fluctuations) a period-doubling cascade appears as the coupling (migration) increases. As $N$ grows an extremely slow dynamic emerges, leading to a flow along a one-dimensional family of almost period 2 solutions. This manifold become a true solutions in the deterministic limit. The degeneracy between different attractors that characterizes the clustering phase of the deterministic system is thus the $N \to \infty$ limit of the slow dynamics manifold.

preprint2006arXiv

Amplitude dependent frequency, desynchronization, and stabilization in noisy metapopulation dynamics

The enigmatic stability of population oscillations within ecological systems is analyzed. The underlying mechanism is presented in the framework of two interacting species free to migrate between two spatial patches. It is shown that that the combined effects of migration and noise cannot account for the stabilization. The missing ingredient is the dependence of the oscillations' frequency upon their amplitude; with that, noise-induced differences between patches are amplified due to the frequency gradient. Migration among desynchronized regions then stabilizes a "soft" limit cycle in the vicinity of the homogenous manifold. A simple model of diffusively coupled oscillators allows the derivation of quantitative results, like the functional dependence of the desynchronization upon diffusion strength and frequency differences. The oscillations' amplitude is shown to be (almost) noise independent. The results are compared with a numerical integration of the marginally stable Lotka-Volterra equations. An unstable system is extinction-prone for small noise, but stabilizes at larger noise intensity.

preprint2000arXiv

"Quantum phase transitions" in classical nonequilibrium processes

Diffusion limited reaction of the Lotka-Volterra type is analyzed taking into account the discrete nature of the reactants. In the continuum approximation, the dynamics is dominated by an elliptic fixed-point. This fixed-point becomes unstable due to discretization effects, a scenario similar to quantum phase transitions. As a result, the long-time asymptotic behavior of the system changes and the dynamics flows into a limit cycle. The results are verified by numerical simulations.