Electric and magnetic axion quark nuggets, their stability and their detection

The present work studies the dynamics of axion quark nuggets introduced in Zhitnitsky (JCAP 0310:010, 2003) and developed further in the works (Zhitnitsky in Phys Rev D 74:043515, 2006; Lawson and Zhitnitsky in Phys Lett B 724, 17, 2013; Lawson and Zhitnitsky in Phys Rev D 95:063521, 2017; Liang and Zhitnitsky in Phys Rev D 94:083502, 2016; Ge et al. in Phys Rev D 97:043008, 2018; Zhitnitsky in Phys Dark Univ 22:1, 2018; Lawson and Zhitnitsky in Phys Dark Univ 100295, 2019; Raza et al. in Phys Rev D 98:103527, 2018; Fischer et al. in Phys Rev D 98:043013, 2018; van Waerbeke and Zhitnitsky in Phys Rev D 99:043535, 2019; Flambaum and Zhitnitsky in Phys Rev D 99:043535, 2019; Lawson and Zhitnitsky in JCAP 02:049, 2017; Ge et al. in Phys Rev D 99:116017, 2019). The new feature considered here is the possibility that these nuggets become ferromagnetic. This possibility was pointed out in Tatsumi (Phys Lett B 489:280 2000) for ordinary quark nuggets, although ferromagnetism may also take place due some anomaly terms found in Son and Zhitnitsky (Phys Rev D 70:074018, 2004), Son and Stephanov (Phys Rev D 77:014021, 2008) and Melitski and Zhitnitsky (Phys Rev D 72:045011, 2005). The purpose of the present letter however, is not to give evidence in favor or against these statements. Instead, it is focused in some direct consequences of this ferromagnetic behavior, if it exists. The first is that the nugget magnetic field induces an electric field due to the axion wall, which may induce pair production by Schwinger effect. Depending on the value of the magnetic field, the pair production can be quite large. A critical value for such magnetic field at the surface of the nugget is obtained, and it is argued that the value of the magnetic field of Tatsumi (2000) is at the verge of stability and may induce large pair production. The consequences of this enhanced pair production may be unclear. It may indicate that the the nugget evaporates, but on the other hand it may be just an indication that the intrinsic magnetic field disappears and the nuggets evolves to a non magnetized state such as in Zhitnitsky (2003), Oaknin and Zhitnitsky (Phys. Rev. D 71:023519, 2005), Zhitnitsky (2006), Lawson and Zhitnitsky (2013), Lawson and Zhitnitsky (2017), Liang and Zhitnitsky (2016), Ge et al. (2018), Zhitnitsky (2018), Lawson and Zhitnitsky (2019), Raza et al. (2018), Fischer et al. (2018), van Waerbeke and Zhitnitsky (2019), Flambaum and Zhitnitsky (2019), Lawson and Zhitnitsky (2017), and Ge et al. (2019). The interaction of such magnetic and electric nugget with the troposphere of the earth is also analyzed. It is suggested that the cross section with the troposphere is enhanced in comparison with a non magnetic nugget but still, it does not violate the dark matter collision bounds. Consequently, these nuggets may be detected by impacts on water or by holes in the mountain craters (Pace VanDevender et al. in Sci Rep 7:8758, 2017). However, if the magnetic field does not decay before the actual universe, then this would lead to high energy electron flux due to its interaction with the electron gases of the Milky Way. This suggests that these magnetized quarks may be a considerably part of dark matter, but only if their hypothetical magnetic and electric fields are evaporated.

a e-mails: firenzecita@hotmail.com; osantil@dm.uba.ar (corresponding author) b e-mail: sempe_100@hotmail.com be unclear. It may indicate that the the nugget evaporates, but on the other hand it may be just an indication that the intrinsic magnetic field disappears and the nuggets evolves to a non magnetized state such as in Zhitnitsky (2003) The interaction of such magnetic and electric nugget with the troposphere of the earth is also analyzed. It is suggested that the cross section with the troposphere is enhanced in comparison with a non magnetic nugget but still, it does not violate the dark matter collision bounds. Consequently, these nuggets may be detected by impacts on water or by holes in the mountain craters (Pace VanDevender et al. in Sci Rep 7:8758, 2017). However, if the magnetic field does not decay before the actual universe, then this would lead to high energy electron flux due to its interaction with the electron gases of the Milky Way. This suggests that these magnetized quarks may be a considerably part of dark matter, but only if their hypothetical magnetic and electric fields are evaporated.

