BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Frank den Hollander

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

ResearcherUnclaimed source record
math.PRmath-phmath.MPcond-mat.dis-nncond-mat.stat-mechmath.COmath.GMnlin.CG

Catalog footprint

What is connected

41works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

41 published item(s)

preprint2026arXiv

The multi-level friendship paradox for sparse random graphs

In Hazra, den Hollander and Parvaneh (2025) we analysed the friendship paradox for sparse random graphs. For four classes of random graphs we characterised the empirical distribution of the friendship biases between vertices and their neighbours at distance $1$, proving convergence as $n\to\infty$ to a limiting distribution, with $n$ the number of vertices, and identifying moments and tail exponents of the limiting distribution. In the present paper we look at the multi-level friendship bias between vertices and their neighbours at distance $k \in \mathbb{N}$ obtained via a $k$-step exploration according to a backtracking or a non-backtracking random walk. We identify the limit of empirical distribution of the multi-level friendship biases as $n\to\infty$ and/or $k\to\infty$. We show that for non-backtracking exploration the two limits commute for a large class of sparse random graphs, including those that locally converge to a rooted Galton-Watson tree. In particular, we show that the same limit arises when $k$ depends on $n$, i.e., $k=k_n$, provided $\lim_{n\to\infty} k_n = \infty$ under some mild conditions. We exhibit cases where the two limits do not commute and show the relevance of the mixing time of the exploration.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Metastability for Glauber dynamics on the complete graph with coupling disorder

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Phase transitions for spatially extended pinning

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)φ(S_i)$ to the interaction Hamiltonian, where $S_i \in \mathbb{Z}$ is the height of monomer $i$ with respect to the interface, $φ: \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 $φ$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Parabolic Anderson Model on a Galton-Watson tree revisited

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The parabolic Anderson model on a Galton-Watson tree

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.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Expansion of percolation critical points for Hamming graphs

The Hamming graph $H(d,n)$ is the Cartesian product of $d$ complete graphs on $n$ vertices. Let $m=d(n-1)$ be the degree and $V = n^d$ be the number of vertices of $H(d,n)$. Let $p_c^{(d)}$ be the critical point for bond percolation on $H(d,n)$. We show that, for $d \in \mathbb N$ fixed and $n \to \infty$, \begin{equation*} p_c^{(d)}= \dfrac{1}{m} + \dfrac{2d^2-1}{2(d-1)^2}\dfrac{1}{m^2} + O(m^{-3}) + O(m^{-1}V^{-1/3}), \end{equation*} which extends the asymptotics found in \cite{BorChaHofSlaSpe05b} by one order. The term $O(m^{-1}V^{-1/3})$ is the width of the critical window. For $d=4,5,6$ we have $m^{-3} = O(m^{-1}V^{-1/3})$, and so the above formula represents the full asymptotic expansion of $p_c^{(d)}$. In \cite{FedHofHolHul16a} \st{we show that} this formula is a crucial ingredient in the study of critical bond percolation on $H(d,n)$ for $d=2,3,4$. The proof uses a lace expansion for the upper bound and a novel comparison with a branching random walk for the lower bound. The proof of the lower bound also yields a refined asymptotics for the susceptibility of a subcritical Erdős-Rényi random graph.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Metastability for Glauber dynamics on random graphs

In this paper we study metastable behaviour at low temperature of Glauber spin-flip dynamics on random graphs. We fix a large number of vertices and randomly allocate edges according to the Configuration Model with a prescribed degree distribution. Each vertex carries a spin that can point either up or down. Each spin interacts with a positive magnetic field, while spins at vertices that are connected by edges also interact with each other via a ferromagnetic pair potential. We start from the configuration where all spins point down, and allow spins to flip up or down according to a Metropolis dynamics at positive temperature. We are interested in the time it takes the system to reach the configuration where all spins point up. In order to achieve this transition, the system needs to create a sufficiently large droplet of up-spins, called critical droplet, which triggers the crossover. In the limit as the temperature tends to zero, and subject to a certain \emph{key hypothesis} implying metastable behaviour, the average crossover time follows the classical \emph{Arrhenius law}, with an exponent and a prefactor that are controlled by the \emph{energy} and the \emph{entropy} of the critical droplet. The crossover time divided by its average is exponentially distributed. We study the scaling behaviour of the exponent as the number of vertices tends to infinity, deriving upper and lower bounds. We also identify a regime for the magnetic field and the pair potential in which the key hypothesis is satisfied. The critical droplets, representing the saddle points for the crossover, have a size that is of the order of the number of vertices. This is because the random graphs generated by the Configuration Model are expander graphs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multi-colony Wright-Fisher with seed-bank

