Born-Infeld magnetars: larger than classical toroidal magnetic fields and implications for gravitational-wave astronomy

# Born-Infeld magnetars: larger than classical toroidal magnetic fields and implications for gravitational-wave astronomy

Received: date / Revised version: date
###### Abstract

Magnetars are neutron stars presenting bursts and outbursts of X- and soft-gamma rays that can be understood with the presence of very large magnetic fields. Thus, nonlinear electrodynamics should be taken into account for a more accurate description of such compact systems. We study that in the context of ideal magnetohydrodynamics and make a realization of our analysis to the case of the well known Born-Infeld (BI) electromagnetism in order to come up with some of its astrophysical consequences. We focus here on toroidal magnetic fields as motivated by already known magnetars with low dipolar magnetic fields and their expected relevance in highly magnetized stars. We show that BI electrodynamics leads to larger toroidal magnetic fields when compared to Maxwell’s electrodynamics. Hence, one should expect higher production of gravitational waves (GWs) and even more energetic giant flares from nonlinear stars. Given current constraints on BI’s scale field, giant flare energetics and magnetic fields in magnetars, we also find that the maximum magnitude of magnetar ellipticities should be . Besides, BI electrodynamics may lead to a maximum increase of order of the GW energy radiated from a magnetar when compared to Maxwell’s, while much larger percentages may arise for other physically motivated scenarios. Thus, nonlinear theories of the electromagnetism might also be probed in the near future with the improvement of GW detectors.

journal: Eur. Phys. J. C

\thankstext

e1e-mail: jonas.pereira@ufabc.edu.br \thankstexte2e-mail: jaziel.coelho@inpe.br \thankstexte3e-mail: rafael.lima@udesc.br

## 1 Introduction

Soft Gamma Repeaters (SGRs) and Anomalous X-Ray pulsars (AXPs) (also known as magnetars (see, e.g., 2014ApJS..212….6O, , and references therein)) are spectacular astrophysical systems which hold the unique possibility of probing yet unknown particle physics under extremely high magnetic fields. Believed to be “lonely wolves”, among other properties such objects are observationally characterized by outbursts of X-ray and soft gamma-ray flares and have rotational periods  s and slowing down rates  s/s (2017ARA&A..55..261K, ). Duncan and Thompson (1992ApJ…392L…9D, ; 1995MNRAS.275..255T, ) proposed that they would be neutron stars (NSs) with huge magnetic fields (magnetars), of the order of G. One of the reasons for this would be their transient activities in the form of giant flares, whose typical luminosities are - erg s (2017ARA&A..55..261K, ), not possible in general to be powered by their rotational energy. 111Recently, though, the authors of Ref. 2017A&A…599A..87C () have discussed the possibility of some SGRs/AXPs being rotation-powered NSs and explored the entire range of NS parameters allowed by the conditions of stability of the star, not only fiducial parameters. Notwithstanding, three counterexamples of magnetars with low surface magnetic fields are already known (2010Sci…330..944R, ; 2012ApJ…754…27R, ; 2014ApJ…781L..17R, ). This apparent “glitch” of the magnetar model has motivated different scenarios for the explanation of SGRs/AXPs, e.g.: the possibility of a fallback disk slowing down a neutron star pulsar up to its current spin period  (see, e.g., 2011ApJ…732L…4A, ; 2013ApJ…764…49T, ); drift waves near the light-cylinder of NSs (see 2010ARep…54..925M, , and references therein); exotic scenarios involving quark stars (2006MNRAS.373L..85X, ); and massive, fast rotating, highly magnetized white dwarfs to explain these types of sources (2012PASJ…64…56M, ; 2014PASJ…66…14C, ; 2017MNRAS.465.4434C, ). Low-B values are related to surface poloidal magnetic fields, obtained by assuming that the spinning down of a star is due to its magnetic dipolar radiation (1986bhwd.book…..S, ), namely (2017ARA&A..55..261K, ). Thus, low-B observations rendered the issue of magnetars much more complex than just outbursts and high dipolar magnetic fields and strengthened the relevance of toroidal magnetic fields in stars.

