Source author record

Fabio Lucio Toninelli

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

ResearcherUnclaimed source record

Catalog footprint

Research overview

22works
4topics
4close collaborators

Published work

22 publications

preprint2020arXiv

Hydrodynamic limit for a 2D interlaced particle process

The Markov dynamics of interlaced particle arrays, introduced by A. Borodin and P. Ferrari in arXiv:0811.0682, is a classical example of (2+1)-dimensional random growth model belonging to the so-called Anisotropic KPZ universality class. In Legras-Toninelli (2017) arXiv:1704.06581, a hydrodynamic limit -- the convergence of the height profile, after space/time rescaling, to the solution of a deterministic Hamilton-Jacobi PDE with non-convex Hamiltonian -- was proven when either the initial profile is convex, or for small times, before the solution develops shocks. In the present work, we give a simpler proof, that works for all times and for all initial profiles for which the limit equation makes sense. In particular, the convexity assumption is dropped. The main new idea is a new viewpoint about "finite speed of propagation" that allows to bypass the need of a-priori control of the interface gradients, or equivalently of inter-particle distances.

preprint2020arXiv

Non-integrable dimers: Universal fluctuations of tilted height profiles

We study a class of close-packed dimer models on the square lattice, in the presence of small but extensive perturbations that make them non-determinantal. Examples include the 6-vertex model close to the free-fermion point, and the dimer model with plaquette interaction previously analyzed in \cite{A,AL,GMT17a,GMT17b}. By tuning the edge weights, we can impose a non-zero average tilt for the height function, so that the considered models are in general not symmetric under discrete rotations and reflections. In the determinantal case, height fluctuations in the massless (or `liquid') phase scale to a Gaussian log-correlated field and their amplitude is a universal constant, independent of the tilt. When the perturbation strength λλ is sufficiently small we prove, by fermionic constructive Renormalization Group methods, that log-correlations survive, with amplitude AA that, generically, depends non-trivially and non-universally on λλ and on the tilt. On the other hand, AA satisfies a universal scaling relation (`Haldane' or `Kadanoff' relation), saying that it equals the anomalous exponent of the dimer-dimer correlation.

preprint2018arXiv

A (2+1)-dimensional Anisotropic KPZ growth model with a smooth phase

