Source author record

Giambattista Giacomin

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

20 published item(s)

preprint2022arXiv

Lyapunov exponent for products of random Ising transfer matrices: the balanced disorder case

We analyze the top Lyapunov exponent of the product of sequences of two by two matrices that appears in the analysis of several statistical mechanics models with disorder: for example these matrices are the transfer matrices for the nearest neighbor Ising chain with random external field, and the free energy density of this Ising chain is the Lyapunov exponent we consider. We obtain the sharp behavior of this exponent in the large interaction limit when the external field is centered: this balanced case turns out to be critical in many respects. From a mathematical standpoint we precisely identify the behavior of the top Lyapunov exponent of a product of two dimensional random matrices close to a diagonal random matrix for which top and bottom Lyapunov exponents coincide. In particular, the Lyapunov exponent is only $\log$-Hölder continuous.

preprint2022arXiv

The zeros of the partition function of the pinning model

We aim at understanding for which (complex) values of the potential the pinning partition function vanishes. The pinning model is a Gibbs measure based on discrete renewal processes with power law inter-arrival distributions. We obtain some results for rather general inter-arrival laws, but we achieve a substantially more complete understanding for a specific one parameter family of inter-arrivals. We show, for such a specific family, that the zeros asymptotically lie on (and densely fill) a closed curve that, unsurprisingly, touches the real axis only in one point (the critical point of the model). We also perform a sharper analysis of the zeros close to the critical point and we exploit this analysis to approach the challenging problem of Griffiths singularities for the disordered pinning model. The techniques we exploit are both probabilistic and analytical. Regarding the first, a central role is played by limit theorems for heavy tail random variables. As for the second, potential theory and singularity analysis of generating functions, along with their interplay, will be at the heart of several of our arguments.

preprint2018arXiv

Continuum limit of random matrix products in statistical mechanics of disordered systems

We consider a particular weak disorder limit ("continuum limit") of matrix products that arise in the analysis of disordered statistical mechanics systems, with a particular focus on random transfer matrices. The limit system is a diffusion model for which the leading Lyapunov exponent can be expressed explicitly in terms of modified Bessel functions, a formula that appears in the physical literature on these disordered systems. We provide an analysis of the diffusion system as well as of the link with the matrix products. We then apply the results to the framework considered by Derrida and Hilhorst [J. Phys. A (1983)], which deals in particular with the strong interaction limit for disordered Ising model in one dimension and that identifies a singular behavior of the Lyapunov exponent (of the transfer matrix), and to the two dimensional Ising model with columnar disorder (McCoy-Wu model). We show that the continuum limit sharply captures the Derrida and Hilhorst singularity. Moreover we revisit the analysis by McCoy and Wu [Phys. Rev. 1968] and remark that it can be interpreted in terms of the continuum limit approximation. We provide a mathematical analysis of the continuum approximation of the free energy of the McCoy-Wu model, clarifying the prediction (by McCoy and Wu) that, in this approximation, the free energy of the two dimensional Ising model with columnar disorder is $C^\infty$ but not analytic at the critical temperature.

preprint2016arXiv

A note on dynamical models on random graphs and Fokker-Planck equations

We address the issue of the proximity of interacting diffusion models on large graphs with a uniform degree property and a corresponding mean field model, i.e. a model on the complete graph with a suitably renormalized interaction parameter. Examples include Erdős-Rényi graphs with edge probability $p_n$, $n$ is the number of vertices, such that $\lim_{n \to \infty}p_n n= \infty$. The purpose of this note it twofold: (1) to establish this proximity on finite time horizon, by exploiting the fact that both systems are accurately described by a Fokker-Planck PDE (or, equivalently, by a nonlinear diffusion process) in the $n=\infty$ limit; (2) to remark that in reality this result is unsatisfactory when it comes to applying it to systems with $N$ large but finite, for example the values of $N$ that can be reached in simulations or that correspond to the typical number of interacting units in a biological system.

preprint2016arXiv

Small noise and long time phase diffusion in stochastic limit cycle oscillators

We study the effect of additive Brownian noise on an ODE system that has a stable hyperbolic limit cycle, for initial data that are attracted to the limit cycle. The analysis is performed in the limit of small noise - that is, we modulate the noise by a factor $\varepsilon \searrow 0$ - and on a long time horizon. We prove explicit estimates on the proximity of the noisy trajectory and the limit cycle up to times $\exp\left(c \varepsilon^{-2}\right)$, $c>0$, and we show both that on the time scale $\varepsilon^{-2}$ the "'dephasing" (i.e., the difference between noiseless and noisy system measured in a natural coordinate system that involves a phase) is close to a Brownian motion with constant drift, and that on longer time scales the dephasing dynamics is dominated, to leading order, by the drift. The natural choice of coordinates, that reduces the dynamics in a neighborhood of the cycle to a rotation, plays a central role and makes the connection with the applied science literature in which noisy limit cycle dynamics are often reduced to a diffusion model for the phase of the limit cycle.

