PRDPRVDAQPhysical Review DPhys. Rev. D2470-00102470-0029American Physical Society10.1103/PhysRevD.97.034505ARTICLESLattice field theories, lattice QCDMeson masses in electromagnetic fields with Wilson fermionsMESON MASSES IN ELECTROMAGNETIC FIELDS WITH …BALI et al.BaliG. S.1,2,*BrandtB. B.3,†EndrődiG.3,‡GläßleB.1,§Institute for Theoretical Physics, Universität Regensburg, D-93040 Regensburg, GermanyTata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005, IndiaInstitute for Theoretical Physics, Goethe University, Max-von-Laue-Strasse 1, 60438 Frankfurt am Main, Germany
gunnar.bali@ur.de
brandt@th.physik.uni-frankfurt.de
endrodi@th.physik.uni-frankfurt.de
benjamin.glaessle@ur.de
12February20181February201897303450525October2017Published by the American Physical Society2018authorsPublished by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Funded by SCOAP3.
We determine the light meson spectrum in QCD in the presence of background magnetic fields using quenched Wilson fermions. Our continuum extrapolated results indicate a monotonous reduction of the connected neutral pion mass as the magnetic field grows. The vector meson mass is found to remain nonzero, a finding relevant for the conjectured ρ-meson condensation at strong magnetic fields. The continuum extrapolation was facilitated by adding a novel magnetic field–dependent improvement term to the additive quark mass renormalization. Without this term, sizable lattice artifacts that would deceptively indicate an unphysical rise of the connected neutral pion mass for strong magnetic fields are present. We also investigate the impact of these lattice artifacts on further observables like magnetic polarizabilities and discuss the magnetic field–induced mixing between ρ-mesons and pions. We also derive Ward-Takashi identities for QCD+QED both in the continuum formulation and for (order a–improved) Wilson fermions.
Deutsche Forschungsgemeinschaft10.13039/501100001659EN 1064/2-1SFB/TRR 55Universität Regensburg10.13039/501100005626INTRODUCTION
Background magnetic fields have a decisive impact on the physics of quarks and gluons and offer a wide range of applications. Strong magnetic fields appears in noncentral heavy-ion collisions [1,2], inside magnetars [3], and might have been generated during the evolution of the early Universe [4]. An important characteristic of the magnetized QCD medium—relevant for each of the above settings—is the effect of the magnetic field on the spectrum of the theory. It has been speculated [5,6] that long-lived magnetic fields might affect hadronization in heavy-ion collisions. In strongly magnetized neutron stars, the magnetic properties of QCD matter, predominantly those of the light hadrons, may have a significant impact on the mass-radius relation [7].
The direct influence of the external field on the phase structure of QCD has received considerable attention recently; see for instance the review [8]. In particular, lattice simulations have revealed that the transition temperature of the chiral restoration/deconfinement crossover monotonously decreases as the magnetic field grows [9,10]; see also Refs. [11,12]. Thinking in terms of a hadronic model, in which the temperature needed to resolve the quark structure of a hadron is of the order of the hadron mass, it is plausible that the magnetic field dependence of the masses of the lightest hadrons do play a relevant role for the structure of the phase diagram. Specifically, pions are expected to be dominant in this respect. It is well known from chiral perturbation theory [13] that charged and neutral pions respond oppositely to the magnetic field B=|B|. While lattice results employing staggered [9] and quenched unimproved Wilson [14] quarks agree that the charged pion mass mπ± is increased by the presence of a magnetic field, there is a discrepancy in the literature about mπ0(B) for intermediate magnetic fields. It monotonously decreases for quenched overlap quarks [15], while for quenched unimproved Wilson fermions, a turning point and subsequent increase was reported in Ref. [14].
In fact, the above mentioned lattice results are based on the connected contribution to the pion correlator, which corresponds to a hypothetic meson with exclusively u¯u or d¯d flavor content. We get back to this point in Sec. II B below.
In light of the above discussion, the two results hint at opposite implications for the QCD phase diagram for strong magnetic fields.
Clarifying this question may also be relevant for yet another aspect of the phase structure of QCD with B>0: It has been conjectured that strong magnetic fields could reduce the vector meson mass to zero and lead to the condensation of ρ-mesons [16]. They may also induce a transition to a superconducting phase [17].
Note that this is the simplest scenario for the superconductivity of the QCD vacuum when exposed to external magnetic fields. In fact, when coupled to the photon field that exhibits gauge fluctuations, the would-be massless ρ mode is absorbed by the photon field in accordance with the Higgs mechanism [18] so that no massless mode remains. For constant background magnetic fields—like the one we are working with here—this mechanism is absent. Also note that the proposed superconducting phase was suggested to exhibit several unusual features like anisotropy and inhomogeneity [17,19], and we do not attempt to capture these subtle characteristics here, which would require dynamical QCD+QED simulations and is much more involved.
On the one hand, QCD inequalities can be used to show that massless ρ-mesons are only allowed if the (connected) neutral pion mass vanishes as well. The existence of a turning point in mπ0(B), above which it grows with increasing B, would thus speak against this scenario. On the other hand, a study using quenched overlap quarks in two-color QCD, based on spatial vector correlators, has reported a critical magnetic field of around 0.9GeV2, beyond which ρ-condensation was found to occur [20].
Resolving the above discrepancy regarding mπ0(B) is thus clearly of high importance. In Ref. [21], we published preliminary results that explain the origin of the difference between the quenched overlap and Wilson results. We have found that Wilson quarks are subject to a B-dependent additive mass renormalization at finite lattice spacings. Although this effect disappears in the continuum limit, it has a drastic impact on the meson masses evaluated at nonzero lattice spacings. We also sketched a method to subtract these artifacts by determining a magnetic field–dependent line of constant physics for the Wilson bare mass parameter.
In this paper, we expand on the ideas introduced in Ref. [21] and discuss the improvement of scaling toward the continuum limit via the B-dependent additive mass renormalization in more detail. Performing this improvement reveals that the neutral pion mass monotonously decreases as B grows and, thus, resolves the existing disagreement between the quenched Wilson and overlap formulations. We also determine the ρ-meson mass and find that it remains nonzero in the whole range of magnetic fields studied, 0<eB<4GeV2. For comparison, we also measure the connected neutral pion mass using dynamical staggered quarks on existing configurations from the study of Ref. [9].
Note that the use of staggered fermions for the spectroscopy of vector mesons at finite B is complicated, since the magnetic field affects the group theoretical construction of vector meson interpolators. Furthermore, the use of the staggered formulation complicates the interpretation of higher-lying states due to taste splitting.
There are further interesting quantities that are related to the spectrum, including magnetic polarizabilities and magnetic moments. These describe the leading-order response of hadrons to the external field and, thus, involve derivatives at zero field; see, e.g., Refs. [15,22,23]. In contrast to observables at finite values of the field, these can also be defined for background electric fields, giving rise to electric polarizabilities. We will show that for Wilson quarks both types of polarizabilities suffer from enhanced lattice artifacts that can (and should) be eliminated using our improvement program.
We emphasize that our results are obtained in the quenched approximation, just as all studies in the literature so far concerning the light meson spectrum at finite B (except for the staggered results of Ref. [9], which, however, focused on pions only). The quenching of virtual quark loops induces a systematic effect that is hard to estimate. At B=0, it is known that the results for the quenched and the dynamical spectrum agree within about 10% (see Refs. [24,25], for instance), indicating that the meson masses are dominantly determined by the properties of valence quarks. This is expected to remain true in the presence of external fields. For example, the mass of the neutral pion, being a pseudo-Goldstone boson even for B>0, depends primarily on the quark masses and not on gluonic properties. For further mesons, one may expect a moderate quenching effect as long as the vacuum exhibits a structure similar to that at B=0. Observables related to the bulk properties of the gluonic field (such as the static quark-antiquark potential considered in Ref. [26]), however, are expected to show stronger quenching effects, since the primary effect of the external field enters only indirectly via the sea quark loops which are not present in the quenched setup.
This paper is organised as follows. We start with a few remarks about the meson spectrum in continuum QCD in the presence of background magnetic fields in Sec. II. This is followed by Sec. III, where our lattice setup is described and our improvement scheme relevant for Wilson fermions is detailed. Section IV contains our results about the spectrum and the analysis of the dependence on quark masses and on the lattice spacing. In Sec. V, we consider the impact of our improvement scheme for magnetic polarizabilities before we conclude in Sec. VI. Three Appendixes contain our results in the free case (Appendix A), the derivation of Ward-Takahashi identities for nonzero electromagnetic fields (Appendix B), and a perturbative check of the B-independence of the multiplicative mass renormalization constant in QCD (Appendix C).
MESON SPECTRUM IN THE CONTINUUMMeson and quark masses for B>0
We begin with a few general comments about the magnetic field dependence of the pion masses in the continuum. For small magnetic fields (and at zero temperature), we may approximate the charged pion as a pointlike free scalar particle with electric charge ±e. In that case, the quantum mechanically allowed energies are the Landau levels, En2=mπ2+(2n+1)|eB|+pz2,n∈Z0+,where mπ is the B=0 pion mass and we have assumed B=Bez. Thus, the magnetic field–dependent mass (the lowest energy with n=pz=0) reads mπ±(B)=mπ2+|eB|.The same argument for the neutral pion gives a magnetic field–independent mass, mπ0(B)=mπ. Thus, the neutral pion remains massless for nonzero magnetic fields if mπ=0, while the charged pions become massive for B≠0.
It turns out that the conclusion of increasing charged pion masses and an almost constant neutral pion mass remains valid if the interaction between pions is taken into account. This can be done consistently within chiral perturbation theory for nonzero magnetic fields [13,27–29]. Indeed, the neutral pion remains a Goldstone boson even for B≠0. In addition, away from the chiral limit, the Gell-Mann-Oakes-Renner relation mf=mπ02Fπ022ψ¯ψ+O(mπ04)still applies and is independent of B[27,28]. Here, Fπ0 is the decay constant of the neutral pion, ψ¯ψ is the (in our convention, positive) chiral condensate, and mf denotes the (degenerate) up and down quark masses. Fπ0 and ψ¯ψ above are defined in the chiral limit, mπ0→0. Equation (3) is consistent with the notion that the quark mass does not depend on the magnetic field.
The independence of the quark masses on the magnetic field also becomes obvious in a completely different limit: at very high temperatures. Here, quarks become quasifree due to asymptotic freedom. Then, the quantum mechanical energy levels are the Landau levels for fermions with charge q, Ef;n2=mf2+2n|qfB|+pz2,n∈Z0+,so that the lowest energy (n=pz=0) is indeed magnetic field independent.
Clearly, if QCD interactions are present, the levels will mix. Nevertheless, it turns out that the lowest Landau level can still be defined unambiguously and remains approximately B independent [30,31].
The free case is of interest for yet another reason: In the strong magnetic field limit B→∞, quarks and gluons decouple, and pure fermionic observables approach their free-case values due to asymptotic freedom (see, for example, Refs. [10,32]). It is therefore instructive to compare our QCD results at large magnetic fields to the free case; see below.
In QCD, the quark mass is also subject to a multiplicative renormalization of the form mfr=Zmmf. Since Zm is related to the ultraviolet behavior of the theory, while the magnetic field is a physical infrared parameter, one expects that Zm is also independent of B. We demonstrate this by an explicit calculation using perturbation theory in Appendix C.
In the following, we will also investigate the effect of the magnetic field on ρ-mesons. Unlike for pions, for vector mesons, the magnetic field interacts directly with the spin of the particle. In the free case, the energies of a pointlike vector meson (with a gyromagnetic ratio g=2 and charge q=±e) are given by [19]En2=mρ2+(2n+1)|eB|-gszqB+pz2,where sz=0, ±1 is the projection of the spin on the magnetic field axis. For sz=1, one may expect the mass of the associated ρ-meson to vanish when eB=mρ2. This has lead to the proposal that the ρ-mesons could condense [16] and initiated the investigation of the possibility of a superconducting vacuum for strong magnetic fields [17], as mentioned in the Introduction. However, QCD inequalities imply [14] that the correlation function in the vector channel is bounded from above by the (connected) neutral correlation function of the pion, so that mρ≥(mπu+mπd)/2.Thus, a vanishing ρ-meson mass is only possible when mπu/πd=0.
Note that the ρ-meson is a resonance in nature, rather than a stable particle, so the extraction of its resonant mass from lattice simulations demands a careful finite size analysis. In our study, however, the ρ-meson can be considered to be a stable particle, since we are working in the quenched approximation.
Mixing of meson states for B>0
When we consider the neutral pion at B=0, where the SUV(2) symmetry is intact, the associated state is a mixture of neutral pions with u¯u and d¯d flavor contents, |π0⟩=12(|πu⟩-|πd⟩).At B>0, however, SUV(2) is broken explicitly by the different quark electric charges, qd=-qu/2=-e/3, where e>0 is the elementary charge. In this case, the neutral pion will be given by the more general relation |π0⟩=α(B)|πu⟩-β(B)|πd⟩withα2(B)+β2(B)=1,where α(B)→1/2 and β(B)→1/2 for B→0. Thus, determining quantities related to the neutral pion in principle involves the computation of the coefficients α and β for each value of the magnetic field. Note that for B>0, the associated correlation matrix contains disconnected diagrams (these terms cancel at B=0 due to isospin symmetry). In the present study, we neglect disconnected diagrams and the associated mixing and work with the individual πu and πd-states instead. This approximation is expected to become valid for strong magnetic fields, where the free case is approached (see the discussion above) and, thus, disconnected diagrams become negligible. Note that the so-defined connected pion states πu and πd are still useful as, for example, they enter in QCD inequalities [14]. We remark furthermore that neutral ρ-mesons are affected by the above issue in the same way.
The presence of a finite external magnetic field enables the mixing between the pion and the ρ-meson with sz=0, B≠0:ρsz=0±,0≡ρ0±,0↔π±,0.Throughout the paper, superscripts of meson states denote electric charge, and subscripts denote the spin. For sz=±1, the above mixing is forbidden due to the conservation of angular momentum. Consequently, pion states contribute to the spectral representation of correlation functions associated with sz=0ρ-mesons, and, in turn, these ρ-meson states also contribute to the spectral representation of pion correlation functions, so both correlation functions show the same leading-order exponential falloff. In fact, once this mixing is enabled, the physical mass eigenstates are mixtures of would-be π and ρ-states which can be written as |(π′)±,0⟩=cos(θ)|π±,0⟩+sin(θ)|ρ0±,0⟩and|(ρ′)0±,0⟩=-sin(θ)|π±,0⟩+cos(θ)|ρ0±,0⟩,where θ is the mixing angle and we denote the state with smaller energy by π′. In principle, the corresponding energies can be extracted from a correlation matrix using the operators for π±,0 and ρ0±,0, while θ will depend on the interpolators in use. On top of these two (and higher excited single-particle) states, there are also states containing more pions (and radial excitations) that may be involved in the mixing. However, their contribution is negligible if mρ′<E1 is satisfied, where E1 is the energy of the lowest-lying multiparticle excitation.
MAGNETIC FIELDS AND QUARK MASSES WITH WILSON FERMIONS ON THE LATTICELattice setup
We consider two flavors of (quenched) unimproved Wilson fermions on the lattice with the fermion action SF=a4∑x[ψ¯(x)(M0+4ra)ψ(x)-12a∑μ=04(ψ¯(x)(r-γμ)Uμ(x)ψ(x+aμ^)+ψ¯(x+aμ^)(r+γμ)Uμ†(x)ψ(x))].Here, a denotes the lattice spacing, ψ=(u,d)⊤ is the Dirac field for flavor u and d, and M0 is the matrix of bare quark masses, M0=12(mu0+md0)1+12(mu0-md0)τ3,where 1 is the unit matrix in flavor space and τ3 is the third Pauli matrix. γμ denote the Euclidean Dirac matrices. r is the coefficient of the Wilson term which we will set to r=1 in the simulations. The quarks couple to gluonic and electromagnetic fields via the link variables Uμ(x), containing the gluonic links UμG(x)∈SU(N) and the electromagnetic link variables uμ(x), Uμ(x)=UμG(x)uμ(x)∈U(N)×SUf(2)withuμ(x)=exp[iea(16+τ32)Aμ(x)].The latter are matrices in flavor space since the electromagnetic charges are different. We are mostly interested in QCD in the presence of an external magnetic field B=Bez. Such a magnetic field can be generated by a vector potential of the form Ax(x)=-ηBy-(1-η)ByLxδ(x-Lx),Ay(x)=(1-η)Bx+ηBxLyδ(y-Ly)andAt(x)=Az(x)=0,where the parameter η cancels in physical observables and thus can be chosen at will. In the following, results are obtained using the “symmetric gauge,” η=1/2, but we have explicitly checked that the results agree when we set η=1. Lμ=Nμa denotes the spatial extent in the μ-direction, and Nμ is the number of lattice points. The boundary δ-terms ensure periodic boundary conditions for the gauge potential and the constancy of the magnetic field in the x, y-plane [33]. Note, that in a finite volume, the magnetic field is quantized according to qdBNxNy=2πNb,withNb∈Zand0≤Nb<NxNy,where the smallest quark charge, qd=-e/3, appears.
Quark mass renormalization
The term proportional to r in Eq. (11) is the well-known Wilson term, which is introduced to remove the doublers from the theory. This term breaks chiral symmetry explicitly so that the quark mass is subject to additive renormalization, besides the multiplicative renormalization that is already necessary in the continuum. From the arguments of Sec. II, we expect the multiplicative renormalization factor Zm to be independent of the external field—or, in other words, a B-independent renormalization scheme can be chosen (see Appendix C). Thus, at finite lattice spacing, the renormalized quark mass of flavor f is given by
Note that a similar equation holds for the case of the coupling to a dynamical electromagnetic field (e.g., Ref. [34]), mfr=Zm(a,e)[mf0-mc;f(a,e)], where both renormalization factors depend on the bare electromagnetic coupling e.
mfr=Zm(a)mf=Zm(a)[mf0-mc;f(a,B)],where amc;f(a,B)=amc(a,0)+dma2|qfB|+O([a2|qfB|]2)and dm is the coefficient of the leading-order B-dependent term. Note that, although dm corrects for a lattice artifact, it is not a Symanzik improvement coefficient. However, it can have a nontrivial value since our theory is not O(a2) [and not even O(a)] improved. In Appendix A, we give an intuitive explanation of why magnetic fields shift the quark mass, i.e., why mc;f becomes B dependent. The argument boils down to the fact that the Wilson term coincides with the action of a scalar particle—multiplied by the lattice spacing so that it vanishes for a→0. In the presence of the magnetic field, the squared mass of this scalar particle is shifted by B according to Eq. (2). Eventually, this induces the B-dependence in the additive mass renormalization. This can be checked explicitly in the free case (where dm=-1/2), for which we show results in Appendix A.
It is worth stressing that the presence of the additive quark mass renormalization is a lattice artifact and mc;f(a,B)→0 for a→0 for any magnetic field. For Wilson fermions, it is customary to introduce the hopping parameter κf and its critical value κc;f, κf=(2amf0+8)-1andκc;f(a,B)=(2amc;f(a,B)+8)-1,so that the B-dependent quark mass renormalization can be translated to the B-dependence of κc,f, 1κc;f(a,B)=1κc(a,0)+2dma2|qfB|+O([a2|qfB|]2).
At finite lattice spacing, the B-dependence of mc;f contaminates physical B-effects. Due to this additional lattice artifact, very fine lattices are necessary for a reliable continuum extrapolation. A better way is to get rid of the contamination by keeping the simulation on the line of constant physics (LCP), characterized by a constant renormalized quark mass mfr. This LCP depends on the magnetic field according to Eq. (16) and will be referred to as LCP(B) from now on.
We remark that, for unimproved Wilson fermions, lines of constant physics are, in general, only defined up to magnetic field–dependent or –independent O(a) lattice artifacts. Here, we choose to correct for the B-dependent artifacts originating from the additive mass renormalization, since these turn out to be numerically large.
Since Zm is, in our definition, B independent, it suffices to tune the hopping parameter κf=κf(a,B) such that the bare (multiplicatively unrenormalized) quark mass mf≡[mf0(a,B)-mc;f(a,B)]=12a[1κ(a,B)-1κc;f(a,B)]remains constant when varying B. We will discuss the impact of this tuning on lattice artifacts in detail in Sec. IV E, where we compare the results obtained when keeping κf constant while varying B to those obtained along the LCP(B). The B-dependence of the quark mass is also inherited by further quantities that are related to derivatives with respect to the magnetic field. Particular examples, namely magnetic moments and polarizabilities, will be discussed in Sec. V.
Before we move on, let us collect some theoretical expectations concerning the behavior of mc;f(a,B). Parity symmetry ensures that mc;f(a,B)=mc;f(a,-B). Below, we show numerical results in the quenched setup, where charged quarks are only included in the valence sector. In this case, |qu/qd|=2 leads to mc;u(a,B)=mc;d(a,2B)or∂mc;u(a,B)∂B|B=0=2∂mc;d(a,B)∂B|B=0.
Determination of the critical mass
The key ingredient to fixing the LCP(B) is the B-dependence of κc;f. To determine it, we can make use of the fact that the neutral pion mass vanishes in the chiral limit, even at finite B (cf. Sec. II). In our quenched setup, this chiral extrapolation is complicated by the appearance of chiral logarithms in quenched chiral perturbation theory [35,36] (see also the right panel of Fig. 1 below). It is more convenient to use the current quark mass m˜f obtained from the axial Ward-Takahashi identity (WI). In Appendix B, we derive the axial and vector WIs in the continuum and for Wilson fermions for B>0. While the charged WIs acquire additional B-dependent terms,
See also Refs. [37–39].
the neutral WIs are found to keep their B=0 form, am˜f(B)=∂0⟨(JA)0f(x0)Pf(0)⟩2⟨Pf(x0)Pf(0)⟩withx0≠0.Here, (JA)f is the (point-split) axial current operator with fermionic fields of flavors f [cf. Eq. (B24)], and Pf is the associated pseudoscalar density. The isosinglet WI suffers from the well-known axial anomaly, see Eq. (B17), which involves only gluonic contributions for background magnetic fields, and obtains contributions from disconnected diagrams, which we neglect. This corresponds to pretending that we have two different flavors carrying the same electric charge, in analogy to replacing π0 by πu, thereby avoiding disconnected diagrams. In the following, we will thus determine κc;f by requiring that the “neutral” current quark mass m˜f, defined via Eq. (22), vanishes at κc;f(B) for each value of B.
110.1103/PhysRevD.97.034505.f1
Results for am˜u(a,B;κ) (left) and amπu(a,B;κ) (right) versus 1/κ-1/κc;u(a,0) for several values of Nb on the β=6 ensemble. The colored bands in the left figure are the results from linear fits, and the diamonds on the am˜u(a,B;κ)=0 line in both panels are the results for 1/κc;u(a,B)-1/κc;u(a,0). The error bars of the extrapolated results include the estimate for the systematic uncertainty explained in the text and are smaller than the symbols.
The correlation functions have been projected to zero momentum in the z-direction. In addition, we have summed the correlation functions over the x, y-plane at the sink. For B=0, this corresponds to projecting to zero momentum in the x- and y-directions. On the level of the electromagnetic vector potential, translational invariance is broken, and Fourier transformation is not well defined in the x, y-plane (see also Ref. [14]). We found that the summation enhances statistics but does not affect our results, as we have checked explicitly by comparing the results from summed correlation functions to those obtained without summation at a fixed position in the x, y-plane. For each value of a and B, we measured m˜f for several values of κ and determined κc;f(a,B) via a linear extrapolation of the form [cf. Eq. (20)] m˜f(a,B;κ)=Cf(a,B)·mf=Cf(a,B)2a[1κ-1κc;f(a,B)].Here, Cf is a proportionality constant, which relates the two estimates for the bare quark mass to each other and, thus, may contain (B-dependent) lattice artifacts. It can be written as Cf(a,B)=Z(a)(1+dZa2|qfB|+O([a2|qfB|]2)),where we have introduced the B=0 proportionality factor (see, e.g., Ref. [40]) Z=ZmZPZAand dZ is the coefficient of the leading-order B-dependent lattice artifact in Cf.
We note in passing that ideally one would like to use the charged axial and vector WIs for the tuning of sums and differences of quark masses. The corresponding identities are worked out including electromagnetic interactions [and also O(a)-improvement] in Appendix B. In this case, no disconnected diagrams contribute. We have also looked at the charged WIs but could not detect plateaus before the signal got lost in the noise. A potential reason for this is that the correlation functions used to construct the charged WIs decay with a mass larger than the smallest mass in the spectrum (since mπ+>mπ0) and thus suffer from the standard signal-to-noise problem. We plan to address this and investigate the charged WIs in more detail in the future.
Lattice setup and lines of constant physics
The simulations are performed in the quenched setup using the Wilson plaquette action [41] and heat bath and over-relaxation updates with a ratio of 1∶4. The generation of configurations and the measurements of correlation functions have been carried out employing a modified version of CHROMA[42]. The run parameters are collected in Table I. We also list the value of the Sommer parameter r0/a obtained from the interpolation provided in Ref. [43]. To convert to physical units, we use r0=0.5fm. The lattice sizes have been tuned to keep the physical volume fixed. This also allows us to achieve the same physical magnetic field—complying with the quantization condition (15)—for the different lattice spacings. The correlation functions have been computed using a Wuppertal smeared [44] source including (three dimensionally) APE smeared [45] link variables and a point sink. For the inversions, we have used the DFL-SAP-GCR solver introduced in Refs. [46,47].
I10.1103/PhysRevD.97.034505.t1
Details of the quenched ensembles. The results for κc(B=0) have been obtained from the linear extrapolation of the data for Nb=0; the first uncertainty is statistical, while the second results from the systematic uncertainty associated with the extrapolation, estimated by using a polynomial of second degree for the extrapolation.
To determine κc,f, we have measured m˜u and m˜d for several values of κ with pion masses between 1 GeV and about 400 MeV on each of these ensembles and for magnetic fields 0<eB<4GeV2. For B=0, the above setup satisfies mπfL≳3 for all measurement points. However, as we will see below, the neutral pion mass decreases significantly with increasing magnetic field, so that finite size effects can potentially increase as B grows. To show that this is not the case, we study finite size effects in Sec. IV F.
We have extracted κc;f(a,B) using a linear chiral extrapolation of the form given in Eq. (23). We show the results for am˜u(a,B;κ) together with the chiral extrapolations for several values of B obtained on the β=6.0 ensemble in Fig. 1 (left). As can be seen from the plot, the value of κc;u(a,B) has a strong dependence on B, while the slope Cu is only mildly affected. In fact, the strong change in κc forced us to measure m˜u for values of κ which would have lead to a negative value for mu at B=0. In Fig. 1 (right), we show the results for the mass of the neutral pion with u¯u flavor content for the smallest values of Nb. The plot shows that a chiral extrapolation using the larger pion masses would lead to similar results. However, in the region of small pion masses, we observe deviations from the straight line, possibly due to the presence of chiral logarithms. We show the behavior of Cu(a,B) in Fig. 2. Toward the continuum limit, Cu approaches unity as is clearly visible in the plot. In addition, for large values of B, we expect Cu to approach the free fermion case for which again Cu→1. This is indicated by the increase in Cu with increasing B. The increase seems to level off at around eB≈2GeV2, which could be a sign for growing lattice artifacts.
210.1103/PhysRevD.97.034505.f2
Results for Cu versus |eB| for different values of the lattice spacing. The dashed line indicates the continuum limit value Cu=1, and the grey solid line is the B=0 line.
The results for κc;u and κc;d are tabulated in Table II. To estimate the systematic uncertainty of the linear ansatz within the present range of quark masses, we have repeated the procedure by including a term quadratic in 1/κ-1/κc;f in Eq. (23). The spread of the results from these two fits is given as the systematic uncertainty in Table II. To visualize the effect of B on mc;f, we introduce the parameter Δmc;f(a,B)=mc;f(a,B)-mc(a,0)mf(a,0;mπ=415MeV).In this ratio and in our quenched setup, the leading B=0 lattice artifacts and renormalization factors cancel, so this is the ideal tool to see to what extent the renormalized quark mass is affected by the shift in mc;f(a,B). Note, that if Δmc;f(a,B)=-1, this means that the B-dependent part of the additive quark mass renormalization is negative and as large in magnitude as the quark mass associated with a pion with a mass of 415 MeV at B=0.
II10.1103/PhysRevD.97.034505.t2
Results for κc;u and κc;d obtained from the linear chiral extrapolation of the current quark masses. The first uncertainty is statistical, and the second is the estimate for the systematic uncertainty obtained from the spread to the result from an extrapolation with a polynomial of second order.
The results for Δmc;u and Δmc;d are displayed for different values of B and the different ensembles in Fig. 3. In the left panel of the figure, we show Δmc;u versus eB for the different lattice spacings [Δmc;d follows from these results and Eq. (21)]. The plot indicates that the mass shift is positive for low magnetic fields, in contrast to the negative shift in the free case [see Eq. (A1)]. However, the effects of higher-order terms in B soon become important, and Δmc;u(a,B) becomes negative and linear in B for strong magnetic fields—again in accordance with the expectation that we approach the free case in this limit. The figure also shows that the magnitude of Δmc;u(a,B) decreases with decreasing lattice spacing—reflecting the fact that the magnetic field dependence of mc;f is a lattice artifact and thus vanishes in the continuum limit. The right panel displays the results in the (Δmu,Δmd)-plane, and the individual data points indicate different values of B, starting from B=0 (solid black square) and proceeding along the colored lines which are included to guide the eye.
310.1103/PhysRevD.97.034505.f3
Left: Results for Δmc;u(a,B) as defined in Eq. (26) versus B for different values of the lattice spacing. Right: Results for Δmc;f(a,B) for different values of B (starting from the black point—B=0—outward) and different lattice spacings. The black dashed line corresponds to the slope of 2, which we expect for B≈0 due to Eq. (21), and the colored lines simply connect the data points and are included to guide the eye. The points marked with black circles correspond to a fixed physical magnetic field of eB=3.3GeV2.
The snail shell–like shape of the curves mapped out in this parameter space demands some discussion. First, all curves start with a slope which is close to 2 (indicated by the dashed black line). This is expected due to the relation given in Eq. (21), which fixes the factor between the derivatives of mc;u and mc;d at B=0. The direction of movement in the plane, however, is determined by the prefactor of the first term appearing in the expansion of mc;f(a,B), Eq. (17), with respect to B, and we have seen above that it is positive. Deviations from the black line indicate the onset of higher-order terms in the expansion of mc;f(a,B). For large magnetic fields, we expect to approach the free theory (see Sec. II), where the slope is again 2, cf. Appendix A. This behavior is also visible in Fig. 3. Next, let us discuss the change in shape for a→0. While the slopes for B→0 and B→∞ are dictated by Eq. (21) and the free theory case, respectively, the extent of the curve depends on the lattice spacing. Remember that the B-dependence of mc;f(a,B) is a lattice artifact, so the coefficients of the expansion of mc;f(a,B) become smaller for a→0, as can be seen in the left panel of Fig. 3. Consequently, the movement along the curve as a function of B becomes slower, and in the continuum limit, the curve will collapse to the starting point. To see this, compare the points of the curves marked with the black circles—these all correspond to the same physical magnetic field eB=3.302GeV2 (Nb=10).
Using the results for κc;u and κc;d tabulated in Table II, we can now tune κu and κd along the LCP(B), so that the additively renormalized quark mass from Eq. (20) remains independent of B.
QUENCHED MESON SPECTRUM IN EXTERNAL MAGNETIC FIELDSComputation of meson masses
The meson masses are extracted from the large Euclidean-time behavior of the correlation functions according to ⟨O(t)O†(0)⟩=A[e-mO(B)t+e-mO(B)(T-t)],where O is either the pseudoscalar density (for the pion) or the vector current (for the ρ-meson), Pa(x)=ψ¯(x)γ5τaψ(x),Vμa(x)=ψ¯(x)γμτaψ(x).The charged mesons a=± involve τ±=(τ1∓iτ2)/2, while our (connected) neutral mesons either have pure u¯u [τu=diag(1,0)] or pure d¯d [τd=diag(0,1)] quark content.
In the presence of the background magnetic field, ρ-mesons can be classified in terms of the projection of their spin onto the magnetic field, sB=szB. In terms of the vector current operator (28), the spin eigenstates with sz=±1 and 0 are given by (see, e.g., Ref. [48]) V±a=12(Vxa∓iVya)andV0a=Vza,respectively. Thus, we can compute the relevant correlation functions for sz=±1 via Cρa,±(t)=12[⟨Vxa(t)(Vxa(0))†⟩+⟨Vya(t)(Vya(0))†⟩±i(⟨Vxa(t)(Vya(0))†⟩-⟨Vya(t)(Vxa(0))†⟩)]and for sz=0 by Cρa,0(t)=⟨Vza(t)(Vza(0))†⟩.We remark that charge conjugation symmetry and parity symmetry ensure that mπ+=mπ- and mρ0+=mρ0- for the modes with spin projection sz=0 and mρ++=mρ-- and mρ-+=mρ+- for the states with nonvanishing spin projection on the B-field axis.
To determine the fit range for the function (27), we have looked for the region of the minimal value of t included in the fit, tmin, where χ2/dof is close to unity and where the extracted meson mass does not depend on the particular choice of tmin. For mesons with masses larger than the neutral pion mass mu/md(B), the correlation functions show the typical exponential decrease of the signal-to-noise ratio, so that the maximal possible value for tmin decreases. In those cases, we have less control over contaminations from excited states, and we note that these effects should be further investigated in the future.
Meson masses from lines of constant physics
We have measured the meson masses along the LCP(B) using the tuned values of κc;u and κc;d described in Sec. III D. We work with quark masses corresponding to pion masses which are above 415 MeV at B=0, so that mπL≳3 at B=0. In this case, finite size effects for the meson spectrum are expected to be small. A study of finite size effects at B≠0 is presented in Sec. IV F. We start the discussion of the results with the quark mass dependence of the meson masses on the β=6.0 ensemble, i.e., at the intermediate lattice spacing a=0.093fm.
In Fig. 4, we show the masses of neutral (left) and charged pions (right) versus B, normalized by the mass at B=0. The neutral pion masses fall off strongly before they level off at intermediate values of B, decreasing monotonously throughout. At around eB≈2GeV2, the mass of the neutral pions has decreased to about 60% to 70% of its B=0 value. The quark mass dependence of this behavior is rather mild, but the effect is enhanced toward lighter quark masses. The charged pions show the increase with the magnetic field [cf. Eq. (2) for noninteracting pions], with the tendency to undershoot the free-case prediction for larger values of B. The same behavior has also been observed for dynamical staggered quarks [9] and perturbatively using one-loop pion-nucleon interactions [49].
410.1103/PhysRevD.97.034505.f4
Results for the neutral (u¯u) pion masses (left) and the masses of π+ (right) versus B for different values of the quark masses on the ensemble with a=0.093fm. The filled curves correspond to the π+ masses in the free case, Eq. (1).
Next, we investigate ρ-mesons. Of particular relevance is the component ρ++ (or equivalently the ρ--), which is the particle conjectured to condense [17] for large values of B. First of all, we have checked that our correlators are well described by the exponential form (27) even for the largest magnetic field eB/mρ2(0)≈6.7.
Note that this is in contrast with the quenched results of Ref. [20]; however, this reference uses much smaller lattice extents. Note furthermore that Ref. [20] employed correlators in the spatial z-direction. Nevertheless, at T=0 and a magnetic field parallel to the z axis, this is expected to be equivalent to correlation functions in the t-direction.
The results for the masses of the ρ++-mesons (left) and the ρ-+ (right) are shown in Fig. 5 versus B, again normalized by the B=0 mass. The plot indicates that the mass of the ρ++ decreases monotonously with B, but it starts to level off at around eB≈2GeV2 and remains in the region of 50% to 70% of its B=0 value. In terms of the inequality (6), this behavior is expected, given the flattening of mπu and mπd at large B. In fact, the ρ++-meson mass does not saturate the bound from Eq. (6) but remains at a multiple of the pion mass larger than 1, as shown in Fig. 6 (left). The plot indicates that the ratio mρ++(B)/mπ0(B) remains constant within 20% with B and increases as the pion becomes lighter. This indicates that one moves further away from the saturation of the bound. The ρ-+ masses show the increase expected from the free case, Eq. (5) with g=2, but tend to undershoot the free-case curve for larger magnetic fields. Once more, this feature has only a very mild quark mass dependence.
510.1103/PhysRevD.97.034505.f5
Results for ρ+-meson masses with sz=+1 (left) and sz=-1 (right) versus B for different values of the quark masses on the ensemble with a=0.093fm. The curves correspond to the masses in the free case, Eq. (5).
610.1103/PhysRevD.97.034505.f6
Left: Results for the ρ+-meson masses with sz=+1 versus B, normalized to the neutral (u¯u) pion mass at finite values of B. Right: Results for ρu-meson masses with sz=±1 versus B for different values of the quark masses on the ensemble with a=0.093fm.
In Fig. 6 (right), we show the masses of the neutral ρ±u-meson with nonzero sz. Since the neutral meson has no magnetic moment, this mass remains unaffected by the magnetic field in the free case, and indeed we find that the mass is insensitive to the spin direction. However, we observe an almost linear rise with B, in this case almost completely independent of the quark mass. We suspect that disconnected diagrams might be important for this channel.
Next, we investigate the dependence of the meson masses on the lattice spacing (the final continuum extrapolation is postponed to Sec. IV D). To this end, we keep the quark mass constant by fixing the pion mass at B=0 to 415 MeV. The results for the neutral (left) and charged pions (right) versus B for different values of the lattice spacing are shown in Fig. 7. The plot indicates that in both cases only mild lattice artifacts are present. The result of the same analysis for the ρ++- and ρ-+-states is shown in Fig. 8 in the left and right panels, respectively. As for the pion masses, we observe only a mild lattice spacing dependence. Concerning the saturation of the bound (6), reducing the lattice spacing tends to drive the results away from saturation, as can be seen from Fig. 9 (left). This means that cutoff effects and quark mass effects act in the same manner in this respect. The mass of the neutral ρ±u-meson with nonzero sz, shown in Fig. 9 (right), is also mainly independent of the lattice spacing with the tendency to show noticeable effects for larger values of B.
710.1103/PhysRevD.97.034505.f7
Results for the neutral (u¯u) pion masses (left) and the masses of π+ (right) versus B for different values of the lattice spacing for mπ(0)=416MeV. The filled curves correspond to the π+ masses in the free case, Eq. (1).
810.1103/PhysRevD.97.034505.f8
Results for ρ+-meson masses with sz=+1 (left) and sz=-1 (right) versus B for different values of the lattice spacing for mπ(0)=416MeV. The curves correspond to the masses in the free case, Eq. (5).
910.1103/PhysRevD.97.034505.f9
Left: Results for the ρ+-meson masses with sz=+1 versus B, normalized to the neutral pion mass at finite values of B for mπ(0)=416MeV. Right: Results for ρu-meson masses with sz=±1 versus B for different values of the lattice spacing for an approximately constant pion mass of 416 MeV at B=0.
Pion–ρ-meson mixing
We will now discuss the mixing between pions and ρ-mesons with spin projection sz=0 along the magnetic field axis, as described in Sec. II B. Since we neglect disconnected diagrams, we will focus on the mixing between the charged meson pair π+ and ρ0+. Note that taking into account this mixing is necessary for the determination of the mass of the heavier of the pair (the ρ0+), but it does not affect the lighter state (the π+).
Due to the mixing described in Eq. (10), for B>0, the naive operators for the states π+ and ρ0+ overlap. The orthogonal basis can be determined via the diagonalization of the correlator matrix (⟨π+π+⟩⟨π+ρ0+⟩⟨ρ0+π+⟩⟨ρ0+ρ0+⟩)(t),and the true masses are the ones associated with the temporal exponential decay of the eigenvalues of this matrix in the t→∞ limit. We remark that the off-diagonal elements need to be multiplied by ±i when we translate from the Euclidean to the Minkowskian Dirac matrices (for the spatial components, these differ by a factor i so that γzE=-iγzM). To extract the eigenvalues of the matrix, we use a generalized eigenvalue problem (GEVP) (see Ref. [50] and references therein) with starting temporal extents of t0=3a, 4a, and 5a for β=5.845, 6.0, and 6.26, respectively. In the construction of the correlator matrix, we have used smeared operators at the source and sink, so the correlator matrix ought to be symmetric. In practice, this only holds within the statistical uncertainty, and we found it beneficial to stabilize the computation of eigenvalues by replacing the off-diagonal elements by (the real part of) their average.
The smaller eigenvalue of the correlator matrix corresponds to the heavier state, i.e., to the ρ0+-meson. To extract the mass of the ρ0+-meson, we perform a fit similar to that in Eq. (27). We show the associated effective masses for the two eigenvalues for the β=6.0 ensemble for eB=0.660 and 2.641GeV2 in Fig. 10. Note, that the GEVP has been setup with t0=4a and that the signal soon becomes lost in noise, so we have left out the results for the effective masses once uncertainties become overly large. For comparison, we have also plotted the result for π+ from the analysis of the pion correlator without the use of the GEVP. The results indicate good agreement for the mass of the ground state in this channel. The final masses for different values of the quark mass are shown in Fig. 11 (left), normalized to the ρ-meson mass at B=0. We see that the mass of the ρ0+ increases and is in agreement with the free-case prediction, the colored lines, even though it shows the tendency to overshoot the prediction for intermediate values of B and to undershoot for B≳3GeV2. This agreement lends support to our method and suggests that the contamination from other excited states is suppressed. The dependence of the mass on the lattice spacing is shown in Fig. 11 (right), revealing that, as before, lattice artifacts are smaller than our statistical errors. Note once more that we are working in the quenched approximation, where sea quarks are absent and the ρ-meson is a stable particle. Unlike in the case with dynamical fermions (e.g., Refs. [51,52]), we thus do not need to consider multipion states in our analysis.
1010.1103/PhysRevD.97.034505.f10
Results for the effective masses for the two different eigenvalues (blue circles for the ground state and red boxes for the first excited state) obtained from the GEVP in the channel relevant for π+ and ρ0+ mixing. The results have been obtained on the ensemble with a=0.093fm with a B=0 pion mass of 416 MeV. The colored areas are the results from a fit to the correlation function. The magenta areas are the ones for the analysis of the π+ correlation function alone, without the use of the GEVP.
1110.1103/PhysRevD.97.034505.f11
Results for ρsz=0± versus B for different values of the quark masses on the ensemble with a=0.093fm (left) and for different values of the lattice spacing at a constant B=0 pion mass of 416 MeV (right). The curves correspond to the masses in the free case, Eq. (5).
Results in the continuum limit
We proceed by performing the continuum extrapolation of the meson masses. As in the previous section, we will work at a fixed B=0 pion mass of 415 MeV. Since we employ unimproved Wilson fermions, meson masses should show lattice artifacts of O(a). Consequently, our continuum extrapolation is carried out linearly in a, and we extrapolate the meson masses for each value of B individually. The results for this continuum extrapolation are shown for two representative cases [the neutral (uu) pion masses and the ρ-+ masses] in Fig. 12. The linear extrapolation works well for most of the cases, giving χ2/dof values in the region between 0.5 and 2. However, we cannot exclude that our data still receive significant contributions by terms of O(a2). To estimate the associated systematic uncertainty, we have performed another linear continuum extrapolation including only the points at the two smallest lattice spacings and use the difference between the two extrapolations as the systematic uncertainty. In Fig. 12, the resulting uncertainty, including the systematic part, is represented by the error bars of the points on the a=0 intersect. The continuum extrapolations at high magnetic field suggest that O(a2)-effects become more pronounced as B grows. In this region, the systematic errors are potentially underestimated and should be checked in future simulations on finer lattices.
1210.1103/PhysRevD.97.034505.f12
Continuum extrapolations of the neutral (uu) pion masses (left) and the ρ-+ masses for different values of B, normalized by the mass at B=0. The filled curves represent linear fits including all data points. The error bar at a=0 includes in addition the systematic error obtained by varying the fit range.
We show the results for the continuum extrapolated meson masses in Fig. 13. The results basically confirm what we already found at finite lattice spacing in the previous section. The masses of the neutral pions (with u¯u and d¯d flavor content) decrease down to about 60%–70% of their B=0 value at B≳2GeV2, while the masses of the charged pions increase, in agreement with the energies in the free case, Eq. (1). The data for mπu can be qualitatively described by a curve of the form mπu(B)mπ(0)=1+a1(eB)21+a2(eB)2.Performing a fit to the data, including also the data for mπd, with B→B/2 following Eq. (21), we obtain a1=3.2(8)GeV-2 and a2=4.8(1.2)GeV-2. The resulting curve is shown in Fig. 13 (top), too.
1310.1103/PhysRevD.97.034505.f13
Continuum results for pion (top) and ρ-meson (bottom) masses versus B for a B=0 pion mass of 415 MeV, normalized to the respective B=0 values. Note, that the scaling of the y axes changes at the dashed horizontal line.
The mass of the ρ++, similarly to that of the neutral pion, decreases down to about 60% of its B=0 value, where it starts to level off. In fact, the ρ++ mass remains at twice the mass of the πu. The masses of the other ρ-mesons with sz=±1 increase, even those of the ρ±u/d-mesons, the B-dependence of which is, as for the neutral pions, only an indirect effect of B, due to the nonvanishing polarizability of the meson.
Comparison to trajectories with B-independent κ
In the previous sections, we have discussed the results for meson masses obtained along the LCP(B), where the bare quark mass has been tuned in order for the renormalized quark mass to remain constant with varying B. In the literature, this tuning has not been considered so far, and we will now show the problems that can arise without it. Of course, results obtained on the LCP(0) should also give the correct continuum limit; however, the lattice artifacts in this case strongly depend on B, as we will demonstrate. In this section, we will focus on the results for the neutral pion with uu flavor content, since these are the ones which are affected most by changes in the quark mass. However, other meson masses are influenced similarly by this effect.
In Fig. 14, we show the results for the masses obtained with constant κ for different values of B compared to those calculated along the LCP(B) for different lattice spacings, corresponding to a B=0 pion mass of 415 MeV. The plot indicates that the results with constant κ suffer from enormous lattice artifacts. These act in a way in which the meson masses are initially below the continuum result, while they overshoot it for strong B. This tendency originates from the nonmonotonous behavior shown in Fig. 3, where the critical mass initially becomes bigger, before it starts to decrease with B. To investigate the different lattice artifacts in more detail, we plot the masses against the lattice spacing for some representative cases in Fig. 15. The plots show that the two sets of results seem to converge toward each other as a→0. However, in most cases, a controlled (linear) continuum extrapolation is not possible for the results at constant κ. The exception is for eB=0.660GeV2, where both linear extrapolations point to the same continuum limit. Summarizing, the results for constant κ show lattice artifacts that are much larger than those of the runs along the LCP(B), which is particularly true for large values of B.
1410.1103/PhysRevD.97.034505.f14
Results for the neutral (u¯u) pion masses versus B for a pion mass of 415 MeV at B=0 and different values of the lattice spacing, normalized to their B=0 value. The open symbols correspond to results obtained with constant κ, while the filled symbols correspond to the LCP(B). We have slightly shifted the results horizontally for better visibility of the different data sets.
1510.1103/PhysRevD.97.034505.f15
Results for the masses of the neutral (u¯u) pions versus the lattice spacing in units of r0, normalized by the mass at B=0. The filled points correspond to the results obtained along the LCP(B), and the filled gray curves represent their continuum extrapolation. The open symbols denote the results obtained without B-dependent improvement. In this case, a continuum extrapolation is only possible for eB=0.66GeV2 (blue curve). The black points on the a=0 line, indicated by the vertical dashed line, are the results from the continuum extrapolations, including the estimate for systematic uncertainties. In the figures, we have slightly shifted some of the results horizontally (including the associated curves) for better visibility of the different data sets.
In Fig. 16, we compare the ρ++-meson mass results obtained from constant κ and from the LCP(B), revealing the same tendency as observed above for the pions.
1610.1103/PhysRevD.97.034505.f16
Results for the ρ++-meson masses versus B for a pion mass of 415 MeV at B=0 and different values of the lattice spacing, normalized to their B=0 value. The open symbols correspond to results obtained with constant κ, while the filled symbols show results along the LCP(B). We have slightly shifted the results horizontally for better visibility of the different data sets.
Finite size effects
In the previous sections, we have seen that the neutral pion mass decreases down to 60% of its B=0 value with increasing B. Consequently, finite size effects will potentially increase, since mπu(B)L becomes smaller. Indeed, for the strongest magnetic field, this combination drops to a value of mπu(B)L≈2, indicating that finite size effects can potentially be sizeable. To investigate this systematic uncertainty, we have generated configurations on a 48×243 lattice at β=6.0, providing a second, larger volume for this β value compared to the 48×163 lattice listed in Table I, which has been used in the previous sections. On these configurations, we have computed the meson masses for eB=1.320, 2.641, and 3.962GeV2—these magnetic fields can be realized with integer flux quanta Nb on both volumes; see Eq. (15). For the computations, we have used the κ values obtained from the LCP(B) determined on the smaller volume.
The results for the two different volumes are shown in Fig. 17 for the neutral pion and the ρ++, which are the two most important mesons for our analysis. The plot indicates that finite size effects are small. For both mesons, the results change by around 2%–3%, which is well below the typical uncertainties of more than 10% of our continuum extrapolated results. The same is true for the other meson masses. Thus, we conclude that increasing B does not enhance finite size effects and that these at present can be neglected in comparison to other systematics.
1710.1103/PhysRevD.97.034505.f17
Results for the masses of the neutral (u¯u) pions (left) and the ρ++-mesons (right) in units of the vacuum Sommer parameter r0=0.5fm. Shown are two different spatial volumes of 163 and 243 at β=6.0, corresponding to mπu(B)L≈2 and 3, respectively, for the largest value of B.
BACKGROUND FIELD METHOD AND POLARIZABILITIES WITH WILSON FERMIONSBackground field method and additive quark mass renormalization
Further interesting features of hadrons are their magnetic moments and the electric and magnetic polarizabilities. The latter describe the indirect response of a bound state to external electromagnetic fields. The background field method gives access to both the magnetic [53] and the electric (see Refs. [54,55], for instance) properties. In fact, in the past few years, a number of groups have started to measure these quantities in lattice QCD using Wilson fermions; see Refs. [23,56–66]. So far, the change of the additive renormalization of the quark mass has been neglected in these studies, each of which were carried out at a single lattice spacing. Before investigating the magnitude of this effect, let us first discuss the implications that the unwanted change in the renormalized quark mass in the presence of an external electric or magnetic field may have in this respect.
The energy of a relativistic particle (hadron) with mass m and charge q in a magnetic field aligned in the z-direction is to leading order given by (e.g., Ref. [67]) EH;n2=m2+(1+2n)|qB|-gHszqB-4πmβH|eB|2+…,n∈Z0+,where gH is the g-factor of the hadron H and βH is the polarizability, which is absent in the free field cases discussed in Sec. II. Note that our convention for the polarizabilities is in agreement with the ones from Refs. [23,65,67], which use the nonrelativistic limit, and it includes a contribution from the tensor polarizability for a particle with spin [23,65]. The ellipses stand for higher-order terms in |eB|. The magnetic moment μ of the particle is related to its g-factor in the nonrelativistic limit. The g-factor and the polarizability represent the leading- and next-to-leading-order responses of the bound state with respect to an external magnetic field. In particular, gH=-1sz[∂EH;02∂|eB|-1]|B=0andβH=18πm∂2EH;02∂|eB|2|B=0.Here, the derivative ∂ is defined as the derivative with respect to |eB| while keeping all other renormalized parameters fixed. When we measure these with Wilson fermions via the background field method, however, we obtain the full derivatives, which receive contributions from the implicit B-dependence of the quark mass, so that dEH;02d|eB|=∂EH;02∂|eB|+∑f∂EH;02∂mf∂mf∂|eB|=∂EH;02∂|eB|-∑f∂EH;02∂mf∂mc;f∂|eB|.Similarly, additional terms emerge for the second derivative as well. Analogous relations hold also for measurements of the electric polarizability. In Eq. (36), the sum is over the quark flavors f present in the measurement. In the quenched setup, the sum includes only contributions from valence quarks. For dynamical quarks exposed to the external field, however, the sum also includes effects from sea quarks.
The additional terms including derivatives of the quark masses with respect to the external field are lattice artifacts and thus vanish in the continuum limit. They are not present for actions where the quark mass is protected by a remnant of chiral symmetry. For Wilson fermions, however, those terms are present and contaminate the results for magnetic moments and polarizabilities. These contaminations can be removed by tuning the quark mass along LCP(B)s as discussed above, resulting in a cancellation of the derivatives of the quark masses with respect to the external field. Even when the tuning is only approximate, it leads to a strong reduction of the additional terms in Eq. (36) and, consequently, a reduction of lattice artifacts. We would also like to emphasize that this effect will always be present (and of the same size), no matter how small the external field is since polarizabilities are derivatives with respect to the field, Eq. (35).
Magnetic moments and polarizabilities of mesons and the impact of LCP(B)s
We proceed by investigating the strength of the effect mentioned above for polarizabilities and g-factors of mesons. Our setup is not optimal for the extraction of these quantities, since the volume in our simulations is rather small and, thus, the smallest available external field following Eq. (15) is relatively large. In fact, we could only use the lowest two values of the external field for the extraction of gH and βH. Obviously, this practice does not lead to a precision determination of these quantities, but it is sufficient for the purpose of demonstrating the improvement achieved.
Let us start with the polarizability of the neutral and the charged pion. Concerning the neutral pion, we can only determine the polarizability of pions with u¯u and d¯d flavor content individually. The results versus the lattice spacing are shown in Fig. 18 (left), both the ones obtained from LCP(B)s (filled symbols) and the ones obtained by keeping κ constant (open symbols). At finite lattice spacing, there is a visible difference between the two sets of results, and lattice artifacts appear to be larger for the results from constant κ, as expected. We did not attempt a continuum extrapolation, but it is clear that a naive linear continuum extrapolation leads to different results for the improved [LCP(B)] case and the unimproved one. However, the data indicate that a linear continuum extrapolation is not valid in the regions of lattice spacings at our disposal. A careful extrapolation to the continuum should, eventually, lead to the same result for both cases. Comparing our results from LCP(B)s to those of Ref. [67], we see that they are in the same ballpark, in particular, those for the pion with dd flavor content. Concerning the charged pion, the rather large uncertainties for the masses at finite B lead to large uncertainties for the polarizability, so we cannot draw any conclusions from them and we will not discuss these results here.
1810.1103/PhysRevD.97.034505.f18
Results for the polarizabilities of the neutral pions with u¯u and d¯d flavor content and for the magnetic moment of the (charged) ρ-meson versus the lattice spacing. The filled symbols correspond to the results obtained from LCP(B)s, and the open symbols are the results from measurements with a constant κ value.
Moving on to the ρ-mesons, the polarizability for the ρu/d with sz=±1 and the ρ0± can be extracted analogously to the ones for the neutral and the charged pions, respectively. However, once more, the uncertainties are too large to draw any definite conclusion, so we exclude these results from the discussion as well. The situation is somewhat better for the g-factor of ρ+ with sz=±1. In this case, the energies in Eq. (34) allow for the separation of the effects from the g-factor and the polarizability by building the combinations Δ=12(Eρ-+2-Eρ++2)andΣ=12(Eρ-+2+Eρ++2).Following Eq. (34), Δ is proportional to g, while Σ only contains the polarizability βρ+ as a free parameter up to O(|eB|2). Note that in the quantity Δ-effects from the additive quark mass renormalization will cancel, since they appear with equivalent prefactors in Eρ-+2 and Eρ++2. Consequently, we expect to obtain similar results for gρ+ from LCP(B)s and with constant κ. The results are shown in Fig. 18 (right). Indeed, the two sets of results show good agreement within uncertainties. The results indicate that gρ+ is close to the free-case value g=2, which is in qualitative agreement with previous lattice computations [68–70] and findings in chiral perturbation theory [71].
Since the relative size of the change in quark mass depends on the bare quark mass we have at B=0 (the change in mc;f is independent of mf), we expect the impact of the additional lattice artifacts to become increasingly important for smaller quark masses. We test this intuition by looking at the polarizabilities of the neutral pions for different values of mπ(B=0) on the lattice with β=6.0. The results are shown in Fig. 19. The increasing impact when going to smaller quark masses is clearly visible in the data. We note that lattice artifacts of this type could potentially also be responsible for the pion mass dependence of the electrical polarizability of neutral pions obtained using the connected part of the correlation function only [59,63,66]. It has been found that the polarizability becomes negative for physical quark masses, which disagrees with expectations from chiral perturbation theory [72–74].
1910.1103/PhysRevD.97.034505.f19
Results for the polarizabilities of the neutral pions with uu and dd flavor content versus mπ at B=0. The filled symbols correspond to the results obtained from LCP(B)s, and the open symbols are the results from measurements with a constant κ value. All results were obtained on the same lattice ensemble with β=6.0.
CONCLUSIONS
In this paper, we investigated the meson spectrum in QCD at zero temperature in the presence of background magnetic fields B with Wilson fermions in the quenched approximation. The new methods introduced in this paper allow us for the first time to perform the continuum limit for the Wilson spectrum at a finite value of the external field. The novelty of our approach is the introduction of a magnetic field–dependent improvement term that ensures that the renormalized quark mass remains independent of B. This requirement defines a magnetic field–dependent line of constant physics LCP(B) for the bare quark mass parameter. This B-dependent tuning of the bare quark mass is absent in the continuum (as we discuss in Appendix C) and only concerns fermion discretizations that suffer from an additive renormalization, such as Wilson quarks, where this tuning is beneficial already in the free case (see Appendix A). Note that a similar tuning along LCP(B)s should also be carried out for nonuniform external fields; see the discussion in Appendix A. The LCP(B) was determined using the lattice Ward-Takahashi identities for nonvanishing electromagnetic fields, which we derived in Appendix B.
We emphasize that the improvement only differs from the naive approach by lattice artifacts, which, however, may be large. In particular, we demonstrated that without the improvement, meson masses suffer from enormous discretization effects for strong magnetic fields. Besides the impact for the strong field limit, we also considered the effect of the improvement for derivatives with respect to the magnetic field at B=0—i.e., magnetic polarizabilities. Also, here, we found that our approach suppresses lattice artifacts and enables a flatter, controlled continuum extrapolation.
Our most important results about the spectrum involve the mass of the neutral (connected) pion, which was found to decrease monotonously as the magnetic field grows, and the mass of the lightest charged ρ-meson, which remains nonzero for the whole range of magnetic fields that we considered. The latter might be relevant for the ρ-meson condensation predicted to set in for strong magnetic fields [16,19]. Nevertheless, we mention again that our quenched analysis of the ρ-meson mass in constant background magnetic fields cannot capture the subtle details of the superconducting vacuum [17,18]. Indeed, there are examples in the literature which show that the interplay between sea and valence quarks can lead to unexpected effects [75]. Eventually, the study of the meson spectrum should thus be repeated including sea quarks. Note, however, that this demands the generation of new configurations for each value of the quark mass and the external field in the (minimal) Nf=1+1 setup at large values of Nt, rendering the study with dynamical fermions at least 2 orders of magnitude more expensive than the present quenched study. In addition to the investigation of lattice artifacts and quark mass effects, we have also checked for finite size effects and found these to be negligible compared to the other uncertainties (see Sec. IV F).
Finally, we elaborate on the connection between the B-dependence of the mass of the lightest hadron and that of the QCD transition temperature. The latter was determined using continuum extrapolated dynamical staggered quarks with physical masses in Refs. [9,10]. Remember that, according to our results, the quantity mπu(B)/mπ(0) has only mild quark mass dependence, cf. Fig. 4. It is thus sensible to compare our continuum extrapolated results to the staggered results using physical pion masses. Employing the same definition involving connected neutral pion correlators, the Wilson and staggered results
The details of the staggered simulation setup are described in Refs. [9,76]. The determination of mπu follows the strategy of Ref. [9] for the measurement of mπ+ and is based on the same lattice ensembles.
for mπu(B)/mπ(0) are shown in Fig. 20. For eB≲1GeV2, where both results are available, the comparison reveals small differences due to the different B=0 pion masses and the inclusion of sea effects in the staggered case.
2010.1103/PhysRevD.97.034505.f20
Comparison of Wilson and staggered continuum extrapolated results for the lightest pion mass (normalized by its B=0 value). Also shown is the QCD transition temperature defined via the inflection point of the average light quark condensate—also in units of its B=0 value. Note the difference in the B=0 pion mass, which is 415 MeV for the Wilson results and 135 MeV for the staggered data.
Regarding the QCD transition temperature, we considered the parametrization of the crossover transition temperature Tc(B) defined using the inflection point of the average light quark condensate [10]. According to Fig. 20, the B-dependence of the curves looks qualitatively very similar—with an initial reduction followed by a saturation to around 60% to 70% of the B=0 value. This supports the picture sketched in the Introduction, where we compared the finite temperature QCD transition to the “melting” of the lightest hadron state. Note that the quantitative difference for intermediate values of B could originate from effects due to the charged pions, which, in this region, are potentially light enough to have a direct influence on the transition temperature. We remark that, while the inclusion of charged sea quarks appears to have a marginal impact on the T=0 pion masses, it is known to make a drastic difference for Tc[8]. A possible explanation is that the transition temperature is driven by bulk effects encoded in the configurations (i.e., the action used to generate them), so an action leading to the correct pion mass is essential to observe the correct features of the transition. For meson masses, measurements are most affected by the properties of valence quarks, which reflects itself in the fact that the quenched spectrum reproduces the QCD spectrum up to about 10% accuracy (see, e.g., Refs. [24,25]). Nevertheless, it would be interesting to understand this subtle difference between measurements of properties of the phase transition and the hadron spectrum in more detail.
ACKNOWLEDGMENTS
This research was funded by the DFG (Emmy Noether Programme EN 1064/2-1 and SFB/TRR 55). The majority of the simulations was performed on the iDataCool cluster of the Institute for Theoretical Physics at the University of Regensburg. The authors are grateful for the useful correspondence with Pavel Buividovich, Maxim Chernodub, Davide Giusti, Rainer Sommer, and Arata Yamamoto. We thank Max Theilig for a careful reading of the Appendixes.
ADDITIVE MASS RENORMALIZATION IN THE FREE CASE
To demonstrate the presence of a B-dependent additive quark mass renormalization, it is instructive to look at the free case. We have already discussed this setting in Ref. [21]; here, we repeat the main findings for the sake of completeness.
The massless Wilson Dirac operator may be written schematically as DW=D+ar·Δ/2, where D is the naive discretization of the (anti-Hermitian) continuum Dirac operator and the second part is the Wilson term. While D describes spin-1/2 particles, Δ is the discretization of the Klein-Gordon operator and, thus, describes scalars. The eigenvalues of the two operators for B>0 can be read off from Eqs. (4) and (1) for massless particles. Therefore (at zero temperature, where the lowest Matsubara frequency is zero), the lowest eigenvalue of D vanishes, while that of the Wilson term equals ar·|qB|/2. Since the two operators commute on the subspace spanned by the lowest eigenmode, the real part of the lowest eigenvalue of DW increases linearly with B. When discretized on the lattice, this conclusion continues to hold for small values of a2qB; see Fig. 21. In fact, for this check, it suffices to diagonalize the Wilson operator on a two-dimensional x, y-plane, since the eigenmodes factorize and form plane waves in the z- and t-directions.
2110.1103/PhysRevD.97.034505.f21
The lowest eigenvalue (in lattice units) of the Wilson Dirac operator with r=1 against the magnetic field (also in lattice units), calculated on a two-dimensional 322 lattice.
This increase of the lowest eigenvalue is equivalent to an additive shift in the quark mass, implying a B-dependent additive renormalization. The critical mass (or the critical hopping parameter κc) is reduced by the magnetic field as amc;f(a,B)=amc(a)-a2|qfB|/2,κc;f-1(a,B)=κc-1(a)-a2|qfB|.Thus, in the free theory, a fixed renormalized quark mass, see Eq. (16), is achieved if the quark mass parameter is shifted by the same amount (note that in the free case, Zm=1). We mention that this additive shift in the spectrum is only present for fermion formulations which break the full chiral symmetry explicitly and does not appear, for example, in the staggered formulation [77].
The effect of this B-dependent tuning can be demonstrated using the free “pion masses,” i.e., the energies associated to the leading decay of the pseudoscalar correlation functions, corresponding to the energies of quark-antiquark states with imposed pion quantum numbers. On the one hand, for neutral correlation functions, both quarks have magnetic moments parallel to B and thus are in the ground state. The associated energy should be E(B)=mu+md. On the other hand, for the charged correlation function, one of the quarks (that with the smaller absolute charge) is forced to have its magnetic moment antiparallel to B. The energy in this case should be E(B)=mu+md2+2|qdB|.
Our numerical results for the energies at fixed κ=0.124 are shown in Fig. 22 (left). The energy of the neutral pion increases with the magnetic field, indicating the unphysical increase of the quark mass by the amount a2|qfB|/2. Tuning the hopping parameters along the trajectory (A1) instead, the neutral pion mass remains constant as it should; see Fig. 22 (right). Our results are also in good agreement with the expectation for the charged “pion.”
2210.1103/PhysRevD.97.034505.f22
Results for the energies associated with free correlation functions in the pseudoscalar channel without (left) and with (right) tuning of κ with B along the trajectory (A1). The colored dashed lines in the left plot show the analytic expectations Ef(B)=E(0)+a|qfB|/2 for energies associated with neutral correlation functions, and the solid line is the expectation for charged pions.
The shift in the bare quark mass basically follows from the magnetic field dependence of the lowest energy state of a noninteracting charged scalar particle. For a homogeneous background magnetic field, this dependence can be found analytically and is given by Eq. (1). We also checked numerically what happens if the magnetic background field is inhomogeneous. In particular, we considered an oscillatory field B(x)∝sin(2πkx/Lx) with k∈Z and a half-half field B(x)∝(-1)xdiv(Lx/2). In both cases, the bare quark mass was observed to increase with growing magnetic field—although not linearly, as for the homogeneous background, but quadratically in a2qB.
AXIAL WARD-TAKAHASHI IDENTITIES FOR QCD+QED
In this Appendix, we will derive axial and vector WIs for QCD in the presence of electromagnetic interactions with two light quark flavors, i.e., for QCD+QED. We will start by (re)deriving the continuum WIs in Appendix B 1, which, to our knowledge, have so far only been discussed in Ref. [37], before we discuss the WIs for (unimproved) Wilson quarks on the lattice in Appendix B 2. Note that an account of the continuum WIs has already appeared in Ref. [21]. The WIs for O(a)-improved Wilson quarks can be derived along the same lines and will be discussed briefly in Appendix B 3. Note that the axial Ward identity for the domain wall fermion discretization has been derived in Ref. [38] and recently for twisted mass fermions at maximal twist in Ref. [39].
We thank Davide Giusti for drawing our attention to this.
Continuum Ward-Takahashi identities
In the continuum, the Euclidean fermion action including u and d quarks in the presence of QED interactions is given by SFconti=∫d4xψ¯(x)(D+M)ψ(x)with the covariant derivative D=DQCD+Dem,DQCD=γμ(∂μ+iGμ(x))1,Dem=iQγμAμ(x);M the matrix of bare quark masses in the continuum, given by Eq. (12), and the charge matrix Q=e(τ06+τ32)withτ0≡1.We have denoted the gluon fields by Gμ(x) and the electromagnetic (photon) field by Aμ(x).
Ward-Takahashi identities can be derived by varying the expectation value of some test operator O under transformations of the fermionic variables in the path integral of the form δψ(x)=i[αVj(x)τj2]ψ(x),δψ¯(x)=-iψ¯(x)[αVj(x)τj2],corresponding to infinitesimal transformations under SUV(2), and δψ(x)=i[αAj(x)τj2γ5]ψ(x),δψ¯(x)=iψ¯(x)[αAj(x)τj2γ5]for infinitesimal transformations under SUA(2). Here, τj is any of the three Pauli matrices for j=1, 2, 3 and the unit matrix for j=0, corresponding to isosinglet transformations. Using these transformations, one obtains (with X=A, V) ∂∂αXj(x)⟨O⟩=⟨∂O∂αXj(x)⟩-⟨∂SF∂αXj(x)O⟩=-⟨∂SF∂αXj(x)O⟩=0in the flavor nonsinglet case, where the second equality is valid as long as the operator O has no support at the point x. In the singlet case, one also has to consider the variation of the measure.
Using Eqs. (B1) and (B2), we obtain ∂SFconti∂αXj(x)=∫d4x∂∂αXj(x)(ψ¯(x)(DQCD+M)ψ(x))+∫d4x∂∂αXj(x)(ψ¯(x)Demψ(x)).The result for the first term is well known for the two types of transformations (see, e.g., Ref. [78]), ∂∂αVj(x)(ψ¯(x)(DQCD+M)ψ(x))=-i∂μ(JV)μj(x)+iψ¯(x)[M,τj/2]ψ(x)and∂∂αAj(x)(ψ¯(x)(DQCD+M)ψ(x))=-i∂μ(JA)μj(x)+iψ¯(x){M,τj/2}γ5ψ(x),where the curly brackets denote the anticommutator, so we only need to compute the additional terms stemming from the second term in Eq. (B7). Here, (JV)μj(x) is the (continuum) vector current from Eq. (28), including the Pauli matrix τj/2, and (JA)μj(x) is the local axial vector current
We use the notation JV and JA instead of V and A in this section to avoid confusion with the electromagnetic vector potential Aμ and to distinguish between point-split and local currents in the lattice regularization.
(JA)μj(x)=ψ¯(x)γμγ5τj2ψ(x).
Considering vector transformations, the additional terms are given by ∂∂αVj(x)(ψ¯(x)Demψ(x))=-e4ψ¯(x)γμAμ(x)[τ3,τj]ψ(x)=-iε3jkeψ¯(x)Aμ(x)γμτk2ψ(x),so together with [M,τj/2]=14(mu-md)[τ3,τj]=iε3jk(mu-md)τk2,we arrive at the continuum vector WI ∂μ(JV)μj(x)=iε3jk(mu-md)Sk(x)-ε3jkeAμ(x)(JV)μk(x).Here, εijk is the totally antisymmetric tensor for i, j, k=1, 2, 3 and zero if either of i, j, k=0, and we have introduced the scalar density Sj(x)=ψ¯(x)τj2ψ(x).
With axial vector transformations, we obtain ∂∂αAj(x)(ψ¯(x)Demψ(x))=-e4ψ¯(x)Aμ(x){γμτ3,γ5τj}ψ(x)=-iε3jkeψ¯(x)Aμ(x)γμγ5τk2ψ(x)and {M,τj/2}=(mu+md)τj2+12(mu-md)δj3,so we get for the continuum axial WI ∂μ(JA)μj(x)=(mu+md)Pj(x)+δj3(mu-md)P0(x)-ε3jkeAμ(x)(JA)μk(x)-e216π2ϵαβμνFαβ(x)Fμν(x)Tr[τjQ2]-δj0116π2ϵαβμνTr[Hαβ(x)Hμν(x)].Here, we have introduced the pseudoscalar densities Pj(x)=ψ¯(x)γ5τj2ψ(x)andP0(x)=ψ¯(x)γ512ψ(x);Hμν and Fμν are the gluonic and electromagnetic field strength tensors, respectively; and Tr denotes the trace over flavor and color indices. The terms in the last line of Eq. (B17) are the ones associated with the Jacobian of the transformation from Eq. (B5)[79], known as the axial anomaly. For QCD+QED, there are two such terms, associated with the topological charge operators for the gluonic and the electromagnetic fields, respectively. Note that for external magnetic fields (and no electric fields), the electromagnetic anomaly is absent.
Ward-Takahashi identities for Wilson fermions
We will now derive similar identities for (unimproved) Wilson fermions using the fermion action from Eq. (11). In particular, we need to compute (∂SF)/(∂αXj(x)) for the transformations from Eqs. (B8) and (B9). The basic idea to separate the resulting equations in parts which are already known and new parts is the same as in the continuum case. In practice, however, the separation is a bit less intuitive since the electromagnetic field enters the action via the electromagnetic link variables uμ(x); see Eq. (13). The separation can be done by bringing the Pauli matrices from the variation of the ψ¯ variables to the right-hand side of the link variables.
Note that this choice corresponds to a particular convention for the definition of point-split bilinear operators. An alternative convention includes the Pauli matrices to the left of the link variables. Both definitions lead to equivalent results as long as one strictly follows one and the same convention analytically and numerically.
When the links uμ(x) are absent, this is possible since the gluonic link variables commute with matrices in flavor space. In the presence of uμ(x), new terms will appear due to the nonvanishing commutator [uμ,τj]=iε3jk[exp(i23eaAμ(x))-exp(-i13eaAμ(x))]τk≡iε3jk(δu)μ(x)τk,where the last equation defines (δu)μ(x).
Let us start again with the vector transformations. Using the results for Wilson fermions from Ref. [80], we arrive at the vector WI a4∑μ∇xμ(J˜V)μj(x)=iε3jka4(mu0-md0)Sk(x)-irε3jka32∑μ[ψ¯(x)UμG(x)(δu)μ(x)τk2ψ(x+μ^)+ψ¯(x)(UG)μ†(x-μ^)(δu)μ†(x-μ^)τk2ψ(x-μ^)]}(1)-iε3jka32∑μ[-ψ¯(x)γμUμG(x)(δu)μ(x)τk2ψ(x+μ^)+ψ¯(x)γμ(UG)μ†(x-μ^)(δu)μ†(x-μ^)τk2ψ(x-μ^)].}(2)In Eq. (B20), we have introduced the asymmetric lattice derivative a∇xμf(x)=f(x)-f(x-μ^), and J˜V is the point-split (conserved) vector current (J˜V)μj(x)=12[ψ¯(x)(γμ-r)Uμ(x)τj2ψ(x+μ^)+ψ¯(x+μ^)(γμ+r)Uμ†(x)τj2ψ(x)].The first two terms are the standard terms appearing in the vector Ward identity, while the terms denoted by (1) and (2) appear due to Eq. (B19).
Similarly, we obtain for the axial WI a4∑μ∇xμ(J˜A)μj(x)=a4(mu0+md0)Pj(x)+δj3(mu0-md0)P0(x)+a4rXj(x)-irε3jka32∑μ[ψ¯(x)UμG(x)(δu)μ(x)γ5τk2ψ(x+μ^)+ψ¯(x)(UG)μ†(x-μ^)(δu)μ†(x-μ^)γ5τk2ψ(x-μ^)]}(1)-iε3jka32∑μ[-ψ¯(x)γ5γμUμG(x)(δu)μ(x)τk2ψ(x+μ^)+ψ¯(x)γ5γμ(UG)μ†(x-μ^)(δu)μ†(x-μ^)τk2ψ(x-μ^)],}(2)where Xj=-12a∑μ{[ψ¯(x)Uμ(x)γ5τj2ψ(x+μ^)+ψ¯(x+μ^)Uμ(x)γ5τj2ψ(x)]+[x→x-μ^]-4ψ¯(x)γ5τj2ψ(x)}is the term associated with the variation of the Wilson term in the standard QCD axial WI (see Ref. [81]) and J˜A is the point-split axial vector current (J˜A)μj(x)=12[ψ¯(x)γμγ5Uμ(x)τj2ψ(x+μ^)+ψ¯(x+μ^)γμγ5Uμ†(x)τj2ψ(x)].Note, that in QCD [without terms (1) and (2)], the presence of Xj is responsible for the additive quark mass renormalization since Xj mixes with the pseudoscalar density. In Eq. (B22), we have neglected the terms associated with the anomaly [see Eq. (B17)].
Let us briefly discuss the properties and the nature of the terms denoted by (1) and (2). These two types of terms are due to the presence of the electromagnetic interactions in the QCD+QED Lagrangian and vanish, as expected, in the pure QCD case, i.e., when Aμ(x)=0 [since (δu)μ(x)=0]. Since the additional terms are proportional to ε3jk, they only contribute to “charged” WIs, i.e., WIs which include combinations of τ1 and τ2. Neutral WIs, i.e., such that are a combination of WIs with τ3 and 1, are unaffected by the presence of the QED interactions.
This is true up to the anomalous terms for the isosinglet axial WI.
The terms denoted by (1) arise from the Wilson term (consequently, they are multiplied by a factor of r), and, as Xj, they are proportional to a dimension-5 operator and vanish in the continuum. However, as Xj, these terms also mix with the scalar (for the vector WI) and pseudoscalar (for the scalar WI) density, respectively, and thus lead to an additional additive quark mass renormalization. In particular, the presence of such a term in the vector WI indicates that the quark mass difference does not vanish when the bare quark masses coincide, meaning that u and d quark masses renormalize with different additive terms. The terms denoted by (2) are operators of dimension 4 and survive the continuum limit. In fact, they exactly resemble the additional terms in the continuum WIs, Eqs. (B13) and (B17), proportional to Aμ(x).
To compute current quark masses, terms of type (2) should be included explicitly in their definition, since these terms are also present in the continuum WI, while terms of type (1) should be left out (similar to Xj). For an example, consider the definition of the sum of u and d current quark masses via the axial WI evaluated for the matrix τ+ from Sec. IV A. Using the pseudoscalar density P+ from Eq. (B18) as the operator O in Eq. (B6), a suitable definition for the sum of current quark masses is given by a(m˜u(B)+m˜d(B))=∂0⟨(J˜A)0+(x0)P+(0)⟩-⟨(2)P+(0)⟩⟨P+(x0)P+(0)⟩,where (2) refers to the insertion of the associated operators from Eq. (B22).
Note that the electromagnetic link variables that need to be included in the point-split currents depend on the particular flavor matrix for which the WI is evaluated. As an example, let us once more consider the WIs for the matrix τ+. Writing uμ(x)=(uμu(x)00uμd(x))withuμf=exp(iqfaAμ(x)),we obtain ψ¯(x)uμ(x)τ+ψ(x)=d¯(x)uμd(x)u(x).
Ward-Takahashi identities for O(a)-improved Wilson fermions
For O(a)-improvement in QCD+QED, we need to include two Sheikholeslami-Wohlert [82] terms in the action, SFSW=-i2cSWa5ψ¯(x)σμνHμν(x)ψ(x)andSFem=-i2cema5ψ¯(x)σμνCμν(x)ψ(x),where Hμν is (a suitable discretization of) the gluonic field strength tensor and Cμν=QFμν is the electromagnetic one, multiplied by the electric charge matrix Q. The terms in the WIs which result from these additional terms depend on the particular choice of discretization. Most commonly used is the clover discretization of the field strength tensor, Hμν(x)=18a2[Qμν(x)-Qνμ(x)]withQμν(x)=Uμ,νG(x)+Uν,-μG(x)+U-μ,-νG(x)+U-ν,μG(x),and similarly for Cμν with UG replaced by u. Here, U±μ,±νG(x) [and similarly u±μ,±ν(x) for Cμν] denotes the multiplication of links around a plaquette in the ±μ, ±ν-direction, starting from point x.
For vector transformations, we only get a contribution from the electromagnetic clover term. In particular, we get terms which are proportional to the commutator [u±μ,±ν(x),τj]=iε3jk[exp(i23eaA±μ,±ν(x))-exp(-i13eaA±μ,±ν(x))]τk≡iε3jk(δu)±μ,±ν(x)τk,where A±μ,±ν(x) denotes the sum of the phase factors appearing in the plaquette u±μ,±ν(x) and the last equation defines (δu)±μ,±ν(x) in analogy to (δu)μ(x). When we define (δC)μν similar to Hμν from Eq. (B29) with U±μ,±νG(x)→(δu)±μ,±ν(x), the additional term on the right-hand side of Eq. (B20) is given by a4∑μ∇xμ(J˜V)μj(x)=(B20)+12ε3jkcema5ψ¯(x)σμν(δC)μντk2ψ(x).
For axial vector transformations, it is the anticommutator of τj and the clover term which is relevant for the WI. For the gluonic clover term, the anticommutator is nonvanishing, leading to the known additional term ∂SFSW∂αAj(x)=cSWa5ψ¯(x)σμνHμνγ5τj2ψ(x).For the electromagnetic clover term, we obtain terms proportional to {u±μ,±ν(x),τj}=[exp(iea23A±μ,±ν(x))+exp(-iea13A±μ,±ν(x))]τj+δ3j(δu)±μ,±ν(x)≡(Σu)±μ,±ν(x)τj+δ3j(δu)±μ,±ν(x),where the last equation defines (Σu)±μ,±ν. We now define (ΣC)μν analogously to (δC)μν with (δu)±μ,±ν→(Σu)±μ,±ν, so that we obtain the extension of Eq. (B22) in the form a4∑μ∇xμ(J˜A)μj(x)=(B22)-icSWa5ψ¯(x)σμνHμνγ5τj2ψ(x)-i2cema5[ψ¯(x)σμν(ΣC)μν(x)γ5τj2ψ(x)+12δ3jψ¯(x)σμν(δC)μν(x)γ5ψ(x)].
MULTIPLICATIVE MASS RENORMALIZATION IN MAGNETIC FIELDS
In this Appendix, we demonstrate that the QCD mass renormalization constant is independent of the background magnetic field. On general grounds, it is expected that the ultraviolet divergent renormalization constants of the theory do not depend on physical parameters like the magnetic field. Here, we show this using one-loop perturbation theory in continuum QCD.
For a similar calculation at nonzero temperature, see Ref. [83].
The mass renormalization constant is obtained from the fermion self-energy, where the double line represents the free fermion propagator in the background magnetic field and the curly line represents the gluon propagator. The SU(3) generators are denoted by ta.
In the Schwinger proper time representation [84], the Feynman-gauge gluon propagator reads iD˜μνab(k)=-δμνδab∫0∞dte-k2t,while the B-dependent fermion propagator is taken from Ref. [85] (analytically continued to imaginary proper times and Wick rotated to Euclidean space-time with Euclidean Dirac matrices), iS˜(p)=∫0∞dse-(m2+p∥2)s-p⊥2/(qB)tanh(qBs)[ip+m-(pxγy-pyγx)tanh(qBs)][1-iγxγytanh(qBs)],where p⊥2=px2+py2 and p∥2=pz2+pt2 and throughout we will assume that qB>0.
In fact, Eq. (C3) is the Fourier transform of the translational invariant part of the propagator and excludes the Schwinger phase. Inverse Fourier transforming Eq. (C1) and multiplying by the Schwinger phase gives the coordinate space representation of both sides [32,85]. Accordingly, the perpendicular p⊥ are not real quantum numbers but referred to as pseudomomenta.
The propagator can be expanded in powers of the magnetic field. The consecutive orders read iS˜0(p)=∫0∞dse-(m2+p2)s[ip+m],iS˜1(p)=-iqBγxγy∫0∞dsse-(m2+p2)s[i(pzγz+ptγt)+m],iS˜2(p)=(qB)2∫0∞dss2e-(m2+p2)s[(ip+m)p⊥2s3-i(pxγx+pyγy)]and can be inserted, order by order in the magnetic field, into Σ˜. In color space, we only have the trivial factor ∑atata=CF1 with CF=4/3. The Dirac structure can be simplified with conventional techniques.
The B=0 contribution to the self-energy reads, after integration over p, iΣ˜0(k)=g28π2CF∫0∞ds∫0∞dt1(s+t)2[2m-ikts+t]e-m2s-k2sts+t.This can be analytically integrated in t. The result at low s behaves like 1/s, such that an UV cutoff is necessary at the lower limit of integration 1/Λ2<s, and the integral diverges as logΛ. To obtain the mass renormalization constant, we need to evaluate Σ˜ at the pole of S˜0. This amounts to the replacement k→-im (compare to Ref. [79], and remember that we are working with Euclidean metric). We arrive at iΣ˜0(k→-im)=3g216π2CFm[logΛ2m2-γE+56],which is the well-known one-loop mass renormalization constant in QCD. Equation (C6) coincides with Schwinger’s result for QED (replacing g→q and CF→1).
We proceed with the O(qB)-term. After integration over p, we get iΣ˜1(k)=g2eB8π2CFγxγy(kzγz+ktγt)∫0∞ds∫0∞dtst(s+t)3e-m2s-k2sts+t,which can be integrated directly over t. The result behaves as s0 at low s and thus involves no ultraviolet divergences. Analytically integrating over s results in iΣ˜1(k)=g2eB8π2CFγxγy(k0γ0+k3γ3)1k2[1+m2k2logm2k2+m2],which is related to the anomalous magnetic moment of the quark [84].
Finally, we consider the O((qB)2)-term. Also, here, we are able to perform the p-integral, followed by the t-integral. The result is linear in s at low s and is thus again ultraviolet finite. In fact, since higher powers of the magnetic field always appear with higher powers of s [see Eq. (C4)], the ultraviolet region s≈0 is getting more and more suppressed as we consider higher and higher orders in the magnetic field. Altogether, this shows that the B-dependent part of the quark self-energy is, to all orders in B, finite. Therefore, the mass renormalization constant contains no B-dependent divergences, and one is free to use a B-independent renormalization scheme to define the renormalized mass.
1D. Kharzeev, Parity violation in hot QCD: Why it can happen, and how to look for it, 633, 260 (2006).PYLBAJ0370-269310.1016/j.physletb.2005.11.0752D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, The effects of topological charge change in heavy ion collisions: ‘Event by event P and CP violation,’A803, 227 (2008).NUPBBO0550-321310.1016/j.nuclphysa.2008.02.2983R. C. Duncan and C. Thompson, Formation of very strongly magnetized neutron stars—Implications for gamma-ray bursts, 392, L9 (1992).ASJOAB1538-435710.1086/1864134T. Vachaspati, Magnetic fields from cosmological phase transitions, 265, 258 (1991).PYLBAJ0370-269310.1016/0370-2693(91)90051-Q5C. Bonati, M. D’Elia, and A. Rucci, Heavy quarkonia in strong magnetic fields, 92, 054014 (2015).PRVDAQ1550-799810.1103/PhysRevD.92.0540146K. Fukushima and Y. Hidaka, Magnetic Shift of the Chemical Freeze-out and Electric Charge Fluctuations, 117, 102301 (2016).PRLTAO0031-900710.1103/PhysRevLett.117.1023017A. Broderick, M. Prakash, and J. M. Lattimer, The equation of state of neutron star matter in strong magnetic fields, 537, 351 (2000).ASJOAB1538-435710.1086/3090108J. O. Andersen, W. R. Naylor, and A. Tranberg, Phase diagram of QCD in a magnetic field: A review, 88, 025001 (2016).RMPHAT0034-686110.1103/RevModPhys.88.0250019G. S. Bali, F. Bruckmann, G. Endrődi, Z. Fodor, S. D. Katz, S. Krieg, A. Schäfer, and K. K. Szabó, The QCD phase diagram for external magnetic fields, 02 (2012) 044.JHEPFG1029-847910.1007/JHEP02(2012)04410G. Endrődi, Critical point in the QCD phase diagram for extremely strong background magnetic fields, 07 (2015) 173.JHEPFG1029-847910.1007/JHEP07(2015)17311V. G. Bornyakov, P. V. Buividovich, N. Cundy, O. A. Kochetkov, and A. Schäfer, Deconfinement transition in two-flavor lattice QCD with dynamical overlap fermions in an external magnetic field, 90, 034501 (2014).PRVDAQ1550-799810.1103/PhysRevD.90.03450112E. M. Ilgenfritz, M. Muller-Preussker, B. Petersson, and A. Schreiber, Magnetic catalysis (and inverse catalysis) at finite temperature in two-color lattice QCD, 89, 054512 (2014).PRVDAQ1550-799810.1103/PhysRevD.89.05451213J. O. Andersen, Thermal pions in a magnetic background, 86, 025020 (2012).PRVDAQ1550-799810.1103/PhysRevD.86.02502014Y. Hidaka and A. Yamamoto, Charged vector mesons in a strong magnetic field, 87, 094502 (2013).PRVDAQ1550-799810.1103/PhysRevD.87.09450215E. V. Luschevskaya, O. E. Solovjeva, O. A. Kochetkov, and O. V. Teryaev, Magnetic polarizabilities of light mesons in SU(3) lattice gauge theory, B898, 627 (2015).NUPBBO0550-321310.1016/j.nuclphysb.2015.07.02316S. Schramm, B. Müller, and A. J. Schramm, Quark-antiquark condensates in strong magnetic fields, 07, 973 (1992).MPLAEQ0217-732310.1142/S021773239200086017M. N. Chernodub, Spontaneous Electromagnetic Superconductivity of Vacuum in Strong Magnetic Field: Evidence from the Nambu–Jona-Lasinio Model, 106, 142003 (2011).PRLTAO0031-900710.1103/PhysRevLett.106.14200318M. N. Chernodub, Comment on “Charged vector mesons in a strong magnetic field”, 89, 018501 (2014).PRVDAQ1550-799810.1103/PhysRevD.89.01850119M. N. Chernodub, Superconductivity of QCD vacuum in strong magnetic field, 82, 085011 (2010).PRVDAQ1550-799810.1103/PhysRevD.82.08501120V. V. Braguta, P. V. Buividovich, M. N. Chernodub, A. Yu. Kotov, and M. I. Polikarpov, Electromagnetic superconductivity of vacuum induced by strong magnetic field: numerical evidence in lattice gauge theory, 718, 667 (2012).PYLBAJ0370-269310.1016/j.physletb.2012.10.08121B. B. Brandt, G. Bali, G. Endrődi, and B. Gläßle, QCD spectroscopy and quark mass renormalisation in external magnetic fields with Wilson fermions, , LATTICE2015 (2016) 265.22E. V. Luschevskaya and O. V. Larina, The ρ and A mesons in a strong abelian magnetic field in SU(2) lattice gauge theory, B884, 1 (2014).NUPBBO0550-321310.1016/j.nuclphysb.2014.04.00323S. R. Beane, E. Chang, S. Cohen, W. Detmold, H. W. Lin, K. Orginos, A. Parreno, M. J. Savage, and B. C. Tiburzi, Magnetic Moments of Light Nuclei from Lattice Quantum Chromodynamics, 113, 252001 (2014).PRLTAO0031-900710.1103/PhysRevLett.113.25200124S. Aoki (CP-PACS Collaboration), Quenched Light Hadron Spectrum, 84, 238 (2000).PRLTAO0031-900710.1103/PhysRevLett.84.23825S. Aoki (JLQCD Collaboration), Light hadron spectroscopy with two flavors of O(a) improved dynamical quarks, 68, 054502 (2003).PRVDAQ0556-282110.1103/PhysRevD.68.05450226C. Bonati, M. D’Elia, M. Mariti, M. Mesiti, F. Negro, A. Rucci, and F. Sanfilippo, Magnetic field effects on the static quark potential at zero and finite temperature, 94, 094007 (2016).PRVDAQ2470-001010.1103/PhysRevD.94.09400727I. A. Shushpanov and A. V. Smilga, Quark condensate in a magnetic field, 402, 351 (1997).PYLBAJ0370-269310.1016/S0370-2693(97)00441-328N. O. Agasian and I. A. Shushpanov, Gell-Mann-Oakes-Renner relation in a magnetic field at finite temperature, 10 (2001) 006.JHEPFG1029-847910.1088/1126-6708/2001/10/00629J. O. Andersen, Chiral perturbation theory in a magnetic background—finite-temperature effects, 10 (2012) 005.JHEPFG1029-847910.1007/JHEP10(2012)00530F. Bruckmann, G. Endrődi, M. Giordano, S. D. Katz, T. G. Kovács, F. Pittler, and J. Wellnhofer, Landau levels in lattice QCD, , LATTICE2016 (2016) 048.31F. Bruckmann, G. Endrődi, M. Giordano, S. D. Katz, T. G. Kovacs, F. Pittler, and J. Wellnhofer, Landau levels in QCD, 96, 074506 (2017).PRVDAQ2470-001010.1103/PhysRevD.96.07450632V. A. Miransky and I. A. Shovkovy, Quantum field theory in a magnetic field: From quantum chromodynamics to graphene and Dirac semimetals, 576, 1 (2015).PRPLCM0370-157310.1016/j.physrep.2015.02.00333M. H. Al-Hashimi and U. J. Wiese, Discrete accidental symmetry for a particle in a constant magnetic field on a torus, 324, 343 (2009).APNYA60003-491610.1016/j.aop.2008.07.00634Sz. Borsanyi, Ab initio calculation of the neutron-proton mass difference, 347, 1452 (2015).SCIEAS0036-807510.1126/science.125705035S. R. Sharpe, Quenched chiral logarithms, 46, 3146 (1992).PRVDAQ0556-282110.1103/PhysRevD.46.314636C. W. Bernard and M. F. L. Golterman, Chiral perturbation theory for the quenched approximation of QCD, 46, 853 (1992).PRVDAQ0556-282110.1103/PhysRevD.46.85337S. Scherer, Introduction to chiral perturbation theory, 27, 277 (2003).ANUPBZ0065-297010.1007/b10051938T. Blum, T. Doi, M. Hayakawa, T. Izubuchi, and N. Yamada, Determination of light quark masses from the electromagnetic splitting of pseudoscalar meson masses computed with two flavors of domain wall fermions, 76, 114508 (2007).PRVDAQ1550-799810.1103/PhysRevD.76.11450839D. Giusti, V. Lubicz, G. Martinelli, F. Sanfilippo, and S. Simula, Strange and charm HVP contributions to the muon (g-2) including QED corrections with twisted-mass fermions, 10 (2017) 157.JHEPFG1029-847910.1007/JHEP10(2017)15740P. Fritzsch, J. Heitger, and N. Tantalo, Non-perturbative improvement of quark mass renormalization in two-flavor lattice QCD, 08 (2010) 074.JHEPFG1029-847910.1007/JHEP08(2010)07441K. G. Wilson, Confinement of quarks, 10, 2445 (1974).PRVDAQ0556-282110.1103/PhysRevD.10.244542R. G. Edwards and B. Joo (SciDAC, LHPC, and UKQCD Collaborations), The Chroma software system for lattice QCD, 140, 832 (2005).NPBSE70920-563210.1016/j.nuclphysbps.2004.11.25443S. Necco and R. Sommer, The N(f)=0 heavy quark potential from short to intermediate distances, B622, 328 (2002).NUPBBO0550-321310.1016/S0550-3213(01)00582-X44S. Güsken, A study of smearing techniques for hadron correlation functions, 17, 361 (1990).NPBSE70920-563210.1016/0920-5632(90)90273-W45M. Falcioni, M. L. Paciello, G. Parisi, and B. Taglienti, Again on SU(3) glueball mass, B251, 624 (1985).NUPBBO0550-321310.1016/0550-3213(85)90280-946M. Lüscher, Solution of the Dirac equation in lattice QCD using a domain decomposition method, 156, 209 (2004).CPHCBZ0010-465510.1016/S0010-4655(03)00486-747M Lüscher, Local coherence and deflation of the low quark modes in lattice QCD, 07 (2007) 081.JHEPFG1029-847910.1088/1126-6708/2007/07/08148B. de Wit and J. Smith, , North-Holland Personal Library (North-Holland, Amsterdam, 1986), Vol. 1.49G. Colucci, E. S. Fraga, and A. Sedrakian, Chiral pions in a magnetic background, 728, 19 (2014).PYLBAJ0370-269310.1016/j.physletb.2013.11.02850B. Blossier, M. Della Morte, G. von Hippel, T. Mendes, and R. Sommer, On the generalized eigenvalue method for energies and matrix elements in lattice field theory, 04 (2009) 094.JHEPFG1029-847910.1088/1126-6708/2009/04/09451C. McNeile and C. Michael (UKQCD Collaboration), Hadronic decay of a vector meson from the lattice, 556, 177 (2003).PYLBAJ0370-269310.1016/S0370-2693(03)00130-852G. S. Bali, S. Collins, A. Cox, G. Donald, M. Göckeler, C. B. Lang, and A. Schäfer (RQCD Collaboration), ρ and K* resonances on the lattice at nearly physical quark masses and Nf=2, 93, 054509 (2016).PRVDAQ2470-001010.1103/PhysRevD.93.05450953G. Martinelli, G. Parisi, R. Petronzio, and F. Rapuano, The proton and neutron magnetic moments in lattice QCD, 116B, 434 (1982).PYLBAJ0370-269310.1016/0370-2693(82)90162-954H. R. Fiebig, W. Wilcox, and R. M. Woloshyn, A study of hadron electric polarizability in quenched lattice QCD, B324, 47 (1989).NUPBBO0550-321310.1016/0550-3213(89)90180-655J. C. Christensen, W. Wilcox, F. X. Lee, and L.-m. Zhou, Electric polarizability of neutral hadrons from lattice QCD, 72, 034503 (2005).PRVDAQ1550-799810.1103/PhysRevD.72.03450356F. X. Lee, L. Zhou, W. Wilcox, and J. C. Christensen, Magnetic polarizability of hadrons from lattice QCD in the background field method, 73, 034503 (2006).PRVDAQ1550-799810.1103/PhysRevD.73.03450357W. Detmold, B. C. Tiburzi, and A. Walker-Loud, Electromagnetic and spin polarizabilities in lattice QCD, 73, 114505 (2006).PRVDAQ1550-799810.1103/PhysRevD.73.11450558F. X. Lee, S. Moerschbacher, and W. Wilcox, Magnetic moments of vector, axial, and tensor mesons in lattice QCD, 78, 094502 (2008).PRVDAQ1550-799810.1103/PhysRevD.78.09450259W. Detmold, B. C. Tiburzi, and A. Walker-Loud, Extracting electric polarizabilities from lattice QCD, 79, 094505 (2009).PRVDAQ1550-799810.1103/PhysRevD.79.09450560W. Detmold, B. C. Tiburzi, and A. Walker-Loud, Extracting nucleon magnetic moments and electric polarizabilities from lattice QCD in background electric fields, 81, 054502 (2010).PRVDAQ1550-799810.1103/PhysRevD.81.05450261T. Primer, W. Kamleh, D. Leinweber, and M. Burkardt, Magnetic properties of the nucleon in a uniform background field, 89, 034508 (2014).PRVDAQ1550-799810.1103/PhysRevD.89.03450862J. M. M. Hall, D. B. Leinweber, and R. D. Young, Finite-volume and partial quenching effects in the magnetic polarizability of the neutron, 89, 054511 (2014).PRVDAQ1550-799810.1103/PhysRevD.89.05451163M. Lujan, A. Alexandru, W. Freeman, and F. Lee, Electric polarizability of neutral hadrons from dynamical lattice QCD ensembles, 89, 074506 (2014).PRVDAQ1550-799810.1103/PhysRevD.89.07450664W. Freeman, A. Alexandru, M. Lujan, and F. X. Lee, Sea quark contributions to the electric polarizability of hadrons, 90, 054507 (2014).PRVDAQ1550-799810.1103/PhysRevD.90.05450765E. Chang, W. Detmold, K. Orginos, A. Parreno, M. J. Savage, B. C. Tiburzi, and S. R. Beane (NPLQCD Collaboration), Magnetic structure of light nuclei from lattice QCD, 92, 114502 (2015).PRVDAQ1550-799810.1103/PhysRevD.92.11450266M. Lujan, A. Alexandru, W. Freeman, and F. X. Lee, Finite volume effects on the electric polarizability of neutral hadrons in lattice QCD, 94, 074506 (2016).PRVDAQ2470-001010.1103/PhysRevD.94.07450667E. V. Luschevskaya, O. E. Solovjeva, and O. V. Teryaev, Magnetic polarizability of pion, 761, 393 (2016).PYLBAJ0370-269310.1016/j.physletb.2016.08.05468J. N. Hedditch, W. Kamleh, B. G. Lasscock, D. B. Leinweber, A. G. Williams, and J. M. Zanotti, Pseudoscalar and vector meson form-factors from lattice QCD, 75, 094504 (2007).PRVDAQ1550-799810.1103/PhysRevD.75.09450469M. Gurtler (QCDSF Collaboration), Vector meson electromagnetic form factors, , LATTICE2008 (2008) 051.70E. V. Luschevskaya, O. E. Solovjeva, and O. V. Teryaev, Lattice Stern-Gerlach experiment, 09 (2017) 142.JHEPFG1029-847910.1007/JHEP09(2017)14271D. Djukanovic, E. Epelbaum, J. Gegelia, and U. G. Meissner, The magnetic moment of the ρ-meson, 730, 115 (2014).PYLBAJ0370-269310.1016/j.physletb.2014.01.00172J. Bijnens and F. Cornet, Two pion production in photon-photon collisions, B296, 557 (1988).NUPBBO0550-321310.1016/0550-3213(88)90032-673J. F. Donoghue, B. R. Holstein, and Y. C. Lin, The reaction gamma Gamma→pi0pi0 and Chiral loops, 37, 2423 (1988).PRVDAQ0556-282110.1103/PhysRevD.37.242374J. Hu, F.-J. Jiang, and B. C. Tiburzi, Pion polarizabilities and volume effects in lattice QCD, 77, 014502 (2008).PRVDAQ1550-799810.1103/PhysRevD.77.01450275F. Bruckmann, G. Endrődi, and T. G. Kovács, Inverse magnetic catalysis and the Polyakov loop, 04 (2013) 112.JHEPFG1029-847910.1007/JHEP04(2013)11276S. Borsanyi, G. Endrődi, Z. Fodor, A. Jakovac, S. D. Katz, S. Krieg, C. Ratti, and K. K. Szabo, The QCD equation of state with dynamical quarks, 11 (2010) 077.JHEPFG1029-847910.1007/JHEP11(2010)07777G. Endrődi, QCD in magnetic fields: from Hofstadter’s butterfly to the phase diagram, , LATTICE2014 (2014) 018.78A. Vladikas, in Modern Perspectives in Lattice QCD: Quantum Field Theory and High Performance Computing, (Oxford University Press, Oxford, 2011), p. 161.79M. E. A. Peskin and D. V. A. Schroeder, , The Advanced Book Program (Addison-Wesley, Reading, USA, 1995).80L. H. Karsten and J. Smit, Lattice fermions: Species doubling, chiral invariance, and the triangle anomaly, B183, 103 (1981).NUPBBO0550-321310.1016/0550-3213(81)90549-681M. Bochicchio, L. Maiani, G. Martinelli, G. C. Rossi, and M. Testa, Chiral symmetry on the lattice with Wilson fermions, B262, 331 (1985).NUPBBO0550-321310.1016/0550-3213(85)90290-182B. Sheikholeslami and R. Wohlert, Improved continuum limit lattice action for QCD with Wilson fermions, B259, 572 (1985).NUPBBO0550-321310.1016/0550-3213(85)90002-183A. Bandyopadhyay, N. Haque, and M. G. Mustafa, The pressure of a weakly magnetized deconfined QCD matter within one-loop Hard-Thermal-Loop perturbation theory, arXiv:1702.02875.84J. S. Schwinger, On gauge invariance and vacuum polarization, 82, 664 (1951).PHRVAO0031-899X10.1103/PhysRev.82.66485I. A. Shovkovy, Magnetic catalysis: A review, 871, 13 (2013).LNPHA40075-845010.1007/978-3-642-37305-3