BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David A. Kessler

David A. Kessler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechPopulations and Evolutionquant-phnlin.PScond-mat.softMolecular Networkscond-mat.mtrl-scimath-phmath.MPSubcellular ProcessesBiological Physicscond-mat.dis-nncond-mat.quant-gasNeurons and Cognitionnlin.CDphysics.atom-phphysics.data-anQuantitative Methods

Catalog footprint

What is connected

42works
18topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

42 published item(s)

preprint2022arXiv

Discrete Sampling of Extreme Events Modifies Their Statistics

Extreme value (EV) statistics of correlated systems are widely investigated in many fields, spanning the spectrum from weather forecasting to earthquake prediction. Does the unavoidable discrete sampling of a continuous correlated stochastic process change its EV distribution? We explore this question for correlated random variables modeled via Langevin dynamics for a particle in a potential field. For potentials growing at infinity faster than linearly and for long measurement times, we find that the EV distribution of the discretely sampled process diverges from that of the full continuous dataset and converges to that of independent and identically distributed random variables drawn from the process's equilibrium measure. However, for processes with sublinear potentials, the long-time limit is the EV statistics of the continuously sampled data. We treat processes whose equilibrium measures belong to the three EV attractors: Gumbel, Fréchet, and Weibull. Our work shows that the EV statistics can be extremely sensitive to the sampling rate of the data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Driving quantum systems with repeated conditional measurements

We investigate the effect of conditional null measurements on a quantum system and find a rich variety of behaviors. Specifically, quantum dynamics with a time independent $H$ in a finite dimensional Hilbert space are considered with repeated strong null measurements of a specified state. We discuss four generic behaviors that emerge in these monitored systems. The first arises in systems without symmetry, along with their associated degeneracies in the energy spectrum, and hence in the absence of dark states as well. In this case, a unique final state can be found which is determined by the largest eigenvalue of the survival operator, the non-unitary operator encoding both the unitary evolution between measurements and the measurement itself. For a three-level system, this is similar to the well known shelving effect. Secondly, for systems with built-in symmetry and correspondingly a degenerate energy spectrum, the null measurements dynamically select the degenerate energy levels, while the non-degenerate levels are effectively wiped out. Thirdly, in the absence of dark states, and for specific choices of parameters, two or more eigenvalues of the survival operator match in magnitude, and this leads to an oscillatory behavior controlled by the measurement rate and not solely by the energy levels. Finally, when the control parameters are tuned, such that the eigenvalues of the survival operator all coalesce to zero, one has exceptional points that corresponds to situations that violate the null measurement condition, making the conditional measurement process impossible.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Non-normalizable quasi-equilibrium solution of the Fokker-Planck equation for nonconfining fields

We investigate the overdamped Langevin motion for particles in a potential well that is asymptotically flat. When the potential well is deep compared to temperature, physical observables like the mean square displacement are essentially time-independent over a long time interval, the stagnation epoch. However the standard Boltzmann-Gibbs (BG) distribution is non-normalizable, given that the usual partition function is divergent. For this regime, we have previously shown that a regularization of BG statistics allows the prediction of the values of dynamical and thermodynamical observables in the non-normalizable quasi-equilibrium state. In this work, based on the eigenfunction expansion of the time-dependent solution of the associated Fokker-Planck equation with free boundary conditions, we obtain an approximate time-independent solution of the BG form, valid for times which are long, but still short compared to the exponentially large escape time. The escaped particles follow a general free-particle statistics, where the solution is a an error function, shifted due to the initial struggle to overcome the potential well. With the eigenfunction solution of the Fokker-Planck equation in hand, we show the validity of the regularized BG statistics and how it perfectly describes the time-independent regime though the quasi-stationary state is non-normalizable.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Infinite ergodic theory meets Boltzmann statistics

We investigate the overdamped stochastic dynamics of a particle in an asymptotically flat external potential field, in contact with a thermal bath. For an infinite system size, the particles may escape the force field and diffuse freely at large length scales. The partition function diverges and hence the standard canonical ensemble fails. This is replaced with tools stemming from infinite ergodic theory. Boltzmann-Gibbs statistics, even though not normalized, still describes integrable observables, like energy and occupation times. The Boltzmann infinite density is derived heuristically using an entropy maximization principle, as well as via a first-principles calculation using an eigenfunction expansion in the continuum of low-energy states. A generalized virial theorem is derived, showing how the virial coefficient describes the delay in the diffusive spreading of the particles, found at large distances. When the process is non-recurrent, e.g. diffusion in three dimensions with a Coulomb-like potential, we use weighted time averages to restore basic canonical relations between time and ensemble averages.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-Hermitian and Zeno limit of quantum systems under rapid measurements