Introduction
The existence of lumps of quark matter has been postulated long ago [17,107]. This is not to be confused with nuclear matter, which is composed by a large number of protons and neutrons whose main interaction is due to nuclear forces. Instead, quark matter is approximately a Fermi gas composed of 3N B quarks constituting a color singlet baryon of baryon number N B , and their interaction is much weak, due the fact of the quotient of the cross section to the mass of the object σ n /M is very small. Some of these states of quark matter are composed only by u and d quarks, and these are known as non strange quark matter. Other form is strange matter, which contains s quarks as well as u and d quarks, in such a way that flavor equilibrium is established by weak interactions d → u +e+ν, s → u +e+ν and s +u → d +u. Strange matter also contain gluons, and a small component of electrons whose role is insure electric charge neutrality [107]. These strange lumps may be formed during an hypothetical phase transition of the early universe from a quark gluon plasma state to a hadron phase. There is debate about the nature of this transition, it can be a first order or second order one, a crossover or even an spinodal type of transition. But there are segregation scenarios which predicts the formation of these lumps, irrespective on the order of the phase transition. An example is given in the reference [11].
Strange form of matter has received special attention, since it is conjectured to have lower energy per baryon number than ordinary nuclei, thus it may constitute a stable state [15,21]. In other words, the energy per baryon is less than 930 MeV, at least at very low temperatures. An heuristic argument for stability comes from the fact that a non strange quark lump posses a Fermi momentum p f ∼ 300-350 MeV, while the strange quark mass is around m s ∼ 80-130 MeV. As the Fermi momentum is very high in comparison, it may be energetically favored for some non strange quarks to become strange, as this conversion may lower the energy and the Fermi momentum of the system. The estimations given in [15,21] suggest that the inclusion of strange matter decrease the energy per baryon by  MeV. The range of baryon numbers allowed for these objects is 10 2 < N B < 10 57 [5]. The upper bound is due to the fact that higher baryon numbers collapse into a black hole. The lower bound arises due to shell effects which raise the energy per baryon.
The bounds 10 2 < N B < 10 57 have been further constrained. By use of detailed balance arguments, it is found that primordial nuggets with N B < 10 52 evaporate. However [69] suggested that nuggets with N B > 10 46 survive hadron emission, and lower charge nuggets also survive due to reabsortion of hadrons. Another constraint comes from nucleosynthesis. The point is that, if quark nuggets exist and are stable, then iron may not be the ground state of nuclear interactions. This is a delicate aspect when studying nucleosynthesis and was initially considered in [10]. For instance, for an inflationary universe with = 1, the prediction of the density of the 3 He and D elements is below the observed value. This problem may be solved if there is some unknown mechanism of photodisintegration of 4 He into this elements. This requires high energy photons which may be produced by decay of massive neutrinos by non standard interactions. Such channels may be sourced by the presence of quark nuggets through a weak interaction of high order, as in this case the decay rate will be close to zero. Despite these observations, it was suggested in [95] that quark nuggets with 10 16 < N B < 10 22 are con-sistent with a universe with = 1 and with the predicted abundances of the light elements. This motivated a special interest for studying the cosmological consequences of these objects, in particular, in dark matter applications.
The mass density of a quark nugget is ρ ∼ 35. 10 16 kg/cm 3 , and its mass is approximately 6m = N B GeV. The surface electron cloud for nuggets with N B < 10 15 extend to a distance r ∼ 10 −8 cm from the surface while for very high baryon numbers this distance can be around 400 fm. Such nuggets therefore possess a Coulomb barrier, with a value close to 10 MeV. Due to this barrier, only neutrons at low energies may be absorbed by the nugget. There is also a probability of emission of a neutron. Further stability issues were considered in [5,6,14,65,68,69,[72][73][74]80,101], and it is believed that stable nuggets have a mass in the range 10 −8 kg< m n < 10 20 kg. Thus, very large quark nuggets may masses of planetary order.
Several years ago it was suggested that quark nuggets may reach a highly magnetized final state [103]. This observation is partially motivated from the Bloch works about the possible ferromagnetism of an electron gas [16]. This hypothesis was supported in [57], who showed that for an electron gas in a fully polarized state the ferromagnetic state is stable. These claims were subsequently tested for a quark liquid in [103] and some evidence supporting a ferromagnetic state was collected for quark nuggets. The presence of an intrinsic magnetic field is an attractive feature, as it may make the nugget more stable if B < 10 16 T [19]. A problem is that such highly magnetized magnetic fields may violate the collision requirements for dark matter. However, it was shown in [81] that due to the fall of the magnetic field as R −3 , the resulting cross section is acceptable for representing dark matter.
The present work is devoted to the effect of such magnetic fields in a different type of quark nuggets, which were introduced in [109], in which the quarks are trapped in the bulk of an axion domain wall. There are several theoretical reasons for which these axion quark nuggets are interesting. First, the axion particle is an attractive candidate for solving the CP problem in QCD [27,51,82,83,96,108]. However, the formation of axion domain walls is problematic, since they generate an energy density which largely overcomes the critical density. Nevertheless, there are some axion models for which the axion domain structure does not possess this problem [97]. These domain walls induce a large value of baryon number due to the fermions living on the wall, which can be shown by bosonization techniques [63]. Due to the surface tension, these wall tend to contract and segregate baryons and anti-baryons. It is likely that the contraction of these objects stops when the surface tension compensate the pressure difference and the resulting bubble enters into the CS phase [7,8]. These nuggets may be formed without relying on the order of the phase transition from the quark gluon plasma to the confined phase.
The organization of the present work is as follows. In Sect. 2.1 the modified Maxwell equations for magnetized quark nuggets in presence of an axion domain wall are considered. It is argued that neither the magnetic field found in [103] nor the axion wall profile are considerably deformed by the interaction between each other. The resulting electric field in the nugget is then estimated. In Sect. 2.2 the pair creation by Schwinger effect in the wall is calculated by use of several methods, and it is argued that the magnetic field [103] is at the verge of stability. In Sect. 3 the fate of a quark nugget impacting in the earth is analyzed and its possible detection is specified. It is also argued that if these nuggets had a density close to dark matter then they would generate a high energy flux due to their interactions with the Milky Way electrons. Section 4 contains the interpretation of the results.

