Local correlation among the chiral condensate, monopoles, and color magnetic fields in Abelian projected QCD

Using the lattice gauge field theory, we study the relation among the local chiral condensate, monopoles, and color magnetic fields in quantum chromodynamics (QCD). First, we investigate idealized Abelian gauge systems of 1) a static monopole-antimonopole pair and 2) a magnetic flux without monopoles, on a four-dimensional Euclidean lattice. In these systems, we calculate the local chiral condensate on quasi-massless fermions coupled to the Abelian gauge field, and find that the chiral condensate is localized in the vicinity of the magnetic field. Second, using SU(3) lattice QCD Monte Carlo calculations, we investigate Abelian projected QCD in the maximally Abelian gauge, and find clear correlation of distribution similarity among the local chiral condensate, monopoles, and color magnetic fields in the Abelianized gauge configuration. As a statistical indicator, we measure the correlation coefficient $r$, and find a strong positive correlation of $r \simeq 0.8$ between the local chiral condensate and an Euclidean color-magnetic quantity ${\cal F}$ in Abelian projected QCD. The correlation is also investigated for the deconfined phase in thermal QCD. As an interesting conjecture, like magnetic catalysis, the chiral condensate is locally enhanced by the strong color-magnetic field around the monopoles in QCD.


I. INTRODUCTION
Quantum chromodynamics (QCD) is an SU(N c ) gauge theory to describe the strong interaction, and has presented many interesting subjects full of variety and difficult problems in physics. Actually, in spite of the simple form of the QCD action, this miracle theory creates hundreds of hadrons and leads to various interesting nonperturbative phenomena, such as color confinement and dynamical chiral-symmetry breaking [1]. This magic is due to the strong coupling of QCD in the low-energy region, and this strong-coupling nature drastically changes the vacuum structure itself. Therefore, perturbative technique is no more workable and analytical treatment of QCD is fairly difficult in the strong-coupling region. As a reliable standard technique, lattice QCD Monte Carlo simulations have been applied to analyze nonperturbative QCD [2,3].
Among the nonperturbative properties of QCD, spontaneous chiral-symmetry breaking is particularly important in our real world. Indeed, chiral symmetry breaking drastically influences the vacuum structure and gives a nontrivial vacuum expectation value of the chiral condensate qq , which plays the role of an order parameter. Also, it is considered that chiral symmetry breaking leads to dynamical quark-mass generation [1,4], and creates most of the matter mass of our Universe, apart from the dark matter, because only small masses of u, d current quarks and electrons are Higgs-origin in atoms [5] and their contribution to the nucleon mass is estimated to be small [6]. In addition, chiral symmetry breaking inevitably accompanies light pions of the Nambu-Goldstone bosons, and their small mass gives range of the nuclear force.
In nonperturbative QCD, color confinement is also one of the most important phenomena in physics, and presents an extremely difficult mathematical problem. Experiments for hadron spectra and lattice QCD studies for various interquark potentials [7][8][9][10] show that the quark confining force is basically characterized by a universal physical quantity of the string tension σ 0.89 GeV/fm. This universal string tension is physically explained by one-dimensional squeezing of the color electric flux, i.e., the color flux-tube formation in hadrons, as is also indicated by lattice QCD for both mesons [3] and baryons [11]. As for the relation between color confinement and chiral symmetry breaking, it is not yet clarified directly from QCD. While almost coincidence between deconfinement and chiral-restoration temperatures [12] suggests their close correlation, a lattice QCD analysis using the Dirac-mode expansion based on the Banks-Casher relation [13] indicates some independence of these phenomena in QCD [14,15].
For the quark confinement mechanism, Nambu [16], 't Hooft [17], and Mandelstam [18] proposed the dual superconductivity scenario, paying attention to analogy with the Abrikosov vortex in the superconductivity, where Cooper-pair condensation leads to the Meissner effect, and the magnetic flux is excluded or squeezed like a one-dimensional tube as the Abrikosov vortex. If the QCD vacuum can be regarded as the dual version of the superconductor, the electric-type color flux is squeezed between (anti)quarks in hadrons, and quark confinement can be physically explained by the dual Meissner effect. Because of the electromagnetic duality, the dual Meissner effect inevitably needs condensation of magnetic objects, i.e., color magnetic monopoles, which correspond to the dual version of the electric Cooper-pair bosonic field.
In the dual-superconductor picture for the QCD vacuum, however, there are two large gaps with QCD.
1. While QCD is a non-Abelian gauge theory, the dual-superconductor picture is based on an Abelian gauge theory subject to the Maxwell-type equations including magnetic currents, where electromagnetic duality is manifest.
2. While QCD includes only color electric variables, i.e., quarks and gluons, as the elementary degrees of freedom, the dual-superconductor picture requires condensation of color magnetic monopoles as a key concept.
Historically, to bridge between QCD and the dual-superconductivity, 't Hooft proposed Abelian gauge fixing [19], partial gauge fixing which only remains Abelian gauge degrees of freedom in QCD. By Abelian gauge fixing, QCD reduces into an Abelian gauge theory, where off-diagonal gluons behave as charged matter fields similar to W ± µ in the Weinberg-Salam model and give the color electric current j µ in terms of the residual Abelian gauge symmetry. As a remarkable fact in the Abelian gauge, color-magnetic monopoles appear as topological objects corresponding to the nontrivial homotopy group Π 2 (SU(N c )/U(1) Nc−1 ) = Z Nc−1 ∞ in a similar manner to appearance of 't Hooft-Polyakov monopoles [20] in the SU(2) non-Abelian Higgs theory. Thus, in the Abelian gauge, QCD is reduced into an Abelian gauge theory including both electric current j µ and magnetic current k µ , which is expected to give a theoretical basis of the dual-superconductor picture for the confinement mechanism, although off-diagonal gluons remain as charged matter fields.
From the viewpoint of Abelianzation of QCD, the maximally Abelian (MA) gauge [21] is an interesting special Abelian gauge. In the MA gauge, off-diagonal gluons have a large effective mass of about 1 GeV in both SU (2) and SU(3) lattice QCD [22][23][24], so that off-diagonal gluons become infrared inactive, and only the Abelian gluon is relevant at distances larger than about 0.2 fm. Also, monopole condensation is suggested from appearance of long entangled monopole worldlines [21,25] and the magnetic screening in lattice QCD [26,27].
In this way, by taking the MA gauge, the QCD vacuum can be regarded as an Abelian dual superconductor at a large scale, and color magnetic monopoles seem to capture essence of nonperturbative QCD. Note however that, even without gauge fixing, there is an evidence of monopole condensation in nonabelian gauge theories [27], and therefore it might be possible to define infrared-relevant monopoles in QCD and to construct the dual superconductor system in more general manner. In fact, MA gauge fixing gives a concrete way to extract infrared-relevant Abelian gauge manifold and monopoles from QCD.
In the context of the dual superconductir picture, close correlation between monopoles and chiral symmetry breaking was pointed out in the dual Ginzburg-Landau theory [28], in SU(2) lattice QCD in the MA gauge [29,30], and in SU(3) lattice QCD [31,32]. Since most of the pioneering lattice studies were done in SU(2) QCD or on a small lattice as 8 3 × 4, we recently investigated SU(3) QCD with a large volume, and find a clear correlation between monopoles and the chiral condensate in SU(3) lattice QCD in the MA gauge [33].
In this paper, as a continuation of Ref. [33], we proceed the lattice works for the relation between chiral symmetry breaking and color magnetic objects including monopoles. In particular, as a new point of this paper, we quantitatively study correlation of the local chiral condensate with color magnetic fields using the lattice gauge theory.
The organization of this paper is as follows. In Section 2, we review the MA gauge and Abelianization of QCD in SU(3) lattice formalism. In Section 3, we prepare magnetic objects in Abelian projected QCD. In Section 4, we consider the local chiral condensate and chiral symmetry breaking in Abelian gauge systems. In Section 5, we present idealized Abelian gauge systems of a static monopole-antimonopole pair on a lattice, and investigate the relation of the local chiral condensate with the magnetic objects. In Section 6, we perform SU(3) lattice QCD Monte Carlo calculations and study the relation among monopoles, magnetic fields, and the local chiral condensate in Abelian projected QCD in the MA gauge. Section 7 is devoted for summary and conclusion.