We consider a multi-colony version of the Wright-Fisher model with seed-bank that was recently introduced by Blath et al. Individuals live in colonies and change type via resampling and mutation. Each colony contains a seed-bank that acts as a genetic reservoir. Individuals can enter the seed-bank and become dormant or can exit the seed-bank and become active. In each colony at each generation a fixed fraction of individuals swap state, drawn randomly from the active and the dormant population. While dormant, individuals suspend their resampling. While active, individuals resample from their own colony, but also from different colonies according to a random walk transition kernel representingmigration. Both active and dormant individuals mutate. We derive a formula for the probability that two individuals drawn randomly from two given colonies are identical by descent, i.e., share a common ancestor. This formula, which is formulated in Fourier language, is valid when the colonies form a discrete torus. We consider the special case of a symmetric slow seed-bank, for which in each colony half of the individuals are in the seed-bank and at each generation the fraction of individuals that swap state is small. This leads to a simpler formula, from which we are able to deduce how the probability to be identical by descent depends on the distance between the two colonies and various relevant parameters. Through an analysis of random walk Green functions, we are able to derive explicit scaling expressions when mutation is slower than migration. We also compute the spatial second moment of the probability to be identical by descent for all parameters when the torus becomes large. For the special case of a symmetric slow seed-bank, we again obtain explicit scaling expressions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phase diagram for a copolymer in a micro-emulsion

In this paper we study a model describing a copolymer in a micro-emulsion. The copolymer consists of a random concatenation of hydrophobic and hydrophilic monomers, the micro-emulsion consists of large blocks of oil and water arranged in a percolation-type fashion. The interaction Hamiltonian assigns energy $-α$ to hydrophobic monomers in oil and energy $-β$ to hydrophilic monomers in water, where $α,β$ are parameters that without loss of generality are taken to lie in the cone $\{(α,β) \in\mathbb{R}^2\colon\,α\geq |β|\}$. Depending on the values of these parameters, the copolymer either stays close to the oil-water interface (localization) or wanders off into the oil and/or the water (delocalization). Based on an assumption about the strict concavity of the free energy of a copolymer near a linear interface, we derive a variational formula for the quenched free energy per monomer that is column-based, i.e., captures what the copolymer does in columns of different type. We subsequently transform this into a variational formula that is slope-based, i.e., captures what the polymer does as it travels at different slopes, and we use the latter to identify the phase diagram in the $(α,β)$-cone. There are two regimes: supercritical (the oil blocks percolate) and subcritical (the oil blocks do not percolate). The supercritical and the subcritical phase diagram each have two localized phases and two delocalized phases, separated by four critical curves meeting at a quadruple critical point. The different phases correspond to the different ways in which the copolymer can move through the micro-emulsion. The analysis of the phase diagram is based on three hypotheses of percolation-type on the blocks. We show that these three hypotheses are plausible, but do not provide a proof.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Breaking of ensemble equivalence in networks

It is generally believed that, in the thermodynamic limit, the microcanonical description as a function of energy coincides with the canonical description as a function of temperature. However, various examples of systems for which the microcanonical and canonical ensembles are not equivalent have been identified. A complete theory of this intriguing phenomenon is still missing. Here we show that ensemble nonequivalence can manifest itself also in random graphs with topological constraints. We find that, while graphs with a given number of links are ensemble-equivalent, graphs with a given degree sequence are not. This result holds irrespective of whether the energy is nonadditive (as in unipartite graphs) or additive (as in bipartite graphs). In contrast with previous expectations, our results show that: (1) physically, nonequivalence can be induced by an extensive number of local constraints, and not necessarily by long-range interactions or nonadditivity; (2) mathematically, nonquivalence is determined by a different large-deviation behaviour of microcanonical and canonical probabilities for a single microstate, and not necessarily for almost all microstates. The latter criterion, which is entirely local, is not restricted to networks and holds in general.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Parabolic Anderson model in a dynamic random environment: random conductances

The parabolic Anderson model is defined as the partial differential equation \partial u(x,t)/\partial t = κΔu(x,t) + ξ(x,t)u(x,t), x\in\Z^d, t\geq 0, where κ\in [0,\infty) is the diffusion constant, Δis the discrete Laplacian, and ξis a dynamic random environment that drives the equation. The initial condition u(x,0)=u_0(x), x\in\Z^d, is taken to be non-negative and bounded. The solution of the parabolic Anderson equation describes the evolution of a field of particles performing independent simple random walks with binary branching: particles jump at rate 2dκ, split into two at rate ξ\vee 0, and die at rate (-ξ) \vee 0. Our focus is on the Lyapunov exponents λ_p(κ) = \lim_{t\to\infty} \frac{1}{t} \log \E([u(0,t)]^p)^{1/p}, p \in \N, and λ_0(κ) = \lim_{t\to\infty} \frac{1}{t}\log u(0,t). We investigate what happens when κΔis replaced by Δ^\cK, where \cK = \{\mathcal{K}(x,y)\colon\,x,y\in\Z^d,\,x \sim y\} is a collection of random conductances between neighbouring sites replacing the constant conductances κin the homogeneous model. We show that the associated annealed Lyapunov exponents are given by the formula λ_p(\cK) = \sup\{λ_p(κ) \colon\,κ\in\Supp(\cK)\}, where \Supp(\cK) is the set of values taken by the \cK-field. We also show that for the associated quenched Lyapunov exponent this formula only provides a lower bound. Our proof is valid for three classes of reversible ξ, and for all \cK satisfying a certain clustering property, namely, there are arbitrarily large balls where \cK is almost constant and close to any value in \Supp(\cK). What our result says is that the Lyapunov exponents are controlled by those pockets of \cK where the conductances are close to the value that maximises the growth in the homogeneous setting.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Random Walk on Random Walks