The equations describing the configuration
If the arguments of [103] about a possible ferromagnetic nature of these nuggets are taken into account, then there should be a magnetic field B due to bulk effects inside the object. For standard nuggets, the reference [103] presents a magnetic field that scales as Strictly speaking, a nugget does not have a spherically symmetric radial magnetic field. In fact it may depend on the polar ψ angle in spherical coordinates. The value (2.1) represents a extrema value say, at certain north or south pole drawn on the nugget. This extreme numerical value at the surface of the nugget r n is estimated as B 0 ∼ 10 12±1 T or B 0 ∼ 10 16±1 Gauss [103]. For ordinary quark nuggets, the work [19] suggests that such large magnetic field increase the stability of the configuration, if it does not reach a critical value of the order B ∼ 10 16 T. Consider now a magnetized quark nugget as above, but including an axion domain wall [35][36][37][38]58,[60][61][62][63]79,90,105,[109][110][111]. The domain wall profile function is denoted by a(r ), and interpolates between two axion vacuums a 0 and a 0 + 2π f a N with f a the axion constant, whose standard values run between 10 9 GeV < f a < 10 12 GeV [97]. The value of N is usually not far from unity. The choice of the window of values for f a has been justified previously as follows. Below the QCD temperature T qcd ∼ 100-150 MeV, there appears an induced periodic potential U (a), and the axion becomes light but massive. A customary assumption is that the axion is at some value of the order of the top of the potential U (a) at the time where this transition occurs. When the Hubble constant is of the same order as the axion mass this pseudo scalar falls to the potential minimum and starts coherent oscillations around it. The initial amplitude, which correspond to a maximum, is A ∼ f a and thus, the energy stored at by these oscillations is of the order E ∼ A 2 m 2 a . The authors of [1,26,86] analyzed the evolution of these oscillations to the present universe and found that the axion energy density today would be larger than the critical one ρ c ∼ 10 −47 GeV 4 unless the bound f a < 10 12 GeV takes place. On the other hand, there are phenomenological observations which fix this scale f a > 10 9 GeV [24]. This lower bound is required for suppressing the power radiated in axions by the helium core of a red giant star to the experimental accuracy level. Current algebra methods [12] estimate the axion mass by the relation This fixes the axion mass in the window 10 −6 eV< m a < 10 −3 eV. The axion potential U (a) is given in these terms by Here n takes integer values, which depends on the particular axion model in consideration. It should be emphasized however, that the previous discussions holds for the standard QCD axions [51,82,83,96]. There exists several alternatives to this picture. An example are axion models which try to avoid the isocurvature problems [66,67,84], examples can be found in [13] and references therein. Nevertheless, the present work is focused in the standard QCD axion [51,82,83,96], as the axion nuggets in [35][36][37][38]58,[60][61][62][63]79,90,105,[109][110][111] seems focused in these type of Goldstone pseudo-bosons.
The explicit form of the function a(r ) describing the axion domain wall is not too relevant in the following discussion. The important point is that a(r ) evolves from one vacuum to the other in a radial distance of the order the last factors take into account the axion mass bound 10 −6 eV< m a < 10 −3 eV, which is assumed to hold in the present work. The baryon number range that the references [35][36][37][38]58,[60][61][62][63]79,90,105,[109][110][111] employ is 10 23 < N B < 10 32 . For the larger baryon number one has 10 −6 cm< R < 10 −3 cm, while for the smallest value 10 −9 cm< R < 10 −6 cm. The mass bound is given by 10 −6 kg< m n < 10 5 kg, thus these objects can be quite heavy. The value N B ∼ 10 25 is attractive in the context of the physics of the 511 KeV line emission from dwarf galaxies [61], but it will be assumed that the whole range is valid.
In presence of an axion field a(x, t), when a magnetic field is turned on, an electric field is induced and viceversa. This follows from an inspection of the Maxwell equations in presence of an axion field, which are given by The axion dependent terms in these equations arise due to the interaction L a = g aγ γ a E · B between the axion and the electromagnetic field. Here g aγ γ = cα/π f a with α = 1/137 the fine structure constant. The constant c is model dependent, but its value is not far from the unity. It is seen from (2.4) that the quantity ρ a = −g aγ γ ∇ · (a B) can be interpreted as a charge source for the electric field, not present in the standard Maxwell theory. Thus, when a magnetic field passes through an axion wall, an electric field is induced due to this source. In addition, U (a) is the axion potential defined in (2.3). From the discussion given above it follows that the magnetic field (2.1) is modified due to the Eq. (2.7). Furthermore, an electric field E is induced, as is seen from equation (2.4). It is explicitly given by (2.9) The presence of the electric field may generate electron positron pairs due to Schwinger effect. The stability of the configuration is then compromised if the resulting electric field (2.9) is very strong. Clearly, the axion profile a(r ) of references [35][36][37][38]58,[60][61][62][63]79,90,105,[109][110][111] should be modified in presence of the magnetic field B, which is also deviated from its original form (2.1). In fact, the Eq. (2.5) shows that ∇ × E = 0, as the magnetic field B is static. In other words, the electric lines for the configuration are never closed. This may be non true if ∇ × (a B) = 0, as seen from (2.4), unless the profile a and the magnetic field B(r ) are adapted for this to happen. Thus, the task of solving the explicit profiles may be a complicated one. For this reason, it will be assumed below that the axion domain wall function a(r ) keep its shape approximately and that the magnetic field B(r ) do not deviate considerably from (2.1). But before considering the decay of the induced electric field by Schwinger effect, it may be convenient to justify these approximations first.
Consider first the approximation that the magnetic field does not deviate considerably from the functional dependence (2.1). The term J a = g aγ γ ∇a × E is the axion induced current in (2.6). If one assumes that the field a make a sudden change of a → a + 2π f a at the nugget surface, then the axion surface magnetization is proportional to M ∼ 2π f a g aγ γ E × n with n the surface unit normal. The induced magnetic field at the surface is of the order |B a | ∼ |M|. But |M| ∼ f a g aγ γ |E| and, by taking into account that E = −g aγ γ a B, it follows that |B a | ∼ |M| ∼ f 2 a g 2 aγ γ |B| ∼ α 2 |B|. This suggest that the magnetic field B a induced by the axion current is around four orders of magnitude smaller that the internal magnetic field (2.1), and thus the axion current does not deviate the value (2.1) considerably.
Another heuristic argument for this statement about B a comes from the study of the magnetic field of an infinite wire with radius r w ∼ m −1 a near its border. This magnetic field is B w ∼ μ 0 J w /r c with r c the standard cylindrical radius. If the field value in the border of the wire is near B = 10 12±1 TeV, then its current should be J w ∼ 10 12±1 TeVm a . On the other hand, the axion current is J a = g aγ γ ∇a × E, this follows from (2.6). If the domain wall a is not considerably modified, it should have a size of length L a ∼ m −1 a and therefore, the current is of order |J a | ∼ g aγ γ am a E. But the Eq. (2.4) shows that E = −g aγ γ a B and therefore Here the fact that a ∼ 2N f a π has been taken into account, with N an integer not far from unity. This current is fourth orders of magnitude larger than J ∼ 10 12±1 TeVm a , which suggest that the induced magnetic field due to the axion wall is of smaller order than the magnetic field (2.1) itself.
The estimations given in the last two paragraphs assume that the profile a(r ) is not strongly modified in presence of a magnetic field. To see that this may be the case, one should check that the term E·B is not considerably larger than ∂ a U in (2.8), otherwise the axion domain wall will be considerably deformed by this new term. Since the potential (2.3) can be approximated for small values of a as ∂ a U ∼ m 2 a a and g aγ γ E · B ∼ g 2 aγ γ a B 2 , with B ∼ 10 11±1 T, one has to check that But m 2 π f 2 π = 10 −4 GeV 4 and 10 22±2 T 2 ∼ 10 −10±2 GeV 4 , and thus the last inequality is true. Therefore, the assumption that the magnetic field B do not deviate from its form (2.1), and that the electric field E does not induce a considerable deformation of the domain wall is, at least, consistent. This is the approximation to be used below.
The field (2.1) is induced by a ferromagnetic behavior of a quark liquid at the bulk of the nugget. On the other hand, the axion wall can be a source of a magnetic field as well. This is understood from an analysis of anomaly terms for axions [76,98,99], in presence of a non zero chemical potential, which is the case for a nugget in the Color Superconducting phase [7,8]. There are two anomaly terms, one is the standard one L 1 ∼ a F μν F μν , which induces the decay a → γ + γ by the ABJ anomaly diagram. However, in presence of a non zero chemical potential μ there is a further WZW type of anomaly term given by [76,98,99] Here B is an external magnetic field, and C a model dependent constant, but not far from unity. From the fact that ∇ · B = 0 it is seen that this term is a total derivative, thus it does not contribute to the equations of motion. However, the presence of this term is not trivial. To see this, consider a flat axion domain wall. Then this term induces an extra energy contribution proportional to μB A with A the area of the wall. This energy is produced by a magnetic moment by unit area located at the wall, given by For a spherical wall, the total magnetic moment may be zero, since the gradient contributions of a at opposite points are equal and of opposite sign. This implies that the leading contribution to the external magnetic field is the quadrupole one. However, for a simple estimation, one may remember that for a magnetized sphere at the south and north poles the magnetic field is |B| = μ 0 |M|, with M its uniform permanent magnetization. If this value is employed as a guide for a magnitude order, then it is mandatory to estimate the value of M. By a naive application of formula (57) of [98] one finds that Here is a gap characterizing the formation Cooper color pairs, by assuming that the evolution of the object is such that the CFL phase is achieved [7,8]. A typical value for the chemical potential is 1 GeV and ∼ 50 MeV. This implies that the magnetic field these references predict may be of the order of B ∼ 10 12 T, which is pretty close to the value of [103]. The physics describing these fields is however different.
In view of the present discussion, a generic magnetic field will be considered below without relying wether it is induced by the axion wall or by the quark liquid inside the nugget. The working assumption is then that an electric field E = −g aγ γ a B is induced, and the task is to estimate its decay probability due to electron positron pair production.