II. MAXIMALLY ABELIAN GAUGE AND ABELIANIZATION OF QCD
To begin with, we briefly review the lattice formalism for maximally Abelian (MA) gauge fixing and Abelianization in QCD.
Continuum QCD is described with the quark field q(x), the gluon field A µ (x) ∈ su(N c ) and the QCD gauge coupling g. In SU(N c ) lattice QCD [3], the gluon field is described as the SU(N c ) link variable U µ (s) ≡ exp (iagA µ (s)) ∈ SU(N c ) on four-dimensional Euclidean lattices with the spacing a and the volume V = L x L y L z L t .
Using the Cartan subalgebra H ≡ (T 3 , T 8 ) in SU(3), MA gauge fixing is defined so as to maximize by the SU(3) gauge transformation, and therefore this gauge fixing strongly suppresses all the off-diagonal fluctuation of the SU(3) gauge field. In the MA gauge, the SU(3) gauge group is partially fixed remaining its maximal torus subgroup U(1) 3 × U(1) 8 with the global Weyl (color permutation) symmetry [34], and QCD is reduced to an Abelian gauge theory. From the SU(3) link variable U MA µ (s) ∈ SU(3) in the MA gauge, we extract the Abelian link variable by maximizing the overlap Note that the distance between u µ (s) and U MA µ (s) becomes the smallest in the SU(3) manifold, and there is a constraint 3 i=1 θ i µ (s) = 0 (mod 2π) reflecting the uni-determinant of u µ (s). Here, θ i µ (s) (i = 1, 2, 3) is taken to be the principal value of −π ≤ θ i µ (s) < π. The Abelian projection is defined by the simple replacement of the SU(3) link variable U µ (s) by the Abelian link variable u µ (s) for each gauge configuration, that is, O[u µ (s)] is called "Abelian dominance" for the operator O [35].

