Showing 1–50 of 116 results for author: Hollander, F d

Self-repellent branching random walk

Authors: Anton Bovier, Lisa Hartung, Frank den Hollander

Abstract: We consider a discrete-time binary branching random walk with independent standard normal increments subject to a penalty $\b$ for every pair of particles that get within distance $\e$ of each other at any time. We give a precise description of the most likely configurations of the particles under this law for $N$ large and $\b,\e$ fixed. Particles spread out over a distance $2^{2N/3}$, essentiall… ▽ More We consider a discrete-time binary branching random walk with independent standard normal increments subject to a penalty $\b$ for every pair of particles that get within distance $\e$ of each other at any time. We give a precise description of the most likely configurations of the particles under this law for $N$ large and $\b,\e$ fixed. Particles spread out over a distance $2^{2N/3}$, essentially in finite time, and subsequently arrange themselves so that at time $2N/3$ they cover a grid of width $\e$ with one particle per site. After time $2N/3$, the bulk of the particles and their descendants do not move anymore, while the particles in a boundary layer of width $2^{N/3}$ form a ``staircase" to the particles in the bulk. At time $N$, each site in the boundary layer is occupied by $2^{N/3}$ particles. △ Less

Submitted 22 July, 2024; originally announced July 2024.

Comments: 20 pages, 4 figures

MSC Class: 60J80; 60G70; 82B44

Tame sparse exponential random graphs

Authors: Suman Chakraborty, Remco van der Hofstad, Frank den Hollander

Abstract: In this paper, we obtain a precise estimate of the probability that the sparse binomial random graph contains a large number of vertices in a triangle. The estimate of log of this probability is correct up to second order, and enables us to propose an exponential random graph model based on the number of vertices in a triangle. Specifically, by tuning a single parameter, we can with high probabili… ▽ More In this paper, we obtain a precise estimate of the probability that the sparse binomial random graph contains a large number of vertices in a triangle. The estimate of log of this probability is correct up to second order, and enables us to propose an exponential random graph model based on the number of vertices in a triangle. Specifically, by tuning a single parameter, we can with high probability induce any given fraction of vertices in a triangle. Moreover, in the proposed exponential random graph model we derive the large deviation principle for the number of edges. As a byproduct, we propose a consistent estimator of the tuning parameter. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Follow up of arXiv:2112.06526v2

MSC Class: 05C80; 60F10; 62F12

Homogeneous nucleation for two-dimensional Kawasaki dynamics

Authors: Simone Baldassarri, Alexandre Gaudillière, Frank den Hollander, Francesca Romana Nardi, Enzo Olivieri, Elisabetta Scoppola

Abstract: This is the third in a series of three papers in which we study a lattice gas subject to Kawasaki dynamics at inverse temperature $β>0$ in a large finite box $Λ_β\subset \mathbb{Z}^2$ whose volume depends on $β$. Each pair of neighbouring particles has a negative binding energy $-U<0$, while each particle has a positive activation energy $Δ>0$. The initial configuration is drawn from the grand-can… ▽ More This is the third in a series of three papers in which we study a lattice gas subject to Kawasaki dynamics at inverse temperature $β>0$ in a large finite box $Λ_β\subset \mathbb{Z}^2$ whose volume depends on $β$. Each pair of neighbouring particles has a negative binding energy $-U<0$, while each particle has a positive activation energy $Δ>0$. The initial configuration is drawn from the grand-canonical ensemble restricted to the set of configurations where all the droplets are subcritical. Our goal is to describe, in the metastable regime $Δ\in (U,2U)$ and in the limit as $β\to\infty$, how and when the system nucleates, i.e., creates a critical droplet somewhere in $Λ_β$ that subsequently grows by absorbing particles from the surrounding gas. In the first paper we showed that subcritical droplets behave as quasi-random walks. In the second paper we used the results in the first paper to analyse how subcritical droplets form and dissolve on multiple space-time scales when the volume is moderately large, namely, $|Λ_β| = \mathrm{e}^{θβ}$ with $Δ< θ< 2Δ-U$. In the present paper we consider the setting where the volume is very large, namely, $|Λ_β| = \mathrm{e}^{Θβ}$ with $Δ< Θ< Γ-(2Δ-U)$, where $Γ$ is the energy of the critical droplet in the local model with fixed volume, and use the results in the first two papers to identify the nucleation time and the tube of typical trajectories towards nucleation. We will see that in a very large volume critical droplets appear more or less independently in boxes of moderate volume, a phenomenon referred to as homogeneous nucleation. One of the key ingredients in the proof is an estimate showing that no information can travel between these boxes on relevant time scales. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 21 pages, 5 figures. arXiv admin note: text overlap with arXiv:2304.14099

MSC Class: 60K35; 82C26; 82C27

Mixing of fast random walks on dynamic random permutations

Authors: Luca Avena, Remco van der Hofstad, Frank den Hollander, Oliver Nagy

Abstract: We analyse the mixing profile of a random walk on a dynamic random permutation, focusing on the regime where the walk evolves much faster than the permutation. Two types of dynamics generated by random transpositions are considered: one allows for coagulation of permutation cycles only, the other allows for both coagulation and fragmentation. We show that for both types, after scaling time by the… ▽ More We analyse the mixing profile of a random walk on a dynamic random permutation, focusing on the regime where the walk evolves much faster than the permutation. Two types of dynamics generated by random transpositions are considered: one allows for coagulation of permutation cycles only, the other allows for both coagulation and fragmentation. We show that for both types, after scaling time by the length of the permutation and letting this length tend to infinity, the total variation distance between the current distribution and the uniform distribution converges to a limit process that drops down in a single jump. This jump is similar to a one-sided cut-off, occurs after a random time whose law we identify, and goes from the value 1 to a value that is a strictly decreasing and deterministic function of the time of the jump, related to the size of the largest component in Erdős-Rényi random graphs. After the jump, the total variation distance follows this function down to 0. △ Less

Submitted 29 February, 2024; originally announced March 2024.

The Moran model with random resampling rates

Authors: Siva Athreya, Frank den Hollander, Adrian Röllin

Abstract: In this paper we consider the two-type Moran model with $N$ individuals. Each individual is assigned a resampling rate, drawn independently from a probability distribution ${\mathbb P}$ on ${\mathbb R}_+$, and a type, either $1$ or $0$. Each individual resamples its type at its assigned rate, by adopting the type of an individual drawn uniformly at random. Let $Y^N(t)$ denote the empirical distrib… ▽ More In this paper we consider the two-type Moran model with $N$ individuals. Each individual is assigned a resampling rate, drawn independently from a probability distribution ${\mathbb P}$ on ${\mathbb R}_+$, and a type, either $1$ or $0$. Each individual resamples its type at its assigned rate, by adopting the type of an individual drawn uniformly at random. Let $Y^N(t)$ denote the empirical distribution of the resampling rates of the individuals with type $1$ at time $Nt$. We show that if ${\mathbb P}$ has countable support and satisfies certain tail and moment conditions, then in the limit as $N\to\infty$ the process $(Y^N(t))_{t \geq 0}$ converges in law to the process $(S(t)\,¶)_{t \geq 0}$, in the so-called Meyer-Zheng topology, where $(S(t))_{t \geq 0}$ is the Fisher-Wright diffusion with diffusion constant $D$ given by $1/D = \int_{{\mathbb R}_+} (1/r)\,{\mathbb P}(\mathrm{d} r)$. △ Less

Submitted 2 February, 2024; originally announced February 2024.

MSC Class: 60J70; 60K35; 92D25

The friendship paradox for sparse random graphs

Authors: Rajat Subhra Hazra, Frank den Hollander, Azadeh Parvaneh

Abstract: Let $G_n$ be an undirected finite graph on $n\in\mathbb{N}$ vertices labelled by $[n] = \{1,\ldots,n\}$. For $i \in [n]$, let $Δ_{i,n}$ be the friendship bias of vertex $i$, defined as the difference between the average degree of the neighbours of vertex $i$ and the degree of vertex $i$ itself when $i$ is not isolated, and zero when $i$ is isolated. Let $μ_n$ denote the friendship-bias empirical d… ▽ More Let $G_n$ be an undirected finite graph on $n\in\mathbb{N}$ vertices labelled by $[n] = \{1,\ldots,n\}$. For $i \in [n]$, let $Δ_{i,n}$ be the friendship bias of vertex $i$, defined as the difference between the average degree of the neighbours of vertex $i$ and the degree of vertex $i$ itself when $i$ is not isolated, and zero when $i$ is isolated. Let $μ_n$ denote the friendship-bias empirical distribution, i.e., the measure that puts mass $\frac{1}{n}$ at each $Δ_{i,n}$, $i \in [n]$. The friendship paradox says that if $G_n$ has no self-loops, then $\int_{\mathbb{R}} xμ_n(\mathrm{d}x) \geq 0$, with equality if and only if in each connected component of $G_n$ all the degrees are the same. We show that if $(G_n)_{n\in\mathbb{N}}$ is a sequence of sparse random graphs that converges to a rooted random tree in the sense of convergence locally in probability, then $μ_n$ converges weakly to a limiting measure $μ$ that is expressible in terms of the law of the rooted random tree. We study $μ$ for four classes of sparse random graphs: the homogeneous Erdős-Rényi random graph, the inhomogeneous Erdős-Rényi random graph, the configuration model and the preferential attachment model. In particular, we compute the first two moments of $μ$, identify the right tail of $μ$, and argue that $μ([0,\infty))\geq\tfrac{1}{2}$, a property we refer to as friendship-paradox significance. △ Less

Submitted 10 May, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Largest Eigenvalue of the Configuration Model and Breaking of Ensemble Equivalence

Authors: Pierfrancesco Dionigi, Diego Garlaschelli, Rajat Subhra Hazra, Frank den Hollander

