Erratum to: Magnetic dominance of axion electrodynamics: photon capture effect and anisotropy of Coulomb potential

A Correction to this paper has been published: 10.1140/epjc/s10052-021-09046-3


Introduction
Upon quantization, the axial chiral U(1) A -invariance is exonerated of representing a quantum symmetry in the theory of the strong interactions, namely Quantum Chromodynamics (QCD). The occurrence of an anomalous contribution in the QCD Lagrangian -linked to a nontrivial topological structure of the QCD vacuum -prevents the conservation of the axial current by dispossessing the theory from the Charge-Parity (CP)-symmetry. However -with astonishing experimental accuracy -QCD manifests itself as a CPpreserving framework and the anomalous "θ −term" seems to be fully irrelevant. This famous controversy -known as the strong CP-problem -is reconciled through the Peccei-Quinn mechanism [1], where a global U(1) PQ -symmetry is promoted to compensate the CP-violating term once it is broken spontaneously. While the described solution is perhaps the most appealing among other possible ideas, it did not come free of a new challenge for the contemporary physics. This is because the associated Nambu-Goldstone boson, i.e. the QCD axion, emerged as a plausible particle candidate [2,3]. Although soon after its introduction this hypothetical degree of freedom was ruled out [4,5], it certainly paved the way for more general class of axions embedded in the DFSZ and KSVZ axion models [6][7][8][9], this way opening a new front linked to their detection. Experimental endeavours toward this task are nowadays being carried out worldwide, based -predominantly -on the axion-diphoton coupling that is present in many effective Lagrangian densities such as axion-electrodynamics (AED) [10]; a framework whose phenomenology is subject to an intense scrutiny owing to its relevance, not just in particle physics, but also in research branches sharing its main features [11][12][13].
At lab scales, direct searches of Axion-Like Particles (ALPs) rely on plausible traces of axion-photon oscillations mediated by a constant magnetic field. Descriptions of the most popular detection methods can be found in Refs. [14][15][16][17][18][19]. The realization of this process could be relevant, indeed, in polarimetric arrangements [20][21][22][23][24] 1 and in Light-Shiningthrough-a-Wall setups [31][32][33][34][35][36][37][38]. Although the outcomes of these ongoing experiments have not verified the existence of the QCD axion yet, valuable bounds in its parameter space have been inferred instead, as well as in ALPs that are predicted in string theory [39][40][41][42] and in various Standard Model extensions which attempt to incorporate the dark matter of our Universe [43][44][45][46][47][48][49]. Constraints on axions -or, more general, ALPs -are also deduced from their potential astrocosmological consequences which are not reflected accordingly by the current observational data of our universe [50][51][52]. A basic assumption underlying this line of argument is that the interplay between the axionic degrees of freedom and the well established Standard Model sector -i.e., photons in first place -is extremely feeble. As a consequence, ALPs produced copiously in the core of stars via the Primakoff effect might escape from there almost freely, constituting a leak of energy that accelerates the cooling of the star and, thus, shortens its lifetime. Therefore, the number of red giants in the helium-burning phase in globular clusters should diminish considerably. That this fact does not take place -at least not significantly -constraints the axion-diphoton coupling g to lie below g 10 −10 GeV −1 for ALP masses m below the keV scale [53][54][55].
Precisely on the surface of stellar objects identified as neutron stars [56][57][58][59] and magnetars, magnetic fields as large as B ∼ O(10 13 -10 15 ) G are predicted to exist. As these strengths are bigger than the characteristic scale B 0 = 4.42 × 10 13 G of Quantum Electrodynamics (QED) [60,[62][63][64], such astrophysical scenarios are propitious for the realization of a variety of yet unobserved quantum processes, which are central in our current understanding on the nature and origin of the pulsar radiation. Notable among them are, the photon capture effect by the magnetic field owing to the resonant behavior [65] of the vacuum polarization in QED [66][67][68] 2 , the one-photon production of electron-positron pairs and the photon splitting effect [72][73][74]. Clearly, the strong-field environments provided by these compact objects can also be favorable for an ALPs phenomenology -as they are in QED -contrary to what is predicted relying on the weak coupling treatment. This occurs, because the aforementioned field strengths could compensate for the weakness of the coupling and significantly stimulate quantum vacuum fluctuations of axion-like fields. This paper is devoted to analyzing this possibility in the ALPs phenomenology. Theoretically, the problems linked to field strengths overpassing the corre-sponding characteristic scale of AED have not yet attracted adequate attention, save Refs. [75,76], where the virtualfermion-induced correction to the ALP mass in an external magnetic field of arbitrary strength was determined via the corresponding self-energy operator. This, however, involves a coupling constant which differs from the one mediating the interaction between an ALP and two photons.
In this paper we provide a further step toward the understanding of this subject by studying the characteristic magnetic field regime resulting from AED for which the polarization operator is responsible. We shall demonstrate that for sufficiently large magnetic fields, the axion-like fields lead to both stronger refraction and screening effects as compared with the predictions of QED [at least when we are kinematically far from the vacuum polarization resonance]. Hence, in the first part of the paper we reveal how virtual ALPs might be responsible for the capture of γ -quanta by the magnetic field lines in the magnetosphere of neutron stars and magnetars before the analogous effect starts acting in QED. We argue that such a phenomenon can prevent the photo-production of an electron-positron pair, and thus, the main channel responsible for the generation of this type of plasma. The suppresion of this production mechanism is incompatible with the existence of the plasma itself, which is understood nowadays as an essential part in the formation of the pulsar radiation. As a consequence of this inconsistency, new bounds on the parameter space of ALPs are inferred by taking as reference the magnetic field associated with two pulsars: SGR 1806-20 [57,58,[77][78][79] and RX J1856.5-3754 [80][81][82][83][84][85][86]. We shall mention at this point that considerable attention is currently paid to the possibility of probing ALPs via observational effects induced by them in neutron star magnetospheres [87][88][89][90].
The second part of our investigation is oriented to study the screening that virtual ALPs induce on the running QED coupling and so on the Coulomb potential of a static pointlike charge, when placed in a strong magnetic field. Contrary to studies developed within QED [91][92][93][94], our investigation reveals that along the field direction, and within the regime where the associated ALPs phenomenology could dominate over the corresponding QED effects, the screening due to ALPs is very weak, no matter how strong the field strength is (within the limits imposed on it by the unitarity). Hence, the modified potential we obtain follows the Coulomb law approximately. As a consequence, the existence of certain ALPs at extremely large magnetic fields might lead to a reestablishment of the formula for the unscreened groundstate energy of a nonrelativistic electron in the hydrogen atom [95]. The longstanding problem linked to this formula, i.e. the fact that the energy spectrum is unbounded from below when the magnetic field grows unlimitedly, was solved when the screening induced by the vacuum polarization of QED was incorporated and the virtuality of ALPs was ignored [91,92,[96][97][98][99]. We will see that in AED, by requir-ing the unitarity of the scattering matrix, a limitation arises for the largest treatable magnetic field, which may cause dramatic changes in judgements about the limiting ground-state energy. Also we show that, in a vicinity of the point-like charge perpendicular to the magnetic field, the Coulomb potential is of Yukawa type. This feature allows for certain energy regimes in which the electromagnetic interaction becomes very weak.
The investigation carried out here relies on a quantized approach of AED developed previously in Ref. [111]. In Sect. 2.1 the leading order contributions to the polarization tensor in the presence of an external electromagnetic field and a homogeneous axion background are established. Later on, in Sect. 2.2 the analisis is restricted to the case in which the latter is absent and the electromagnetic background is a constant magnetic field. Parallely, an approximationindependent eigenmode expansion of the polarization tensor is given in terms of its eigenvectors and eigenvalues, valid in a magnetic field both in QED and AED. We determine the polarization tensor with an accuracy of the second order in the axion-diphoton coupling with the magnetic field. Although the neutral nature of quantum vacuum fluctuations due to axion-like fields contrasts with those associated with the charged Dirac ones, the outcomes resulting from the axion-photon theory pretty much resemble those obtained in QED. This similarity is stressed in Sect. 2.2, where various asymptotes of the eigenvalues of the polarization tensor are established in parallel. Also in that section the plausible regions of parameters of undiscarded ALPs, where the contribution of the latter could dominate over the QED effects, are established.
In Sect. 3.1 we define the mass shell by solving the dispersion equation for a photon of the extraordinary polarization mode interacting with an ALP in the presence of a magnetic field to find two branches of the spectrum: the massless one and the massive one for the axion-like particle. Parallely, we find that the flattening of the massless branch near small values of the transverse momentum results in vanishing of the group velocity component across the magnetic field. This phenomenon forms a basis for the photon capture considered in the subsequent sections 3.2 and 3.3. In the former we recall the effect of γ -quantum capture in QED and argue that in AED the same effect might destroy the accepted mechanism of electron-positron plasma production in pulsars and hence of their radio-emission, if it was not eliminated by the already existing restrictions on the ALP mass and coupling constant. In Sect. 3.3 we deal with the capture of soft photons that occurs in AED. We consider the refraction of radiation, when it leaves the region of strong field, and assure ourselves that the heat radiation from the pulsar surface belonging to the extraordinary mode is canalized parallel to the magnetic field, and it may leave it, whereas the radiation polarized in the orthogonal direction is not canalized and is subject to the angular distribution determined by the refraction law. We believe that the effect that the radiation of a certain polarization condenses close to the magnetic field direction, while the orthogonally polarized radiation is distributed over the angle, might become a valuable criterion for probing ALPs among the plausible signatures induced in the thermal spectra of neutron stars [86,100].
Finally, in Sect. 4.1 we determine the photon propagator in order to use it subsequently in Sect. 4.2 in the establishment of the running QED coupling. We show that, in the regime where AED dominates over QED, the screening caused by the magnetized vacuum can screen a point-like charge to zero when the magnetic field grows unlimitely at a fixed momentum. Properties of the electromagnetic interaction in this axion-dominating scenario are further inferred via the axion-modified Coulomb potential. Its anisotropization, short-ranging and dimensional reduction are revealed manifestly through numerical evaluations. We show, in particular, that for field strengths larger than the characteristic scale associated with this effective theory, the modified potential is short-ranged in the direction perpendicular to the magnetic field while it follows approximately the behavior of Coulomb's law along B B B. Section 4.3 is devoted to support analytically these results. There, asymptotic expressions for the modified Coulomb potential in the limits of strong magnetic fields are determined. In addition, the imprints that virtual ALPs leave within the ground state energy of a nonrelativistic electron in a hydrogen atom are evaluated. Our conclusions are exposed finally in Sect. 5, whereas some important details of our calculations are presented in Appendices A and B.

General considerations in a constant magnetic field and a homogeneous axion background
To leading order in the axion-diphoton coupling g, the Lagrangian density describing AED combines the standard Maxwell Lagrangian with the free Lagrangian density of the pseudoscalar field φ(x) and with the scalar interaction Lagrangian density involving the dual of the electromagnetic field tensorf μν = 1 2 μναβ f αβ where f μν = ∂ μ a ν − ∂ ν a μ and 0123 = 1. Explicitly, Hereafter, we use a metric with signature diag(g μν ) = (1, −1, −1, −1), and Heaviside−Lorentz units with the speed of light and the Planck constant set to unity c =h = 1.
To abbreviate, the following notations have also been used: In contrast to the latter, the former is mediated by quantum fluctuations of pseudoscalar fields. Here, the horizontal dashed line represents the ALP propagator Δ F (x,x), whereas the horizontal wavy lines must be understood as amputated photon legs. The lowest diagram is the contribution resulting from the original axion-diphoton vertex Here m refers to the ALP mass, whereas g denotes its coupling strength to two photons. As the interacting term has a mass dimension −5, this framework belongs to the class of perturbatively nonrenormalizable theories. In this connection, L can be regarded as a Wilsonian effective Lagrangian parametrizing the leading order contribution of its ultraviolet completion, i.e. a renormalizable theory linked to physics beyond the Standard Model which might contain new heavy particles at the energy scale Λ UV ∼ 1 g . If this hypothetical theory satisfies fundamental principles such as Lorentz and gauge invariance, and preserves -in addition -the unitarity, quantum processes due to fluctuations of modes of the fields in L, taking place at energies substantially below Λ UV , should also be described in a unitary way.
In the following we evaluate consequences of the vacuum polarization tensor in AED under the influence of a con- (1) by f μν +F μν , and due to this operation, the Lagrangian density L acquires an additional interaction L F = − 1 2 gφF f , which is responsible for the appearance of the upper Feynman diagram shown in Fig. 1. Since cold relic axions -resulting from vacuum misalignment in the early universe [45][46][47]55] -are a viable candidate for dark matter, quantum fluctuations of the axion-like field φ(x) might occur not on an empty vacuum but rather over a background field which is supposed to fit the problem of the missing mass of the universe. This coherent classical field is identified as a solution of the equation of motion in the presence of the magnetic field. Ignoring the effect due to the local virial dark matter velocity v ∼ 10 −3 , it can be approximated by a homogeneous sinusoidally time-varying field: 3 with DM ≈ 0.3 GeV/cm 3 referring to the local dark matter density [102,103], Ω ≡ 1 + g 2 B 2 /m 2 and τ ≡ mt. However, the existence of Eq. (2) comes accompanied with a time-dependent electric field parallel to B B B As a consequence, the propagation of photons actually takes place on an electromagnetic field background whose invariant [104]. Clearly, the inclusion of Φ(t) in L is carried out via the replacement φ → φ + Φ(t). This operation leads to the emergence of an interaction vertex which is exhibited in the lower panel of Fig. 1. At this point we find convenient to remark that the inclusion of the axion background is only meaningful if a certain ALP is identified as the main constituent of the dark matter. ALPs which are not related to the latter scenario are by themselves well motivated particle candidates and, indeed, most of the so-called direct searches are oriented to look for ALPs belonging to this last group. For details we refer the reader to Refs. [14][15][16][17][18].
As the graphs mentioned so far drive the dispersive and electrostatic effects in AED, it turns out to be convenient to have a theoretical formulation in which its analytical expression appears explicitly. To this end, we rewrite the effective action of the theory S = d 4 xL in the following form: Notice that the last term in the second line of this expression is responsible for the vertex involved in the one-loop diagram [lowest graph depicted in Fig. 1]. The calculation of this Feynman diagram was developed in Ref. [111], and unless necessary, its contribution will be ignored from now on. To bring the theory into a tractable form, the pseudoscalar ALP field φ(x) is integrated out from the generating functional of the photon Green's functions where j ≡ j μ (x) is the external source and the term proportional to ζ −1 guarantees a covariant quantization. The effective action that results from this functional integration is Here p 2 −m 2 +i0 stands for the ALP propagator with the i0−prescription allowing to elude the pole at p 2 = m 2 . The source term contained within Γ [a], i.e. J μ (t) = gη ν (t)F νμ , depends upon the strong magnetic field and the homogeneous time-dependent axion dark matter background The time dependence of this four-pseudovector field leads to a violation of the Lotentz-invariance, in addition to the one induced by the background magnetic field. Besides, the bilinear term in the expression above defines the inverse photon Green function, where = ∂ 2 /∂t 2 − ∇ 2 is the d'Alembert operator and the leading-order term of the polarization tensor Π μν (x,x) combines the two contributions represented in Fig. 1: While the first line in the expression above is a symmetric nonlocal tensor depending quadratically on the external electromagnetic field, the term contained in the second line is antisymmetric translation-non-invariant and local, depending linearly on η μ (t). The presence of this contribution within the polarization tensor resembles the Chern-Simons term found in the context of the standard model plasma [105], with the particularity that η α plays the role of the axial vector responsible for the generation of long-range magnetic fields in the ground state. Its occurrence induces properties on the quantum vacuum akin to a medium in which -in addition to the breaking of the Lorentz and time-translational symmetries -the parity invariance is no longer preserved. It is worth pointing out that the propagation of electromagnetic waves in parity violating scenarios has been studied previously [see for instance Refs. [106][107][108][109]]. In contrast to Eq. (8), the strongest nonlocal interacting term -second line in Eq. (6) -depends linearly on the external background and quadratically on small-amplitude electromagnetic waves. This feature guarantees that, for very strong magnetic fields, a controlled perturbative expansion in the coupling g 2 can be carried out. Consequently, corrections to Π μν (x,x) will always remain smaller as compared to the expression above, independently of how strong the external field is [for details see Appendix A]. Hence, Eq. (8) should be understood as the starting point of further considerations. Observe that the treatment followed so far provides the following equations of motion for the small-amplitude electromagnetic wave [∂a = 0]: We anticipate that this dispersion equation is essential to establish the capture effect in AED (see Sect. 3).
At this point it is worth mentioning that when both η α (t) and E (t) vary slowly, the problem under consideration absorbs some features linked to QED embedded into a combination of constant magnetic B and electric E fields studied earlier in the works published under the participation of one of us [see [113][114][115]]. However, in order to facilitate the understanding the effect of photon capture, we shall refer hereafter to the case in which ALPs and dark matter are unrelated quantities. An analysis including both the effect of the cold axion background and the external magnetic field is beyond the scope of this manuscript. Nevertheless, the reader interested in this subject may find attractive the predictions reported in Ref. [110], where -contrary to our case -the magnetic field has been taken B m/g.

Polarization tensor, eigenvalues and dominance over QED counterpart
When the external field is a constant magnetic-like background characterized by A μ (x) = − 1 2 F μ ν x ν and the invari- for ALPs which do not represent the main constituent of dark matter, the Fourier transform of the polarization tensor can be written in the following form iΠ μν (q 1 , q 2 ) = (2π) 4 δ 4 (q 1 − q 2 )P μν (q). Because of the Dirac delta δ 4 (q 1 −q 2 ), the tensorial portion P μν (q) does not depend on which choice q = q 1 or q = q 2 is taken. Besides, it turns out to be symmetric P μν (q) = P νμ (q), transverse q μ P μν (q) = 0 and can be diagonalized in a similar way as it occurs with the vacuum polarization tensor established for QED with magnetic field [114][115][116] The scalar eigenvalues j in the approximation considered are where the form factor linked to the one-loop correction satisfies π(q 2 ) 1 provided q 2 is not exponentially large. where π(q 2 ) ∼ g 2 m 2 1 is the form factor linked to the one-loop correction [lowest graph in Fig. 1]. As we are interested in the strong field regime B m/g, the contribution due to the latter can be safely ignored. Note that the external field and hence the anisotropy is introduced through the interaction with an ALP. In the special Lorentz frame, where the external field is purely magnetic and directed along the axis denoted as , the Lorentz scalar qF 2 q/(2F) becomes q 2 0 −q 2 where the photon energy is q 0 and the momentum component along the magnetic field is q q q . There is also an extra Lorentz scalar qF 2 q/(2F) = q 2 ⊥ , with q q q ⊥ denoting the projection of the momentum q q q onto the plane perpendicular to B B B, so that 2Fq 2 = qF 2 q − qF 2 q. At this point it is worth remarking that the set of eigenvectors μ j with j = 1, 2, 3 is built from first principles, independently of any approximation used in the calculations [114,115]. Explicitly, The eigenvectors contained in this decomposition are mutually orthogonal jμ μ = −δ j 2 and transverse q μ μ j = 0 with j = 1, 2, 3. Besides, they satisfy the completeness rela- j . When AED is embedded in the QED action, the eigenvalues of the polarization tensor [see Eq. (11)] add to the corresponding eigenvalues of the QED polarization tensor. In this subsection we exploit this property to establish regions of ALPs parameters for which a dominance over the QED effects could take place. This seems to be particularly accessible because the nontrivial eigenvalue of the vacuum polarization tensor, i.e. 2 [see Eq. (11)] exhibits a quadratic growth in the external field while the corresponding eigenvalue in QED grows linearly [118][119][120][121] Here α = 1/137 denotes the fine structure constant, b = B/B 0 and B 0 = m 2 e /e ≈ 4.42 × 10 13 G the characteristic QED scale. In this context, m e and e refer to the electron mass and the absolute value of its charge. It is worth mentioning that the application of these asymptotes requires q 2 ⊥ m 2 e b. In order to establish undiscarded ALPs parameters which might lead to strong refraction as well as to a more pronounced screening property than the one predicted from QED, we first consider the situation in which q 2 0 − q 2 min{m 2 , m 2 e } and q 2 ⊥ min{m 2 , m 2 e b}. In this case the nontrivial eigenvalue of the polarization tensor [see Eq. (11)] is real: and a direct comparison with the infrared behavior of QED 2 [first line in Eq. (13)] can be established. Here, we have introduced the parameter b = B/B c with B c = m/g denoting the characteristic field scale associated with AED. We will see very shortly that this scale plays in AED the role that B 0 does in QED, in the sense that for magnetic field strengths larger than it [B B c ], the effects induced by the vacuum polarization mediated by ALPs are expected to become much more pronounced than in the contrary situation [B B c ]. Now, as the asymptotic expression given in Eq. (14) has the same momentum dependence as its QED counterpart [first line in Eq. (13)], one can establish a region of dominance of AED over QED directly on the parameter space of ALPs. Indeed, by demanding that Eq. (14) exceeds the corresponding QED expression we find b 2 > αb/(3π). This condition leads to the relation To evaluate the consequences of this inequality we shall refer to various existing estimates of the ALP parameters g and m.
Let us star with considering the best laboratory limit obtained by the OSCAR collaboration: g < 4 × 10 −8 GeV −1 for m 100 µeV [38]. For the boundary values in this range, the characteristic magnetic field B c ≈ 1.2 × 10 14 G is a little above the value of characteristic magnetic field in QED, B 0 = 4.4 × 10 13 G. The latter is comparable with the magnetic fields that might exist in the magnetosphere of neutron stars and magnetars. For a magnetic field exceeding this characteristic scale by two orders of magnitude -which could be realized at the surface of cosmological gamma-ray bursters if they are rotation-powered neutron stars [126][127][128] -we obtain that b 2 ∼ 10 4 is four orders of magnitude larger than αb/(3π) ∼ 1. Under such a circumstance the vacuum polarization linked to virtual ALPs would prevail over the one associated with virtual electron-positron pairs. We remark that much more stringent bounds -resulting from the nonobservation of astro-cosmological consequence linked to ALPs -are available [see Fig. 3]. For instance, the modified solar sound-speed profile can be monitored by helioseismology, providing a conservative upper limit: g < 5 × 10 −10 GeV −1 for m ∼ (10 −5 -10 −3 ) eV [55]. With these estimates the characteristic value B c ∼ (10 15 -10 17 ) G remains within the limits admitted in the core of neutron stars 10 18 G, 5 and the dominance of AED over QED could also take place. To gain a more extended idea, we have inserted in Fig 11)] is also real: and can be compared with the corresponding QED expression given in the second line in Eq. (13). Hence, the vacuum polarization of AED would dominate over the QED contribution if mb (2αb/π ) 1 /2 m e corresponding to g m e [2α/(B B 0 )] 1 /2 . Observe that, if the largest magnetic field admitted in neutron stars and magnetars is 10 18 G, the previous condition translates into g 0.47 GeV −1 . We find opportune to emphasize that this value is out of the region exhibited in Fig. 3, and so discarded for any ALP mass m embedded within the interval (10 −8 -10 9 ) eV. Therefore, in the energy regime analyzed above no dominance of AED over QED occurs for B 10 18 G. However, for a mag- 5 Studies relying on the Virial theorem corroborate this upper limit. For details, we refer the reader to Ref. [112].
netic field B ∼ 10 22 G the dominance could take place if g 10 −3 GeV −1 . For ALPs masses m < 10 MeV, values of g overpassing the previous limitation have been ruled out [see Fig. 3]. A small window 10 −2 GeV < m < 1 GeV with g < 10 −2 GeV −1 remains for which B c ∼ (10 19 -10 21 ) G. While fields of this order of magnitude are not within the reach of neutron stars and magnetars, it is likely that they have existed right after the electroweak phase transition [122,123].
We want to conclude this section indicating that the regions of dominance found here have to be considered as exploratory. This is because, in an axion-dominating regime the fine-structure constant used so far could undergo a substantial screening. Indeed, a more elaborate procedure will require to replace α by the value of the running QED coupling α scr (q q q) which depends upon the magnetic field strength and the typical momentum transfer of the physical scenario. This problem is further investigated in Sect. 4.2.