III. MAGNETIC OBJECTS IN ABELIAN PROJECTED QCD
In this section, we prepare magnetic objects in Abelian projected QCD in four-dimensional Euclidean space-time.

A. Monopoles in Abelian projected QCD
In this subsection, we define monopoles in Abelian projected QCD in lattice formalism [40]. Like the ordinary SU(3) plaquette, the Abelian plaquette variable is defined as whereμ is the µ-directed unit vector in the lattice unit. The Abelian field strength θ i µν (s) (i = 1, 2, 3) is the principal value of the exponent in u µν (s), and is defined as with the forward derivative ∂ µ . Note that θ i µν (s) is U(1) 2 gauge invariant and corresponds to the regular Abelian field strength in the continuum limit of a → 0, while n i µν (s) corresponds to the singular gauge-variant Dirac string [40]. The electric current j i µ and the monopole current k i µ are defined from the Abelian field strength θ i µν as where ∂ µ is the backward derivative andθ µν is the dual tensor ofθ µν ≡ 1 2 µναβ θ αβ . Both electric and monopole currents are U(1) 2 gauge invariant, according to U(1) 2 gauge invariance of θ i µν (s). In the lattice formalism, k i µ (s) is located at the dual lattice L 4 dual of s α + 1 2 with flowing in µ direction [41]. In this way, Abelian projected QCD includes both electric current j i µ and monopole current k i µ . Remarkably, lattice QCD shows monopole dominance, i.e., dominant role of monopoles for quark confinement in the MA gauge [42]. Also, lattice QCD shows monopole dominance for chiral symmetry breaking, that is, monopoles in the MA gauge crucially contribute to spontaneous chiral-symmetry breaking in both SU (2) [29,31] and SU(3) lattice QCD [33].
In the lattice formalism, the monopole current k i µ appears on the dual lattice L 4 dual of s α + 1 2 , and therefore we define the local monopole density where P (s) denotes the dual lattices in the vicinity of s, i.e., P (s) = s ∈ L 4 dual |s − s | = 1 . Note here that the distance between the site s and its closest dual site s is |s − s | = 4 1 ( 1 2 ) 2 = 1 in the four-dimensional Euclidean space-time.

