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

Eyal Lubetzky

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

ResearcherUnclaimed source record
math.PRmath.COmath-phmath.MPInformation Theorymath.ITData Structures and AlgorithmsComputer Science and Game TheoryDiscrete Mathematics

Catalog footprint

What is connected

47works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

47 published item(s)

preprint2022arXiv

On the limiting law of line ensembles of Brownian polymers with geometric area tilts

We study the line ensembles of non-crossing Brownian bridges above a hard wall, each tilted by the area of the region below it with geometrically growing pre-factors. This model, which mimics the level lines of the $(2+1)$D SOS model above a hard wall, was studied in two works from 2019 by Caputo, Ioffe and Wachtel. In those works, the tightness of the law of the top $k$ paths, for any fixed $k$, was established under either zero or free boundary conditions, which in the former setting implied the existence of a limit via a monotonicity argument. Here we address the open problem of a limit under free boundary conditions: we prove that as the interval length, followed by the number of paths, go to $\infty$, the top $k$ paths converge to the same limit as in the free boundary case, as conjectured by Caputo, Ioffe and Wachtel.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The threshold for stacked triangulations

A \emph{stacked triangulation} of a $d$-simplex $\mathbf{o}=\{1,\ldots,d+1\}$ ($d\geq 2$) is a triangulation obtained by repeatedly subdividing a $d$-simplex into $d+1$ new ones via a new vertex (the case $d=2$ is known as an Appolonian network). We study the occurrence of such a triangulation in the Linial--Meshulam model, i.e., for which $p$ does the random simplicial complex $Y\sim \mathcal{Y}_d(n,p)$ contain the faces of a stacked triangulation of the $d$-simplex $\mathbf{o}$, with its internal vertices labeled in $[n]$. In the language of bootstrap percolation in hypergraphs, it pertains to the threshold for $K_{d+2}^{d+1}$, the $(d+1)$-uniform clique on $d+2$ vertices. Our main result identifies this threshold for every $d\geq 2$, showing it is asymptotically $(α_d n)^{-1/d}$, where $α_d$ is the growth rate of the Fuss--Catalan numbers of order $d$. The proof hinges on a second moment argument in the supercritical regime, and on Kalai's algebraic shifting in the subcritical regime.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotics in percolation on high-girth expanders

We consider supercritical bond percolation on a family of high-girth $d$-regular expanders. Alon, Benjamini and Stacey (2004) established that its critical probability for the appearance of a linear-sized ("giant'') component is $p_c=1/(d-1)$. Our main result recovers the sharp asymptotics of the size and degree distribution of the vertices in the giant and its 2-core at any $p>p_c$. It was further shown in [ABS04] that the second largest component, at any $0<p<1$, has size at most $n^ω$ for some $ω<1$. We show that, unlike the situation in the classical Erdős-Rényi random graph, the second largest component in bond percolation on a regular expander, even with an arbitrarily large girth, can have size $n^{ω'}$ for $ω'$ arbitrarily close to $1$. Moreover, as a by-product of that construction, we answer negatively a question of Benjamini (2013) on the relation between the diameter of a component in percolation on expanders and the existence of a giant component. Finally, we establish other typical features of the giant component, e.g., the existence of a linear path.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cycle lengths in sparse random graphs

We study the set ${\cal L}(G)$ of lengths of all cycles that appear in a random $d$-regular $G$ on $n$ vertices for a fixed $d\geq 3$, as well as in Erdős--Rényi random graphs on $n$ vertices with a fixed average degree $c>1$. Fundamental results on the distribution of cycle counts in these models were established in the 1980's and early 1990's, with a focus on the extreme lengths: cycles of fixed length, and cycles of length linear in $n$. Here we derive, for a random $d$-regular graph, the limiting probability that ${\cal L}(G)$ simultaneously contains the entire range $\{\ell,\ldots,n\}$ for $\ell\geq 3$, as an explicit expression $θ_\ell=θ_\ell(d)\in(0,1)$ which goes to $1$ as $\ell\to\infty$. For the random graph ${\cal G}(n,p)$ with $p=c/n$, where $c\geq C_0$ for some absolute constant $C_0$, we show the analogous result for the range $\{\ell,\ldots,(1-o(1))L_{\max}(G)\}$, where $L_{\max}$ is the length of a longest cycle in $G$. The limiting probability for ${\cal G}(n,p)$ coincides with $θ_\ell$ from the $d$-regular case when $c$ is the integer $d-1$. In addition, for the directed random graph ${\cal D}(n,p)$ we show results analogous to those on ${\cal G}(n,p)$, and for both models we find an interval of $c ε^2 n$ consecutive cycle lengths in the slightly supercritical regime $p=\frac{1+ε}n$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Empirical spectral distributions of sparse random graphs

We study the spectrum of a random multigraph with a degree sequence ${\bf D}_n=(D_i)_{i=1}^n$ and average degree $1 \ll ω_n \ll n$, generated by the configuration model, and also the spectrum of the analogous random simple graph. We show that, when the empirical spectral distribution (ESD) of $ω_n^{-1} {\bf D}_n $ converges weakly to a limit $ν$, under mild moment assumptions (e.g., $D_i/ω_n$ are i.i.d. with a finite second moment), the ESD of the normalized adjacency matrix converges in probability to $ν\boxtimes σ_{\rm sc}$, the free multiplicative convolution of $ν$ with the semicircle law. Relating this limit with a variant of the Marchenko--Pastur law yields the continuity of its density (away from zero), and an effective procedure for determining its support. Our proof of convergence is based on a coupling between the random simple graph and multigraph with the same degrees, which might be of independent interest. We further construct and rely on a coupling of the multigraph to an inhomogeneous Erdős-Rényi graph with the target ESD, using three intermediate random graphs, with a negligible fraction of edges modified in each step.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximum and shape of interfaces in 3D Ising crystals

Dobrushin (1972) showed that the interface of a 3D Ising model with minus boundary conditions above the $xy$-plane and plus below is rigid (has $O(1)$-fluctuations) at every sufficiently low temperature. Since then, basic features of this interface -- such as the asymptotics of its maximum -- were only identified in more tractable random surface models that approximate the Ising interface at low temperatures, e.g., for the (2+1)D Solid-On-Solid model. Here we study the large deviations of the interface of the 3D Ising model in a cube of side-length $n$ with Dobrushin's boundary conditions, and in particular obtain a law of large numbers for $M_n$, its maximum: if the inverse-temperature $β$ is large enough, then $M_n / \log n \to 2/α_β$ as $n\to\infty$, in probability, where $α_β$ is given by a large deviation rate in infinite volume. We further show that, on the large deviation event that the interface connects the origin to height $h$, it consists of a 1D spine that behaves like a random walk, in that it decomposes into a linear (in $h$) number of asymptotically-stationary weakly-dependent increments that have exponential tails. As the number $T$ of increments diverges, properties of the interface such as its surface area, volume, and the location of its tip, all obey CLTs with variances linear in $T$. These results generalize to every dimension $d\geq 3$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximum of Branching Brownian motion in a periodic environment

We study the maximum of Branching Brownian motion (BBM) with branching rates that vary in space, via a periodic function of a particle's location. This corresponds to a variant of the F-KPP equation in a periodic medium, extensively studied in the last 15 years, admitting pulsating fronts as solutions. Recent progress on this PDE due to Hamel, Nolen, Roquejoffre and Ryzhik ('16) implies tightness for the centered maximum of BBM in a periodic environment. Here we establish the convergence in distribution of specific subsequences of this centered maximum, and identify the limiting distribution. Consequently, we find the asymptotic shift between the solution to the corresponding F-KPP equation with Heavyside initial data and the pulsating wave, thereby answering a question of Hamel et al. Analogous results are given for the cases where the Brownian motion is replaced by an Ito diffusion with periodic coefficients, as well as for nearest-neighbor branching random walks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Noise sensitivity of critical random graphs