We investigate in depth the relation between the first detection time of an isolated quantum system that is repeatedly perturbed by strong local measurements with a large fixed frequency $1/τ$, determining whether it is in some given state $| ψ_\text{d} \rangle$, and the time of absorption to the same state of the same system with the added imaginary potential $2i\hbar | ψ_\text{d} \rangle \langle ψ_\text{d} | / τ$. As opposed to previous works, we compare directly the solutions of both problems in the small $τ$, i.e., Zeno, limit. We find a scaling collapse in $F(t)$ with respect to $τ$ and compute the total detection probability as well as the moments of the first detection time probability density $F(t)$ in the Zeno limit. We show that both solutions approach the same result in this small $τ$ limit, as long as the initial state $| ψ_\text{in} \rangle$ is not parallel to the detection state, i.e. as long as $| \langle ψ_\text{d} | ψ_\text{in} \rangle | < 1$. However, when this condition is violated, the small probability density to detect the state on time scales much larger than $τ$ is precisely a factor of four different for all such times. We express the solution of the Zeno limit of both problems formally in terms of an electrostatic analogy. Our results are corroborated with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ordered Hexagonal Patterns via Notch-Delta Signaling

Many developmental processes in biology utilize Notch-Delta signaling to construct an ordered pattern of cellular differentiation. This signaling modality is based on nearest-neighbor contact, as opposed to the more familiar mechanism driven by the release of diffusible ligands. Here, exploiting this "juxtacrine" property, we present an exact treatment of the pattern formation problem via a system of nine coupled ordinary differential equations. The possible patterns that are realized for realistic parameters can be analyzed by considering a co-dimension 2 pitchfork bifurcation of this system. This analysis explains the observed prevalence of hexagonal patterns with high Delta at their center, as opposed to those with central high Notch levels. We show that outside this range of parameters, in particular for low cis-coupling, a novel kind of pattern is produced, where high Delta cells have high Notch as well. It also suggests that the biological system is only weakly first order, so that an additional mechanism is required to generate the observed defect-free patterns. We construct a simple strategy for producing such defect-free patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantization of the mean decay time for non-Hermitian quantum systems

We show that the mean time, which a quantum particle needs to escape from a system to the environment, is quantized and independent from most dynamical details of the system. In particular, we consider a quantum system with a general Hermitian Hamiltonian $\hat{H}$ and one decay channel, through which probability dissipates to the environment with rate $Γ$. When the system is initially prepared exactly in the decay state, the mean decay time $\langle T \rangle$ is quantized and equal to $w/(2Γ)$. $w$ is the number of distinct energy levels, i.e. eigenvalues of $\hat{H}$, that have overlap with the decay state, and is also the winding number of a transform of the resolvent in the complex plane. Apart from the integer $w$, $\langle T \rangle$ is completely independent of the system's dynamics. The complete decay time distribution can be obtained from an electrostatic analogy and features rare events of very large dissipation times for parameter choices close to critical points, where $w$ changes, e.g. when a degeneracy is lifted. Experiments of insufficient observation time may thus measure a too small value of $w$. We discuss our findings in a disordered tight-binding model and in the two-level atom in a continuous-wave field.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Saffman-Taylor Fingers at Intermediate Noise

We study Saffman-Taylor flow in the presence of intermediate noise numerically by using both a boundary-integral approach as well as the Kadanoff-Liang modified Diffusion-Limited Aggregation model that incorporates surface tension and reduced noise. For little to no noise, both models result reproduce the well-known Saffman-Taylor finger. We compare both models in the region of intermediate noise where we get occasional tip-splitting events, focusing on the ensemble-average. We show that as the noise in the system is increased, the mean behavior in both models approaches the $\cos^2(πy/W)$ transverse density profile far behind the leading front. We also investigate how the noise scales and affects both models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The first detection time of a quantum state under random probing