Schwinger effect and estimation of the critical magnetic field of the configuration
In order to study the electron positron creation by the induced electric field E, one may avoid the complication of the inhomogeneity of the electric field E and simply assume that E ∼ g aγ γ a 0 B 0 with B 0 the surface magnetic field of the bubble. This is the roughest possible approximation. In this situation, the vacuum persistence probability is given by Here the rate of pair creation w by a constant electric field E is expressed by the well known Schwinger formula This formula implies that the vacuum transition probability can be written as follows The approximate volume occupied by the field is roughly (2.11) If the critical field is reached, then evaporation in electron positron pairs may be relevant. This will be the case when Taking into account that |E| ∼ g aγ γ f a |B| happens when the numerical value of the magnetic field at the surface is B ∼ 10 11 T. The conclusion given above is avoiding some subtle point. From the relation E = −g aγ γ a B it follows that the magnetic field and the electric field of the bubble are collinear. In this situation, it is known that the magnetic field B suppresses pair creation for scalar fields if it is large enough, but enhance it for fermions [100,102]. Thus, the pair production corresponding to the value B found above may be too large. To see that this is not the case, recall that the rate of electron positron pair creation when the collinear magnetic field is turned on is given by [102] From this formula it is seen that a large magnetic field increases w. In the present case, the magnetic field is large, but it is related to the electric field by E = −g aγ γ a B with a taking values close to 2π . By taking into account that g aγ γ = cα/π f a it follows that E ∼ α B, with α the fine structure constant. The rate given above is then 1 n e − m 2 e π n eE coth(137π n). (2.12) The leading term of (2.10) and (2.12) represent the pair production rate by unit time and unit volume. These are Now, a rough approach for estimating the value of the field E for which pair creation is significant may be to evaluate the first rate (2.13) at the value E c = m 2 e /e and find the electric field E such that the second rate has the same numerical value. The cotangent factor in (2.12) is close to the unity and can be neglected. Then the numerical relation defining E is thus given by The last equation can be cast in the Lambert form From here, the solution is given in terms in the Lambert function W (x) as follows This means that E ∼ E c /4.5. Thus, the magnetic field does not change considerably the value of the electrical field for which pair creation is appreciable. For this reason, the presence of the magnetic field when studying pair creation can be neglected in this specific case, even though the nugget magnetization is quite large. The estimation made above assumes that the electric field is uniform in the volume V f ∼ 4π R 3 /3. But this estimation may be not accurate, as the field E is varying inside the wall. Thus, it is of interest to study the effect of the inhomogeneity of the field inside the wall. A possible approach is to assume that the radius of the bubble is large enough and to analyze what happens at the pole. There are several works that derive results when such inhomogeneities are present [20,25,[30][31][32][33]41,[43][44][45][52][53][54]77,78], and in the following these references will be followed closely. In order to study the role of the inhomogeneities, the electric can be approximated by a one dimensional one E = E(z) where the axisẑ connects the pole to the center of the sphere. The field to be considered is of the Sauter form This choice is simply for convenience, as Schwinger pair creation is understood for this types of potentials [20,25,[30][31][32][33]41,[43][44][45][52][53][54]77,78]. We have investigated other type of shapes, as shown in the appendix. However, the results are mathematically more complicated and we do not believe that they change qualitatively the present result. One of the reason for believing that is that the Sauter electric field (2.14) is localized in a region of width R, and this imitates the field living on the axion wall bulk. Of course the functional form postulated above correspond to a flat situation and not to a spherical one. But for the sough estimation it may be enough. An important point is that, if the instanton method is to be used for estimating the pair production, the two conditions eE R >> m e and m e R >> 1 should be fulfilled. The second is immediately satisfied since l e ∼ m −1 e ∼ 10 −13 cm, which is much smaller that the nugget radius estimated above. The first condition can be rewritten as eE >> 1/Rl e . But for electric fields near to the critical value one has eE ∼ 1/l 2 e and therefore eE R ∼ R/l 2 e >> m e , as l e << R. Thus, both conditions are satisfied for fields with values close or larger than the Schwinger critical field. In this situation, there exists a formula for the pair production per area [77,78] which, applied to the present case, gives The region R of integration is given by In order to understand if the spatial inhomogeneity decreases pair production, assume that the electric field is such that δ = eE 0 /π m 2 e >> 1. This is equivalent to say that the pair production becomes significant when the value of E 0 is considerably larger the critical field (2.11). Then the last formula can be approximated by In these terms, one may calculate a value E 0 such that the last expression is equal to the value of N e in Eq. (2.13) evaluated at E c . This calculation yields the following result But this equation should be applied only for δ >> 1, as follows from the assumption that E 0 >> π E c . This condition contradicts the last equation. This implies that δ ∼ 1 and therefore the sough field can not be much larger than E c . This reasoning suggest that the spatial inhomogeneity does not alter considerably mean value of pair production on the axion wall. The result given are somehow similar to one related to scalar pair production for Sauter fields and found in [31], but generalized to fermions. In that reference, it is found that the quotient of the rate of pair creation for a Sauter field w s and for a uniform electric field w e is given by This quotient shows that the rate is suppressed due to the inhomogeneity parameter γ , which tends to zero when R → ∞. But the inhomogeneity parameter for the quark nugget with size R considered here is much less than unity, and there is no significant numerical difference between w s and w e . The results of the previous paragraph generalize partially this result to the case in which there is electron positron pair production instead of scalar fields. The above discussion suggest hat a magnetic field of B ∼ 10 11 T may induce pair production considerably. The value of the reference [103] is B ∼ 10 12±1 T, which shows that these magnetic field is at the very verge for stability. The magnetic field of the references [76,98,99] are also of this order. Thus, there may be considerably pair production, unless the combined physics axion quark gluon plasma does considerably change the value of the magnetic field.
It is not clear for us what is the real fate of such nugget. These apparent instability may indicate that the nugget evaporates, or simply that the ferromagnetic state is not effective due to pair decay of the electric, and thus the magnetic field disappears. In the following, the second possibility will be considered. In this case, the nugget will evolve to the state described in [35][36][37][38]58,[60][61][62][63]79,90,105,[109][110][111]