We study noise sensitivity of properties of the largest components $({\cal C}_j)_{j\geq 1}$ of the random graph ${\cal G}(n,p)$ in its critical window $p=(1+λn^{-1/3})/n$. For instance, is the property "$|{\cal C}_1|$ exceeds its median size" noise sensitive? Roberts and Şengül (2018) proved that the answer to this is yes if the noise $ε$ is such that $ε\gg n^{-1/6}$, and conjectured the correct threshold is $ε\gg n^{-1/3}$. That is, the threshold for sensitivity should coincide with the critical window---as shown for the existence of long cycles by the first author and Steif (2015). We prove that for $ε\gg n^{-1/3}$ the pair of vectors $ n^{-2/3}(|{\cal C}_j|)_{j\geq 1}$ before and after the noise converges in distribution to a pair of i.i.d. random variables, whereas for $ε\ll n^{-1/3}$ the $\ell^2$-distance between the two goes to 0 in probability. This confirms the above conjecture: any Boolean function of the vector of rescaled component sizes is sensitive in the former case and stable in the latter. We also look at the effect of the noise on the metric space $n^{-1/3}({\cal C}_j)_{j\geq 1}$. E.g., for $ε\geq n^{-1/3+o(1)}$, we show that the joint law of the spaces before and after the noise converges to a product measure, implying noise sensitivity of any property seen in the limit, e.g., "the diameter of ${\cal C}_1$ exceeds its median."

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tightness and tails of the maximum in 3D Ising interfaces

Consider the 3D Ising model on a box of side length $n$ with minus boundary conditions above the $xy$-plane and plus boundary conditions below it. At low temperatures, Dobrushin (1972) showed that the interface separating the predominantly plus and predominantly minus regions is localized: its height above a fixed point has exponential tails. Recently, the authors proved a law of large numbers for the maximum height $M_n$ of this interface: for every $β$ large, $M_n/ \log n\to c_β$ in probability as $n\to\infty$. Here we show that the laws of the centered maxima $(M_n - \mathbb{E}[M_n])_{n\geq 1}$ are uniformly tight. Moreover, even though this sequence does not converge, we prove that it has uniform upper and lower Gumbel tails (exponential right tails and doubly exponential left tails). Key to the proof is a sharp (up to $O(1)$ precision) understanding of the surface large deviations. This includes, in particular, the shape of a pillar that reaches near-maximum height, even at its base, where the interactions with neighboring pillars are dominant.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universality for Langevin-like spin glass dynamics

We study dynamics for asymmetric spin glass models, proposed by Hertz et al. and Sompolinsky et al. in the 1980's in the context of neural networks: particles evolve via a modified Langevin dynamics for the Sherrington--Kirkpatrick model with soft spins, whereby the disorder is i.i.d. standard Gaussian rather than symmetric. Ben Arous and Guionnet (1995), followed by Guionnet (1997), proved for Gaussian interactions that as the number of particles grows, the short-term empirical law of this dynamics converges a.s. to a non-random law $μ_\star$ of a ``self-consistent single spin dynamics,'' as predicted by physicists. Here we obtain universality of this fact: For asymmetric disorder given by i.i.d. variables of zero mean, unit variance and exponential or better tail decay, at every temperature, the empirical law of sample paths of the Langevin-like dynamics in a fixed time interval has the same a.s. limit $μ_\star$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cutoff on all Ramanujan graphs

We show that on every Ramanujan graph $G$, the simple random walk exhibits cutoff: when $G$ has $n$ vertices and degree $d$, the total-variation distance of the walk from the uniform distribution at time $t=\frac{d}{d-2}\log_{d-1} n + s\sqrt{\log n}$ is asymptotically $\mathbb{P}(Z > c\, s)$ where $Z$ is a standard normal variable and $c=c(d)$ is an explicit constant. Furthermore, for all $1 \leq p \leq \infty$, $d$-regular Ramanujan graphs minimize the asymptotic $L^p$-mixing time for SRW among all $d$-regular graphs. Our proof also shows that, for every vertex $x$ in $G$ as above, its distance from $n-o(n)$ of the vertices is asymptotically $\log_{d-1} n$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On replica symmetry of large deviations in random graphs

The following question is due to Chatterjee and Varadhan (2011). Fix $0<p<r<1$ and take $G\sim G(n,p)$, the Erdős-Rényi random graph with edge density $p$, conditioned to have at least as many triangles as the typical $G(n,r)$. Is $G$ close in cut-distance to a typical $G(n,r)$? Via a beautiful new framework for large deviation principles in $G(n,p)$, Chatterjee and Varadhan gave bounds on the replica symmetric phase, the region of $(p,r)$ where the answer is positive. They further showed that for any small enough $p$ there are at least two phase transitions as $r$ varies. We settle this question by identifying the replica symmetric phase for triangles and more generally for any fixed $d$-regular graph. By analyzing the variational problem arising from the framework of Chatterjee and Varadhan we show that the replica symmetry phase consists of all $(p,r)$ such that $(r^d,h_p(r))$ lies on the convex minorant of $x\mapsto h_p(x^{1/d})$ where $h_p$ is the rate function of a binomial with parameter $p$. In particular, the answer for triangles involves $h_p(\sqrt{x})$ rather than the natural guess of $h_p(x^{1/3})$ where symmetry was previously known. Analogous results are obtained for linear hypergraphs as well as the setting where the largest eigenvalue of $G\sim G(n,p)$ is conditioned to exceed the typical value of the largest eigenvalue of $G(n,r)$. Building on the work of Chatterjee and Diaconis (2012) we obtain additional results on a class of exponential random graphs including a new range of parameters where symmetry breaking occurs. En route we give a short alternative proof of a graph homomorphism inequality due to Kahn (2001) and Galvin and Tetali (2004).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Random walks on the random graph