B. General argument for magnetic instability and magnetic objects in the QCD vacuum
In QCD in the MA gauge, color magnetic monopoles generally appear, and play an important role in nonperturbative properties, which migth looks curious, since the original QCD action does not have monopoles.
However, some active roles of magnetic objects would be natural in QCD, because QCD itself has color magnetic instability, and spontaneous generation of color magnetic fields generally takes place, as Savvidy first pointed out in 1977 [43,44].
In fact, in the QCD vacuum in the Minkowski space-time, the gluon condensate G a µν G a µν takes a large positive value, which physically means that the QCD vacuum is filled with color magnetic fields. Since the gluon condensate is expressed with color magnetic fields H a and color electric fields E a as its large positivity means inevitable significant generation of color magnetic fields. Thus, some superior role of magnetic objects is expected instead of electric objects in the Minkowski QCD vacuum.
In the Euclidean space-time, because of the space-time SO(4) symmetry, the roles of magnetic and electric fields become similar. Actually, the gluon condensate is written as G a µν G a µν = 2( H 2 a + E 2 a ), where the electromagnetic duality is manifest. Then, in Euclidean QCD, the electric field often behaves as a magnetic field, and therefore we regard the Euclidean electric field as a sort of the magnetic field in this paper.

C. Lorentz invariant quantities in Abelian projected QCD
In Abelian projected QCD, there are two Lorentz invariant quantities F and G in the Euclidean space-time: with the color magnetic field ( H i ) j ≡ 1 2 jkl F kl i and the color electric field ( E i ) j ≡ F j4 i . These quantities are also invariant under the residual U(1) 2 gauge transformation and global Weyl transformation [34], i.e., permutation of the color index, in the MA gauge.
Here, F is parity-even and expresses total magnitude of magnetic fields in the Euclidean space-time, since the electric field behaves as a magnetic field there. In this paper, we simply call F "magnetic quantity" in Euclidean gauge theories. Note that G is parity-odd and is just the Abelian projected quantity of the topological charge density on instantons in QCD, which might relate to chiral symmetry breaking.
In the lattice formalism, the field strength tensor is a plaquette variable spanning at s, s +μ, s +ν, and s +μ +ν, so that we define the Abelian field strength F i µν (s) as the local average of clover-type four plaquettes, and consider F and G as local quantities in each Abelian gauge configuration.

IV. LOCAL CHIRAL CONDENSATE AND CHIRAL SYMMETRY BREAKING IN GAUGE THEORIES
In this section, we consider the chiral condensate and chiral symmetry breaking in the gauge theory in terms of the quark propagator.

