The stochastic gravitational wave background from magnetars

Magnetars have already been a potential candidate as gravitational wave sources that could be detected by current and future terrestrial as well as ground based gravitational wave detectors. In this article, we focus on the gravitational wave emission from the distorted rotating neutron stars. The deformation is assumed to be symmetric around an axis that is perpendicular to the rotation axis. The form is applied in the context of a neutron star whose magnetic field has been deformed on its own. By introducing the effects from all magnetars in the Universe, based on various proposed magnetic field configurations, such as poloidal, toroidal; stochastic gravitational wave background (SGWB) can be generated. We choose to figure out exactly how the observations of the stochastic gravitational-wave background should be used to understand much more about physics correlated with the magnetar behaviour, based on the restriction on the eccentricity of the magnetar.


Introduction
The detection of gravitational waves (GWs) was achieved by the Advanced LIGO and VIRGO team [1,2,3,4,5,6,7,8], and the presence of GWs is estimated by Einstein's general relativity.This significant step forward in GW research provides an opportunity to investigate gravity in intensely dynamic and robust regions.A few events, such as black hole-black hole mergers, neutron star-neutron star mergers, or black hole-neutron star mergers [9,10,11], are considered to have a different class of results.When used in conjunction with some electromagnetic contemporaries, it may be possible to better understand physics, particularly in the most extreme regimes of gravitational fields, densities, magnetic fields, and other variables.
The contributions of superposition from many independent and unresolved GW sources are expected to produce the stochastic gravitational wave background (SGWB).The SGWB could be cosmological, arising in various inflationary models [12,13,14] or cosmic string models [15,16,17]; it could be astrophysical because of the superposition of waves generated by many astrophysical sources, e.g., compact binary coalescences (CBCs) [19,18,20,21]; or it could be from neutron stars [22,23] (quadrupole emission, including magnetars) or initial instabilities [24,25,26].Binary neutron star mergers with multi messengers, such as GWs and electromagnetic counterparts of the signals, are studied to estimate the Hubble constant, which measures the rate of expansion of the Universe, and to study the behaviour of the matter which are denser than an atomic nucleus [27].The Hubble constant is estimated to be around 70 km/sec/Mpc.This figure is in line with previous estimates of Planck 2015 results [28,29].To explain the observed features of soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs), magnetars, neutron stars with a strong magnetic field were first proposed [30].The strong magnetic field introduced induced a quadrupolar deformation in the magnetar structure, generating GWs during its rapid spinning, in addition to driving the powerful electromagnetic radiation that enabled observation of these objects.
In the frequency bandwidth of LIGO, VIRGO, and KARGA [2,31,32] interferometric detectors, rapidly rotating neutron stars could be a promising source of consecutive gravitational waves.Second-and third-generation GW detectors could be more beneficial to detect the SGWB caused by magnetars.Extensive studies are in the following references [33,34,35].In Advanced LIGO's second observing run, the LIGO investigated short-and intermediate-duration GW signals from four magnetar bursts.They found no proof of the signal and ascertained a constraint on the root sum square of the maximum dimensionless strain (from incoming intermediate-duration GWs of a span of the spectrum).Based on optimal orientation, the collaborations place upper bounds on the isotropic GW energy for a magnetar with a known distance (SGR 1806220 (8.7 kpc)) [3].Macquet et al. [36] recommended a new GW search algorithm.
A stationary rotating body that is perfectly symmetric does not emit GWs.A pulsar accordingly must deviate from axisymmetry to radiate gravitationally.If an axisymmetric, rotating body radiates gravitationally, its axis of rotation cannot coincide with the axis attributed to the symmetry [37].The rotation axis of a neutron star might not even coincide with a principal axis of its moment of inertia due to its violent formation (supernova) or its environment, and the star may precess [38].If the magnetic axis and the rotational axis do not line up, it is assumed that the magnetic pressure can be the origin of the contortion of the star.The magnetic fields of neutron stars are known to be substantial.Another reason for considering the mechanism of asymmetries is the advancement of non-axisymmetric instabilities.The gravitational radiation reaction [39] or nuclear matter viscosity [40,41] drive rapid rotation in the neutron stars.An internal magnetic field that causes the deformation is also associated with the period derivative ( Ṗ ) of the pulsar [42].SGWB generated from magnetic deformation of newly born magnetars was studied by Cheng et al. [43] in the context of the growing tilt angle between the star's spin and magnetic axis due to the internal viscosity effects.The magnetars with ultra-strong toroidal magnetic fields contribute to the GWs' background spectra.The background radiation above 1000 Hz is suppressed, unlike the background spectra estimated by assuming constant tilt angles at π/2.
Gravitational radiation is produced at two frequencies in the general case: Ω and 2Ω, where Ω is the rotational frequency.Instead of considering analytical expressions of GWs, some have provided the amplitudes h + and h × directly from the Riemann tensor [41].
In addition, to accomplish powerful electromagnetic radiation, the observations of such strongly magnetized compact objects are encouraged.SGWB is formed by the superposition of GWs produced by all the magnetars in the Universe.We classified our study in the following manners: in Section 2, the model of the SGWB has been modified with sophistication.Different magnetic field configurations (poloidal and toroidal field configuration) within the magnetars have been studied based on ellipticities of the objects in Section 3, and we conclude and discuss in Section 4.