preprint2015arXiv

Generalized Poland-Scheraga denaturation model and two-dimensional renewal processes

The Poland-Scheraga model describes the denaturation transition of two complementary - in particular, equally long - strands of DNA, and it has enjoyed a remarkable success both for quantitative modeling purposes and at a more theoretical level. The solvable character of the homogeneous version of the model is one of features to which its success is due. In the bio-physical literature a generalization of the model, allowing different length and non complementarity of the strands, has been considered and the solvable character extends to this substantial generalization. We present a mathematical analysis of the homogeneous generalized Poland-Scheraga model. Our approach is based on the fact that such a model is a homogeneous pinning model based on a bivariate renewal process, much like the basic Poland-Scheraga model is a pinning model based on a univariate, i.e. standard, renewal. We present a complete analysis of the free energy singularities, which include the localization-delocalization critical point and (in general) other critical points that have been only partially captured in the physical literature. We obtain also precise estimates on the path properties of the model.

preprint2015arXiv

Pinning and disorder relevance for the lattice Gaussian free field

This paper provides a rigorous study of the localization transition for a Gaussian free field on $\mathbb{Z}^d$ interacting with a quenched disordered substrate that acts on the interface when the interface height is close to zero. The substrate has the tendency to localize or repel the interface at different sites and one can show that a localization-delocalization transition takes place when varying the average pinning potential $h$: the free energy density is zero in the delocalized regime, that is for $h$ smaller than a threshold $h_c$, and it is positive for $h>h_c$. For $d\ge 3$ we compute $h_c$ and we show that the transition happens at the same value as for the annealed model. However we can show that the critical behavior of the quenched model differs from the one of the annealed one. While the phase transition of the annealed model is of first order, we show that the quenched free energy is bounded above by $ (h-h_c)_+^2$ times a positive constant and that, for Gaussian disorder, the quadratic behavior is sharp. Therefore this provides an example in which a {\sl relevant disorder critical exponent} can be made explicit: in theoretical physics disorder is said to be {\sl relevant} when the disorder changes the critical behavior of a system and, while there are cases in which it is known that disorder is relevant, the exact critical behavior is typically unknown. For $d=2$ we are not able to decide whether the quenched and annealed critical points coincide, but we provide an upper bound for the difference between them.

preprint2015arXiv

Weak noise and non-hyperbolic unstable fixed points: Sharp estimates on transit and exit times

We consider certain one dimensional ordinary stochastic differential equations driven by additive Brownian motion of variance $\varepsilon ^2$. When $\varepsilon =0$ such equations have an unstable non-hyperbolic fixed point and the drift near such a point has a power law behavior. For $\varepsilon >0$ small, the fixed point property disappears, but it is replaced by a random escape or transit time which diverges as $\varepsilon \searrow0$. We show that this random time, under suitable (easily guessed) rescaling, converges to a limit random variable that essentially depends only on the power exponent associated to the fixed point. Such random variables, or laws, have therefore a universal character and they arise of course in a variety of contexts. We then obtain quantitative sharp estimates, notably tail properties, on these universal laws.

preprint2014arXiv

Landau damping in the Kuramoto model

We consider the Kuramoto model of globally coupled phase oscillators in its continuum limit, with individual frequencies drawn from a distribution with density of class $C^n$ ($n\geq 4$). A criterion for linear stability of the uniform stationary state is established which, for basic examples of frequency distributions, is equivalent to the standard condition on the coupling strength in the literature. We prove that, under this criterion, the Kuramoto order parameter, when evolved under the full nonlinear dynamics, asymptotically vanishes (with polynomial rate $n$) for every trajectory issued from sufficiently small $C^n$ perturbation. The proof uses techniques from the Analysis of PDEs and closely follows recent proofs of the nonlinear Landau damping in the Vlasov equation and Vlasov-HMF model.

preprint2013arXiv

Log-periodic critical amplitudes: a perturbative approach

Log-periodic amplitudes appear in the critical behavior of a large class of systems, in particular when a discrete scale invariance is present. Here we show how to compute these critical amplitudes perturbatively when they originate from a renormalization map which is close to a monomial. In this case, the log-periodic amplitudes of the subdominant corrections to the leading critical behavior can also be calculated.