Abstract: We analyse the largest eigenvalue of the adjacency matrix of the configuration model with large degrees, where the latter are treated as hard constraints. In particular, we compute the expectation of the largest eigenvalue for degrees that diverge as the number of vertices $n$ tends to infinity, uniformly on a scale between $1$ and $\sqrt{n}$, and show that a weak law of large numbers holds. We co… ▽ More We analyse the largest eigenvalue of the adjacency matrix of the configuration model with large degrees, where the latter are treated as hard constraints. In particular, we compute the expectation of the largest eigenvalue for degrees that diverge as the number of vertices $n$ tends to infinity, uniformly on a scale between $1$ and $\sqrt{n}$, and show that a weak law of large numbers holds. We compare with what was derived in our earlier paper "Central limit theorem for the principal eigenvalue and eigenvector of Chung-Lu random graphs" for the Chung-Lu model, which in the regime considered represents the corresponding configuration model with soft constraints, and show that the expectation is shifted down by $1$ asymptotically. This shift is a signature of breaking of ensemble equivalence between the hard and soft (also known as micro-canonical and canonical) versions of the configuration model. The latter result generalizes a previous finding in "A spectral signature of breaking of ensemble equivalence for constrained random graphs" obtained in the case when all degrees are equal. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: 21 pages, 3 figures

Continuum graph dynamics via population dynamics: well-posedness, duality and equilibria

Authors: Andreas Greven, Frank den Hollander, Anton Klimovsky, Anita Winter

Abstract: In this paper we consider stochastic processes taking values in a set of continuum graphs we call graphemes, defined as equivalence classes of sequences of vertices labelled by N embedded in an uncountable Polish space (with the cardinality of the continuum), together with an N x N connection matrix with entries 0 or 1 specifying the absence or presence of edges between pairs of vertices. In parti… ▽ More In this paper we consider stochastic processes taking values in a set of continuum graphs we call graphemes, defined as equivalence classes of sequences of vertices labelled by N embedded in an uncountable Polish space (with the cardinality of the continuum), together with an N x N connection matrix with entries 0 or 1 specifying the absence or presence of edges between pairs of vertices. In particular, we construct a Markov process on a Polish state space G of graphemes suitable to describe the time-space path of countable graphs. The class of dynamics we propose arises by specifying simple rules for the evolution of finite graphs and passing to the limit of infinite graphs. The evolution of graphemes is characterised by well-posed martingale problems, and leads to strong Markov processes with the Feller property. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: 78 pages, 2 figures

MSC Class: 05C80; 60J68; 60J70; 92D25

Quenched large deviations in renewal theory

Authors: Frank den Hollander, Marco Zamparo

Abstract: In this paper we introduce and study renewal-reward processes in random environments where each renewal involves a reward taking values in a Banach space. We derive quenched large deviation principles and identify the associated rate functions in terms of variational formulas involving correctors. We illustrate the theory with three examples: compound Poisson processes in random environments, pinn… ▽ More In this paper we introduce and study renewal-reward processes in random environments where each renewal involves a reward taking values in a Banach space. We derive quenched large deviation principles and identify the associated rate functions in terms of variational formulas involving correctors. We illustrate the theory with three examples: compound Poisson processes in random environments, pinning of polymers at interfaces with disorder, and returns of Markov chains in dynamic random environments. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: Submitted to the special issue of Stochastic Processes and their Applications in honor of Francis Comets

The grapheme-valued Wright-Fisher diffusion with mutation

Authors: Andreas Greven, Frank den Hollander, Anton Klimovsky, Anita Winter

Abstract: In [Athreya, den Hollander, Röllin; 2021, arXiv:1908.06241] models from population genetics were used to define stochastic dynamics in the space of graphons arising as continuum limits of dense graphs. In the present paper we exhibit an example of a simple neutral population genetics model for which this dynamics is a Markovian diffusion that can be characterised as the solution of a martingale pr… ▽ More In [Athreya, den Hollander, Röllin; 2021, arXiv:1908.06241] models from population genetics were used to define stochastic dynamics in the space of graphons arising as continuum limits of dense graphs. In the present paper we exhibit an example of a simple neutral population genetics model for which this dynamics is a Markovian diffusion that can be characterised as the solution of a martingale problem. In particular, we consider a Markov chain in the space of finite graphs that resembles a Moran model with resampling and mutation. We encode the finite graphs as graphemes, which can be represented as a triple consisting of a vertex set, an adjacency matrix and a sampling measure. We equip the space of graphons with convergence of sample subgraph densities and show that the grapheme-valued Markov chain converges to a grapheme-valued diffusion as the number of vertices goes to infinity. We show that the grapheme-valued diffusion has a stationary distribution that is linked to the Poisson-Dirichlet distribution. In a companion paper [Greven, den Hollander, Klimovsky, Winter; 2023], we build up a general theory for obtaining grapheme-valued diffusions via genealogies of models in population genetics. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 25 pages

MSC Class: 05C80; 60J68; 60J70; 92D25

The annealed parabolic Anderson model on a regular tree

Authors: Frank den Hollander, Daoyi Wang

Abstract: We study the total mass of the solution to the parabolic Anderson model on a regular tree with an i.i.d. random potential whose marginal distribution is double-exponential. In earlier work we identified two terms in the asymptotic expansion for large time of the total mass under the quenched law, i.e., conditional on the realisation of the random potential. In the present paper we do the same for… ▽ More We study the total mass of the solution to the parabolic Anderson model on a regular tree with an i.i.d. random potential whose marginal distribution is double-exponential. In earlier work we identified two terms in the asymptotic expansion for large time of the total mass under the quenched law, i.e., conditional on the realisation of the random potential. In the present paper we do the same for the annealed law, i.e., averaged over the random potential. It turns out that the annealed expansion differs from the quenched expansion. The derivation of the annealed expansion is based on a new approach to control the local times of the random walk appearing in the Feynman-Kac formula for the total mass. In particular, we condition on the backbone to infinity of the random walk, truncate and periodise the infinite tree relative to the backbone to obtain a random walk on a finite subtree with a specific boundary condition, employ the large deviation principle for the empirical distribution of Markov renewal processes on finite graphs, and afterwards let the truncation level tend to infinity to obtain an asymptotically sharp asymptotic expansion. △ Less

Submitted 8 July, 2023; originally announced July 2023.

Large deviation principle for the norm of the Laplacian matrix of inhomogeneous Erdős-Rényi random graphs

Authors: Rajat Subhra Hazra, Frank den Hollander, Maarten Markering

Abstract: We consider an inhomogeneous Erdős-Rényi random graph $G_N$ with vertex set $[N] = \{1,\dots,N\}$ for which the pair of vertices $i,j \in [N]$, $i\neq j$, is connected by an edge with probability $r_N(\tfrac{i}{N},\tfrac{j}{N})$, independently of other pairs of vertices. Here, $r_N\colon\,[0,1]^2 \to (0,1)$ is a symmetric function that plays the role of a reference graphon. Let $λ_N$ be the maxima… ▽ More We consider an inhomogeneous Erdős-Rényi random graph $G_N$ with vertex set $[N] = \{1,\dots,N\}$ for which the pair of vertices $i,j \in [N]$, $i\neq j$, is connected by an edge with probability $r_N(\tfrac{i}{N},\tfrac{j}{N})$, independently of other pairs of vertices. Here, $r_N\colon\,[0,1]^2 \to (0,1)$ is a symmetric function that plays the role of a reference graphon. Let $λ_N$ be the maximal eigenvalue of the Laplacian matrix of $G_N$. We show that if $\lim_{N\to\infty} \|r_N-r\|_\infty = 0$ for some limiting graphon $r\colon\,[0,1]^2 \to (0,1)$, then $λ_N/N$ satisfies a downward LDP with rate $\binom{N}{2}$ and an upward LDP with rate $N$. We identify the associated rate functions $ψ_r$ and $\widehatψ_r$, and derive their basic properties. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 22 pages

Droplet dynamics in a two-dimensional rarefied gas under Kawasaki dynamics

Authors: Simone Baldassarri, Alexandre Gaudillière, Frank den Hollander, Francesca Romana Nardi, Enzo Olivieri, Elisabetta Scoppola

Abstract: This is the second in a series of three papers in which we study a lattice gas subject to Kawasaki conservative dynamics at inverse temperature $β>0$ in a large finite box $Λ_β\subset\mathbb Z^2$ whose volume depends on $β$. Each pair of neighbouring particles has a negative binding energy $-U<0$, while each particle has a positive activation energy $Δ>0$. The initial configuration is drawn from t… ▽ More This is the second in a series of three papers in which we study a lattice gas subject to Kawasaki conservative dynamics at inverse temperature $β>0$ in a large finite box $Λ_β\subset\mathbb Z^2$ whose volume depends on $β$. Each pair of neighbouring particles has a negative binding energy $-U<0$, while each particle has a positive activation energy $Δ>0$. The initial configuration is drawn from the grand-canonical ensemble restricted to the set of configurations where all the droplets are subcritical. Our goal is to describe, in the metastable regime $Δ\in(U,2U)$ and in the limit as $β\to\infty$, how and when the system nucleates. In the first paper we showed that subcritical droplets behave as quasi-random walks. In the present paper we use the results in the first paper to analyse how subcritical droplets form and dissolve on multiple space-time scales when the volume is moderately large, i.e., $|Λ_β|=\mathrm e^{Θβ}$ with $Δ<Θ<2Δ-U$. In the third paper we consider the setting where the volume is very large, namely, $|Λ_β|=\mathrm e^{Θβ}$ with $Δ<Θ<Γ-(2Δ-U)$, where $Γ$ is the energy of the critical droplet in the local model with fixed volume, and use the results in the first two papers to identify the nucleation time. We will see that in a very large volume critical droplets appear more or less independently in boxes of moderate volume, a phenomenon referred to as homogeneous nucleation. Since Kawasaki dynamics is conservative, i.e., particles are preserved, we need to control non-local effects in the way droplets are formed and dissolved. This is done via a deductive approach: the tube of typical trajectories leading to nucleation is described via a series of events on which the evolution of the gas consists of droplets wandering around on multiple space-time scales. △ Less

Submitted 4 June, 2024; v1 submitted 27 April, 2023; originally announced April 2023.

Comments: 46 pages, 11 figures

MSC Class: 60K35; 82C26; 82C27

Spatial populations with seed-bank: finite-systems scheme

Authors: Andreas Greven, Frank den Hollander