In this paper we study a random walk in a one-dimensional dynamic random environment consisting of a collection of independent particles performing simple symmetric random walks in a Poisson equilibrium with density $ρ\in (0,\infty)$. At each step the random walk performs a nearest-neighbour jump, moving to the right with probability $p_{\circ}$ when it is on a vacant site and probability $p_{\bullet}$ when it is on an occupied site. Assuming that $p_\circ \in (0,1)$ and $p_\bullet \neq \tfrac12$, we show that the position of the random walk satisfies a strong law of large numbers, a functional central limit theorem and a large deviation bound, provided $ρ$ is large enough. The proof is based on the construction of a renewal structure together with a multiscale renormalisation argument.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gibbs-non-Gibbs dynamical transitions for mean-field interacting Brownian motions

We consider a system of real-valued spins interacting with each other through a mean-field Hamiltonian that depends on the empirical magnetization of the spins via a general potential. The system is subjected to a stochastic dynamics where the spins perform independent Brownian motions. As in \cite{FedHoMa13}, which considers the Curie-Weiss model with Ising spins interacting via a quadratic potential and subjected to independent spins flips, we follow the program outlined in \cite{vEFedHoRe10}. We show that in the thermodynamic limit the system is non-Gibbs at time $t \in (0,\infty)$ if and only if there exists an $α\in \mathbb{R}$ such that the large deviation rate function for the trajectory of the magnetization conditional on hitting the value $α$ at time $t$ has multiple global minimizers. We further show that different minimizing trajectories are different at time $t=0$. We give conditions on the potential under which the system is Gibbs at time $t=0$, classify the possible scenarios of being Gibbs at time $t \in (0,\infty)$ in terms of the second difference quotient of the potential, and show that the system cannot become Gibbs once it has become non-Gibbs, i.e., there is a unique and explicitly computable crossover time $t_c \in [0,\infty]$ from Gibbs to non-Gibbs. We give examples of immediate loss of Gibbsianness ($t_c=0$), short-time conservation of Gibbsianness, large-time loss of Gibbsianness ($t_c\in (0,\infty)$), and preservation of Gibbsianness ($t_c=\infty$). Depending on the potential, the system can be Gibbs or non-Gibbs at the cross-over time time $t=t_c$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Renormalisation of hierarchically interacting Cannings processes

The present paper brings a new class of interacting jump processes into focus. We start from a single-colony $C^Λ$-process, which arises as the continuum-mass limit of a $Λ$-Cannings individual-based population model, where $Λ$ is a finite non-negative measure that describes the offspring mechanism. After that we introduce a system of hierarchically interacting $C^Λ$-processes, where the interaction comes from migration and reshuffling-resampling based on measures $(Λ_k)_{k}$ both acting in $k$-blocks of the hierarchical group. We refer to this system as the $C_N^{c,Λ}$-process. The dual process of the $C_N^{c,Λ}$-process is a spatial coalescent with multi-level block coalescence. For the above system we carry out a full renormalisation analysis in the hierarchical mean-field limit $N\to\infty$. Our main result is that, in the limit as $N\to\infty$, on each scale $k\in\mathbb{N}_0$ the $k$-block averages of the $C_N^{c,Λ}$-process converge to a random process that is a superposition of a $C^{Λ_k}$-process and a Fleming-Viot process, the latter with a volatility $d_k$ and with a drift of strength $c_k$ towards the limiting $(k+1)$-block average. It turns out that $d_k$ is a function of $c_l$ and $Λ_l$ for all $0\leq l<k$. Thus, it is through the volatility that the renormalisation manifests itself. We discuss the implications of the scaling of $d_k$ for the behaviour on large space-time scales of the $C_N^{c,Λ}$-process. We compare the outcome with what is known from the renormalisation analysis of hierarchically interacting Fleming-Viot diffusions, pointing out several new features. We obtain a new classification for when the process exhibits clustering, respectively, exhibits local coexistence. Finally, we show that for finite $N$ the same dichotomy between clustering and local coexistence holds as for $N\to\infty$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A general smoothing inequality for disordered polymers