Numerical simulations of ordinary twisted magnetic fields in magnetars suggest they should be poloidal field dominated, which is at odds with some inferences from observations (see details in 2010MNRAS.406.2540C, ; 2013MNRAS.435L..43C, , and references therein), such as the three above mentioned low-poloidal-B magnetars. This has been addressed in the context of a modified twisted magnetic field, where ratios of the toroidal energy to the total magnetic energy in stars could be up to (2013MNRAS.435L..43C, ). This has been done in the context of Maxwellian electromagnetism and toroidal fields as high as G could be obtained. However, for such large magnetic fields, nonlinearities of the electromagnetism should play a significant role in the observational properties of the systems bearing them. Motivations for nonlinear electrodynamics in the context of high fields naturally come from QED (2010PhR…487….1R, ) but in this paper we focus on an older approach, due to Born and Infeld (1934RSPSA.144..425B, ), interesting even nowadays due to the reasons that follow.

Born-Infeld’s electromagnetism was conceived with the purpose of healing the energy singularities of point-like charged particles such as the electron and is motivated by the finiteness of physical observables in special relativity (1934RSPSA.144..425B, ). It was showed in the 1980s that it is consistent with the low-energy limit of string theory and this has, since then, drawn very much attention to it (see 1997hep.th….2087R, , and references therein). It already has an exact black hole solution (see for instance 2015PhRvD..91f4048P, , and references therein) though in this context singularities are unavoidable (2001PhRvD..63d4005B, ). When applied to the hydrogen atom, Born-Infeld theory as the description for electromagnetic interactions meets the observational spectrum of this atom only if its scale factor is larger than the one inferred by Born and Infeld themselves within the unitarian viewpoint (2006PhRvL..96c0402C, ; 2011PhLA..375.1391F, ; 2018EPJC…78..143A, ), approximately statvolt cm (or also  G) (1934RSPSA.144..425B, ) (around 100 times QED’s critical fields 2010PhR…487….1R ()). There are also constraints on coming from particle accelerators. It has been recently shown that LHC light-by-light scattering in Pb-Pb collisions would restrict to be larger than statvolt cm G [ G statvolt cm] (2017PhRvL.118z1802E, ). However, it is important to point out that the kinematic cuts made in ATLAS and the self-consistency of the linear analysis used for cross section calculations are such that this experiment does not allow for the test of smaller values of (2017PhRvL.118z1802E, ). All the above means that the window G  G cannot be assessed by LHC experiments so far and hence they leave a “hole” in the probe of the parameter. Given that magnetars are believed to have surface fields as high as G (even larger fields in their interiors), they could be the ideal systems for assessing scale fields to the Born-Infeld theory in the region the LHC cannot. This is important and is a motivation our analysis because the Born-Infeld electrodynamics can only be disregarded when its whole space of parameters is excluded. Besides, magnetars could also be ideal testbeds for other nonlinear theories of the electromagnetism. Impressively enough, perhaps due to the intrinsic difficulties already present in the Maxwell theory, such analysis seems scarce in the literature of magnetars. We try to partially fill this gap in here.

We structure this work in the following way. Next section is devoted to the deduction of the field equations in nonlinear electromagnetism for the case of ideal magnetohydrodynamics. In Sec. 3.1 we work out a realization of nonlinear toroidal fields, by focusing on some consequences of it for the case of Born-Infeld electrodynamics. Some astrophysical consequences thereof are investigated in Sec. 4, especially related to the increase of the magnetic ellipticity and gravitational-wave energy budget when compared to their Maxwellian counterparts. Finally, discussions and conclusions are given in Sec. 5. We work with Gaussian and geometric units.

## 2 Nonlinear electrodynamics in ideal magnetohydrodynamics

For a Lagrangian density depending upon both invariants of the electromagnetism, and , the dual of (1975ctf..book…..L, ), the appropriate electromagnetic action is

 Sem=a∫d4x√−gL(F,G)−∫d4x√−gjμAμ, (1)

where is a numerical factor that depends upon the system of units used, is the determinant of the spacetime metric (assumed in this work to be given), is the current four-vector of the system and is the four-potential of the electromagnetic fields. We work with Gaussian units, which means we take . By assuming that and varying Eq. (1) with respect to , we obtain the field equations

 ∂μ(√−gLFFμν+√−gLG\raisebox0.86pt$\lx@stackrel∗F$μν) = 4π√−gjνor (LFFμν+LG\raisebox0.86pt$\lx@stackrel∗F$μν);μ = 4πjν, (2)

complemented with (a natural consequence of the definition of ) (1975ctf..book…..L, )

 ∂μ(√−g\raisebox0.86pt$\lx@stackrel∗F$μν)=0or(\raisebox0.86pt$\lx@stackrel∗F$μν);μ=0 (3)

