Gluon Propagators in $2+1$ Lattice QCD with Nonzero Isospin Chemical Potential

The static longitudinal and transverse gluon propagators in the Landau gauge are studied in $2+1$ lattice QCD with nonzero isospin chemical potential $\mu_I$. Parameterization of the momentum dependence of the propagators is provided for all values of the chemical potential under study. We find that the longitudinal propagator is infrared suppressed at nonzero $\mu_I$ with suppression increasing with increasing $\mu_I$. It is found, respectively, that the electric screening mass is increasing with increasing $\mu_I$. Additionally, we analyze the difference between two propagators as a function of the momentum and thus compare interactions in chromoelectric and chromomagnetic sectors.

It is also believed that the study of QCD at nonzero µ I might help, in one way or another, to solve more difficult problem of QCD at nonzero baryon chemical potential. First, QCD with nonzero µ I is an ideal test system for some lattice approaches to QCD with nonzero baryon chemical potential such as Taylor expansion, analytical continuation and others. Second, lattice results obtained at nonzero µ I by means of first principle calculations can be compared with respective results of the effective field theories of QCD and phenomenological models. Such comparison allows to test the effectiveness of QCD models and thus to estimate their applicability to the study of QCD at nonzero µ B .
The gluon propagators are among important quantities to study. Contemporary perturbative computations of the equation of state of high-density strong-interacting matter employ infrared behavior of gluon propagators described in terms of screening (see, e.g. [32] and references therein). Nonperturbative estimates of such screening effects on the basis of lattice simulations may be used for testing the methods of resummation of perturbation series.
In early calculations of pressure at nonzero quark chemical potentials [33], the full gluon propagators appear in the expressions for the thermodynamical potential Ω(T, V, µ). If the Green functions that appear in such expressions were evaluated nonperturbatively, the dynamics of pressure build-up and other thermodynamic quantities could be investigated in detail.
The gluon propagator results also will be useful in testing future results obtained at nonzero µ I by the functional methods like Dyson-Schwinger equation and functional renormalization group.
Landau gauge gluon propagators in the non-Abelian gauge theories at zero and nonzero temperature were extensively studied in the infrared range of momenta by various methods. We shall note lattice gauge theory, Dyson-Schwinger equations, Gribov-Zwanziger approach. At the same time the studies at nonzero quark chemical potential are restricted to a few papers only.
Here we close this gap for the case of nonzero isospin chemical potential. Let us note that for the lattice QCD at nonzero baryon chemical potential the results can be obtained only for small values of µ B because of the sign problem [34].
It is known that in QCD with nonzero µ I there is a phase transition from the normal to the superfluid phase at µ I = µ π /2. Recently results were obtained [10] indicating that at large µ I the theory undergoes smooth transition to the deconfinement phase. It is interesting to check if the gluon propagators change at respective values of µ I .
In this paper we make the first study of the influence of µ I on the Landau gauge gluon propagators using lattice QCD approach. We use the same lattice action as in Ref. [10] and in fact the same set of the lattice configurations. Our goal is to study how the gluon propagators change with an increase of µ I . In particular, we find strong suppression of the (color)-electric (longitudinal) gluon propagator with increasing µ I which is reflected in increasing of the respective screening mass. We parameterize the propagators in the wide range of momentum values using the fit function which is the tree level prediction of the Refined Gribov-Zwanziger approach. We also compute an observable introduced recently in [35], a difference between the (color-)electric and (color-)magnetic propagators, and obtain its dependence on the momentum and isospin chemical potential.
The paper is organized as follows. In Section II we specify details of the lattice setup to be used: lattice action, definition of the propagators and details of the simulation. In Section III we consider the propagators at low momenta to compute respective screening masses and determine their dependence on µ I . In Section IV the fit over a wide range of momenta is performed using Gribov-Stingl fit function motivated by the Refined Gribov-Zwanziger effective action. We demonstrate that for the longitudinal propagator this fit works well over the whole momenta domain under study for all values of µ I considered, whereas for the transverse propagator it works only at small µ I .
In Section V the dependence of the difference between the electric and magnetic propagators on the momentum and isospin chemical potential is discussed. We show that this difference decreases exponentially fast with momenta similarly to the case of SU (2) QCD [35]. The last section is devoted to the discussion of the results and to conclusions to be drawn.

