Type of dual superconductivity for the $SU(2)$ Yang--Mills theory

TL;DR

This paper determines the type of dual superconductivity responsible for quark confinement in SU(2) Yang--Mills theory by solving field equations and fitting lattice simulation data, confirming it as type I and analyzing flux tube forces.

Contribution

It provides a comprehensive solution for the vortex in the U(1) gauge-scalar model and refines the classification of dual superconductivity in SU(2) Yang--Mills theory as type I.

Findings

Confirms type I dual superconductivity for quark confinement.

Improves accuracy of the Ginzburg--Landau parameter estimation.

Shows attractive forces among chromoelectric flux tubes.

Abstract

We investigate the type of dual superconductivity responsible for quark confinement. For this purpose, we solve the field equations of the $U(1)$ gauge-scalar model to obtain the static vortex solution in the whole range without restricting to the long-distance region. Then we use the resulting magnetic field of the vortex to fit the gauge-invariant chromoelectric field connecting a pair of quark and antiquark which was measured by numerical simulations for $SU(2)$ Yang--Mills theory on a lattice. This result improves the accuracy of the fitted value for the Ginzburg--Landau parameter to reconfirm the type I dual superconductivity for quark confinement which was claimed by preceding works based on a fitting using the Clem ansatz. Moreover, we calculate the Maxwell stress tensor to obtain the distribution of the force around the flux tube. This result suggests that the attractive force acts among chromoelectric flux tubes, in agreement with the type I dual superconductivity.

Equations

$$\rho [U]:=$$

$$-\frac{1}{\mathrm{tr}(1)}\frac{⟨\mathrm{tr}(U_P)\mathrm{tr}(W[U])⟩}{⟨\mathrm{tr}(W[U])⟩},$$

$$\rho [U]=$$

$$\simeq$$

$$F_{\mu \nu }[U]=\frac{\sqrt{\beta }}{2}\rho [U],\beta =\frac{4}{g^2}.$$

$$V_{x,\mu }:={\tilde{V}}_{x,\mu }/\sqrt{\frac{1}{2}\mathrm{tr}[{\tilde{V}}_{x,\mu }^{†}{\tilde{V}}_{x,\mu }]},$$

$${\tilde{V}}_{x,\mu }:=U_{x,\mu }+{\mathbf{n}}_x{U}_{x,\mu }{\mathbf{n}}_{x+\mu }.$$

$$\rho [V]:=$$

$$-\frac{1}{\mathrm{tr}(1)}\frac{⟨\mathrm{tr}(V_P)\mathrm{tr}(W[V])⟩}{⟨\mathrm{tr}(W[V])⟩},$$

$$\rho [V]\simeq$$

$$F_{\mu \nu }[V]=$$

$$k_{\mu }:=\frac{1}{2}{ϵ}_{\mu \nu \rho \sigma }{\nabla }_{\nu }{F}_{\rho \sigma }[V],$$

$$\mathcal{L}=-\frac{1}{4}{F}_{\mu \nu }{F}^{\mu \nu }+{\left(D_{\mu }\varphi \right)}^{\ast }{D}^{\mu }\varphi -\frac{{\lambda }^2}{2}{\left({\varphi }^{\ast }\varphi -v^2\right)}^2,$$

$$F_{\mu \nu }(x):=$$

$$D_{\mu }\varphi (x):=$$

$$D^{\mu }{D}_{\mu }\varphi =$$

$${\partial }^{\mu }{F}_{\mu \nu }=$$

$$j_{\nu }:=$$

$$A_0(x)=0,\mathbf{A}(x)=A(\rho ){\mathbf{e}}_{\phi },\varphi (x)=vf(\rho )e^{in\phi },$$

$$\displaystyle-\frac{1}{\rho}\frac{d}{d\rho}\biggl{[}\rho\frac{d}{d\rho}f(\rho)\biggr{]}+\biggl{[}\frac{n}{\rho}-qA(\rho)\biggr{]}^{2}f(\rho)$$

$$\displaystyle\hskip 56.9055pt=\lambda^{2}v^{2}\bigl{[}1-f^{2}(\rho)\bigr{]}f(\rho),$$

$$\displaystyle\frac{d}{d\rho}\biggl{[}\frac{1}{\rho}\frac{d}{d\rho}\left(\rho A(\rho)\right)\biggr{]}=j_{\varphi}(\rho),$$

$$\displaystyle j_{\varphi}(\rho)=2q^{2}v^{2}\biggl{[}A(\rho)-\frac{n}{q\rho}\biggr{]}f^{2}(\rho).$$

$$\mathbf{B}(x)=\nabla \times \mathbf{A}(x)=\frac{1}{\rho }\frac{d}{d\rho }\left(\rho A(\rho )\right){\mathbf{e}}_z.$$

$$T^{\mu \nu }=$$

$$+{\left(D^{\mu }\varphi \right)}^{\ast }\left(D^{\nu }\varphi \right)-g^{\mu \nu }\left(D_{\rho }\varphi \right){\left(D^{\rho }\varphi \right)}^{\ast }$$

$$+\frac{{\lambda }^2}{2}{g}^{\mu \nu }{\left(v^2-{\varphi }^{\ast }\varphi \right)}^2.$$

$$E=$$

$$=$$

$$\displaystyle+v^{2}\biggl{[}\frac{d}{d\rho}f(\rho)\biggr{]}^{2}+v^{2}\biggl{[}\frac{n}{\rho}-qA(\rho)\biggr{]}^{2}f^{2}(\rho)$$

$$\displaystyle+\left.\frac{\lambda^{2}v^{4}}{2}\bigl{[}1-f^{2}(\rho)\bigr{]}^{2}\right\}.$$

$$f(\rho )=1,A(\rho )=\frac{n}{q\rho },$$

$$f(\rho )\stackrel{\rho \to \infty }{\to }1,A(\rho )\stackrel{\rho \to \infty }{\to }\frac{n}{q\rho },$$

$$f(\rho )\stackrel{\rho \to 0}{\to }0,A(\rho )\stackrel{\rho \to 0}{\to }0,$$

$$\Phi :=$$

$$=$$

$$A(\rho )=\frac{n}{q\rho }a(\rho ).$$

Full text

††preprint: CHIBA-EP-236, KEK Preprint 2018-82, 2019.09.02

I Introduction