We study random walks on the giant component of the Erdős-Rényi random graph ${\cal G}(n,p)$ where $p=λ/n$ for $λ>1$ fixed. The mixing time from a worst starting point was shown by Fountoulakis and Reed, and independently by Benjamini, Kozma and Wormald, to have order $\log^2 n$. We prove that starting from a uniform vertex (equivalently, from a fixed vertex conditioned to belong to the giant) both accelerates mixing to $O(\log n)$ and concentrates it (the cutoff phenomenon occurs): the typical mixing is at $(ν{\bf d})^{-1}\log n \pm (\log n)^{1/2+o(1)}$, where $ν$ and ${\bf d}$ are the speed of random walk and dimension of harmonic measure on a ${\rm Poisson}(λ)$-Galton-Watson tree. Analogous results are given for graphs with prescribed degree sequences, where cutoff is shown both for the simple and for the non-backtracking random walk.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Information percolation and cutoff for the stochastic Ising model

We introduce a new framework for analyzing Glauber dynamics for the Ising model. The traditional approach for obtaining sharp mixing results has been to appeal to estimates on spatial properties of the stationary measure from within a multi-scale analysis of the dynamics. Here we propose to study these simultaneously by examining "information percolation" clusters in the space-time slab. Using this framework, we obtain new results for the Ising model on $(\mathbb{Z}/n\mathbb{Z})^d$ throughout the high temperature regime: total-variation mixing exhibits cutoff with an $O(1)$-window around the time at which the magnetization is the square-root of the volume. (Previously, cutoff in the full high temperature regime was only known for $d\leq 2$, and only with an $O(\log\log n)$-window.) Furthermore, the new framework opens the door to understanding the effect of the initial state on the mixing time. We demonstrate this on the 1D Ising model, showing that starting from the uniform ("disordered") initial distribution asymptotically halves the mixing time, whereas almost every deterministic starting state is asymptotically as bad as starting from the ("ordered") all-plus state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strong noise sensitivity and random graphs

The noise sensitivity of a Boolean function describes its likelihood to flip under small perturbations of its input. Introduced in the seminal work of Benjamini, Kalai and Schramm [Inst. Hautes Études Sci. Publ. Math. 90 (1999) 5-43], it was there shown to be governed by the first level of Fourier coefficients in the central case of monotone functions at a constant critical probability $p_c$. Here we study noise sensitivity and a natural stronger version of it, addressing the effect of noise given a specific witness in the original input. Our main context is the Erdős-Rényi random graph, where already the property of containing a given graph is sufficiently rich to separate these notions. In particular, our analysis implies (strong) noise sensitivity in settings where the BKS criterion involving the first Fourier level does not apply, for example, when $p_c\to0$ polynomially fast in the number of variables.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An exposition to information percolation for the Ising model

Information percolation is a new method for analyzing stochastic spin systems through classifying and controlling the clusters of information-flow in the space-time slab. It yielded sharp mixing estimates (cutoff with an $O(1)$-window) for the Ising model on $Z^d$ up to the critical temperature, as well as results on the effect of initial conditions on mixing. In this expository note we demonstrate the method on lattices (more generally, on any locally-finite transitive graph) at very high temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cover time of a random graph with a degree sequence II: Allowing vertices of degree two

We study the cover time of a random graph chosen uniformly at random from the set of graphs with vertex set $[n]$ and degree sequence $\mathbf{d}=(d_i)_{i=1}^n$. In a previous work, the asymptotic cover time was obtained under a number of assumptions on $\mathbf{d}$, the most significant being that $d_i\geq 3$ for all $i$. Here we replace this assumption by $d_i\geq 2$. As a corollary, we establish the asymptotic cover time for the 2-core of the emerging giant component of $\mathcal{G}(n,p)$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cutoff for the East process

The East process is a 1D kinetically constrained interacting particle system, introduced in the physics literature in the early 90's to model liquid-glass transitions. Spectral gap estimates of Aldous and Diaconis in 2002 imply that its mixing time on $L$ sites has order $L$. We complement that result and show cutoff with an $O(\sqrt{L})$-window. The main ingredient is an analysis of the front of the process (its rightmost zero in the setup where zeros facilitate updates to their right). One expects the front to advance as a biased random walk, whose normal fluctuations would imply cutoff with an $O(\sqrt{L})$-window. The law of the process behind the front plays a crucial role: Blondel showed that it converges to an invariant measure $ν$, on which very little is known. Here we obtain quantitative bounds on the speed of convergence to $ν$, finding that it is exponentially fast. We then derive that the increments of the front behave as a stationary mixing sequence of random variables, and a Stein-method based argument of Bolthausen ('82) implies a CLT for the location of the front, yielding the cutoff result. Finally, we supplement these results by a study of analogous kinetically constrained models on trees, again establishing cutoff, yet this time with an $O(1)$-window.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cycle factors and renewal theory

For which values of $k$ does a uniformly chosen $3$-regular graph $G$ on $n$ vertices typically contain $ n/k$ vertex-disjoint $k$-cycles (a $k$-cycle factor)? To date, this has been answered for $k=n$ and for $k \ll \log n$; the former, the Hamiltonicity problem, was finally answered in the affirmative by Robinson and Wormald in 1992, while the answer in the latter case is negative since with high probability most vertices do not lie on $k$-cycles. Here we settle the problem completely: the threshold for a $k$-cycle factor in $G$ as above is $κ_0 \log_2 n$ with $κ_0=[1-\frac12\log_2 3]^{-1}\approx 4.82$. Precisely, we prove a 2-point concentration result: if $k \geq κ_0 \log_2(2n/e)$ divides $n$ then $G$ contains a $k$-cycle factor w.h.p., whereas if $k<κ_0\log_2(2n/e)-\frac{\log^2 n}n$ then w.h.p. it does not. As a byproduct, we confirm the "Comb Conjecture," an old problem concerning the embedding of certain spanning trees in the random graph $G(n,p)$. The proof follows the small subgraph conditioning framework, but the associated second moment analysis here is far more delicate than in any earlier use of this method and involves several novel features, among them a sharp estimate for tail probabilities in renewal processes without replacement which may be of independent interest.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of $(2+1)$-dimensional SOS surfaces above a wall: Slow mixing induced by entropic repulsion