A. Local chiral condensate in lattice QCD
In this subsection, we briefly review the local chiral condensate in lattice formalism. The local chiral condensate can be calculated with the quark propagator for each gauge configuration U = {U µ (s)} generated with the Monte Carlo method.
As the lattice fermion, we here adopt the Kogut-Susskind (KS) fermion [3]. For the KS fermion, the Dirac operator with , and the KS Dirac eigenvalue equation takes the form of Here, the quark field q α (s) is described by a spinless Grassmann variable χ(s) [3], and the chiral condensate per flavor is evaluated as qq = χχ /4 in the continuum limit. The local chiral condensate can be calculated using the quark propagator of the KS fermion with a small quark mass m. The chiral-limit value is estimated by the chiral extrapolation of m → 0. As a technical caution, the chiral and continuum limits do not commute for the KS fermion at the quenched level, although this problem would be absent in full QCD [45].
For the gauge configuration U = {U µ (s)}, the Euclidean KS fermion propagator is given by with the color index i and j. This propagator is numerically obtained by solving the large-scale linear equation with a point source. The local chiral condensate for the gauge configuration {U µ (s)} is expressed with the propagator as Here, we consider the net chiral condensate by subtracting the contribution from the trivial vacuum U = 1 as where the subtraction term is exactly zero in the chiral limit m = 0. The global chiral condensate is obtained by taking its average over the space-time x and the gauge ensembles U 1 , U 2 , ..., U N , B. Chiral symmetry breaking in Abelian gauge system In this subsection, we analytically investigate relation between chiral symmetry breaking and the field strength in Euclidean Abelian gauge systems. For the simple argument, we consider Euclidean U(1) gauge systems with quasimassless Dirac fermions coupled to the U(1) gauge field, although it is straightforward to generalize this argument to Abelian projected QCD with U(1) 2 gauge symmetry.
In the U(1) gauge system, the chiral condensate is proportional to the functional trace of the fermion propagator, with the covariant derivative D µ ≡ ∂ µ + igA µ , the field strength F µν , σ · F ≡ σ µν F µν and σ µν ≡ i 2 [γ µ , γ ν ]. Note that D 2 − m 2 is a negative-definite operator, and all of its eigenvalues are negative. Since the trace of any odd-number product of γ-matrices is zero, we find with Because of the overall factor m in I, I goes to zero in the chiral limit m → 0, unless the denominator becomes zero in this limit.
Since the operator (D 2 − m 2 ) 2 in the denominator is positive definite, to realize the zero denominator in I in Eq. (20), we need a significant negative contribution from the other three terms including F, G, or [D 2 , σ · F ]. For instance, in the absence of the field strength, i.e., F µν ≡ 0, one finds near m 0 with the UV cutoff Λ and the degeneracy γ. According to the positive denominator in the integrand, I Fµν ≡0 has no IR singularity to cancel m of the numerator, and therefore I Fµν ≡0 goes to zero in the chiral limit of m → 0.
To cancel m in the numerator of I, we need a significantly large amount of the field strength so as to present zero mode in the denominator of I and to keep I non-zero in the chiral limit. Note here that F(≥ 0) always gives a negative (non-positive) contribution in the denominator of I, while the contribution from G or [D 2 , σ · F ] can be positive and negative. In fact, the magnetic quantity F can give the zero mode in the denominator of I, even without the contribution from G and [D 2 , σ · F ]. In contrast, in Euclidean Abelian gauge systems, F ≡ 1 4 F 2 µν = 0 means F µν = 0, and then G = [D 2 , σ · F ] = 0.
To conclude, the magnetic quantity F is expected to be significantly important to realize chiral symmetry breaking in Euclidean Abelian gauge theories, although, in some cases, the contribution from G and [D 2 , σ · F ] can assist the realization of chiral symmetry breaking.
In a special case of constant F µν , one finds [D 2 , σ · F ] = 0 for the Abelian system, and obtains because of trγ 5 = trσ µν = trγ 5 σ µν = 0. For more special case of a constant magnetic field, there occurs the Landaulevel quantization, and the spatial degrees of freedom perpendicular to the magnetic field is frozen in the lowest Landau level. This infrared effective low-dimensionalization of the charged spinor dynamics induces chiral symmetry breaking in the chiral limit [46][47][48], which is known as magnetic catalysis [49].