Abstract: We consider a system of interacting Fisher-Wright diffusions with seed-bank. Individuals carry type one of two types, live in colonies, and are subject to resampling and migration as long as they are active. Each colony has a structured seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active again. As geographic space label… ▽ More We consider a system of interacting Fisher-Wright diffusions with seed-bank. Individuals carry type one of two types, live in colonies, and are subject to resampling and migration as long as they are active. Each colony has a structured seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active again. As geographic space labelling the colonies we consider a countable Abelian group endowed with the discrete topology. In earlier work we showed that the system has a one-parameter family of equilibria controlled by the relative density of the two types. Moreover, these equilibria exhibit a dichotomy of coexistence (= locally multi-type equilibrium) versus clustering (= locally mono-type equilibrium). We identified the parameter regimes for which these two phases occur, and found that these regimes are different when the mean wake-up time of a dormant individual is finite or infinite. The goal of the present paper is to establish the finite-systems scheme, i.e., identify how a finite truncation of the system (both in the geographic space and in the seed-bank) behaves as both the time and the truncation level tend to infinity, properly tuned together. If the wake-up time has finite mean, then there is a single universality class for the scaling limit. On the other hand, if the wake-up time has infinite mean, then there are two universality classes depending on how fast the truncation level of the seed-bank grows compared to the truncation level of the geographic space. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 56 pages, 7 figures

MSC Class: 60J70; 60K35; 92D325

Metastability of Glauber dynamics with inhomogeneous coupling disorder

Authors: Anton Bovier, Frank den Hollander, Saeda Marello, Elena Pulvirenti, Martin Slowik

Abstract: We introduce a general class of mean-field-like spin systems with random couplings that comprises both the Ising model on inhomogeneous dense random graphs and the randomly diluted Hopfield model. We are interested in quantitative estimates of metastability in large volumes at fixed temperatures when these systems evolve according to a Glauber dynamics, i.e.\ where spins flip with Metropolis trans… ▽ More We introduce a general class of mean-field-like spin systems with random couplings that comprises both the Ising model on inhomogeneous dense random graphs and the randomly diluted Hopfield model. We are interested in quantitative estimates of metastability in large volumes at fixed temperatures when these systems evolve according to a Glauber dynamics, i.e.\ where spins flip with Metropolis transition probabilities at inverse temperature $β$. We identify conditions ensuring that with high probability the system behaves like the corresponding system where the random couplings are replaced by their averages. More precisely, we prove that the metastability of the former system is implied with high probability by the metastability of the latter. Moreover, we consider relevant metastable hitting times of the two systems and find the asymptotic tail behaviour and the moments of their ratio. This work provides an extension of the results known for the Ising model on the Erdős--R{é}nyi random graph. The proofs use the potential-theoretic approach to metastability in combination with concentration inequalities. △ Less

Submitted 9 July, 2024; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 30 pages

MSC Class: 60K35; 60K37; 82B20; 82B44; 82C44

Graphon-valued processes with vertex-level fluctuations

Authors: Peter Braunsteins, Frank den Hollander, Michel Mandjes

Abstract: We consider a class of graph-valued stochastic processes in which each vertex has a type that fluctuates randomly over time. Collectively, the paths of the vertex types up to a given time determine the probabilities that the edges are active or inactive at that time. Our focus is on the evolution of the associated empirical graphon in the limit as the number of vertices tends to infinity, in the s… ▽ More We consider a class of graph-valued stochastic processes in which each vertex has a type that fluctuates randomly over time. Collectively, the paths of the vertex types up to a given time determine the probabilities that the edges are active or inactive at that time. Our focus is on the evolution of the associated empirical graphon in the limit as the number of vertices tends to infinity, in the setting where fluctuations in the graph-valued process are more likely to be caused by fluctuations in the vertex types than by fluctuations in the states of the edges given these types. We derive both sample-path large deviation principles and convergence of stochastic processes. We demonstrate the flexibility of our approach by treating a class of stochastic processes where the edge probabilities depend not only on the fluctuations in the vertex types but also on the state of the graph itself. △ Less

Submitted 4 September, 2022; originally announced September 2022.

Comments: 41 pages, 4 figures

MSC Class: 05C80; 60C05; 60F10

Discordant edges for the voter model on regular random graphs

Authors: Luca Avena, Rangel Baldasso, Rajat Subhra Hazra, Frank den Hollander, Matteo Quattropani

Abstract: We consider the two-opinion voter model on a regular random graph with n vertices and degree $d \geq 3$. It is known that consensus is reached on time scale n and that on this time scale the volume of the set of vertices with one opinion evolves as a Fisher-Wright diffusion. We are interested in the evolution of the number of discordant edges (i.e., edges linking vertices with different opinions),… ▽ More We consider the two-opinion voter model on a regular random graph with n vertices and degree $d \geq 3$. It is known that consensus is reached on time scale n and that on this time scale the volume of the set of vertices with one opinion evolves as a Fisher-Wright diffusion. We are interested in the evolution of the number of discordant edges (i.e., edges linking vertices with different opinions), which can be thought as the perimeter of the set of vertices with one opinion, and is the key observable capturing how consensus is reached. We show that if initially the two opinions are drawn independently from a Bernoulli distribution with parameter $u \in (0, 1)$, then on time scale 1 the fraction of discordant edges decreases and stabilises to a value that depends on d and u, and is related to the meeting time of two random walks on an infinite tree of degree d starting from two neighbouring vertices. Moreover, we show that on time scale n the fraction of discordant edges moves away from the constant plateau and converges to zero in an exponential fashion. Our proofs exploit the classical dual system of coalescing random walks and use ideas from Cooper et al. (2010) built on the so-called First Visit Time Lemma. We further introduce a novel technique to derive concentration properties from weak-dependence of coalescing random walks on moderate time scales. △ Less

Submitted 10 April, 2024; v1 submitted 2 September, 2022; originally announced September 2022.

Comments: Manuscript matching published version

Report number: ALEA, 21, 431-464 MSC Class: 05C81; 60K35

Central limit theorem for the principal eigenvalue and eigenvector of Chung-Lu random graphs

Authors: Pierfrancesco Dionigi, Diego Garlaschelli, Rajat Subhra Hazra, Frank den Hollander, Michel Mandjes

Abstract: A Chung-Lu random graph is an inhomogeneous Erdős-Rényi random graph in which vertices are assigned average degrees, and pairs of vertices are connected by an edge with a probability that is proportional to the product of their average degrees, independently for different edges. We derive a central limit theorem for the principal eigenvalue and the components of the principal eigenvector of the ad… ▽ More A Chung-Lu random graph is an inhomogeneous Erdős-Rényi random graph in which vertices are assigned average degrees, and pairs of vertices are connected by an edge with a probability that is proportional to the product of their average degrees, independently for different edges. We derive a central limit theorem for the principal eigenvalue and the components of the principal eigenvector of the adjacency matrix of a Chung-Lu random graph. Our derivation requires certain assumptions on the average degrees that guarantee connectivity, sparsity and bounded inhomogeneity of the graph. △ Less

Submitted 7 July, 2022; originally announced July 2022.

MSC Class: 60B20; 60C05; 60K35

Sparse random graphs with many triangles

Authors: Suman Chakraborty, Remco van der Hofstad, Frank den Hollander

Abstract: In this paper we consider the Erdős-Rényi random graph in the sparse regime in the limit as the number of vertices $n$ tends to infinity. We are interested in what this graph looks like when it contains many triangles, in two settings. First, we derive asymptotically sharp bounds on the probability that the graph contains a large number of triangles. We show that conditionally on this event, with… ▽ More In this paper we consider the Erdős-Rényi random graph in the sparse regime in the limit as the number of vertices $n$ tends to infinity. We are interested in what this graph looks like when it contains many triangles, in two settings. First, we derive asymptotically sharp bounds on the probability that the graph contains a large number of triangles. We show that conditionally on this event, with high probability the graph contains an almost complete subgraph, i.e., the triangles form a near-clique, and has the same local limit as the original Erdős-Rényi random graph. Second, we derive asymptotically sharp bounds on the probability that the graph contains a large number of vertices that are part of a triangle. If order $n$ vertices are in triangles, then the local limit (provided it exists) is different from that of the Erdős-Rényi random graph. Our results shed light on the challenges that arise in the description of real-world networks, which often are sparse, yet highly clustered, and on exponential random graphs, which often are used to model such networks. △ Less

Submitted 17 February, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: 31 pages, updated references

MSC Class: 05C80; 60F10

Spatial Populations with seed-bank: renormalisation on the hierarchical group

Authors: Andreas Greven, Frank den Hollander, Margriet Oomen

Abstract: We consider a system of interacting diffusions labeled by a geographic space that is given by the hierarchical group $Ω_N$ of order $N\in\mathbb{N}$. Individuals live in colonies and are subject to resampling and migration as long as they are active. Each colony has a seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active… ▽ More We consider a system of interacting diffusions labeled by a geographic space that is given by the hierarchical group $Ω_N$ of order $N\in\mathbb{N}$. Individuals live in colonies and are subject to resampling and migration as long as they are active. Each colony has a seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active again. The migration kernel has a hierarchical structure: individuals hop between colonies at a rate that depends on the hierarchical distance between the colonies. The seed-bank has a layered structure: when individuals become dormant they acquire a colour that determines the rate at which they become active again. The latter allows us to model seed-banks whose wake-up times have a fat tail. We analyse a system of coupled stochastic differential equations that describes the population in the large-colony-size limit. For fixed $N\in\mathbb{N}$, the system exhibits a dichotomy between coexistence and clustering. We identify the range of parameters controlling the migration and the seed-bank for which clustering prevails. We carry out a multi-scale renormalisation analysis in the hierarchical mean-field limit $N\to\infty$. We show that block averages on hierarchical space-time scale $k \in \mathbb{N}$ perform a diffusion with a renormalised diffusion function that depends on $k$. In the clustering regime, after an appropriate scaling with $k$, this diffusion function converges to the Fisher-Wright diffusion function as $k\to\infty$, irrespective of the diffusion function controlling the resampling. For several subclasses of parameters we identify the speed at which the scaled renormalised diffusion function converges to the Fisher-Wright diffusion as $k\to\infty$. We show that the seed-bank reduces the speed compared to the model without seed-bank. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: 167 pages

MSC Class: Primary 60J70; 60K35; Secondary 92D25

