Magnetic polarizability of pion
Abstract
We explore the energy dependence of mesons off the background Abelian magnetic field on the base of quenched SU(3) lattice gauge theory and calculate the magnetic dipole polarizability of charged and neutral pions for various lattice volumes and lattice spacings. The contribution of the magnetic hyperpolarizability to the neutral pion energy has been also found.
keywords:
lattice QCD, SU(3) gluodynamics, magnetic field, magnetic polarizability, pseudoscalar meson1 Introduction
Quantum Chromodynamics in strong magnetic fields is a promising topic for research. Magnetic fields of hadronic scale could exist in the Early Universe Hector:00 () and could be formed in cosmic objects like magnetars and neutron stars. They can also be achieved in terrestrial laboratories (RHIC, LHC, FAIR, NICA) Skokov:2009 (). The external electromagnetic fields can be utilized as a probe of QCD properties, the recent progress obtained on this way in lattice gauge theories is discussed in Massimo:2015 (). The energy levels of hadrons in external magnetic field can be usefull for the calculation of the cross sections Savage:2015 (). The energies of mesons at the nonzero magnetic field were calculated in various phenomenological approaches Simonov:2013 (); Liu:2015 (); Taya:2015 (); Kawaguchi:2016 (); Hattori:2016 (), within the QCD sum rules Cho:2015 (); Gubler:2015 () and in the lattice gauge theories Bali:2015 (); Luschevskaya:2015a (); Luschevskaya:2015b ().
Background magnetic Abelian fields also enable to calculate the magnetic polarizabilities of hadrons. In order to get the dipole magnetic polarizability and the hyperpolarizability we measure the energy of a meson as a function of the magnetic field. The magnetic polarizabilities are important physical characteristics describing the distibution of quark currents inside a hadron in an external field. For the first time the concept of polarizability for the nuclear matter was used by A.B. Migdal in the analysis of the scattering of low energy gamma quanta by atomic nuclei Migdal (). For hadrons the notion of the polarizability was discussed in papers klein (); baldin ().
There were some discrepancies between the experimental obtained value of the magnetic and the electric polarizabilities of the charged mesons and some theoretical predictions based on the chiral perturbation theory Gasser:2005 (); Aleksejevs:2013 ().
Measurement of the electrical and the magnetic polarizabilities of mesons was performed on the spectrometer SIGMAAJAX in Serpukhov, on the electron synchrotron Pakhra (LPI) in Moscow, on the MarkII detector at SLAC, at COMPASS (CERN) and other experiments.
According to the obtained data from these experiments, the value of the polarizability of the charged mesons is positive. The most precise value of the charged pion electric polarisability has been obtained experimentally by the COMPASS experiment Adolph:2015 () under the assumption . Comparison with the successful predictions of the experiments and the chiral perturbation theory Gasser:2005 (); Aleksejevs:2013 () is interesting for fundamental science.
In this work we consider a behaviour of the ground state energy of pions in the strong magnetic fields on the base of the SU(3) lattice gauge theory. The details of the calculations are briefly sketched in section 2. We discuss the dipole magnetic polarizabilities of and mesons in sections 3 and 4 accordingly. The magnetic hyperpolarizability of the neutral pion is calculated in section 5.
2 Details of calculations
2.1 The improved gauge action
Ensemble  
8.20  0.115  245  
8.10  0.125  285  

8.20  0.115  200  
8.30  0.105  235  
8.45  0.095  195  
8.60  0.084  180  
8.20  0.115  275  

