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

Rick Durrett

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

ResearcherUnclaimed source record
math.PRphysics.soc-phPopulations and EvolutionSocial and Information NetworksTissues and Organscond-mat.dis-nncond-mat.stat-mechmath.COnlin.AO

Catalog footprint

What is connected

19works
9topics
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

19 published item(s)

preprint2020arXiv

Controlling the spread of COVID-19 on college campuses

This research was done during the DOMath program at Duke University from May 18 to July 10, 2020. At the time, Duke and other universities across the country were wrestling with the question of how to safely welcome students back to campus in the Fall. Because of this, our project focused on using mathematical models to evaluate strategies to suppress the spread of the virus on campus, specifically in dorms and in classrooms. For dorms, we show that giving students single rooms rather than double rooms can substantially reduce virus spread. For classrooms, we show that moving classes with size above some cutoff online can make the basic reproduction number $R_0<1$, preventing a wide spread epidemic. The cutoff will depend on the contagiousness of the disease in classrooms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The q-voter model on the torus

In the $q$-voter model, the voter at $x$ changes its opinion at rate $f_x^q$, where $f_x$ is the fraction of neighbors with the opposite opinion. Mean-field calculations suggest that there should be coexistence between opinions if $q<1$ and clustering if $q>1$. This model has been extensively studied by physicists, but we do not know of any rigorous results. In this paper, we use the machinery of voter model perturbations to show that the conjectured behavior holds for $q$ close to 1. More precisely, we show that if $q<1$, then for any $m<\infty$ the process on the three-dimensional torus with $n$ points survives for time $n^m$, and after an initial transient phase has a density that it is always close to 1/2. If $q>1$, then the process rapidly reaches fixation on one opinion. It is interesting to note that in the second case the limiting ODE (on its sped up time scale) reaches 0 at time $\log n$ but the stochastic process on the same time scale dies out at time $(1/3)\log n$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Genealogies in Expanding Populations

The goal of this paper is to prove rigorous results for the behavior of genealogies in a one-dimensional long range biased voter model introduced by Hallatschek and Nelson [25]. The first step, which is easily accomplished using results of Mueller and Tribe [38], is to show that when space and time are rescaled correctly, our biased voter model converges to a Wright-Fisher SPDE. A simple extension of a result of Durrett and Restrepo [18] then shows that the dual branching coalescing random walk converges to a branching Brownian motion in which particles coalesce after an exponentially distributed amount of intersection local time. Brunet et al. [8] have conjectured that genealogies in models of this type are described by the Bolthausen-Sznitman coalescent, see [39]. However, in the model we study there are no simultaneous coalescences. Our third and most significant result concerns "tracer dynamics" in which some of the initial particles in the biased voter model are labeled. We show that the joint distribution of the labeled and unlabeled particles converges to the solution of a system of stochastic partial differential equations. A new duality equation that generalizes the one Shiga [44] developed for the Wright-Fisher SPDE is the key to the proof of that result.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Latent Voter Model on Locally Tree-Like Random Graphs

In the latent voter model, which models the spread of a technology through a social network, individuals who have just changed their choice have a latent period, which is exponential with rate $λ$, during which they will not buy a new device. We study site and edge versions of this model on random graphs generated by a configuration model in which the degrees $d(x)$ have $3 \le d(x) \le M$. We show that if the number of vertices $n \to\infty$ and $\log n \ll λ_n \ll n$ then the latent voter model has a quasi-stationary state in which each opinion has probability $\approx 1/2$ and persists in this state for a time that is $\ge n^m$ for any $m<\infty$. Thus, even a very small latent period drastically changes the behavior of the voter model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Evolving Voter Model on Thick Graphs