Dispersion laws and group velocities
Although resulting from different nature, the singularity exhibited in the second eigenvalue at q 2 = m 2 [see Eq. (13)] somewhat resembles the known situation in QED, where there appears an inverse square root singularity in QED 2 at the border to the continuum of free electron and positron states produced by a single photon: the cyclotronic resonance of the vacuum [65,116], or the pole where the electron and positron are bound to form a positronium atom [69][70][71]. The genuine dispersion relations associated with the mode-2 are neither q 0 = |q q q| nor q 0 = q q q 2 + m 2 but rather Observe that the dispersion law labeled with a negative sign describes a massless particle as q 0− vanishes when both components q ⊥ and q vanish. Conversely, the dispersion relation q 0+ corresponds to a massive branch because lim q ,⊥ →0 q 2 0+ = m 2 * with m * = m √ 1 + b 2 the effective mass, which reduces to the ALP mass m only when the external field is zero. Therefore, a small-amplitude electromagnetic wave characterized by the second propagation mode is actually a state in which massive and massless particles coexist. Both branches of the dispersion law [see Eq. (17)] are depicted in Fig. 2 for various values of the b parameter. As b grows, the massless dispersion curves tend to stick to the horizontal axis. This behavior can be understood from Eq. (17) when the case q ⊥ m * is considered. To obtain the corresponding asymptotic formula one has to take into account that, for any value of b, the factor 4b 2 /(1 + b 2 ) 4. Hence, if q ⊥ m * /2, the second term within the squared root is very small. Consequently, This result describes the trend that q 0− exhibits in Fig. 2 when the magnetic field grows gradually and applies whenever m * q ⊥ . Hence, the pattern in this figure obtained is to be interpreted in such a way that an ALP and mode-2 photon interact in the magnetic field to produce the spectrum of two disconnected branches: one massless and one massive. The massless branch may be referred to as belonging to a photon modified by its interaction with an ALP, while the massive branch corresponds to an axion-like particle modified by its interaction with the mode-2 photon. No mixed state between the axion-like branch and the photon is formed. Here we encounter a vast difference with QED in this type of background, where the free photon dispersion curve quasi-intersects with the dispersion curve of an electron-positron pair -mutually free or bound into a positronium atom -before the interaction between the photon and the pair is taken into account. This quasi-intersection results in forming the photon-pair mixed state, the polariton, after this interaction comes into play [66][67][68][69][70][71]. In such a context, photon-positronium oscillations occur [130,131].
When calculating the components of the group velocity v ⊥, = ∂q 0− /∂q ⊥, linked to Eq. (18), we find v ⊥ ≈ As a consequence, the angle θ between the direction of propagation of the electromagnetic energy and the external mag- does not coincide with the one [tan ϑ = q ⊥ /q ] between the photon momentum q q q and B B B, [tan θ < tan ϑ]. Furthermore, we observe that it tends to vanish when q ⊥ m * ≈ mb the faster, the stronger the field is, since one has [see Eq. (19) This implies that the group velocity of massless excitations tends to be parallel to the magnetic field for hard, as well as for soft photons with q ⊥ m * . The described situation resembles the main feature on which the capture effect of photons with energy much lower than the first pair creation threshold relies in QED [68]. Following our discussion below Eq. (14), we expect that this phenomenon could also be realized via a certain class of virtual ALPs rather than by virtual electron-positron pairs, provided b 1. We will see very shortly that, in addition to γ −quanta [66,67], and as a direct consequence of Eq. (19), X-rays, optical light [68], infrared radiation and electromagnetic microwaves could undergo the capture due to the magnetized vacuum of AED near the surface of neutron stars.
Let us investigate the refractive properties linked to the massive branch [lower sign in Eq. (17)]. When b 1 and q ⊥ m * ≈ mb the associated dispersion relation approaches q 2 0+ ≈ 2q 2 ⊥ + q 2 + m 2 * . In this case the com- As a result, the cosine of the angle η between the group velocity vector v v v and the momentum is In terms of the angle ϑ between the photon momentum and the magnetic field [see below Eq. (19)], this expression reads cos η ≈ 1 + 2 tan 2 ϑ √ 1 + 4 tan 2 ϑ cos ϑ.
Within the range 0 < ϑ < π/2 the angle η remains very small and achieves its maximum value of η ≈ arccos[0.943] = 0.339 rad at ϑ ≈ 36 • . This is the deflection that the direction of propagation of the wave undergoes from its momentum vector due to the anisotropy of the medium.
3.2 Implications of plausible captures of γ γ γ quanta mediated by virtual axion-like particles Let us now consider the photon capture effect in different terms. In understanding that the group velocity specifies the direction of propagation of the wave envelope, the suppression of its perpendicular component implies the effect of capture of the mode-2 photon by the strong magnetic field, as described previously in QED. We will see, however, that this phenomenon has quite different features in the present context. Observe that Eqs. (19) and (22) apply for q ⊥ m * with b 1. Beyond this limit, i.e. for q ⊥ m * , the components of the group velocity approach and the angle θ = arctan(v ⊥ /v ) between the direction of propagation of the electromagnetic energy and B B B does not deviate substantially from the one formed by the momentum of the small-amplitude wave and the external field [θ ≈ ϑ = arctan(q ⊥ /q )]. Suppose that we have a magnetic field with its lines of force curved as it is adopted in models describing the magnetosphere of pulsars and neutron stars. Whenever the curvature radius of the field is much larger than the photon wavelength and the Compton wavelength of ALPs, the analysis developed so far will also apply. When propagating in such a field a photon changes its momentum components q ⊥ and q , while its energy q 0 is kept constant, because the external field does not depend on time. Thereby the photon evolves along its dispersion curve. Let primarily the so-called curvature γ -quantum with q 0 > 2m e be injected in the direction parallel to the magnetic field, as it is the case within the widely accepted scenario of the events that take place in the pulsar magnetosphere. The curvature gamma-quantum is emitted by an electron propagating along a curved magnetic line of force. As the photon propagates inside the dipoleshaped magnetic field it gains a crucial pitch angle with the direction of this field, kinematically sufficient for the photon to create an electron-positron pair. This takes place when q 2 0 − q 2 = 4m 2 e . A "classical" on-shell photon is characterized by the vacuum dispersion law q 2 0 − q 2 = q 2 ⊥ [dashed line in Fig. 2] unless the magnetic field exceeds too much the characteristic scale B = 0.1 m 2 e /e ≈ 4.42 × 10 12 G, in which case the resonant phenomena [66] must be taken into account. 6 Therefore, before reaching the necessary pitch angle -which happens for mode-2 photons at q ⊥ = 2m ethe photon propagates straightforwardly in the coordinate space and then it produces the pair. This is the mechanism of the formation of the electron-positron plasma by mode-2 photons responsible for the further directed electromagnetic emission of the pulsar. Now, in accordance with Eq. (21) the admission of the axion-photon interaction makes the curvature photon of mode-2 be captured by the magnetic line of force practically immediately after the photon is emitted, provided q ⊥ m * . This condition might lead to the capture of massless mode-2 photons with energy larger than the characteristic energy scale associated with the electron mass, i.e. for q 0 > m e . However, if this was the case, a single sufficiently hardbelonging to the γ -range -photon in the magnetic field has no chance to produce a pair for long, such that the established mechanism of the pulsar radiation would be destroyed, provided B c B 0 . 7 Indeed, following our preceding discussion, in the magnetic field the production of an electronpositron pair by a γ -photon is possible if q 2 0− − q 2 4m 2 e . 8 Substituting q 2 0− − q 2 = 4m 2 e into Eq. (17) we obtain the equation for the threshold value of q 2 ⊥ that borders this process. Its solution is ⎠ . 6 When taking into account the resonant phenomena owing to the singularity of the polarization tensor in momentum space in the point, where the photon creates a free or bound pair (the positronium atom), Footnote 6 continued it was shown that the capture effect is present in QED. Its role in QED was reduced to a certain increase of the interval where he photon freely propagates before it turns into the pair. We shall not take this correction into account in the consideration that follows. 7 As for mode-3 photons being not affected by the presence of ALPs these could contribute to the plasma creation only if their energy is larger than m e 1 + (1 + b) 1/2 > 2m e . For sufficiently large magnetic field they cannot save the necessary plasma creation mechanism resting on the less energetic mode-2 photons, once the latter ones are excluded from the game by the presence of ALPs, as described above. 8 Kinematically, the process of pair creation by a "photon" belonging to the massive branch is allowed provided q 2 0 − q 2 4m 2 e . However, the curvature γ -quantum, i.e. the one emitted in the pulsar magnetosphere by an electron moving along a curved magnetic line of force, is presumably the standard zero-mass photon. The option of creating excitations belonging to the masive branch is not seen. Thereby, the plasma creation by the massive-branch "photon" is left beyond discussion.
The accepted mechanism of plasma generation in neutron stars is destroyed if it may take place at the values q 2 ⊥ thr exceeding substantially the threshold value of q 2 ⊥ thr = 4m 2 e characteristic of QED in a magnetic field. This occurs if the three conditions 4bm e m, mb 2m e and b 1 (27) are met simultaneously, along with B c B 0 . In such a case the photon is captured by the magnetic line of force, and it cannot create a pair before it has propagated far enough along the line to leave the region of the strong dipole field of the pulsar where, however, the one-photon pair creation is strongly suppressed. Clearly, the occurrence of this phenomenon would suppress the main channel responsible for supplying an electron-positron plasma in the neutron star magnetospheres. This suppression would forbid the existence of the plasma itself, and therefore, the so far accepted source of pulsar radiation. Other recognized radiation emission mechanisms such as the splitting of photons turn out to be suppressed as well because -in the regime under consideration -QED undergoes a sort of quantum triviality [read discussion below Eq. (38)]. The incompatibility between both models demands, either to look for alternative mechanisms that would explain the pulsar emissions in the presence of virtual ALPs, or to exclude the (g, m) region dictated by Eq. (27) provided the field strength is known from observations. While the determination of the magnetic fields linked to neutron stars are made indirectly and turn out to be modeldependent, it is natural to investigate the upper bounds that result from the latter possibility.
Let us consider the scenario provided by RX J1856.5-3754 -a representative pulsar belonging to the seven nearby X-ray thermal isolated neutron stars known as "Magnificent Seven"-to which considerable observational resources have been devoted [80][81][82][83][84][85][86]. In the surface of RX J1856.5-3754 the magnetic field is B ≈ 1.5 × 10 13 G [83]. The area excluded by the existence of an electron-positron plasma in the surface of this neutron star is colored in light-blue [see Fig. 3]. Our result is expected to be reliable whenever the ALP mass m exceeds substantially the energy scale linked to the curvature radius ∼ R −1 curv , a restriction introduced by the adiabatic approximation. Hence, if we identify the curvature radius R curv with the star radius R ∼ 12 km, the previous condition will allow us to probe ALPs with masses m 10 −11 eV. More stringent constraints can be established from neutron stars holding a more intense magnetic field. At the time of writing, the strongest magnetic field belong to the soft γ −repeater SGR 1806-20 [57,58]. Although a field as strong as B ≈ 7×10 15 G has been inferred for this magnetar [77], in our numerical assessment we used a more conservative value B ≈ 2.4×10 15 G [78,79]. It is worth mentioning that, for this field strength, the effect of the QED vacuum polarization [see  [133,134]. Colored in green and light red appear the regions ruled out by particle decay experiments. Furthermore, the portion discarded by investigating the energy loss in horizontal branch (HB) stars are shown in dark red, whereas the excluded area resulting from the solar monitoring of a plausible ALP flux (CAST+SUMICO) has been added in olive [124,125]). In light brown the portion discarded by measuring the duration of the neutrino signal of the supernova SN1987A is depicted. The light blue area results from cosmological investigations Eq. (13)] still can be treated as a perturbation. The respective region ruled out by the noncapture of a γ −quanta mediated by ALPs, prior to the creation of an electron-positron pair is exhibited in Fig. 3 [light red colored]. Noteworthy, the area discarded in this way turns out to be much wider than the one linked to RX J1856.5-3754, and improves existing bounds in the region m ∼ 10 −10 -10 −5 eV.