Further comments about the approximations made
There are some approximations made during the previous analysis, such as (2.9), which were justified by use of the classical Maxwell equations (2.4)-(2.8). However, due to the strong nature of the electromagnetic field considered, there are quantum corrections to be added to Maxwell equations. This is due to the fact that the influence on the quantum vacuum of these fields may generate non linearities, which should be properly included in the scenario. These linearities are an active research area, and have several interesting effects. First, the photon polarization tensor in presence of homogeneous electromagnetic external fields was studied in [48,50]. The effect of slowly varying inhomogeneities was also considered for instance in [49], and the magnetic response in presence of Coloumb sources was studied in [3]. In addition, the presence of magnetic fields produce an enhancement at the light shinning through wall experiments designed in [28,29], which may be seen in a near future. Further applications of these non linearities for laser fields were considered in [40]. In addition, the interaction between micro-bubbles with ultra intense laser pulses, by taking into account the QED vacuum polarization, was studied in [56].
A research of QED vacuum polarization in presence of axion fields was initiated in [34]. These result may have interesting cosmological applications.
The fact that the physics of QED non linearities is so rich raises an issue about the affirmation made above that the nugget electric field E n > E c , as it may be the case that the non linearities change its numerical value and render it subcritical E n < E c . We will argue below that this is not the case. In order to see these effects, recall that the exact one loop QED effective lagrangian is given by This is the well known Euler-Heisenberg lagrangian. Here L m is the Maxwell lagrangian, L aγ is the axion photon coupling defined above and the following quantities have been introduced As is well known, the Euler-Heisenberg lagrangian contains a real part and an imaginary one. The imaginary part describes pair creation, which was analyzed in previous sections. However, the real part describe corrections to the Maxwell equations due to quantum effects. It is difficult to evaluate the integral (2.15) explicitly. However, in some limits, its approximate form is known. For instance, the fields satisfying the condition << E c and β << E c are known as weak fields. In this limit, the real part of the lagrangian (2.15) may be approximated by (see [55,92,93] and references therein) The modified Maxwell equations that follow from (2.16) are Here τ n = nπ e , τ m = mπ eβ , and the following function has been introduced, with E(z) is the exponential integral function. Now, if it were the case that B >> B c and E << E c , then 2S ∼ −|B| 2 and P << −S. Thus The denominator in (2.19) grows with n and m. The functions J (z) in the numerator tend to zero for such large values.
Based on this, one may consider the terms with n, m = 0, 1 as the leading ones. This leads to the following approximated lagrangian 1 Still, if << β, this may be approximated by By taking into account that ∼ E and β ∼ B, and that This lagrangian should be added to the standard Maxwell one, together with the photon axion coupling. In similar fashion as above, we have played with the resulting equations of motion. We find that, if it is assumed that E << E c , the corrections to the Maxwell equations are small. On the other hand, if the corrections are small, this implies that the classical picture is still valid and E >> E c . This contradiction suggest that the assumption of small electric field is not quite right. The arguments given above are not complete proofs, and they do not have into account higher loop calculations [39,47]. But we conjecture that the effect of the non linearities will not correct the electric field from a classical value E >> E c to a corrected value E << E c . It is worthy to mention that there exists a discussion in the literature about the possible role of the non linearities in QED for rotationally powered pulsars, with large magnetic fieldsr, as shown in [91] and references therein. The reference [91] argues that these effect are negligible, at least for certain types of neutron stars. Based on these heuristic argumetns, we assume that the non linearities of QED do not change the supercritical value of E n and that the relation (2.9) and the estimated values obtained in previous section still apply for axion quark nuggets.