In high energy physics, quark confinement is a long-standing problem to be solved in the framework of quantum field theories, especially quantum chromodynamics (QCD). The dual superconductivity picture dsp1 ; dsp2 ; dsp3 for the QCD vacuum is known as one of the most promising scenarios for quark confinement. For a review of the dual superconductivity picture, see, e.g., PhysRept . For this hypothesis to be realized, we must show the existence of some magnetic objects which can cause the dual Meissner effect. Then, the resulting chromofields are squeezed into the flux tube by the dual Meissner effect. This situation should be compared with the Abrikosov–Nielsen–Olesen (ANO) vortex NO1 ; NO2 in the $U(1)$ gauge-scalar model as a model describing the superconductor. In the context of the superconductor, in type II the repulsive force acts among the vortices, while in type I the attractive force acts. The boundary of the type I and type II is called the Bogomol’nyi–Prasad–Sommerfield (BPS) limit and no forces act among the vortices. From the viewpoint of the dual superconductivity picture, the type of dual superconductor characterizes the vacuum of the Yang–Mills theory or QCD for quark confinement.

The type of dual superconductor has been investigated for a long time by fitting the chromoelectric flux obtained by lattice simulations to the magnetic field of the ANO vortex. The preceding studies type2-1 ; type2-2 ; type2-3 ; type2-4 done in 1990’s concluded that the vacuum of the Yang–Mills theory is of type II or the border of type I and type II as a dual superconductor. In these studies, however, the fitting range was restricted to a long-distance region from a flux tube. The improved studies Koma1 ; Koma2 concluded that the vacuum of the Yang–Mills theory can be classified as weakly type I dual superconductor. Recent studies Cea-Cosmai1 ; Cea-Cosmai2 ; type1-1 ; type1-2 based on the standard framework of lattice gauge theory, and studies Kato-Kondo-Shibata ; Shibata-Kondo-Kato-Shinohara based on the new formulation KSSMKI08 ; SKS10 , on the other hand, show that the vacua of the $SU(2)$ and $SU(3)$ Yang–Mills theories are strictly type I dual superconductor. In these works Cea-Cosmai1 ; Cea-Cosmai2 ; type1-1 ; type1-2 ; Kato-Kondo-Shibata ; Shibata-Kondo-Kato-Shinohara , the Clem ansatz Clem was used to incorporate also the short distance behavior of a flux tube. The Clem ansatz assumes an analytical form for the behavior of the complex scalar field (as the order parameter of a condensation of the Cooper pairs), which means that it still uses an approximation. In this work, we shall fit the chromoelectric flux tube to the magnetic field of the ANO vortex in the $U(1)$ gauge-scalar model without any approximations to examine the type of dual superconductor. Indeed, we determine the Ginzburg–Landau (GL) parameter by fitting the lattice data of the chromoelectric flux to the numerical solution of the ANO vortex in the whole range. The resulting value of the GL parameter reconfirms that the dual superconductivity of $SU(2)$ Yang–Mills theory is of type I.

In addition, in order to estimate the force acting among the flux tubes, we investigate the Maxwell stress force carried by a single vortex configuration. Recently, the Maxwell stress force distribution around a quark-antiquark pair was directly measured on a lattice via the gradient flow method EMT . Our results should be compared with theirs. For this purpose, we shall calculate the energy-momentum tensor originating from a single ANO vortex solution to obtain the distribution of the Maxwell stress force corresponding to the obtained value of the GL parameter.

This paper is organized as follows. In Section $\rm{II}$, we introduce an operator to measure chromofields produced by a pair of quark and antiquark on a lattice. We review the results of lattice measurements in Kato-Kondo-Shibata . In Section $\rm{III}$, we give a brief review of the ANO vortex in the $U(1)$ gauge-scalar model. Then, we discuss the type of superconductor characterized by the GL parameter. In Section $\rm{IV}$, we explain a new method of fitting after giving a brief review of the fitting method based on the Clem ansatz adopted in the previous study Kato-Kondo-Shibata in order to compare our new result with the previous one. In Section V, we study the distribution of the force around a single flux tube by considering the Maxwell stress tensor. In Section $\rm{VI}$, we summarize our results. In Appendix A, we explain the advantage of the operator which we propose based on the new formulation to measure the gauge-invariant field strength on a lattice.

II Operator on a lattice to measure the flux tube

In order to measure the chromofield strength $\mathscr{F}_{\mu\nu}$ generated by a pair of static quark and antiquark belonging to the fundamental representation of the gauge group $G=SU(2)$, we use the gauge-invariant operator proposed by Di Giacomo, Maggiore, and Olejnik Giacomo using the Wilson loop operator $W[U]$ defined by the Yang–Mills link variable $U\in SU(2)$ along a path $C$ ($L\times T$ rectangular):111We use the notation $\rho[U]$ to indicate the average coming from the operators defined in terms of the original link variable $U$, since we define the similar operator defined in terms of the different variable later.

[TABLE]

where $U_{P}$ is a single plaquette constructed by $U$ and $L[U]$ is called the Schwinger line operator connecting the Wilson loop operator $W[U]$ and the plaquette $U_{P}$. See FIG.1 for the setup of the operator $W[U]L[U]U_{P}L^{\dagger}[U]$. In the continuum limit where the lattice spacing $\epsilon$ vanishes $\epsilon\to 0$, $\rho[U]$ reduces to

[TABLE]

where $\mathscr{A}\in\mathfrak{su}(2)$ stands for the gauge field of the continuum $SU(2)$ Yang–Mills theory, which is related to the link variable $U$ as $U_{x,\mu}=\exp\left(-ig\epsilon\mathscr{A}_{\mu}(x)\right)$. Thus, the field strength $F_{\mu\nu}[U]$ generated by a quark-antiquark pair can be obtained by

[TABLE]

FIG.2 shows the chromofield strength $F_{\mu\nu}[U]$ measured at the midpoint of the $q\bar{q}$ pair for the $8\times 8$ Wilson loop on the $24^{4}$ lattice at $\beta=2.5$ Kato-Kondo-Shibata . In this paper, we have fixed the physical scale of the lattice spacing $\epsilon=0.08320\ \mathrm{fm}$ at $\beta=2.5$ for $SU(2)$ by fixing the physical string tension $\sigma_{\rm phys}=(440\ \mathrm{MeV})^{2}$ according to the relation $\sigma_{\rm lat}=\sigma_{\rm phys}\epsilon^{2}$ Kato-Kondo-Shibata . Our results are consistent with the preceding studies Cea-Cosmai1 ; Cea-Cosmai2 .