In the evolving voter model, when an individual interacts with a neighbor having an opinion different from theirs, they will with probability $1-α$ imitate the neighbor but with probability $ α$ will sever the connection and choose a new neighbor at random (i) from the graph or (ii) from those with the same opinion. Durrett et al. used simulation and heuristics to study these dynamics on sparse graphs. Recently Basu and Sly have studied this system with $1-α= ν/N$ on a dense Erdős-Rényi graph $G(N,1/2)$ and rigorously proved that there is a phase transition from rapid disconnection into components with a single opinion to prolonged persistence of discordant edges as $ν$ increases. In this paper, we consider the intermediate situation of Erdős-Rényi random graphs with average degree $L=N^a$ where $0 < a < 1$. Most of the paper is devoted to a rigorous analysis of an approximation of the dynamics called the approximate master equation. Using ideas of \cite{LMR} and \cite{Silk} we are able to analyze these dynamics in great detail.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coexistence of grass, saplings and trees in the Staver-Levin forest model

In this paper, we consider two attractive stochastic spatial models in which each site can be in state 0, 1 or 2: Krone's model in which 0${}={}$vacant, 1${}={}$juvenile and 2${}={}$a mature individual capable of giving birth, and the Staver-Levin forest model in which 0${}={}$grass, 1${}={}$sapling and 2${}={}$tree. Our first result shows that if $(0,0)$ is an unstable fixed point of the mean-field ODE for densities of 1's and 2's then when the range of interaction is large, there is positive probability of survival starting from a finite set and a stationary distribution in which all three types are present. The result we obtain in this way is asymptotically sharp for Krone's model. However, in the Staver-Levin forest model, if $(0,0)$ is attracting then there may also be another stable fixed point for the ODE, and in some of these cases there is a nontrivial stationary distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evolutionary Games on the Torus with Weak Selection

We study evolutionary games on the torus with $N$ points in dimensions $d\ge 3$. The matrices have the form $\bar G = {\bf 1} + w G$, where ${\bf 1}$ is a matrix that consists of all 1's, and $w$ is small. As in Cox Durrett and Perkins \cite{CDP} we rescale time and space and take a limit as $N\to\infty$ and $w\to 0$. If (i) $w \gg N^{-2/d}$ then the limit is a PDE on ${\bf R}^d$. If (ii) $N^{-2/d} \gg w \gg N^{-1}$, then the limit is an ODE. If (iii) $w \ll N^{-1}$ then the effect of selection vanishes in the limit. In regime (ii) if we introduce a mutation $μ$ so that $μ/w \to \infty$ slowly enough then we arrive at Tarnita's formula that describes how the equilibrium frequencies are shifted due to selection.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spatial Evolutionary Games with small selection coefficients

Here we will use results of Cox, Durrett, and Perkins for voter model perturbations to study spatial evolutionary games on $Z^d$, $d\ge 3$ when the interaction kernel is finite range, symmetric, and has covariance matrix $σ^2I$. The games we consider have payoff matrices of the form ${\bf 1} + wG$ where ${\bf 1}$ is matrix of all 1's and $w$ is small and positive. Since our population size $N=\infty$, we call our selection small rather than weak which usually means $w =O(1/N)$. The key to studying these games is the fact that when the dynamics are suitably rescaled in space and time they convergence to solutions of a reaction diffusion equation (RDE). Inspired by work of Ohtsuki and Nowak and Tarnita et al we show that the reaction term is the replicator equation for a modified game matrix and the modifications of the game matrix depend on the interaction kernel only through the values of two or three simple probabilities for an associated coalescing random walk. Two strategy games lead to an RDE with a cubic nonlinearity, so we can describe the phase diagram completely. Three strategy games lead to a pair of coupled RDE, but using an idea from our earlier work, we are able to show that if there is a repelling function for the replicator equation for the modified game, then there is coexistence in the spatial game when selection is small. This enables us to prove coexistence in the spatial model in a wide variety of examples where the replicator equation of the modified game has an attracting equilibrium with all components positive. Using this result we are able to analyze the behavior of four evolutionary games that have recently been used in cancer modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spatial networks evolving to reduce length

