BZPEER
More

Discover

About BZPEERSearchGraph

Research tools

HomeCollectionsResearchCollaborationPublish

Network

ExpertsCommunities

Opportunities

Listings

Account

Log inCreate account
Log inCreate account

Source author record

Jason Schweinsberg

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record
math.PRPopulations and Evolutionmath-phmath.MPTissues and Organs

Catalog footprint

Research overview

17works
5topics
4close collaborators

Actions

Is this your work?

Log in to claim
Open graphBrowse papers

Published work

17 publications

preprint2022arXiv

Particle configurations for branching Brownian motion with an inhomogeneous branching rate

Aiming to understand the distribution of fitness levels of individuals in a large population undergoing selection, we study the particle configurations of branching Brownian motion where each particle independently moves as Brownian motion with negative drift, particles can die or undergo dyadic fission, and the difference between the birth rate and the death rate is proportional to the particle's location. Under some assumptions, we obtain the limit in probability of the number of particles in any given interval and an explicit formula for the asymptotic empirical density of the fitness distribution. We show that after a sufficiently long time, the fitness distribution from the lowest to the highest fitness levels approximately evolves as a traveling wave with a profile which is asymptotically related the the Airy function. Our work complements the results in Roberts and Schweinsberg (2021), giving a fuller picture of the fitness distribution.

Open access
preprint2021arXiv

A Gaussian particle distribution for branching Brownian motion with an inhomogeneous branching rate

Motivated by the goal of understanding the evolution of populations undergoing selection, we consider branching Brownian motion in which particles independently move according to one-dimensional Brownian motion with drift, each particle may either split into two or die, and the difference between the birth and death rates is a linear function of the position of the particle. We show that, under certain assumptions, after a sufficiently long time, the empirical distribution of the positions of the particles is approximately Gaussian. This provides mathematically rigorous justification for results in the biology literature indicating that the distribution of the fitness levels of individuals in a population over time evolves like a Gaussian traveling wave.

Open access
preprint2020arXiv

ΛΛ-coalescents arising in populations with dormancy

Consider a population evolving from year to year through three seasons: spring, summer and winter. Every spring starts with NN dormant individuals waking up independently of each other according to a given distribution. Once an individual is awake, it starts reproducing at a constant rate. By the end of spring, all individuals are awake and continue reproducing independently as Yule processes during the whole summer. In the winter, NN individuals chosen uniformly at random go to sleep until the next spring, and the other individuals die. We show that because an individual that wakes up unusually early can have a large number of surviving descendants, for some choices of model parameters the genealogy of the population will be described by a ΛΛ-coalescent. In particular, the beta coalescent can describe the genealogy when the rate at which individuals wake up increases exponentially over time. We also characterize the set of all ΛΛ-coalescents that can arise in this framework.

Open access
preprint2020arXiv

Mutation timing in a spatial model of evolution

Motivated by models of cancer formation in which cells need to acquire kk mutations to become cancerous, we consider a spatial population model in which the population is represented by the dd-dimensional torus of side length LL. Initially, no sites have mutations, but sites with i−1i-1 mutations acquire an iith mutation at rate μiμ_i per unit area. Mutations spread to neighboring sites at rate αα, so that tt time units after a mutation, the region of individuals that have acquired the mutation will be a ball of radius αtαt. We calculate, for some ranges of the parameter values, the asymptotic distribution of the time required for some individual to acquire kk mutations. Our results, which build on previous work of Durrett, Foo, and Leder, are essentially complete when k=2k = 2 and when μi=μμ_i = μ for all ii.

Open access
preprint2015arXiv

Rigorous results for a population model with selection I: evolution of the fitness distribution

We consider a model of a population of fixed size NN undergoing selection. Each individual acquires beneficial mutations at rate μNμ_N, and each beneficial mutation increases the individual's fitness by sNs_N. Each individual dies at rate one, and when a death occurs, an individual is chosen with probability proportional to the individual's fitness to give birth. Under certain conditions on the parameters μNμ_N and sNs_N, we obtain rigorous results for the rate at which mutations accumulate in the population and the distribution of the fitnesses of individuals in the population at a given time. Our results confirm predictions of Desai and Fisher (2007).