SGWB Model
Stochastic gravitational radiation background is a spontaneous GW signal accumulated by many weak, unresolved gravity wave sources.It resembles the cosmic microwave background radiation (CMBR), which appears to be a stochastic background of electromagnetic radiation [44] in several ways.The spectral features of the gravitational background can be characterized in the same fashion as the CMBR by defining how well the energy is distributed in frequency.The normalized GW energy spectrum is commonly used to describe the SGWB [45], as follows.
where dρ gw is the energy density in the range of frequency band f → f + df .ρ c is the critical energy density (today), related with the Hubble constant (H 0 ) as follows: where G and c are the Newtonian gravitational constants and light speed, respectively.Energy (e) propagating over an area (A) in time (t) can be related by energy flux (F ) = e/At.Now, from the definition of the SGWB, the energy density spectrum becomes as follows: Again, the total flux can be written as the integral of fluxes emanating from all sources: Here, R(z) is the rate of magnetars as a function of redshift (z).The gravitational wave energy spectrum is dEgw df and the luminosity distance d L (z) is correlated with the proper distance r(z) as d L (z) = (1 + z)r(z).
The (1 + z) factor accounts for cosmic advancement and transforms the time in the source picture to the detector picture.Eventually, the redshift influence of the comoving volume is captured by E(Ω m , Ω λ , z) as follows: where Ω m and Ω λ are the energy density in matter and in dark energy.R(z) is the rate of magnetars or the GW sources in terms of red shift (z) as observed on Earth, measured in local time.It is correlated with the magnetars' rate per unit comoving volume R V (z), as follows: where R V (z) is measured in the source frame time.Triaxial-shaped rotating neutron stars (NSs) have a time-varying quadrupole moment, which causes them to emit GWs at the double rate of rotational frequency.The complete spectral gravitational energy emitted by an NS [35] with a rotational period of P 0 that slows down due to magnetic dipole torques and GW emission is given by where R is the magnetar radius, B is the magnetic field, I its moment of inertia, and is the ellipticity of the magnetar.The first term in the bracket comes from the electromagnetic radiation's angular acceleration, while the second term comes from the contribution of GW emission acceleration.

Configuration of Magnetars
The magnetic field and the ellipticity of the magnetars are related separately, depending on their field configurations.SGWB searches can be used to examine these various field structures and the physics that also prompt them-in particular, a poloidal magnetic field expanded from the interior to the exterior of the magnetar.A toroidal magnetic field aspect has also been suggested to occur only within the torus-shaped region of the magnetar.

Poloidal Field Configuration
The variables of self-consistent models of magnetized neutron stars can be used to estimate the value of B for distinct models of magnetic field distributions.The magnetic field distribution has a significant impression on the disturbance at a specified magnetic dipole moment.In comparison to the uniformly magnetized ideal conductor case, a stochastic magnetic field or a stellar superconductor interior varies significantly with different parameters, which could lead to GWs.Bonazzola and Gourgoulhon [41] proposed that the ellipticity of the poloidal field configuration is Here, β is a dimensionless quantity.This quantity is taken into consideration in the equation of state and magnetic field geometry.R is the radius of the magnetar, B is the magnetic field associated with the corresponding magnetar, and the moment of inertia of the body is I.A slightly higher ellipticity leads to an adequately high distortion parameter, and the values are almost stagnant for a range of ellipticity, although they increase rapidly for the lower range.The energy spectrum corresponding to the configuration thus depends on the distortion parameters and the associated magnetic field.

Twisted-Torus Magnetic Field Configuration
Braithwaite and Spruit [46] developed the twisted-torus magnetic field configuration.It is claimed that instead of the pure poloidal field configuration, the twisted-torus magnetic field configuration has a universal equilibrium structure of the magnetic field, which is potentially stable.A toroidal magnetic field is closed in the interior of the magnetars in the twisted-torus field structure.However, the deformation is mainly dependent on the magnetic field.Huge deformations can be expected if the magnetic field is intense enough during the early stellar stages.
The poloidal field becomes stabilized and twisted by the toroidal field.The ellipticity can be modeled with the help of a dimensionless parameter, which depends only on the equation of state and field geometry, as follows: According to Ciolfi et al. [47] for the twisted-torus field configuration, for compact stars in the realistic regime, the dimensionless parameter (k) must lie within the range 4 → 9.However, from the Akmal-Pandharipande-Ravenhall equation of state, it is clear k = 4 provides small compactness [48], while Glendnning [49] claimed that the significant compactness could arise around k = 9.The emission of GWs from a distorted star is mainly based on its quadrupole ellipticity.Highly deformed, eccentric stellar bodies offer a lower gravitational wave energy spectrum.Even with a higher radius, they can reduce rapidly.

