Charting the scaling region of the Ising universality class in two and three dimensions

We study the behaviour of a universal combination of susceptibility and correlation length in the Ising model in two and three dimensions, in presence of both magnetic and thermal perturbations, in the neighbourhood of the critical point. In three dimensions we address the problem using a parametric representation of the equation of state. In two dimensions we make use of the exact integrability of the model along the thermal and the magnetic axes. Our results can be used as a sort of “reference frame” to chart the critical region of the model. While our results can be applied in principle to any possible realization of the Ising universality class, we address in particular, as speciﬁc examples, three instances of Ising behaviour in ﬁnite temperature QCD related in various ways to the deconﬁnement transition. In particular, in the last of these examples, we study the critical ending point in the ﬁnite density, ﬁnite temperature phase diagram of QCD. In this ﬁnite density framework, due to well know sign problem, Montecarlo simulations are not possible and thus a direct comparison of experimental results with QFT/Statmech predictions like the one we discuss in this paper may be important. Moreover in this example it is particularly diﬃcult to disentangle “magnetic-like” from “thermal-like” observables and thus an explicit charting of the neighbourhood of the critical point can be particularly useful.


Introduction
Despite its apparent simplicity the Ising model is one of the cornerstones of modern Statistical Mechanics. Over the years it has become a theoretical laboratory to test new ideas, ranging from Symmetry Breaking to Conformal Field Theories. Moreover, thanks to its exact solvability in two dimensions [1,2] and the ease with which it can be simulated in three dimensions, it has been widely used as a benchmark to test new numerical approaches and innovative approximations.
The corresponding universality class, in the Renormalization Group sense [3], is of central importance in theoretical physics due to its many experimental realizations in different physical contexts, ranging from condensed matter to high energy physics. At the same time it describes the critical behaviour of a lot of different spin models and, in the limit of high temperatures, also of gauge theories.
From a Statistical Mechanics point of view it represents the simplest way to describe systems with short range interactions and a scalar order parameter (density or uniaxial magnetization) which undergo a symmetry breaking phase transition. From a Quantum Field Theory (QFT) point of view it is the simplest example of a unitary Conformal Field Theory (CFT) [4] perturbed by only two relevant operators: the "spin" operator (which is Z 2 odd) and the "energy" operator (Z 2 even) [5].
Thanks to integrability, conformal perturbation and bootstrap [6,7] lots of results are known, both in two and in three dimensions, on the behaviour of the model at the critical point, or when only one of the two perturbing operators is present. However, typically, the interesting regime for most of the experimental realizations of the model is when both the perturbing operators are present and much less is known in this situation.
The aim of this paper is to partially fill this gap by studying a suitable universal combination of thermodynamic quantities (see below for the precise definition) in presence of both perturbing operators. In three dimensions we shall address the problem using a parametric representation of the equation of state [8], while in two dimensions we shall make use of the exact integrability of the model in presence of a single perturbation [5]. Using these tools we shall be able to predict the value of this quantity in the whole phase space of the model in the neighbourhood of the critical point. These values can be used as a sort of "reference frame" to chart the critical region of the model. The universal combination that we shall study involves the magnetic susceptibility and thus our proposal is particularly effective when the model is characterized by an explicit Z 2 symmetry. When this is not the case, like for the liquid-vapour transition or for the finite density QCD example that we shall discuss below, the explicit knowledge of our universal combination may help to identify the exact directions in the phase space of the model with respect to which the magnetic susceptibility must be evaluated.
Thanks to universality, our results hold not only for the standard nearest neighbour Ising model, but also for any possible realization of the Ising universality class and in fact we shall use the high precision Montecarlo estimates obtained from an improved version of the Ising model to benchmark and test our results [9,10,11,12,13,14,15].
In particular, we shall concentrate in the second part of the paper on realizations in the context of high energy physics, suggested by the lattice regularization of QCD. We shall discuss three instances of Ising behaviour in finite temperature QCD related in various ways to the deconfinement transition. In the last of these examples, we shall address the critical ending point of finite density QCD. In this case, due to well know sign problem, Montecarlo simulations are not possible and thus a direct comparison of experimental results with QFT/Statmech predictions like the one we discuss in this paper may be important. This paper is organized as follows. Sect. 2 is devoted to a general introduction to the model and to the universal combination of thermodynamic quantities which is the main subject of the paper. In sect. 3 we shall address the problem in three dimensions using a suitable parametric representation of the equation of state of the model. We shall also show that the same approach cannot be used in two dimensions. In sect. 4 we shall then address the two dimensional case using appropriate expansions around the exact solutions of the model. Finally sect. 5 will be devoted to the discussion of a set of examples in high temperature QCD. We collected in the appendices some additional material which may be useful to reproduce our numerical analysis.