For the generation of quenched lattice configurations we used the tadpole improved LüscherWeisz action Luscher:1985 (), which reduces the ultraviolet lattice artifacts. The action has the form
(1) 
where is the plaquette (denoted by pl) or 12 rectangular loop term (rt), is the input tadpole factor computed at zero temperature Bornyakov:2005 (). Our simulations have been carried out on the symmetrical lattices. The parameters of the lattice ensembles and number of the configurations are listed in the Table 1.
2.2 Fermionic spectrum
We solve the Dirac equation numerically
(2) 
and found the eigenfunctions and the eigenvalues for a test quark in an external gauge field .
For this goal we use the massive overlap operator Neuberger:1997 (). It has the following form
(3) 
where is the quark mass, is the parameter in our calculations, is the massless overlap Dirac operator, which preserves chiral invariance even at the finite lattice spacing . It may be written as
(4) 
where is the WilsonDirac operator with the negative mass term , is the Wilson hopping term, is the Hermitian WilsonDirac operator.
We construct the polynomial approximation for the function
(5) 
This approximation should be valid on the entire spectrum of matrix . Since
(6) 
and , then . The MinMax polynomial approximation is used for the function at , where .
The polynomial of a degree can be the best approximation, if it minimizes the maximal relative error
(7) 
where . The polynom is represented by the series
(8) 
where are the Chebyshev polynomials defined in the range . The detailed description of the algorithm used for the calculation of can be found in Giusti:2003 (). The resulting polynomial of the matrix has the same eigenfunctions as the original matrix. The eigenvalues can be found from the eigenfunctions using the formula for a some operator Neff:2001 (). From the polynomial approximation for the sign function we get the approximation for the overlap Dirac operator. We find its eigenfunctions and eigenvalues, which are used for the calculation of the propagators and the correlators.
As we consider pure lattice gauge theory, the Abelian magnetic field is introduced only into the overlap Dirac operator. In the symmetric gauge the magnetic field parallel to ’z’ axis has the form
(9) 
The total gauge field is the sum of the nonAbelian gluonic field and the Abelian field of magnetic photons
(10) 
where are the colour indices, are the Lorentz indices. Quark fields obey periodic boundary conditions in space and antiperiodic boundary conditions in time. In order to match (10) with the periodic boundary conditions we apply the additional xdependent boundary twist for fermions AlHashimi:2009 ().
In the finite lattice volume the magnetic flux trough any twodimensional face of the hypecube is quantized. So the magnetic field value is
(11) 
where is the elementary charge and is the numbers of lattice sites in spatial directions.
2.3 Calculation of correlation functions
To observe the ground state energy for a meson we construct the interpolating operator creating the state with the corresponding quantum numbers. In case the pseudoscalar charged meson the interpolating operator is described by the equations
(12) 
The interpolating operator for the neutral pion is
(13) 
It should be mentioned that in Eucledean space . We are interested in the 2point lattice correlation function of the interpolating operators
(14) 
where are Dirac gamma matrices, and are lattice coordinates. The spatial lattice coordinate and are the numbers of lattice sites in time direction.
The correlation function can be represented as the sum of connected and disconnected parts
(15) 
where is the Dirac propagator. For the isovector currents the disconnected part of the correlation function has to be zero due to cancelation of and quarks contributions, and we verify this on the lattice (see Sec. 4). The massive Dirac propagator is calculated using the lowest WilsonDirac eigenmodes
(16) 
We perform a discrete Fourier transformation of (15) numerically and set because we are interested in the ground state energy. To obtain the masses we expand the correlation function to the exponential series
(17) 
where and are the corresponding Hilbert space operators.
For the large the main contribution to the correlator (17) comes from the term corresponding to the ground state. Due to the periodic boundary conditions the correlator in the leading order has the following form
(18) 
where is a constant, is the ground state energy, is the lattice spacing. We found the ground state energies fitting the lattice correlators by the equation (18) at .
3 Dipole polarizability of meson
The energy levels of a free charged pointlike particle in a constant Abelian magnetic field parallel to z axis are described by the following equation:
(19) 
where is the particle momentum along ’z’ axis, is the principal quantum number, is the particle charge, is its spin projection, is the dimensionless quantity, that characterizes the magnetic moment of the particle and is the energy at zero magnetic field (the mass).
For the charged pion in the ground state in the rest frame it is necessary to utilize , , and . In the strong external magnetic field the charged pion is not a pointlike particle anymore and its internal structure can be described by the magnetic polarizabilities. If we take into account the dipole magnetic polarizability and the first order magnetic hyperpolarizability , then in the relativistic case the pion energy squared has the following form
(20) 
The energy of a charged pion was calculated from the correlation function
(21) 
The energies squared are shown in Fig.1 for the lattice volumes , , , lattice spacings , , , quark masses and . We fit our lattice data by formula (20) at , where , and are the fit parameters. The fitting curves are also depicted in Fig.1. At the magnetic fields we observe a linear dependence of the energy squared versus the magnetic field and it increases with the field for the all sets of lattice data. At higher fields the nonlinear terms in the magnetic field contribute to the pion energy.
Fitting the lattice data we do not consider the term with the second order magnetic hyperpolarizability () in (20), because from the data fits it follows that its relative contribution is small in comparison with the terms proportional to and .
In Table 2 we show the values of the magnetic dipole polarizability and hyperpolarizability for several lattice data sets. The best values obtained are for the lattice volume , lattice spacing and for the lattice volume and lattice spacing .
In two loops of the chiral perturbation theory it was predicted that the magnetic polarizability of equals to Aleksejevs:2013 (); Ivanov:2015 (), that is close to our result. The polarizabilities for the and lattice at quark mass and are not represented in Table 2 because of the poor relative accuracy. It is a very subtle effect but the work in this direction is carried out, including the smaller bare quark masses.
In 2015 the COMPASS collaboration at CERN has investigated the pion Compton scattering Adolph:2015 (). They have found the pion electric polarizability equal to under the assumption of the ChPT, which is true in the exact chiral limit Donoghue:1989 (). We have observed agreement of this value with our lattice results. This is also in accordance with the earlier analysis of MARK II group data of the cross section of the process Boyer:1990 () which was made in Babusci:1992 () and gave the value .
The Serpukhov group have found the value exploring the radiative pion nucleon scattering in 1983 Antipov:1983 (). Without using the relation they obtained .
The analysis of the experimental data on helicity amplitudes of the reaction gave the value Filkov:2006 ().
4 Dipole polarizability of meson
We calculate the energy of neutral pion from the correlation function
(22) 
It is taken into account that the and quarks interact with the external magnetic field differently because of their nonequal charges.
We have checked on the lattice that the disconnected part of the correlator (15) is zero within the error range. In Fig.3 it is shown for three values of the magnetic field. In Fig.4 the connected correlators are represented with the fits (18) for comparison. In the case of exact isospin symmetry the disconnected part has to be zero so in what follows we neglect it at all.
Fig.2 shows the ground state energy of the meson depending on the field value squared for ensembles at the quark mass and at the quark masses , (see Table 1). If the neutral pion would consist of one type of quark, (), then its energy decreases slower (faster) than for the real pion, see Fig.2.
For the pion each lattice data set is fitted at using the function
(23) 
where we determine and as the fit parameters. is the energy at zero magnetic field and is the magnetic polarizability of neutral pion which is presented in Table 3. We observe the linear energy dependence on the magnetic field squared at small values of the field.
In Fig.5 we represent the values of the meson as a function of the quark mass for ensemble , the values for ensembles , at quark mass are also depicted.
The value of diminishes with the quark mass and extrapolation to the chiral limit gives the number for the lattice volume and lattice spacing . It is a bit higher than the value of magnetic polarizability obtained in our previous work Luschevskaya:2015a (), because we considered meson consisting of one type of quarks . In the framework of the chiral perturbation theory the value of the is equal to in one (two) loops Gasser:2005 ().
Our result coincides in sign with the prediction of the ChPT and differs in value, and there can be various reasons for that, including the contributions of higher chiral loops. In turn, from the lattice side, the dynamical quark loops may be relevant.
5 Magnetic hyperpolarizability of neutral pion
The terms of the higher degrees of the magnetic field give contribution to the energy of the pseudoscalar meson at very high magnetic field. Fig.6 shows the ground state energy depending on the magnetic field squared at . The lattice data fits are described by the formula
(24) 
where is the magnetic hyperpolarizability and is the hyperpolarizability of higher order, , and are the fit parameters. The term proportional to is parity forbidden: the physical interpretation is that can’t decay to three (magnetic) photons.
The decay is allowed, so the term presents in formula (24). The values are negative for the all sets of the lattice data and showed in Table 4.
We found the value at lattice spacing and at for the lattice volume , and quark mass .
At the lowest quark mass we obtained using for fitting formula (24).
Note also that the energy shows the qualitative tendency of flattening at large magnetic field so that it should not turn to zero and cause the phase transition, in complete similarity to meson case studied earlier Luschevskaya:2015b (). The decrease at very large fields is likely to be compensated by higher hyperpolarizabilities. Still, the quantitative analysis of this behaviour requires more investigations.
6 Conclusions
We have calculate the magnetic dipole polarizability of the charged pion for several lattice spacings, lattice volumes and bare lattice quark mass . For the lattice volume and the lattice spacing we obtain the value . At finest lattice with the lattice spacing and the volume we get . We found the agreement of this result with the observation of COMPASS collaboration.
We have explored the energy dependence of neutral pion off the external Abelian magnetic field and have found its dipole magnetic polarizability and hyperpolarizability. The dipole magnetic polarizability depends on the lattice quark mass and the chiral extrapolation was performed. In the chiral limit we have obtained the magnetic dipole polarizability of the is equal to for the lattice volume and the lattice spacing , that is close to the prediction of the ChPT Aleksejevs:2013 ().
The contribution of the hyperpolarizability have been also revealed to the energy at considered strong magnetic field. This is very tiny effect and for the all sets . For the finest lattice with spacing we found . We have observed qualitative behaviour of energy which is not favouring the appearance of tachyonic mode and phase transition.
7 Acknowledgements
O.T. is indebted to A. Guskov and M. Ivanov for useful discussions.
The authors are grateful to FAIRITEP supercomputer center where these numerical calculations were performed. This work is completely supported by a grant from the Russian Science Foundation (project number 161210059).
/d.o.f.  
/d.o.f.  

