1. Introduction
Dual superconductor is known to be one of the earliest and well established models of quark confinement in QCD [
1,
2] (for a review, see [
3]). It suggests a model for the Yang–Mills vacuum which is based on the condensate of a magnetically charged Higgs field. The corresponding dual Abelian Higgs model (DAHM) represents a 4D relativistic generalization of the Landau-Ginzburg theory of dual superconductivity. This scenario of confinement is fully supported by the dedicated lattice simulations [
4,
5]. One has to realize, though, that the monopole condensation provides confinement of only the diagonal degrees of freedom, through the violation of the Abelian Bianchi identities. In this sense, the dual-superconductor model differs from the Yang–Mills theory, where confinement (of both diagonal and off-diagonal degrees of freedom) does not require violation of either Abelian or non-Abelian Bianchi identities. That is, in the dual superconductor, confinement is described in terms of a field dual to that of the dual vector boson, i.e., an antisymmetric-tensor field
which couples to a world sheet bounded by the contour of the ’t Hooft loop (for a review, see [
6]). Unlike the gauge-field strength tensor, the field
contains an additional term
, i.e.,
, which violates the Abelian Bianchi identities. This term yields the tensor structure
in the correlation function
, which, in turn, yields confinement. (The same is true for compact QED [
7].) Rather, the non-Abelian field-strength tensor
does not contain a term
, and yet the correlation function
contains the tensor structure
[
8,
9] (for a review, see [
10]). Consequently, in the non-Abelian case, one can still fix a certain gauge, e.g., the Fock–Schwinger gauge which expresses
in terms of
(to obtain the propagator
with the confinement effects, which can further be used in the non-local Nambu–Jona-Lasinio model [
11] corresponding to QCD with light quarks in the Gaussian approximation), whereas this is no longer possible in a confining Abelian theory like the dual superconductor or compact QED.
The fundamentally important problem for the description of confinement is the construction of a string represenation of the Wilson loop in the Yang–Mills theory, or of the ’t Hooft loop in DAHM (for a review, see e.g., [
6]). The partition function of DAHM, in the London limit, admits such a representation in terms of closed dual Abrikosov–Nielsen–Olesen strings. Indeed, as the dual Higgs field has singularities at the cores of those strings, one can unambiguously reformulate the integration over the phase of the Higgs field in terms of the integration over the strings’ world sheets. Moreover, it turns out that the Jacobian, corresponding to this change of variables in the functional integral, yields the Polchinski–Strominger term [
12] with such a coefficient that the string conformal anomaly (for a review, see [
13]) gets cancelled out in 4D rather than 26D [
14]. While this is certainly a successful approach to the construction of a fully quantum string representation of the DAHM partition function, it remains unclear whether, in the Yang–Mills theory, a representation of the Wilson loop in terms of an integral over open-string world sheets can be constructed in a similar way. The main obstacle in this way is the identification, in the Yang–Mills integration measure, of those degrees of freedom that can be related to the string’s world-sheet coordinates. In a recent paper [
15], an alternative approach has been proposed, where the Polchinski–Strominger term emerges from the quartic cumulant in the cumulant expansion of the Wilson loop. That has been done in the 3D Yang–Mills theory at finite temperature, with the correlations between the elements of the string’s world sheet modeled by those of the London limit of DAHM, which can be calculated analytically. The corresponding finite-temperature compactification of the theory is important for the generation of the Polchinski–Strominger term and the cancellation of the string conformal anomaly at a certain temperature, which turns out to be approximately twice smaller than the deconfinement critical temperature. In the first part of the present paper, this approach will be generalized to the case of the 4D Yang–Mills theory compactified on a torus, so that it can, in particular, become applicable to the zero-temperature theory in a space with two compact spatial dimensions.
In the second part of the paper, another DAHM-related issue will be addressed, which is the calculation of the correlation length in that theory away from the London limit. In this way, an analogy between the dual representations of the London limit of DAHM and the low-energy limit of 4D compact QED turns out to be helpful, because of the known Higgs-induced corrections to the 3D compact QED, which take place when the latter theory emerges from the 3D Georgi–Glashow model. In the Summary section, some concluding remarks will be presented.
2. Quantum-String Representation of the Wilson Loop in the Yang–Mills Theory Compactified on a Torus
Let us consider the Yang–Mills theory compactified on a torus, which is, e.g., the case if this theory is compactified on
[
16] and further taken at finite temperature [
17], or just in the continuum limit of a lattice version of the theory with periodic boundary conditions (for a review, see e.g., [
18]). In such a compactified version of the Yang–Mills theory, we are interested in the string representation of a Wilson loop corresponding to a static quark and an antiquark separated along one of the non-compact dimensions (as opposed to the so-called torelon, which is a closed string emerging in case when that dimension is compact [
19]). Let us further model the confining part of the two-point correlation function of gluonic field strengths by the corresponding expression for the dual field strengths in the London limit of DAHM, which reads [
20,
21] (for a review, see [
22])
, where
. Here,
is the mass of the dual vector boson, with
being the magnetic coupling constant, and
being the v.e.v. of the dual Higgs field. Note that
, where:
is just the massive-boson propagator, which obeys the equation
, and
henceforth denotes a Macdonald function, which is a modified Bessel function of the second kind. Due to the compactification, we further have
. Here,
is a 2D vector in the non-compact dimensions,
and
are the circumferences of the two circles which are forming the torus, and
n and
are the winding numbers. We further use integral representation
, along with the following generalization of the Poisson sum formula, which can be considered as a discrete version of the Gaussian integral [
23]:
The subsequent
t-integration yields
where
and
. Applying now the same generalized Poisson sum formula once again, in the form
we obtain the following double-sum representation for the function
:
As long as
, this expression can be approximated as:
since all the terms with a non-zero
k or
become exponentially suppressed.
By equating to each other the string tensions in the DAHM in the London limit and in the Yang–Mills theory in the Gaussian approximation, we will now establish a correspondence between these two theories. The string tension in the DAHM in the London limit, corresponding to the string interconnecting two unit electric charges, stems from the non-local string action
and reads [
20,
21,
22]
with
M standing for the Higgs mass. The last equality here is valid within the leading logarithmic approximation to the London limit [
24,
25], which assumes that not only
, but also
. On the other hand, in the Yang–Mills theory in the Gaussian approximation, the non-local string action, with the exponential parametrization for the function
, reads [
8,
9,
10]
In this expression,
N is the number of colors,
g is the Yang–Mills coupling, and we have directly set the Yang–Mills vacuum correlation length equal to
. Furthermore, parameter
, which determines the mixing between the confining- and the non-perturbative non-confining contributions to the correlation function
, is numerically equal to
in the SU(3)-case of interest [
26], so we set
. The resulting string tension reads
, yielding the following correspondence:
Let us now replace
by
in the parametrization of the correlation function
. That yields for the normalization coefficient
A:
Following the approach of Ref. [
15], we will further make the same replacement in the confining structure of the quartic cumulant [
27],
where
are the relative coordinates of the points
, and
Owing to the proximity of
to
and of
to
, ensured by the factor
, we have:
By using Equations (
2) and (
4), we further have for this expression:
Setting now
and
, we obtain the following Polchinski–Strominger term [
12] in the string representation of a Wilson loop (cf. Reference [
15]):
Here,
and
are the expressions for the Laplacian
and the scalar curvature of the world sheet in the conformal gauge, where the induced metric
has the diagonal form
, and
with
. Furthermore,
is the Green’s function of the Laplacian, and we have denoted
and
. We have also used the facts that, for
,
at
, and that
.
For the Polchinski–Strominger term (
6) to cancel the bosonic-string conformal anomaly in 4D, the coefficient in front of the integral in Equation (
6) should be equal to
(cf. Refs. [
6,
15]). This yields the following relation for the circumferences of the compactification circles, whose fulfillment provides the cancellation of the conformal anomaly:
The gluon condensate,
, calculated from the aforementioned expression for the string tension,
, by means of the phenomenological SU(3) values
and [
28]
, reads
. That yields the numerical relation
. Accordingly, we have
, which justifies approximation (
2).
Note also that there is no contradiction with the fact that the compactification radii,
and
, can be smaller than the vacuum correlation length
, since the Wilson loop is extending in the non-compact dimensions. As for the string, its existence requires to have the size of at least one of the two compact dimensions larger than the vacuum correlation length. For example, in the case when one of the two compact dimensions is
, which corresponds to finite temperature
, the spatial dimension with compactification radius
shrinks in the zero-temperature limit. Nevertheless, the string is still well defined in the resulting (2+1)D Yang–Mills theory. Note also that the limit of the small radius of spatial compactification is precisely the one where the Yang–Mills theory compactified on
admits an analytic description of confinement [
16]. In the opposite limit of large temperatures,
, with
standing for the deconfinement critical temperature, where the Yang–Mills theory effectively becomes purely spatial, the obtained condition (
7) can also be fulfilled, making the string representation of the corresponding spatial Wilson loop fully quantum. Indeed, in that limit, the chromomagnetic condensate is proportional to the fourth power of the inverse chromomagnetic-vacuum correlation length,
, i.e., [
29,
30]
, where
g now stands for the finite-temperature Yang–Mills coupling, and
i and
j are the purely spatial indices. Therefore, in this high-temperature limit, condition (
7) parametrically becomes
. As follows from this relation, both
and
vanish at
, but
stays parametrically larger than the vacuum correlation length
, since
. We notice that the thus emerged ultrasoft scale
defines the so-called sound attenuation length in the gluon plasma (see e.g., [
31,
32,
33]). In particular, it can be shown [
34,
35] that the correlation function
with the correlation length defined by that scale, yields the known perturbative expression for the shear viscosity of the plasma [
36],
. The unexpected emergence of this ultrasoft scale, within the present context of the string representation of a spatial Wilson loop in the high-temperature limit, can be viewed as a remarkable example of possible interrelations between
a priori different areas of physics of the Yang–Mills fields.
3. DAHM away from the London Limit: Vacuum Correlation Length and the Distribution of Monopole Currents
The partition function of 4D compact QED in the continuum limit has the form [
22]:
where
is the dual gauge field,
is its field-strength tensor,
,
is the UV cutoff,
is the monopole fugacity of dimensionality
, the factor of 2 in front of it is due to the contributions of both monopoles and antimonopoles, and
. The Debye mass of the field
, stemming from this expression, is
. In the low-energy limit, which corresponds to the limit of low monopole densities, dualization of Equation (
8) yields the following representation of the theory in terms of an antisymmetric-tensor field
[
22]:
where
is the strength tensor of the field
.
At the same time, the partition function of the London limit of DAHM (with closed dual Abrikosov–Nielsen–Olesen strings disregarded for this comparison, due to the absence of closed strings in compact QED) has the form [
22]:
Thus, the two theories are equivalent to each other, with the equality
yielding an expression for the fugacity,
where we have naturally set
.
Now, by an analogy with the 3D compact QED, we will proceed beyond the London limit of the 4D DAHM. In the case when the 3D compact QED is stemming from the 3D Georgi–Glashow model, for a not-infinitely-heavy Higgs field in the latter model (as compared to the mass of the W-boson), monopoles start interacting with each other via the Higgs field, so that their Coulomb interaction acquires a correction [
37],
. Here,
and
are the charges (in the units of
) of two monopoles located at the points
and
,
is the dual-photon propagator, and
is the Higgs-boson propagator. We therefore expect that, in the 4D DAHM away from the London limit, the interaction between two monopoles, propagating along the world lines
and
, becomes
, where
,
is given by Equation (
1), and we have taken into account that the (dual) Higgs is a scalar particle. Hence, we can write down the following partition function of the DAHM away from the London limit (cf. the corresponding partition function stemming from the 3D Georgi–Glashow model with a not-infinitely-heavy Higgs, considered in Ref. [
38]):
where we have once again set
, in
. The latter choice is natural, since the typical length of a monopole trajectory,
, is
.
We further perform the
-integration in Equation (
9), by using the cumulant expansion in the Gaussian approximation (cf. Refs. [
39,
40]; for a review, see [
41]), i.e.,
where
and
. Averages entering this expression read
and
. As we are always considering constant
m, i.e., a constant Higgs v.e.v.
, by staying outside the cores of the dual Abrikosov–Nielsen–Olesen strings,
, which corresponds to having
, we should replace 0 by
in
. This replacement yields
. For the same reason of considering only the distances
, we can, with the exponential accuracy, approximate
. This yields for Equation (
9):
where
, since
. Expanding now the product of the two cosines up to the quadratic term,
, we obtain the following corrected squared Debye mass (cf. Refs. [
39,
40,
41] for the corresponding expression in the 3D Georgi–Glashow model):
Thus, with a deviation from the London limit, the DAHM vacuum correlation length
, equal to
, acquires a negative correction
. Note that this decrease of the vacuum correlation length is in line with the attractive force, which is known to start acting onto two parallel same-oriented Abrikosov–Nielsen–Olesen strings away from the London limit due to the Higgs-boson interactions (in addition to the repulsive force, which is produced by the gauge-boson interactions) [
42]. Such a competition between the forces which the gauge-boson and the Higgs-boson interactions excert onto two interacting topological defects is similar for Abrikosov–Nielsen–Olesen strings and for monopoles (cf. above).
Note also that the obtained decrease of the vacuum correlation length in the vicinity of the London limit parallels its strong decrease in the Bogomolny limit, where the gauge- and the Higgs-boson masses are equal to the same quantity
[
43]. Indeed, the string tension in the Bogomolny limit differs from the one in the London limit, Equation (
3), by the absence of the factor
[
25], i.e., it reads:
Since, within the dual superconductor model, both Equations (
3) and (
11) should have the same aforementioned phenomenological value
, we obtain:
Accordingly, the vacuum correlation length in the Bogomolny limit,
, is smaller than the vacuum correlation length in the London limit,
, by the same large factor of
.
Let us now dualize theory (
10), which will yield a potential of monopole currents away from the London limit. To this end, we write
, where the first term in the exponential on the r.h.s. represents Coulomb interaction of monopole currents, while the dual field
has no kinetic term anymore, and can thus be integrated over in the saddle-point approximation. For a rescaled field
, the action has the form
Assuming
of the form:
, we obtain the following saddle-point equation:
Seeking its solution in the form
, where
, we obtain in the leading approximation:
,
. This solution to the saddle-point equation yields the following potential of monopole currents:
where
. Retaining further the linear-in-
terms in the saddle-point equation, we obtain its following solution:
. We explicitly see that the inequality
holds at low monopole densities,
, in which case
. Accordingly, in the corresponding expression for the potential,
only the quadratic term should be retained in this low-density approximation. That yields
Here,
consistently substitutes
which enters the same expression for
in the case where the correction is not taken into account, i.e.,
.
4. Summary
This paper has been devoted to the two aspects of the dual superconductor scenario of confinement. The first one was the derivation of the string representation of a Wilson loop in the 4D Yang–Mills theory compactified on a torus. For the Gaussian approximation to the Yang–Mills theory, we have used the explicit form of the two-point correlation function of gluonic field strengths, which can be obtained in the London limit of the dual superconductor. With a deviation from the Gaussian approximation, the quartic cumulant of gluonic field strengths starts contributing to the Wilson loop. That cumulant leads to the appearance, in the string representation, of a term corresponding to the two-point correlation of the Euler characteristic of the string world sheet. At short distances, this correlation gets reduced to a massless 2D propagator, yielding the Polchinski–Strominger term (
6) in the string representation. When relation (
7) between the circumferences of the compactification circles holds, that term cancels string conformal anomaly in 4D, thereby making string representation fully quantum. While we find this result interesting, we notice that its derivation relies on the mentioned explicit form of the two-point correlation function of gluonic field strengths which corresponds to the London limit of the dual superconductor. Indeed, only that correlation function yields, upon the compactification of the theory on the torus, the correlation function (
2), whose short-distance logarithmic asymptotics eventually yields the Polchinski–Strominger term. We notice here one more remarkable fact stemming from the logarithmic short-distance asymptotics of the Macdonald function
which describes also the electric field of a straight-line dual Abrikosov–Nielsen–Olesen string, i.e., a flux tube, in the London limit of the dual superconductor. Namely, pairs produced in the field of such a flux tube have an exponential distribution, in contrast to the standard Schwinger’s Gaussian distribution, which takes place away from the London limit [
44,
45].
The second aspect of the dual superconductor scenario, addressed in this paper, was the deviation from the London limit. We have used the correspondence between the dual formulations of DAHM in the London limit and of the low-energy limit of 4D compact QED, to write down partition function (
9) of DAHM away from the London limit. It is similar to the partition function of the grand canonical ensemble of monopoles in the 3D Georgi–Glashow model where the Higgs field is considered not infinitely heavy. Having integrated, in the Gaussian approximation, over the scalar field which represents in this formulation the dual Higgs field, we have obtained a decrease of the vacuum correlation length of DAHM away from the London limit, by a factor of
, where
m and
M are the masses of the dual vector boson and the dual Higgs field, respectively. Such a decrease parallels the known fact about an attractive force caused by the Higgs-boson exchanges between two parallel same-oriented straight-line dual Abrikosov–Nielsen–Olesen strings, as well as the strong decrease of the vacuum correlation length in the Bogomolny limit. Finally, by dualizing partition function (
9), we have obtained its representation in terms of monopole currents, which includes Coulomb interaction of those currents along with their corresponding potential.