This note sharpens the smoothing inequality of Giacomin and Toninelli for disordered polymers. This inequality is shown to be valid for any disorder distribution with locally finite exponential moments, and to provide an asymptotically sharp constant for weak disorder. A key tool in the proof is an estimate that compares the effect on the free energy of tilting, respectively, shifting the disorder distribution. This estimate holds in large generality (way beyond disordered polymers) and is of independent interest.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extremal geometry of a Brownian porous medium

The path W[0,t] of a Brownian motion on a d-dimensional torus T^d run for time t is a random compact subset of T^d. We study the geometric properties of the complement T^d \ W[0,t] for t large and d >= 3. In particular, we show that the largest regions in this complement have a linear scale phi = [(d log t)/(d-2)kt]^{1/(d-2)}, where k is the capacity of the unit ball. More specifically, we identify the sets E for which T^d \ W[0,t] contains a translate of phi E, and we count the number of disjoint such translates. Furthermore, we derive large deviation principles for the largest inradius of T^d \ W[0,t] for t large and the epsilon-cover time of T^d for epsilon small. Our results, which generalise laws of large numbers proved by Dembo, Peres and Rosen, are based on a large deviation principle for the shape of the component with largest capacity in T^d \ W_rho[0,t], where W_rho[0,t] is the Wiener sausage of radius rho = rho(t), with rho(t) chosen much smaller than phi but not too small. The idea behind this choice is that T^d \ W[0,t] consists of "lakes", whose linear size is of order phi, connected by narrow "channels". We also derive large deviation principles for the principal Dirichlet eigenvalue and for the maximal volume of the components of T^d \ W_rho[0,t] for t large. Our results give a complete picture of the extremal geometry of T^d \ W[0,t] and of the optimal strategy for W[0,t] to realise the extremes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large deviation principles for words drawn from correlated letter sequences

When an i.i.d.\ sequence of letters is cut into words according to i.i.d.\ renewal times, an i.i.d.\ sequence of words is obtained. In the \emph{annealed} LDP (large deviation principle) for the empirical process of words, the rate function is the specific relative entropy of the observed law of words w.r.t.\ the reference law of words. In Birkner, Greven and den Hollander \cite{BGdH10} the \emph{quenched} LDP (= conditional on a typical letter sequence) was derived for the case where the renewal times have an \emph{algebraic} tail. The rate function turned out to be a sum of two terms, one being the annealed rate function, the other being proportional to the specific relative entropy of the observed law of letters w.r.t.\ the reference law of letters, obtained by concatenating the words and randomising the location of the origin. The proportionality constant equals the tail exponent of the renewal process. The purpose of the present paper is to extend both LDP's to letter sequences that are not i.i.d. It is shown that both LDP's carry over when the letter sequence satisfies a mixing condition called \emph{summable variation}. The rate functions are again given by specific relative entropies w.r.t.\ the reference law of words, respectively, letters. But since neither of these reference laws is i.i.d., several approximation arguments are needed to obtain the extension.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Random walk in a high density dynamic random environment

The goal of this note is to prove a law of large numbers for the empirical speed of a green particle that performs a random walk on top of a field of red particles which themselves perform independent simple random walks on $\Z^d$, $d \geq 1$. The red particles jump at rate 1 and are in a Poisson equilibrium with density $μ$. The green particle also jumps at rate 1, but uses different transition kernels $p'$ and $p''$ depending on whether it sees a red particle or not. It is shown that, in the limit as $μ\to\infty$, the speed of the green particle tends to the average jump under $p'$. This result is far from surprising, but it is non-trivial to prove. The proof that is given in this note is based on techniques that were developed in \cite{KeSi} to deal with spread-of-infection models. The main difficulty is that, due to particle conservation, space-time correlations in the field of red particles decay slowly. This places the problem in a class of random walks in dynamic random environments for which scaling laws are hard to obtain.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scaling of a random walk on a supercritical contact process

A proof is provided of a strong law of large numbers for a one-dimensional random walk in a dynamic random environment given by a supercritical contact process in equilibrium. The proof is based on a coupling argument that traces the space-time cones containing the infection clusters generated by single infections and uses that the random walk eventually gets trapped inside the union of these cones. For the case where the local drifts of the random walk are smaller than the speed at which infection clusters grow, the random walk eventually gets trapped inside a single cone. This in turn leads to the existence of regeneration times at which the random walk forgets its past. The latter are used to prove a functional central limit theorem and a large deviation principle. The qualitative dependence of the speed, the volatility and the rate function on the infection parameter is investigated, and some open problems are mentioned.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The parabolic Anderson model in a dynamic random environment: basic properties of the quenched Lyapunov exponent