Magnetically-induced canalization of linearly polarized radiation as a signature of axion-like particles
The direct thermal radiation emitted by the hot surface of neutron stars is directly observable in the X-ray range [see, e.g., Ref. [132]]. In this subsection we inspect the consequence associated with the capture of their mode-2 part by the strong magnetic field due to the presence of ALPs, as pointed out below Eq. (21). In this connection, we will start demanding that q 2 < m 2 e and q 2 ⊥ m 2 * < m 2 e , so that the photon energy q 0 m e . Observe that the previous conditions guarantee the applicability of Eq. (18), and therefore, the capture of heat radiation due to a plausible dominance of AED over QED.
Let us consider the refraction law of electromagnetic waves -with energy belonging to the thermal range -when they leave the region of strong magnetic field and enter the space free of it. This consideration will be independent of what mechanisms are responsible for formation of the vacuum polarization and dispersion laws, in other words of whether there are ALPs or we just face the usual QED with magnetic field. Suppose that an electromagnetic wave is created at the surface of a neutron star inside a strong magnetic field, which gradually disappears. Conventionally, the field B B B changes only along its own direction, but it is transversely homogeneous. We assume that both the photon wave length λ and the Compton wave length of the ALP are much smaller than the length scale over which the magnetic field varies substantially. Under such a restriction, the adiabatic approximation is valid, wherein the dielectric properties of the medium may be considered as depending on the coordinate. Then the photon energy q 0 and the photon momentum component q ⊥ across the magnetic field do not change, when the wave leaves the region occupied by the magnetic field, due to the corresponding space-time translational invariance. On the contrary, the momentum component along the magnetic field q does change, because the magnetic field varies in this direction. We shall supply this component with the subscript "∞", when it relates to the region in which the magnetic field is not present.
There are two different dispersion laws: The last relation is a consequence of relativistic covariance, the function f (q 2 ⊥ ) is model-dependent and it is different for each propagating mode. According to Eq. (18) one has provided the conditions q ⊥ m * and b 1 are fulfilled. 9 The refraction index at infinity is n ∞ = |q q q ∞ |/q 0 = 1, whereas at the magnetic cap of the star n = |q q q| The adiabatic approximation allows us to divide the region between the star surface and the observatory into a set of equally thin portions with length so that -within the boundaries of each of them -the corresponding refraction index can be assumed constant [see Fig. 4], and the Snell's law can be applied to each interface: n sin ϑ = n 1 sin ϑ 1 = n 2 sin ϑ 2 = . . .
where ϑ ∞ (ϑ) is the angle between the outgoing [incoming] momentum q q q ∞ (q q q) and the magnetic field direction. With the use of the relation sin ϑ ∞ = q ⊥ /|q q q ∞ | = q ⊥ /q 0 and the substitution of Eq. (30) into the expression above we find 9 Note that in mode-3 the function f (q 2 ⊥ ) is proportional to q 2 ⊥ at small q ⊥ . Therefore, this region of momentum is an arena for utmost manifestation of birefringence. This can also be written as provided that q 0 sin ϑ ∞ m * and b 1. The critical angle ϑ cr for total internal reflection is achieved when ϑ ∞ = π/2. It depends on the energy: Only approximately one half of the isotropic long-wave radiation has a chance to leave the magnetic field. The capture effect further modifies the situation. We have found [see Eq. (19)] that on the massless photon branch, for q 2 ⊥ as small as q 2 ⊥ = q 2 0 sin 2 ϑ ∞ m 2 * all photons are canalized towards the direction of the magnetic field. The fulfillment of this inequality is guaranteed for all incident angles provided that the photon energy q 0 m * < m e . As the previous condition is stronger than the one required at the beginning of this subsection, one is forced to limit q 2 m 2 * < m 2 e . On the other hand, the dispersion curve passes through the point ϑ ∞ = ϑ = 0, which means that the wave parallel to the magnetic field is not deflected, when it goes out of it. Therefore the whole mode-2 part of the radiation would gather in a beam going out of the magnetic field parallel to it. Another possibility is that q 0 m * . Under this condition the canalization of the second propagation mode takes place with the particularization that the refractive effect becomes very weak: 0 < sin ϑ ∞ m * /q 0− 1. In this case, the propagation of the gathered radiation would not deviate considerably from the magnetic field direction and the condition of dominance [see Eq. (18)] still holds.
As the mode-2, in the ALP-dominating regime, is canalized parallelly to the background magnetic field, most of its energy is concentrated at momenta |q ⊥ | m * . Following the analysis developed in Appendix B [see Eq. (B.3)], mode-2 is a transverse wave with its electric field lying in the plane spanned by the vectors B B B and q q q. In accord with the analysis above [see Eq. (33)], it leaves the background field with the refraction angle close to the incident angle, when both are small. On the other hand, the mode-3 component of the heat radiation, whose electric field is normal to the plane formed by the background magnetic field and the direction of propagation, is not canalized. (Recall, that it is not influenced by the presence of ALPs and remains as calculated via QED [see Ref. [68]]). Its original angular distribution is only modulated by the refraction law Eq. (33), and it is not concentrated along any direction. Therefore, if ALPs exist and the vacuum polarization induced by them dominates over the one of QED, the heat photons with energies q 0 m * < m e polarized in the plane to which the vectors B B B and q q q belong, peak along the direction of the magnetic field against the smooth background of photons polarized orthogonal to this plane. Hence, if the direction of the observed photon flux deviates from being parallel to the magnetic field we would face the smooth background of photons polarized orthogonal to this plane-that are not subject to the capture. The registered radiation would be scant of photons polarized in the said plane.
The peaking (or not peaking) of photons in the parallel direction constitutes a suitable way to make sure that ALPs exist (or do not). In the latter case, constraints on the parameter space of ALPs could also be inferred directly from the observed photon flux. Clearly, this procedure would require to derive an expression for the flux of canalized photons to be detected by including both the capture and the refraction effect. While the establishment of such a formula is beyond the scope of this investigation, we plan to treat this problem in a forthcoming publication.
It is also worth mentioning at this point that there exists already significant progress in high-purity polarimetric techniques for X-ray probes [135,136] which are expected to be exploited in the envisaged experiment at HIBEF collaboration [137,138]. These techniques are planned to be used in polarimeters designed to investigate the soft X-ray emission in the strong magnetic fields of neutron stars and white dwarfs [139][140][141][142]. 10 Eventually, these polarimetric measurements could allow us as well to probe the dominance of AED over QED, and provide new constraints on the ALP parameter space [143]. Indeed, if the proposed effect could be extinguished from the future observations one will be able to ban the values of the ALP coupling g > q 0 /B, where the pulsar value is taken for the magnetic field B and the x-ray value for the photon energy q 0 .