Open access
preprint2015arXiv

Rigorous results for a population model with selection II: genealogy of the population

We consider a model of a population of fixed size NN undergoing selection. Each individual acquires beneficial mutations at rate μNμ_N, and each beneficial mutation increases the individual's fitness by sNs_N. Each individual dies at rate one, and when a death occurs, an individual is chosen with probability proportional to the individual's fitness to give birth. Under certain conditions on the parameters μNμ_N and sNs_N, we show that the genealogy of the population can be described by the Bolthausen-Sznitman coalescent. This result confirms predictions of Desai, Walczak, and Fisher (2013), and Neher and Hallatschek (2013).

Open access
preprint2014arXiv

The evolving beta coalescent

In mathematical population genetics, it is well known that one can represent the genealogy of a population by a tree, which indicates how the ancestral lines of individuals in the population coalesce as they are traced back in time. As the population evolves over time, the tree that represents the genealogy of the population also changes, leading to a tree-valued stochastic process known as the evolving coalescent. Here we will consider the evolving coalescent for populations whose genealogy can be described by a beta coalescent, which is known to give the genealogy of populations with very large family sizes. We show that as the size of the population tends to infinity, the evolution of certain functionals of the beta coalescent, such as the total number of mergers, the total branch length, and the total length of external branches, converges to a stationary stable process. Our methods also lead to new proofs of known asymptotic results for certain functionals of the non-evolving beta coalescent.

Open access
preprint2013arXiv

Critical branching Brownian motion with absorption: particle configurations

We consider critical branching Brownian motion with absorption, in which there is initially a single particle at x>0x > 0, particles move according to independent one-dimensional Brownian motions with the critical drift of −2-\sqrt{2}, and particles are absorbed when they reach zero. Here we obtain asymptotic results concerning the behavior of the process before the extinction time, as the position xx of the initial particle tends to infinity. We estimate the number of particles in the system at a given time and the position of the right-most particle. We also obtain asymptotic results for the configuration of particles at a typical time.

Open access
preprint2013arXiv

The genealogy of branching Brownian motion with absorption

We consider a system of particles which perform branching Brownian motion with negative drift and are killed upon reaching zero, in the near-critical regime where the total population stays roughly constant with approximately N particles. We show that the characteristic time scale for the evolution of this population is of order (log⁡N)3(\log N)^3, in the sense that when time is measured in these units, the scaled number of particles converges to a variant of Neveu's continuous-state branching process. Furthermore, the genealogy of the particles is then governed by a coalescent process known as the Bolthausen-Sznitman coalescent. This validates the nonrigorous predictions by Brunet, Derrida, Muller and Munier for a closely related model.

Open access
preprint2012arXiv

Critical branching Brownian motion with absorption: survival probability

We consider branching Brownian motion on the real line with absorption at zero, in which particles move according to independent Brownian motions with the critical drift of −2-\sqrt{2}. Kesten (1978) showed that almost surely this process eventually dies out. Here we obtain upper and lower bounds on the probability that the process survives until some large time tt. These bounds improve upon results of Kesten (1978), and partially confirm nonrigorous predictions of Derrida and Simon (2007).

Open access
preprint2011arXiv

Consensus in the two-state Axelrod model

The Axelrod model is a spatial stochastic model for the dynamics of cultures which, similarly to the voter model, includes social influence, but differs from the latter by also accounting for another social factor called homophily, the tendency to interact more frequently with individuals who are more similar. Each individual is characterized by its opinions about a finite number of cultural features, each of which can assume the same finite number of states. Pairs of adjacent individuals interact at a rate equal to the fraction of features they have in common, thus modeling homophily, which results in the interacting pair having one more cultural feature in common, thus modeling social influence. It has been conjectured based on numerical simulations that the one-dimensional Axelrod model clusters when the number of features exceeds the number of states per feature. In this article, we prove this conjecture for the two-state model with an arbitrary number of features.

Open access
preprint2011arXiv

Dynamics of the evolving Bolthausen-Sznitman coalescent