/d.o.f.  
References
 D. Grasso and H.R. Rubinstein, Phys. Rept. 348, 163 (2001), arXiv: astroph/0009061.
 V. Skokov, A. Illarionov and V. Toneev, Int. J. Mod. Phys. A 24, 5925 (2009), arXiv: 0907.1396[nuclth].
 Massimo D’Elia, PoS LATTICE2014 020 (2015), IFUPTH201503, arXiv:1502.06047 [heplat].
 S.R. Beane et al., NPLQCD Collaboration, Phys.Rev.Lett. 115 no.13, 132001 (2015), arXiv:1505.02422 [heplat].
 M.A. Andreichikov, B.O. Kerbikov, V.D. Orlovsky, Yu.A. Simonov, Phys.Rev. D 87 (2013) 094029 (2013), arXiv: 1304.2533 [hepph].
 H. Liu, L. Yu and M. Huang, Phys. Rev. D 91, 014017 (2015), arXiv:1408.1318 [hepph].
 H. Taya, Phys.Rev. D 92, no.1, 014038 (2015), arXiv:1412.6877 [hepph].
 M. Kawaguchi, S. Matsuzaki, Phys.Rev. D 93, no.12, 125027 (2016), arXiv:1511.06990.
 K. Hattori, T. Kojo and N. Su, Nucl.Phys. A 951 1 (2016), arXiv:1512.07361 [hepph].
 S. Cho, K. Hattori, S.H. Lee, K. Morita and Sho Ozaki, Phys.Rev. D 91, no.4, 045025 (2015), arXiv:1411.7675 [hepph].
 Ph. Gubler et. al., Phys.Rev. D 93, no.5, 054026 (2016), arXiv:1512.08864.
 G. Bali, B.B. Brandt, G. Endrodi and B. Glaessle, arXiv:1510.03899 [heplat].
 E.V. Luschevskaya, O.E. Solovjeva, O.E. Kochetkov and O.V. Teryaev, Nucl. Phys. B 898, 627 (2015).
 E.V. Luschevskaya, O.E. Kochetkov, O.V. Teryaev and O.E. Solovjeva, JETP Letters 101, Issue 10, 674 (2015).
 A.B. Migdal, J.Phys. USSR 8, 331 (1944)
 A. Klein, Phys. Rev. 99, 998 (1955)
 A.M. Baldin, Nucl. Phys.18, 310 (1960)
 J. Gasser, M.A. Ivanov and M.E. Sainio, Nucl. Phys. B 728, 31 (2005).
 A. Aleksejevs and S. Barkanova, Nucl. Phys. Proc. Suppl. 245, 17 (2013), arXiv:1309.3313 [hepph].
 C. Adolph et al. (COMPASS Collaboration), Phys. Rev. Lett. 114, 062002 (2015).
 M. Lüscher and P. Weisz, Commun. Math. Phys. 97 59 (1985).
 H. Neuberger, Phys. Lett. B 417, 141 (1998), arXiv: heplat/9707022.
 M.A. Ivanov, Int.J.Mod.Phys.Conf.Ser. 39 1560104 (2015), ePrint: arXiv:1509.05225.
 J.F. Donoghue, B.R. Holstein, Phys.Rev.D 40, 2378 (1989).
 J. Boyer et al., Phys.Rev.D 42, 1350 (1990)
 D. Babusci, S. Bellucci, G. Giordano, G. Matone, A. M. Sandorfi, M.A. Moinester, Phys.Lett.B 277, 158 (1992)
 Yu. M. Antipov et al., Phys.Lett.B 121, 445 (1983).
 V.G. Bornyakov, E.M. Ilgenfritz, and M. MüllerPreussker, Phys. Rev. D 72, 054511 (2005), heplat/0507021.
 L. Giusti, C. Hoelbling, M. Lüscher, and H. Wittig, Comp. Phys. Commun. 153, 31 (2003).
 H. Neff, N. Eicker, Th. Lippert, J.W. Negele, and K. Schilling, Phys. Rev. D 64, 114509 (2001), arXiv: heplat/0106016.
 M.H. AlHashimi and U.J. Wiese, Annals Phys. 324 343 (2009), arXiv: 0807.0630 [quantph].
 L.V.Fil’kov, V.L.Kashevarov, Phys.Rev.C 73 035210 (2006), arXiv:nuclth/0512047.