where , , , , is a totally antisymmetric tensor. Note that since the Lagrangian must be an even function of the invariant (due to symmetry requirements from electrodynamics), is an odd function of , which means it is zero when the fields are orthogonal.

In order to simplify our description and evidence the physical picture involved in our analysis, we assume now the case where the electromagnetic fields are orthogonal and take the spacetime metric to be the Minkowski metric. Crossed fields naturally arise in the context of ideal magnetohydrodynamics (MHD) [see e.g., (1986bhwd.book…..S, ; 2009lema.book…..S, )], which we will consider throughout this work, and the use of a flat spacetime metric, though being an idealized case, will allow us to find the maximal changes of physical quantities regarding Maxwell’s electrodynamics. 222Indeed, when general relativistic corrections are taken into account, quantities such as the magnetic fields should decrease with respect to their flat spacetime counterparts by a function depending on the compactness factor (2004MNRAS.352.1161R, ; 2015ApJ…799…23B, ; 2017A&A…599A..87C, ). We plan to investigate the case of nonlinear electrodynamics in curved geometries as due to magnetars in a forthcoming work. As commented previously, for the case under interest and Eqs. (2) and (3) can be cast as (we restore Gaussian units here)

 ∇⋅(LF→E)=4πρ, (4)
 ∇⋅→B=0, (5)
 ∇×→E=−1c∂→B∂t (6)

and

 ∇×(LF→B)=1c∂LF→E∂t+4πc→j, (7)

where we have defined as the system’s charge density and its current vector.

Since and are given aspects of a system, one could take as reference Maxwell’s electromagnetism. There, and , which allows us to recast Eqs. (4) and (7) as

 ∇⋅(LF→E−→EMa)=0 (8)

and

 ∇×(LF→B−→BMa)=1c∂∂t(∂LF→E−→EMa). (9)

Care must be taken at this point due to Eqs. (5) and (6), which are formally identical to their Maxwellian counterparts but are related to the fields and , instead of and . This means, for instance, that a dipolar magnetic field in Maxwell’s electromagnetism cannot in general have the same functional form in nonlinear electrodynamics.

Nonetheless, for a very good conductive region of a rigidly rotating star with an angular frequency in the regime of small velocities () , one has that (ideal MHD (1986bhwd.book…..S, ; 2009lema.book…..S, ))

 →E=−→vc×→B=−→ω×→rc×→B, (10)

which comes from the assumption that Ohm’s law also holds for nonlinear electrodynamics, taken in this work as a first approach, and the consideration the current density is finite in the limit of infinite conductivity. 333Let us elaborate on this point. If one assumes the validity of Ohm’s law also in nonlinear electrodynamics, then in a locally comoving reference frame (or rest-frame) (with respect to a conducting fluid), one has that , where and are the current density and the (nonlinear) electric field there, respectively, and is the conductivity of the fluid. Take now the limit (ideal conductor). If one assumes that is finite (could be any), then, necessarily, . Make use now of another inertial coordinate system (usually called the laboratory frame) such that moves with respect to it with velocity . When , it follows that (1975ctf..book…..L, ). (Due to the smallness of the charge density in a conductor after a small characteristic time (2001imhd.book…..D, ), from the transformation laws of four-vectors in the limit of , it also follows that .) Thus, in the limit of ideal MHD, Eq. (10) follows as a kinematic constraint, and hence would be valid for any theory of the electromagnetism.

If one assumes that the magnetic field lines are dragged along the motion of the star (frozen), then it follows that (Alfvén’s theorem) (2009lema.book…..S, )

 ∂→B∂t=∇×(→v×→B)=∇×(→ω×→r×→B). (11)

However, in nonlinear electrodynamics this point might be subtler in principle. If one inserts (from Ohm’s law in first order in ; see footnote) into Eq. (6), takes to be constant, assumes that is negligible when compared to the current density (which could be justified due to the fact that in a given limit the equations of nonlinear electrodynamics tend to the Maxwell ones where that holds true for large (2001imhd.book…..D, )), one has that

 ∂→B∂t=∇×(→v×→B)−c24πσ∇×∇×(LF→B). (12)

One sees from the above equation that Eq. (11) is recovered when (ideal MHD) and that Eq. (10) arises in this limit. The above means that the magnetic flux through any closed loop in a perfectly conducting medium is constant also in nonlinear electrodynamics if the conditions leading to Eq. (12) hold. We leave analyses of other scenarios (for instance finite , violation of Ohm’s law, etc.) to be investigated elsewhere.

From Eq. (8), one has

 LF→E=→EMa+∇×→C, (13)