General information on the Ising universality class
The Ising model has a global Z 2 symmetry and is characterized by two relevant operators which encode the Z 2 odd (σ) and Z 2 even ( ) perturbations of the critical point. From a QFT point of view, the model in the vicinity of the critical point can be written as a perturbed Conformal Field Theory where (x) and σ(x) are the energy and spin operators and represent the continuum limits of the lattice operators ij σ i σ j and i σ i respectively. These operators are conjugated to the reduced temperature t = 1 Tc (T − T c ) and magnetic field H, which measure the deviation from the critical point. The action S CF T is the conformal-invariant action of the model at the critical point. In two dimensions this is the action of a free massless Majorana fermion with central charge c = 1/2.
Thanks to the exact integrability of the model for H = 0 (pure thermal perturbation) and for t = 0 (pure magnetic perturbation) much is known of this QFT in two dimensions. In particular all the critical exponents and the universal ratios are know exactly and, as we shall see below, reliable expansions around the integrable lines can be constructed for several observables. In three dimensions there are not exact results, but from the recent progress of the bootstrap approach and the improvement of Montecarlo methods several universal quantities can be evaluated with very high precision.
The most important realization of this QFT is the spin Ising model on a cubic (in d = 3) or square (in d = 2) lattice, which we shall use in the following to fix notations. As it is well known, the model is defined by the following energy function where the spins σ i can take the values σ i = ±1, the index i labels the sites of the lattice, the symbol ij means that the sum is performed over pairs of nearest neighbour sites, J is the coupling strength between spins (we assume a positive isotropic interaction so that for all pairs of nearest neighbour spins J ij = J > 0), andĤ is the external magnetic field. The partition function of the model is let us define β = J/k B T and H =Ĥ/k B T . For H = 0 the model is explicitly Z 2 symmetric and is characterized by two phases, a low temperature phase in which the Z 2 symmetry is spontaneously broken and a spontaneous magnetization is present and a high temperature phase in which the Z 2 symmetry is restored. The two phases are separated by a critical point. If one switches on the magnetic field it becomes apparent that the low T phase is actually a line of first order phase transitions which ends with the critical point. In the following we shall be interested in the scaling region in the vicinity of this critical point. Standard Renormalization Group arguments tell us that in this limit the irrelevant operators of the model can be neglected and the behaviour is completely described only by the two relevant operators , σ and one can perform a continuum limit of the model which leads exactly to the QFT described by eq. 1. From the partition function defined above it is easy to obtain all the thermodynamic observables. In particular, following the standard notation we have for the magnetization and the magnetic susceptibility The exponential correlation length can be extracted from the large distance decay of the spin-spin connected correlator as In several practical applications it is also useful the so called "secondmoment" correlation length which is defined through the second moment of the spin-spin correlation function as and it is simpler to evaluate than ξ both in numerical simulations and in experiments.
In the scaling limit the critical behaviour of all thermodynamic quantities is controlled by the two "scaling exponents" x and x σ which are universal and are shared by all physical realizations of the Ising universality class. The corresponding amplitudes are not universal, but one can construct suitable combinations in which the non-universal features of the model cancel out (see appendix A) and represent testable predictions of the Ising QFT to be compared with any possible realization of the Ising universality class.
While this is a well studied subject when only a single perturbation is present, its extension to the whole scaling region of the model, where both the H and t perturbations are present is not straightforward.
The main goal of this paper is to show that such an extension can be easily obtained making use of a parametric representation of the model and that the resulting universal quantities can be used as a natural "reference frame" to chart the scaling region of the Ising universality class.
While the parametric approach is completely general and could be applied in principle to any universal combination of thermodynamic quantities, in this paper we shall study in particular the following ratio and its natural extension to the second moment correlation length where Γ − , ξ − and ξ 2,− denote the amplitudes of χ, ξ and ξ 2 along the t < 0, H = 0 axis (see appendix A for detailed definitions and normalizations). This choice is motivated by the fact that the two observables which appear in the ratio are rather easy to evaluate, both in numerical simulations and in experiments, since they only involve derivatives or correlations of the order parameter and are normalized with respect to the values they have along the critical line of first order phase transitions, which is easy to identify (again, both numerically and experimentally).
The main drawback of this choice is that it assumes an explicit realization of the Z 2 symmetry. While in many interesting applications, like for the liquid-vapour transition in which the role of the perturbing parameter is played by the density, this is not the case and the Z 2 symmetry is just an "emergent" symmetry. The typical approach in these cases, following Rehr and Mermin [16], is to realize the t, H perturbations as suitable linear combinations of the actual variables of the model.
In the scaling region, when both the relevant perturbations are present, all the thermodynamic observables depend on the scaling combination 1 The three limits in which only one of the two perturbations is present (H = 0, t < 0), (H = 0, t = 0) and (H = 0, t > 0) correspond respectively to η = −∞, η = 0 and η = +∞. In these limits Ω can be written in terms of the standard universal amplitude ratios Q 2 , Γ + /Γ − and ξ − /ξ + (see appendix A) as follows These values can be used as benchmarks to test the reliability of our estimates and as "anchors" of the reference frame we are constructing.