V. ABELIAN GAUGE SYSTEM WITH A STATIC MONOPOLE-ANTIMONOPOLE PAIR ON A LATTICE
In the QCD vacuum, complicated monopole world-lines generally emerge in the MA gauge [21,25], and therefore it is difficult to clarify the primary correlation with the chiral condensate among the magnetic objects such as monopoles, F and G.
In this section, to seek for the primary correlation with the chiral condensate, we create idealized Abelian gauge system with a monopole-antimonopole pair on a lattice, and investigate the relation among the local chiral condensate, monoples, and magnetic fields. Also, we consider a magnetic flux system without monopoles.
A. Static monopole-antimonopole pair systems To begin with, we deal with an idealized Abelian gauge system of a static monopole-antimonopole pair on a periodic lattice of the four-dimensional Euclidean space-time.
In the three-dimensional space R 3 , let us consider a static monopole-antimonopole pair with the distance of l in z-direction. To realize such a lattice gauge system, we set the Abelian link-variable u µ (s) to be otherwise u µ (s) = 1. Figure 1 shows the building-block plaquette to realize a static monopole-antimonopole pair on the lattice. Here, only the red link-variables take a nontrivial value of i.
As for the phase variable θ µ (s), which corresponds to the Abelian gluon, one finds otherwise θ µ (s) = 0. For the all-red plaqutte with s x = s y = 0 and 1 ≤ s z ≤ l, one gets which leads to the Dirac string of n xy (s) = −1 and zero field strength θ xy (s) = 0, because of the definition of the field strength θ µν (s) and the Dirac string n µν (s), Thus, the all-red plaquette induces the singular Dirac string at its center on the dual lattice. In fact, for the idealized system in Fig. 1(b), a Dirac string appears inside the all-red plaquette. At the terminal of the Dirac string, a monopole or an anti-monopole appears on the dual lattice, as shown in Fig. 2. Actually, the three-dimensional spatial cube including only one all-red plauette has a static (anti)monopole at its center (on the dual lattice), because only one n kl (s) has nonzero value of ±1 among the six independent plaquettes composing the cube, Thus, this idealized system includes a static monopole at ( 1 2 , 1 2 , 1 2 ) and a static anti-monopole at ( 1 2 , 1 2 , l + 1 2 ) in spatial R 3 .
This monopole and anti-monopole system has also physical magnetic flux around the line segment connecting the monopole pair. In fact, a physical magnetic field is created in the neighboring plaquette u xy (s) of the all-red plaquette in Fig. 1. In this idealized system, only the plaquette u xy (s) including one red link takes a nontrivial value as u xy (s) = −i = e −iπ/2 in case with one nontrivial link, (30) otherwise u µν (s) = 1. Note here that, by gauge transformation, the location of the Dirac string is generally changed, but the physical field strength is never changed. Figure 3 shows the local chiral condensate χχ(s) u and the magnetic quantity F ≡ 1 4 F 2 µν = 1 2 H 2 for l = 4 in the three dimensional space R 3 . In this demonstration, the quark mass is taken to be m = 0.01 in the lattice unit.
For this idealized static system, there actually appears a magnetic field H, i.e., nonzero flux of F = 1 2 H 2 > 0, in space between the monopole and the anti-monopole, and the local chiral condensate takes a significant value in the vicinity of the magnetic field. In contrast, one finds G = H · E = 0 everywhere, since only spatial plaquettes take a nontrivial value and E = 0. Thus, in this system, it is likely that the magnetic field stemming from monopoles has the primary correlation with the local chiral condensate.

B. Static magnetic flux system
Next, let us investigate a static magnetic flux system without monopoles. Owing to the spatial periodicity, the special case of l = L z in the static monopole-antimonopole system has no (anti)monopoles, because of the magneticcharge cancellation. In this special case of l = L z , there only exists a physical static magnetic flux along z-direction. Figure 4 shows the local chiral condensate χχ(s) u and the magnetic quantity F ≡ 1 4 F 2 µν = 1 2 H 2 for l = L z in spatial R 3 , taking the quark mass of m = 0.01 in the lattice unit.
Again, the local chiral condensate takes a significant value in the vicinity of the magnetic field H, i.e., nonzero flux of F = 1 2 H 2 > 0, even without (anti)monopoles. Note also that this system has G = H · E = 0 everywhere, because of E = 0. Therefore, in this idealized system, we conclude that the magnetic field or F has the primary correlation with the local chiral condensate.

