Removing Phase Transitions from Gibbs Measures

Total Page:16

File Type:pdf, Size:1020Kb

Removing Phase Transitions from Gibbs Measures Removing Phase Transitions from Gibbs Measures Ian E. Fellows Mark S. Handcock Fellows Statistics University of California, Los Angeles Abstract While the Gibbs measure has shown wide applicability as a general modeling framework, the presence of mul- Gibbs measures are a fundamental class of tiphase behavior has caused some practical difficulty. distributions for the analysis of high dimen- Nowhere is this more apparent than in the network sional data. Phase transitions, which are also science field, where the popular Exponential-family known as degeneracy in the network science Random Graph Model (ERGM) has been repeatedly literature, are an emergent property of these shown to suffer from model frailty due to the exis- models that well describe many physical sys- tence of multiphase parameter configurations. This is tems. However, the reach of the Gibbs mea- known as \model degeneracy" in the network science sure is now far outside the realm of physi- literature, and is perhaps the most fundamental chal- cal systems, and in many of these domains lenge facing ERGMs both in theory [Schweinberger, multiphase behavior is a nuisance. This nui- 2011, Handcock, 2003, Bhamidi et al., 2008, Chatter- sance often makes distribution fitting impos- jee et al., 2013] and practice [Robins et al., 2007, Gold- sible due to failure of the MCMC sampler, enberg et al., 2010]. and even when an MLE fit is possible, if the The contribution of this work is to introduce a new solution is near a phase transition point, the class of statistical models for high dimensional data plausibility of the fit can be highly question- that retains the flexibility of the Gibbs measure, with- able. We introduce a modification to the out the undesirable multiphase properties. The form Gibbs distribution that reduces the effects of this model class allows us to ensure that it is not of phase transitions, and with properly cho- multiphase at any parameter values, and that the like- sen hyper-parameters, provably removes all lihood puts significant probability mass around the multiphase behavior. We show that this new configuration that is being modeled. distribution is just as easy to fit via MCM- CMLE as the Gibbs measure, and provide Let X be a set of random variables with realization examples in the Ising model from statistical x and N be the sample space in which X may exist. physics and ERGMs from network science. A Gibbs measure is defined by the probability mass function 1 Introduction 1 P θ g (x) p (xjθ) = e k k k x 2 N gibbs Z(θ) The Gibbs measure is a class of probability distribu- where θ is a vector of d parameters, g is a vector valued P tions describing the configuration of a (usually high P i θigi(x) set of sufficient statistics, and Z(θ) = x e is dimensional) set of random variables. Initially it was the partition function assuring that the probabilities motivated primarily by statistical mechanics [Gibbs sum to unity. For notational convenience we will be et al., 1902], though in more recent decades it has assuming throughout that x is discrete. The continu- found widespread usage in many areas of statistics and ous case may be handled with ease by replacing sums machine learning, including computer vision [Li, 2012], over x by integrals. neuroscience [Stephens et al., 2011] and the modeling of social networks [Frank and Strauss, 1986] to name Since the likelihood is only a function of the sufficient a few. statistics, the probability of observing a set of statistics is th Proceedings of the 20 International Conference on Artifi- 1 P θkgk(x) cial Intelligence and Statistics (AISTATS) 2017, Fort Laud- pgibbs(g(x)jθ) = N(g(x)) e k x 2 N ; erdale, Florida, USA. JMLR: W&CP volume 54. Copy- Z(θ) right 2017 by the author(s). where N(g(x)) is the count of configurations with Removing Phase Transitions from Gibbs Measures statistics equal to g(x). disparate types of Gibbs measures [Geyer and Thomp- son, 1992, Hunter and Handcock, 2012, Descombes p is an exponential family distribution [Barndorff- gibbs et al., 1997]. Nielsen, 1978] with natural parameters θ. An alternate way to define the model is in terms of the mean value parameters, which are defined as the expected values 1.1 Example: The Ising Model of the sufficient statistics µ(θ) = E (g(X)). There θ A canonical example of a Gibbs measure is the Ising is a one-to-one relationship between the natural and [1925] model of ferromagnetism, which in this case we mean value parameters, so a model may be expressed will define over an n-by-n toroidal lattice. Each x in either terms [Barndorff-Nielsen, 1978]. ij may either have magnetic spin up or down (xij 2 One of the most attractive features of pgibbs as a class {−1; 1g 8 i; j 2 f1; :::; ng), and is influenced by its of general modeling distributions is that it satisfies the four neighbors x(i−1)j, x(i+1)j, xi(j−1), and xi(j+1). maximum entropy principle, in that among all possi- Here we allow the subscripts to wrap, such that a sub- ble probability distributions (p0) for X with identical script of 0 refers to index n and a subscript of n + 1 mean values (i.e. Ep0 (g(X)) = µ(θ)), pgibbs has the refers to index 1. The Ising probability model has largest entropy [Jaynes, 1957]. Because entropy is a two sufficient statistics. The first models the gen- measure of the informativeness of the distribution, the eral tendency for up or down spins and is defined P Gibbs measure may be thought of as the least infor- as g1(x) = i;j xij. The second is the number of mative distribution consistent with the mean value pa- neighbors sharing the same spin and is defined as P rameters. g2(x) = i;j xij(x(i+1)j + xi(j+1)). Performing likelihood-based inference on Gibbs mea- While the Ising model was developed in the context sures is a non-trivial task because for even moderately of statistical physics, the tendency of like to be con- multivariate X, the sum in the normalizing constant nected to like is a more universal concept, and thus the Z is typically computationally intractable, and thus it model has found extensive use in domains far afield is impossible to evaluate the likelihood directly. How- from magnetism. To name a few examples, it is used ever, indirect sampling methods may be used by ex- in computer vision for image denoising [Cohen et al., pressing the likelihood and its derivatives in terms of 2015], in network modeling the classic two-star model expectations and covariances. The log likelihood of can be viewed as an Ising model [Palla et al., 2004, the distribution is Park and Newman, 2004], and in sociology it is use- ful in understanding group choice [Galam, 1997] and X `gibbs(θjx) = θkgk(x) − log(Z(θ)); urban segregation [Stauffer, 2008]. k and the first and second derivatives may then be found 2 Phases In Gibbs Measures to be @`gibbs The Gibbs distribution captures many of the charac- = gi(x) − µi(θ) @θi teristics of typical physical systems, even those which in other domains may lead to computational difficul- and 2 ties and implausible models. In particular, Gibbs mea- @ `gibbs = −cov(gi(X); gj(X)): sures often exhibit the phenomena of a phase transi- @θ @θ i j tion point. This has been extensively studied in the physics literature [Georgii, 2011]. In the context of It is clear from setting the first derivative to zero that social network modeling, multiphase behavior has re- at the maximum likelihood estimate (MLE) of θ, the sulted in significant difficulty in finding good models mean value parameters must match the observed suf- for large complex networks [Handcock, 2003, Schwein- ficient statistics (i.e. g(x) = µ(θ^ )). mle berger, 2011, Chatterjee et al., 2013, Horv´atet al., Since pgibbs is known up to a constant, we may use 2015]. Markov Chain Monte Carlo (MCMC) methods to gen- Consider an Ising model where θ = 0 and θ is very erate a sample from it. Using this sample, the above 1 2 large. If one were to start an MCMC chain at a con- equations may be approximated using weighted sam- figuration where most vertices are set to 1. At equi- ple means and covariances in place of the population librium we would find almost perfect alignment to 1. expectations and covariances. The precise mechanics Similarly, if the starting point for the chain was mostly of using these approximations to find the maximum -1 then the equilibrium state would be a global -1 align- likelihood estimate of θ is known as Markov Chain ment. Monte Carlo Maximum Likelihood Estimation (MCM- CMLE), and has been used successfully across many Such behavior has important implications for prac- Ian E. Fellows, Mark S. Handcock by Horv´atet al. [2015], is a mathematical convenience removing the effects of small local fluctuations and as- suring the existence of a Hessian. While the choice of N* could result in a probability distribution, this is not a hard requirement of its construction, and in many cases, N* could be made to be arbitrarily close to N while still maintaining 2nd order differentiability. Following Horv´atet al. [2015] we define uniphase and multiphase as: Definition 1 A distribution p(g(x)) is uniphase with respect to p∗ if it contains no local minima or saddle points, and is multiphase otherwise. Horv´atet al. [2015] developed an insightful connec- tion between multiphase behavior and the configura- tion count function N.
Recommended publications
  • Ruelle Operator for Continuous Potentials and DLR-Gibbs Measures
    Ruelle Operator for Continuous Potentials and DLR-Gibbs Measures Leandro Cioletti Artur O. Lopes Departamento de Matem´atica - UnB Departamento de Matem´atica - UFRGS 70910-900, Bras´ılia, Brazil 91509-900, Porto Alegre, Brazil [email protected] [email protected] Manuel Stadlbauer Departamento de Matem´atica - UFRJ 21941-909, Rio de Janeiro, Brazil [email protected] October 25, 2019 Abstract In this work we study the Ruelle Operator associated to a continuous potential defined on a countable product of a compact metric space. We prove a generaliza- tion of Bowen’s criterion for the uniqueness of the eigenmeasures. One of the main results of the article is to show that a probability is DLR-Gibbs (associated to a continuous translation invariant specification), if and only if, is an eigenprobability for the transpose of the Ruelle operator. Bounded extensions of the Ruelle operator to the Lebesgue space of integrable functions, with respect to the eigenmeasures, are studied and the problem of exis- tence of maximal positive eigenfunctions for them is considered. One of our main results in this direction is the existence of such positive eigenfunctions for Bowen’s potential in the setting of a compact and metric alphabet. We also present a version of Dobrushin’s Theorem in the setting of Thermodynamic Formalism. Keywords: Thermodynamic Formalism, Ruelle operator, continuous potentials, Eigenfunc- tions, Equilibrium states, DLR-Gibbs Measures, uncountable alphabet. arXiv:1608.03881v6 [math.DS] 24 Oct 2019 MSC2010: 37D35, 28Dxx, 37C30. 1 Introduction The classical Ruelle operator needs no introduction and nowadays is a key concept of Thermodynamic Formalism.
    [Show full text]
  • The Dobrushin Comparison Theorem Is a Powerful Tool to Bound The
    COMPARISON THEOREMS FOR GIBBS MEASURES∗ BY PATRICK REBESCHINI AND RAMON VAN HANDEL Princeton University The Dobrushin comparison theorem is a powerful tool to bound the dif- ference between the marginals of high-dimensional probability distributions in terms of their local specifications. Originally introduced to prove unique- ness and decay of correlations of Gibbs measures, it has been widely used in statistical mechanics as well as in the analysis of algorithms on random fields and interacting Markov chains. However, the classical comparison theorem requires validity of the Dobrushin uniqueness criterion, essentially restricting its applicability in most models to a small subset of the natural parameter space. In this paper we develop generalized Dobrushin comparison theorems in terms of influences between blocks of sites, in the spirit of Dobrushin- Shlosman and Weitz, that substantially extend the range of applicability of the classical comparison theorem. Our proofs are based on the analysis of an associated family of Markov chains. We develop in detail an application of our main results to the analysis of sequential Monte Carlo algorithms for filtering in high dimension. CONTENTS 1 Introduction . .2 2 Main Results . .4 2.1 Setting and notation . .4 2.2 Main result . .6 2.3 The classical comparison theorem . .8 2.4 Alternative assumptions . .9 2.5 A one-sided comparison theorem . 11 3 Proofs . 12 3.1 General comparison principle . 13 3.2 Gibbs samplers . 14 3.3 Proof of Theorem 2.4................................. 18 3.4 Proof of Corollary 2.8................................. 23 3.5 Proof of Theorem 2.12................................ 25 4 Application: Particle Filters .
    [Show full text]
  • Continuum Percolation for Gibbs Point Processes Dominated by a Poisson Process, and Then We Use the Percolation Results Available for Poisson Processes
    Electron. Commun. Probab. 18 (2013), no. 67, 1–10. ELECTRONIC DOI: 10.1214/ECP.v18-2837 COMMUNICATIONS ISSN: 1083-589X in PROBABILITY Continuum percolation for Gibbs point processes∗ Kaspar Stuckiy Abstract We consider percolation properties of the Boolean model generated by a Gibbs point process and balls with deterministic radius. We show that for a large class of Gibbs point processes there exists a critical activity, such that percolation occurs a.s. above criticality. Keywords: Gibbs point process; Percolation; Boolean model; Conditional intensity. AMS MSC 2010: 60G55; 60K35. Submitted to ECP on May 29, 2013, final version accepted on August 1, 2013. Supersedes arXiv:1305.0492v1. 1 Introduction Let Ξ be a point process in Rd, d ≥ 2, and fix a R > 0. Consider the random S B B set ZR(Ξ) = x2Ξ R(x), where R(x) denotes the open ball with radius R around x. Each connected component of ZR(Ξ) is called a cluster. We say that Ξ percolates (or R-percolates) if ZR(Ξ) contains with positive probability an infinite cluster. In the terminology of [9] this is a Boolean percolation model driven by Ξ and deterministic radius distribution R. It is well-known that for Poisson processes there exists a critical intensity βc such that a Poisson processes with intensity β > βc percolates a.s. and if β < βc there is a.s. no percolation, see e.g. [17], or [14] for a more general Poisson percolation model. For Gibbs point processes the situation is less clear. In [13] it is shown that for some | downloaded: 28.9.2021 two-dimensional pairwise interacting models with an attractive tail, percolation occurs if the activity parameter is large enough, see also [1] for a similar result concerning the Strauss hard core process in two dimensions.
    [Show full text]
  • Mean-Field Monomer-Dimer Models. a Review. 3
    Mean-Field Monomer-Dimer models. A review. Diego Alberici, Pierluigi Contucci, Emanuele Mingione University of Bologna - Department of Mathematics Piazza di Porta San Donato 5, Bologna (Italy) E-mail: [email protected], [email protected], [email protected] To Chuck Newman, on his 70th birthday Abstract: A collection of rigorous results for a class of mean-field monomer- dimer models is presented. It includes a Gaussian representation for the partition function that is shown to considerably simplify the proofs. The solutions of the quenched diluted case and the random monomer case are explained. The presence of the attractive component of the Van der Waals potential is considered and the coexistence phase coexistence transition analysed. In particular the breakdown of the central limit theorem is illustrated at the critical point where a non Gaussian, quartic exponential distribution is found for the number of monomers centered and rescaled with the volume to the power 3/4. 1. Introduction The monomer-dimer models, an instance in the wide set of interacting particle systems, have a relevant role in equilibrium statistical mechanics. They were introduced to describe, in a simplified yet effective way, the process of absorption of monoatomic or diatomic molecules in condensed-matter physics [15,16,39] or the behaviour of liquid solutions composed by molecules of different sizes [25]. Exact solutions in specific cases (e.g. the perfect matching problem) have been obtained on planar lattices [24,26,34,36,41] and the problem on regular lattices is arXiv:1612.09181v1 [math-ph] 29 Dec 2016 also interesting for the liquid crystals modelling [1,20,27,31,37].
    [Show full text]
  • Markov Random Fields and Gibbs Measures
    Markov Random Fields and Gibbs Measures Oskar Sandberg December 15, 2004 1 Introduction A Markov random field is a name given to a natural generalization of the well known concept of a Markov chain. It arrises by looking at the chain itself as a very simple graph, and ignoring the directionality implied by “time”. A Markov chain can then be seen as a chain graph of stochastic variables, where each variable has the property that it is independent of all the others (the future and past) given its two neighbors. With this view of a Markov chain in mind, a Markov random field is the same thing, only that rather than a chain graph, we allow any graph structure to define the relationship between the variables. So we define a set of stochastic variables, such that each is independent of all the others given its neighbors in a graph. Markov random fields can be defined both for discrete and more complicated valued random variables. They can also be defined for continuous index sets, in which case more complicated neighboring relationships take the place of the graph. In what follows, however, we will look only at the most appro- achable cases with discrete index sets, and random variables with finite state spaces (for the most part, the state space will simply be {0, 1}). For a more general treatment see for instance [Rozanov]. 2 Definitions 2.1 Markov Random Fields Let X1, ...Xn be random variables taking values in some finite set S, and let G = (N, E) be a finite graph such that N = {1, ..., n}, whose elements will 1 sometime be called sites.
    [Show full text]
  • A Note on Gibbs and Markov Random Fields with Constraints and Their Moments
    NISSUNA UMANA INVESTIGAZIONE SI PUO DIMANDARE VERA SCIENZIA S’ESSA NON PASSA PER LE MATEMATICHE DIMOSTRAZIONI LEONARDO DA VINCI vol. 4 no. 3-4 2016 Mathematics and Mechanics of Complex Systems ALBERTO GANDOLFI AND PIETRO LENARDA A NOTE ON GIBBS AND MARKOV RANDOM FIELDS WITH CONSTRAINTS AND THEIR MOMENTS msp MATHEMATICS AND MECHANICS OF COMPLEX SYSTEMS Vol. 4, No. 3-4, 2016 dx.doi.org/10.2140/memocs.2016.4.407 ∩ MM A NOTE ON GIBBS AND MARKOV RANDOM FIELDS WITH CONSTRAINTS AND THEIR MOMENTS ALBERTO GANDOLFI AND PIETRO LENARDA This paper focuses on the relation between Gibbs and Markov random fields, one instance of the close relation between abstract and applied mathematics so often stressed by Lucio Russo in his scientific work. We start by proving a more explicit version, based on spin products, of the Hammersley–Clifford theorem, a classic result which identifies Gibbs and Markov fields under finite energy. Then we argue that the celebrated counterexample of Moussouris, intended to show that there is no complete coincidence between Markov and Gibbs random fields in the presence of hard-core constraints, is not really such. In fact, the notion of a constrained Gibbs random field used in the example and in the subsequent literature makes the unnatural assumption that the constraints are infinite energy Gibbs interactions on the same graph. Here we consider the more natural extended version of the equivalence problem, in which constraints are more generally based on a possibly larger graph, and solve it. The bearing of the more natural approach is shown by considering identifi- ability of discrete random fields from support, conditional independencies and corresponding moments.
    [Show full text]
  • Large Deviations and Applications 29. 11. Bis 5. 12. 1992
    MATHEMATISCHES FORSCHUNGSINSTITUT OBERWOLFACH Tagungsbericht 51/1992 Large Deviations and Applications 29. 11. bis 5. 12. 1992 Zu dieser Tagung unter der gemeinsamen Leitung von E. Bolthausen (Zurich), J. Gartner (Berlin) und S. R. S. Varadhan (New York) trafen sich Mathematiker und mathematische Physiker aus den verschieden- sten Landern mit einem breiten Spektrum von Interessen. Die Theorie vom asymptotischen Verhalten der Wahrscheinlichkeiten groer Abweichungen ist einer der Schwerpunkte der jungeren wahr- scheinlichkeitstheoretischen Forschung. Es handelt sich um eine Pra- zisierung von Gesetzen groer Zahlen. Gegenstand der Untersuchun- gen sind sowohl die Skala als auch die Rate des exponentiellen Abfalls der kleinen Wahrscheinlichkeiten, mit denen ein untypisches Verhal- ten eines stochastischen Prozesses auftritt. Die Untersuchung dieser kleinen Wahrscheinlichkeiten ist fur viele Fragestellungen interessant. Wichtige Themen der Tagung waren: – Anwendungen groer Abweichungen in der Statistik – Verschiedene Zugange zur Theorie groer Abweichungen – Stochastische Prozesse in zufalligen Medien – Wechselwirkende Teilchensysteme und ihre Dynamik – Statistische Mechanik, Thermodynamik – Hydrodynamischer Grenzubergang – Verhalten von Grenzachen, monomolekulare Schichten – Dynamische Systeme und zufallige Storungen – Langreichweitige Wechselwirkung, Polymere – Stochastische Netzwerke Die Tagung hatte 47 Teilnehmer, es wurden 42 Vortrage gehalten. 1 Abstracts G. Ben Arous (joint work with A. Guionnet) Langevin dynamics for spin glasses We study the dynamics for the Sherrington-Kirkpatrick model of spin glasses. More precisely, we take a soft spin model suggested by Som- polinski and Zippelius. Thus we consider diusions interacting with a Gaussian random potential of S-K-type. (1) We give an “annealed” large deviation principle for the empirical measure at the process level. We deduce from this an annealed law of large numbers and a central limit theorem.
    [Show full text]
  • Stein's Method for Discrete Gibbs Measures
    The Annals of Applied Probability 2008, Vol. 18, No. 4, 1588–1618 DOI: 10.1214/07-AAP0498 c Institute of Mathematical Statistics, 2008 STEIN’S METHOD FOR DISCRETE GIBBS MEASURES1 By Peter Eichelsbacher and Gesine Reinert Ruhr-Universit¨at Bochum and University of Oxford Stein’s method provides a way of bounding the distance of a prob- ability distribution to a target distribution µ. Here we develop Stein’s method for the class of discrete Gibbs measures with a density eV , where V is the energy function. Using size bias couplings, we treat an example of Gibbs convergence for strongly correlated random vari- ables due to Chayes and Klein [Helv. Phys. Acta 67 (1994) 30–42]. We obtain estimates of the approximation to a grand-canonical Gibbs ensemble. As side results, we slightly improve on the Barbour, Holst and Janson [Poisson Approximation (1992)] bounds for Poisson ap- proximation to the sum of independent indicators, and in the case of the geometric distribution we derive better nonuniform Stein bounds than Brown and Xia [Ann. Probab. 29 (2001) 1373–1403]. 0. Introduction. Stein [17] introduced an elegant method for proving convergence of random variables toward a standard normal variable. Bar- bour [2, 3] and G¨otze [10] developed a dynamical point of view of Stein’s method using time-reversible Markov processes. If µ is the stationary dis- tribution of a homogeneous Markov process with generator , then X µ if and only if E g(X) = 0 for all functions g in the domainA ( ) of∼ the operator . ForA any random variable W and for any suitable functionD A f, to assess theA distance Ef(W ) f dµ we first find a solution g of the equation | − | g(Rx)= f(x) f dµ.
    [Show full text]
  • Boundary Conditions for Translation- Invariant Gibbs Measures of the Potts Model on Cayley Trees Daniel Gandolfo, M
    Boundary Conditions for Translation- Invariant Gibbs Measures of the Potts Model on Cayley Trees Daniel Gandolfo, M. M. Rakhmatullaev, U.A. Rozikov To cite this version: Daniel Gandolfo, M. M. Rakhmatullaev, U.A. Rozikov. Boundary Conditions for Translation- In- variant Gibbs Measures of the Potts Model on Cayley Trees. Journal of Statistical Physics, Springer Verlag, 2017, 167 (5), pp.1164. 10.1007/s10955-017-1771-5. hal-01591298 HAL Id: hal-01591298 https://hal.archives-ouvertes.fr/hal-01591298 Submitted on 26 Apr 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Boundary conditions for translation-invariant Gibbs measures of the Potts model on Cayley trees D. Gandolfo Aix Marseille Univ, Universit´ede Toulon, CNRS, CPT, Marseille, France. [email protected] M. M. Rahmatullaev and U. A. Rozikov Institute of mathematics, 29, Do'rmon Yo'li str., 100125, Tashkent, Uzbekistan. [email protected] [email protected] Abstract We consider translation-invariant splitting Gibbs measures (TISGMs) for the q- state Potts model on a Cayley tree of order two. Recently a full description of the TISGMs was obtained, and it was shown in particular that at sufficiently low temper- atures their number is 2q − 1.
    [Show full text]
  • Arxiv:2103.07391V1 [Math.PR] 12 Mar 2021 Studied
    GIBBS MEASURES OF POTTS MODEL ON CAYLEY TREES: A SURVEY AND APPLICATIONS U.A. ROZIKOV Abstract. In this paper we give a systematic review of the theory of Gibbs mea- sures of Potts model on Cayley trees (developed since 2013) and discuss many applications of the Potts model to real world situations: mainly biology, physics, and some examples of alloy behavior, cell sorting, financial engineering, flocking birds, flowing foams, image segmentation, medicine, sociology etc. 1. Introduction The Potts model is defined by a Hamiltonian (energy) of configurations of spins which take one of q possible values on vertices of a lattice. The model is used to study the behavior of systems having multiple states (spins, colors, alleles, etc). Since the model has a rich mathematical formulation it has been studied extensively. Usually the results of investigations of a system (in particular, the Potts model) are presented by a measure assigning a number to each (interesting) suitable properties of the system. The Gibbs measure is one of such important measures in many problems of probability theory and statistical mechanics. It is the measure associated with the Hamiltonian of a (biological, physical) system. Each Gibbs measure gives a state of the system. The main problem for a given Hamiltonian on a countable lattice is to describe its all possible Gibbs measures. In case of uniqueness of the Gibbs measure (for all values of the parameters, in particular a parameter can be temperature), the system does not change its state. The existence of some values of parameters at which the uniqueness of Gibbs measure switches to non-uniqueness is interpreted as phase transition (the system will change its state).
    [Show full text]
  • Lecture Notes Gibbs Measures and Phase Transitions. Part 2
    Lecture notes Gibbs measures and phase transitions. Part 2 Anton Bovier Rheinische Friedrich-Wilhelms Universit¨at Bonn Institut f¨ur Angewandte Mathematik Endenicher Allee 60 53115 Bonn Contents 1 Gibbsian formalism and metastates page 1 1.1 Introduction 1 1.2 Random Gibbs measures and metastates 3 1.3 Remarks on uniqueness conditions 13 1.4 Phase transitions 13 1.5 The Edwards–Anderson model. 17 2 The random-field Ising model 19 2.1 The Imry–Ma argument 19 2.2 The Aizenman–Wehr method 28 2.2.1 Translation covariant states 28 2.2.2 Order parameters and generating functions 30 2.3 The Bricmont-Kupiainen renormalization group 36 2.3.1 Renormalization group and contour models 38 2.3.2 The ground-states 44 2.3.3 The Gibbs states at finite temperature 58 Bibliography 77 i 1 Gibbsian formalism and metastates Longtemps les objets dont s’occupent les math´ematiciens ´etaient pour la pluspart mal d´efinis; on croyait les connaˆıtre, parce qu’on se les repr´esentait avec le sens ou l’imagination; mais on n’en avait qu’une image grossi`ere et non une id´ee pr´ecise sur laquelle le raisonnement pˆut avoir prise1. Henri Poincar´e, La valeur de la science We now turn to the main topic of this book, disordered systems. We split this into two parts, treating in turn lattice models and mean-field models. From the physical point of view, the former should be more realistic and hence more relevant, so it is natural that we present the general mathemat- ical framework in this context.
    [Show full text]
  • Distribution Dependent Stochastic Differential Equations
    Distribution Dependent Stochastic Differential Equations∗ Xing Huanga), Panpan Renb), Feng-Yu Wanga) a) Center for Applied Mathematics, Tianjin University, Tianjin 300072, China b) Department of Mathematics, City University of HongKong, HongKong, China [email protected], [email protected], [email protected] December 29, 2020 Abstract Due to their intrinsic link with nonlinear Fokker-Planck equations and many other applications, distribution dependent stochastic differential equations (DDSDEs for short) have been intensively investigated. In this paper we summarize some recent progresses in the study of DDSDEs, which include the correspondence of weak solutions and nonlinear Fokker-Planck equations, the well-posedness, regularity estimates, exponential ergodicity, long time large deviations, and comparison theorems. AMS subject Classification: 60B05, 60B10. Keywords: DDSDE, nonlinear Fokker-Planck equation, Bismut formula, Wasserstein distance, gradient estimate. arXiv:2012.13656v1 [math.PR] 26 Dec 2020 1 Introduction To characterize nonlinear PDEs in Vlasov’s kinetic theory, Kac [27, 28] proposed the “propagation of chaos” of mean field particle systems, which stimulated McKean [33] to study nonlinear Fokker-Planck equations using stochastic differential equations with distribution dependent drifts, see [45] for a theory on mean field particle systems and applications. In general, a nonlinear Fokker-Planck equation can be characterized by the following distri- bution dependent stochastic differential equations (DDSDEs for short): L L (1.1) dXt = b(t, Xt, Xt )dt + σ(t, Xt, Xt )dWt, ∗Supported in part by NNSFC (11771326, 11831014, 11921001, 11801406). 1 where Wt is an m-dimensional Brownian motion on a complete filtration probability space (Ω, F , P), L is the distribution (i.e.
    [Show full text]