We solve for the statistics of the first detection of a quantum system in a particular desired state, when the system is subject to a projective measurement at independent identically distributed random time intervals. We present formulas for the probability of detection in the $n$th attempt. We calculate as well the mean and mean square both of the number of the first successful detection attempt and the time till first detection. We present explicit results for a particle initially localized at a site on a ring of size $L$, probed at some arbitrary given site, in the case when the detection intervals are distributed exponentially. We prove that, for all interval distributions and finite-dimensional Hamiltonians, the mean detection time is equal to the mean attempt number times the mean time interval between attempts. We further prove that for the return problem when the initial and target state are identical, the total detection probability is unity and the mean attempts till detection is an integer, which is the size of the Hilbert space (symmetrized about the target state). We study an interpolation between the fixed time interval case to an exponential distribution of time intervals via the Gamma distribution with constant mean and varying width. The mean arrival time as a function of the mean interval changes qualitatively as we tune the inter-arrival time distribution from very narrow (delta peaked) to exponential, as resonances are wiped out by the randomness of the sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Biological Regulatory Networks are Minimally Frustrated

Characterization of the differences between biological and random networks can reveal the design principles that enable the robust realization of crucial biological functions including the establishment of different cell types. Previous studies, focusing on identifying topological features that are present in biological networks but not in random networks, have, however, provided few functional insights. We use a Boolean modeling framework and ideas from spin glass literature to identify functional differences between five real biological networks and random networks with similar topological features. We show that minimal frustration is a fundamental property that allows biological networks to robustly establish cell types and regulate cell fate choice, and this property can emerge in complex networks via Darwinian evolution. The study also provides clues regarding how the regulation of cell fate choice can go awry in a disease like cancer and lead to the emergence of aberrant cell types.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Uncertainty and symmetry bounds for the quantum total detection probability

We investigate a generic discrete quantum system prepared in state $|ψ_\text{in}\rangle$, under repeated detection attempts aimed to find the particle in state $|d\rangle$, for example a quantum walker on a finite graph searching for a node. For the corresponding classical random walk, the total detection probability $P_\text{det}$ is unity. Due to destructive interference, one may find initial states $|ψ_\text{in}\rangle$ with $P_\text{det}<1$. We first obtain an uncertainty relation which yields insight on this deviation from classical behavior, showing the relation between $P_\text{det}$ and energy fluctuations: $ ΔP \,\mathrm{Var}[\hat{H}]_d \ge | \langle d| [\hat{H}, \hat{D}] | ψ_\text{in} \rangle |^2$ where $ΔP = P_\text{det} - |\langleψ_\text{in}|d\rangle |^2$, and $\hat{D} = |d\rangle\langle d|$ is the measurement projector. Secondly, exploiting symmetry we show that $P_\text{det}\le 1/ν$ where the integer $ν$ is the number of states equivalent to the initial state. These bounds are compared with the exact solution for small systems, obtained from an analysis of the dark and bright subspaces, showing the usefulness of the approach. The upper bounds works well even in large systems, and we show how to tighten the lower bound in this case.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heavy-tailed phase-space distributions beyond Boltzmann-Gibbs and equipartition: Statistics of confined cold atoms

The Boltzmann-Gibbs density, a central result of equilibrium statistical mechanics, relates the energy of a system in contact with a thermal bath to its equilibrium statistics. This relation is lost for non-thermal systems such as cold atoms in optical lattices, where the heat bath is replaced by the laser beams of the lattice. We investigate in detail the stationary phase-space probability for Sisyphus cooling under harmonic confinement. In particular, we elucidate whether the total energy of the system still describes its stationary state statistics. We find that this is true for the center part of the phase-space density for deep lattices, where the Boltzmann-Gibbs density provides an approximate description. The relation between energy and statistics also persists for strong confinement and in the limit of high energies, where the system becomes underdamped. However, the phase-space density now exhibits heavy power-law tails. In all three cases we find expressions for the leading order phase-space density and corrections which break the equivalence of probability and energy and violate energy equipartition. The non-equilibrium nature of the steady state is confounded by explicit violations of detailed balance. We complement these analytical results with numerical simulations to map out the intricate structure of the phase-space density.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear self-adapting wave patterns

We propose a new type of traveling wave pattern, one that can adapt to the size of physical system in which it is embedded. Such a system arises when the initial state has an instability that extends down to zero wavevector, connecting at that point to two symmetry modes of the underlying dynamical system. The Min system of proteins in E. coli is such as system with the symmetry emerging from the global conservation of two proteins, MinD and MinE. For this and related systems, traveling waves can adiabatically deform as the system is increased in size without the increase in node number that would be expected for an oscillatory version of a Turing instability containing an allowed wavenumber band with a finite minimum.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generic Criticality in Ecological and Neuronal Networks