VI. LATTICE QCD STUDY FOR LOCAL CHIRAL CONDENSATE, MONOPOLES, AND MAGNETIC FIELDS
In our previous study with lattice QCD, we observed a strong correlation between the local chiral condensate and monopoles in Abelian projected QCD [33]. As a possible reason of this correlation, we conjectured that the strong magnetic field around monopoles is responsible to chiral symmetry breaking in QCD, similarly to the magnetic catalysis [46][47][48][49].
In this section, using lattice QCD Monte Carlo simulations, we investigate the relation among the local chiral condensate, monopoles, and magnetic fields in Abelian projected QCD. In this paper, the SU(3) lattice QCD simulation is performed using the standard plaquette action at the quenched level. In each space-time direction, we impose the periodic boundary condition for link variables, and the anti-periodic boundary condition for quarks in order to describe also thermal QCD.
For the numerical Monte Carlo calculation, we basically adopt the lattice parameter of β ≡ 2N c /g 2 = 6.0 and the size V = 24 4 . The lattice spacing a 0.1 fm is obtained from the string tension σ = 0.89 GeV/fm [37]. Also, we adopt β = 6.0 and V = 24 3 × 6 for the high-temperature deconfined phase at T 330 MeV above the critical temperature. Using the pseudo-heat-bath algorithm, we generate 100 and 200 gauge configurations for V = 24 4 and 24 3 × 6, respectively. All the gauge configurations are taken every 500 sweeps after thermalization of 5,000 sweeps. MA gauge fixing is performed with the stopping criterion that the deviation ∆R MA / (4V ) becomes smaller than 10 −5 in 100 iterations. For the calculation of the local chiral condensate, we use the quark propagator of the KS fermion with the quark mass of m = 0.01, 0.015, 0.02 in the lattice unit, Here, the quark mass is taken to be finite, since the chiral and continuum limits do not commute for the KS fermion at the quenched level [45]. The jackknife method is used for statistical error estimates.
For each lattice gauge configuration of Abelian projected QCD in the MA gauge, we calculate the local monopole density ρ L (s), the local chiral condensate, and the Lorentz invariants F and G, defined in Section III

A. Distribution similarity between local chiral condensate and magnetic variables
To begin with, we pick up a gauge configuration generated in lattice QCD on V = 24 4 at β = 6.0, and investigate correlation between the local chiral condensate and magnetic variables. Figure 5 shows the local chiral condensate χχ(s) u with the quark mass of m = 0.02, the local monopole density ρ L (s), and the Lorentz invariants F(s) and |G(s)|, respectively, as well as the monopole location in the space R 3 at a time slice, for a typical gauge configuration of Abelian projected QCD.
From Fig. 5(a), one finds that the local chiral condensate tends to take a large value near the monopole location [33]. Since monopoles appear on the dual lattice, we show the local monopole density ρ L (s), as the average on closest dual sites. Of course, ρ L (s) takes a large value near the monopole. The distribution of the the local monopole density resembles that of the local chiral condensate, as was pointed out in Ref. [33]. Figures 5(c) and (d) show the Lorentz invariants F and G, respectively. As a new result in this paper, we find that the distributions of F and G also resemble that of the local chiral condensate.
The close relation of monopoles with F and G might be understood, since the field strength tensor relates to monopoles as ∂ µF i µν = k i ν . Roughly speaking, the monopole can be a kind of source of F and G. In contrast, their similarity with the local chiral condensate is fairly nontrivial.
In any case, we find clear correlation of distribution similarity among the local chiral condensate, the local monopole density, and the Lorentz invariants F and G in Abelian projected QCD in the MA gauge.