In this paper we study the parabolic Anderson equation \partial u(x,t)/\partial t=κΔu(x,t)+ξ(x,t)u(x,t), x\in\Z^d, t\geq 0, where the u-field and the ξ-field are \R-valued, κ\in [0,\infty) is the diffusion constant, and $Δ$ is the discrete Laplacian. The initial condition u(x,0)=u_0(x), x\in\Z^d, is taken to be non-negative and bounded. The solution of the parabolic Anderson equation describes the evolution of a field of particles performing independent simple random walks with binary branching: particles jump at rate 2dκ, split into two at rate ξ\vee 0, and die at rate (-ξ)\vee 0. Our goal is to prove a number of basic properties of the solution u under assumptions on $ξ$ that are as weak as possible. Throughout the paper we assume that $ξ$ is stationary and ergodic under translations in space and time, is not constant and satisfies \E(|ξ(0,0)|)<\infty, where \E denotes expectation w.r.t. ξ. Under a mild assumption on the tails of the distribution of ξ, we show that the solution to the parabolic Anderson equation exists and is unique for all κ\in [0,\infty). Our main object of interest is the quenched Lyapunov exponent λ_0(κ)=\lim_{t\to\infty}\frac{1}{t}\log u(0,t). Under certain weak space-time mixing conditions on ξ, we show the following properties: (1)λ_0(κ) does not depend on the initial condition u_0; (2)λ_0(κ)<\infty for all κ\in [0,\infty); (3)κ\mapsto λ_0(κ) is continuous on [0,\infty) but not Lipschitz at 0. We further conjecture: (4)\lim_{κ\to\infty}[λ_p(κ)-λ_0(κ)]=0 for all p\in\N, where λ_p (κ)=\lim_{t\to\infty}\frac{1}{pt}\log\E([u(0,t)]^p) is the p-th annealed Lyapunov exponent. Finally, we prove that our weak space-time mixing conditions on ξare satisfied for several classes of interacting particle systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The parabolic Anderson model in a dynamic random environment: space-time ergodicity for the quenched Lyapunov exponent

We continue our study of the parabolic Anderson equation $\partial u(x,t)/\partial t = κΔu(x,t) + ξ(x,t)u(x,t)$, $x\in\Z^d$, $t\geq 0$, where $κ\in [0,\infty)$ is the diffusion constant, $Δ$ is the discrete Laplacian, and $ξ$ plays the role of a \emph{dynamic random environment} that drives the equation. The initial condition $u(x,0)=u_0(x)$, $x\in\Z^d$, is taken to be non-negative and bounded. The solution of the parabolic Anderson equation describes the evolution of a field of particles performing independent simple random walks with binary branching: particles jump at rate $2dκ$, split into two at rate $ξ\vee 0$, and die at rate $(-ξ) \vee 0$. We assume that $ξ$ is stationary and ergodic under translations in space and time, is not constant and satisfies $\E(|ξ(0,0)|)<\infty$, where $\E$ denotes expectation w.r.t.\ $ξ$. Our main object of interest is the quenched Lyapunov exponent $λ_0 (κ) = \lim_{t\to\infty} \frac{1}{t}\log u(0,t)$. In earlier work we showed that under certain mild space-time mixing assumptions the limit exists $ξ$-a.s., is finite and continuous on $[0,\infty)$, is globally Lipschitz on $(0,\infty)$, is not Lipschitz at 0, and satisfies $λ_0(0) = \E(ξ(0,0))$ and $λ_0(κ) > \E(ξ(0,0))$ for $κ\in (0,\infty)$.In the present paper we show that $\lim_{κ\to\infty} λ_0(κ) =\E(ξ(0,0))$ under an additional space-time mixing condition on $ξ$. This result shows that the parabolic Anderson model exhibits space-time ergodicity in the limit of large diffusivity. This fact is interesting because there are choices of $ξ$ that fulfill our assumption for which the annealed Lyapunov exponent $λ_1(κ) = \lim_{t\to\infty} \frac{1}{t}\log \E(u(0,t))$ is infinite on $[0,\infty)$, a situation that is referred to as strongly catalytic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variational characterization of the critical curve for pinning of random polymers