where is an arbitrary vector. If one assumes that the time derivative of is negligible [which should be justified by the assumption of large conductivity and would be equivalent to the disregard of the time derivative term of Eq. (7) w.r.t. , as happens in Maxwell’s ideal magnetohydrodynamics (2009lema.book…..S, )], then it follows from Eq. (9) that

 LF→B=→BMa+∇f, (14)

where is also an arbitrary function. One can solve Eq. (14) for (given a nonlinear electrodynamics), which, after replaced in Eq. (5), will fix the function .

From Eqs. (13), (14) and (10), it immediately follows that

 →EMa+∇×→C=−→ω×→rc×(→BMa+∇f), (15)

which from the Maxwellian relationship of the fields [] implies

 ∇×→C=−→ω×→rc×∇f. (16)

Thus, any non-null induces a non-null and this is done through the assumption of the system’s large conductivity, exactly as suggested by previous intuitive arguments. We stress that Eqs. (13), (14) and (16) are simple just because they relate fields in different theories. Without the solutions for and , which are not easy in general, they are just indicative.

The electromagnetic energy density in nonlinear electrodynamics can be easily obtained with the mixed 00 component of the electromagnetic energy momentum tensor, which to the case of Eq. (1) for orthogonal fields is given by (2017PhRvD..95b5011M, )

 16πTνμ=−4LFFμαFγβgνγgαβ+Lδνμ, (17)

where is the Kronecker delta function (1975ctf..book…..L, ). When the hypothesis of very conductive media is taken into account (), Eq. (17) implies that

 16πT00≈L. (18)

## 3 Nonlinear toroidal fields

In order to have some insights into the influence of nonlinear electrodynamics on magnetars and motivated by the expected dominance of azimuthal fields inside such systems (2013MNRAS.435L..43C, ), here we focus on the nonlinear field equations for purely toroidal fields with axial symmetry (thought of as a rough approximation to the real scenario), such that , where is the azimuthal unit vector. This field is such that Eq. (5) is automatically satisfied. From Eq. (10), the case of toroidal fields in rigidly rotating stars (we define the -axis such that , which means , being any) result in , which automatically satisfies Eq. (6) [as well as Eq. (11)] and from Eqs. (13) and (16) implies . Thus, Eqs. (14) (the only remaining nonlinear equation) leads to

 LFBϕ=BMaϕ. (19)

Since in the small field regime of nonlinear electrodynamics (2016ApJ…823…97C, ), it follows generically from the above equation that the magnitudes of nonlinear toroidal fields in ideal MHD are in general larger than their Maxwellian counterparts. Given that, we investigate an interesting case with an exact solution in the next subsection.

### 3.1 Born-Infeld toroidal field analysis

Born-Infeld’s Lagrangian is defined as (see for instance 2014PhRvA..89d3822D, ; 2017PhRvD..95b5011M, )

 LB.I≐4b2(√1+F2b2−G216b4−1), (20)

where is the scale field of the theory and from the definition of the invariants of the electromagnetism, and . Thus, for orthogonal fields, our main interest in this work, the only relevant derivative to the Lagrangian is

 LF=b√b2+→B2−→E2≈b√b2+→B2, (21)

since for very good conducting stars it follows that . To the case of Born-Infeld toroidal magnetic fields, Eqs. (19) and (21) lead to a very simple solution, namely,

 Bϕ=bBMaϕ√b2−(BMaϕ)2, (22)

which clearly evidences the increase nonlinear electrodynamics impinges on .

For toroidal-dominated fields, it follows from Eqs. (18), (20) and (22) that the energy densities in Born-Infeld and Maxwell’s theories are related by

 (23)

which is always larger than the unit for , needed for the consistency of Born-Infeld toroidal fields. This means that toroidal fields in nonlinear electrodynamics imply larger energy reservoirs when compared to Maxwell’s predictions.

### 3.2 Ellipticities for Born-Infeld toroidal fields

Strong toroidal magnetic fields also have important implications for the ellipticity, as we discuss now. Assume that the magnetic ellipticity of a star is proportional to the mean value of , i.e.,

 ϵ=c1⟨B2ϕ⟩, (24)

where is a constant that depends upon aspects of the star such as its radius, compactness, etc. Indeed, when toroidal magnetic fields are dominant, one has that (1969ApJ…157.1395O, )

 ϵ∝−1Eg∫dVB2ϕ, (25)