II. SIMULATION DETAILS
We carry out our study of the 2 + 1 lattice QCD with light quark mass m u = m d ≡ m ud and strange quark mass m s using 28 4 lattices for a set of the chemical potentials in the range aµ I ∈ (0, 0.42). The tree level improved Symanzik gauge action [36] and the staggered fermion action with a pion source term were used to generate the lattice gauge field configurations [10].
The partition function (after integration over the quark fields) is of the form [7][8][9][10] where λ is coupling of the pionic source term, light quarks lattice operator M ud and strange quark lattice operator M s are / D(µ I ) is the staggered lattice Dirac operator with quark chemical potential µ I = µ u = −µ d . The unphysical pionic source term which breaks U (1) symmetry explicitly is necessary to observe the spontaneous breaking of this symmetry in a finite volume in the limit λ → 0 at high enough µ I [1]. The lattice configurations were generated [10] at β = 4.036, am ud = 0.0077, am s = 0.0271, λ = 0.5am ud . The lattice spacing of the ensemble in physical units was determined via the Sommer scale value r 0 = 0.468 fm [37] to be a = 0.0687 fm. Then the pion mass in physical units is m π = 380 MeV [10]. Our lattice spacing is small enough to be sure that usual lattice discretization effects are small. However, at nonzero isospin chemical potential new lattice artifacts may appear. To reach large quark densities without lattice artifacts one needs sufficiently small lattice spacing to satisfy condition aµ I 1. Actual values of aµ I under consideration are given in Table I of Appendix A, one can see that aµ I < 0.42. Thus we can hope that respective lattice artifacts are small.
At the same time to study the gluon propagators in the infrared region it is necessary to employ large physical volume. As a result of a compromise between these two requirements our lattice size is rather moderate: L = 28a = 1.92 fm. This implies a potential problem of large finite volume effects at small momenta (the minimal momentum p min ≈ 650 MeV).
Our lattice with equal number of temporal and spatial sites N s = N t = 28 corresponds to the temperature T = 1 28a ≈ 100 MeV. However, the energy of the lightest excitation at µ I = 0 is estimated as E exc ≥ p min = 650 MeV. Since E exc T , this temperature can be considered as approximately zero. At the greatest isospin chemical potential under study, µ I = 1.2 GeV, the energy of the lightest excitation can be estimated as follows. Using the model of free quarks, we notice that in the respective ground state all modes with | p| < µ I are occupied by ud quark states (27 modes with | p| ≤ p min √ 3 in the case under consideration). The lightest excitation occurs when one light quark jumps from the mode with | p| = p min √ 3 to the mode with | p| = 2p min . The energy needed for this purpose is E exc ∼ 170 MeV, therefore, in a rough approximation, the temperature T ∼ 100 MeV can be thought of as zero even at µ I = 1.2 GeV.
In our study of the gluon propagators we employ the standard definition of the lattice gauge vector potential A x,µ [38]: The lattice Landau gauge fixing condition is which is equivalent to finding an extremum of the gauge functional with respect to gauge transformations ω x . To fix the Landau gauge we use the simulated annealing (SA) algorithm with finalizing overrelaxation [39,40]. The gluon propagator D ab µν (p) is defined as follows: ] and the physical momenta p µ are defined by the relations ap i = 2 sin (πq i /N s ), ap 4 = 2 sin (πq 4 /N t ); p = 4 µ=1 p 2 µ . At nonzero µ I the O(4) symmetry is broken and there are two tensor structures for the gluon propagator [41] : where (symmetric) orthogonal projectors P T ;L µν (p) are defined for p = ( p = 0; p 4 = 0) as follows Therefore, two scalar propagators -longitudinal D L (p) and transverse D T (p) -are given by For p = 0 they are defined as follows: D T (p) is associated with the magnetic sector, D L (p)with the electric sector. We consider the soft modes p 4 = 0 and use the notation D L,T (p) = D L,T (0, | p|).