In this paper we look at the pinning of a directed polymer by a one-dimensional linear interface carrying random charges. There are two phases, localized and delocalized, depending on the inverse temperature and on the disorder bias. Using quenched and annealed large deviation principles for the empirical process of words drawn from a random letter sequence according to a random renewal process [Birkner, Greven and den Hollander, Probab. Theory Related Fields 148 (2010) 403-456], we derive variational formulas for the quenched, respectively, annealed critical curve separating the two phases. These variational formulas are used to obtain a necessary and sufficient criterion, stated in terms of relative entropies, for the two critical curves to be different at a given inverse temperature, a property referred to as relevance of the disorder. This criterion in turn is used to show that the regimes of relevant and irrelevant disorder are separated by a unique inverse critical temperature. Subsequently, upper and lower bounds are derived for the inverse critical temperature, from which sufficient conditions under which it is strictly positive, respectively, finite are obtained. The former condition is believed to be necessary as well, a problem that we will address in a forthcoming paper. Random pinning has been studied extensively in the literature. The present paper opens up a window with a variational view. Our variational formulas for the quenched and the annealed critical curve are new and provide valuable insight into the nature of the phase transition. Our results on the inverse critical temperature drawn from these variational formulas are not new, but they offer an alternative approach, that is, flexible enough to be extended to other models of random polymers with disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variational description of Gibbs-non-Gibbs dynamical transitions for spin-flip systems with a Kac-type interaction

We continue our study of Gibbs-non-Gibbs dynamical transitions. In the present paper we consider a system of Ising spins on a large discrete torus with a Kac-type interaction subject to an independent spin-flip dynamics (infinite-temperature Glauber dynamics). We show that, in accordance with the program outlined in \cite{vEFedHoRe10}, in the thermodynamic limit Gibbs-non-Gibbs dynamical transitions are \emph{equivalent} to bifurcations in the set of global minima of the large-deviation rate function for the trajectories of the empirical density \emph{conditional} on their endpoint. More precisely, the time-evolved measure is non-Gibbs if and only if this set is not a singleton for \emph{some} value of the endpoint. A partial description of the possible scenarios of bifurcation is given, leading to a characterization of passages from Gibbs to non-Gibbs and vice versa, with sharp transition times. Our analysis provides a conceptual step-up from our earlier work on Gibbs-non-Gibbs dynamical transitions for the Curie-Weiss model, where the mean-field interaction allowed us to focus on trajectories of the empirical magnetization rather than the empirical density.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Copolymer with pinning: variational characterization of the phase diagram

This paper studies a polymer chain in the vicinity of a linear interface separating two immiscible solvents. The polymer consists of random monomer types, while the interface carries random charges. Both the monomer types and the charges are given by i.i.d. sequences of random variables. The configurations of the polymer are directed paths that can make i.i.d. excursions of finite length above and below the interface. The Hamiltonian has two parts: a monomer-solvent interaction ("copolymer") and a monomer-interface interaction ("pinning"). The quenched and the annealed version of the model each undergo a transition from a localized phase (where the polymer stays close to the interface) to a delocalized phase (where the polymer wanders away from the interface). We exploit the approach developed in [5] and [3] to derive variational formulas for the quenched and the annealed free energy per monomer. These variational formulas are analyzed to obtain detailed information on the critical curves separating the two phases and on the typical behavior of the polymer in each of the two phases. Our main results settle a number of open questions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Free energy of a copolymer in a micro-emulsion

In this paper we consider a two-dimensional model of a copolymer consisting of a random concatenation of hydrophilic and hydrophobic monomers, immersed in a micro-emulsion of random droplets of oil and water. The copolymer interacts with the micro-emulsion through an interaction Hamiltonian that favors matches and disfavors mismatches between the monomers and the solvents, in such a way that the interaction with the oil is stronger than with the water. The configurations of the copolymers are directed self-avoiding paths in which only steps up, down and right are allowed. The configurations of the micro-emulsion are square blocks with oil and water arranged in percolation-type fashion. The only restriction imposed on the path is that in every column of blocks its vertical displacement on the block scale is bounded. The way in which the copolymer enters and exits successive columns of blocks is a directed self-avoiding path as well, but on the block scale. We refer to this path as the coarse-grained self-avoiding path. We are interested in the limit as the copolymer and the blocks become large, in such a way that the copolymer spends a long time in each block yet visits many blocks. This is a coarse-graining limit in which the space-time scales of the copolymer and of the micro-emulsion become separated. We derive a variational formula for the quenched free energy per monomer, where quenched means that the disorder in the copolymer and the disorder in the micro-emulsion are both frozen. In a sequel paper we will analyze this variational formula and identify the phase diagram. It turns out that there are two regimes, supercritical and subcritical, depending on whether the oil blocks percolate or not along the coarse-grained self-avoiding path. The phase diagrams in the two regimes turn out to be completely different.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kawasaki dynamics with two types of particles: critical droplets