Parametric Representation
It is useful to introduce a parametric representation of the critical equation of state, that not only satisfies the scaling hypothesis but additionally allows a simpler implementation of the analytic properties of the equation of state itself. Following [17], we express the thermodynamic variables t, H in terms of a couple of parameters R, θ both positive. Intuitively, the first measures the distance form the critical point in the (t, H) plane while the latter corresponds to the angular displacement along lines of constant R around the critical point. The parametrization of t and H results in a parametric expression of M as well, explicitly Calling θ 0 > 1 the smallest positive zero of the function h(θ), we see from the system (eq. 11) that the domain of interest in the (R, θ) plane is 0 ≤ θ ≤ θ 0 for every R ≥ 0 and that θ = θ 0 corresponds to the H = 0, t < 0 axis, θ = 1 to the H = 0, t = 0 axis and θ = 0 to the H = 0, t > 0 one.
The key point of the whole analysis is the determination of the function h(θ). There are some general properties which h(θ) must satisfy. It must be analytic in this physical domain in order to satisfy the regularity properties of the critical equation of state, i.e. the so-called Griffiths' analyticity [18].
Moreover, it must be an odd function of θ because of the Z 2 symmetry of the system. The most general choice is thus a polynomial of the type Using standard QFT methods [19,8,20,21] one can extract the coefficients h 2n+1 from a high temperature expansion of the free energy of the model. Using a variational method it is possible to obtain reliable and stable estimates of the coefficients up to h 7 in three dimensions and to h 11 in two dimensions. We can use the known amplitude ratios as benchmarks to evaluate the reliability of these parametric representations. As we will see h 7 will be enough to obtain estimates which in three dimensions agree within the errors with the amplitude ratios. The situation is worse in two dimensions, and this will prompt us to address the 2d case with a different approach.
A similar parametric representation can be introduced also for the correlation length. Following [22] we parametrize the square mass of the underlying QFT i.e. ξ −2 as follows and similarly, for the second moment correlation length The constants c and c 2 can be fixed using the universal ratios ξ + /ξ − and ξ 2,+ /ξ 2,− respectively and we use then the Q 2 and (Q 2 ) 2nd ratios to test the parametric representation.
Truncating at the quadratic order eq.s 13 and 14 could seem a too drastic approximation, but we will see below that in d = 3 it gives quite good results. The same is not true in d = 2 where however, as we anticipated above, we shall use a different approach to evaluate Ω(η).
Using the above results we may construct a parametric representation of Ω as a function of θ with and a similar expression for Ω 2 (θ) with c → c 2 .
There are a few universal features of Ω that we may deduce from this expression and hold for any possible realization of the Ising universality class • Ω is a monotonic decreasing function of θ (see fig. 1 for a plot of Ω(θ) in d = 3), therefore it can be inverted. From any given experimental estimate of Ω a precise value of θ may be extracted.
In practical applications one may be interested in the expression of Ω as a function of η = t/|H| 1 βδ . It is useful to introduce a "universal" version of the scaling variable defined asη Using the parametric representation of eq. 11 it is possible to write the expansion of θ as a function ofη in the neighbourhood of the three singular pointsη = (±∞, 0) and then Ω as a function ofη, which we plot in fig. 2 below in the d = 3 case (see appendix B for more details on this expansion).