We investigate the dynamics of two models of biological networks with purely suppressive interactions between the units; species interacting via niche competition and neurons via inhibitory synaptic coupling. In both of these cases, power-law scaling of the density of states with probability arises without any fine-tuning of the model parameters. These results argue against the increasingly popular notion that non-equilibrium living systems operate at special critical points, driven by there by evolution so as to enable adaptive processing of input data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Microbranching in mode-I fracture using large scale simulations of amorphous and perturbed lattice models

We study the high-velocity regime mode-I fracture instability when small microbranches start to appear near the main crack, using large scale simulations. Some of the features of those microbranches have been reproduced qualitatively in smaller scale studies (using ${\cal O}(10^4)$ atoms) on both a model of an amorphous materials (via the continuous random network model) and using perturbed lattice models. In this study, larger scale simulations (${\cal O}(10^6)$ atoms) were performed using multi-threading computing on a GPU device, in order to achieve more physically realistic results. First, we find that the microbranching pattern appears to be converging with the lattice width. Second, the simulations reproduce the growth of the size of a microbranch as a function of the crack velocity, as well as the increase of the amplitude of the derivative of the electrical resistance RMS with respect to the time as a function of the crack velocity. In addition, the simulations yield the correct branching angle of the microbranches, and the power law governing the shape of the microbranches seems to be lower than one, so that the side cracks turn over in the direction of propagation of the main crack as seen in experiment.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deviations from Boltzmann-Gibbs equilibrium in confined optical lattices

Cold atoms in dissipative optical lattices have long been known to exhibit anomalous kinetics due to an effective nonlinear friction force. Here we show that confining the spatial motion of the atoms will lead to an anomalous non-Boltzmann-Gibbs equilibrium state, with a power law tail at large energies. Only in the limit of deep optical lattices, do we regain the Boltzmann-Gibbs state. For strong confinement relative to the damping, we find an explicit expression for the equilibrium phase-space distribution, which generally differs from the canonical Boltzmann-Gibbs state at all energies. Both in the low and high energy limits, the equilibrium distribution is a function of the system's Hamiltonian. At intermediate energies, however, the distribution is not a function of energy only and equipartition is violated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fractional Edgeworth Expansion: Corrections to the Gaussian-Lévy Central Limit Theorem

In this article we generalize the classical Edgeworth expansion for the probability density function (PDF) of sums of a finite number of symmetric independent identically distributed random variables with a finite variance to sums of variables with an infinite variance which converge by the generalized central limit theorem to a Lévy $α$-stable density function. Our correction may be written by means of a series of fractional derivatives of the Lévy and the conjugate Lévy PDFs. This series expansion is general and applies also to the Gaussian regime. To describe the terms in the series expansion, we introduce a new family of special functions and briefly discuss their properties. We implement our generalization to the distribution of the momentum for atoms undergoing Sisyphus cooling, and show the improvement of our leading order approximation compared to previous approximations. In vicinity of the transition between Lévy and Gauss behaviors, convergence to asymptotic results slows down.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scaling solution in the large population limit of the general asymmetric stochastic Luria-Delbrück evolution process

One of the most popular models for quantitatively understanding the emergence of drug resistance both in bacterial colonies and in malignant tumors was introduced long ago by Luria and Delbrück. Here, individual resistant mutants emerge randomly during the birth events of an exponentially growing sensitive population. A most interesting limit of this process occurs when the population size $N$ is large and mutation rates are low, but not necessarily small compared to $1/N$. Here we provide a scaling solution valid in this limit, making contact with the theory of Levy $α$-stable distributions, in particular one discussed long ago by Landau. One consequence of this association is that moments of the distribution are highly misleading as far as characterizing typical behavior. A key insight that enables our solution is that working in the fixed population size ensemble is not the same as working in a fixed time ensemble. Some of our results have been presented previously in shortened form.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Distribution of the Area under a Bessel Excursion and its Moments