B. Correlation coefficients between local chiral condensate and magnetic variables
In this subsection, we quantify the similarity between the local chiral condensate χχ(s) u and magnetic variables, i.e., ρ L (s), F(s) and G(s), defined in Section III. To this end, we use all the generated 100 gauge configurations in lattice QCD on V = 24 4 at β = 6.0, and calculate the local chiral condensate at 2 4 distant space-time points for each gauge configuration, resulting 1600 data points at each quark mass. Figure 6 shows the scatter plot between the local chiral condensate χχ(s) u and magnetic variables, i.e., the local monopole density ρ L (s), Lorentz invariants F(s) and |G(s)|, respectively, using 100 gauge configurations of Abelian projected QCD in the MA gauge, with the quark mass of m = 0.01, 0.015, 0.02 in the lattice unit. In Fig. 6, positive correlation is qualitatively found between the local chiral condensate and the magnetic variables, ρ L , F and |G|, respectively.
Next, we consider a quantitative analysis using correlation coefficients between the local chiral condensate and the magnetic variables, as a statistical indicator of correlation. In general, for arbitrary two statistical ensembles {A i } and {B i }, their correlation coefficient r is defined as We measure correlation coefficients between the local chiral condensate | χχ(s) u | and three magnetic variables, ρ L (s) α , F(s) α and |G(s)| α , at various exponent α, using 100 gauge configurations of Abelian projected QCD in SU(3) lattice QCD at β = 6.0 on V = 24 4 , for the quark mass m = 0.01 in the lattice unit. Table I shows the result for the correlation coefficients. Quantitatively, the magnetic quantity F has the strongest correlation with the chiral condensate rather than ρ L and G. As a conclusion of this paper, we find a strong positive correlation of r 0.8 between the local chiral condensate | χχ(s) u | and the magnetic quantity F(s) in the confined vacuum of Abelian projected QCD.
C. High-temperature deconfined phase Finally, we investigate also a high-temperature deconfined phase in lattice QCD on V = 24 3 × 6 at β = 6.0, where the temperature is T 330 MeV above the critical temperature. We generate 200 gauge configurations, and calculate the local chiral condensate at 2 3 distant space points at a time slice for each gauge configuration, resulting 1600 data points at each quark mass. Figure 7 shows the scatter plot between the local chiral condensate | χχ(s) u | and magnetic variables, i.e., the local monopole density ρ L (s), and Lorentz invariants F(s) and |G(s)|, respectively, using 200 gauge configurations of Abelian projected QCD in the MA gauge, with the quark mass of m = 0.01, 0.015, 0.02 in the lattice unit.
We show in Table II correlation  From Fig. 7 and Table II, all the correlations between the local chiral condensate and the three magnetic variables, ρ L , F and G, become weaker in the deconfined phase, where the chiral condensate itself goes to zero in the chiral limit.

VII. SUMMARY AND CONCLUSION
We have studied the relation among the local chiral condensate, monopoles, and magnetic fields, using the lattice gauge theory, as a continuation of Ref. [33].
First, we have created idealized Abelian gauge systems of 1) a static monopole-antimonopole pair and 2) a magnetic flux without monopoles, on a four-dimensional Euclidean lattice. In these systems, we have calculated the local chiral condensate on quasi-massless fermions coupled to the Abelian gauge field, and have found that the chiral condensate is localized in the vicinity of the magnetic field.
Second, performing SU(3) lattice QCD Monte Carlo simulations, we have investigated Abelian projected QCD in the maximally Abelian gauge, and have found clear correlation of distribution similarity among the local chiral condensate, color monopoles, and color magnetic fields in the Abelianized gauge configuration.
As a statistical indicator, we have measured the correlation coefficient r, and have found a strong positive correlation of r 0.8 between the local chiral condensate and the Euclidean color-magnetic quantity F.
We have also examined the local correlation in the deconfined phase of thermal QCD, and have found that the correlation between the local chiral condensate and magnetic variables becomes weaker.
Thus, in this paper, we have observed a strong correlation between the local chiral condensate and magnetic fields in both idealized Abelian gauge systems and Abelian projected QCD. From these results, we conjecture that the chiral condensate is locally enhanced by the strong color-magnetic field around the monopoles in Abelian projected QCD, like magnetic catalysis.
Note however that this correlation does not necessarily mean that chiral symmetry breaking is caused by the nonuniform magnetic field. To realize spontaneous chiral-symmetry breaking, as was discussed in subsection IV B, we need some zero mode in the denominator of I in the chiral limit. In the context of the dual superconductor picture, this might be realized by condensation of monopoles, as was suggested in the dual Ginzburg-Landau theory [28].
To conclude, once chiral symmetry is spontaneously broken, the local chiral condensate is expected to have a strong correlation with the color magnetic field.