preprint2012arXiv

Global attractor and asymptotic dynamics in the Kuramoto model for coupled noisy phase oscillators

We study the dynamics of the large N limit of the Kuramoto model of coupled phase oscillators, subject to white noise. We introduce the notion of shadow inertial manifold and we prove their existence for this model, supporting the fact that the long term dynamics of this model is finite dimensional. Following this, we prove that the global attractor of this model takes one of two forms. When coupling strength is below a critical value, the global attractor is a single equilibrium point corresponding to an incoherent state. Conversely, when coupling strength is beyond this critical value, the global attractor is a two-dimensional disk composed of radial trajectories connecting a saddle equilibrium (the incoherent state) to an invariant closed curve of locally stable equilibria (partially synchronized state). Our analysis hinges, on the one hand, upon sharp existence and uniqueness results and their consequence for the existence of a global attractor, and, on the other hand, on the study of the dynamics in the vicinity of the incoherent and synchronized equilibria. We prove in particular non-linear stability of each synchronized equilibrium, and normal hyperbolicity of the set of such equilibria. We explore mathematically and numerically several properties of the global attractor, in particular we discuss the limit of this attractor as noise intensity decreases to zero.

preprint2012arXiv

Oscillatory critical amplitudes in hierarchical models and the tail of the Harris random variable

Oscillatory critical amplitudes have been repeatedly observed in hierarchical models and, in the cases that have been taken into consideration, these oscillations are so small to be hardly detectable. Hierarchical models are tightly related to iteration of maps and, in fact, very similar phenomena have been repeatedly reported in many fields of mathematics, like combinatorial evaluations and discrete branching processes. It is precisely in the context of branching processes with bounded off-spring that T. Harris, in 1948, first set forth the possibility that the logarithm of the moment generating function of the rescaled population size, in the super-critical regime, does not grow near infinity as a power, but it has an oscillatory prefactor. These oscillations have been observed numerically only much later and, while the origin is clearly tied to the discrete character of the iteration, the amplitude size is not so well understood. The purpose of this note is to reconsider the issue for hierarchical models and in what is arguably the most elementary setting -- the pinning model -- that actually just boils down to iteration of polynomial maps (and, notably, quadratic maps). In this note we show that the oscillatory critical amplitude for pinning models and the oscillating pre factor connected to the Harris random variable coincide. Moreover we make explicit the link between these oscillatory functions and the geometry of the Julia set of the map, making thus rigorous and quantitative some ideas set forth by B. Derrida, C. Itzykson and J. M. Luck (CMP 1984).

preprint2012arXiv

Synchronization and random long time dynamics for mean-field plane rotators

We consider the natural Langevin dynamics which is reversible with respect to the mean-field plane rotator (or classical spin XY) measure. It is well known that this model exhibits a phase transition at a critical value of the interaction strength parameter K, in the limit of the number N of rotators going to infinity. A Fokker-Planck PDE captures the evolution of the empirical measure of the system as N goes to infinity, at least for finite times and when the empirical measure of the system at time zero satisfies a law of large numbers. The phase transition is reflected in the fact that the PDE for K above the critical value has a stable manifold of stationary solutions, that are equivalent up to rotations. These stationary solutions are actually unimodal densities parametrized by the position of their maximum (the synchronization phase or center). We characterize the dynamics on times of order N and we show substantial deviations from the behavior of the solutions of the PDE. In fact, if the empirical measure at time zero converges as N goes to infinity to a probability measure (which is away from a thin set that we characterize) and if time is speeded up by N, the empirical measure reaches almost instantaneously a small neighborhood of the stable manifold, to which it then sticks and on which a non-trivial random dynamics takes place. In fact the synchronization center performs a Brownian motion with a diffusion coefficient that we compute. Our approach therefore provides, for one of the basic statistical mechanics systems with continuum symmetry, a detailed characterization of the macroscopic deviations from the large scale limit -- or law of large numbers -- due to finite size effects. But the interest for this model goes beyond statistical mechanics, since it plays a central role in a variety of scientific domains in which one aims at understanding synchronization phenomena.

preprint2011arXiv

Coherence stability and effect of random natural frequencies in populations of coupled oscillators