Consider a population of fixed size that evolves over time. At each time, the genealogical structure of the population can be described by a coalescent tree whose branches are traced back to the most recent common ancestor of the population. As time goes forward, the genealogy of the population evolves, leading to what is known as an evolving coalescent. We will study the evolving coalescent for populations whose genealogy can be described by the Bolthausen-Sznitman coalescent. We obtain the limiting behavior of the evolution of the time back to the most recent common ancestor and the total length of the branches in the tree. By similar methods, we also obtain a new result concerning the number of blocks in the Bolthausen-Sznitman coalescent.

Open access
preprint2010arXiv

Survival of near-critical branching Brownian motion

Consider a system of particles performing branching Brownian motion with negative drift μ=2−εμ= \sqrt{2 - ε} and killed upon hitting zero. Initially there is one particle at x>0x>0. Kesten showed that the process survives with positive probability if and only if ε>0ε>0. Here we are interested in the asymptotics as $\eps\to 0$ of the survival probability Qμ(x)Q_μ(x). It is proved that if L=π/εL= π/\sqrtε then for all x∈Rx \in \R, lim⁡ε→0Qμ(L+x)=θ(x)∈(0,1)\lim_{ε\to 0} Q_μ(L+x) = θ(x) \in (0,1) exists and is a travelling wave solution of the Fisher-KPP equation. Furthermore, we obtain sharp asymptotics of the survival probability when x<Lx<L and L−x→∞L-x \to \infty. The proofs rely on probabilistic methods developed by the authors in a previous work. This completes earlier work by Harris, Harris and Kyprianou and confirms predictions made by Derrida and Simon, which were obtained using nonrigorous PDE methods.

Open access
preprint2010arXiv

The number of small blocks in exchangeable random partitions

Suppose ΠΠ is an exchangeable random partition of the positive integers and ΠnΠ_n is its restriction to {1,...,n}\{1, ..., n\}. Let KnK_n denote the number of blocks of ΠnΠ_n, and let Kn,rK_{n,r} denote the number of blocks of ΠnΠ_n containing rr integers. We show that if 0<α<10 < α< 1 and Kn/(nαℓ(n))K_n/(n^α \ell(n)) converges in probability to Γ(1−α)Γ(1-α), where ℓ\ell is a slowly varying function, then Kn,r/(nαℓ(n))K_{n,r}/(n^α \ell(n)) converges in probability to αΓ(r−α)/r!αΓ(r - α)/r!. This result was previously known when the convergence of Kn/(nαℓ(n))K_n/(n^α \ell(n)) holds almost surely, but the result under the hypothesis of convergence in probability has significant implications for coalescent theory. We also show that a related conjecture for the case when KnK_n grows only slightly slower than nn fails to be true.

Open access
preprint2008arXiv

The waiting time for m mutations

We consider a model of a population of fixed size N in which each individual gets replaced at rate one and each individual experiences a mutation at rate μ. We calculate the asymptotic distribution of the time that it takes before there is an individual in the population with m mutations. Several different behaviors are possible, depending on how μchanges with N. These results have applications to the problem of determining the waiting time for regulatory sequences to appear and to models of cancer development.

Open access
preprint2007arXiv

Beta-coalescents and continuous stable random trees

Coalescents with multiple collisions, also known as ΛΛ-coalescents, were introduced by Pitman and Sagitov in 1999. These processes describe the evolution of particles that undergo stochastic coagulation in such a way that several blocks can merge at the same time to form a single block. In the case that the measure ΛΛ is the Beta⁡(2−α,α)\operatorname {Beta}(2-α,α) distribution, they are also known to describe the genealogies of large populations where a single individual can produce a large number of offspring. Here, we use a recent result of Birkner et al. to prove that Beta-coalescents can be embedded in continuous stable random trees, about which much is known due to the recent progress of Duquesne and Le Gall. Our proof is based on a construction of the Donnelly--Kurtz lookdown process using continuous random trees, which is of independent interest. This produces a number of results concerning the small-time behavior of Beta-coalescents. Most notably, we recover an almost sure limit theorem of the present authors for the number of blocks at small times and give the multifractal spectrum corresponding to the emergence of blocks with atypical size. Also, we are able to find exact asymptotics for sampling formulae corresponding to the site frequency spectrum and the allele frequency spectrum associated with mutations in the context of population genetics.

Open access

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.