Spatially Inhomogeneous Populations with Seed-banks: II. Clustering Regime

Authors: Frank den Hollander, Shubhamoy Nandan

Abstract: We consider a spatial version of the classical Moran model with seed-banks where the constituent populations have finite sizes. Individuals live in colonies labelled by $\mathbb{Z}^d$, $d\geq 1$, playing the role of a geographic space, carry one of the two $types$: $\heartsuit$ or $\spadesuit$, and change type via $resampling$ as long as they are $active$. Each colony contains a seed-bank into whi… ▽ More We consider a spatial version of the classical Moran model with seed-banks where the constituent populations have finite sizes. Individuals live in colonies labelled by $\mathbb{Z}^d$, $d\geq 1$, playing the role of a geographic space, carry one of the two $types$: $\heartsuit$ or $\spadesuit$, and change type via $resampling$ as long as they are $active$. Each colony contains a seed-bank into which individuals can enter to become $dormant$, suspending their resampling until they exit the seed-bank and become active again. Individuals resample not only from their own colony, but also from other colonies according to a symmetric random walk transition kernel. The latter is referred to as $migration$. The sizes of the populations vary across colonies and remain constant in time. It was shown in Hollander and Nandan (2021) that the system is well-defined, admits a family of equilibria parametrized by the initial density of type $\heartsuit$, and exhibits a dichotomy between $clustering$ (mono-type equilibrium) and $coexistence$ (multi-type equilibrium). This dichotomy is determined by a clustering criterion that is given in terms of a dual of the system, which consists of a system of $interacting$ coalescing random walks. In this paper we provide an alternative clustering criterion, given in terms of an auxiliary dual that is simpler than the original dual, and identify a range of parameters for which the criterion is met, which we refer to as the $clustering$ $regime$. It turns out that if the sizes of the active populations are non-clumping (i.e., do not take arbitrarily large values in finite regions of the geographic space) and the relative strengths of the seed-banks (i.e., the ratio of the sizes of dormant and active population in each colony) are bounded uniformly over the geographic space, then clustering prevails if and only if the (symmetrised) migration kernel is recurrent. △ Less

Submitted 20 December, 2021; v1 submitted 31 July, 2021; originally announced August 2021.

Comments: Version 2, added intuitive explainations, remarks etc. The old assumption 2.9 on migration kernel is now weakened to Assumption 2.11

MSC Class: Primary 60K35; Secondary 92D25

Switching interacting particle systems: scaling limits, uphill diffusion and boundary layer

Authors: Simone Floreani, Cristian Giardinà, Frank den Hollander, Shubhamoy Nandan, Frank Redig

Abstract: In this paper we consider three classes of interacting particle systems on $\mathbb Z$: independent random walks, the exclusion process, and the inclusion process. We allow particles to switch their jump rate (the rate identifies the type of particle) between $1$ (fast particles) and $ε\in[0,1]$ (slow particles). The switch between the two jump rates happens at rate $γ\in(0,\infty)$. In the exclus… ▽ More In this paper we consider three classes of interacting particle systems on $\mathbb Z$: independent random walks, the exclusion process, and the inclusion process. We allow particles to switch their jump rate (the rate identifies the type of particle) between $1$ (fast particles) and $ε\in[0,1]$ (slow particles). The switch between the two jump rates happens at rate $γ\in(0,\infty)$. In the exclusion process, the interaction is such that each site can be occupied by at most one particle of each type. In the inclusion process, the interaction takes places between particles of the same type at different sites and between particles of different type at the same site. We derive the macroscopic limit equations for the three systems, obtained after scaling space by $N^{-1}$, time by $N^2$, the switching rate by $N^{-2}$, and letting $N\to\infty$. The limit equations for the macroscopic densities associated to the fast and slow particles is the well-studied double diffusivity model. This system of reaction-diffusion equations was introduced to model polycrystal diffusion and dislocation pipe diffusion, with the goal to overcome the limitations imposed by Fick's law. In order to investigate the microscopic out-of-equilibrium properties, we analyse the system on $[N]=\{1,\ldots,N\}$, adding boundary reservoirs at sites $1$ and $N$ of fast and slow particles, respectively. Inside $[N]$ particles move as before, but now particles are injected and absorbed at sites $1$ and $N$ with prescribed rates that depend on the particle type. We compute the steady-state density profile and the steady-state current. It turns out that uphill diffusion is possible, i.e., the total flow can be in the direction of increasing total density. This phenomenon, which cannot occur in a single-type particle system, is a violation of Fick's law made possible by the switching between types. △ Less

Submitted 14 July, 2021; originally announced July 2021.

MSC Class: 60J70; 60K35 (Primary); 82C26; 92D25 (Secondary)

Metastability for Glauber dynamics on the complete graph with coupling disorder

Authors: Anton Bovier, Frank den Hollander, Saeda Marello

Abstract: Consider the complete graph on $n$ vertices. To each vertex assign an Ising spin that can take the values $-1$ or $+1$. Each spin $i \in [n]=\{1,2,\dots, n\}$ interacts with a magnetic field $h \in [0,\infty)$, while each pair of spins $i,j \in [n]$ interact with each other at coupling strength $n^{-1} J(i)J(j)$, where $J=(J(i))_{i \in [n]}$ are i.i.d. non-negative random variables drawn from a pr… ▽ More Consider the complete graph on $n$ vertices. To each vertex assign an Ising spin that can take the values $-1$ or $+1$. Each spin $i \in [n]=\{1,2,\dots, n\}$ interacts with a magnetic field $h \in [0,\infty)$, while each pair of spins $i,j \in [n]$ interact with each other at coupling strength $n^{-1} J(i)J(j)$, where $J=(J(i))_{i \in [n]}$ are i.i.d. non-negative random variables drawn from a probability distribution with finite support. Spins flip according to a Metropolis dynamics at inverse temperature $β\in (0,\infty)$. We show that there are critical thresholds $β_c$ and $h_c(β)$ such that, in the limit as $n\to\infty$, the system exhibits metastable behaviour if and only if $β\in (β_c, \infty)$ and $h \in [0,h_c(β))$. Our main result is a sharp asymptotics, up to a multiplicative error $1+o_n(1)$, of the average crossover time from any metastable state to the set of states with lower free energy. We use standard techniques of the potential-theoretic approach to metastability. The leading order term in the asymptotics does not depend on the realisation of $J$, while the correction terms do. The leading order of the correction term is $\sqrt{n}$ times a centred Gaussian random variable with a complicated variance depending on $β,h$, on the law of $J$ and on the metastable state. The critical thresholds $β_c$ and $h_c(β)$ depend on the law of $J$, and so does the number of metastable states. We derive an explicit formula for $β_c$ and identify some properties of $β\mapsto h_c(β)$. Interestingly, the latter is not necessarily monotone, meaning that the metastable crossover may be re-entrant. △ Less

Submitted 27 April, 2022; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 41 pages, 6 figures

MSC Class: 60K35; 60K37; 82B20; 82B44; 82C44

Journal ref: Commun. Math. Phys. 392 (2022), 307-345

Breaking of ensemble equivalence for dense random graphs under a single constraint

Authors: Frank den Hollander, Maarten Markering

Abstract: Two ensembles are often used to model random graphs subject to constraints: the microcanonical ensemble (= hard constraint) and the canonical ensemble (= soft constraint). It is said that breaking of ensemble equivalence (BEE) occurs when the specific relative entropy of the two ensembles does not vanish as the size of the graph tends to infinity. The latter means that it matters for the scaling p… ▽ More Two ensembles are often used to model random graphs subject to constraints: the microcanonical ensemble (= hard constraint) and the canonical ensemble (= soft constraint). It is said that breaking of ensemble equivalence (BEE) occurs when the specific relative entropy of the two ensembles does not vanish as the size of the graph tends to infinity. The latter means that it matters for the scaling properties of the graph whether the constraint is met for every single realisation of the graph or only holds as an ensemble average. In the literature, it was found that BEE is the rule rather than the exception for two classes: sparse random graphs when the number of constraints is of the order of the number of vertices and dense random graphs when there are two or more constraints that are frustrated. In the present paper we establish BEE for a third class: dense random graphs with a single constraint, namely, on the density of a given finite simple graph. We show that BEE occurs only in a certain range of choices for the density and the number of edges of the simple graph, which we refer to as the BEE-phase. We show that, in part of the BEE-phase, there is a gap between the scaling limits of the averages of the maximal eigenvalue of the adjacency matrix of the random graph under the two ensembles, a property that is referred to as spectral signature of BEE. Proofs are based on an analysis of the variational formula on the space of graphons for the limiting specific relative entropy derived by Den Hollander et al. (2018), in combination with an identification of the minimising graphons and replica symmetry arguments. We show that in the replica symmetric region of the BEE-phase, as the size of the graph tends to infinity, the microcanonical ensemble behaves like an Erdős-Rényi random graph, while the canonical ensemble behaves like a mixture of two Erdős-Rényi random graphs. △ Less

Submitted 4 October, 2022; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 21 pages, 3 figures. arXiv admin note: text overlap with arXiv:1807.07750

MSC Class: 05C80; 60C05; 60F10; 82B20

The Parabolic Anderson Model on a Galton-Watson tree revisited

Authors: Frank den Hollander, Daoyi Wang

Abstract: In [1] a detailed analysis was given of the large-time asymptotics of the total mass of the solution to the parabolic Anderson model on a supercritical Galton-Watson random tree with an i.i.d. random potential whose marginal distribution is double-exponential. Under the assumption that the degree distribution has bounded support, two terms in the asymptotic expansion were identified under the quen… ▽ More In [1] a detailed analysis was given of the large-time asymptotics of the total mass of the solution to the parabolic Anderson model on a supercritical Galton-Watson random tree with an i.i.d. random potential whose marginal distribution is double-exponential. Under the assumption that the degree distribution has bounded support, two terms in the asymptotic expansion were identified under the quenched law, i.e., conditional on the realisation of the random tree and the random potential. The second term contains a variational formula indicating that the solution concentrates on a subtree with minimal degree according to a computable profile. The present paper extends the analysis to degree distributions with unbounded support. We identify the weakest condition on the tail of the degree distribution under which the arguments in [1] can be pushed through. To do so we need to control the occurrence of large degrees uniformly in large subtrees of the Galton-Watson tree. △ Less