In the previous study Kato-Kondo-Shibata , we used the new formulation KSSMKI08 ; SKS10 of the lattice Yang–Mills theory by decomposing the gauge field $U_{x,\mu}$ into $V_{x,\mu}$ and $X_{x,\mu}$, $U_{x,\mu}=X_{x,\mu}V_{x,\mu}$, where $V_{x,\mu}\in SU(2)$ called the restricted link variable is supposed to have the same transformation law as the original link variable $U$ under the gauge transformation, and a remaining part $X_{x,\mu}\in SU(2)$ called the remaining site variable transforms in an adjoint way under the gauge transformation. The restricted link variable $V_{x,\mu}$ plays a very important role for realizing the dual superconductor picture, since the dominant mode for quark confinement is extracted from it, for example, $V_{x,\mu}$ induces naturally the magnetic current. See, e.g., PhysRept for more details.

In the new formulation, the key ingredient is the color direction field $\bm{n}_{x}$ which takes the value in the $SU(2)$ Lie algebra with a constraint of a unit length $\bm{n}_{x}\cdot\bm{n}_{x}=1$ and transforms in an adjoint way under the gauge transformation. The color direction field $\bm{n}_{x}$ is in advance obtained as a functional of the original link variable $U_{x,\mu}$ by solving the reduction condition IKKMSS06 ; PhysRept . Then the restricted link variable $V_{x,\mu}$ is expressed in terms of the original link variable $U_{x,\mu}$ and the color direction field $\bm{n}_{x}$ as

[TABLE]

See Appendix A for more details.

In view of these, we propose to use the operator $\rho[V]$ similar to (1) by replacing the full link variable $U$ by the restricted link variable $V$:

[TABLE]

where $W[V]$ is the restricted Wilson loop operator obtained by replacing the link variable $U$ by $V$. In the continuum limit $\epsilon\to 0$, $\rho[V]$ reduces to

[TABLE]

and therefore, we can define the chromofield strength $F_{\mu\nu}[V]$ generated by $q\bar{q}$ pair for the restricted link variable $V$ by

[TABLE]

FIG.3 shows the restricted chromofield strength $F_{\mu\nu}[V]$ measured in the same settings as $F_{\mu\nu}[U]$ Kato-Kondo-Shibata .

In Appendix A, we demonstrate advantages of using $\rho[V]$ constructed from the restricted link variable $V$ based on the new formulation, in sharp contrast to the preceding operator $\rho[U]$ defined in terms of the original link variable $U$ based on the ordinary framework of lattice gauge theory: (i) The operator $\rho[V]$ enables us to extract the non-trivial gauge-invariant and Abelian-like field strength which is used to measure the chromoelectrix flux, in sharp contrast to the gauge-covariant non-Abelian field strength. (ii) The operator $\rho[V]$ does not depend on the choice of the Schwinger lines $L,L^{\dagger}$, namely, the shape of $L,L^{\dagger}$ and the position $z$ at which the Schwinger lines are inserted.

Recent study renormalize suggests that the operator $\rho[U]$ undergoes nontrivial renormalizations, which depend on the length and on the number of cusps in the Schwinger lines. The study Baker suggests that the extended smearing behaves like an effective renormalization of the operator $\rho[U]$. For $\rho[V]$, however, such renormalizations are not necessary since $\rho[V]$ does not depend on the Schwinger lines. On the other hand, renormalization or smearing for the restricted Wilson loop operator and the probe should be taken into account. In the previous study Kato-Kondo-Shibata , we used the hypercubic blocking (HYP) method HYP once for the link variables on the Wilson loop to reduce high-energy noises for both $U$ and $V$. However, we find numerically that for the restricted field $V$, the measured expectation value hardly differs from the unsmeared case.

First of all, we observe that the $z$-component of the restricted chromoelectric field $E_{z}[V]$ forms a uniform flux tube compared with a non-uniform one $E_{z}[U]$ Kato-Kondo-Shibata ; Shibata-Kondo-Kato-Shinohara , since the effect due to the static sources placed at a finite distance in $E_{z}[V]$ is smaller than $E_{z}[U]$. Therefore, the restricted chromoelectric flux $E_{z}[V]$ can be well approximated by the ANO vortex with an infinite length. Moreover, it was shown in the previous studies Kato-Kondo-Shibata ; Shibata-Kondo-Kato-Shinohara that the type of dual superconductor determined only by the flux tube does not change irrespective of whether we use $E_{z}[U]$ or $E_{z}[V]$. By these reasons, we shall use the data of $E_{z}[V]$ for fitting.

It should be noticed that we can define the magnetic current $k_{\mu}$ induced by the chromofield $F_{\mu\nu}[V]$ as

[TABLE]

with the lattice derivative $\nabla_{\nu}$ so that the conservation law $\nabla_{\mu}k_{\mu}=0$ holds Kato-Kondo-Shibata ; Shibata-Kondo-Kato-Shinohara . Since the nontrivial component of the chromofield $F_{\mu\nu}[V]$ is only the $z$-component $E_{z}[V]$ of the chromoelectric field (see the left panel of FIG.3), the induced magnetic current $k_{\mu}$ has only the component $k_{\varphi}$ circulating around a flux tube. The left panel of FIG.4 is an illustration of the relation between the chromoelectric field $\bm{E}$ and the induced magnetic current $k_{\mu}$. The right panel of FIG.4 is a plot of the chromoelectric field $E_{z}[V]$ and the magnetic current $k_{\varphi}$ induced around a single chromoelectric flux tube.

III The gauge-scalar model and type of superconductor

III.1 The Abrikosov–Nielsen–Olesen vortex

In this subsection, we give a brief review of the $U(1)$ gauge-scalar model with the Lagrangian density given by

[TABLE]

where $\lambda$ is a coupling constant of the scalar self-interaction, and $v$ is a value of the magnitude $|\phi(x)|$ of the complex scalar field $\phi(x)$ at the vacuum $|x|=\infty$. The asterisk $({}^{\ast})$ denotes the complex conjugation. The field strength $F_{\mu\nu}$ of the $U(1)$ gauge field $A_{\mu}$ and the covariant derivative $D_{\mu}\phi$ of the scalar field $\phi$ are defined by

[TABLE]

where $q$ is the electric charge of the scalar field $\phi(x)$. The Euler–Lagrange equations are given as

