Mesons in ultraintense magnetic field: an evaded collapse
Abstract
Spectra of mesons are investigated in the framework of the Hamiltonian obtained from the relativistic path integral in external homogeneous magnetic field. The spectra of all 12 spinisospin wave states generated by  and mesons with different spin projections, are studied analytically as functions of the field strength. Three types of behavior with characteristic splittings are found. The results are in agreement with recent lattice calculations.
Mesons in ultraintense magnetic field: an evaded collapse
B.O. Kerbikov^{a}^{a}aemail: borisk@itep.ru \@andauthorM.A. Andreichikov^{b}^{b}bemail: andreichicov@mail.ru \@andlastauthorYu.A. Simonov^{c}^{c}cemail: simonov@itep.ru
1A.I. Alikhanov Institute for Theoretical and Experimental Physics,
Moscow 117218, Russia
\@textsuperscript2Lebedev Physical Institute, Moscow 119991, Russia
\@textsuperscript3Moscow Institute of Physics and Technology, Dolgoprudny 141700,
Moscow Region, Russia

Abstract.
1 Introduction
The interest to the behavior of quarks, hadrons and atoms in strong magnetic field (MF) has been very high during the last decade. The outbrake of the research activity in this field was inspired by the fact that MF up to ^{1}^{1}1 We use the relativistic system of units . Then 1 GeV is generated during the early stages of peripheral heavyion collisions at RHIC and LHC. The field about four orders of magnitudes less is anticipated to operate in magnetars. The immediate question is what happens to the mass and the wave function of a meson embedded in such a strong MF. The answer to this question has been searched for in various approaches (see [1] for a list of references) including lattice simulations. In the present work the problem is investigated in the framework of the relativistic path integral Hamiltonian (PIH) formalism [2, 3, 4]. For pion this method has to be supplemented by the elements of chiral dynamics [5]. The analytical results will be compared with the lattice calculations presented recently in [1]. Before getting involved with the details of calculations it makes sense to relate the MF strength to some characteristic physical parameter which defines the spectrum of quarkantiquark meson states. From the textbooks we know that for the hydrogen atom the critical, or the socalled “atomic field”, is . This value corresponds to the situation when the magnetic, or Landau radius is equal to the Bohr radius. The QCD coupling constant , the meson radius at is determined by the QCD string tension GeV [4]. It is therefore natural to define for the hadron spectra the critical MF as which yields fm. This value is approximately equal or smaller than the typical hadron size.
The determination of the hadron spectrum in MF is not an easy task. The first problem is to separate the centerofmass (c.m.) motion. For the neutral nonrelativistic system in MF this can be done (with some qualifications) making use of the Pseudomomentum [6, 7, 8, 9]. This approach was extended to the relativistic sector within the PIH framework in [4]. For a charged meson Pseudomementum method is applicable only for an unphysical model of a meson with two equally charged quarks [4]. In this contribution we present the results on the meson spectrum in MF both within the Pseudomomentum approach and in a new analytical method of constituents separation (CS). It will be argued that its accuracy is within 15% for ultrastrong MF , and within 20% for . The method allows to study neutral and charged mesons in the same way. The results for the neutral mesons will be obtained both in Pseudomomentum and CS approaches. In this way the accuracy of the CS will be tested.
The most important question we have to answer is whether the meson spectrum in MF is bounded from below. In other words, does the meson mass reaches zero value at some MF strength. We shall point out the two dynamical mechanisms that might have led to such a collapse and explain why this does not happen. The paper is organized as follows. In section 2 the relativistic Hamiltonian based on path integral FeynmanFockSchwinger representation is written down and the spectral problem is formulated. In section 3 we discuss the possible types of meson mass trajectories in MF. Section 4 contains the analysis of perturbative corrections and potential reasons for the collapse of meson state in MF. In section 5 we present the main results in comparison with lattice calculations.
2 The relativistic Hamiltonian and the spectral problem
To find the meson masses in MF we use the path integral Hamiltonian (PIH) method based on the FeynmanFockSchwinger representation [2, 3, 4]. It allows with the help of Wilson loop to treat the interaction of quarks with external Abelian and nonAbelian fields in a gaugeinvariant way. As it was shown in [2, 3, 4, 10] the quarkantiquark spectral problem in MF in PIH formalism is reduced to the bound states problem for the relativistic Hamiltonian which includes all the nonperturbative dynamics
(1) 
Here is the th quark dynamical mass, or the einbein variable [2, 3, 4]. The MF is convenient to take in the symmetric gauge since this gauge allows to define the angular momentum projection of each quark as a quantum number. The next step is to perform minimization with respect to which yields the physical spectrum
(2) 
(3) 
The total meson mass is a sum of the nonperturbative (dynamical) one obtained from (1)(3) and the firstorder perturbative contributions
(4) 
where is the onegluon exchange potential, and and are the spinspin and self energy contributions.
For the neutral hadrons (mesons and baryons) the eigenvalue problem (2) admits the exact solution which is obtained by the separation of the c.m. motion. To this end the pseudomomentum operator is introduced [6, 7, 8, 9]
(5) 
In MF the pseudomomentum takes the role of the mechanical momentum, commutes with the Hamiltonian, and is therefore a constant of motion. Physically, is conserved since it takes into account the Lorentz force acting on particles in MF.
In [1] we have proposed a more general approach which allows to investigate the mass spectra of both neutral and charged mesons. This is the constituent separation (CS) method. The c.m. position is fixed at the origin and an effective string tension is attributed to each quark
(6) 
In this picture quarks may be considered as quasiindependent of the nonperturbative part of the interaction.
3 Meson trajectories in strong MF
In strong MF it is convenient to use for the spin degrees of freedom in (1) a basis in which the operator is diagonal. The four vectors forming this basis are . One can easily obtain three types of asymptotic meson trajectories at . The character of the trajectory is determined by the signs of the quark charges and the spin directions. According to the terminology adopted in atomic physics the trajectory is called lowfield seeking (LFS) if the energy decreases as the MF decreases. The state which at is MFindependent may be called zerofield seeking (ZFS). The types of asymptotic trajectories are
(7) 
The trajectory LFS2 exhibits stronger MF dependence than LFS1. From (7) it follows that the family which contains and quarks is distributed among the three above classes in the following way: belongs to ZFS, and belong to LFS1, rests in LFS2. The same situation up to the sign change holds for . The states and contain and components. The charges of and are different and this results in an additional double splitting.
4 Perturbative corrections
As it was shown in [10], the color Coulomb interaction in presence of the MF contains a potential danger of a collapse for ZFS states. This reminds the ”falltothecenter” phenomenon in the hydrogen atom enbedded in strong MF. The onegluon exchange matrix element has a negative sign, its absolute value grows with , and as shown on Fig.1 ”falltothecenter” takes place at . To prevent the collapse the screening by the quark loops at the lowest Landau level (LLL) was introduced in [10].
As one can see from Fig.1, the color Coulomb collapse is really evaded. We remind that in superstrong MF radiative corrections screen the Coulomb potential in the hydrogen atom thus leading to the freezing of the ground state energy at the value [11, 12]. It is interesting to note that asymptotically at the matrix element vanishes.
Another threat of a collapse comes from the hyperfine spinspin interaction . In the firstorder perturbation theory in PIH formalism it corresponds to the colormagnetic interaction of the form
(8) 
In strong MF the ground state wave function acquires the form of an ellipsoid elongated in the direction of MF. At the transverse and longitudinal radii are , . This means that the focusing of the wave function at the origin and a divergent factor in the matrix element of . Not that the problem of singularity due to function interaction exists without MF as well. It is cured by smearing function [13, 14]. In PIH formalism there is a natural cutoff parameter . It corresponds to the correlation length of the stochastic vacuum gluonic field. The function is replaced by
(9) 
In this way the ”falltothecenter”is prevented for all ZFS states except for the meson. The trajectory is stabilized if one takes its chiral degrees of freedom into account. We also note that in [3] a general theorem was proven according to which the eigenvalues of the relativistic Hamiltonian in MF are positive. The explicit account of pion chiral dynamics [5] confirm this result. The main point is that GMOR relations remain valid for neutral pions in arbitrarystrong MF, while charged pions loose their chiral properties at .
5 Results and conclusions
Below we present the results of our analytic calculations in comparison with the recent lattice results from [1]. In Fig.2 the meson mass evolution in MF is shown.
In Fig.3 and Fig.4 results for and are exposed. One should keep in mind that and components give rise to their own trajectories. The growing trajectories belong to the LFS2 class and the splitting is equal to . In Fig.5 we present the mass evolution of chiral and nonchiral in comparison with the lattice data. The chiral effects provide the decrease of the mass to its physical value at .
In this work we have evaluated the trajectories of and meson masses as functions of the external MF. The meson quark content and pion chiral dynamics were thoroughly taken into account. The most interesting problem was whether the mass remains finite in arbitrary strong MF. The collapse might have happened either due to color Coulomb interaction, or due to spinspin potential proportional to function. We have shown that in both cases there are physical reasons why collapse is evaded. The analytic calculations of meson mass trajectories give the results which are in agreement with recent lattice simulations.
Acknowledgement
The authors are grateful to E.V.Luschevskaya and O.E.Solovjeva for kindly providings with their lattice results and to O.V.Teryaev for discussions. The authors are supported by the Russian Science Foundation grant number 161210414.
References
 [1] M.A.Andreichikov, B.O.Kerbikov, E.V.Luschevskaya, Yu.A.Simonov, O.E.Solovjeva, arXiv:1610.06887 (2016).
 [2] Yu.A.Simonov, Phys.Rev. D88, 025028 (2013), arXiv:1303.4952.
 [3] Yu.A.Simonov, Phys.Rev. D88, 053004 (2013), arXiv:1304.0365.
 [4] M.A.Andreichikov, B.O.Kerbikov, V.D.Orlovsky, Yu.A.Simonov, Phys.Rev. D87, 094029 (2013), arXiv:1304.2533.
 [5] V.D.Orlovsky, Yu.A.Simonov, JHEP 1309, 135 (2013), arXiv:1306.2322.
 [6] W.E.Lamb, Phys.Rev. 85, 59 (1952).
 [7] L.P.Gor’kov, I.E.Dzyaloshinskii, Soviet Phys. JETP 26, 449 (1968).
 [8] J.E.Avron, I.W.Herbst, B.Simon, Ann.Phys. (NY) 114, 431 (1978).
 [9] H.Grotsch, R.A.Hegstrom, Phys.Rev. A4, 59 (1976).
 [10] M.A.Andreichikov, V.D.Orlovsky, Yu.A.Simonov, Phys.Rev.Lett. 110, 162002 (2013).
 [11] A.E.Shabad, V.V.Usov, Phys.Rev.Lett. 98 180403 (2007).
 [12] M.I.Vysotsky, S.I.Godunov, Phys.Uspekhi 57, 194 (2014).
 [13] B.O.Kerbikov, M.I.Polikarpov, I.V.Shevchenko, Nucl.Phys. B331, 19 (1990).
 [14] P.Brovetto, V.Maxia, M.Salis, Nuovo Cim. 14, 795 (1992).
 [15] G.Bali, B.B.Brandt, G.Endrodi, B.Glassle, 33rd International Symposium on Lattice field theory, 1418 Jul 2015, Kobe, Japan, arXiv:1510.03899.
 [16] Y.Hidaka, A.Yamamoto, Phys.Rev. D87, 094502 (2013).