Motivated by results of Henry, Pralat and Zhang (PNAS 108.21 (2011): 8605-8610), we propose a general scheme for evolving spatial networks in order to reduce their total edge lengths. We study the properties of the equilbria of two networks from this class, which interpolate between three well studied objects: the Erdős-Rényi random graph, the random geometric graph, and the minimum spanning tree. The first of our two evolutions can be used as a model for a social network where individuals have fixed opinions about a number of issues and adjust their ties to be connected to people with similar views. The second evolution which preserves the connectivity of the network has potential applications in the design of transportation networks and other distribution systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A first order phase transition in the threshold-$θ\ge 2$ contact process on random $r$-regular graphs and $r$-trees

We consider the discrete-time threshold-$θ\ge 2$ contact process on a random r-regular graph on n vertices. In this process, a vertex with at least θoccupied neighbors at time t will be occupied at time t+1 with probability p, and vacant otherwise. We show that if $θ\ge 2$ and $r \ge θ+2$, $ε_1$ is small and p is at least $p_1(ε_1)$, then starting from all vertices occupied the fraction of occupied vertices stays above $1-2ε_1$ up to time $\exp(γ_1(r)n)$ with probability at least $1 - \exp(-γ_1(r)n)$. In the other direction, we show that for $p_2 < 1$ there is an $ε_2(p_2)>0$ so that if $p \le p_2$ and the number of occupied vertices in the initial configuration is at most $ε_2(p_2)n$, then with high probability all vertices are vacant at time $C_2(p_2) \log(n)$. These two conclusions imply that on the random r-regular graph there cannot be a quasi-stationary distribution with density of occupied vertices between 0 and $ε_2(p_1)$, and allow us to conclude that the process on the r-tree has a first order phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A multi-opinion evolving voter model with infinitely many phase transitions

We consider an idealized model in which individuals' changing opinions and their social network coevolve, with disagreements between neighbors in the network resolved either through one imitating the opinion of the other or by reassignment of the discordant edge. Specifically, an interaction between $x$ and one of its neighbors $y$ leads to $x$ imitating $y$ with probability $(1-α)$ and otherwise (i.e., with probability $α$) $x$ cutting its tie to $y$ in order to instead connect to a randomly chosen individual. Building on previous work about the two-opinion case, we study the multiple-opinion situation, finding that the model has infinitely many phase transitions. Moreover, the formulas describing the end states of these processes are remarkably simple when expressed as a function of $β= α/(1-α)$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase Transitions in the Quadratic Contact Process on Complex Networks

The quadratic contact process (QCP) is a natural extension of the well studied linear contact process where infected (1) individuals infect susceptible (0) neighbors at rate $λ$ and infected individuals recover ($1 \longrightarrow 0$) at rate 1. In the QCP, a combination of two 1's is required to effect a $0 \longrightarrow 1$ change. We extend the study of the QCP, which so far has been limited to lattices, to complex networks. \comment{as a model for the change in a population through sexual reproduction and death.} We define two versions of the QCP -- vertex centered (VQCP) and edge centered (EQCP) with birth events $1-0-1 \longrightarrow 1-1-1$ and $1-1-0 \longrightarrow 1-1-1$ respectively, where `$-$' represents an edge. We investigate the effects of network topology by considering the QCP on random regular, Erdős-Rényi and power law random graphs. We perform mean field calculations as well as simulations to find the steady state fraction of occupied vertices as a function of the birth rate. We find that on the random regular and Erdős-Rényi graphs, there is a discontinuous phase transition with a region of bistability, whereas on the heavy tailed power law graph, the transition is continuous. The critical birth rate is found to be positive in the former but zero in the latter.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Population genetics of neutral mutations in exponentially growing cancer cell populations

In order to analyze data from cancer genome sequencing projects, we need to be able to distinguish causative, or "driver," mutations from "passenger" mutations that have no selective effect. Toward this end, we prove results concerning the frequency of neutural mutations in exponentially growing multitype branching processes that have been widely used in cancer modeling. Our results yield a simple new population genetics result for the site frequency spectrum of a sample from an exponentially growing population.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotic behavior of Aldous' gossip process