[TABLE]

where we define the electric current $j_{\mu}$ by

[TABLE]

In order to describe the vortex solution, we introduce the cylindrical coordinate system $(\rho,\varphi,z)$ for the spatial coordinates with the associated unit vectors $\bm{e}_{\rho},\bm{e}_{\varphi},$ and $\bm{e}_{z}$, and adopt a static and axisymmetric ansatz:

[TABLE]

where $n$ is an integer. Under this ansatz, the field equations (13) and (12) are cast into

[TABLE]

where a non-vanishing component $j_{\varphi}$ of the electric current is written as

[TABLE]

Moreover, the magnetic field $\bm{B}$ is given in the present ansatz by

[TABLE]

To determine the boundary conditions, let us consider the static energy $E$. The energy-momentum tensor $T^{\mu\nu}$ is obtained from the Lagrangian density (9) as

[TABLE]

Notice that this energy-momentum tensor is symmetric, i.e., $T^{\mu\nu}=T^{\nu\mu}$. Then, the static energy $E$ is obtained as

[TABLE]

In what follows, we consider the energy per unit length of a vortex to avoid the divergence, since the energy density $T^{00}$ does not depend on $z$.

The static energy $E$ given by (21) is nonnegative, $E\geq 0$. The equality $E=0$ holds if and only if

[TABLE]

are satisfied. Since the equation (22) is the solution of the field equations (16) and (17), we call it the vacuum solution.

Therefore, we require the solution to satisfy the boundary conditions for $\rho\to\infty$:

[TABLE]

so that the static energy $E$ does not diverge in the long-distance region $\rho\gg 1$. Indeed, these boundary conditions describe that in the long-distance region, the scalar field $\phi(x)$ goes to its vacuum value $|\phi(\infty)|=v$ and the gauge field $A_{\mu}(x)$ becomes the pure gauge configuration.

In the limit $\rho\to 0$, we assume

[TABLE]

so that the energy $E$ does not have a short-distance divergence.

Now we can clarify the meaning of the integer $n$ by using the boundary conditions. Let us consider the magnetic flux $\Phi$ passing through the surface $S$ bounded by a circle $C$ with the center at the origin and the large radius $\rho\to\infty$,

[TABLE]

which implies that the integer $n$ corresponds to the quantization of the magnetic flux. By this reason, we call the integer $n$ the topological charge, especially the winding number of a vortex.

Motivated by the vacuum solution (22), we modify the ansatz for the gauge field $A(\rho)$ as

[TABLE]

Moreover, in order to make the field equations dimensionless, we introduce the dimensionless variable:

[TABLE]

and redefine the profile functions as $f(\rho)=f(R)$ and $a(\rho)=a(R)$. Thus, the field equations (16), (17), and (18) are rewritten into

[TABLE]

where the prime $({}^{\prime})$ stands for the derivative with respect to $R$. The boundary conditions are also modified as

[TABLE]

We have simultaneously solved the field equations (28) and (29) in a numerical way under the boundary conditions (30) and (31). The left panel of FIG.5 shows the solutions $f$ and $a$ of the field equations (28) and (29) as functions of $R$ with a unit winding number $n=1$ for various values $\kappa$, $\kappa=\frac{1}{2},\frac{1}{\sqrt{2}}$, and $1$ of the Ginzburg–Landau (GL) parameter, which is defined by

[TABLE]

For the physical meaning of the GL parameter, see the next section. This solution is called the Abrikosov–Nielsen–Olesen (ANO) vortex NO1 ; NO2 . When we introduce the dimensionless magnetic field $b(R)$ and electric current $j(R)$, (18) and (19) are rewritten as follows:

[TABLE]

The right panel of FIG.5 shows the dimensionless magnetic field $b(R)$ corresponding to (34). Notice that the magnetic field $b(R)$ has no short-distance divergences, which is supported by the boundary condition (24). This means that the boundary condition (24) implies the regularity of the magnetic field $b(R)$ and the finiteness of the energy $E$ for a short distance.

III.2 Type of the superconductor

In order to investigate the asymptotic forms of the profile functions in the long-distance region $R\gg 1$, we introduce $g$ and $w$ in place of $f$ and $a$ as functions of $R$ by

[TABLE]

where $|g(R)|,|w(R)|\ll 1$ for $R\gg 1$. Then, the field equations for $g$ and $w$ read

[TABLE]

The second equation (37) can be solved by using the modified Bessel function of the second kind $K_{\nu}(x)$ as

[TABLE]

which behaves for $R\gg 1$ as

[TABLE]

Therefore, the magnetic field $B_{z}(R)$ has the asymptotic form for $R\gg 1$:

[TABLE]

where we have used the formula $zK_{\nu}^{\prime}(z)+\nu K_{\nu}(z)=-zK_{\nu-1}(z)$.

Inserting the asymptotic form (39) of $w(R)$ into the first equation (36), we have the closed equation for $g(R)$

[TABLE]

The solution of this inhomogeneous equation is given by

[TABLE]

where the first term is the general solution of the homogeneous equation obtained by ignoring the right hand side of (41) and the second term is a particular solution of (41). In terms of the dimensionful variable $\rho$, $g(R)$ behaves as

[TABLE]

which means that the fall-off factor of the scalar field must be distinguished by the value of $\lambda/q$.

We can define two typical lengths $\delta$ and $\xi$ by

[TABLE]

and the ratio by

[TABLE]

The length $\delta$ is called the penetration length (or depth), at which the magnitude of the magnetic field $B_{z}$ falls to $1/e\simeq 37\%$ of its original value at the origin $\rho=0$. The length $\xi$ is called the coherent length because the magnitude of the scalar field $|\phi(x)|$ grows to $1-1/e\simeq 63\%$ of its vacuum value $v$. (See FIG.6.) Taking into account the fall-off rates (or the masses) of the gauge and scalar fields, the mass of the gauge field $m_{V}=\sqrt{2}qv$ is larger than that of the scalar field $m_{S}=\sqrt{2}\lambda v$ for $\kappa<\frac{1}{\sqrt{2}}$, while for $\kappa>\frac{1}{\sqrt{2}}$ the opposite situation occurs. At the critical value $\kappa=\frac{1}{\sqrt{2}}$, the two masses $m_{V}$ and $m_{S}$ become equal: $m_{V}=m_{S}$. Therefore, the superconductor is classified by the value of the ratio $\kappa$ as