where is the magnitude of the gravitational energy of the star and is a volume element. Thus, by defining , Eq. (24) ensues with depending on the compactness of the star and its radius, all byproducts of its microphysics.

From Eq. (22), it follows that the Born-Infeld to the Max-wellian ellipticity ratio is

 ϵBIϵMa=⎛⎝˙EBIGW˙EMaGW⎞⎠12≈(ΔEBIGWΔEMaGW)12≈⎡⎢⎣1−⎛⎝BMaϕb⎞⎠2⎤⎥⎦−1, (26)

where is the energy loss due to gravitational waves. We have assumed that is slowly varying within the star (where it takes place), which might be taken for maximal estimates. Since from consistency of Eq. (22), one has from the above equation that the production of gravitational waves is larger in the Born-Infeld theory than in Maxwell’s.

Assume now that is given, so as an upper limit to the norm of the magnetic ellipticity of a star and the Maxwellian toroidal magnetic field. (One could easily estimate by means of magnetar’s observables such as flare luminosities and upper limits to the ellipticities could be inferred from gravitational-wave analysis.) In this case, from Eq. (24), it follows that Maxwell’s ellipticity is known. If one takes the true theory of the electromagnetism as Born-Infeld’s, thus

 b≥BMaϕ(1−1C)−12,C≐|ϵul||ϵMa|, (27)

where stands for the observational upper limit of the magnetic ellipticity. From Eqs. (22) and (24) one sees that when the measured ellipticity is the one connected with the Born-Infeld theory, then it follows that automatically. Besides, Eq. (27) is totally equivalent to the consistency of Eq. (22) since the minimum value of is , which is larger than . Therefore, Born-Infeld electrodynamics is self-consistent.

## 4 Possible Astrophysical Implications of Born-Infeld magnetars

LIGO measurements already constrain ellipticities in ordinary pulsars and their magnitude should be smaller than approximately (2017ApJ…839…12A, ). The values of the ellipticities they found should be taken as upper limits since no continuous gravitational-wave signals have been detected from neutron stars. In what follows we do not assume the above value holds true for magnetars, but we rather find an upper limit to it by means of outcomes of Born-Infeld theory applied to the hydrogen atom. We will focus on the systems supposed to have dominant toroidal fields such that the analysis of the previous section could be taken into account. This is actually believed to be the case in all magnetars and even possibly in several pulsars (2013MNRAS.435L..43C, ).

In this context, mean toroidal Maxwellian fields could be estimated with giant flare events by dint of

 (BMaϕ)2=6EflR3, (28)

where is the radius of the magnetar and is the total energy released in the giant flare. Precise calculations for the case of a homogeneous ellipsoid show that (1969ApJ…157.1395O, )

 |c1|≈R4GM2, (29)

where is Newton’s constant and is the magnetar’s mass. From Eq. (27) [or Eq. (26)], could be constrained if a minimum value for , , was given, implying that

 |ϵul|=⎡⎢⎣1−⎛⎝BMaϕbmin⎞⎠2⎤⎥⎦−1|ϵMa|. (30)

Maximum values for are estimated by taking the maximum energy released during giant flares, around erg to the magnetar SGR 1806–20 (2015RPPh…78k6901T, ; 2017ARA&A..55..261K, ). In this case, toroidal magnetic fields of order G arise and from Eqs.  (24) and (29), by making use of fiducial magnetar parameters [ km and M], it follows that . Taking into account hydrogen experiment outcomes (2018EPJC…78..143A, ), one learns that the absolute minimum value for Born-Infeld’s scale field is statvolt cm (or G) (1934RSPSA.144..425B, ) [most recent data for the mass and charge of the electron have been used]. Thus, from the above and Eqs. (30) and (26), we have that and .

Naturally, if larger values of are the case (for instance associated with uncertainties in distance measurements or even newly-born magnetars), then larger ellipticities would emerge. If erg, say erg, which might be possible to giant flare events when they are associated with nonlinear electrodynamics [which increases the magnetic energy reservoir of highly magnetized stars, see Eq. (23)], or are associated with smaller than usual short GRBs 444The authors of Ref. 2005Natur.434.1098H () have estimated that an appreciable quantity of short GRBs might be related to extragalactic unstable magnetars. However, Ref. 2017arXiv171106593T () has cast doubt on whether or not the recent subluminous short gamma-ray burst GRB 170817 A (2017ApJ…848L..13A, ), associated with GW 170817 (2017PhRvL.119p1101A, ), might be related to a magnetar flare. [typical isotropic energy of ordinary short-GRBs are in the range erg (2015JHEAp…7…73D, ; 2017arXiv171106593T, )], then [ G] and . If, instead, only magnetic fields as high as G are possible in magnetars, from Eq. (26) one has that , which in turn would imply in . Thus, upper limits to the magnitude of magnetic ellipticities of magnetars would be around () and the maximum increase of gravitational-wave energy released in Born-Infeld theory with respect to Maxwell’s theory would be around (even larger than ) for flare energy up to () erg.