III. SCREENING MASSES
Various definitions of the screening masses were discussed in great detail in our previous studies [35,42].
Here we use the definition in terms of the Taylor expansion of D −1 L,T (p) in p 2 at low momenta: This method was used in [43] in the studies of lattice QCD at finite temperature and was applied to QC 2 D in [42]. One can note that the Yukawa-type fit-function successfully employed in the previous studies of lattice gluodynamics at zero and finite temperatures in Refs. [44][45][46] does not work in the case under consideration because the minimal nonzero momentum on our lattices is rather large: p min ≈ 640 MeV. Thus we have no data points in deep infrared domain, where higher-order terms in formula (15) can be neglected. As explained below, we have to take into account up to the O (p 2 ) 3 terms of the Taylor expansion of the inverse propagators.
It should be noted that the screening massesm 2 E,M are related to the correlation lengths: where ξ E,M are usually defined in terms of propagators as follows [47]: To extract the screening masses from the data, we employ the fit function (15) keeping terms up to p 2 or (p 2 ) 2 or (p 2 ) 3 , depending on the fit range. We will call below respective fit functions as fit function with two terms, with three terms or with four terms. We perform the fit over the range p low ≤ p < p cut , where p low is equal to either zero or p min and p cut varies from 1.3 GeV to 2 GeV.
We find that, in the general case, both m E and m M depend substantially on the choice of p low , p cut and the fit function. To improve stability of the results, we consider the propagators as functions of 2 variables, D −1 L,T (p, µ I ). In so doing, we assume smooth dependence of a propagator on both variables and fit the overall set of data by a polynomial in p 2 and µ 2 I thus reducing drastically the  total number of estimated parameters. We choose the polynomial by selecting a finite number of terms in the series we omit the superscript L or T of the coefficients c i,j .
Odd powers of µ I are prohibited by isospin symmetry of strong interactions and, in particular, gluon propagators are invariant under the change µ I → − µ I . Such polynomials providing a reliable quality of the fit are found by trial-and-error method in the longitudinal case, the resulting expressions for each fit range are given below. The results of this (that is, combined) and the direct (that is, by formulas like (15)) fits for the inverse prop-agator at a given value of µ I are compared in Fig. 1. It should be remembered that an increase of the number of estimated parameters gives rise to increasing of uncertainties in their values, whereas a decrease of this number results in either a dramatic decrease of the fit p-value or in a substantial shrinking of the fit range.
In the case of transverse propagator both fitting procedures are not stable with respect to elimination of the zero momentum; data on larger lattices and increased statistics are needed to arrive at reliable results. In this study, screening in the magnetic sector is discussed in terms of the dressing function, see Section V.

A. Chromoelectric screening mass
We fit our data by the function (15) with three or four terms using p cut = 1.3, 1.5 or 2.0 GeV and p low = 0 or p min . The results of the fits are stable with respect to exclusion of the zero momentum only for p cut = 2 GeV. In this case, the fit function (15) with three terms works well at µ I ≥ 0.4 GeV; at 0 ≤ µ I ≤ 0.2 GeV we employ the fit function (15) with four terms. These results are shown in the right panel of Fig. 2 by squares. In performing this fitting procedure we determine in total 23 fit parameters (when all values of µ I are taken into account).
As we explained above one can substantially reduce the number of fit parameters by considering the inverse propagator as function of p and µ I . This approach also makes it possible to derive m 2 E as function of µ I and estimate respective confidence interval for each value of µ I ∈ (0, 1.2 GeV).
An appropriate goodness of fit χ 2 /34 = 1.66, pvalue ≈ 0.009 can be achieved for the momentum range 0 ≤ p < 1.5 GeV with the following 15-parameter fit function The error in this quantity which takes into account the correlations between the fit parameters was calculated with the use of the REDUCE computer algebra system.
The respective corridor of errors is shown by the error bars in Fig.2.
It should be emphasized that the fit (20) (fit B) at both values of p cut and the fit (15) (fit A) at p cut = 2.0 GeV are stable with respect to an exclusion of p = 0 data. This indicates that the dynamics below p min has only a little effect on the screening in the chromoelectric sector.
We show the results of fit A (15) with huge errors (for p cut = 1.5 GeV) in the left panel of Fig.2 in order to illustrate the advantages of fit B (20). Some disagreement between the results of the fit A and B is observed, especially in the right panel of Fig.2 at µ I ∼ 1.0 GeV (p cut = 2.0 GeV). It should be emphasized that, in performing fit A, the term ∼ p 8 is not taken into account when p cut = 2.0 GeV and the terms ∼ p 6 and ∼ p 8 are not taken into account when p cut = 1.5 GeV, whereas all these terms are included in fit function (20) in both cases. The mentioned disagreement indicates that these terms play a significant role at p > 1 GeV. In view of this facts we conclude that the results obtained with fit B are more reliable provided that m 2 E is a smooth function of µ I . The systematic error of fit B can be estimated by a comparison of the error corridors for the cutoff momenta p cut = 1.5 GeV and p cut = 2.0 GeV as well as the error corridor for the case when the zero momentum is eliminated. It is smaller than the statistical error.
One can see in Fig.2 that the chromoelectric screening mass remains approximately constant at 0 ≤ µ I ≤ 0.4 GeV and increases with µ I at its greater values. The dependence of m E on µ I found here is qualitatively similar to its dependence on the temperature at T > T c both in pure gluodynamics and QCD as was demonstrated by lattice simulations in [44,46]. It is also similar to the dependence of m E on the baryon chemical potential in QC 2 D [35,42].
Recollecting the dependence of m E and m M on µ q and T in QC 2 D obtained in [48], one is tempted to assume that the screening in the magnetic sector will increase with T starting from T ∼ 200 ÷ 300 MeV at all values of µ I , whereas m E will show a pronounced increase only at small values of µ I . Testing this hypothesis should be the subject of future studies.