Results in d = 3
In three dimensions we have [8,20] and the resulting value for the smallest positive root is θ 2 0 = 1.37861. Inserting these values in eq. 33 we find Γ + /Γ − ∼ 4.76 which is in good agreement (with a difference of the order of 1%) with the Montecarlo estimate Γ + /Γ − = 4.714(4) reported in [11,13].
Using the known values of ξ + /ξ − and ξ 2,+ /ξ 2,− we then obtain c = 0.0416 and c 2 = 0.0776. To check the reliability of this parametric representation we may use them to predict Q 2 and (Q 2 ) 2nd : we end up with Q 2 ∼ 1.250 and (Q 2 ) 2nd ∼ 1.209 which, again, are in agreement (with a difference of the order of 3%) with the best numerical estimates reported in eq.s 28, 29. This tells us that we can trust the parametric representation of Ω(θ) discussed above. We plot the result in fig. 1. We find in particular Ω iso = 1.08376... and Ω max = 1.35502... In fig. 2 we plot the expression of Ω as a function ofη.

Exact expression for Ω(η) in the d = 2 case
In d = 2, one can obtain much more precise results performing a perturbative expansion around the exact solutions along the two axes H = 0 and t = 0. A powerful tool to study the free energy of a perturbed CFT is the well known Truncated Conformal Space Approach (TCSA) [24,25]. In our case, thanks to the exact mapping of the H = 0 model to the QFT of a free Majorana fermion it is possible to construct a more effective version of the TCSA which uses the free fermions as a basis, the "Truncated Free-Fermion Space Approach" [5]. With this approach it is possible to evaluate the free energy for almost all values of η [5,26,27,28] and from that of the susceptibility. With similar methods it is also possible to study the perturbed mass spectrum of the theory [26,27,28] and hence the correlation length ξ which is the inverse of the lowest mass of the spectrum M 1 . The resulting expansions for χ and M 1 , in the three singular limits, are reported in appendix C.
Combining these quantities we obtain a precise estimate for Ω, which we plot in fig. 3. The resulting polynomial expansions in the three regions of interest are The coefficients are reported in tab. 1.  Table 1: Expansion coefficients of Ω(η) in the three regimes of interest, according to eq.s 20.
As expected, the three limiting cases Ω(±∞), Ω(0) (which correspond to the n = 0 values in the table) agree with the universal values obtained plugging in eq. 10 the universal ratios quoted in eq. 25.

Three examples of Ising behaviour in finite temperature QCD
Among the many physical realizations of the Ising universality class in this paper we decided to focus on three examples taken from high energy physics and in particular from the lattice regularization of QCD at finite temperature.