The modified Coulomb potential in axion-electrodynamics
In this section we study how the screening caused by a plausible dominance of AED over QED modifies the running QED coupling and so the Coulomb potential.

Photon propagator
Aside from influencing the propagation of photons, the polarization tensor [see Eq. (11)] also induces modifications of Coulomb's law. This distortion is determined here through the temporal component of the electromagnetic fourpotential Here this is accomplished by using the identity is given in Eq. (8) which corresponds to summing up the geometric series of one-photon-reducible diagrams. However, the close analogy between our problem and the one previously analyzed in QED allows us to apply the model-and approximation-independent expression for D μν (x,x) found in Ref. [91] directly. This general formula links the previous object with the corresponding eigenvalues and eigenvectors of the polarization tensor in QED. It only remains to carry out an appropriate replacement of these eigenvalues by those resulting from the photon-ALP oscillations. Hence, up to an inessential longitudinal contribution, the photon propagator reads where the shorthand notationd 4 q ≡ d 4 q/(2π) 4 has been used. We insert this formula into Eq. (35) bearing in mind the details explained below. As a consequence, the axionmodified Coulomb potential is given by When establishing this expression we used the fact that 0 1 = 0 3 = 0, i.e. that the first and third eigenmodes do not contribute to the electrostatic interaction. This feature could be anticipated because, in the stationary limit q 0 = 0, "virtual" photons associated with μ 1,3 carry magnetic fields only, while off-shell photons of mode-2 are purely electric.