IV. DEPENDENCE ON THE MOMENTUM AND ISOSPIN CHEMICAL POTENTIAL
First we consider the momentum dependence of the gluon propagators for various values of µ I in more detail. The propagators are renormalized according to the MOM scheme to satisfy the condition D L,T (p = κ) = 1/κ 2 at κ = 3 GeV. In Fig 3 (left) we present the momentum dependence for the longitudinal propagator D L (p) for some values of µ I . In the infrared domain it clearly decreases with an increase of µ I . This dependence on µ I is similar to the dependence of D L (p) on T at T > T c in both gluodynamics and QCD [49][50][51] as well as on µ q in QC 2 D [35].
In Fig. 3(right) the momentum dependence for the transverse propagator D T (p) for the same values of µ I is shown. It is clearly seen that D T (p) is substantially less sensitive to an increase of µ I . In the deep infrared region µ I dependence of the transverse propagator differs substantially from that of the longitudinal propagator. It increases with µ I at small µ I until it reaches a weakly pronounced peak (at µ I 0.4 GeV for 0 ≤ p ≤ 0.6 GeV and at µ I 0.8 GeV at p = 0.9 GeV, see also  Fig.4 (right panel) by purple triangles; this behavior is also in contrast to that of the longitudinal propagator, which decreases with µ I at all momenta as is seen in the left panel of Fig. 4 It is known that at a finite temperature the propagator D T (p) has a clear maximum at the value of momentum increasing with temperature. Our data give no evidence for such maximum at a small momentum, however, we do not rule out its existence. Now we proceed to an interpolation formula for our data. It was demonstrated many times [51][52][53][54][55] that the infrared behavior of the gluon propagators at zero and finite temperature can be well described by the fit function proposed in [57] as well as the tree level prediction of the Refined Gribov-Zwanziger approach [58], This fit function is in fact identical to that used in our previous study [35], the relations between the two sets of fit parameters are obvious.
Since the minimal nonzero momentum p min for the lattice under study is rather large our fit over the infrared domain may suffer from finite-volume effects. Nevertheless, we believe that we find qualitatively correct dependence of the propagators on µ I , at least for D L at large µ I .
To lower uncertainties in the fit parameters M i , we perform fit for the ratio where p 0 should lie in the fit domain; our choice is p 0 = 2 GeV. In so doing, the normalization factors Z L,T are not determined and thus should be computed by the formula With this Z factor, formula (23) gives the propagator renormalized at p = κ. The fit domain in the longitudinal case ranges up to p cut = 5 GeV; in the transverse case the cutoff momentum for each value of µ I is indicated in Table II. A few words about the fit stability should be said. In the longitudinal case, the fit parameters are unaffected by elimination of zero momentum, in the transverse case the fit parameters change significantly at 4 of 7 values of µ I when zero momentum is excluded.
The fit parameters for D L (p) depend only weakly on the cutoff momentum, as is shown, as an example, in Table III for µ I = 0; the range 2.2GeV < p cut < 5.0 GeV is considered (at lower value of p cut the fit parameters are poorly determined). At the other values of µ I the p cut -dependence of M i is also insignificant, however, the minimum value of p cut at which the fit parameters can be determined increases with µ I and reaches 4.5 GeV at µ I = 1.2 GeV.
The fit parameters and goodness-of-fit for D T (p) depend on p cut as follows. At µ I = 0.0 and 0.2 GeV we arrive at a reasonable p-value for 1.9 < p cut < 5.0 GeV, M i depend on p cut only weakly. At µ I = 0.4 GeV, a reasonable p-value for 1.9 < p cut < 5.0 GeV also holds, however, M i depend significantly on p cut at p cut > 3 GeV. At 0.6 ≤ µ I ≤ 1.2 GeV we obtain p-value ≥ 0.05 only at p cut ≤ 2.2 GeV (this one can see in Fig. 7(right) for the transverse dressing function). At 1.9 ≤ p cut ≤ 2.2 GeV M i vary within their errors and at p cut ≤ 1.9 GeV they are poorly determined. The parameters of this fit are presented in tables I and II for the longitudinal and transverse propagators, respectively. The parameters M 2 1 , M 2 2 , and M 4 3 are given in the respective powers of GeV for the sake of comparison with [55]. Results of this fit are also shown in Fig. 3 together with our lattice data.
The dependence of these fit parameters on the isospin chemical potential is presented in Figs.5 and 6, where our results are compared with those obtained in [55] in SU (3) gluodynamics at T = 0 on a large lattice with high statistical precision.
The parameters M It is instructive to look also at the respective dressing functions J L,T (p) defined as It is seen in Figure 7 (left) that with increasing µ I the maximum of the longitudinal dressing function goes down and shifts to the right, thus approaching dressing function of a massive scalar particle. We note once more that this dependence on µ I is very similar to dependence on the temperature, see e.g. Ref. [50,56]. The transverse dressing function increases at µ I < 0.8 GeV so that it approaches its peak at p ∼ 0.7 ÷ 0.9 GeV, at larger µ I the height of the peak decreases with increasing µ I and the position shifts to p ∼ 1.0 ÷ 1.1 GeV.
The difference between the transverse and longitudinal dressing functions ∆J = J T − J L is shown in Fig. 8, right panel. It equals zero at µ I = 0 and then it rapidly increases over the range 0 < µ I < 0.4 GeV. Then, as µ I increases from 0.4 to 1.2 GeV, the height of its peak fluctuates over 0.8 < ∆J < 1.1, the width increases and the position shifts to higher values of p.
It should be emphasized that, at 0 < µ I < 0.4 GeV, when both screening masses remain constant, ∆J shows a rapid change.

V. FEATURES OF DL − DT BEHAVIOR
In the previous two sections we demonstrated that the difference between the longitudinal and transverse propagators increases with an increase of the isospin chemical potential. Therewith, they come close to each other at high momenta for a fixed value of µ I . In this section we study the rate of the decrease of ∆(p) = D T (p) − D L (p) with increasing momentum and how the picture changes with increasing µ I . A similar comparison of these two propagators was made in [35] in QC 2 D as well as in [46] in SU (3) gluodynamics at finite temperatures where the ratio D L (p)/D T (p) was considered.
Here we show that, in the theory under study, the difference between transverse and longitudinal propagators, ∆(p) = D T (p) − D L (p) at p 4 = 0 shows clear exponential dependence on p.
Our numerical results for ∆(p) are presented in Fig. 8  (left), together with the fit function The exponential decreasing is well established starting from some momentum p 0 depending on µ I . We found that 0 < p 0 < p min for 0.2 ≤ µ I ≤ 1.0 GeV and p 0 ≈ 1.2 GeV for µ I = 1.2 GeV. The upper limit of the fit domain p 0 < p < p cut is determined by the requirement that ∆(p) differs from zero by more than 3 statistical errors (see Table IV). As a check we compare the fit function (27) with the power fit function motivated by a power-like behavior of both gluon propagators when p → ∞. At 0.2 ≤ µ I ≤ 0.4 both fit functions work well because ∆(p) does not vanish for only a few momenta. At µ I > 0.4 GeV only the exponential fit function (27) works. Moreover, we strengthen this conclusion by considering three-parameter power-like fit function as follows: it also works only for 0.2 ≤ µ I ≤ 0.4, the respective pvalues as well as the fitting parameters in formula (27) are presented in Table IV. We also tried several other power-like fit functions which also have failed.
Since our data can be well approximated by the exponential function and poorly -by a power-like fit function even with 3 parameters, we arrive at the conclusion that the Gribov-Stingl fit cannot work for D L (p) and D T (p) simultaneously at large µ I . Assuming that they do work simultaneously, we obtain the power-like dependence of the difference D T (p) − D L (p), whereas, as we found, it cannot be fitted by a power-like func-tion. Of course, the above reasoning is valid provided that the fit domain is sufficiently large and statistical precision is sufficiently high to discriminate between powerlike and exponential behavior. Given statistical precision as in the case under consideration, power-like and exponential behavior are differentiated when the fit domain is larger than p min ≤ p < 2.2 GeV as one can conclude from a comparison of Table IV and Fig. 8. This is the case for µ I ≥ 0.6 GeV. Thus our conclusion that ∆(p; µ I ) c exp( − ν ·p) cannot be fitted by a power-like fit function at µ I ≥ 0.6 GeV agrees with the statement in Section IV that the Gribov-Stingl fit (23) for D T (p) does not work at these values of µ I when p cut > 2.2 GeV, whereas the same function can be fitted well to D L (p) even at p cut = 5 GeV for all µ I under consideration.
We show the curves resulting from our fit (27) to ∆(p) for some values of µ I in Fig. 8 (left panel) as well as the difference between the dressing functions (right panel). The dependence of the parameters c and ν on the isospin chemical potential is shown in Fig. 9. The position and the height of the peak of the difference between the dressing functions (30) is given by the formulas the respective values are plotted in Fig. 10.
It is interesting to notice that, at µ I ≤ 0.4 GeV, ∆J max rapidly increases, whereas p max does not change. A completely different situation occurs at µ I ≥ 0.4 GeV: ∆J max changes only slightly, whereas p max increases substantially and shows a clear linear dependence on µ I at µ I ≥ 0.6 GeV. That is, the behavior of p max is similar to the behavior of the electric screening mass (which is associated with the interaction radius), whereas ∆J max (which is associated with the strength of interaction) changes in an opposite manner.

VI. CONCLUSIONS
We presented first results of the study of the static longitudinal and transverse propagators in the Landau gauge of lattice QCD at nonzero isospin chemical potential.
Our main observations are as follows. We found that the longitudinal propagator D L (p) is more and more suppressed in the infrared with increasing µ I . This is reflected, in particular, in the increase of the chromoelectric screening mass. Such dependence of D L (p) on µ I is analogous to its dependence on temperature at T > T c [49][50][51]. It is also similar to dependence of D L (p) on the quark chemical potential µ q in QC 2 D [35]. We found much weaker dependence on µ I for the transverse propagator D T (p). The question of whether the transverse propagator in the infrared domain can be adequately described in terms of screening mass may be the subject of future studies.
We determined the parameters of the Gribov-Stingl fit function (23) associated with the Gribov-Zwanziger confinement scenario. The parameters for the transverse propagator depend only weakly on µ I and agree well with those in pure gluodynamics at T = 0, the parameters M 2 and M 3 for the longitudinal propagator increase with increasing µ I . The difference between the parameters of the longitudinal and transverse propagators approaches a local peak at µ I = 0.4 GeV. This might be caused by the transition from the normal to the superfluid phase which is expected at µ I ∼ 0.2 GeV.
We investigated the behavior of the dressing functions, the excess of which over one in the infrared region indicates the presence of long-range forces. At large values of µ I the transverse dressing function is greater than the longitudinal one indicating that the chromoelectric interactions are screened at shorter distances than the chromomagnetic interactions.
Finally, we studied the difference ∆(p) = D L (p) − D T (p) and found that it decreases exponentially with the momentum at large p similarly to the case of the QC 2 D theory [35]. We also investigated the µ I dependence of the parameters of the exponential fit function (27) as well as the parameters p max and δJ max characterizing the behavior of the difference between the longitudinal and transverse dressing functions. This dependence indicates that with increasing µ I the domain where chromomagnetic forces dominate over chromoelecric ones extends to ever shorter distances.
Our results for the gluon propagators obtained using the first principles approach of the lattice QCD should be definitely useful to test other approaches to QCD at finite isospin chemical potential like perturbation theory [32] or functional methods.     (27) to our data for ∆(p). In the last two columns the values of the goodness-of-fit parameter corresponding to (27) and (29)