The deconfinement transition in the SU (2) pure gauge theory
The most direct realization of the Ising universality class in Lattice Gauge Theories (LGTs) is given by the deconfinement transition of pure gauge theories with a symmetry group G which has Z 2 as center. This result is based on the Svetitsky Yaffe approach [29] to the study of finite temperature (d+1)-dimensional pure gauge theories. The main idea of [29] is to construct a d-dimensional effective theory from the original one by integrating out the spacelike links and keeping as only remaining degrees of freedom the Polyakov loops. These loops are then treated as spins of an effective ddimensional model whose global symmetry must coincide with the center of the original gauge group. If both the phase transitions of the original gauge theory and that of the effective spin model are continuous, they must belong to the same universality class and one can use the effective model (which is usually much simpler than the original gauge theory) to extract informations on the deconfinement transition of the original model. These results are very general: if in particular we focus on gauge theories with a gauge group G whose center is Z 2 (like for instance, Z 2 itself, SU (2) or Sp(2N )), the deconfinement transition will belong to the Ising universality class. A well studied example of this correspondence is the SU (2) model [30,31,32]. Even if the gauge group is not SU (3) and the model only contains gluonic degrees of freedom this simplified model shares a lot of properties with QCD: the presence of a confining flux tube at low temperatures, a deconfined phase at high temperature, a rich glueball spectrum, asymptotic freedom. For these reasons it has been studied a lot in the past both in (2+1) [33,34,35] and in (3+1) [36] dimensions. Thanks to the Svetitsky-Yaffe construction, several gauge invariant observables of the SU (2) model can be mapped to equivalent Ising observables • The Polyakov loop is mapped to the spin (Z 2 odd) operator and thus the Polyakov loop susceptibility is mapped to the magnetic susceptibility χ of the Ising model.
• The deconfining transition of the gauge model corresponds to the magnetization transition of the Ising model. In particular, the confining phase (low temperature of the gauge model) is mapped to the Z 2 symmetric phase (high temperature phase) of the Ising model, while the deconfined phase is mapped to the broken symmetry phase of the Ising model.
• The Wilson action (i.e. the trace of the ordered product of the gauge field along the links of a plaquette) is mapped to the energy operator 2 of the Ising model.
• The screening mass of the gauge model in the deconfined phase is mapped to the mass of the Ising model in the low temperature phase, while N t σ, where σ is the string tension and N t is the inverse temperature of the gauge model (i.e. the size N t of the lattice in the compactified direction which defines the finite temperature setting in LGTs) is mapped to the inverse of the high temperature correlation length of the Ising model. This mapping has been widely used in the past years to predict the behaviour of various physical observables of the gauge model near the deconfinement transition, like for instance the short distance behaviour of the Polyakov loop correlator [37], the width of the flux tube [38], the Hagedorn-like behaviour of the glueball spectrum [39] or the behaviour of the universal ξ/ξ 2 ratio [40]. More generally the mapping allows one to relate all the thermodynamic observables of the two models and in particular also Ω(η), which can thus be evaluated in the SU (2) LGT and then compared with the QFT prediction discussed in the previous section.

QCD with dynamical quarks and the Columbia plot
The situation is different if we study full QCD, i.e. if we include dynamical quarks in the model. In this case the center symmetry is explicitly broken by the Dirac operator and all the above considerations do not hold anymore.
For physical values of the quark masses there is no phase transition between the high temperature quark-gluon plasma phase and the low temperature confined phase which are only separated by a smooth crossover [41,42]. However in the phase diagram of the model one finds a rich structure of phase transitions as the values of the masses of the quarks are varied [43]. This pattern of phase transitions is summarized in the well known Columbia plot, which we report here in fig. 4. The plot describes the nature of the phase transitions as a function of the quark masses. On the two axes are respectively the masses of up-down quarks (x axis) and the mass of the strange quark (y axis). We see in the central part of the plot a wide region, where the physical point lies, in which there is no phase transition but only a crossover between the two phases. In the top right and in the bottom left corners there are two regions where the transition is of the first order. These regions end with two lines of second order transitions which are expected to both be of the Ising type. The phase diagram reported in the Columbia Plot can be studied with standard Montecarlo simulations and in the vicinity of the critical lines the results of these simulations could be mapped using our tools to the Ising phase diagram and then compared with the Ising predictions as we discussed above for the simpler case of the pure gauge SU (2) model.
It is important to notice that the two critical regions are of different nature. The one in the top-right corner is a deconfinement transition similar to the one discussed in the previous section. In fact in the limit of infinite mass quarks the model becomes a pure gauge theory. In this limit the SU (3) gauge model has a first order deconfinement transition (differently from the SU (2) one discussed above which is of second order). As the mass of the quarks decreases the gap in the order parameter decreases and the first order region ends into a critical line of Ising-like phase transitions.
The critical region in the bottom left portion of the Columbia Plot has a completely different origin. It describes the restoration of the chiral symmetry in QCD at finite temperature and small quark masses. In QCD with three massless quark flavours the chiral phase transition is expected to be first order and to remain of first order even for small but non-zero values of the quark masses. As the quark masses increase, the gap in the order parameter decreases and the first order region terminates in a critical line of the Ising type. Above this line chiral symmetry is restored through a smooth crossover. In this case, the Z 2 symmetry is an "emergent" symmetry and the identification of the H and t axes of the equivalent 3d Ising model is not trivial (see a discussion on this issue in [44,45] and in sect. 5.3 below) and thus a universal charting of the scaling region could indeed be useful.
Even if the precise location of the critical line is still debated, it seems that the physical point is not too far from this bottom left Ising line. If this is the case, then our analysis could be applied, maybe with some degree of approximation, also to the physical point.