[TABLE]

The ratio $\kappa$ is called the Ginzburg–Landau (GL) parameter. The limit $\kappa\to\infty$, which is realized by $\xi\to 0$ or $m_{S}\to\infty$, is called the London limit.

IV Type of dual superconductor

To determine the type of dual superconductivity for $SU(2)$ Yang–Mills theory, we simultaneously fit the chromoelectric field and the induced magnetic current obtained by the lattice simulation Kato-Kondo-Shibata (see FIG.3 and FIG.4) to the magnetic field and electric current of the $n=1$ ANO vortex.

IV.1 The previous study using the Clem ansatz

In this subsection, we give a review of the approximated method of fitting with the Clem ansatz Clem . The previous studies Cea-Cosmai1 ; Cea-Cosmai2 ; type1-1 ; type1-2 ; Kato-Kondo-Shibata ; Shibata-Kondo-Kato-Shinohara considered only the regression of the chromoelectric flux, however in this paper, we also take into account the regression of the induced magnetic current to compare with our new method. In the Clem ansatz adopted to the $U(1)$ gauge-scalar model, the scalar profile function $f(\rho)$ is assumed to be

[TABLE]

where $\zeta$ is a variational parameter for the core radius of the ANO vortex and $\rho$ is a dimensionful variable $\rho=R/(qv)$. For the profile function of the gauge field $a(\rho)$, we introduce a new function $w(\rho)$ by

[TABLE]

which satisfies the boundary condition (30), i.e., $a(\rho=0)=0$. Then, the field equation (29) for the gauge field is now written as the differential equation for $w$:

[TABLE]

where we have defined a variable $x:=\sqrt{\rho^{2}+\zeta^{2}}$. The solution is given by the modified Bessel function of the second kind $K_{\nu}(z)$ as

[TABLE]

and hence

[TABLE]

Therefore, the magnetic field $B(\rho)$ is given by

[TABLE]

where we have defined

[TABLE]

with the external flux $\Phi=2\pi n/q$. The electric current $\bm{J}(\rho)=J(\rho)\bm{e}_{\varphi}$ is also written as:

[TABLE]

In the present setting, the energy per unit length $E$ can be calculated by restricting ourselves to the unit vortex with $n=1$ as

[TABLE]

where we have introduced the parameter $s=\sqrt{2}qv\zeta$. Since the vortex solution is obtained by minimizing the energy with respect to the parameter $s$, or $\zeta$, for a given GL parameter $\kappa$, the energy (57) must satisfy

[TABLE]

Therefore, the GL parameter $\kappa$ is given by

[TABLE]

In the previous study Kato-Kondo-Shibata , we adopted the fitting only for the chromoelectric flux. In this paper, we adopt the fitting for the chromoelectric flux and magnetic current simultaneously. In what follows, we use values measured in the lattice unit, e.g., the distance $\hat{y}=y/\epsilon$ with a lattice spacing $\epsilon$, the chromoelectric flux $E_{z}(\hat{y})=F_{34}[V](\hat{y})$ in (7), and the magnetic current $k_{\varphi}(\hat{y})$ in (8). Then, we denote the set of data as $(\hat{y}_{i},E_{z}(\hat{y}_{i}),\delta E_{z}(\hat{y}_{i}))$ for the chromoelectric field, and $(\hat{y}_{j},k_{\varphi}(\hat{y}_{j}),\delta k_{\varphi}(\hat{y}_{j}))$ for the induced magnetic current, where $\delta\mathcal{O}$ represents the error of the measurement $\mathcal{O}$.

To define dimensionless regression functions, let us rescale the parameters $\beta$ and $\zeta$ to be dimensionless by using the lattice spacing $\epsilon$ as

[TABLE]

and hence a parameter $\alpha$ is rescaled as

[TABLE]

We also rescale the magnetic field $B$ and the electric current $J$ as

[TABLE]

Then, we can define regression functions by

[TABLE]

with the dimensionless variable $\hat{\rho}:=\rho/\epsilon$ in the lattice unit. Then, the error functions of the regression with the weights are given by

[TABLE]

When we assume that these errors obey independent standard normal distributions, the parameters $\hat{\alpha},\hat{\beta}$, and $\hat{\zeta}$ can be estimated by maximizing the log-likelihood function $\ell(\hat{\alpha},\hat{\beta},\hat{\zeta})$ for (65) and (66) defined by

[TABLE]

The GL parameter $\kappa$ is determined according to (59) in terms of the estimated values $\hat{\beta}_{\star}$ and $\hat{\zeta}_{\star}$ by

[TABLE]

The obtained values in the previous work Kato-Kondo-Shibata , which can be achieved by ignoring the second term in (67) and restricting the fitting range to $2\leq\hat{\rho}\leq 8$, are given by

[TABLE]

where ${\rm MSR}_{\rm flux}$ is the sum of squared residuals for the regression of (65) divided by the degrees of freedom (d.o.f.) for fitting: (the number of data points) minus (the number of independent variational parameters), i.e., ${\rm d.o.f.}=7-3=4$.

By incorporating also the regression of the electric current $J$, the fitting result is in good agreement with (69):

[TABLE]

It should be noticed that the fitting range is restricted to $2\leq\hat{\rho}\leq 8$ as well as (69).

We further investigate the fitting by using the whole range $0\leq\hat{\rho}\leq 8$. The result obtained by using only the flux is given by

[TABLE]

By including the magnetic current, the result is given by

[TABLE]

The fitting result (72) is shown in the left panel of FIG.7. We find that the inclusion of the short range modifies the value of the GL parameter $\kappa$ to a smaller one. We also find that the inclusion of the regression for the magnetic current indeed improves the accuracy of the fitting in both cases of the fitting range $2\leq\hat{\rho}\leq 8$ and $0\leq\hat{\rho}\leq 8$.

IV.2 The new method

In this subsection, we shall fit the chromoelectric flux and the magnetic current to the magnetic field and the electric current of the ANO vortex simultaneously without any approximations. The advantage of the new method could be that the value of the GL parameter $\kappa$ is a direct fitting parameter unlike the case in the Clem ansatz.