Aldous [(2007) Preprint] defined a gossip process in which space is a discrete $N\times N$ torus, and the state of the process at time $t$ is the set of individuals who know the information. Information spreads from a site to its nearest neighbors at rate 1/4 each and at rate $N^{-α}$ to a site chosen at random from the torus. We will be interested in the case in which $α<3$, where the long range transmission significantly accelerates the time at which everyone knows the information. We prove three results that precisely describe the spread of information in a slightly simplified model on the real torus. The time until everyone knows the information is asymptotically $T=(2-2α/3)N^{α/3}\log N$. If $ρ_s$ is the fraction of the population who know the information at time $s$ and $\varepsilon$ is small then, for large $N$, the time until $ρ_s$ reaches $\varepsilon$ is $T(\varepsilon)\approx T+N^{α/3}\log (3\varepsilon /M)$, where $M$ is a random variable determined by the early spread of the information. The value of $ρ_s$ at time $s=T(1/3)+tN^{α/3}$ is almost a deterministic function $h(t)$ which satisfies an odd looking integro-differential equation. The last result confirms a heuristic calculation of Aldous.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Traveling waves of selective sweeps

The goal of cancer genome sequencing projects is to determine the genetic alterations that cause common cancers. Many malignancies arise during the clonal expansion of a benign tumor which motivates the study of recurrent selective sweeps in an exponentially growing population. To better understand this process, Beerenwinkel et al. [PLoS Comput. Biol. 3 (2007) 2239--2246] consider a Wright--Fisher model in which cells from an exponentially growing population accumulate advantageous mutations. Simulations show a traveling wave in which the time of the first $k$-fold mutant, $T_k$, is approximately linear in $k$ and heuristics are used to obtain formulas for $ET_k$. Here, we consider the analogous problem for the Moran model and prove that as the mutation rate $μ\to0$, $T_k\sim c_k\log(1/μ)$, where the $c_k$ can be computed explicitly. In addition, we derive a limiting result on a log scale for the size of $X_k(t)={}$the number of cells with $k$ mutations at time $t$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evolution in predator-prey systems

We study the adaptive dynamics of predator-prey systems modeled by a dynamical system in which the traits of predators and prey are allowed to evolve by small mutations. When only the prey are allowed to evolve, and the size of the mutational change tends to 0, the system does not exhibit long term prey coexistence and the trait of the resident prey type converges to the solution of an ODE. When only the predators are allowed to evolve, coexistence of predators occurs. In this case, depending on the parameters being varied, we see (i) the number of coexisting predators remains tight and the differences in traits from a reference species converge in distribution to a limit, or (ii) the number of coexisting predators tends to infinity, and we calculate the asymptotic rate at which the traits of the least and most "fit" predators in the population increase. This last result is obtained by comparison with a branching random walk killed to the left of a linear boundary and a finite branching-selection particle system.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evolutionary dynamics of tumor progression with random fitness values

Most human tumors result from the accumulation of multiple genetic and epigenetic alterations in a single cell. Mutations that confer a fitness advantage to the cell are known as driver mutations and are causally related to tumorigenesis. Other mutations, however, do not change the phenotype of the cell or even decrease cellular fitness. While much experimental effort is being devoted to the identification of the different functional effects of individual mutations, mathematical modeling of tumor progression generally considers constant fitness increments as mutations are accumulated. In this paper we study a mathematical model of tumor progression with random fitness increments. We analyze a multi-type branching process in which cells accumulate mutations whose fitness effects are chosen from a distribution. We determine the effect of the fitness distribution on the growth kinetics of the tumor. This work contributes to a quantitative understanding of the accumulation of mutations leading to cancer phenotypes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A waiting time problem arising from the study of multi-stage carcinogenesis

We consider the population genetics problem: how long does it take before some member of the population has $m$ specified mutations? The case $m=2$ is relevant to onset of cancer due to the inactivation of both copies of a tumor suppressor gene. Models for larger $m$ are needed for colon cancer and other diseases where a sequence of mutations leads to cells with uncontrolled growth.

0 citations0 reviewsNo code linkedOpen access