This is the third in a series of three papers in which we study a two-dimensional lattice gas consisting of two types of particles subject to Kawasaki dynamics at low temperature in a large finite box with an open boundary. Each pair of particles occupying neighboring sites has a negative binding energy provided their types are different, while each particle has a positive activation energy that depends on its type. There is no binding energy between particles of the same type. At the boundary of the box particles are created and annihilated in a way that represents the presence of an infinite gas reservoir. We start the dynamics from the empty box and are interested in the transition time to the full box. This transition is triggered by a critical droplet appearing somewhere in the box. In the first paper we identified the parameter range for which the system is metastable, showed that the first entrance distribution on the set of critical droplets is uniform, computed the expected transition time up to and including a multiplicative factor of order one, and proved that the nucleation time divided by its expectation is exponentially distributed, all in the limit of low temperature. These results were proved under three hypotheses, and involved three model-dependent quantities: the energy, the shape and the number of critical droplets. Here prove the third hypothesis and identify the shape and the number of critical droplets. The geometric properties of subcritical, critical and supercritical droplets, which are crucial in determining the metastable behavior of the system are identified. The geometry turns out to be considerably more complex than for Kawasaki dynamics with one type of particle, for which an extensive literature exists. The main motivation behind our work is to understand metastability of multi-type particle systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Law of large numbers for non-elliptic random walks in dynamic random environments

We prove a law of large numbers for a class of $\Z^d$-valued random walks in dynamic random environments, including non-elliptic examples. We assume for the random environment a mixing property called \emph{conditional cone-mixing} and that the random walk tends to stay inside wide enough space-time cones. The proof is based on a generalization of a regeneration scheme developed by Comets and Zeitouni for static random environments and adapted by Avena, den Hollander and Redig to dynamic random environments. A number of one-dimensional examples are given. In some cases, the sign of the speed can be determined.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lectures on random polymers

These lecture notes are a guided tour through the fascinating world of polymer chains interacting with themselves and/or with their environment. The focus is on the mathematical description of a number of physical and chemical phenomena, with particular emphasis on phase transitions and space-time scaling. The topics covered, though only a selection, are typical for the area. Sections 1-3 describe models of polymers without disorder, Sections 4-6 models of polymers with disorder. Appendices A-E contain tutorials in which a number of key techniques are explained in more detail.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pinning of a random walk by a random walk: proof of a conjecture

In [3] the radius of convergence of the generating function of the collision local time of two independent copies of an irreducible, symmetric and transient random walk on Zd, d \geq 1, was studied. Two versions were considered: z1, the radius of convergence when one walk is averaged over and the other is kept fixed; z2, the radius of convergence when both walks are averaged over. While z2 can be easily computed, no explicit expression is available for z1. In [3] it was conjectured that z1 > z2 under a weak regularity assumption on the random walk. However, this gap was only proved for strongly transient random walk. Subsequently, in [1], [4] and [5] the gap was proved for a subclass of random walks that are transient but not strongly transient. In the present note we settle the full conjecture. The proof is based on a comparison of variational formulas for z1 and z2 derived in [3]. Under additional weak regularity assumptions on the random walk, we derive an upper bound for z1. We further show that the gap persists when the random walk has exponential waiting times. The gaps imply the presence of an intermediate phase for a model of two random walks with a pinning interaction. As explained in [3], they also imply the presence of intermediate phases in three classes of interacting stochastic systems: coupled branching processes, interacting diffusions, and directed polymers with bulk disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Variational description of Gibbs-non-Gibbs dynamical transitions for the Curie-Weiss model

We perform a detailed study of Gibbs-non-Gibbs transitions for the Curie-Weiss model subject to independent spin-flip dynamics ("infinite-temperature" dynamics). We show that, in this setup, the program outlined in van Enter, Fernández, den Hollander and Redig can be fully completed, namely that Gibbs-non-Gibbs transitions are equivalent to bifurcations in the set of global minima of the large-deviation rate function for the trajectories of the magnetization conditioned on their endpoint. As a consequence, we show that the time-evolved model is non-Gibbs if and only if this set is not a singleton for some value of the final magnetization. A detailed description of the possible scenarios of bifurcation is given, leading to a full characterization of passages from Gibbs to non-Gibbs -and vice versa- with sharp transition times (under the dynamics Gibbsianness can be lost and can be recovered). Our analysis expands the work of Ermolaev and Kulske who considered zero magnetic field and finite-temperature spin-flip dynamics. We consider both zero and non-zero magnetic field but restricted to infinite-temperature spin-flip dynamics. Our results reveal an interesting dependence on the interaction parameters, including the presence of forbidden regions for the optimal trajectories and the possible occurrence of overshoots and undershoots in the optimal trajectories. The numerical plots provided are obtained with the help of MATHEMATICA.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A crossover for the bad configurations of random walk in random scenery