Stochastic growth processes in dimension (2+1)(2+1) were conjectured by D. Wolf, on the basis of renormalization-group arguments, to fall into two distinct universality classes, according to whether the Hessian HρH_ρ of the speed of growth v(ρ)v(ρ) as a function of the average slope ρρ satisfies det⁡Hρ>0\det H_ρ>0 ("isotropic KPZ class") or det⁡Hρ≤0\det H_ρ\le 0 ("anisotropic KPZ (AKPZ)" class). The former is characterized by strictly positive growth and roughness exponents, while in the AKPZ class fluctuations are logarithmic in time and space. It is natural to ask (a) if one can exhibit interesting growth models with "smooth" stationary states, i.e., with O(1)O(1) fluctuations (instead of logarithmically or power-like growing, as in Wolf's picture) and (b) what new phenomena arise when v(⋅)v(\cdot) is not smooth, so that HρH_ρ is not defined. The two questions are actually related and here we provide an answer to both, in a specific framework. We define a (2+1)(2+1)-dimensional interface growth process, based on the so-called shuffling algorithm for domino tilings. The stationary, non-reversible measures are translation-invariant Gibbs measures on perfect matchings of Z2\mathbb Z^2, with 22-periodic weights. If ρ≠0ρ\ne0, fluctuations are known to grow logarithmically in space and to behave like a two-dimensional GFF. We prove that fluctuations grow at most logarithmically in time and that det⁡Hρ<0\det H_ρ<0: the model belongs to the AKPZ class. When ρ=0ρ=0, instead, the stationary state is "smooth", with correlations uniformly bounded in space and time; correspondingly, v(⋅)v(\cdot) is not differentiable at ρ=0ρ=0 and we extract the singularity of the eigenvalues of HρH_ρ for ρ∼0ρ\sim 0.

preprint2015arXiv

Height fluctuations in non-integrable classical dimers

We rigorously establish the asymptotic equivalence between the height function of interacting dimers on the square lattice and the massless Gaussian free field. Our theorem explains the microscopic origin of the sine-Gordon field theory description away from the free fermion point, which has previously been elusive. We use a novel technique, based on the combination of discrete holomorphicity with exact, constructive, renormalization group methods, which has the potential of being applicable to a variety of other non-integrable models at or close to criticality.

preprint2015arXiv

On the probability of staying above a wall for the (2+1)-dimensional SOS model at low temperature

We obtain sharp asymptotics for the probability that the (2+1)-dimensional discrete SOS interface at low temperature is positive in a large region. For a square region ΛΛ, both under the infinite volume measure and under the measure with zero boundary conditions around ΛΛ, this probability turns out to behave like exp⁡(−τβ(0)Llog⁡L)\exp(-τ_β(0) L \log L ), with τβ(0)τ_β(0) the surface tension at zero tilt, also called step free energy, and LL the box side. This behavior is qualitatively different from the one found for continuous height massless gradient interface models.

preprint2015arXiv

Stationary measure of the driven two-dimensional q-Whittaker particle system on the torus

We consider a q-deformed version of the uniform Gibbs measure on dimers on the periodized hexagonal lattice (equivalently, on interlacing particle configurations, if vertical dimers are seen as particles) and show that it is invariant under a certain irreversible q-Whittaker dynamic. Thereby we provide a new non-trivial example of driven interacting two-dimensional particle system, or of (2+1)-dimensional stochastic growth model, with explicit stationary measure. We emphasize that this measure is far from being a product Bernoulli measure. These Gibbs measures and dynamics both arose earlier in the theory of Macdonald processes. The q=0 degeneration of the Gibbs measures reduce to the usual uniform dimer measures with given tilt, the degeneration of the dynamics originate in the study of Schur processes and the degeneration of the results contained herein were recently treated in work of the second author.

preprint2014arXiv

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

We study the Glauber dynamics for the (2+1)D(2+1)\mathrm{D} Solid-On-Solid model above a hard wall and below a far away ceiling, on an L×LL\times L box of Z2\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≍(1/β)log⁡LH\asymp(1/β)\log L. As an essential step in understanding the effect of entropic repulsion on the Glauber dynamics we determine the equilibrium height HH to within an additive constant: H=(1/4β)log⁡L+O(1)H=(1/4β)\log L+O(1). We then show that starting from zero initial conditions the surface rises to its final height HH through a sequence of metastable transitions between consecutive levels. The time for a transition from height h=aHh=aH, a∈(0,1)a\in(0,1), to height h+1h+1 is roughly exp⁡(cLa)\exp(cL^a) for some constant c>0c>0. In particular, the mixing time of the dynamics is exponentially large in LL, that is, TMIX≥ecLT_{\mathrm{MIX}}\geq e^{cL}. We also provide the matching upper bound TMIX≤ec′LT_{\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 LL.

preprint2014arXiv

Interaction versus entropic repulsion for low temperature Ising polymers

Contours associated to many interesting low-temperature statistical mechanics models (2D Ising model, (2+1)D SOS interface model, etc) can be described as self-interacting and self-avoiding walks on Z2\mathbb Z^2. When the model is defined in a finite box, the presence of the boundary induces an interaction, that can turn out to be attractive, between the contour and the boundary of the box. On the other hand, the contour cannot cross the boundary, so it feels entropic repulsion from it. In various situations of interest a crucial technical problem is to prove that entropic repulsion prevails over the pinning interaction: in particular, the contour-boundary interaction should not modify significantly the contour partition function and the related surface tension should be unchanged. Here we prove that this is indeed the case, at least at sufficiently low temperature, in a quite general framework that applies in particular to the models of interest mentioned above.

preprint2014arXiv

Multi-level pinning problems for random walks and self-avoiding lattice paths

We consider a generalization of the classical pinning problem for integer-valued random walks conditioned to stay non-negative. More specifically, we take pinning potentials of the form ∑j≥0εjNj\sum_{j\geq 0}ε_j N_j, where NjN_j is the number of visits to the state jj and {εj}\{ε_j\} is a non-negative sequence. Partly motivated by similar problems for low-temperature contour models in statistical physics, we aim at finding a sharp characterization of the threshold of the wetting transition, especially in the regime where the variance σ2σ^2 of the single step of the random walk is small. Our main result says that, for natural choices of the pinning sequence {εj}\{ε_j\}, localization (respectively delocalization) occurs if σ−2∑j≥0(j+1)εj≥δ−1σ^{-2}\sum_{ j\geq0}(j+1)ε_j\geqδ^{-1} (respectively ≤δ\le δ), for some universal δ<1δ<1. Our finding is reminiscent of the classical Bargmann-Jost-Pais criteria for the absence of bound states for the radial Schrödinger equation. The core of the proof is a recursive argument to bound the free energy of the model. Our approach is rather robust, which allows us to obtain similar results in the case where the random walk trajectory is replaced by a self-avoiding path γγ in Z2\mathbb Z^2 with weight exp⁡(−β∣γ∣)\exp(-β|γ|), ∣γ∣|γ| being the length of the path and β>0β>0 a large enough parameter. This generalization is directly relevant for applications to the above mentioned contour models.

preprint2013arXiv

A one-dimensional coagulation-fragmentation process with a dynamical phase transition

We introduce a reversible Markovian coagulation-fragmentation process on the set of partitions of {1,…,L}\{1,\ldots,L\} into disjoint intervals. Each interval can either split or merge with one of its two neighbors. The invariant measure can be seen as the Gibbs measure for a homogeneous pinning model \cite{cf:GBbook}. Depending on a parameter λλ, the typical configuration can be either dominated by a single big interval (delocalized phase), or be composed of many intervals of order 11 (localized phase), or the interval length can have a power law distribution (critical regime). In the three cases, the time required to approach equilibrium (in total variation) scales very differently with LL. In the localized phase, when the initial condition is a single interval of size LL, the equilibration mechanism is due to the propagation of two "fragmentation fronts" which start from the two boundaries and proceed by power-law jumps.

preprint2013arXiv

Lozenge tilings, Glauber dynamics and macroscopic shape

We study the Glauber dynamics on the set of tilings of a finite domain of the plane with lozenges of side 1/L. Under the invariant measure of the process (the uniform measure over all tilings), it is well known that the random height function associated to the tiling converges in probability, in the scaling limit L→∞L\to\infty, to a non-trivial macroscopic shape minimizing a certain surface tension functional. According to the boundary conditions the macroscopic shape can be either analytic or contain "frozen regions" (Arctic Circle phenomenon). It is widely conjectured, on the basis of theoretical considerations, partial mathematical results and numerical simulations for similar models, that the Glauber dynamics approaches the equilibrium macroscopic shape in a time of order L2+o(1)L^{2+o(1)}. In this work we prove this conjecture, under the assumption that the macroscopic equilibrium shape contains no "frozen region".

preprint2013arXiv

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

Consider the classical (2+1)(2+1)-dimensional Solid-On-Solid model above a hard wall on an L×LL\times L box of $\bbZ^2$. The model describes a crystal surface by assigning a non-negative integer height ηxη_x to each site xx in the box and 0 heights to its boundary. The probability of a surface configuration ηη is proportional to exp⁡(−βH(η))\exp(-β\mathcal{H}(η)), where ββ is the inverse-temperature and H(η)\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)=⌊(1/4β)log⁡L⌋H(L)=\lfloor (1/4β)\log L\rfloor and H(L)−1H(L)-1. Moreover, for most values of LL the height is concentrated on the single value H(L)H(L). Next, we study the ensemble of level lines corresponding to the heights (H(L),H(L)−1,...)(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 LkL_k, the rescaled macroscopic level line at height H(Lk)−nH(L_k)-n has a limiting shape if the fractional parts of (1/4β)log⁡Lk(1/4β)\log L_k converge to a noncritical value. The scaling limit is an explicit convex subset of the unit square QQ and its boundary has a flat component on the boundary of QQ. Finally, the highest macroscopic level line has Lk1/3+o(1)L_k^{1/3+o(1)} fluctuations along the flat part of the boundary of its limiting shape.

preprint2012arXiv

The shape of the (2+1)(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×LL\times L box at a large inverse-temperature ββ the height of most sites concentrates on a single level h=⌊(1/4β)log⁡L⌋h = \lfloor (1/4β)\log L\rfloor for most values of LL. For a sequence of diverging boxes the ensemble of level lines of heights (h,h−1,...)(h,h-1,...) has a scaling limit in Hausdorff distance iff the fractional parts of (1/4β)log⁡L(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 hh-level lines feature L1/3+o(1)L^{1/3+o(1)} fluctuations from the side boundaries.

preprint2011arXiv

Mixing times of monotone surfaces and SOS interfaces: a mean curvature approach

We consider stochastic spin-flip dynamics for: (i) monotone discrete surfaces in Z^3 with planar boundary height and (ii) the one-dimensional discrete Solid-on-Solid (SOS) model confined to a box. In both cases we show almost optimal bounds O(L^2polylog(L)) for the mixing time of the chain, where L is the natural size of the system. The dynamics at a macroscopic scale should be described by a deterministic mean curvature motion such that each point of the surface feels a drift which tends to minimize the local surface tension. Inspired by this heuristics, our approach consists in bounding the dynamics with an auxiliary one which, with very high probability, follows quite closely the deterministic mean curvature evolution. Key technical ingredients are monotonicity, coupling and an argument due to D.B. Wilson in the framework of lozenge tiling Markov Chains. Our approach works equally well for both models despite the fact that their equilibrium maximal height fluctuations occur on very different scales (logarithmic for monotone surfaces and L^{1/2} for the SOS model). Finally, combining techniques from kinetically constrained spin systems together with the above mixing time result, we prove an almost diffusive lower bound of order 1/L^2 up to logarithmic corrections for the spectral gap of the SOS model with horizontal size L and unbounded heights.

preprint2010arXiv

"Zero" temperature stochastic 3D Ising model and dimer covering fluctuations: a first step towards interface mean curvature motion

We consider the Glauber dynamics for the Ising model with "+" boundary conditions, at zero temperature or at temperature which goes to zero with the system size (hence the quotation marks in the title). In dimension d=3 we prove that an initial domain of linear size L of "-" spins disappears within a time τ_+ which is at most L^2(\log L)^c and at least L^2/(c\log L), for some c>0. The proof of the upper bound proceeds via comparison with an auxiliary dynamics which mimics the motion by mean curvature that is expected to describe, on large time-scales, the evolution of the interface between "+" and "-" domains. The analysis of the auxiliary dynamics requires recent results on the fluctuations of the height function associated to dimer coverings of the infinite honeycomb lattice. Our result, apart from the spurious logarithmic factor, is the first rigorous confirmation of the expected behavior τ_+\simeq const\times L^2, conjectured on heuristic grounds. In dimension d=2, τ_+ can be shown to be of order L^2 without logarithmic corrections: the upper bound was proven in [Fontes, Schonmann, Sidoravicius, 2002] and here we provide the lower bound. For d=2, we also prove that the spectral gap of the generator behaves like c/L for L large, as conjectured in [Bodineau-Martinelli, 2002].

preprint2010arXiv

Convergence to equilibrium of biased plane partitions

We study a single-flip dynamics for the monotone surface in (2+1) dimensions obtained from a boxed plane partition. The surface is analyzed as a system of non-intersecting simple paths. When the flips have a non-zero bias we prove that there is a positive spectral gap uniformly in the boundary conditions and in the size of the system. Under the same assumptions, for a system of size M, the mixing time is shown to be of order M up to logarithmic corrections.

preprint2010arXiv

Polymer dynamics in the depinned phase: metastability with logarithmic barriers

We consider the stochastic evolution of a (1 + 1)-dimensional polymer in the depinned regime. At equilibrium the system exhibits a double well structure: the polymer lies(essentially) either above or below the repulsive line. As a consequence one expects a metastable behavior with rare jumps between the two phases combined with a fast thermalization inside each phase. However the energy barrier between these two phases is only logarithmic in the system size L and therefore the two relevant time scales are only polynomial in L with no clear-cut separation between them. The whole evolution is governed by a subtle competition between the diffusive behavior inside one phase and the jumps across the energy barriers. In particular the usual scenario in which the tunneling time coincides with the exponential of the energy barrier breaks down. Our main results are: (i) a proof that the mixing time of the system lies between L^{5/2} and L^{5/2+2}; (ii) the identification of two regions associated with the positive and negative phase of the polymer together with the proof of the asymptotic exponentiality of the tunneling time between them with rate equal to a half of the spectral gap.

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 LL at low temperature and with random boundary conditions whose distribution PP stochastically dominates the extremal plus phase. An important special case is when PP is concentrated on the homogeneous all-plus configuration, where the mixing time TmixT_{mix} is conjectured to be polynomial in LL. In [MT] it was shown that for a large enough inverse-temperature ββ and any ε>0ε>0 there exists c=c(β,ε)c=c(β,ε) such that lim⁡L→∞P(Tmix≥exp⁡(cLε))=0\lim_{L\to\infty}P(T_{mix}\geq \exp({c L^ε}))=0. In particular, for the all-plus boundary conditions and ββ large enough Tmix≤exp⁡(cLε)T_{mix} \leq \exp({c L^ε}). Here we show that the same conclusions hold for all ββ larger than the critical value βcβ_c and with exp⁡(cLε)\exp({c L^ε}) replaced by Lclog⁡LL^{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)].

preprint2009arXiv

Hierarchical pinning models, quadratic maps and quenched disorder

We consider a hierarchical model of polymer pinning in presence of quenched disorder, introduced by B. Derrida, V. Hakim and J. Vannimenius in 1992, which can be re-interpreted as an infinite dimensional dynamical system with random initial condition (the disorder). It is defined through a recurrence relation for the law of a random variable {R_n}_{n=1,2,...}, which in absence of disorder (i.e., when the initial condition is degenerate) reduces to a particular case of the well-known Logistic Map. The large-n limit of the sequence of random variables 2^{-n} log R_n, a non-random quantity which is naturally interpreted as a free energy, plays a central role in our analysis. The model depends on a parameter alpha>0, related to the geometry of the hierarchical lattice, and has a phase transition in the sense that the free energy is positive if the expectation of R_0 is larger than a certain threshold value, and it is zero otherwise. It was conjectured by Derrida et al. (1992) that disorder is relevant (respectively, irrelevant or marginally relevant) if 1/2<alpha<1 (respectively, alpha<1/2 or alpha=1/2), in the sense that an arbitrarily small amount of randomness in the initial condition modifies the critical point with respect to that of the pure (i.e., non-disordered) model if alpha is larger or equal to 1/2, but not if alpha is smaller than 1/2. Our main result is a proof of these conjectures for the case alpha different from 1/2. We emphasize that for alpha>1/2 we find the correct scaling form (for weak disorder) of the critical point shift.

preprint2009arXiv

Marginal relevance of disorder for pinning models

The effect of disorder on pinning and wetting models has attracted much attention in theoretical physics. In particular, it has been predicted on the basis of the Harris criterion that disorder is relevant (annealed and quenched model have different critical points and critical exponents) if the return probability exponent alpha, a positive number that characterizes the model, is larger than 1/2. Weak disorder has been predicted to be irrelevant (i.e. coinciding critical points and exponents) if alpha < 1/2. Recent mathematical work has put these predictions on firm grounds. In renormalization group terms, the case alpha = 1/2 is a 'marginal case' and there is no agreement in the literature as to whether one should expect disorder relevance or irrelevance at marginality. The question is particularly intriguing also because the case alpha = 1/2 includes the classical models of two-dimensional wetting of a rough substrate, of pinning of directed polymers on a defect line in dimension (3+1) or (1+1) and of pinning of an heteropolymer by a point potential in three-dimensional space. Here we prove disorder relevance both for the general alpha = 1/2 pinning model and for the hierarchical version of the model proposed by B. Derrida, V. Hakim and J. Vannimenus (JSP, 1992), in the sense that we prove a shift of the quenched critical point with respect to the annealed one. In both cases we work with Gaussian disorder and we show that the shift is at least of order exp(-1/β^4) for beta small, if beta is the standard deviation of the disorder.

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.