Submitted 30 March, 2021; originally announced March 2021.

MSC Class: 60H25; 82B44; 05C80

Linking the mixing times of random walks on static and dynamic random graphs

Authors: Luca Avena, Hakan Güldaş, Remco van der Hofstad, Frank den Hollander, Oliver Nagy

Abstract: This paper considers non-backtracking random walks on random graphs generated according to the configuration model. The quantity of interest is the scaling of the mixing time of the random walk as the number of vertices of the random graph tends to infinity. Subject to mild general conditions, we link two mixing times: one for a static version of the random graph, the other for a class of dynamic… ▽ More This paper considers non-backtracking random walks on random graphs generated according to the configuration model. The quantity of interest is the scaling of the mixing time of the random walk as the number of vertices of the random graph tends to infinity. Subject to mild general conditions, we link two mixing times: one for a static version of the random graph, the other for a class of dynamic versions of the random graph in which the edges are randomly rewired but the degrees are preserved. The link is provided by the probability that the random walk has not yet stepped along a previously rewired edge. We use this link to compute the scaling of the mixing time for three specific classes of random rewirings. Depending on the speed and the range of the rewiring relative to the current location of the random walk, the mixing time may exhibit no cut-off, one-sided cut-off or two-sided cut-off, a trichotomy that was also found in earlier work. Interestingly, for a class of dynamics that are `mesoscopic', i.e., non-local and non-global, we find new behaviour with six subregimes. Proofs are built on a new and flexible coupling scheme, in combination with sharp estimates on the degrees encountered by the random walk in the static and the dynamic version of the random graph. Some of these estimates require sharp control on possible short-cuts in the graph between the edges that are traversed by the random walk. △ Less

Submitted 7 February, 2021; v1 submitted 20 December, 2020; originally announced December 2020.

Comments: 38 pages, 9 figures

MSC Class: 05C81; 37A25; 60K37; 82C27

Principles of seed banks: complexity emerging from dormancy

Authors: Jay T. Lennon, Frank den Hollander, Maite Wilke-Berenguer, Jochen Blath

Abstract: Across the tree of life, populations have evolved the capacity to contend with suboptimal conditions by engaging in dormancy, whereby individuals enter a reversible state of reduced metabolic activity. The resulting seed banks are complex, storing information and imparting memory that gives rise to multi-scale structures and networks spanning collections of cells to entire ecosystems. We outline t… ▽ More Across the tree of life, populations have evolved the capacity to contend with suboptimal conditions by engaging in dormancy, whereby individuals enter a reversible state of reduced metabolic activity. The resulting seed banks are complex, storing information and imparting memory that gives rise to multi-scale structures and networks spanning collections of cells to entire ecosystems. We outline the fundamental attributes and emergent phenomena associated with dormancy and seed banks, with the vision for a unifying and mathematically based framework that can address problems in the life sciences, ranging from global change to cancer biology. △ Less

Submitted 6 May, 2021; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 60 pages, 2 figures, 6 boxes

Spatially Inhomogeneous Populations with Seed-banks: Duality, Existence and Clustering

Authors: Frank den Hollander, Shubhamoy Nandan

Abstract: We consider a system of interacting Moran models with seed-banks. Individuals live in colonies and are subject to resampling and migration as long as they are $active$. Each colony has a seed-bank into which individuals can retreat to become $dormant$, suspending their resampling and migration until they become active again. The colonies are labelled by $\mathbb{Z}^d$, $d \geq 1$, playing the role… ▽ More We consider a system of interacting Moran models with seed-banks. Individuals live in colonies and are subject to resampling and migration as long as they are $active$. Each colony has a seed-bank into which individuals can retreat to become $dormant$, suspending their resampling and migration until they become active again. The colonies are labelled by $\mathbb{Z}^d$, $d \geq 1$, playing the role of a $geographic\, space$. The sizes of the active and the dormant population are $finite$ and depend on the $location$ of the colony. Migration is driven by a random walk transition kernel. Our goal is to study the equilibrium behaviour of the system as a function of the underlying model parameters. In the present paper we show that, under mild condition on the sizes of the active population, the system is well-defined and has a dual. The dual consists of a system of $interacting$ coalescing random walks in an $inhomogeneous$ environment that switch between active and dormant. We analyse the dichotomy of $coexistence$ (= multi-type equilibria) versus $clustering$ (= mono-type equilibria), and show that clustering occurs if and only if two random walks in the dual starting from arbitrary states eventually coalesce with probability one. The presence of the seed-bank $enhances\, genetic\, diversity$. In the dual this is reflected by the presence of time lapses during which the random walks are dormant and do not move. △ Less

Submitted 6 July, 2021; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: Changed title, added references and fixed typos. A different proof of the original Theorem 3.16 (Clustering Criterion) now stated as Theorem 3.17, is presented. Added a new lemma (Lemma 3.8) to justification of well-posedness of the spatial dual

MSC Class: 60J70 (Primary) 60K35; 92D25 (Secondary)

A sample-path large deviation principle for dynamic Erdős-Rényi random graphs

Authors: Peter Braunsteins, Frank den Hollander, Michel Mandjes

Abstract: We consider a dynamic Erdős-Rényi random graph (ERRG) on $n$ vertices in which each edge switches on at rate $λ$ and switches off at rate $μ$, independently of other edges. The focus is on the analysis of the evolution of the associated empirical graphon in the limit as $n\to\infty$. Our main result is a large deviation principle (LDP) for the sample path of the empirical graphon observed until a… ▽ More We consider a dynamic Erdős-Rényi random graph (ERRG) on $n$ vertices in which each edge switches on at rate $λ$ and switches off at rate $μ$, independently of other edges. The focus is on the analysis of the evolution of the associated empirical graphon in the limit as $n\to\infty$. Our main result is a large deviation principle (LDP) for the sample path of the empirical graphon observed until a fixed time horizon. The rate is $\binom{n}{2}$, the rate function is a specific action integral on the space of graphon trajectories. We apply the LDP to identify (i) the most likely path that starting from a constant graphon creates a graphon with an atypically large density of $d$-regular subgraphs, and (ii) the mostly likely path between two given graphons. It turns out that bifurcations may occur in the solutions of associated variational problems. △ Less

Submitted 27 September, 2020; originally announced September 2020.

A spectral signature of breaking of ensemble equivalence for constrained random graphs

Authors: Pierfrancesco Dionigi, Diego Garlaschelli, Frank den Hollander, Michel Mandjes

Abstract: For random systems subject to a constraint, the microcanonical ensemble requires the constraint to be met by every realisation ("hard constraint"), while the canonical ensemble requires the constraint to be met only on average ("soft constraint"). It is known that for random graphs subject to topological constraints breaking of ensemble equivalence may occur when the size of the graph tends to inf… ▽ More For random systems subject to a constraint, the microcanonical ensemble requires the constraint to be met by every realisation ("hard constraint"), while the canonical ensemble requires the constraint to be met only on average ("soft constraint"). It is known that for random graphs subject to topological constraints breaking of ensemble equivalence may occur when the size of the graph tends to infinity, signalled by a non-vanishing specific relative entropy of the two ensembles. We investigate to what extent breaking of ensemble equivalence is manifested through the largest eigenvalue of the adjacency matrix of the graph. We consider two examples of constraints in the dense regime: (1) fix the degrees of the vertices (= the degree sequence); (2) fix the sum of the degrees of the vertices (= twice the number of edges). Example (1) imposes an extensive number of local constraints and is known to lead to breaking of ensemble equivalence. Example (2) imposes a single global constraint and is known to lead to ensemble equivalence. Our working hypothesis is that breaking of ensemble equivalence corresponds to a non-vanishing difference of the expected values of the largest eigenvalue under the two ensembles. We verify that, in the limit as the size of the graph tends to infinity, the difference between the expected values of the largest eigenvalue in the two ensembles does not vanish for (1) and vanishes for (2). A key tool in our analysis is a transfer method that uses relative entropy to determine whether probabilistic estimates can be carried over from the canonical ensemble to the microcanonical ensemble, and illustrates how breaking of ensemble equivalence may prevent this from being possible. △ Less

Submitted 1 November, 2021; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 16 pages. Updated version of the previous work: the result is generalized to broader density regimes for the considered random graph models

MSC Class: 60B20; 60C05; 60K35

Journal ref: Electron. Commun. Probab. 26: 1-15 (2021)

Large deviation principle for the maximal eigenvalue of inhomogeneous Erdős-Rényi random graphs

Authors: Arijit Chakrabarty, Rajat Subhra Hazra, Frank den Hollander, Matteo Sfragara

Abstract: We consider an inhomogeneous Erdős-Rényi random graph $G_N$ with vertex set $[N] = \{1,\dots,N\}$ for which the pair of vertices $i,j \in [N]$, $i\neq j$, is connected by an edge with probability $r(\tfrac{i}{N},\tfrac{j}{N})$, independently of other pairs of vertices. Here, $r\colon\,[0,1]^2 \to (0,1)$ is a symmetric function that plays the role of a reference graphon. Let $λ_N$ be the maximal ei… ▽ More We consider an inhomogeneous Erdős-Rényi random graph $G_N$ with vertex set $[N] = \{1,\dots,N\}$ for which the pair of vertices $i,j \in [N]$, $i\neq j$, is connected by an edge with probability $r(\tfrac{i}{N},\tfrac{j}{N})$, independently of other pairs of vertices. Here, $r\colon\,[0,1]^2 \to (0,1)$ is a symmetric function that plays the role of a reference graphon. Let $λ_N$ be the maximal eigenvalue of the adjacency matrix of $G_N$. It is known that $λ_N/N$ satisfies a large deviation principle as $N \to \infty$. The associated rate function $ψ_r$ is given by a variational formula that involves the rate function $I_r$ of a large deviation principle on graphon space. We analyse this variational formula in order to identify the properties of $ψ_r$, specially when the reference graphon is of rank 1. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: 21 pages

Phase transitions for spatially extended pinning

Authors: Francesco Caravenna, Frank den Hollander

Abstract: We consider a directed polymer of length $N$ interacting with a linear interface. The monomers carry i.i.d. random charges $(ω_i)_{i=1}^N$ taking values in $\mathbb{R}$ with mean zero and variance one. Each monomer $i$ contributes an energy $(βω_i-h)\varphi(S_i)$ to the interaction Hamiltonian, where $S_i \in \mathbb{Z}$ is the height of monomer $i$ with respect to the interface,… ▽ More We consider a directed polymer of length $N$ interacting with a linear interface. The monomers carry i.i.d. random charges $(ω_i)_{i=1}^N$ taking values in $\mathbb{R}$ with mean zero and variance one. Each monomer $i$ contributes an energy $(βω_i-h)\varphi(S_i)$ to the interaction Hamiltonian, where $S_i \in \mathbb{Z}$ is the height of monomer $i$ with respect to the interface, $\varphi: \mathbb{Z} \to [0,\infty)$ is the interaction potential, $β\in [0,\infty)$ is the inverse temperature, and $h \in \mathbb{R}$ is the charge bias parameter. The configurations of the polymer are weighted according to the Gibbs measure associated with the interaction Hamiltonian, where the reference measure is given by a Markov chain on $\mathbb{Z}$. We study both the quenched and the annealed free energy per monomer in the limit as $N\to\infty$. We show that each exhibits a phase transition along a critical curve in the $(β, h)$-plane, separating a localized phase (where the polymer stays close to the interface) from a delocalized phase (where the polymer wanders away from the interface). We derive variational formulas for the critical curves and we obtain upper and lower bounds on the quenched critical curve in terms of the annealed critical curve. In addition, for the special case where the reference measure is given by a Bessel random walk, we derive the scaling limit of the annealed free energy as $β, h \downarrow 0$ in three different regimes for the tail exponent of $\varphi$. △ Less

Submitted 8 March, 2021; v1 submitted 31 May, 2020; originally announced June 2020.

Comments: 36 pages, 4 figures. Improved presentation, final version to appear in PTRF

MSC Class: 60K37; 82B44; 82B41

Spatial populations with seed-bank: well-posedness, duality and equilibrium

Authors: Andreas Greven, Frank den Hollander, Margriet Oomen

Abstract: We consider a system of interacting Fisher-Wright diffusions with seed-bank. Individuals live in colonies and are subject to resampling and migration as long as they are active. Each colony has a structured seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active again. As geographic space labelling the colonies we consider… ▽ More We consider a system of interacting Fisher-Wright diffusions with seed-bank. Individuals live in colonies and are subject to resampling and migration as long as they are active. Each colony has a structured seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active again. As geographic space labelling the colonies we consider a countable Abelian group $\mathbb{G}$ endowed with the discrete topology. The key example of interest is the Euclidean lattice $\mathbb{G}=\mathbb{Z}^d$. Our goal is to classify the long-time behaviour of the system in terms of the underlying model parameters. In particular, we want to understand in what way the seed-bank enhances genetic diversity. We introduce three models of increasing generality, namely, individuals become dormant: (1) in the seed-bank of their colony; (2) in the seed-bank of their colony while adopting a random colour that determines their wake-up time; (3) in the seed-bank of a random colony while adopting a random colour. The extension in (2) allows us to model wake-up times with fat tails while preserving the Markov property of the evolution. For each of the three models we show that the system converges to a unique equilibrium depending on a single density parameter that is determined by the initial state, and exhibits a dichotomy of coexistence (= locally multi-type equilibrium) versus clustering (= locally mono-type equilibrium) depending on the parameters controlling the migration and the seed-bank. The dichotomy between clustering and coexistence in model 1 is determined by migration only. In models (2) and (3), when the wake-up time has infinite mean, the dichotomy is determined by both the exchange with the seed-bank and migration. It turns out that the seed-bank affects the long-time behaviour both quantitatively and qualitatively. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: 74 pages

MSC Class: 60J70 (Primary); 60K35 (Primary); 92D25 (Secondary)

The parabolic Anderson model on a Galton-Watson tree

Authors: Frank den Hollander, Wolfgang König, Renato S. dos Santos

Abstract: We study the long-time asymptotics of the total mass of the solution to the parabolic Anderson model (PAM) on a supercritical Galton-Watson random tree with bounded degrees. We identify the second-order contribution to this asymptotics in terms of a variational formula that gives information about the local structure of the region where the solution is concentrated. The analysis behind this formul… ▽ More We study the long-time asymptotics of the total mass of the solution to the parabolic Anderson model (PAM) on a supercritical Galton-Watson random tree with bounded degrees. We identify the second-order contribution to this asymptotics in terms of a variational formula that gives information about the local structure of the region where the solution is concentrated. The analysis behind this formula suggests that, under mild conditions on the model parameters, concentration takes place on a tree with minimal degree. Our approach can be applied to finite locally tree-like random graphs, in a coupled limit where both time and graph size tend to infinity. As an example, we consider the configuration model or, more precisely, the uniform simple random graph with a prescribed degree sequence. △ Less

Submitted 27 July, 2020; v1 submitted 14 January, 2020; originally announced January 2020.

Comments: 32 pages

MSC Class: 60H25; 82B44; 05C80

Wireless random-access networks with bipartite interference graphs

Authors: Sem C. Borst, Frank den Hollander, Francesca R. Nardi, Matteo Sfragara

Abstract: We consider random-access networks where nodes represent servers with a queue and can be either active or inactive. A node deactivates at unit rate, while it activates at a rate that depends on its queue length, provided none of its neighbors is active. We consider arbitrary bipartite graphs in the limit as the initial queue lengths become large and identify the transition time between the two sta… ▽ More We consider random-access networks where nodes represent servers with a queue and can be either active or inactive. A node deactivates at unit rate, while it activates at a rate that depends on its queue length, provided none of its neighbors is active. We consider arbitrary bipartite graphs in the limit as the initial queue lengths become large and identify the transition time between the two states where one half of the network is active and the other half is inactive. The transition path is decomposed into a succession of transitions on complete bipartite subgraphs. We formulate a randomized greedy algorithm that takes the graph as input and gives as output the set of transition paths the network is most likely to follow. Along each path we determine the mean transition time and its law on the scale of its mean. Depending on the activation rates, we identify three regimes of behavior. △ Less

Submitted 12 September, 2023; v1 submitted 9 January, 2020; originally announced January 2020.

Comments: 42 pages

Crossover times in bipartite networks with activity constraints and time-varying switching rates

Authors: Sem Borst, Frank den Hollander, Francesca Nardi, Siamak Taati

Abstract: In this paper we study the performance of a bipartite network in which customers arrive at the nodes of the network, but not all nodes are able to serve their customers at all times. Each node can be either active or inactive, and two nodes connected by a bond cannot be active simultaneously. This situation arises in wireless random-access networks where, due to destructive interference, stations… ▽ More In this paper we study the performance of a bipartite network in which customers arrive at the nodes of the network, but not all nodes are able to serve their customers at all times. Each node can be either active or inactive, and two nodes connected by a bond cannot be active simultaneously. This situation arises in wireless random-access networks where, due to destructive interference, stations that are close to each other cannot use the same frequency band. We consider a model where the network is bipartite, the active nodes switch themselves off at rate 1, and the inactive nodes switch themselves on at a rate that depends on time and on which half of the bipartite network they are in. An inactive node cannot become active when one of the nodes it is connected to by a bond is active. The switching protocol allows the nodes to share activity among each other. In the limit as the activation rate becomes large, we compute the crossover time between the two states where one half of the network is active and the other half is inactive. This allows us to assess the overall activity of the network depending on the switching protocol. Our results make use of the metastability analysis for hard-core interacting particle models on finite bipartite graphs derived in an earlier paper. They are valid for a large class of bipartite networks, subject to certain assumptions. Proofs rely on a comparison with switching protocols that are not time-varying, through coupling techniques. △ Less

Submitted 12 February, 2022; v1 submitted 30 December, 2019; originally announced December 2019.

Comments: 32 pages, 2 figure

MSC Class: 60K25; 60K30; 60K35; 90B15; 90B18

Glauber dynamics on the Erdős-Rényi random graph

Authors: Frank den Hollander, Oliver Jovanovski

Abstract: We investigate the effect of disorder on the Curie-Weiss model with Glauber dynamics. In particular, we study metastability for spin-flip dynamics on the Erdős-Rényi random graph $ER_n(p)$ with $n$ vertices and with edge retention probability $p \in (0,1)$. Each vertex carries an Ising spin that can take the values $-1$ or $+1$. Single spins interact with an external magnetic field… ▽ More We investigate the effect of disorder on the Curie-Weiss model with Glauber dynamics. In particular, we study metastability for spin-flip dynamics on the Erdős-Rényi random graph $ER_n(p)$ with $n$ vertices and with edge retention probability $p \in (0,1)$. Each vertex carries an Ising spin that can take the values $-1$ or $+1$. Single spins interact with an external magnetic field $h \in (0,\infty)$, while pairs of spins at vertices connected by an edge interact with each other with ferromagnetic interaction strength $1/n$. Spins flip according to a Metropolis dynamics at inverse temperature $β$. The standard Curie-Weiss model corresponds to the case $p=1$, because $ER_n(1) = K_n$ is the complete graph on $n$ vertices. For $β>β_c$ and $h \in (0,p χ(βp))$ the system exhibits \emph{metastable behaviour} in the limit as $n\to\infty$, where $β_c=1/p$ is the \emph{critical inverse temperature} and $χ$ is a certain \emph{threshold function} satisfying $\lim_{λ\to\infty} χ(λ) =1$ and $\lim_{λ\downarrow 1} χ(λ)=0$. We compute the average crossover time from the \emph{metastable set} (with magnetization corresponding to the `minus-phase') to the \emph{stable set} (with magnetization corresponding to the `plus-phase'). We show that the average crossover time grows exponentially fast with $n$, with an exponent that is the same as for the Curie-Weiss model with external magnetic field $h$ and with ferromagnetic interaction strength $p/n$. We show that the correction term to the exponential asymptotics is a multiplicative error term that is \emph{at most polynomial} in $n$. For the complete graph $K_n$ the correction term is known to be a multiplicative constant. △ Less