The critical ending point of the QCD phase diagram at finite chemical potential
The most interesting application of our results is for QCD at finite baryon density, which is realized by adding a finite chemical potential µ to the QCD lagrangian. This regime is particularly interesting since it can be explored experimentally in heavy-ion collisions [46,47,48,49,50,51,52,53] and at the same time it cannot be studied using Montecarlo simulations due to the well known sign problem (see for instance [54,55] for a discussion of the sign problem in this context).
In this regime the QCD phase diagram is expected to reveal interesting novel phases [56,57]. In particular it is widely expected that the hadronic phase (low T , low µ) should be separated from the quark-gluon plasma phase (high T , high µ) by a line of first order transitions with a critical endpoint at finite critical values of T and µ (see fig. 5) which should again belong to the Ising universality class [58,59,60,61,62,63,64,65]. Also for this model, as for the liquid-vapour transition (or the chiral transition discussed above), the Z 2 symmetry is not realized explicitly but it is just an "emergent" symmetry. This class of models is typically addressed with the "mixing-of-coordinates" scheme proposed in [16]. The approach was pursued for the finite density QCD case in [61] and [63], leading to very interesting results. In both cases the mapping between Ising and QCD coordinates was performed via the parametric representation of the equation of state. In this respect, the explicit expression of Ω in terms of θ that we discussed in this paper could be used as a shortcut in the process and could facilitate the identification of Ising-like behaviours in the experimental data.
As more and more experimental results are obtained, it will become possible to directly test them with universal predictions from the Ising model and it will be important to have a precise charting of the Ising phase diagram to organize results and drive our understanding of strongly coupled QCD in this regime. Our paper is a first step in this direction. We proposed and studied one particular combination, chosen for its simplicity from a theoretical point of view and its accessibility from an experimental and numerical point of view, but other combinations are possible and could be studied using, as we suggest here, parametric representations in d = 3 and/or expansion around the exactly integrable solutions in d = 2.  x = 1 A.2 Numerical estimates of the amplitude ratios in d = 3 In three dimensions there are no exact results but, thanks to the recent improvement of the bootstrap approach [6,7,66] and to the remarkable precision of recent Montecarlo simulations [11,13,15], reliable numerical estimates for all these quantities exist x = 1.412625 (10) x σ = 0.5181489(10) It is also possible to choose realizations of the model in which the amplitude of the first irrelevant operator is tuned toward zero, thus allowing a more efficient approach to the fixed point. This is for instance the idea followed in [9,10,11,12,13,14,15] to improve the numerical estimates of universal quantities using Montecarlo simulations. The scaling parameter η is defined modulo a non-universal constant h 0 which sets its scale and depends on the specific model at hand, i.e. on the specific realization of the Ising universality class in which one is interested. However, given such a realization, it is rather easy to fix h 0 . The simplest option is to measure (numerically or experimentally) the magnetization M and the susceptibility χ along two directions (or in the same direction) in the (t, H) plane, for instance along the critical line of first order phase transitions.
Then from the ratio of the two amplitudes one can easily extract h 0 . We report here for completeness the corresponding expressions in the case in which one measures, besides the amplitude B of the magnetization, the value of Γ + or that of Γ − or For instance, in the case of the 3d Ising model defined by eq.s 2, 3 we have [67] B = 1.6920(5) and Γ − = 0.2394 (5) to which corresponds h 0 ∼ 0.923 while in the case of the model tuned so as to eliminate the first irrelevant operator discussed above we find [13] B = 1.9875(3) and Γ + = 0.14300 (5) to which corresponds h 0 ∼ 8.62.

B.2 Parametric representation of the magnetic susceptibility
From the parametric representation of the critical equation of state (eq. 11) we may obtain the magnetic susceptibility as follows χ(R, θ) = m 0 h 0 R −γ 1 − θ 2 + 2βθ 2 2βδθh(θ) + (1 − θ 2 )h (θ) (32) and from this expression it is easy to extract the amplitude ratio B.3 Expansion of θ as a function ofη in the neighbourhood of the three singular pointsη = (±∞, 0) We report here the first few terms: in the θ → θ 0 limit in the θ → 1 limit and finally in the θ → 0 limit

C Exact results in d = 2
We report here the expansions, in the three regimes of interest, for χ and for M 1 ≡ 1/ξ, obtained using the results reported in [5,26,27,28].    Table 3: Expansion coefficients of the correlation length in the three regimes of interest, according to eq.s 41.