Such a fitting can be done by using the regression functions $B$ and $J$ constructed by the solutions, $f(R)$ and $a(R)$, of the field equations (28) and (29) through the dimensionless magnetic field $b(R)$ in (34) and the electric current $j(R)$ in (33). However, there are difficulties to estimate the model parameters, when we flow the same procedure as in the previous subsection. When we construct the regression functions $B$ and $J$ from the numerical solutions $f(R)$ and $a(R)$ by solving the field equations (28) and (29), we also calculate the regression functions numerically. Indeed, it is necessary to numerically calculate the derivative in (34) separately, and this causes a large numerical error even if one obtains the solutions $f(R)$ and $a(R)$ with small errors. To avoid these difficulties, we reorganize the field equations to include both $b(R)$ and $j(R)$ as independent unknown functions by

[TABLE]

where we have decomposed the second order differential equation (29) for the gauge profile function $a(R)$ into two independent first order differential equations (74) and (75) and one algebraic equation (76). We solve these coupled equations simultaneously. We impose the following boundary conditions for four unknown functions $f(R),a(R),b(R)$, and $j(R)$:

[TABLE]

From (33) and (34) , we obtain the regression functions with the dimensionless variational parameters $\hat{\eta}:=qv^{2}\epsilon^{2}$ $\hat{\tau}:=q^{2}v^{3}\epsilon^{3}$ in the lattice unit by

[TABLE]

where $\hat{\rho}:=\rho/\epsilon$ is a dimensionless variable, and $\kappa$ is the GL parameter.

By numerically solving (73)–(76) simultaneously and maximizing the log-likelihood function (67) with the regression functions (79) by varying the parameters $\hat{\eta},\hat{\tau}$, and $\kappa$, we estimate the model parameters $\hat{\eta},\hat{\tau}$, and $\kappa$. Note that since the coupled differential equations (73)–(76) with respect to $R=\hat{\tau}\hat{\rho}$ depends on only the GL parameter $\kappa$, the variation of the parameters $\hat{\tau}$ and $\hat{\eta}$ does not deform the functions $b(R)$ and $j(R)$. Thus we obtain the results:

[TABLE]

The fitting result is shown in the right panel of FIG.7.

We further obtain the penetration length $\delta$ and coherent length $\xi$ defined in (46) by using the fitted values (80) as

[TABLE]

FIG.8 shows the penetration and coherent lengths for the fitted value of the GL parameter $\kappa$ with corresponding functions $|\phi|/v$ and $B_{z}/(qv)^{2}$. See also FIG.6.

This new result shows that the vacuum of $SU(2)$ Yang–Mills theory is of type I, $\kappa=0.565\pm 0.053<1/\sqrt{2}\approx 0.707$, which is consistent with the results based on the Clem ansatz (69) and (70) within errors. We find that the inclusion of the regression for the magnetic current (70), (72), and (80) give small errors of the GL parameter $\kappa$ than the excluded ones (69) and (71). We also observe that the sums of squared residuals for both the flux and current in the new method become smaller than the fitting method based on the Clem ansatz. Therefore, the inclusion of the fitting for the magnetic current is important to improve the accuracy.

We should be aware of the effect of changing the fitting range. If we choose the fitting range $2\leq\hat{\rho}\leq 8$, the value of the GL parameter $\kappa$ agrees with the result of $0\leq\hat{\rho}\leq 8$ within their errorbars. This fact quite differs from the previous result based on the Clem ansatz represented in (69)–(72).

V Distribution of the stress force around a vortex

In what follows, to clarify the difference between type I and II of dual superconductors in view of force among the chromoelectric fluxes, we investigate the Maxwell stress tensor according to the proposal EMT2-1 ; EMT2-2 ; EMT2-3 . We find that the components (20) of the energy-momentum-stress tensor $T^{\mu\nu}$ around an ANO vortex are written under the ansatz (15), (26) and (34) as

[TABLE]

and all the off-diagonal components vanish.222 Here, we change the sign of $T^{jk}$ defined in (20) by using the ambiguity of the overall sign of the Noether current in order to reproduce the conventional Maxwell stress tensor.

FIG.9 shows $T^{\rho\rho},$ and $T^{\varphi\varphi}$, and $T^{zz}$ for various GL parameter $\kappa=\frac{1}{5},0.565,\frac{1}{\sqrt{2}},1$, and $\infty$ with a unit winding number.

One finds that $T^{\rho\rho}$ is always positive in type I, while always negative in type II:

[TABLE]

[TABLE]

At the boundary between them, i.e., the BPS limit $\kappa=\frac{1}{\sqrt{2}}$, $T^{\rho\rho}$ is identically zero:

[TABLE]

It should be noticed that the components $T^{\rho\rho}$ and $T^{\varphi\varphi}$ are not independent, since the conservation law of the Noether current $\partial^{\mu}T_{\mu\nu}=0$ leads to

[TABLE]

This yields that the sign of $T^{\varphi\varphi}(R)$ flips and hence there is a critical value $R=R_{\ast}$ where $T^{\varphi\varphi}(R_{\ast})=0$. See the middle panel of FIG.9.

Next, we investigate the force acting on the area element of the flux tube. By using the Maxwell stress tensor, the stress force $\bm{F}$ acting on the infinitesimal area element $d\bm{S}$ is given by

[TABLE]

where $\bm{n}$ is a normal vector perpendicular to the area element $dS$, and $\Delta S$ stands for the area of $dS$. See FIG.10. The left and mid panels show the situations for the ANO vortex, while the right panel shows the corresponding situation in the electromagnetism, where a pair of electric charges $\pm q$ is located at $\mp\infty$ on the $z$-axis.

If we choose $\bm{n}$ to be equal to the normal vector pointing the $\rho$-direction, i.e., $\bm{n}=\bm{e}_{\rho}$, the corresponding stress force $\bm{F}^{(\rho)}$ reads

[TABLE]

Since $T^{\rho\rho}$ obeys (86) and (87), we observe that $\bm{F}^{(\rho)}\cdot\bm{e}_{\rho}=T^{\rho\rho}\Delta S$ is always positive in type I, while always negative in type II. Therefore, we find that $\bm{F}^{(\rho)}$ represents the attractive force for type I, while the repulsive force for type II.

If we choose $\bm{n}$ as the unit vector for the $\varphi$-direction, $\bm{n}=\bm{e}_{\varphi}$, the corresponding stress force $\bm{F}^{(\varphi)}$ is written as

[TABLE]

The sign of $\bm{F}^{(\varphi)}\cdot\bm{e}_{\varphi}=T^{\varphi\varphi}\Delta S$ changes, since the sign of $T^{\varphi\varphi}$ flips at some critical value $R=R_{\ast}$. This feature could be an artifact due to the infinite length of the ANO vortex and should be investigated in a more realistic situation.