A Bessel excursion is a Bessel process that begins at the origin and first returns there at some given time $T$. We study the distribution of the area under such an excursion, which recently found application in the context of laser cooling. The area $A$ scales with the time as $A \sim T^{3/2}$, independent of the dimension, $d$, but the functional form of the distribution does depend on $d$. We demonstrate that for $d=1$, the distribution reduces as expected to the distribution for the area under a Brownian excursion, known as the Airy distribution, deriving a new expression for the Airy distribution in the process. We show that the distribution is symmetric in $d-2$, with nonanalytic behavior at $d=2$. We calculate the first and second moments of the distribution, as well as a particular fractional moment. We also analyze the analytic continuation from $d<2$ to $d>2$. In the limit where $d\to 4$ from below, this analytically continued distribution is described by a one-sided Lévy $α$-stable distribution with index $2/3$ and a scale factor proportional to $[(4-d)T]^{3/2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coexistence in an inhomogeneous environment

We examine the two-dimensional extension of the model of Kessler and Sander of competition between two species identical except for dispersion rates. In this class of models, the spatial inhomogeneity of reproduction rates gives rise to an implicit cost of dispersal, due to the tendency to leave favorable locations. Then, as in the Hamilton-May model with its explicit dispersal cost, the tradeoff between dispersal case and the beneficial role of dispersal in limiting fluctuations, leads to an advantage of one dispersal rate over another, and the eventual extinction of the disadvantaged species. In two dimensions we find that while the competition leads to the elimination of one species at high and low population density, at intermediate densities the two species can coexist essentially indefinitely. This is a new phenomenon not present in either the one-dimensional form of the Kessler-Sander model nor in the totally connected Hamilton-May model, and points to the importance of geometry in the question of dispersal.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Micro-branching in mode-I fracture in a randomly perturbed lattice

We study mode-I fracture in lattices with noisy bonds. In contrast to previous attempts, by using a small parameter that perturbs the force-law between the atoms in perfect lattices and using a 3-body force law, simulations reproduce the qualitative behavior of the beyond steady-state cracks in the high velocity regime, including reasonable micro-branching. As far as the physical properties such as the structure factor $g(r)$, the radial or angular distributions, these lattices share the physical properties of perfect lattices rather than that of an amorphous material (e.g., the continuous random network model). A clear transition can be seen between steady-state cracks, where a single crack propagates in the midline of the sample and the regime of unstable cracks, where micro-branches start to appear near the main crack, in line with previous experimental results. This is seen both in a honeycomb lattice and a fully hexagonal lattice. This model reproduces the main physical features of propagating cracks in brittle materials, including the behavior of velocity as a function of driving displacement and the increasing amplitude of oscillations of the electrical resistance. In addition, preliminary indications of power-law behavior of the micro-branch shapes can be seen, potentially reproducing one of the most intriguing experimental results of brittle fracture.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transport and the first passage time problem with application to cold atoms in optical traps

Measurements of spatial diffusion of cold atoms in optical lattices have revealed anomalous super-diffusion, which is controlled by the depth of the optical lattice. We use first passage time statistics to derive the diffusion front of the atoms. In particular, the distributions of areas swept under the first passage curve till its first arrival, and of areas under the Bessel excursion are shown to be powerful tools in the analysis of the atomic cloud. A rather general relation between first passage time statistics and diffusivity is discussed, showing that first passage time analysis is a useful tool in the calculation of transport coefficients. A brief introduction to the semi-classical description of Sisyphus cooling is provided which yields a rich phase diagram for the dynamics.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How input fluctuations reshape the dynamics of a biological switching system

An important task in quantitative biology is to understand the role of stochasticity in biochemical regulation. Here, as an extension of our recent work [Phys. Rev. Lett. 107, 148101 (2011)], we study how input fluctuations affect the stochastic dynamics of a simple biological switch. In our model, the on transition rate of the switch is directly regulated by a noisy input signal, which is described as a nonnegative mean-reverting diffusion process. This continuous process can be a good approximation of the discrete birth-death process and is much more analytically tractable. Within this new setup, we apply the Feynman-Kac theorem to investigate the statistical features of the output switching dynamics. Consistent with our previous findings, the input noise is found to effectively suppress the input-dependent transitions. We show analytically that this effect becomes significant when the input signal fluctuates greatly in amplitude and reverts slowly to its mean.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mass Dependence of Instabilities of an Oscillator with Multiplicative and Additive Noise

We study the instabilities of a harmonic oscillator subject to additive and dichotomous multiplicative noise, focussing on the dependance of the instability threshold on the mass. For multiplicative noise in the damping, the instability threshold is crossed as the mass is decreased, as long as the smaller damping is in fact negative. For multiplicative noise in the stiffness, the situation is more complicated and in fact the transition is reentrant for intermediate noise strength and damping. For multiplicative noise in the mass, the results depend on the implementation of the noise. One can take the velocity or the momentum to be conserved as the mass is changed. In these cases increasing the mass destabilizes the system. Alternatively, if the change in mass is caused by the accretion/loss of particles to the Brownian particle, these processes are asymmetric with momentum conserved upon accretion and velocity upon loss. In this case, there is no instability, as opposed to the other two implementations. We also study the distribution of the energy, finding a power-law cutoff at a value which increases with time.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Theory of fractional-Lévy kinetics for cold atoms diffusing in optical lattices

Recently, anomalous superdiffusion of ultra cold 87Rb atoms in an optical lattice has been observed along with a fat-tailed, Lévy type, spatial distribution. The anomalous exponents were found to depend on the depth of the optical potential. We find, within the framework of the semiclassical theory of Sisyphus cooling, three distinct phases of the dynamics, as the optical potential depth is lowered: normal diffusion; Lévy diffusion; and x ~ t^3/2 scaling, the latter related to Obukhov's model (1959) of turbulence. The process can be formulated as a Lévy walk, with strong correlations between the length and duration of the excursions. We derive a fractional diffusion equation describing the atomic cloud, and the corresponding anomalous diffusion coefficient.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Universal Dimer in a Collisionally Opaque Medium: Experimental Observables and Efimov Resonances

A universal dimer is subject to secondary collisions with atoms when formed in a cloud of ultracold atoms via three-body recombination. We show that in a collisionally opaque medium, the value of the scattering length that results in the maximum number of secondary collisions may not correspond to the Efimov resonance at the atom-dimer threshold and thus can not be automatically associated with it. This result explains a number of controversies in recent experimental results on universal three-body states and supports the emerging evidence for the significant finite range corrections to the first excited Efimov energy level.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Propagating mode-I fracture in amorphous materials using the continuous random network (CRN) model

We study propagating mode-I fracture in two dimensional amorphous materials using atomistic simulations. We used the continuous random network (CRN) model of an amorphous material, creating samples using a two dimensional analogue of the WWW (Wooten, Winer & Weaire) Monte-Carlo algorithm. For modeling fracture, molecular-dynamics simulations were run on the resulting samples. The results of our simulations reproduce the main experimental features. In addition to achieving a steady-state crack under a constant driving displacement (which had not yet been achieved by other atomistic models for amorphous materials), the runs show micro-branching, which increases with driving, transitioning to macro-branching for the largest drivings. Beside the qualitative visual similarity of the simulated cracks to experiment, the simulation also succeeds in explaining the experimentally observed oscillations of the crack velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Infinite invariant densities for anomalous diffusion in optical lattices and other logarithmic potentials

We solve the Fokker-Planck equation for Brownian motion in a logarithmic potential. When the diffusion constant is below a critical value the solution approaches a non-normalizable scaling state, reminiscent of an infinite invariant density. With this non-normalizable density we obtain the phase diagram of anomalous diffusion for this important process. We briefly discuss the consequence for a range of physical systems including atoms in optical lattices and charges in vicinity of long polyelectrolytes. Our work explains in what sense the infinite invariant density and not Boltzmann's equilibrium describes the long time limit of these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-Universal Extinction Transition for Boundary Active Site

We present a generalized model of a diffusion-reaction system where the reaction occurs only on the boundary. This model reduces to that of Barato and Hinrichsen when the occupancy of the boundary site is restricted to zero or one. In the limit when there is no restriction on the occupancy of the boundary site, the model reduces to an age dependent Galton-Watson branching process and admits an analytic solution. The model displays a boundary-induced phase transition into an absorbing state with rational critical exponents and exhibits aging at criticality below a certain fractal dimension of the diffusion process. Surprisingly the behavior in the critical regime for intermediate occupancy restriction $N$ varies with $N$. In fact, by varying the lifetime of the active boundary particle or the diffusion coefficient in the bulk, the critical exponents can be continuously modified.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Frenet algorithm for simulations of fluctuating continuous elastic filaments

We present a new algorithm for generating the equilibrium configurations of fluctuating continuous elastic filaments, based on a combination of statistical mechanics and differential geometry. We use this to calculate the distribution function of the end-to-end distance of filaments with nonvanishing spontaneous curvature and show that for small twist and large bending rigidities, there is an intermediate temperature range in which the filament becomes nearly completely stretched. We show that volume interactions can be incorporated into our algorithm, demonstrating this through the calculation of the effect of excluded volume on the end-to-end distance of the filament.

0 citations0 reviewsNo code linkedOpen access