We study the Glauber dynamics for the $(2+1)\mathrm{D}$ Solid-On-Solid model above a hard wall and below a far away ceiling, on an $L\times L$ box of $\mathbb{Z}^2$ with zero boundary conditions, at large inverse-temperature $β$. It was shown by Bricmont, El Mellouki and Fröhlich [J. Stat. Phys. 42 (1986) 743-798] that the floor constraint induces an entropic repulsion effect which lifts the surface to an average height $H\asymp(1/β)\log L$. As an essential step in understanding the effect of entropic repulsion on the Glauber dynamics we determine the equilibrium height $H$ to within an additive constant: $H=(1/4β)\log L+O(1)$. We then show that starting from zero initial conditions the surface rises to its final height $H$ through a sequence of metastable transitions between consecutive levels. The time for a transition from height $h=aH$, $a\in(0,1)$, to height $h+1$ is roughly $\exp(cL^a)$ for some constant $c>0$. In particular, the mixing time of the dynamics is exponentially large in $L$, that is, $T_{\mathrm{MIX}}\geq e^{cL}$. We also provide the matching upper bound $T_{\mathrm{MIX}}\leq e^{c'L}$, requiring a challenging analysis of the statistics of height contours at low temperature and new coupling ideas and techniques. Finally, to emphasize the role of entropic repulsion we show that without a floor constraint at height zero the mixing time is no longer exponentially large in $L$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Harmonic pinnacles in the Discrete Gaussian model

The 2D Discrete Gaussian model gives each height function $η: \mathbb{Z}^2\to\mathbb{Z}$ a probability proportional to $\exp(-β\mathcal{H}(η))$, where $β$ is the inverse-temperature and $\mathcal{H}(η) = \sum_{x\sim y}(η_x-η_y)^2$ sums over nearest-neighbor bonds. We consider the model at large fixed $β$, where it is flat unlike its continuous analog (the Gaussian Free Field). We first establish that the maximum height in an $L\times L$ box with 0 boundary conditions concentrates on two integers $M,M+1$ with $M\sim \sqrt{(1/2πβ)\log L\log\log L}$. The key is a large deviation estimate for the height at the origin in $\mathbb{Z}^2$, dominated by "harmonic pinnacles", integer approximations of a harmonic variational problem. Second, in this model conditioned on $η\geq 0$ (a floor), the average height rises, and in fact the height of almost all sites concentrates on levels $H,H+1$ where $H\sim M/\sqrt{2}$. This in particular pins down the asymptotics, and corrects the order, in results of Bricmont, El-Mellouki and Fröhlich (1986), where it was argued that the maximum and the height of the surface above a floor are both of order $\sqrt{\log L}$. Finally, our methods extend to other classical surface models (e.g., restricted SOS), featuring connections to $p$-harmonic analysis and alternating sign matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The threshold for combs in random graphs

For $k\mid n$ let $Comb_{n,k}$ denote the tree consisting of an $(n/k)$-vertex path with disjoint $k$-vertex paths beginning at each of its vertices. An old conjecture says that for any $k=k(n)$ the threshold for the random graph $G(n,p)$ to contain $Comb_{n,k}$ is at $p\asymp \frac{\log n}n$. Here we verify this for $k \leq C\log n$ with any fixed $C>0$. In a companion paper, using very different methods, we treat the complementary range, proving the conjecture for $k\geq κ_0 \log n$ (with $κ_0\approx 4.82$).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universality of cutoff for the Ising model

On any locally-finite geometry, the stochastic Ising model is known to be contractive when the inverse-temperature $β$ is small enough, via classical results of Dobrushin and of Holley in the 1970's. By a general principle proposed by Peres, the dynamics is then expected to exhibit cutoff. However, so far cutoff for the Ising model has been confirmed mainly for lattices, heavily relying on amenability and log Sobolev inequalities. Without these, cutoff was unknown at any fixed $β>0$, no matter how small, even in basic examples such as the Ising model on a binary tree or a random regular graph. We use the new framework of information percolation to show that, in any geometry, there is cutoff for the Ising model at high enough temperatures. Precisely, on any sequence of graphs with maximum degree $d$, the Ising model has cutoff provided that $β<κ/d$ for some absolute constant $κ$ (a result which, up to the value of $κ$, is best possible). Moreover, the cutoff location is established as the time at which the sum of squared magnetizations drops to 1, and the cutoff window is $O(1)$, just as when $β=0$. Finally, the mixing time from almost every initial state is not more than a factor of $1+ε_β$ faster then the worst one (with $ε_β\to0$ as $β\to 0$), whereas the uniform starting state is at least $2-ε_β$ times faster.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scaling limit and cube-root fluctuations in SOS surfaces above a wall

Consider the classical $(2+1)$-dimensional Solid-On-Solid model above a hard wall on an $L\times L$ box of $\bbZ^2$. The model describes a crystal surface by assigning a non-negative integer height $η_x$ to each site $x$ in the box and 0 heights to its boundary. The probability of a surface configuration $η$ is proportional to $\exp(-β\mathcal{H}(η))$, where $β$ is the inverse-temperature and $\mathcal{H}(η)$ sums the absolute values of height differences between neighboring sites. We give a full description of the shape of the SOS surface for low enough temperatures. First we show that with high probability the height of almost all sites is concentrated on two levels, $H(L)=\lfloor (1/4β)\log L\rfloor$ and $H(L)-1$. Moreover, for most values of $L$ the height is concentrated on the single value $H(L)$. Next, we study the ensemble of level lines corresponding to the heights $(H(L),H(L)-1,...)$. We prove that w.h.p. there is a unique macroscopic level line for each height. Furthermore, when taking a diverging sequence of system sizes $L_k$, the rescaled macroscopic level line at height $H(L_k)-n$ has a limiting shape if the fractional parts of $(1/4β)\log L_k$ converge to a noncritical value. The scaling limit is an explicit convex subset of the unit square $Q$ and its boundary has a flat component on the boundary of $Q$. Finally, the highest macroscopic level line has $L_k^{1/3+o(1)}$ fluctuations along the flat part of the boundary of its limiting shape.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stochastic coalescence in logarithmic time

The following distributed coalescence protocol was introduced by Dahlia Malkhi in 2006 motivated by applications in social networking. Initially there are n agents wishing to coalesce into one cluster via a decentralized stochastic process, where each round is as follows: every cluster flips a fair coin to dictate whether it is to issue or accept requests in this round. Issuing a request amounts to contacting a cluster randomly chosen proportionally to its size. A cluster accepting requests is to select an incoming one uniformly (if there are such) and merge with that cluster. Empirical results by Fernandess and Malkhi suggested the protocol concludes in O(log n) rounds with high probability, whereas numerical estimates by Oded Schramm, based on an ingenious analytic approximation, suggested that the coalescence time should be super-logarithmic. Our contribution is a rigorous study of the stochastic coalescence process with two consequences. First, we confirm that the above process indeed requires super-logarithmic time w.h.p., where the inefficient rounds are due to oversized clusters that occasionally develop. Second, we remedy this by showing that a simple modification produces an essentially optimal distributed protocol; if clusters favor their smallest incoming merge request then the process does terminate in O(log n) rounds w.h.p., and simulations show that the new protocol readily outperforms the original one. Our upper bound hinges on a potential function involving the logarithm of the number of clusters and the cluster-susceptibility, carefully chosen to form a supermartingale. The analysis of the lower bound builds upon the novel approach of Schramm which may find additional applications: rather than seeking a single parameter that controls the system behavior, instead one approximates the system by the Laplace transform of the entire cluster-size distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anatomy of the giant component: The strictly supercritical regime

In a recent work of the authors and Kim, we derived a complete description of the largest component of the Erdős-Rényi random graph $G(n,p)$ as it emerges from the critical window, i.e. for $p = (1+ε)/n$ where $ε^3 n \to\infty$ and $ε=o(1)$, in terms of a tractable contiguous model. Here we provide the analogous description for the supercritical giant component, i.e., the largest component of $G(n,p)$ for $p = λ/n$ where $λ>1$ is fixed. The contiguous model is roughly as follows: Take a random degree sequence and sample a random multigraph with these degrees to arrive at the kernel; Replace the edges by paths whose lengths are i.i.d. geometric variables to arrive at the 2-core; Attach i.i.d. Poisson Galton-Watson trees to the vertices for the final giant component. As in the case of the emerging giant, we obtain this result via a sequence of contiguity arguments at the heart of which are Kim's Poisson-cloning method and the Pittel-Wormald local limit theorems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cutoff for general spin systems with arbitrary boundary conditions

The cutoff phenomenon describes a sharp transition in the convergence of a Markov chain to equilibrium. In recent work, the authors established cutoff and its location for the stochastic Ising model on the $d$-dimensional torus $(Z/nZ)^d$ for any $d\geq 1$. The proof used the symmetric structure of the torus and monotonicity in an essential way. Here we enhance the framework and extend it to general geometries, boundary conditions and external fields to derive a cutoff criterion that involves the growth rate of balls and the log-Sobolev constant of the Glauber dynamics. In particular, we show there is cutoff for stochastic Ising on any sequence of bounded-degree graphs with sub-exponential growth under arbitrary external fields provided the inverse log-Sobolev constant is bounded. For lattices with homogenous boundary, such as all-plus, we identify the cutoff location explicitly in terms of spectral gaps of infinite-volume dynamics on half-plane intersections. Analogous results establishing cutoff are obtained for non-monotone spin-systems at high temperatures, including the gas hard-core model, the Potts model, the anti-ferromagnetic Potts model and the coloring model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Glauber Dynamics for the mean-field Potts Model

We study Glauber dynamics for the mean-field (Curie-Weiss) Potts model with $q\geq 3$ states and show that it undergoes a critical slowdown at an inverse-temperature $β_s(q)$ strictly lower than the critical $β_c(q)$ for uniqueness of the thermodynamic limit. The dynamical critical $β_s(q)$ is the spinodal point marking the onset of metastability. We prove that when $β<β_s(q)$ the mixing time is asymptotically $C(β, q) n \log n$ and the dynamics exhibits the cutoff phenomena, a sharp transition in mixing, with a window of order $n$. At $β=β_s(q)$ the dynamics no longer exhibits cutoff and its mixing obeys a power-law of order $n^{4/3}$. For $β>β_s(q)$ the mixing time is exponentially large in $n$. Furthermore, as $β\uparrow β_s$ with $n$, the mixing time interpolates smoothly from subcritical to critical behavior, with the latter reached at a scaling window of $O(n^{-2/3})$ around $β_s$. These results form the first complete analysis of mixing around the critical dynamical temperature --- including the critical power law --- for a model with a first order phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mixing time of near-critical random graphs

Let $\mathcal{C}_1$ be the largest component of the Erdős--Rényi random graph $\mathcal{G}(n,p)$. The mixing time of random walk on $\mathcal {C}_1$ in the strictly supercritical regime, $p=c/n$ with fixed $c>1$, was shown to have order $\log^2n$ by Fountoulakis and Reed, and independently by Benjamini, Kozma and Wormald. In the critical window, $p=(1+\varepsilon)/n$ where $λ=\varepsilon^3n$ is bounded, Nachmias and Peres proved that the mixing time on $\mathcal{C}_1$ is of order $n$. However, it was unclear how to interpolate between these results, and estimate the mixing time as the giant component emerges from the critical window. Indeed, even the asymptotics of the diameter of $\mathcal{C}_1$ in this regime were only recently obtained by Riordan and Wormald, as well as the present authors and Kim. In this paper, we show that for $p=(1+\varepsilon)/n$ with $λ=\varepsilon^3n\to\infty$ and $λ=o(n)$, the mixing time on $\mathcal{C}_1$ is with high probability of order $(n/λ)\log^2λ$. In addition, we show that this is the order of the largest mixing time over all components, both in the slightly supercritical and in the slightly subcritical regime [i.e., $p=(1-\varepsilon)/n$ with $λ$ as above].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Random greedy triangle-packing beyond the 7/4 barrier

The random greedy algorithm for constructing a large partial Steiner-Triple-System is defined as follows. Begin with a complete graph on $n$ vertices and proceed to remove the edges of triangles one at a time, where each triangle removed is chosen uniformly at random out of all remaining triangles. This stochastic process terminates once it arrives at a triangle-free graph, and a longstanding open problem is to estimate the final number of edges, or equivalently the time it takes the process to conclude. The intuition that the edge distribution is roughly uniform at all times led to a folklore conjecture that the final number of edges is $n^{3/2+o(1)}$ with high probability, whereas the best known upper bound is $n^{7/4+o(1)}$. It is no coincidence that various methods break precisely at the exponent 7/4 as it corresponds to the inherent barrier where co-degrees become comparable to the variations in their values that arose earlier in the process. In this work we significantly improve upon the previous bounds by establishing that w.h.p. the number of edges in the final graph is at most $ n^{5/3+o(1)} $. Our approach relies on a system of martingales used to control key graph parameters, where the crucial new idea is to harness the self-correcting nature of the process in order to control these parameters well beyond the point where their early variation matches the order of their expectation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Random triangle removal

Starting from a complete graph on $n$ vertices, repeatedly delete the edges of a uniformly chosen triangle. This stochastic process terminates once it arrives at a triangle-free graph, and the fundamental question is to estimate the final number of edges (equivalently, the time it takes the process to finish, or how many edge-disjoint triangles are packed via the random greedy algorithm). Bollobás and Erdős (1990) conjectured that the expected final number of edges has order $n^{3/2}$, motivated by the study of the Ramsey number $R(3,t)$. An upper bound of $o(n^2)$ was shown by Spencer (1995) and independently by Rödl and Thoma (1996). Several bounds were given for variants and generalizations (e.g., Alon, Kim and Spencer (1997) and Wormald (1999)), while the best known upper bound for the original question of Bollobás and Erdős was $n^{7/4+o(1)}$ due to Grable (1997). No nontrivial lower bound was available. Here we prove that with high probability the final number of edges in random triangle removal is equal to $n^{3/2+o(1)}$, thus confirming the 3/2 exponent conjectured by Bollobás and Erdős and matching the predictions of Spencer et al. For the upper bound, for any fixed $ε>0$ we construct a family of $\exp(O(1/ε))$ graphs by gluing $O(1/ε)$ triangles sequentially in a prescribed manner, and dynamically track all homomorphisms from them, rooted at any two vertices, up to the point where $n^{3/2+ε}$ edges remain. A system of martingales establishes concentration for these random variables around their analogous means in a random graph with corresponding edge density, and a key role is played by the self-correcting nature of the process. The lower bound builds on the estimates at that very point to show that the process will typically terminate with at least $n^{3/2-o(1)}$ edges left.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The shape of the $(2+1)$D SOS surface above a wall

We give a full description for the shape of the classical (2+1)\Dim Solid-On-Solid model above a wall, introduced by Temperley (1952). On an $L\times L$ box at a large inverse-temperature $β$ the height of most sites concentrates on a single level $h = \lfloor (1/4β)\log L\rfloor$ for most values of $L$. For a sequence of diverging boxes the ensemble of level lines of heights $(h,h-1,...)$ has a scaling limit in Hausdorff distance iff the fractional parts of $(1/4β)\log L$ converge to a noncritical value. The scaling limit is explicitly given by nested distinct loops formed via translates of Wulff shapes. Finally, the $h$-level lines feature $L^{1/3+o(1)}$ fluctuations from the side boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Index coding via linear programming

Index Coding has received considerable attention recently motivated in part by real-world applications and in part by its connection to Network Coding. The basic setting of Index Coding encodes the problem input as an undirected graph and the fundamental parameter is the broadcast rate $β$, the average communication cost per bit for sufficiently long messages (i.e. the non-linear vector capacity). Recent nontrivial bounds on $β$ were derived from the study of other Index Coding capacities (e.g. the scalar capacity $β_1$) by Bar-Yossef et al (2006), Lubetzky and Stav (2007) and Alon et al (2008). However, these indirect bounds shed little light on the behavior of $β$: there was no known polynomial-time algorithm for approximating $β$ in a general network to within a nontrivial (i.e. $o(n)$) factor, and the exact value of $β$ remained unknown for any graph where Index Coding is nontrivial. Our main contribution is a direct information-theoretic analysis of the broadcast rate $β$ using linear programs, in contrast to previous approaches that compared $β$ with graph-theoretic parameters. This allows us to resolve the aforementioned two open questions. We provide a polynomial-time algorithm with a nontrivial approximation ratio for computing $β$ in a general network along with a polynomial-time decision procedure for recognizing instances with $β=2$. In addition, we pinpoint $β$ precisely for various classes of graphs (e.g. for various Cayley graphs of cyclic groups) thereby simultaneously improving the previously known upper and lower bounds for these graphs. Via this approach we construct graphs where the difference between $β$ and its trivial lower bound is linear in the number of vertices and ones where $β$ is uniformly bounded while its upper bound derived from the naive encoding scheme is polynomially worse.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lexicographic products and the power of non-linear network coding

We introduce a technique for establishing and amplifying gaps between parameters of network coding and index coding. The technique uses linear programs to establish separations between combinatorial and coding-theoretic parameters and applies hypergraph lexicographic products to amplify these separations. This entails combining the dual solutions of the lexicographic multiplicands and proving that they are a valid dual of the product. Our result is general enough to apply to a large family of linear programs. This blend of linear programs and lexicographic products gives a recipe for constructing hard instances in which the gap between combinatorial or coding-theoretic parameters is polynomially large. We find polynomial gaps in cases in which the largest previously known gaps were only small constant factors or entirely unknown. Most notably, we show a polynomial separation between linear and non-linear network coding rates. This involves exploiting a connection between matroids and index coding to establish a previously unknown separation between linear and non-linear index coding rates. We also construct index coding problems with a polynomial gap between the broadcast rate and the trivial lower bound for which no gap was previously known.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A note on the random greedy triangle-packing algorithm

The random greedy algorithm for constructing a large partial Steiner-Triple-System is defined as follows. We begin with a complete graph on $n$ vertices and proceed to remove the edges of triangles one at a time, where each triangle removed is chosen uniformly at random from the collection of all remaining triangles. This stochastic process terminates once it arrives at a triangle-free graph. In this note we show that with high probability the number of edges in the final graph is at most $ O\big( n^{7/4}\log^{5/4}n \big) $.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Choice-memory tradeoff in allocations

In the classical balls-and-bins paradigm, where $n$ balls are placed independently and uniformly in $n$ bins, typically the number of bins with at least two balls in them is $Θ(n)$ and the maximum number of balls in a bin is $Θ(\frac{\log n}{\log \log n})$. It is well known that when each round offers $k$ independent uniform options for bins, it is possible to typically achieve a constant maximal load if and only if $k=Ω(\log n)$. Moreover, it is possible w.h.p. to avoid any collisions between $n/2$ balls if $k>\log_2n$. In this work, we extend this into the setting where only $m$ bits of memory are available. We establish a tradeoff between the number of choices $k$ and the memory $m$, dictated by the quantity $km/n$. Roughly put, we show that for $km\gg n$ one can achieve a constant maximal load, while for $km\ll n$ no substantial improvement can be gained over the case $k=1$ (i.e., a random allocation). For any $k=Ω(\log n)$ and $m=Ω(\log^2n)$, one can achieve a constant load w.h.p. if $km=Ω(n)$, yet the load is unbounded if $km=o(n)$. Similarly, if $km>Cn$ then $n/2$ balls can be allocated without any collisions w.h.p., whereas for $km<εn$ there are typically $Ω(n)$ collisions. Furthermore, we show that the load is w.h.p. at least $\frac{\log(n/m)}{\log k+\log\log(n/m)}$. In particular, for $k\leq\operatorname {polylog}(n)$, if $m=n^{1-δ}$ the optimal maximal load is $Θ(\frac{\log n}{\log\log n})$ (the same as in the case $k=1$), while $m=2n$ suffices to ensure a constant load. Finally, we analyze nonadaptive allocation algorithms and give tight upper and lower bounds for their performance.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical Ising on the square lattice mixes in polynomial time

The Ising model is widely regarded as the most studied model of spin-systems in statistical physics. The focus of this paper is its dynamic (stochastic) version, the Glauber dynamics, introduced in 1963 and by now the most popular means of sampling the Ising measure. Intensive study throughout the last three decades has yielded a rigorous understanding of the spectral-gap of the dynamics on $\Z^2$ everywhere except at criticality. While the critical behavior of the Ising model has long been the focus for physicists, mathematicians have only recently developed an understanding of its critical geometry with the advent of SLE, CLE and new tools to study conformally invariant systems. A rich interplay exists between the static and dynamic models. At the static phase-transition for Ising, the dynamics is conjectured to undergo a critical slowdown: At high temperature the inverse-gap is O(1), at the critical $β_c$ it is polynomial in the side-length and at low temperature it is exponential in it. A seminal series of papers verified this on $\Z^2$ except at $β=β_c$ where the behavior remained a challenging open problem. Here we establish the first rigorous polynomial upper bound for the critical mixing, thus confirming the critical slowdown for the Ising model in $\Z^2$. Namely, we show that on a finite box with arbitrary (e.g. fixed, free, periodic) boundary conditions, the inverse-gap at $β=β_c$ is polynomial in the side-length. The proof harnesses recent understanding of the scaling limit of critical Fortuin-Kasteleyn representation of the Ising model together with classical tools from the analysis of Markov chains.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Diameters in supercritical random graphs via first passage percolation

We study the diameter of $C_1$, the largest component of the Erdős-Rényi random graph $G(n,p)$ in the emerging supercritical phase, i.e., for $p = \frac{1+ε}n$ where $ε^3 n \to \infty$ and $ε=o(1)$. This parameter was extensively studied for fixed $ε> 0$, yet results for $ε=o(1)$ outside the critical window were only obtained very recently. Prior to this work, Riordan and Wormald gave precise estimates on the diameter, however these did not cover the entire supercritical regime (namely, when $ε^3 n\to\infty$ arbitrarily slowly). Łuczak and Seierstad estimated its order throughout this regime, yet their upper and lower bounds differed by a factor of $1000/7$. We show that throughout the emerging supercritical phase, i.e. for any $ε=o(1)$ with $ε^3 n \to \infty$, the diameter of $C_1$ is with high probability asymptotic to $D(ε,n)=(3/ε)\log(ε^3 n)$. This constitutes the first proof of the asymptotics of the diameter valid throughout this phase. The proof relies on a recent structure result for the supercritical giant component, which reduces the problem of estimating distances between its vertices to the study of passage times in first-passage percolation. The main advantage of our method is its flexibility. It also implies that in the emerging supercritical phase the diameter of the 2-core of $C_1$ is w.h.p. asymptotic to $(2/3)D(ε,n)$, and the maximal distance in $C_1$ between any pair of kernel vertices is w.h.p. asymptotic to $(5/9)D(ε,n)$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Explicit expanders with cutoff phenomena

The cutoff phenomenon describes a sharp transition in the convergence of an ergodic finite Markov chain to equilibrium. Of particular interest is understanding this convergence for the simple random walk on a bounded-degree expander graph. The first example of a family of bounded-degree graphs where the random walk exhibits cutoff in total-variation was provided only very recently, when the authors showed this for a typical random regular graph. However, no example was known for an explicit (deterministic) family of expanders with this phenomenon. Here we construct a family of cubic expanders where the random walk from a worst case initial position exhibits total-variation cutoff. Variants of this construction give cubic expanders without cutoff, as well as cubic graphs with cutoff at any prescribed time-point.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimal whitespace synchronization strategies

The whitespace-discovery problem describes two parties, Alice and Bob, trying to establish a communication channel over one of a given large segment of whitespace channels. Subsets of the channels are occupied in each of the local environments surrounding Alice and Bob, as well as in the global environment between them (Eve). In the absence of a common clock for the two parties, the goal is to devise time-invariant (stationary) strategies minimizing the synchronization time. This emerged from recent applications in discovery of wireless devices. We model the problem as follows. There are $N$ channels, each of which is open (unoccupied) with probability $p_1,p_2,q$ independently for Alice, Bob and Eve respectively. Further assume that $N \gg 1/(p_1 p_2 q)$ to allow for sufficiently many open channels. Both Alice and Bob can detect which channels are locally open and every time-slot each of them chooses one such channel for an attempted sync. One aims for strategies that, with high probability over the environments, guarantee a shortest possible expected sync time depending only on the $p_i$'s and $q$. Here we provide a stationary strategy for Alice and Bob with a guaranteed expected sync time of $O(1 / (p_1 p_2 q^2))$ given that each party also has knowledge of $p_1,p_2,q$. When the parties are oblivious of these probabilities, analogous strategies incur a cost of a poly-log factor, i.e.\ $\tilde{O}(1 / (p_1 p_2 q^2))$. Furthermore, this performance guarantee is essentially optimal as we show that any stationary strategies of Alice and Bob have an expected sync time of at least $Ω(1/(p_1 p_2 q^2))$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quasi-polynomial mixing of the 2D stochastic Ising model with "plus" boundary up to criticality

We considerably improve upon the recent result of Martinelli and Toninelli on the mixing time of Glauber dynamics for the 2D Ising model in a box of side $L$ at low temperature and with random boundary conditions whose distribution $P$ stochastically dominates the extremal plus phase. An important special case is when $P$ is concentrated on the homogeneous all-plus configuration, where the mixing time $T_{mix}$ is conjectured to be polynomial in $L$. In [MT] it was shown that for a large enough inverse-temperature $β$ and any $ε>0$ there exists $c=c(β,ε)$ such that $\lim_{L\to\infty}P(T_{mix}\geq \exp({c L^ε}))=0$. In particular, for the all-plus boundary conditions and $β$ large enough $T_{mix} \leq \exp({c L^ε})$. Here we show that the same conclusions hold for all $β$ larger than the critical value $β_c$ and with $\exp({c L^ε})$ replaced by $L^{c \log L}$ (i.e. quasi-polynomial mixing). The key point is a modification of the inductive scheme of [MT] together with refined equilibrium estimates that hold up to criticality, obtained via duality and random-line representation tools for the Ising model. In particular, we establish new precise bounds on the law of Peierls contours which quantitatively sharpen the Brownian bridge picture established e.g. in [Greenberg-Ioffe (2005)],[Higuchi (1979)],[Hryniv (1998)].

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cutoff for the Ising model on the lattice

Introduced in 1963, Glauber dynamics is one of the most practiced and extensively studied methods for sampling the Ising model on lattices. It is well known that at high temperatures, the time it takes this chain to mix in $L^1$ on a system of size $n$ is $O(\log n)$. Whether in this regime there is cutoff, i.e. a sharp transition in the $L^1$-convergence to equilibrium, is a fundamental open problem: If so, as conjectured by Peres, it would imply that mixing occurs abruptly at $(c+o(1))\log n$ for some fixed $c>0$, thus providing a rigorous stopping rule for this MCMC sampler. However, obtaining the precise asymptotics of the mixing and proving cutoff can be extremely challenging even for fairly simple Markov chains. Already for the one-dimensional Ising model, showing cutoff is a longstanding open problem. We settle the above by establishing cutoff and its location at the high temperature regime of the Ising model on the lattice with periodic boundary conditions. Our results hold for any dimension and at any temperature where there is strong spatial mixing: For $\Z^2$ this carries all the way to the critical temperature. Specifically, for fixed $d\geq 1$, the continuous-time Glauber dynamics for the Ising model on $(\Z/n\Z)^d$ with periodic boundary conditions has cutoff at $(d/2λ_\infty)\log n$, where $λ_\infty$ is the spectral gap of the dynamics on the infinite-volume lattice. To our knowledge, this is the first time where cutoff is shown for a Markov chain where even understanding its stationary distribution is limited. The proof hinges on a new technique for translating $L^1$ to $L^2$ mixing which enables the application of log-Sobolev inequalities. The technique is general and carries to other monotone and anti-monotone spin-systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Mixing time of critical Ising model on trees is polynomial in the height

In the heat-bath Glauber dynamics for the Ising model on the lattice, physicists believe that the spectral gap of the continuous-time chain exhibits the following behavior. For some critical inverse-temperature $β_c$, the inverse-gap is bounded for $β< β_c$, polynomial in the surface area for $β= β_c$ and exponential in it for $β> β_c$. This has been proved for $\Z^2$ except at criticality. So far, the only underlying geometry where the critical behavior has been confirmed is the complete graph. Recently, the dynamics for the Ising model on a regular tree, also known as the Bethe lattice, has been intensively studied. The facts that the inverse-gap is bounded for $β< β_c$ and exponential for $β> β_c$ were established, where $β_c$ is the critical spin-glass parameter, and the tree-height $h$ plays the role of the surface area. In this work, we complete the picture for the inverse-gap of the Ising model on the $b$-ary tree, by showing that it is indeed polynomial in $h$ at criticality. The degree of our polynomial bound does not depend on $b$, and furthermore, this result holds under any boundary condition. We also obtain analogous bounds for the mixing-time of the chain. In addition, we study the near critical behavior, and show that for $β> β_c$, the inverse-gap and mixing-time are both $\exp[Θ((β-β_c) h)]$.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Censored Glauber Dynamics for the mean field Ising Model

We study Glauber dynamics for the Ising model on the complete graph on $n$ vertices, known as the Curie-Weiss Model. It is well known that at high temperature ($β< 1$) the mixing time is $Θ(n\log n)$, whereas at low temperature ($β> 1$) it is $\exp(Θ(n))$. Recently, Levin, Luczak and Peres considered a censored version of this dynamics, which is restricted to non-negative magnetization. They proved that for fixed $β> 1$, the mixing-time of this model is $Θ(n\log n)$, analogous to the high-temperature regime of the original dynamics. Furthermore, they showed \emph{cutoff} for the original dynamics for fixed $β<1$. The question whether the censored dynamics also exhibits cutoff remained unsettled. In a companion paper, we extended the results of Levin et al. into a complete characterization of the mixing-time for the Currie-Weiss model. Namely, we found a scaling window of order $1/\sqrt{n}$ around the critical temperature $β_c=1$, beyond which there is cutoff at high temperature. However, determining the behavior of the censored dynamics outside this critical window seemed significantly more challenging. In this work we answer the above question in the affirmative, and establish the cutoff point and its window for the censored dynamics beyond the critical window, thus completing its analogy to the original dynamics at high temperature. Namely, if $β= 1 + δ$ for some $δ> 0$ with $δ^2 n \to \infty$, then the mixing-time has order $(n / δ)\log(δ^2 n)$. The cutoff constant is $(1/2+[2(ζ^2 β/ δ- 1)]^{-1})$, where $ζ$ is the unique positive root of $g(x)=\tanh(βx)-x$, and the cutoff window has order $n / δ$.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Non-backtracking random walks mix faster

We compute the mixing rate of a non-backtracking random walk on a regular expander. Using some properties of Chebyshev polynomials of the second kind, we show that this rate may be up to twice as fast as the mixing rate of the simple random walk. The closer the expander is to a Ramanujan graph, the higher the ratio between the above two mixing rates is. As an application, we show that if $G$ is a high-girth regular expander on $n$ vertices, then a typical non-backtracking random walk of length $n$ on $G$ does not visit a vertex more than $(1+o(1))\frac{\log n}{\log\log n}$ times, and this result is tight. In this sense, the multi-set of visited vertices is analogous to the result of throwing $n$ balls to $n$ bins uniformly, in contrast to the simple random walk on $G$, which almost surely visits some vertex $Ω(\log n)$ times.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Shannon capacity of a graph and the independence numbers of its powers

The independence numbers of powers of graphs have been long studied, under several definitions of graph products, and in particular, under the strong graph product. We show that the series of independence numbers in strong powers of a fixed graph can exhibit a complex structure, implying that the Shannon Capacity of a graph cannot be approximated (up to a sub-polynomial factor of the number of vertices) by any arbitrarily large, yet fixed, prefix of the series. This is true even if this prefix shows a significant increase of the independence number at a given power, after which it stabilizes for a while.

0 citations0 reviewsNo code linkedOpen access