This thesis presents and applies the M-layer construction, a recent diagrammatic method for studying the critical behavior of lattice models in finite dimensions. The construction provides a unified tool to compute correlation functions and critical exponents through a perturbative expansion in topological loops, whose leading order corresponds to the behavior on a Bethe lattice. The first application of the method is on well-known problems, such as the Ising model and the percolation problem, where it successfully reproduces standard results, thus validating its efficiency. It is then applied to the Edwards–Anderson spin glass model in a magnetic field at zero temperature, a problem whose finite-dimensional physics remains widely debated. The remarkable new result of the M-layer construction is the identification, for the first time through analytical computations, of a phase transition in a field below the upper critical dimension. The results, obtained through analytical computations and partially checked numerically, show that this approach provides a promising alternative to standard Renormalization Group techniques for investigating phase transitions in disordered and frustrated systems.
This thesis presents and applies the M-layer construction, a recent diagrammatic method for studying the critical behavior of lattice models in finite dimensions. The construction provides a unified tool to compute correlation functions and critical exponents through a perturbative expansion in topological loops, whose leading order corresponds to the behavior on a Bethe lattice. The first application of the method is on well-known problems, such as the Ising model and the percolation problem, where it successfully reproduces standard results, thus validating its efficiency. It is then applied to the Edwards–Anderson spin glass model in a magnetic field at zero temperature, a problem whose finite-dimensional physics remains widely debated. The remarkable new result of the M-layer construction is the identification, for the first time through analytical computations, of a phase transition in a field below the upper critical dimension. The results, obtained through analytical computations and partially checked numerically, show that this approach provides a promising alternative to standard Renormalization Group techniques for investigating phase transitions in disordered and frustrated systems.
Abstract In statistical physics, one of the standard methods to study second order phase transitions is the renormalization group that usually leads to an expansion around the corresponding fully connected solution. Unfortunately, often in disordered models, some important finite dimensional second-order phase transitions are qualitatively different or absent in the corresponding fully connected model: in such cases the standard expansion fails. Recently, a new method, the M-layer one, has been introduced that performs an expansion around a different soluble mean field model: the Bethe lattice one. This new method has been already used to compute the upper critical dimension D U of different disordered systems such as the Random Field Ising model or the Spin glass model with field. If then one wants to go beyond and construct an expansion around D U to understand how critical quantities get renormalized, the actual computation of all the numerical factors is needed. This next step has still not been performed, being technically more involved. In this paper we perform this computation for the ferromagnetic Ising model without quenched disorder, in finite dimensions: we show that, at one-loop order inside the M-layer approach, we recover the continuum quartic field theory and we are able to identify the coupling constant g and the other parameters of the theory, as a function of macroscopic and microscopic details of the model such as the lattice spacing, the physical lattice dimension and the temperature. This is a fundamental step that will help in applying in the future the same techniques to more complicated systems, for which the standard field theoretical approach is impracticable.
Abstract For every physical model defined on a generic graph or factor graph, the Bethe M -layer construction allows building a different model for which the Bethe approximation is exact in the large M limit, and coincides with the original model for M = 1 . The 1 / M perturbative series is then expressed by a diagrammatic loop expansion in terms of so-called fat diagrams. Our motivation is to study some important second-order phase transitions that do exist on the Bethe lattice, but are either qualitatively different or absent in the corresponding fully connected case. In this case, the standard approach based on a perturbative expansion around the naive mean field theory (essentially a fully connected model) fails. On physical grounds, we expect that when the construction is applied to a lattice in finite dimension there is a small region of the external parameters, close to the Bethe critical point, where strong deviations from mean-field behavior will be observed. In this region, the 1 / M expansion for the corrections diverges, and can be the starting point for determining the correct non-mean-field critical exponents using renormalization group arguments. In the end, we will show that the critical series for the generic observable can be expressed as a sum of Feynman diagrams with the same numerical prefactors of field theories. However, the contribution of a given diagram is not evaluated by associating Gaussian propagators to its lines, as in field theories: one has to consider the graph as a portion of the original lattice, replacing the internal lines with appropriate one-dimensional chains, and attaching to the internal points the appropriate number of infinite-size Bethe trees to restore the correct local connectivity of the original model. The actual contribution of each (fat) diagram is the so-called line-connected observable, which also includes contributions from sub-diagrams with appropriate prefactors. In order to compute the corrections near to the critical point, Feynman diagrams (with their symmetry factors) can be read directly from the appropriate field-theoretical literature; the computation of momentum integrals is also quite similar; the extra work consists of computing the line-connected observable of the associated fat diagram in the limit of all lines becoming infinitely long.
In a recent work [M. Castellana and G. Parisi, Phys. Rev. E 82, 040105(R) (2010)], the large-scale behavior of the simplest non-mean-field spin-glass system has been analyzed, and the critical exponent related to the divergence of the correlation length has been computed at two loops within the ε-expansion technique by two independent methods. By performing the explicit calculation of the critical exponents at two loops, one obtains that the two methods yield the same result. This shows that the underlying renormalization group ideas apply consistently in this disordered model, in such a way that an ε-expansion can be set up. The question of the extension to high orders of this ε-expansion is particularly interesting from the physical point of view. Indeed, once high orders of the series in ε for the critical exponents are known, one could check the convergence properties of the series, and find out if the ordinary series resummation techniques, yielding very accurate predictions for the Ising model, work also for this model. If this is the case, a consistent and predictive non-mean-field theory for such a disordered system could be established. In that regard, in this work we expose the underlying techniques of such a two-loop computation. We show with an explicit example that such a computation could be quite easily automatized, i.e., performed by a computer program, in order to compute high orders of the ε-expansion, and so eventually make this theory physically predictive. Moreover, all the underlying renormalization group ideas implemented in such a computation are widely discussed and exposed.
"Renorma1ization-group methods have been applied in the study of quenched random magnetic systems in recent years. We begin with a brief review of second-order phase transitions in pure, homogeneous systems and also of the .. renorma1ization group framework. Then we provide an introduction to quenched random magnetic systems. Next, momentum-space methods and position-space techniques as applied to quenched random magnets are outlined and compared. Grinstein and Luther applied the Wilson-Fisher E-expansion to random n-vector models; Khme1'nitsky discovered that the random Ising model (n = 1) possessed a ""random"" fixed point of 0(squareroot(e)). This fixed point was found to have one marginal and one irrelevant operator. We have investigated the stability of this fixed point using Ca11an-Symanzik equations and renorma1ized perturbation theory. We find the fixed point stable in the next order; we have also obtained critical exponents to one higher order. Next, position-space techniques are used to study some simple model systems. In addition to critical exponents, global thermodynamic properties are determined. These calculations are based on the Migda1-Kadanoff approximate recursion relations suitably generalized to the inhomogeneous case. Firstly we study the randomly bond-dilute two-dimensional nearest- neighbor Ising model on a square lattice. Calculations give both thermal and magnetic exponents associated with the percolative fixed point. Differential recursion relations yield a phase diagram which is in quantitative agreement with all known results. Curves for the specific heat, percolation probability, and magnetization are displayed. The critical region of the specific heat becomes unobservably narrow well above the percolation threshold Pc. This provides a possible explanation for the apparent specific-heat rounding in certain experiments. We then study the Edwards-Anderson model of a spin glass. The current theoretical situation, which is far from satisfactory at present, is briefly reviewed. We treat the spin-l/2 Ising model with independently random nearestneighbor interactions in dimensionalities d = 2, 3, and 4. The phase diagram, which is in qualitative agreement with mean-field results, exhibits paramagnetic, ferromagnetic, antiferromagnetic, and spin-glass phases. The spinglass and paramagnetic phases meet along an extended second-order phase boundary, which terminates in two tricritical points. Critical and tricritical exponents are calculated. The spin-glass specific-heat exponent turns out to be large and negative, compatibly with recent experiments which show a rounded specific heat anomaly. Global specific-heat curves are also displayed for d = 2."
One way to perform field theory computations for the bond percolation problem is through the Kasteleyn and Fortuin mapping to the n+1 n + 1 states Potts model in the limit of n \to 0 n → 0 . In this paper, we show that it is possible to recover the \epsilon ϵ -expansion for critical exponents in finite dimension directly using the M M -layer expansion, without the need to perform any analytical continuation. Moreover, we also show explicitly that the critical exponents for site and bond percolation are the same. This computation provides a reference for applications of the M M -layer method to systems where the underlying field theory is unknown or disputed.
The analysis of the error of perturbative series by comparing it to the exact solution is an important tool to understand the non-perturbative physics of statistical models. For some toy models, a new method can be used to calculate higher order weak coupling expansion and modified perturbation theory can be constructed. However, it is nontrivial to generalize the new method to understand the critical behavior of high dimensional spin and gauge models. Actually, it is a big challenge in both high energy physics and condensed matter physics to develop accurate and efficient numerical algorithms to solve these problems. In this thesis, one systematic way named tensor renormalization group method is discussed. The applications of the method to several spin and gauge models on a lattice are investigated. theoretically, the new method allows one to write an exact representation of the partition function of models with local interactions. E.g. O(N) models, Z2 gauge models and U(1) gauge models. Practically, by using controllable approximations, results in both finite volume and the thermodynamic limit can be obtained. Another advantage of the new method is that it is insensitive to sign problems for models with complex coupling and chemical potential. Through the new approach, the Fisher's zeros of the 2D O(2) model in the complex coupling plane can be calculated and the finite size scaling of the results agrees well with the Kosterlitz-Thouless assumption. Applying the method to the O(2) model with a chemical potential, new phase diagram of the models can be obtained. The structure of the tensor language may provide a new tool to understand phase transition properties in general.