We consider the (noisy) Kuramoto model, that is a population of N oscillators, or rotators, with mean-field interaction. Each oscillator has its own randomly chosen natural frequency (quenched disorder) and it is stirred by Brownian motion. In the limit N goes to infty this model is accurately described by a (deterministic) Fokker-Planck equation. We study this equation and obtain quantitatively sharp results in the limit of weak disorder. We show that, in general, even when the natural frequencies have zero mean the oscillators synchronize (for sufficiently strong interaction) around a common rotating phase, whose frequency is sharply estimated. We also establish the stability properties of these solutions (in fact, limit cycles). These results are obtained by identifying the stable hyperbolic manifold of stationary solutions of an associated non disordered model and by exploiting the robustness of hyperbolic structures under suitable perturbations. When the disorder distribution is symmetric the speed vanishes and there is a one parameter family of stationary solutions : in this case we provide more precise stability estimates. The methods we use apply beyond the Kuramoto model and we develop here the case of active rotator models, that is the case in which the dynamics of each rotator in absence of interaction and noise is not simply a rotation.

preprint2011arXiv

Transitions in active rotator systems: invariant hyperbolic manifold approach

Our main focus is on a general class of active rotators with mean field interactions, that is globally coupled large families of dynamical systems on the unit circle with non-trivial stochastic dynamics. Each isolated system is a diffusion process on a circle, with drift -delta V', where V' is a periodic function and delta is an intensity parameter. It is well known that the interacting dynamics is accurately described, in the limit of infinitely many interacting components, by a Fokker-Planck PDE and the model reduces for delta=0 to a particular case of the Kuramoto synchronization model, for which one can show the existence of a stable normally hyperbolic manifold of stationary solutions for the corresponding Fokker-Planck equation (we are interested in the case in which this manifold is non-trivial, that happens when the interaction is sufficiently strong, that is in the synchronized regime of the Kuramoto model). We use the robustness of normally hyperbolic structures to infer qualitative and quantitative results on the |delta|< delta0 cases, with delta0 a suitable threshold: as a matter of fact, we obtain an accurate description of the dynamics on the invariant manifold for delta=0 and we link it explicitly to the potential V . This approach allows to have a complete description of the phase diagram of the active rotators model, at least for |delta|< delta0, thus identifying for which values of the parameters (notably, noise intensity and/or coupling strength) the system exhibits periodic pulse waves or stabilizes at a quiescent resting state. Moreover, some of our results are very explicit and this brings a new insight into the combined effect of active rotator dynamics, noise and interaction. The links with the literature on specific systems, notably neuronal models, are discussed in detail.

preprint2010arXiv

The weak coupling limit of disordered copolymer models

A copolymer is a chain of repetitive units (monomers) that are almost identical, but they differ in their degree of affinity for certain solvents. This difference leads to striking phenomena when the polymer fluctuates in a nonhomogeneous medium, for example, made of two solvents separated by an interface. One may observe, for instance, the localization of the polymer at the interface between the two solvents. A discrete model of such system, based on the simple symmetric random walk on $\mathbb {Z}$, has been investigated in [Bolthausen and den Hollander, Ann. Probab. 25 (1997), 1334-1366], notably in the weak polymer-solvent coupling limit, where the convergence of the discrete model toward a continuum model, based on Brownian motion, has been established. This result is remarkable because it strongly suggests a universal feature of copolymer models. In this work, we prove that this is indeed the case. More precisely, we determine the weak coupling limit for a general class of discrete copolymer models, obtaining as limits a one-parameter [$α\in(0,1)$] family of continuum models, based on $α$-stable regenerative sets.

preprint2009arXiv

Dynamical aspects of mean field plane rotators and the Kuramoto model

The Kuramoto model has been introduced in order to describe synchronization phenomena observed in groups of cells, individuals, circuits, etc... We look at the Kuramoto model with white noise forces: in mathematical terms it is a set of N oscillators, each driven by an independent Brownian motion with a constant drift, that is each oscillator has its own frequency, which, in general, changes from one oscillator to another (these frequencies are usually taken to be random and they may be viewed as a quenched disorder). The interactions between oscillators are of long range type (mean field). We review some results on the Kuramoto model from a statistical mechanics standpoint: we give in particular necessary and sufficient conditions for reversibility and we point out a formal analogy, in the N to infinity limit, with local mean field models with conservative dynamics (an analogy that is exploited to identify in particular a Lyapunov functional in the reversible set-up). We then focus on the reversible Kuramoto model with sinusoidal interactions in the N to infinity limit and analyze the stability of the non-trivial stationary profiles arising when the interaction parameter K is larger than its critical value K_c. We provide an analysis of the linear operator describing the time evolution in a neighborhood of the synchronized profile: we exhibit a Hilbert space in which this operator has a self-adjoint extension and we establish, as our main result, a spectral gap inequality for every K>K_c.

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.