Conclusions
In this article, the magnetars are categorized based on their ellipticity.We focused on the finite values (not in a shallow regime) of eccentricities and the corresponding parameters in the range of realistic stellar models.The strong magnetic field anticipates the deformation of the quadrupolar distortion in the magnetar structure, developing GWs during its rapid spinning.The sign of the quadrupole ellipticity is essential.< 0 (corresponds to a prolate structure) could be the spin-flip process.It can introduce viscous forces.The angle between the magnetic axis and the rotation axis developed on a dissipation time-scale becomes orthogonal.This procedure would be directly connected to the GW emission [47].The poloidal field appears to produce the oblate stellar configuration, whereas the twisted-torus field seems to make prolate.The ellipticity has always been positive for the poloidal field and dominates over the toroidal field.As a result, the spin-flip formalism is incompatible with twisted-torus structures.
If the interior of the magnetars is considered to be a type I superconductor, emitting a magnetic field from a part of the star, the distortion parameters are much higher.In such cases, more significant distortions are expected, with an overall low magnetic dipole moment.Extremely high stresses in the star's crust could also contribute to higher β values in the type-II superconducting interior.In the context of a poloidally dominated magnetic field structure, it is thus essential to incorporate SGWB constraints on B and thus obtain information about the equation of state of the magnetars.
Variations of the dimensionless parameter (β and k, respectively) of the poloidal field configuration and twisted-torus field configuration with the ellipticity are shown in Figure 1 and Figure 3, respectively.In both configurations, the values of β and k are effectively very low for very small values of ellipticity (< 10 −3 ) and become almost saturated with higher values of .The blue band in Figure 3 defines the realistic regime of k for the compact stars.Based on Ciolfi et al.'s [47] claim that 4 ≤ k ≤ 9, the variation of ellipticity is shown for different values of magnetic field (B); Figure 4.For B = 1.0 × 10 17 G, the of the realistic compact stars has a range 0.0402 → 0.0905, whereas the range is 0.0102 → 0.0227 for B = 0.5 × 10 17 G, and 0.0037 ≤ ≤ 0.0082 for B = 0.3 × 10 17 G.From the above analysis, it is clear that higher deformation requires a hierarchically higher magnetic field.
Variation of Ω gw (f ) based on the free parameters and dimensionless parameter (β and k, respectively) for the poloidal field configuration and twisted-torus field configuration are shown in Figure 2 and Figure 4, respectively.It is found that the normalized energy spectrum for the twisted-torus configuration is higher compared to the poloidal configuration.It makes the possibility of getting twisted-torus configured magnetars higher than the poloidal configuration.It is also interesting to note that the Ω gw decreases sufficiently for a higher radius.
Figure 5 also recommends that the energy spectrum for twisted-torus field configuration be higher than the poloidal field configuration.It registers that for the same structural parameter, twisted-torus field configured magnetars can easily be found.
The predicted energy spectrum is less than the so far detected energy spectrum.Nevertheless, for both the configurations, higher ellipticity makes the system more plausible to find.Recently, the collaborations have been searching for the quasi- The expected evolution of the GW energy spectra emitted by the magnetars for different field configuration with frequency for a specific parameter's space.Overall, the energy spectrum tends to increase rapidly as frequency increases.The orange and blue lines represent the results calculated by involving parameters I = 10 43 kg km 2 , R = 10 km, = 0.02.monochromatic gravitational wave signals from the energetic young X-ray pulsar PSR J0537-6910 [1].The upcoming presentations of the Advanced LIGO, VIRGO, or KAGRA are awaited to detect the GWs from magnetars.This study can finally affirm that the SGWB energy spectrum for the magnetars with toroidal configuration will emit more detectable GWs than the poloidal configuration.

Figure 1 .
Figure 1.Expected evolution of the dimensionless parameter (β) with ellipticity ( ) for the free parameters radius (R), magnetic field (B), and moment of inertia (I).The parametric valuesare exhibited in the figure.

Figure 2 .
Figure2.Ω gw (f ) for different parametric preferences for the magnetar, produced SGWB high frequency model of the poloidal field configuration.The free parameters are and R, respectively.Furthermore, R is in km, as shown in the figure.

Figure 3 .
Figure 3.The dimensionless parameter (k) versus the ellipticity ( ) for the magnetars with different strengths.The only free parameter here is magnetic field (B), same as in Figure 1.The parametric values are shown in the figure.The range of for being a realistic stellar model for each strength is marked by dashed lines.

Figure 4 .
Figure 4. Variation of Ω gw (f ) for different parametric preferences for the magnetar producing SGWB from the high-frequency model of the twisted-torus magnetic field configuration.Here, and R are the free parameters, respectively, and R are measured in km.Values of are chosen within the realist stellar range referred in Figure 3.
Figure 5.The expected evolution of the GW energy spectra emitted by the magnetars for different field configuration with frequency for a specific parameter's space.Overall, the energy spectrum tends to increase rapidly as frequency increases.The orange and blue lines represent the results calculated by involving parameters I = 10 43 kg km 2 , R = 10 km, = 0.02.