The other choice of $\bm{n}$ is to be parallel to the ANO vortex, i.e., $\bm{n}=\bm{e}_{z}$. The corresponding stress force $\bm{F}^{(z)}$ can be written as

[TABLE]

FIG.11 shows the distribution of the stress forces $\bm{F}^{(\rho)}$ and $\bm{F}^{(z)}$ in $y-z$ plane. Therefore, $\bm{F}^{(z)}$ represents the attractive force. Since $T^{zz}$ is always positive $T^{zz}>0$ due to (83), $\bm{F}^{(z)}$ points the same direction regardless of the value of the GL parameter $\kappa$.

It should be noted that the situation of the type II superconductor is similar to the electromagnetism, see the mid and right panels of FIG.10.

Using the parameters obtained by fitting to the ANO vortex, we can reproduce the distribution of the Maxwell stress force around the flux tube, which is shown in FIG.12. This result indeed supports the type I dual superconductor for quark confinement.

Our analysis on the Maxwell stress tensor around an ANO vortex agrees with the result obtained by the preceding work EMT ; EMT2-1 ; EMT2-2 ; EMT2-3 .

VI Conclusion

In this paper, we have studied the type of dual superconductivity for the $SU(2)$ Yang–Mills theory based on the new method:

• •

We have introduced the restricted field $V$ to extract the dominant mode for quark confinement and define the induced magnetic current in a gauge-invariant way.

• •

We have solved the field equations of the ANO vortex in the $U(1)$ gauge-scalar model numerically without any approximations. The previous method is based on the Clem ansatz which assumes an analytic form of the complex scalar field without solving the field equations.

• •

We have used the resulting magnetic field and the electric current to fit respectively the chromoelectric flux tube and the induced magnetic current obtained by lattice simulations. In the previous method, only the regression for the chromoelectric flux tube was considered.

We have reconfirmed that the vacuum of the $SU(2)$ Yang–Mills theory is of type I as a dual superconductor with the GL parameter $\kappa=0.565\pm 0.053$ by using the new method.

This result obtained by the new method should be compared with the previous method based on the Clem ansatz:

• •

We found that the result of type I agrees with Kato-Kondo-Shibata reproduced and supplemented by (69)–(72). In the new method, we determined the GL parameter with good accuracy.

• •

We have investigated the sensitivity for the fitting range. We found that the inclusion of the short range modifies the value of the GL parameter $\kappa$ to a smaller one under the Clem ansatz. For the new method, on the other hand, we found that the inclusion or exclusion of the short range does not effect the GL parameter $\kappa$.

• •

We also found that the new method proposed in this paper improves the accuracy of the fitting as seen from the error of the GL parameter, or the mean of squared residuals in both methods. Therefore, the inclusion of the regression of the magnetic current is important.

Moreover, we have obtained the distribution of the Maxwell stress force around the flux tube by using the obtained GL parameter. It was observed that there exists the attractive force among the chromoelectric flux tubes, which also supports the type I dual superconductor.

Acknowledgement

The authors would like to thank Hideo Suganuma for valuable discussions, especially suggestions on error estimations. They would like to express sincere thanks to Ryosuke Yanagihara, Takumi Iritani, Masakiyo Kitazawa, and Tetsuo Hatsuda for very helpful and illuminating discussions on the Maxwell stress tensor in the early stage of their investigations, on which a part of the result presented in section V is based. This work was supported by Grant-in-Aid for Scientific Research, JSPS KAKENHI Grant Number (C) No.15K05042 and No.19K03840. S.N. thanks Nakamura Sekizen-kai for a scholarship.

Appendix A Restricted field variable in the new formulation

In this paper, we have used the new formulation KSSMKI08 ; SKS10 of the lattice Yang–Mills theory: the gauge field variable $U_{x,\mu}$ is decomposed into $V_{x,\mu}$ and $X_{x,\mu}$, $U_{x,\mu}=X_{x,\mu}V_{x,\mu}$. Here $V_{x,\mu}\in SU(2)$ is called the restricted link variable which has the same transformation law as the original link variable $U$ under the gauge transformation, while $X_{x,\mu}\in SU(2)$ is a remaining part called the remaining site variable which transforms in an adjoint way under the gauge transformation.

In this decomposition we introduce the so-called color direction field $\bm{n}_{x}$ which is a Lie-algebra valued field written using the Pauli matrices $\sigma_{A}$ ($A=1,2,3$) as

[TABLE]

and is subject to the condition of unit length

[TABLE]

Therefore, the color direction field $\bm{n}_{x}$ does not cover the whole $SU(2)$ and takes the value only in the coset $G/H=SU(2)/U(1)$. An expression exhibiting manifestly this fact is given using $\Omega_{x}\in SU(2)$ as

[TABLE]

which is supposed to transform in an adjoint way under the gauge transformation. The color direction field $\bm{n}_{x}$ is obtained by solving the reduction condition for a given set of the original link variables $U_{x,\mu}$. Therefore, the color direction field $\bm{n}_{x}$ is understood as a functional of the original link variables $U_{x,\mu}$. See, e.g., PhysRept for a review.

First, the expression (4) for the restricted link variable $V_{x,\mu}$ is obtained IKKMSS06 by solving the defining equations:

[TABLE]

By solving the defining equations (97) and (98), indeed, the link variable $V_{x,\mu}$ is determined up to an overall normalization constant in terms of the site variable $\bm{n}_{x}$ and the original link variable $U_{x,\mu}$ as $\tilde{V}_{x,\mu}$:

[TABLE]

However, the defining equation (97) is linear in $V_{x,\mu}$. Therefore, the normalization of $V_{x,\mu}$ cannot be determined by the defining equation alone. Consequently, unitarity is not guaranteed for the general solution of the defining equation and hence a unitarity condition must be imposed afterwards. Fortunately, this issue is easily solved at least for the $SU(2)$ group, since the speciality condition $\det V_{x,\mu}=1$ determines the normalization. Then the special unitary link variable $V_{x,\mu}$ is obtained after the normalization of $\tilde{V}_{x,\mu}$ as

[TABLE]

It is also shown IKKMSS06 that the naive continuum limit $\epsilon\rightarrow 0$ of the link variable $V_{x,\mu}=\exp(-i\epsilon g\mathscr{V}_{\mu}(x))$ reduces to the continuum expression:

[TABLE]

which agrees with the expression of the restricted field in the Cho-Duan-Ge decomposition in the continuum Cho80-1 ; Cho80-2 ; DG79 . This is indeed the case for the remaining variable $X_{x,\mu}=\exp(-i\epsilon g\mathscr{X}_{\mu}(x))$. The $SU(N)$ group is treated in KSSMKI08 ; SKS10 ; PhysRept .

Second, we focus on the fact that the color direction field $\bm{n}_{x}$ is covariantly constant under the restricted link variable $V_{x,\mu}$ by construction (97). The meaning of “covariantly constant” is that one can perform the parallel transport of a vector along a path from point $x$ to point $y$, so that the result is independent of the path chosen. In particular, parallel transport along a closed loop should leave the vector unchanged. Suppose we parallel transport $\bm{n}_{x}$ from a point $x$ to a point $x+\mu+\nu$ via two different paths: (I) from $x$ to $x+\mu$ to $x+\mu+\nu$, and (II) from $x$ to $x+\nu$ to $x+\mu+\nu$. We quickly see from (97) that

[TABLE]

Equating the right-hand side of the first and second lines, and using the unitarity of the $V$’s, we find that the color direction field $\bm{n}_{x}$ at site $x$ satisfies the relation

[TABLE]

where $V_{P_{x}}$ is the plaquette variable at $x$, namely, the product of link variables starting at $x$ along the plaquette $P$. The relation (103) is equivalent to

[TABLE]

due to the unitarity of $V_{P_{x}}$, $V_{P_{x}}V_{P_{x}}^{\dagger}=1$. This consideration can be generalized: Let $C_{x}$ be any contour on the lattice beginning and ending at site $x$, and let $V_{C_{x}}$ be the holonomy equal to the product of $V$ link variables around the loop $C_{x}$. It is not hard to see that (97) implies

[TABLE]

This equation holds at every site $x$ for every possible contour $C_{x}$. A naive inspection of these equations would yield the view that the equation (105) following from (103) or (104) is possible only when the $V$ field is a pure gauge with vanishing field strength. However, this is not true, because an $SU(2)$ group element $V\in SU(2)$ obeys the identity,

[TABLE]

See e.g., SKS10 for a proof of this identity. Indeed, this form of $V$, a linear combination of a unit matrix $\bm{1}$ and the color direction field $\bm{n}$ is consistent with (104) and (103).

Recall that we have already given the explicit expression for the restricted link variable $V_{x,\mu}$ as (4):

[TABLE]

This result shows that the product $\prod_{<x,\mu>\in C}V_{x,\mu}$ of $V_{x,\mu}$ along the loop $C$ does not agree with the pure gauge. However, it should be remarked that this conclusion is due to the special property of the color direction field $\bm{n}_{x}$ which we have adopted in this construction as mentioned above: The color direction field $\bm{n}_{x}$ takes the value only in the coset $G/H=SU(2)/U(1)$,

[TABLE]

In contrast, if we start from a field $\hat{\phi}_{x}$ taking the value in the whole $SU(2)$,

[TABLE]

and impose the covariant-constant condition

[TABLE]

then the restricted variable just agrees with the pure gauge

[TABLE]

This case was indeed shown in the continuum formulation Kondo18 .

The above remark is also understood in the continuum version of the new formulation PhysRept which agrees with the naive continuum limit obtained by taking the limit of vanishing lattice spacing $\epsilon\to 0$ in the lattice version. The restricted plaquette variable is expanded as

[TABLE]

The continuum field strength $\mathscr{F}_{\mu\nu}[\mathscr{V}](x)$ of the restricted field $V_{\mu}(x)$ is proportional to the color direction field $\bm{n}(x)$ by construction:

[TABLE]

This implies that the gauge transformation $\mathscr{F}_{\mu\nu}[\mathscr{V}](x)\to\Omega(x)\mathscr{F}_{\mu\nu}[\mathscr{V}](x)\Omega(x)^{\dagger}$ of the field strength $\mathscr{F}_{\mu\nu}[\mathscr{V}](x)$ is carried by the color direction field $\bm{n}(x)$ which transforms in the adjoint way $\bm{n}(x)\to\Omega(x)\bm{n}(x)\Omega(x)^{\dagger}$ under the gauge transformation $\Omega(x)\in SU(2)$ so that the field strength $f_{\mu\nu}(x)$ defined using the color direction field $\bm{n}(x)$ with a unit length $\bm{n}(x)\cdot\bm{n}(x)=1$ by

[TABLE]

is invariant under the gauge transformation. Note that (112) and (113) are consistent with (106), since in the continuum limit, $2\mathrm{tr}\left(\bm{n}_{x}V_{P}\right)$ reduces to the gauge-invariant field strength $\bm{n}(x)\cdot\mathscr{F}_{\mu\nu}[\mathscr{V}](x)$

[TABLE]

Therefore, (104) does not mean that the $V$ field is a pure gauge with vanishing field strength $\mathscr{F}_{\mu\nu}[\mathscr{V}]$. Thus we have the non-trivial gauge-invariant field strength $f_{\mu\nu}(x)$ which is used to measure the chromoelectrix flux. This feature is more clearly seen using the non-Abelian Stokes theorem for the Wilson loop operator, see e.g., PhysRept .

Finally, in oder to see the physical meaning of the restricted operator $\rho[V]$ and the independence from the Schwinger lines to be inserted, we cast it into another form. The covariant constantness (97) of the color direction field $\bm{n}_{x}$ under the restricted link variable $V_{x,\mu}$ yields another expression for $\rho[V]$. In fact, we apply (106) to $V_{P}$ in (5) and use the property (97) to parallel transport the color direction field $\bm{n}_{x}$ on the plaquette $P$ to $\bm{n}_{z}$ at the point $z$ on the line connecting a pair of quark and antiquark, see Fig. 1:

[TABLE]

where we have used the unitarity of $L_{V}$, $L_{V}L_{V}^{\dagger}=1$ in the last step.

Thus we find that $\rho[V]$ does not depend on both the choice of the Schwinger lines $L,L^{\dagger}$ and the position $z$ at which the color direction field is inserted. This is not the case for $\rho[U]$ constructed from the original gauge variable $U$. These facts demonstrate advantages of using quantities like $\rho[V]$ constructed from the restricted variable based on the new formulation.