Submitted 22 July, 2020; v1 submitted 22 December, 2019; originally announced December 2019.

Graphon-valued stochastic processes from population genetics

Authors: Siva Athreya, Frank den Hollander, Adrian Röllin

Abstract: The goal of this paper is to develop a theory of graphon-valued stochastic processes, and to construct and analyse a natural class of such processes arising from population genetics. We consider finite populations where individuals change type according to Wright-Fisher resampling. At any time, each pair of individuals is linked by an edge with a probability that is given by a type-connection matr… ▽ More The goal of this paper is to develop a theory of graphon-valued stochastic processes, and to construct and analyse a natural class of such processes arising from population genetics. We consider finite populations where individuals change type according to Wright-Fisher resampling. At any time, each pair of individuals is linked by an edge with a probability that is given by a type-connection matrix, whose entries depend on the current empirical type distribution of the entire population via a fitness function. We show that, in the large-population-size limit and with an appropriate scaling of time, the evolution of the associated adjacency matrix converges to a random process in the space of graphons, driven by the type-connection matrix and the underlying Wright-Fisher diffusion on the multi-type simplex. In the limit as the number of types tends to infinity, the limiting process is driven by the type-connection kernel and the underlying Fleming-Viot diffusion. △ Less

Submitted 17 August, 2019; originally announced August 2019.