## 5 Discussions and conclusions

In order to explain magnetars’ observational properties, it seems accepted today in the literature that very large resultant magnetic fields should take place in the vicinities of their surfaces and also interiors. Fields as high as  G should be present in order to function as energy reservoirs to the flare activities magnetars display. With such high fields, it seems reasonable to assume nonlinear electrodynamics could give a more precise description of magnetars and concomitantly they could be natural candidates for probing nonlinear electrodynamics astrophysically. Given that there are known magnetars with low-poloidal magnetic fields, toroidal fields should be particularly relevant for highly magnetized stars. Fields of order  G are a natural scale for magnetic fields in magnetars and this value is of the same order as the scale field in Born-Infeld’s original theory. This was one of our motivations to our analysis.

Axially symmetric toroidal nonlinear fields have a very simple solution in the context of magnetars’ ideal magnetohydrodynamics because electric fields are not induced. Note this would be the case only where toroidal fields would exist and be much larger than poloidal ones, believed to be the case just in the interior of magnetars (2013MNRAS.435L..43C, ). Besides, nonlinear electrodynamics leads to toroidal magnetic field increase when compared to their Maxwellian counterparts. This would imply the increase of the magnetic energy of a magnetar, exactly as we have showed, which could lead to more energetic flare activities. Besides, more magnetic energy could be converted into gravitational-wave energy, as also explicitly showed. The kinematic reason for so is the increase in magnitude of the Born-Infeld magnetic ellipticity when compared to its Maxwellian counterpart.

Hydrogen atom byproducts show that Born-Infeld’s scale field should be larger than the one obtained by Born and Infeld themselves. Actually, such a value defines an absolutely lower limit to the scale field and has to be taken into account for physically relevant constraints. (For the upper limit on which could not have been assessed by the ATLAS experiment, the changes introduced by the Born-Infeld Lagrangian are negligible when compared to the Maxwell theory.) When that is done, Born-Infeld theory predicts an upper limit to the magnetic ellipticity of a magnetar, which also lead to an upper limit to the gravitational-wave energy it may emit. Giant flares are important in this case because they result in the largest energy budgets of the system, which according to the magnetar theory should have exclusive magnetic origins. We have showed that the upper limit to the magnitude of the ellipticity should be within the range for current observations of the magnetar SRG 1806-20. Born-Infeld’s gravitational-wave energy for giant flare events of around erg would be at most larger than their classical counterparts, while for flare energy up to erg the maximum percentage could be much higher than . Naturally, the previous incredibly large Born-Infeld gravitational-wave energy increase should be seen only as an indicative value, given that should be larger than . When LIGO/VIRGO are able to constrain magnetar ellipticities, we see from Eq. (27) that minimum values for the Born-Infeld’s scale field could be inferred astrophysically. This seems very interesting since it would work as a possible cross-check to the already available constraints on and be able to probe the region of parameters current experiments cannot.

It is argued that a fiducial upper limit to the norm of the ellipticity should be around when asymmetries supported by anisotropic stresses built up during the crystallization period of the crust are taken into account (2008PhLB..668….1K, ). It is worth mentioning that Refs. 2016ApJ…831…35D (); 2017EPJC…77..350D () have predicted GW amplitudes for all known pulsars and, when an extremely optimistic case is considered, ellipticities should be at most (for PSR J1846-0258). Thus, since the predicted GW amplitudes are extremely small, observation times of thousands of years would be needed even for advanced detectors such as aLIGO and AdVirgo, and the planned Einstein Telescope might not be able to detect these pulsars. All of the above actually evidences the relevance of finding the maximum increase of gravitational-wave energy in Born-Infeld theory with respect to Maxwell’s theory. The conclusion in the context of Born-Infeld theory is that magnetars with giant flare energy erg would be the most promising candidates for gravitational-wave detections and by consequence potential tests of nonlinear electrodynamics.