Enhanced cross section at the troposphere
The study of the cross section of a quark nugget passing through the troposphere is important for understanding which methods are suitable for their detection.
It should be recalled that for a non magnetized quark nugget [23] it is usually assumed that its cross section σ n is given by the cross sectional area of its core mass density. In this terms the following formula for the energy loss of a quark nugget d E ds = −σ n ρv 2 , σ n = π 3m 4πρ n 2 3 , is found. One may consider velocities of the order v ∼ 250 km/s, which are characteristics of the Sun galaxy rotation.
Here ρ is the medium density the nugget is passing through. The formula given above is inspired by the physics of meteorites. The quark nugget density is approximately ρ n ∼ 10 18 kg/m 3 . This formula allows to analyze if the nuclearity passing through the earth will accumulate at the crust or if it will pass through the earth. Besides, a nuclearity that pass through the air or water deposit part of its energy in terms of visible light. The authors of [23] calculate the fraction of the dissipated energy in terms of visible light, by assuming that the light is emitted as a black body radiation from an expanding cylindrical shock wave. In this terms, they derive a formula for the luminosity of these objects and discuss possible detection mechanisms based on these formulas. This picture changes when the nugget entering into the earth has a magnetic field B turned on. A first attempt for studying their energy deposition is to consider some known formulas for magnetic objects such as monopoles, which in principle may be also applied to a magnetized nugget [18]. However, it may be convenient to generalize the results obtained in [81] since they take into account that the nugget is surrounded by a plasma when it enters in the troposphere. In this approach, the physics of a magnetized nugget is as follows. When the nugget enters into the 20 km region covered by the troposphere, it finds a neutral medium. There is oxygen and nitrogen in this region. The oxygen binding energy is around 8 KeV. The Zeeman term that arises from the interaction of these atoms and the nugget magnetic field is given by It is seen by passing to natural units that L z ∼ 6.10 −5 eV T −1 B. Therefore the Zeeman energy will have a numerical value close to the keV scale when B ∼ 10 7±1 T. By taking into account that the magnetic field depends on the radius by a law of the form B ∼ r −3 in the model introduced in [103] and by taking into account that at the surface B ∼ 10 11±1 T, one may assume that, inside a radius r ∼ (10 − 10 2 )r n , there is enough Zeeman energy for ionizing oxygen and hydrogen. Thus, in this region, the surrounding medium is approximately an ionized plasma. Assume that the nugget has a velocity v, which is not relativistic. In the system in which the bubble is at rest there is an incoming plasma with a speed −v. The particles tend to make a round trajectory, of Larmor type, when approaching the region where the magnetic field is strong, the plasma pressure does not allows this to fully happen. Thus, there is a surface formed in which the plasma pressure equals the magnetic pressure. The plasma ram pressure is P r ∼ ρv 2 up to a model dependent constant which is not far from unity. The magnetic pressure is P m ∼ B 2 /μ 0 . By taking into account the dependence B = B 0 r 3 n /r 3 is considered, then the condition P m = P r shows that the plasma particles are allowed to approach the nugget up to a radius The cross section then can be approximated in natural units by [81] Q n ∼ πr 2 e ∼ In other words, this cross section is simply the two dimensional area of a disc of radius r e , which is enhanced due to the magnetic field. At this point, the calculations given above are the ones already presented in [81]. This effect was employed in that reference to interpret the results of the experiments described in [85], which did not take into account the effect of the possible magnetic field of the nugget. These calculations assume that the nugget electric field is zero. However, as discussed in previous sections, in presence of an axion domain wall surrounding the object plus a magnetic field, an induced electric field E n is induced. Thus, one needs to evaluate the effect of the electric field in those formulas. Our estimation is that the formula (3.2) is not significantly changed, and the reasoning goes as follows. The plasma may be considered a fairly good conductor, and it will try to expel somehow this electric field by forming a charged surface, a shield, which distorts the value of r e . There are two electric field contributions, the nugget original field E n and one the induced plasma E p . This electric field repels charges of one sign and attracts the others. On the other hand there is a tendency of the plasma to screen the electric field, as the plasma is a good conductor. It is not clear if the presence of the electric field deforms considerably the magnetic surface r = r e found above or not. However, some estimation can be made. Consider the new surface r = r e , the border to which the charges are allowed to approach the electrified and magnetized nugget. If this surface is though as the boundary of a conductor, the force on a surface element with area S is F = σ E c S where E c is the surface electric field. The electric field of an ideal conductor is normal to the surface and E c ∼ σ/ 0 with σ the surface charge density. Thus, the resulting force can be expressed as F = 0 E 2 c S, from where it follows the pressure P = 0 E 2 . By extrapolating this picture to the present situation, one may consider a pressure given by P = 0 E 2 with E = E n + E p , with both the induced plasma electric field E p . It may be fairly reasonable to assume that both fields are of the same order on the sought surface r = r e . As shown in previous subsections, the electric nugget field is given by By taking into account that g aγ γ = cα/ f p π and that a = f a θ with θ an angular variable, it follows that The constant c is not far from unity and the angular variable θ can take values between 0 ≤ θ < 2π N with N an integer, it is seen that the electric pressure P e ∼ 0 E 2 is less or of the same order of magnitude than the magnetic one P m ∼ B 2 for controlled values of N . Thus, it is likely that the region determined by (3.1) is not significantly deformed, but the charge distribution is modified in order to screen the nugget electric field. Therefore, it may be considered that the formula (3.2) applies even for this type of electrified and magnetized nuggets. It should be remarked that the formula (3.2) is obtained purely by a classical considerations. It is true if the quantum corrections described in previous section do not alter the order of magnitude of r e . In the following, we assume that this is the case.
There are observables consequences of the cross section (3.2), which follows from reference [81]. This can be briefly described as follows. The formula (3.2) allows to calculate the energy loss of a nugget of mass M and a given velocity of the order of 250 km/s through the atmosphere and three meters of material with density ρ ∼ 1120 kg/m 3 . This simulation is presented in [81] and it is found that the mass range for nuggets is such that they may to penetrate the atmosphere and deposit 10 3 -10 6 J/length. For instance, this curve shows that nuggets with masses between 0.1-1 kg deposit 00.3-50 MJ/m. This energy deposition is enough to form a shock wave in water which decays fast into an acoustic pulse. The pressure of this pulse can be monitored by three or more time synchronized sensors, which may determine the impact point by interpolation. From this information, one may deduce the energy deposition and, by use of the curves found in [81], obtain the mass of the nugget. In order to decrease the risk of obtaining a false positive, it is desirable to use three sensors in coincidence.
The measurements requires to understand the energy deposition on water to pressure pulses. For simulations of this, one may use the CTH shock physics code [75] and the SESAME equation of state in water [22]. Some of these simulation are presented in [81]. A good location for detection is the Great Salt Lake in USA, since its salinity does not allow animal life which may contaminate acoustic pulses. There are further technical problems to be solved. The incident flux as a function of the mass is not known, as the quark nuggets are assumed to have a continuum mass distribution. The Auger [2] and Ice cube [42] experiments may through light about this distribution. But, under the assumption that the events by mass unity varies inversely with the mass, the authors of [81] used a 250 h sensor which collected events at a rate 300/year, thus obtaining a graphic of the quantity of events by mass unity. Their results suggest that detectability is possible. In addition, the conversion of the deposited energy into acoustic energy and the attenuation of the of the acoustic pressure at the detectors as a function of the distance point should also be characterized, for the Great Salt Lake environment. The authors of [81] have studied the effects of the minerals of the lake in the attenuation of higher frequencies of the acoustic pulse and have found that the conversion into acoustic energy has some similarities with TNT explosions. These preliminary experiments were performed by C55 hydrophones, and further details can be found in the original reference.
The result given above is interesting, since it opens the possibility that these axion nuggets may be found in near future experiments. However, as discussed below, these electrified and magnetized axion quark nuggets may generate some high energy rays that are not observed due to their interaction with electron gases in the Milky Way. In this sense, the fact that the electric field is larger than the critical one E c may give an explanation of this absence, since this field may have decayed before the actual universe.

High energy ray generation
If there were electrified nuggets composing a large part of dark matter, then these nuggets would interact with electron gases in the Milky Way galaxy. The density of these electron gases is n e ∼ 1 cm 3 . These electrons, after colliding with the nuggets, will accelerate and acquire a energy eEm −1 a . A rough estimation about the number of collision per unit time that a nugget experiences is Due to the oscillations of the axion field, the electrons acquire an energy E e ∼ eE/m a . These acceleration mechanisms were discussed in other context for instance in [46]. This implies that when a astronomical observation of 1 kpc 3 is done, the number of nuggets would be where m n denotes the mass of the nugget and ρ d denotes the local density of the dark matter ρ d 0.3GeV cm −3 . The electrons are scattered by the collisions with the nuggets. The number N of the electrons scattered per unit time in the area is given by These electrons possess the very high energy eE/m a when the electric field is the strong one discussed in the present paper. For instance, when R is 10 −9 cm and m a = 10 −3 eV, then, the electron energy eE/m a ∼ 10 4 GeV and the number of the electrons N is 10 20 /s for baryon number N B = 10 28 . This is for a velocity v e ∼ 100 km/s. If such a huge number of the electrons with high energy 10 4 GeV were emitted in the area with the volume 1 kpc 3 , they would be already detected. Thus, it is likely that these nuggets may exist today only in a non magnetized state or, if they are magnetic and electric, their density should be considerably small than dark matter.
Therefore, if the nuggets reach a non electrified neither magnetized state, the detection of water impacts discussed in the previous subsection become more involved, as the cross section becomes considerably smaller. However, there are for now several searches of standard quark nuggets and it may be interesting to see if they apply to the present case. An interesting possibility is the detection at the Alpha Magnetic Spectrometer [4,70,71,94], which is sensitive to a wide range of mass and charge for these objects. In fact, if these objects reach the non magnetized neither electrified state, they may be similar to strangelets, which have been intensively searched at this spectrometer. However, there is an important detail to be remarked. One assumption for detection for strangelets is that neutron stars are strange stars. If this is so, when two neutron stars form a binare system which loss energy due to gravitational radiation, they finally collide, and this results in a strangelet cosmic ray. The galactic coalescence of our galaxy is approximately one even every 10,000 years. Simulations estimate that the release of 10 −10 solar masses per year. On the other hand, it is not clear if a considerable axion wall quark nugget flux may be obtained from such binare collisions, as it is not clear if such axion walls would develop. If the axion quark nuggets are assumed to be the main component of dark matter, then the flux is given by [60] This flux is below the sensitivity of the standard searches for dark matter, but is close to the flux of cosmic rays above the GKZ bound. Thus, a further possibility is to consider air showers [59] that these objects may generate. The shower that is generated by these objects have some similarities with the ones produced by a single primary with ultra high energy. But it is precisely their subtle differences which may play a significant role in detection. When an anti-nugget hit the atmosphere, the dominant interaction is strong mediated matter anti-matter annihilation. The resulting shower will be mainly composed by light mesons and their decay products. There are electrons and positrons produced, which are unable to escape the nugget. However, the muons are able to escape, loose energy to the atmosphere, generate fluorescence light before decay into electron products. This leads of a number of charged particles as a function of the height, which grown to a maximum particle and then decreasing rapidly, and is very similar to the one produced by a single primary. However, due to the fact that the nugget is not relativistic in comparison, the nugget shower fluorescence will elapse more time. The reason is that the matter produced by the nugget is confined to a region few kilometers around the object, and the combined system moves slowly. The time scales involving the lateral surface profile are also larger. A further characteristic feature is that the lower velocity of the primary particles will lead to a correlation between the arrival direction and the earths motion direction with respect to the galaxy. This arrival direction is related to dark matter distribution and as such, does not have any correlation with the galaxy objects. In addition, the muon spectroscopy may distinguish both showers since the muon spectrum for the nuggets have a cutoff at energies of the order of the GeV while for ultra high energy primaries this is not the case. Further characterization of these showers can be found in [59].
The differences between ultra high energy may play a very important role in future detection experiments. These include the AUGER observatory [2] and the Telescope Array [104] which may detect these showers by their fluorescence detectors. The reference [88] is able to put flux limits that may apply to strangelets, axion wall quark nuggets and other objects. The macro detector at Gran Sasso may also be sensitive to these objects [9]. It is interesting to remark that the results presented in [106] suggest that the quark nuggets are not excluded when their mass is somewhat larger than m > 10 22 kg. The study of the constraints that follows from ancient mica exposure, of ages of the order of 10 9 years, is also an interesting lead [87]. Another possibility is that these objects may be detected by impact craters [89]. This craters should be compatible with energetic events without the presence of meteorites. The group of authors [81] is searching for such events in County Donegal in Ireland, and have analyzed several reported events in India or Nicaragua, but where unable to find unambiguous results. However, further experiments may be performed and the crater may be distinguished from usual meteorites by the absence of CUDO and meteorite materials.

Discussion
In the present work a ferromagnetic axionic quark nugget was considered. It was shown that the internal magnetic field induces an electric field on the axionic wall which decays into electron positron pairs if the magnetic field is above the critical value B c = 10 11 T, which close but below the range given in [103]. The possible role of the inhomogeneities of the electric field was also taken into account, but it was shown that the deviation from the standard Schwinger formula is small and the critical value remains also unaltered. The fact that the critical value is below the value of [103] may be in fact healthy. It may indicate not that the nugget is unstable, but that in fact the ferromagnetic state does not take place. The calculation done contains several approximations, and to study the real evolution of the nugget is for sure more complicated. But we conjecture that the ferromagnetic state is not achieved when an axion wall is present. Consequently, the interaction of these nuggets with the electron gas in the Milky Way does not give a huge number of electrons with high energy contradicting the observations. In other words, we suggest that ordinary nuggets may reach a ferromagnetic state, but the ones with axionic walls evolve to a non ferromagnetic state of the form [35][36][37][38]58,[60][61][62][63]79,90,105,[109][110][111], due to the influence of the axion wall during its the evolution.
On the other hand, the energy deposition on the earth due to these objects was studied. For the electrified and magnetized state, this was investigated by generalizing the magnetopause model considered in [81]. This model follows an analogy between the plasma surrounding the nugget and the solar wind interacting with the magnetic field of the earth. In the present work, the effect of the electric field was included and it was shown, by use of the modified Maxwell equations in presence of an axion, that the electric field pressure on the plasma is less than the magnetic one. Consequently, the cross section for a nugget interacting with a plasma is not considerably modified. A direct application of the results given in [81] implies that these nuggets, even in the electrified and magnetized state, do not violate the collision bounds for dark matter. However, as we conjectured that the electrified and magnetized state do not take place, alternatives for non electrified neither magnetized nuggets were also briefly discussed. Perhaps the study of air showers is promising, since the products of the decay can be qualitatively different from a single primary with ultra high energy.
Acknowledgements O. S is supported by CONICET, Argentina. We gratefully acknowledge a discussion with A. Ilderton for providing us some details about his work on pair production, and to F. Karbstein for orientation about effects of non linearities in QED.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical study and no experimental data has been listed. Experimental sources has been cited in the references.] 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: Electric field from a spherical domain wall with constant magnetic field
In the present paper, the role of inhomogeneities were estimated in terms of the Sauter field (2.14). We have worked out some other possible forms of the electric field, but the results are more complicated and do not change much the presented picture. We collect however in the present appendix some of the formulas that we have obtained which do not correspond to the Sauter field.
The axion domain wall profile the reference [64] employs is given by (1.23) The formula (1.22) corresponds to the regime R trans > r > R 0 and (1.23) to the regime r > R trans , where the constants γ = 1/2m a y A = π/2 f a have been introduced. The value R trans is given in [64]. For simplicity, assume that the magnetic field inside a nugget has a constant value and direction. In spherical coordinates, such magnetic field is given by Given a magnetic field and an axion wall, an electric field E is induced. The Poisson equation for the potential V such that ∇V = −E that arises from the term ∇(a B) in (2.4) is expressed in the form (1.26) the first expression corresponds to R trans > r > R 0 and the second r > R trans . We can solve both equations using the Green function method, but due to the fact that the effective charge is complicated multipole expansion is required. Due to the azimuthal symmetry, the Green function becomes where r > and r < denotes the biggest and smallest respectively between r and r . One of the boundary conditions is that the potential goes to zero at infinity. The monopole therm is zero for both equations. In this section, the dipole case is worked out explicitly, the correction terms for the general case are given in the next one. In the region R 0 < r < R trans the solution is of the form where Q 1 , Q 2 y Q 3 denotes integrals that come from the moments of the field given by In the region R trans < r < R the resulting potential is In this case the integrals Q 4 , Q 5 y Q 6 are described by (1.35) The integrals described in each solution can't be solved by analytic means in general. Therefore, some reasonable approximations are needed for solving them. For the first integral (1.29), one may use the expansion and that, for γ (R trans − R 0 ) << 1 By use of this, it is obtained that For the second equation we are gonna use the Laurent series of the hyperbolic tangent and truncate the series to get and, like the first equation, we are gonna use sech 2 (x − b) ∼ sech 2 b, both approximations apply in the regime a(R trans − R 0 ) << 1. With these approximations, it is found (1.37) For the third integral, one may take into account that tanh(x − b)/x 1/x, for a sufficiently large nugget radius, γ R > 1. A direct integration leads to (1.38) The next integral, can be solved using the same tricks we used for the first one, the result is The fifth one may be found by the approximation tanh(x − b) ∼ 1 − 2e −2(x−b) , and by the expansion sech 2 (x − b) = 4e 2(b−x) (1 − e 2(b−x) ). Both approximations lead to the expression − π 2 e 2(γ R 0 +δ−x) 6(γ x) 2 +6γ x+3 +2x 3 r Rtrans + e 2γ (R 0 +δ−2x) 24γ 3 × e 2γ (R 0 +δ) 32(γ x) 3 + 24(γ x) 2 + 12(γ x) + 3 +8e 2γ x 4(γ x) 3 + 6(γ x) 2 + 6(γ x) + 3 r

Rtrans
(1.40) The last one can be integrated using the same trick as the ones used in the third one. The result of this procedure yields (1.41)

A.1 Higher order corrections
Having the functional form of the potential for the dipolar case is not the end of the story, due to the fact that there are gonna be higher order correction due to the fact that cot θ can be expanded in an infinite series of Legendre functions. The solution of the potential in the regime r > R trans at all orders is given by (1.42) This expansion includes the terms for the dipolar case found in the previous section, plus higher orders. First, let notice that the dipolar order is important due to the fact that cos θ is proportional to Y 0 1 and will be orthogonal to the rest of the spherical harmonics. Second, let's also notice that cot θ goes to ∞ at 0 and to −∞ at π , thus infinite terms in the expansion are needed.
The higher order corrections, when the approximation tanh(x − b) = 1 − 2e −2(x−b) was taken into account, are given by Here (a, b) denotes the incomplete gamma function. The corrections for the other region R 0 < r < R trans are similar to the ones here and we will not discuss here the functional form here.

A.2 Schwinger effect for the non ferromagnetic axion quark nugget
After finding approximate expression for the electric field induced due to the axion nugget, the next step is to discuss the decay due to Schwinger effect. The following discussion relies particularly on the work [54], and it will be applied to a dipole type of electric field outside the nugget and a constant field inside. For simplicity, scalar field production will be considered. The Klein-Gordon equation for a pure electric field, by taking into account that A 0 = V and A i = 0 reads as follows In the following, a reduced one dimensional problem in one direction will be considered. This reduction is valid if the three dimensional electric fields E(x, y, z) and the potential V (x, y, z) are such that the tunneling length a = z + − z − to be defined below is smaller than the variation δx ⊥ of the potential V the x y plane. In other words, Taking only the radial tunneling has two implications. First the angular part of the potential is discarded and only a radial part V (r ) is taken into account. In this situation, the ansatz employed in the literature [54] is the following one = e i t ψ l (r ) r (θ)e imϕ . (1.45) In the previous section, some approximated expressions for the electric potential of the nugget were found. It is difficult to obtain results with these expressions. For simplicity the potential to be employed is a dipole type one which, in one dimension, reads as follows Here a = p b ot 2 c 2 +m 2 e c 4 − 2 (qv 0 R 3 ) 2 and b = 2 qv 0 R 3 . In these terms the event density in four space is given by with D s the spin degeneration. The formula (1.58) takes into account only an electric field. But these methods can be generalized when a magnetic field is present, if one assumes for simplicity that the magnetic field is constant [54]. The magnetic field in a magnetic nugget is gonna be parallel to the electric field. The resulting event density is given by e −π(E c /E 0 )G(0,E(z)) . (1.59) Here the function sinh(π x) σ = 1/2 (1.60) has been introduced. The limit B → 0 the formula (1.58) is recovered.