MSC Class: 05C80; 60J68; 60K35

The Widom-Rowlinson model: Mesoscopic fluctuations for the critical droplet

Authors: Frank den Hollander, Sabine Jansen, Roman Kotecký, Elena Pulvirenti

Abstract: We study the critical droplet for a close-to-equilibrium Widom-Rowlinson model of interacting particles in the two-dimensional continuum at low temperatures. The critical droplet is the set of macroscopic states that correspond to saddle points for the passage from a low-density vapour to a high-density liquid. We first show that the critical droplet is close to a disc of a certain deterministic r… ▽ More We study the critical droplet for a close-to-equilibrium Widom-Rowlinson model of interacting particles in the two-dimensional continuum at low temperatures. The critical droplet is the set of macroscopic states that correspond to saddle points for the passage from a low-density vapour to a high-density liquid. We first show that the critical droplet is close to a disc of a certain deterministic radius. After that we analyse the mesoscopic fluctuations of the surface of the critical droplet. These are built on microscopic fluctuations of the particles in the boundary layer, and provide a rigorous foundation for what in physics literature is known as capillary waves. Our results represent the first detailed analysis of the surface fluctuations down to mesoscopic and microscopic precision. As such they constitute a fundamental step in the study of phase separation in continuum interacting particle systems from the perspective of stochastic geometry. At the same time, they serve as a basis for the study of a non-equilibrium version of the Widom-Rowlinson model, to be analysed elsewhere, where they lead to a correction term in the Arrhenius formula for the average vapour-liquid crossover time. △ Less

Submitted 2 July, 2019; v1 submitted 30 June, 2019; originally announced July 2019.

Comments: 64 pages, 10 figures

MSC Class: 60J45; 60J60; 60K35; 82C21; 82C22; 82C27

Random walk in cooling random environment: recurrence versus transience and mixed fluctuations

Authors: Luca Avena, Yuki Chino, Conrado da Costa, Frank den Hollander

Abstract: This is the third in a series of papers in which we consider one-dimensional Random Walk in Cooling Random Environment (RWCRE). The latter is obtained by starting from one-dimensional Random Walk in Random Environment (RWRE) and resampling the environment along a sequence of deterministic times, called refreshing times. In the present paper we explore two questions for general refreshing times. Fi… ▽ More This is the third in a series of papers in which we consider one-dimensional Random Walk in Cooling Random Environment (RWCRE). The latter is obtained by starting from one-dimensional Random Walk in Random Environment (RWRE) and resampling the environment along a sequence of deterministic times, called refreshing times. In the present paper we explore two questions for general refreshing times. First, we investigate how the recurrence versus transience criterion known for RWRE changes for RWCRE. Second, we explore the fluctuations for RWCRE when RWRE is either recurrent or satisfies a classical central limit theorem. We show that the answer depends in a delicate way on the choice of the refreshing times. An overarching goal of our paper is to investigate how the behaviour of a random process with a rich correlation structure can be affected by resettings. △ Less

Submitted 24 July, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: 50 pages, 4 figures

MSC Class: 60F05; 60G50; 60K37

Spectra of Adjacency and Laplacian Matrices of Inhomogeneous Erdős-Rényi Random Graphs

Authors: Arijit Chakrabarty, Rajat Subhra Hazra, Frank den Hollander, Matteo Sfragara

Abstract: Inhomogeneous Erdős-Rényi random graphs $\mathbb G_N$ on $N$ vertices in the non-dense regime are considered in this paper. The edge between the pair of vertices $\{i,j\}$ is retained with probability $\varepsilon_N\,f(\frac{i}{N},\frac{j}{N})$, $1 \leq i \neq j \leq N$, independently of other edges, where $f\colon\,[0,1] \times [0,1] \to [0,\infty)$ is a continuous function such that… ▽ More Inhomogeneous Erdős-Rényi random graphs $\mathbb G_N$ on $N$ vertices in the non-dense regime are considered in this paper. The edge between the pair of vertices $\{i,j\}$ is retained with probability $\varepsilon_N\,f(\frac{i}{N},\frac{j}{N})$, $1 \leq i \neq j \leq N$, independently of other edges, where $f\colon\,[0,1] \times [0,1] \to [0,\infty)$ is a continuous function such that $f(x,y)=f(y,x)$ for all $x,y \in [0,1]$. We study the empirical distribution of both the adjacency matrix $A_N$ and the Laplacian matrix $Δ_N$ associated with $\mathbb G_N$ in the limit as $N \to \infty$ when $\lim_{N\to\infty} \varepsilon_N = 0$ and $\lim_{N\to\infty} N\varepsilon_N = \infty$. In particular, it is shown that the empirical spectral distributions of $A_N$ and $Δ_N$, after appropriate scaling and centering, converge to deterministic limits weakly in probability. For the special case where $f(x,y) = r(x)r(y)$ with $r\colon\,[0,1] \to [0,\infty)$ a continuous function, we give an explicit characterization of the limiting distributions. Furthermore, applications of the results to constrained random graphs, Chung-Lu random graphs and social networks are shown. △ Less

Submitted 15 October, 2019; v1 submitted 26 July, 2018; originally announced July 2018.

Comments: Revised version. To appear in Random Matrices: theory and applications

Breaking of ensemble equivalence for perturbed Erdős-Rényi random graphs

Authors: F. den Hollander, M. Mandjes, A. Roccaverde, N. J. Starreveld

Abstract: In [18] we analysed a simple undirected random graph subject to constraints on the total number of edges and the total number of triangles. We considered the dense regime in which the number of edges per vertex is proportional to the number of vertices. We showed that, as soon as the constraints are \emph{frustrated}, i.e., do not lie on the Erdős-Rényi line, there is breaking of ensemble equivale… ▽ More In [18] we analysed a simple undirected random graph subject to constraints on the total number of edges and the total number of triangles. We considered the dense regime in which the number of edges per vertex is proportional to the number of vertices. We showed that, as soon as the constraints are \emph{frustrated}, i.e., do not lie on the Erdős-Rényi line, there is breaking of ensemble equivalence, in the sense that the specific relative entropy per edge of the \emph{microcanonical ensemble} with respect to the \emph{canonical ensemble} is strictly positive in the limit as the number of vertices tends to infinity. In the present paper we analyse what happens near the Erdős-Rényi line. It turns out that the way in which the specific relative entropy tends to zero depends on whether the total number of triangles is slightly larger or slightly smaller than typical. We investigate what the constrained random graph looks like asymptotically in the microcanonical ensemble. △ Less

Submitted 12 June, 2021; v1 submitted 20 July, 2018; originally announced July 2018.

Comments: arXiv admin note: text overlap with arXiv:1703.08058

Transition time asymptotics of queue-based activation protocols in random-access networks

Authors: Sem Borst, Frank den Hollander, Francesca R. Nardi, Matteo Sfragara

Abstract: We consider networks where each node represents a server with a queue. An active node deactivates at unit rate. An inactive node activates at a rate that depends on its queue length, provided none of its neighbors is active. For complete bipartite networks, in the limit as the queues become large, we compute the average transition time between the two states where one half of the network is active… ▽ More We consider networks where each node represents a server with a queue. An active node deactivates at unit rate. An inactive node activates at a rate that depends on its queue length, provided none of its neighbors is active. For complete bipartite networks, in the limit as the queues become large, we compute the average transition time between the two states where one half of the network is active and the other half is inactive. We show that the law of the transition time divided by its mean exhibits a trichotomy, depending on the activation rate functions. △ Less

Submitted 9 February, 2021; v1 submitted 16 July, 2018; originally announced July 2018.

Comments: 32 pages

MSC Class: 60K25; 60K30; 90B15; 90B18

Journal ref: Stochastic Processes and their Applications, Volume 130, Issue 12, December 2020, Pages 7483-7517

Random walks on dynamic configuration models: a trichotomy

Authors: Luca Avena, Hakan Guldas, Remco van der Hofstad, Frank den Hollander

Abstract: We consider a dynamic random graph on $n$ vertices that is obtained by starting from a random graph generated according to the configuration model with a prescribed degree sequence and at each unit of time randomly rewiring a fraction $α_n$ of the edges. We are interested in the mixing time of a random walk without backtracking on this dynamic random graph in the limit as $n\to\infty$, when $α_n$… ▽ More We consider a dynamic random graph on $n$ vertices that is obtained by starting from a random graph generated according to the configuration model with a prescribed degree sequence and at each unit of time randomly rewiring a fraction $α_n$ of the edges. We are interested in the mixing time of a random walk without backtracking on this dynamic random graph in the limit as $n\to\infty$, when $α_n$ is chosen such that $\lim_{n\to\infty} α_n (\log n)^2 = β\in [0,\infty]$. In [1] we found that, under mild regularity conditions on the degree sequence, the mixing time is of order $1/\sqrt{α_n}$ when $β=\infty$. In the present paper we investigate what happens when $β\in [0,\infty)$. It turns out that the mixing time is of order $\log n$, with the scaled mixing time exhibiting a one-sided cutoff when $β\in (0,\infty)$ and a two-sided cutoff when $β=0$. The occurrence of a one-sided cutoff is a rare phenomenon. In our setting it comes from a competition between the time scales of mixing on the static graph, as identified by Ben-Hamou and Salez [4], and the regeneration time of first stepping across a rewired edge. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Comments: 14 pages, 5 figures

MSC Class: 60K37; 82C27

Random walk in cooling random environment: ergodic limits and concentration inequalities

Authors: L. Avena, Y. Chino, C. da Costa, F. den Hollander

Abstract: In previous work by Avena and den Hollander, a model of a one-dimensional random walk in a dynamic random environment was proposed where the random environment is resampled from a given law along a growing sequence of deterministic times. In the regime where the increments of the resampling times diverge, which is referred to as the cooling regime, a weak law of large numbers and certain fluctuati… ▽ More In previous work by Avena and den Hollander, a model of a one-dimensional random walk in a dynamic random environment was proposed where the random environment is resampled from a given law along a growing sequence of deterministic times. In the regime where the increments of the resampling times diverge, which is referred to as the cooling regime, a weak law of large numbers and certain fluctuation properties were derived under the annealed measure. In the present paper we show that a strong law of large numbers and a quenched large deviation principle hold as well. In the cooling regime, the random walk can be represented as a sum of independent variables, distributed as the increments of a random walk in a static random environment over increasing periods of time. Our proofs require suitable multi-layer decompositions of sums of random variables controlled by moments bounds and concentration estimates. Along the way we derive two results of independent interest, namely, a concentration inequality for the cumulants of the displacement in the static random environment and an ergodic theorem that deals with limits of sums of triangular arrays representing the structure of the cooling regime. We close by discussing our present understanding of homogenisation effects as a function of the speed of divergence of the increments of the resampling times. △ Less

Submitted 8 March, 2018; originally announced March 2018.

Comments: 33 pages, 8 figures

Properties of additive functionals of Brownian motion with resetting

Authors: Frank den Hollander, Satya N. Majumdar, Janusz M. Meylahn, Hugo Touchette

Abstract: We study the distribution of additive functionals of reset Brownian motion, a variation of normal Brownian motion in which the path is interrupted at a given rate and placed back to a given reset position. Our goal is two-fold: (1) For general functionals, we derive a large deviation principle in the presence of resetting and identify the large deviation rate function in terms of a variational for… ▽ More We study the distribution of additive functionals of reset Brownian motion, a variation of normal Brownian motion in which the path is interrupted at a given rate and placed back to a given reset position. Our goal is two-fold: (1) For general functionals, we derive a large deviation principle in the presence of resetting and identify the large deviation rate function in terms of a variational formula involving large deviation rate functions without resetting. (2) For three examples of functionals (positive occupation time, area and absolute area), we investigate the effect of resetting by computing distributions and moments, using a formula that links the generating function with resetting to the generating function without resetting. △ Less

Submitted 25 March, 2019; v1 submitted 30 January, 2018; originally announced January 2018.

Comments: 23 pages, 4 figures

Journal ref: J. Phys. A: Math. Theor. 52, 175001, 2019

Torsional rigidity for cylinders with a Brownian fracture

Authors: M. van den Berg, F. den Hollander

Abstract: We obtain bounds for the expected loss of torsional rigidity of a cylinder $Ω_L=(-L/2,L/2) \times Ω\subset \R^3$ of length $L$ due to a Brownian fracture that starts at a random point in $Ω_L,$ and runs until the first time it exits $Ω_L$. These bounds are expressed in terms of the geometry of the cross-section $Ω\subset \R^2$. It is shown that if $Ω$ is a disc with radius $R$, then in the limit a… ▽ More We obtain bounds for the expected loss of torsional rigidity of a cylinder $Ω_L=(-L/2,L/2) \times Ω\subset \R^3$ of length $L$ due to a Brownian fracture that starts at a random point in $Ω_L,$ and runs until the first time it exits $Ω_L$. These bounds are expressed in terms of the geometry of the cross-section $Ω\subset \R^2$. It is shown that if $Ω$ is a disc with radius $R$, then in the limit as $L \rightarrow \infty$ the expected loss of torsional rigidity equals $cR^5$ for some $c\in (0,\infty)$. We derive bounds for $c$ in terms of the expected Newtonian capacity of the trace of a Brownian path that starts at the centre of a ball in $\R^3$ with radius $1,$ and runs until the first time it exits this ball. △ Less

Submitted 27 November, 2017; originally announced November 2017.

Comments: 18 pages

MSC Class: 35J20; 60G50

Covariance structure behind breaking of ensemble equivalence in random graphs

Authors: Diego Garlaschelli, Frank den Hollander, Andrea Roccaverde

Abstract: For a random graph subject to a topological constraint, the microcanonical ensemble requires the constraint to be met by every realisation of the graph (`hard constraint'), while the canonical ensemble requires the constraint to be met only on average (`soft constraint'). It is known that breaking of ensemble equivalence may occur when the size of the graph tends to infinity, signalled by a non-ze… ▽ More For a random graph subject to a topological constraint, the microcanonical ensemble requires the constraint to be met by every realisation of the graph (`hard constraint'), while the canonical ensemble requires the constraint to be met only on average (`soft constraint'). It is known that breaking of ensemble equivalence may occur when the size of the graph tends to infinity, signalled by a non-zero specific relative entropy of the two ensembles. In this paper we analyse a formula for the relative entropy of generic discrete random structures recently conjectured by Squartini and Garlaschelli. We consider the case of a random graph with a given degree sequence (configuration model), and show that in the dense regime this formula correctly predicts that the specific relative entropy is determined by the scaling of the determinant of the matrix of canonical covariances of the constraints. The formula also correctly predicts that an extra correction term is required in the sparse regime and in the ultra-dense regime. We further show that the different expressions correspond to the degrees in the canonical ensemble being asymptotically Gaussian in the dense regime and asymptotically Poisson in the sparse regime (the latter confirms what we found in earlier work), and the dual degrees in the canonical ensemble being asymptotically Poisson in the ultra-dense regime. In general, we show that the degrees follow a multivariate version of the Poisson-Binomial distribution in the canonical ensemble. △ Less

Submitted 13 August, 2018; v1 submitted 12 November, 2017; originally announced November 2017.

Comments: Journal of Statistical Physics, (), 1-19

Journal ref: Journal of Statistical Physics (2018) 173:644-662

Metastability of hard-core dynamics on bipartite graphs

Authors: Frank den Hollander, Francesca R. Nardi, Siamak Taati

Abstract: We study the metastable behaviour of a stochastic system of particles with hard-core interactions in a high-density regime. Particles sit on the vertices of a bipartite graph. New particles appear subject to a neighbourhood exclusion constraint, while existing particles disappear, all according to independent Poisson clocks. We consider the regime in which the appearance rates are much larger than… ▽ More We study the metastable behaviour of a stochastic system of particles with hard-core interactions in a high-density regime. Particles sit on the vertices of a bipartite graph. New particles appear subject to a neighbourhood exclusion constraint, while existing particles disappear, all according to independent Poisson clocks. We consider the regime in which the appearance rates are much larger than the disappearance rates, and there is a slight imbalance between the appearance rates on the two parts of the graph. Starting from the configuration in which the weak part is covered with particles, the system takes a long time before it reaches the configuration in which the strong part is covered with particles. We obtain a sharp asymptotic estimate for the expected transition time, show that the transition time is asymptotically exponentially distributed, and identify the size and shape of the critical droplet representing the bottleneck for the crossover. For various types of bipartite graphs the computations are made explicit. Proofs rely on potential theory for reversible Markov chains, and on isoperimetric results. In a follow-up paper we will use our results to study the performance of random-access wireless networks. △ Less

Submitted 12 January, 2018; v1 submitted 27 October, 2017; originally announced October 2017.

Comments: 56 pages, 17 figures; small corrections (hypothesis H0, example 4.3, typos)

MSC Class: 60C05; 60K35; 60K37; 82C27

Journal ref: Electronic Journal of Probability 2018, Vol. 23, paper no. 97, 1-65

Annealed scaling for a charged polymer in dimensions two and higher

Authors: Quentin Berger, Frank den Hollander, Julien Poisat

Abstract: This paper considers an undirected polymer chain on $\mathbb{Z}^d$, $d \geq 2$, with i.i.d.\ random charges attached to its constituent monomers. Each self-intersection of the polymer chain contributes an energy to the interaction Hamiltonian that is equal to the product of the charges of the two monomers that meet. The joint probability distribution for the polymer chain and the charges is given… ▽ More This paper considers an undirected polymer chain on $\mathbb{Z}^d$, $d \geq 2$, with i.i.d.\ random charges attached to its constituent monomers. Each self-intersection of the polymer chain contributes an energy to the interaction Hamiltonian that is equal to the product of the charges of the two monomers that meet. The joint probability distribution for the polymer chain and the charges is given by the Gibbs distribution associated with the interaction Hamiltonian. The object of interest is the \emph{annealed free energy} per monomer in the limit as the length $n$ of the polymer chain tends to infinity. We show that there is a critical curve in the parameter plane spanned by the charge bias and the inverse temperature separating an \emph{extended phase} from a \emph{collapsed phase}. We derive the scaling of the critical curve for small and for large charge bias and the scaling of the annealed free energy for small inverse temperature, which are both anomalous. We show that in a subset of the collapsed phase the polymer chain is \emph{subdiffusive}, namely, on scale $(n/\log n)^{1/(d+2)}$ it moves like a Brownian motion conditioned to stay inside a ball with a deterministic radius and a randomly shifted center. We expect this scaling to hold throughout the collapsed phase. We further expect that in the extended phase the polymer chain scales like a weakly self-avoiding walk. Proofs are based on a detailed analysis for simple random walk of the downward large deviations of the self-intersection local time and the upward large deviations of the range. Part of our scaling results are rough, and we formulate conjectures under which they can be sharpened. The existence of the free energy remains an open problem, which we are able to settle in a subset of the collapsed phase for a subclass of charge distributions. △ Less

Submitted 22 August, 2017; originally announced August 2017.

Comments: 35 pages, 6 figures

MSC Class: 82D60; 60K35; 82B27