In a realistic stellar model, magnetic fields should present both poloidal and toroidal components and should be dependent upon the azimuthal coordinate and time too (besides and ). However, for slowly rotating highly magnetized nonlinear stars it seems reasonable to start with simplified stationary and axially symmetric models, as suggested by some MHD simulations (2013MNRAS.435L..43C, ). Therefore, one could conceive models in the form . Notwithstanding, in the context of nonlinear theories, finding and is expected to be even harder than in Maxwell’s electrodynamics due to the natural coupling of these components. When toroidal fields are very large, though, one could uncouple the system of equations. Besides, it is also pending stability analysis for magnetic fields in the context of nonlinear electrodynamics and when Ohm’s law does not hold. It is already known that Maxwell’s theory in the context of ideal MHD leads to instabilities of purely poloidal or toroidal fields (see e.g., 2011MNRAS.412.1394L (); 2011MNRAS.412.1730L () and references therein). These analyses are interesting on their own and are left for future works.

Additionally, for describing more realistically nonlinear electrodynamics in magnetars and the effects thereof, curved geometries should also be taken into account. However, in a first approach, it seems justifiable to neglect them such that maximal field strength (and derived quantities) and a better understanding of the physics taking place in the above scenario could be obtained. This would allow one to find which aspects are more important to be focused on in more precise analyses. We have seen that magnetic ellipticities are clear candidates. Curved geometries may significantly affect them by changing the spatial distributions of magnetic fields. Besides, the presence of both poloidal and toroidal magnetic fields and their couplings might also cancel out one another effects on the magnetic ellipticity. Due to the connections of the above effects with magnetar GWs, we plan to investigate these issues elsewhere.

Summing up, we have showed that when ideal magnetohydrodynamics and nonlinear electrodynamics are taken into account in magnetars, axially symmetric toroidal fields should be larger than their Maxwellian counterparts. This implies larger magnetic energy and ellipticities, which could thus increase the emission of gravitational waves from nonlinear highly magnetized stars. Current constraints on the Born-Infeld theory, giant-flare energetics and magnetic fields in magnetars point to a maximum increase of of the energy emitted in the form of gravitational waves by Born-Infeld magnetars when compared to Maxwellian ones. When larger giant-flare energy are taken into account, in principle plausible due to the nonlinearities of the electromagnetism which could increase the magnetic energy reservoirs of magnetars, much higher percentages may appear. Thus, the possibility may arise for also probing nonlinear electrodynamics with the advancement of gravitational-wave detectors.

###### Acknowledgements.
J.P.P. acknowledges the financial support given by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) under grants No. 2015/04174-9 and 2017/21384-2. J.G.C. is likewise grateful to the support of FAPESP (2013/15088-0 and 2013/26258-4). R.C.R.L. acknowledges the support of Fundação de Amparo à Pesquisa e Inovação do Estado de Santa Catarina (FAPESC) under grant No. 2017TR1761.

## References