In this paper, we consider a random walk and a random color scenery on Z. The increments of the walk and the colors of the scenery are assumed to be i.i.d. and to be independent of each other. We are interested in the random process of colors seen by the walk in the course of time. Bad configurations for this random process are the discontinuity points of the conditional probability distribution for the color seen at time zero given the colors seen at all later times. We focus on the case where the random walk has increments 0, +1 or -1 with probability epsilon, (1-epsilon)p and (1-epsilon)(1-p), respectively, with p in [1/2,1] and epsilon in [0,1), and where the scenery assigns the color black or white to the sites of Z with probability 1/2 each. We show that, remarkably, the set of bad configurations exhibits a crossover: for epsilon=0 and p in (1/2,4/5) all configurations are bad, while for (p,epsilon) in an open neighborhood of (1,0) all configurations are good. In addition, we show that for epsilon=0 and p=1/2 both bad and good configurations exist. We conjecture that for all epsilon in [0,1) the crossover value is unique and equals 4/5. Finally, we suggest an approach to handle the seemingly more difficult case where epsilon>0 and p in [1/2,4/5), which will be pursued in future work.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kawasaki dynamics with two types of particles: stable/metastable configurations and communication heights

This is the second in a series of three papers in which we study a two-dimensional lattice gas consisting of two types of particles subject to Kawasaki dynamics at low temperature in a large finite box with an open boundary. Each pair of particles occupying neighboring sites has a negative binding energy provided their types are different, while each particle has a positive activation energy that depends on its type. There is no binding energy between particles of the same type. At the boundary of the box particles are created and annihilated in a way that represents the presence of an infinite gas reservoir. We start the dynamics from the empty box and are interested in the transition time to the full box. This transition is triggered by a critical droplet appearing somewhere in the box. In the first paper we identified the parameter range for which the system is metastable, showed that the first entrance distribution on the set of critical droplets is uniform, computed the expected transition time up to and including a multiplicative factor of order one, and proved that the nucleation time divided by its expectation is exponentially distributed, all in the limit of low temperature. These results were proved under three hypotheses, and involve three model-dependent quantities: the energy, the shape and the number of critical droplets. In this second paper we prove the first and the second hypothesis and identify the energy of critical droplets. The paper deals with understanding the geometric properties of subcritical, critical and supercritical droplets, which are crucial in determining the metastable behavior of the system. The geometry turns out to be considerably more complex than for Kawasaki dynamics with one type of particle, for which an extensive literature exists. The main motivation behind our work is to understand metastability of multi- type particle systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quenched Lyapunov exponent for the parabolic Anderson model in a dynamic random environment

We continue our study of the parabolic Anderson equation $\partial u/\partial t = κΔu + γξu$ for the space-time field $u\colon\,\Z^d\times [0,\infty)\to\R$, where $κ\in [0,\infty)$ is the diffusion constant, $Δ$ is the discrete Laplacian, $γ\in (0,\infty)$ is the coupling constant, and $ξ\colon\,\Z^d\times [0,\infty)\to\R$ is a space-time random environment that drives the equation. The solution of this equation describes the evolution of a "reactant" $u$ under the influence of a "catalyst" $ξ$, both living on $\Z^d$. In earlier work we considered three choices for $ξ$: independent simple random walks, the symmetric exclusion process, and the symmetric voter model, all in equilibrium at a given density. We analyzed the \emph{annealed} Lyapunov exponents, i.e., the exponential growth rates of the successive moments of $u$ w.r.t.\ $ξ$, and showed that these exponents display an interesting dependence on the diffusion constant $κ$, with qualitatively different behavior in different dimensions $d$. In the present paper we focus on the \emph{quenched} Lyapunov exponent, i.e., the exponential growth rate of $u$ conditional on $ξ$. We first prove existence and derive some qualitative properties of the quenched Lyapunov exponent for a general $ξ$ that is stationary and ergodic w.r.t.\ translations in $\Z^d$ and satisfies certain noisiness conditions. After that we focus on the three particular choices for $ξ$ mentioned above and derive some more detailed properties. We close by formulating a number of open problems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A large-deviation view on dynamical Gibbs-non-Gibbs transitions

We develop a space-time large-deviation point of view on Gibbs-non-Gibbs transitions in spin systems subject to a stochastic spin-flip dynamics. Using the general theory for large deviations of functionals of Markov processes outlined in Feng and Kurtz [11], we show that the trajectory under the spin-flip dynamics of the empirical measure of the spins in a large block in Z^d satisfies a large deviation principle in the limit as the block size tends to infinity. The associated rate function can be computed as the action functional of a Lagrangian that is the Legendre transform of a certain non-linear generator, playing a role analogous to the moment-generating function in the Gartner-Ellis theorem of large deviation theory when this is applied to finite-dimensional Markov processes. This rate function is used to define the notion of "bad empirical measures", which are the discontinuity points of the optimal trajectories (i.e., the trajectories minimizing the rate function) given the empirical measure at the end of the trajectory. The dynamical Gibbs-non-Gibbs transitions are linked to the occurrence of bad empirical measures: for short times no bad empirical measures occur, while for intermediate and large times bad empirical mea- sures are possible. A future research program is proposed to classify the various possible scenarios behind this crossover, which we refer to as a "nature-versus-nurture" transition.

0 citations0 reviewsNo code linkedOpen access