Running QED coupling and Coulomb potential
Let us start this subsection by analyzing general features linked to the running QED coupling α scr (q q q) when the dominance of AED over QED polarization tensor takes place. To this end, we shall refer ourselves to the electrostatic energy U(x x x) = −ea 0 (x x x) between two electrons. The insertion of Eq. (37) into this formula with q → e defines the running QED coupling. Indeed, by expressing U( with α = e 2 /(4π) defined when the magnetic field background vanishes, i.e. at b → 0 and |q q q| → 0 [60,61].
Notably, if the dominance of AED over QED takes place, α scr (q q q) is free of the characteristic Landau pole arising in QED that rises the coupling constant to an infinite value for large enough momentum. 11 Furthermore, its dependence on the magnetic field leads to an anisotropic behavior along and transverse to B B B. This dependence decreases the value of the running QED coupling as the magnetic field B grows while keeping the momentum fixed. Indeed, for q 2 ⊥ q 2 m 2 , one finds that α scr (q q q) ∼ αb −2 1 tends to be screened to zero. This phenomenon makes QED an almost "trivial" theory in the aforementioned regime, meaning that the corresponding matter sector and the small-amplitude electromagnetic waves turn out to be very weakly interacting between each other.
The described situation becomes particularly dramatic if the axion mass turns out to be larger than the first pair creation threshold, i.e. m > 2m e because -in line with the discussion of Sect. 3.2 -this would imply a complete suppression of radiation emission mechanisms such as the recombination of pairs and the photon splitting effect as well as the main production channels in highly-magnetized pulsars. This quantum triviality [146][147][148] persists if, for instance, both q 2 q 2 ⊥ m 2 and 1 q 2 ⊥ /q 2 b 2 . Indeed, in this case the running QED coupling also becomes extremely small α scr (q q q) ∼ αq 2 /(q 2 ⊥ b 2 ) 1. Also, electrons and positrons with energy larger than m > 2m e , i.e. with q 2 ⊥ m 2 q 2 and q 2 /m 2 b 2 might undergo the quantum triviality since α scr (q q q) ∼ αq 2 /m 2 * 1. Conversely, for q 2 ⊥ /m 2 1 b 2 q 2 /m 2 or 1 b 2 q 2 ⊥ /q 2 m 2 /q 2 the screening effects are negligibly small and the running QED coupling approximates to it standard value α scr (q q q) ∼ α.
Now we proceed to investigate how the running QED coupling modifies the Coulomb potential. To this end, we write Eq. (37) in an equivalent form where the original integration variables have been changed to cylindrical ones, and the integral representation of the Bessel function of order zero has been used [149]. Here ϕ (φ) is the polar angle associated with q q q ⊥ (x x x ⊥ ) and the square of the momentum q q q must be understood as q q q 2 = q 2 ⊥ + q 2 . The integration over q can be performed by using the residues theorem. For this, we point out that the integrand in Eq. (39) has four poles on the imaginary axis: with the positive roots taken for x < 0 and the negative ones for x > 0. As before, m * = m √ 1 + b 2 with b = g B/m is the ALP mass dressed by the external field [see below Eq. (17)]. With these details in mind, and after developing the change of variable u = 2bq ⊥ /(m * √ 1 + b 2 ), it turns out that In this formula, a λ * = q/(4πλ * ) is the Coulomb potential evaluated at the effective wavelength λ * = m −1 * of an ALP.
Depending on whether the magnetic field strength is smaller or larger than the characteristic scale of the theory B c = m/g, this quantity will reduce either to the Compton wavelength of the axion λ = m −1 or to a sort of "Larmour" length , respectively. As the framework under consideration allows for exploring distances larger than the natural length scale of the theory g ∼ Λ −1 UV [111], the effective wavelengths of interest are those for which the condition λ * g is fullfilled. Observe that, in the strong field limit [b 1], the latter translates into an upper bound for the largest magnetic field that can be considered in the theory B max = g −2 ∼ Λ 2 UV , corresponding to a parameter b max = (mg) −1 with 1 mg. Equation (41) has been written in terms of the dimensionless coordinatesx ⊥, = |x x x ⊥, | 2λ * 1 + 1 b 2 and functions which grow monotonically with the increasing of u: For large values of their argument [u 1], they behave as whereas for small values of u 1 they approach The formula found above [Eq. (41) with (42)  In the left panel [|x x x | = 0], the curves tend to stick closer to each other, falling down rather sharply as the magnetic field grows gradually and the distance from the particle location is much smaller than the Compton wavelength of an ALP [|x x x ⊥ | λ]. Observe that the red dotted curve exhibits a region in which the Coulomb interaction vanishes practically. This flattening is a manifestation of QED triviality due to the dominance of AED at asymptotically large B [read also the discussion below Eq. (38)]. Conversely, when |x x x ⊥ | = 0, In such a situation, the integrand that remains combines J 0 (x ⊥ u) with a function that tends asymptotically to 1 with the growing of u, whatever be the value of b. Hence, the whole integrand reaches the maximum value as the argument of the Bessel function J 0 (x ⊥ u) satifies the conditioñ x ⊥ u 1. Following this analysis we infer that, when |x x x ⊥ | is small on the scale of λ * /(1 + 1/b 2 ) 1 /2 , the main contribution to the integral in Eq. (41) results from the region in which 1 u x −1 ⊥ , provided b |x x x ⊥ |/λ. As a consequence, one can exploit Eq. (43) and integrate over u by using formula (6.616.1) of Ref. [150]. 12 With these details in mind, 12 The integral tabulated in Ref. [150] reads the modified potential approaches This formula looks like a pure Yukawa potential in which the scale ∼ λ * (1 + 1/b 2 ) − 1 /2 characterizes the short-range behavior. This kind of dependence is somewhat expected, since the second propagation mode coexists with a massive branch determined by the ALP sector.
is independent of the ALP mass and 1 + 1/b 2 ∼ 1. The resulting expression can be read off directly from Eq. (45), and tends to be tiny as b grows progressively for any nonzero distance |x x x ⊥ | λ * . The sharp trend exhibited by the curves in the left panel of Fig. 5, nearby the particle location [|x x x ⊥ | = 0], can be understood qualitatively by considering the generalized limit of Eq. (45) when b → ∞: Footnote 12 continued Changing the variable to t = (x 2 + 2γ x) 1/2 , it becomes ∞ 0 dt t which is the desired representation for our calculations.
Here, Ei(z) = − ∞ −z dt t −1 e −t is the exponential integral function [149]. With this result at hand, Eq. (45) becomes a magnetic field independent delta potential Although this result explains why the curves tend to collapse toward the vertical axis with the growing of b [see Fig. 5], and so the trivial interplay for |x ⊥ | = 0, it has to be considered with caution. This is because the b parameter is actually limited from above by b max and x ⊥ must be larger than g [see discussion below Eq. (41)]. It is worth pointing out that, in the limit of b → b max and |x x x ⊥ | → g the axion-Coulomb potential is regular and reads lim b→b max When x x x ⊥ = 0 0 0, the Bessel function J 0 (x ⊥ u) = 1, and the exponentials which remain in Eq. (41) decrease very fast for large values ofx u. Hence, for distances |x x x | λ * (1 + 1/b 2 ) − 1 /2 , the main contribution to the integral in Eq. (41) results from the region in which 1 u x −1 . Once again, we exploit Eq. (43) and keep the factor u/ √ 1 + u 2 as it stands. These steps reduce the integral to one in which (6.616.1) of Ref. [150] can be used. Consequently, The outcome above applies whenever b |x x x |/λ. In the limit of b 1, this expression approaches We note that, as the exponent involved in this formula is much smaller than unity [m * |x x x | 1], the screening induced by the vacuum polarization does not deviate Eq. (49) from the Coulomb potential a 0 (0, This fact agrees with the behavior shown in the right panel of Fig. 5, which exhibits that -for distances smaller than the ALP wavelength -the behavior of the curves is very similar to the one associated with Coulomb's law. The limit of |x x x | → g is readable from Eq. (49) and gives a 0 (0, g) ≈ q/(4πg).
Observe that this formula differs from the standard Coulomb law by a factor 1/2. Thus, the behavior of the axion-modified Coulomb potential in the longitudinal direction saturates in the sense that -unlike the direction perpendicular to B B B [see Eq. (45)] -it reaches a universal shape, independent of the external field strength. This fact helps to understand why the curves in the right panel of Fig. 5 converge to the red one, the behavior of which differs from the case x x x = 0 [left panel] as b 1. The procedure for obtaining Eq. (48) can be extended easily to the case in which |x x x ⊥ | λ * /(1 + 1/b 2 ) 1 /2 . This gives rise to the following short-range anisotropic potential where |x x x| = (|x x x ⊥ | 2 + |x x x 2 ) 1 /2 and m * ≈ mb. In the limit of After applying once again the integral (6.616.1) of Ref. [150] the modified Coulomb potential reads This formula is an anisotropic Coulomb's law which tends to vanish as the external field grows unless |x x x ⊥ | = 0. It resembles closely the longe-range behavior found within a pure QED context, when the external magnetic field exceeds the characteristic scale associated with this framework B 0 = m 2 e /e ≈ 4.42 × 10 13 G [91,92,94]. However, in QED, b 2 is replaced by the factor αb/(3π) [see first line in Eq. (13)]. Observe that, at |x x x ⊥ | = 0, Eq. (52) reduces to a pure Coulomb 1. In such a case, the part of the integrand which does not include the Bessel function can be Taylor expanded in u. As a consequence, the long-range behavior of the axion-Coulomb potential, perpendicular to the direction of the field, reads [b 1] The axion-modified Coulomb law should affect, first of all, the field of an atomic nucleus, placed in a magnetic field. We now show that, although the potential along the magnetic field [b 1] resembles the ∼ 1/|x x x | law approximately, the ground-state energy of a nonrelativistic electron of a hydrogen atom: 13 remains bounded from below. Observe that, the singularity is here cut off by the Larmor length L B = (m e √ b) −1 , whereas a B = (m e α) −1 is the Bohr radius. Both length scales depend on α which is not affected substantially in the considered regime. As a consequence of integrating out x , we end up with the expression of Elliott and Loudon [95] [q = e]: Prior to the inclusion of the QED screening due to the vacuum polarization tensor [91,92], it was commonly accepted that Eq. (54) is unbounded from below when the magnetic field grows unlimitedly. With its reappearance one could think that a possible dominance of AED over QED restores this problem. However, as long as AED is combined with QED, the resulting framework is nonrenormalizable. This means that the above expression would provide reliable predictions whenever the Larmor length exceeds the characteristic scale of AED [L B g]. This bounds the field parameter b = B/B 0 (m e g) −2 and simultaneously limits the ground state energy ε 0 > −2α 2 m e ln 2 1 m e αg , ruling out the reinstatement of the aforementioned problem.

Conclusion
We have considered the interaction between axion-like particles and photons inside the background of a superstrong constant magnetic field. Only the extraordinary photon mode -mode-2 -whose electric field lies in the plane spanned by the external magnetic field and the photon momentum interacts with the ALP , while the ordinary mode -mode-3remains the same as it was in strong-magnetic-field QED. We 13 Valid in a shallow well potential have obtained that the axion contribution to the polarization operator grows quadratically with the magnetic field, while the QED part of the polarization operator of the extraordinary mode is known [118][119][120][121] to grow only linearly (while the ordinary-mode part grows much more slowly). In this ALP-dominating regime we established the common ALPextraordinary-photon dispersion law determining their mass shell. It consists of two nonintersecting branches, with no mixed state forming -in contrast to QED, where the photon mixes with the electron-positron state to form a polariton. One branch is massless in the sense that the frequency as a function of spatial momentum turns to zero when all components of the latter are zero. This implies that the mass defined as the rest energy is zero. This branch is associated with the photon modified by its interaction with an ALP. The other branch is massive and it is attributed to an ALP, modified by its interaction with the photon. The massless branch exhibits the important feature of flattening for components perpendicular to the B B B direction which are smaller than the ALP mass. This feature leads to strong birefringence, because the dispersion curve of the mode-3 photon, which remains the same as in QED, does not show such feature. The flattening results in bending the direction of the wave envelope propagation (the group velocity) towards the external field and in the capture of mode-2 photons by curved lines of force of the magnetic field of pulsars. Contrary to the analogous effect known in QED, in the present scenario not only γ -quanta of curvature radiation undergo the capture, but also heat photons in the X-ray range. Under certain circumstances the capture of gamma-quanta is more efficient in AED and its potential consequence might be the destruction of the well established mechanism of formation of pulsar radiation. We have used this incompatibility to exclude a certain region in the axion parameter space. Furthermore, when applied to the softer thermal radiation of the pulsar surface the capture causes that the mode-2 X-ray quanta all gather together to peak near the direction of the magnetic field, while beyond this direction all photons are those of mode 3, i.e. they are polarized orthogonal to the plane spanned by the magnetic field and the observation direction. Their angular distribution is determined by the refraction at the border of the strong magnetic field area, studied in the paper, and by the mechanisms of formation of the soft radiation (see their review in Ref. [132]).
Despite the charge-neutrality of quantum vacuum fluctuations of axion like fields, we have seen that the presence of a magnetic field whose strength overpasses the characteristic scale of AED could induce a strong screening of static electric charges. This screening depends on the field strength, and for certain energy regimes the running coupling could be screened almost to zero, making the QED building blocks very weakly interacting between each other. Moreover, the dependence on the magnetic background promotes an aniso-tropic behavior in the Coulomb potential. Indeed, at distances from the source much smaller than the Compton wavelength of the ALP, and for magnetic field strengths larger than m/g, the Coulomb potential of a point-like charge is replaced by a Yukawa-type potential in the direction perpendicular to the magnetic field, while along the field the Coulomb law is preserved in a modified form at any length scale. We find that at unlimitedly large magnetic fields the longstanding problemovercome in QED -that the ground-state energy of a hydrogen atom is unbounded from below, seems to be reinstated in AED. However, the nonrenormalizable feature of this theory rules out this possibility.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The calculations in this study have been carried out by analytical means; the data shown in the figures can be reproduced from the formulas given in the corresponding sections.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: Higher order contributions in strong fields [b 1]
This appendix is devoted to precise that, in the strong field regime [B B c ], the polarization tensor of AED [P μν ∝ b 2 = (B/B c ) 2 ] is not affected substantially by higher-order corrections to the diagram depicted in Fig. 1. To this end, we will refer to the nonlocal theory described by Γ [a] in Eq. (6). Bearing in mind that only one-particle-irreducible diagrams may contribute to the inverse two-point Green's function of photons, we find the two next-to-leading order diagrams depicted in Fig. 6. The precise structure of the involved vertices is not necessary to understand the following arguments. The only important property of them to keep track of is that the three-point nonlocal vertex is proportional to b, while the one connecting four points does not depend on the external field. Another important feature of this framework is that the photon propagator involved in each of these diagrams is the one modified by the leading order contribution of the polarization tensor in Eq. (8). This axion-modified photon propagator has been determined explicitly in Eqs. (36), and behaves as D μν (q) ∼ b −2 when b → ∞. Hence, at very strong field strengths [b 1], the two one-loop pieces in Fig. 6 turn out to be suppressed by a factor ∝ b −2 as compared with P μν ∝ b 2 [see Eqs. (10), (11) and (14)].
We note that the insertion of the nonlocal three point vertex comes accompanied with an axion-modified photon propagator at least. This feature prevents the growing of higher order contributions with the field strength, and so, protects the perturbative character of each contribution in the expansion. For instance, at two-loop level, there are eight Feynman diagrams which contribute to the irreducible two-point function of photons [see Fig. 7]. With the exception of the last graph -which is ∝ b −8 -the remaining six contributions are suppressed by a factor ∝ b −6 when b 1. Clearly, this pattern also manifests in irreducible vertices involving more than two points, a fact that has been exploited implicitly in our two-loop analysis since the establishment of some of these diagrams relies already on the one-loop correction of the three-point nonlocal vertex, for instances. This consideration justifies the applicability of the pursued approach within pure AED up to the unitarity limit of B max = Λ 2 UV ∼ g −2 , or b max = (mg) −1 , imposed by the nonrenornalizability.

Appendix B: Degrees of freedom in AED
In QED, there are two photon degrees of freedom that propagate in the magnetized vacuum with different dispersion laws [116]. This property of birefringence is a clear manifestation of the anisotropy induced by the external field via the polarization tensor. Let us investigate a similar issue within the framework of AED. Primarily, this theory contains three degrees of freedom: two (massless) photon degrees and one where q (2) 0 = q 0− [see Eq. (20)]. The other two values are independent of the ALPs presence and are determined by the dispersion laws of mode-1, which in the loop approximation of QED coincides with the vacuum dispersion law q (1) 0 = |q q q|, and of mode-3 whose dispersion law is of the general form q 1 2 , the same as q 0− [see Eq. (20)]. The eigenvectors μ j are chosen to differ -in the rest frame of the magnetic field -from the vectors (12) by a factor: μ j = μ j /(q ⊥ B 2 ). The division by q ⊥ is crucial, as we shall see, to exclude mode-1 and to gather conveniently the two massless degrees of freedom in modes-2 and -3. Note that the set of eigenvectors¯ j remains an orthogonal basis also when the background field is absent. Therefore, these formulae may be used, too.
Then the electric e e e j ∼ q 0¯ j − q q q¯ 0 j and magnetic h h h j ∼ q q q ×¯ j fields of each mode behave as: e e e 1 ∼ −n n n ⊥ q 0 q 2 , h h h 1 ∼ − n n n ⊥ × q q q q 2 , e e e 2⊥ ∼ n n n ⊥ q 2 , e e e 2 ∼ q q q q ⊥ (q 2 − q 2 0 ), h h h 2 ∼ q 0 q q q × n n n ⊥ , e e e 3 ∼ −q 0 n n n ⊥ × q q q , h h h 3⊥ ∼ −n n n ⊥ q 2 , h h h 3 ∼ q q q q ⊥ , (B.3) Here n n n ⊥, = q q q ⊥ /|q q q ⊥ | is a unit vector along the perpendicular direction of B B B. It may seem that e e e (2) in Eq. (B.3) is singular in q ⊥ = 0. This is not the case, however, because if q ⊥ = 0 also q 2 − q 2 0 = 0 for a massless branch [see Eq. (18)]. Now, in the vacuum with no background magnetic field there are two polarizational degrees of freedom, and these are modes 3 and 2. The mode-3 is identically transverse, e e e 3 · q q q = h h h 3 · q q q = e e e 3 · h h h 3 = 0, whereas the mode-2 becomes transverse under the vacuum dispersion law: e e e 2 · q q q = q 2 q q q 2 − q 2 0 /q ⊥ = 0 if q 2 0 = q 2 + q 2 ⊥ . In this case, the relations h h h 2 · q q q = e e e 2 · h h h 2 = 0 are fulfilled identically. As for mode-1, it is seen from Eq. (B.3) that it nullifies in the vacuum as long as the dispersion law is q 2 = 0.
In the magnetic field, mode-2 is the so-called extraordinary wave, whose longitudinal electric component is generally nonzero, e e e 2 · q q q = 0, since now the dispersion law implies q 2 = 0. The exclusion is provided in the background magnetic field by the momentum parallel to the latter: q q q ⊥ = 0 0 0. In this case the longitudinal component of mode-2, e e e 2 · q q q = q 2 (q 2 − q 2 0 )/q ⊥ → 0 due to the dispersion law q 2 0 − q 2 = f (q 2 ⊥ ) [see Eq. (28)], since for the massless branch the function f (q 2 ⊥ ) turns to zero as fast as q 2 ⊥ in QED or faster than q 2 ⊥ in the ALPs theory according to Eq. (21).
The identically transverse mode-3 continues to represent the second degree of freedom in the magnetic field, too. As for mode-1, in the background magnetic field, at least as long as the one-loop approximation of QED is used, its dispersion law remains q 2 = 0 [116]. Therefore the electric and magnetic fields in mode-1 vanish according to Eq. (B.3), leaving again only two polarizational degrees of freedom on the mass-shell q 2 = j . The vanishing of mode-1 may more directly be seen already from the expression for μ 1 in Eq. (12). The four-gradient term ∼ q μ in it gives no contribution to the field strength, whereas the q 2 -term nullifies thanks to the dispersion law q 2 = 0. Therefore, with the choice of the amplitude of the free waves fixed as in Eq. (B.2), mode-1 remains unphysical for every direction of propagation. Because of this reason, this mode is not considered throughout the paper.