• (1) S.A. Olausen, V.M. Kaspi, ApJS212, 6 (2014).
• (2) V.M. Kaspi, A.M. Beloborodov, Ann. Rev. Astron. Astrophys. 55, 261 (2017).
• (3) R.C. Duncan, C. Thompson, ApJL392, L9 (1992).
• (4) C. Thompson, R.C. Duncan, MNRAS275, 255 (1995).
• (5) J.G. Coelho, D.L. Cáceres, R.C.R. de Lima, M. Malheiro, J.A. Rueda, R. Ruffini, A&A599, A87 (2017).
• (6) N. Rea, P. Esposito, R. Turolla, G.L. Israel, S. Zane, L. Stella, S. Mereghetti, A. Tiengo, D. Götz, E. Göğüş, C. Kouveliotou, Science 330, 944 (2010).
• (7) N. Rea, G.L. Israel, P. Esposito, J.A. Pons, A. Camero-Arranz, R.P. Mignani, R. Turolla, S. Zane, M. Burgay, A. Possenti, S. Campana, T. Enoto, N. Gehrels, E. Göǧüş, D. Götz, C. Kouveliotou, K. Makishima, S. Mereghetti, S.R. Oates, D.M. Palmer, R. Perna, L. Stella, A. Tiengo, ApJ754, 27 (2012).
• (8) N. Rea, D. Viganò, G.L. Israel, J.A. Pons, D.F. Torres, ApJL781, L17 (2014).
• (9) M.A. Alpar, Ü. Ertan, Ş. Çalışkan, ApJL732, L4 (2011).
• (10) J.E. Trümper, K. Dennerl, N.D. Kylafis, Ü. Ertan, A. Zezas, ApJ764, 49 (2013).
• (11) I.F. Malov, Astronomy Reports 54, 925 (2010).
• (12) R.X. Xu, D.J. Tao, Y. Yang, MNRAS373, L85 (2006).
• (13) M. Malheiro, J.A. Rueda, R. Ruffini, PASJ64, 56 (2012)
• (14) J.G. Coelho, M. Malheiro, PASJ66, 14 (2014).
• (15) D.L. Cáceres, S.M. de Carvalho, J.G. Coelho, R.C.R. de Lima, J.A. Rueda, MNRAS465, 4434 (2017).
• (16) S.L. Shapiro, S.A. Teukolsky, Black Holes, White Dwarfs and Neutron Stars: The Physics of Compact Objects (Wiley-VCH, Weinheim, 1986)
• (17) R. Ciolfi, V. Ferrari, L. Gualtieri, MNRAS406, 2540 (2010).
• (18) R. Ciolfi, L. Rezzolla, MNRAS435, L43 (2013).
• (19) R. Ruffini, G. Vereshchagin, S. Xue, Phys. Rep. 487, 1 (2010).
• (20) M. Born, L. Infeld, Royal Society of London Proceedings Series A 144, 425 (1934).
• (21) D.A. Rasheed, ArXiv High Energy Physics - Theory e-prints (1997)
• (22) J.P. Pereira, J.A. Rueda, PRD91(6), 064048 (2015).
• (23) K.A. Bronnikov, PRD63(4), 044005 (2001).
• (24) H. Carley, M.K.H. Kiessling, Physical Review Letters 96(3), 030402 (2006).
• (25) J. Franklin, T. Garon, Physics Letters A 375, 1391 (2011).
• (26) P.N. Akmansoy, L.G. Medeiros, European Physical Journal C 78, 143 (2018).
• (27) J. Ellis, N.E. Mavromatos, T. You, Physical Review Letters 118(26), 261802 (2017).
• (28) L.D. Landau, E.M. Lifshitz, The classical theory of fields (1975)
• (29) D.D. Schnack, Lectures in Magnetohydrodynamics (Springer, Berlin, 2009)
• (30) L. Rezzolla, B.J. Ahmedov, MNRAS352, 1161 (2004).
• (31) R. Belvedere, J.A. Rueda, R. Ruffini, ApJ799, 23 (2015).
• (32) P.A. Davidson, An introduction to magnetohydrodynamics (2001)
• (33) H.J. Mosquera Cuesta, G. Lambiase, J.P. Pereira, PRD95(2), 025011 (2017).
• (34) J.G. Coelho, J.P. Pereira, J.C.N. de Araujo, ApJ823, 97 (2016).
• (35) V.A. De Lorenci, J.P. Pereira, PRA89(4), 043822 (2014).
• (36) J.P. Ostriker, J.E. Gunn, ApJ157, 1395 (1969).
• (37) B.P. Abbott, R. Abbott, T.D. Abbott, M.R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R.X. Adhikari, et al., ApJ839, 12 (2017).
• (38) R. Turolla, S. Zane, A.L. Watts, Reports on Progress in Physics 78(11), 116901 (2015).
• (39) K. Hurley, S.E. Boggs, D.M. Smith, R.C. Duncan, R. Lin, A. Zoglauer, S. Krucker, G. Hurford, H. Hudson, C. Wigger, W. Hajdas, C. Thompson, I. Mitrofanov, A. Sanin, W. Boynton, C. Fellows, A. von Kienlin, G. Lichti, A. Rau, T. Cline, Nature434, 1098 (2005).
• (40) H. Tong, C. Yu, ArXiv e-prints (2017)
• (42) B.P. Abbott, R. Abbott, T.D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R.X. Adhikari, V.B. Adya, et al., Physical Review Letters 119(16), 161101 (2017).
• (43) P. D’Avanzo, Journal of High Energy Astrophysics 7, 73 (2015).
• (44) P.G. Krastev, B.A. Li, A. Worley, Physics Letters B 668, 1 (2008).
• (45) J.C.N. de Araujo, J.G. Coelho, C.A. Costa, ApJ831, 35 (2016).
• (46) J.C.N. de Araujo, J.G. Coelho, C.A. Costa, European Physical Journal C 77, 350 (2017).
• (47) S.K. Lander, D.I. Jones, MNRAS412, 1394 (2011).
• (48) S.K. Lander, D.I. Jones, MNRAS412, 1730 (2011).
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters