Search

Magnetar vibrational modes are theorized to be associated with energetic X-ray flares. Regular searches for gravitational waves from these modes have been performed by Advanced LIGO (Laser Interferometer Gravitational-wave Observatory) and Advanced Virgo, with no detections so far. Presently, search results are given in upper limits on the root sum square of the integrated gravitational-wave strain. However, the increased sensitivity of current detectors and the promise of future detectors invite the consideration of more astrophysically motivated methods. We present a framework for augmenting gravitational-wave searches to measure or place direct limits on magnetar astrophysical properties in various search scenarios using a set of phenomenological and analytical models. [ABSTRACT FROM AUTHOR]

While the appeal of their extraordinary radio luminosity to our curiosity is undiminished, the nature of fast radio bursts (FRBs) has remained unclear. The challenge has been due in part to small sample sizes and limited understanding of telescope selection effects. We here present the first inclusion of the entire set of one-off FRBs from CHIME/FRB Catalog 1 in frbpoppy. Where previous work had to curate this data set, and fit for few model parameters, we have developed full multi-dimensional Markov chain Monte Carlo (MCMC) capabilities for frbpoppy – the comprehensive, open-science FRB population synthesis code – that allow us to include all one-off CHIME bursts. Through the combination of these two advances we now find the best description of the real, underlying FRB population, with higher confidence than before. We show that 4 ± 3 × 103 one-off FRBs go off every second between Earth and z = 1; and we provide a mock catalog based on our best model, for straightforward inclusion in other studies. We investigate CHIME side-lobe detection fractions, and FRB luminosity characteristics, to show that some bright, local FRBs are still being missed. We find strong evidence that FRB birth rates evolve with the star formation rate of the Universe, even with a hint of a short (0.1−1 Gyr) delay time. The preferred contribution of the hosts to the FRB dispersion agrees with a progenitor birth location in the host disk. This population-based evidence solidly aligns with magnetar-like burst sources, and we conclude FRBs are emitted by neutron stars. [ABSTRACT FROM AUTHOR]

The thermal evolution of isolated neutron stars is a key element in unravelling their internal structure and composition and establishing evolutionary connections among different observational subclasses. Previous studies have predominantly focused on one-dimensional or axisymmetric two-dimensional models. In this study, we present the thermal evolution component of the novel three-dimensional magnetothermal code MATINS (MAgneto-Thermal evolution of Isolated Neutron Star). MATINS employs a finite volume scheme and integrates a realistic background structure, along with state-of-the-art microphysical calculations for the conductivities, neutrino emissivities, heat capacity, and superfluid gap models. This paper outlines the methodology employed to solve the thermal evolution equations in MATINS , along with the microphysical implementation that is essential for the thermal component. We test the accuracy of the code and present simulations with non-evolving magnetic fields of different configurations (all with electrical currents confined to the crust and a magnetic field that does not thread the core), to produce temperature maps of the neutron star surface. Additionally, for a specific magnetic field configuration, we show one fully coupled evolution of magnetic field and temperature. Subsequently, we use a ray-tracing code to link the neutron star surface temperature maps obtained by MATINS with the phase-resolved spectra and pulsed profiles that would be detected by distant observers. This study, together with our previous article focused on the magnetic formalism, presents in detail the most advanced evolutionary code for isolated neutron stars, with the aim of comparison with their timing properties, thermal luminosities and the associated X-ray light curves. [ABSTRACT FROM AUTHOR]

Utilizing four archival X-ray data sets taken with the Hard X-ray Detector onboard Suzaku , timing studies were performed on three magnetars, 1E 1841−045 (observed in 2006), SGR 0501+4516 (2008), and 1RXS J170849.0−400910 (2009 and 2010). Their pulsations were reconfirmed, typically in an energy range of 12–50 keV. The 11.783 s pulses of 1E 1841−045 and those of SGR 0501+4516 at 5.762 s were periodically phase modulated, with a long period of |$\approx 23.4$| and |$\approx 16.4$|  ks, respectively. The pulse-phase modulation was also observed, at |$\approx 46.5$|  ks, from two data sets of 1RXS J170849.0−400910. In all these cases, the modulation amplitude was 6 per cent to 16 per cent of the pulse cycle. Including previously confirmed four objects, this characteristic timing behaviour is now detected from seven magnetars in total, and interpreted as a result of free precession of neutron stars that are deformed to an asphericity of |$\sim 10^{-4}$|⁠. Assuming that the deformation is due to magnetic stress, these magnetars are inferred to harbour toroidal magnetic fields of |$B_{\rm t}\sim 10^{16}$|  G. By comparing the estimated |$B_{\rm t}$| of these objects with their poloidal dipole field |$B_{\rm d}$|⁠ , the |$B_{\rm t}/B_{\rm d}$| ratio is found to increase with their characteristic age. Therefore, the toroidal fields of magnetars are likely to last longer than their poloidal fields. This explains the presence of some classes of neutron stars that have relatively weak |$B_{\rm d}$| but are suspected to hide strong |$B_{\rm t}$| inside them. [ABSTRACT FROM AUTHOR]

Four X-ray datasets of the soft gamma repeater SGR 1806−20, taken with the Gas Imaging Spectrometer (GIS) onboad ASCA, were analyzed. Three of them were acquired over 1993 October 9–20, the last one in 1995 October. Epoch-folding analysis of the 2.8–12 keV signals confirmed the ∼7.6 s pulses in these data, which Kouveliotou et al. (1998, Nature, 393, 235) reported as one of the earliest pulse detections from this object. In the 1995 observation, 3–12 keV pulses were phase modulated with a period of T  = 16.4 ± 0.4 ks, and an amplitude of ∼1 s. This makes a fourth example of the behavior observed from magnetars. As in the previous three sources, the pulse-phase modulation of SGR 1806−20 disappeared at ≲2.5 keV, where the soft X-ray component dominates. In the 1993 datasets, this periodic modulation was reconfirmed, and successfully phase-connected coherently across the 11 d interval. As a result, the modulation period was refined to T  = 16.435 ± 0.024 ks. The implied high stability of the phenomenon strengthens its interpretation in terms of free precession of the neutron star, which is deformed to an asphericity of ∼10−4, presumably by the stress of toroidal magnetic fields reaching ∼1016 G. Toroidal fields of this level can be common among magnetars. [ABSTRACT FROM AUTHOR]

Magnetars are slowly rotating neutron stars that possess the strongest magnetic fields known in the cosmos (1014 − 1015 G). They display a range of transient high-energy electromagnetic activity. The brightest and most energetic of these events are the gamma-ray bursts (GRBs) known as magnetar giant flares (MGFs), with isotropic energies Eiso ≈ 1044 − 1046 erg. Only seven MGF detections have been made to date: three unambiguous events occurred in our Galaxy and the Magellanic Clouds, and the other four MGF candidates are associated with nearby star-forming galaxies. As all seven identified MGFs are bright at Earth, additional weaker events likely remain unidentified in archival data. We conducted a search of the Fermi Gamma-ray Burst Monitor database for candidate extragalactic MGFs and, when possible, collected localization data from the Interplanetary Network (IPN) satellites. Our search yielded one convincing event, GRB 180128A. IPN localizes this burst within NGC 253, commonly known as the Sculptor Galaxy. The event is the second MGF in modern astronomy to be associated with this galaxy and the first time two bursts have been associated with a single galaxy outside our own. Here we detail the archival search criteria that uncovered this event and its spectral and temporal properties, which are consistent with expectations for a MGF. We also discuss the theoretical implications and finer burst structures resolved from various binning methods. Our analysis provides observational evidence of an eighth identified MGF. [ABSTRACT FROM AUTHOR]

Core-collapse supernovae (CCSNe) are considered the primary magnetar formation channel, with 15 magnetars associated with supernova remnants (SNRs). A large fraction of these should occur in massive stellar binaries that are disrupted by the explosion, meaning that |$\sim 45~{{\ \rm per\ cent}}$| of magnetars should be nearby high-velocity stars. Here, we conduct a multiwavelength search for unbound stars, magnetar binaries, and SNR shells using public optical (uvgrizy bands), infrared (J, H, K , and Ks bands), and radio (888 MHz, 1.4 GHz, and 3 GHz) catalogues. We use Monte Carlo analyses of candidates to estimate the probability of association with a given magnetar based on their proximity, distance, proper motion, and magnitude. In addition to recovering a proposed magnetar binary, a proposed unbound binary, and 13 of 15 magnetar SNRs, we identify two new candidate unbound systems: an OB star from the Gaia catalogue we associate with SGR J1822.3−1606, and an X-ray pulsar we associate with 3XMM J185246.6 + 003317. Using a Markov Chain Monte Carlo simulation that assumes all magnetars descend from CCSNe, we constrain the fraction of magnetars with unbound companions to |$5\lesssim f_u \lesssim 24~{{\ \rm per\ cent}}$|⁠ , which disagrees with neutron star population synthesis results. Alternate formation channels are unlikely to wholly account for the lack of unbound binaries as this would require |$31\lesssim f_{nc} \lesssim 66~{{\ \rm per\ cent}}$| of magnetars to descend from such channels. Our results support a high fraction (⁠|$48\lesssim f_m \lesssim 86~{{\ \rm per\ cent}}$|⁠) of pre-CCSN mergers, which can amplify fossil magnetic fields to preferentially form magnetars. [ABSTRACT FROM AUTHOR]

It was recently discovered that the time correlations of repeating fast radio bursts (FRBs) are similar to earthquake aftershocks. Motivated by the association between FRBs and magnetars, here we report correlation function analyses in the time-energy space for the 563 periodic radio pulses and the 579 X-ray short bursts from the magnetar SGR 1935+2154, which is known to have generated FRBs. Although radio pulses are concentrated near the fixed phase of the rotational cycle, we find that when multiple pulses occur within a single cycle, their correlation properties (aftershock production probability, aftershock rate decaying in power of time, and more) are similar to those of extragalactic FRBs and earthquakes. A possible interpretation is that the radio pulses are produced by rupture of the neutron star crust, and the first pulse within one cycle is triggered by external force periodically exerted on the crust. The periodic external force may be from the interaction of the magnetosphere with material ejected in an outburst. For X-ray bursts, we found no significant correlation signal, though correlation on the same time-scale as radio pulses may be hidden due to the long event duration. The aftershock similarity between the periodic radio pulsation and FRBs is surprising, given that the two are energetically very different, and therefore the energy sources would be different. This suggests that the essence of FRB-like phenomena is starquakes, regardless of the energy source, and it is important to search for FRB-like bursts from neutron stars with various properties or environments. [ABSTRACT FROM AUTHOR]

Recently, there has been renewed interest in the Be X-ray binary (Be/XRB) SAX J1324.4−6200 because of its spatial coincidence with a variable γ-ray source detected by Fermi/LAT. To explore more thoroughly its properties, new X-ray observations were carried out in 2023 by NuSTAR, XMM-Newton, and Swift satellites, jointly covering the energy range from 0.2 − 79 keV. SAX J1324.4−6200 was caught at an X-ray flux of ∼10−11 erg cm−2 s−1. The X-ray spectrum fits well with an absorbed power law with a high energy cutoff. Other acceptable fits require an additional blackbody component (kTbb ≈ 1.1 keV) or a Gaussian in absorption (Egabs ≈ 6.9 keV). We measured a NuSTAR spin period of 175.8127 ± 0.0036 s and an XMM-Newton spin period of 175.862 ± 0.025 s. We show that all the available spin period measurements of SAX J1324.4−6200, spanning 29 yr, are highly correlated with time, resulting in a remarkably stable spin-down of Ṗ = 6.09 ± 0.06 × 10−9 s s−1. We find that if SAX J1324.4−6200 hosts an accretion-powered pulsar, accretion torque models indicate a surface magnetic field of ∼1012 − 13 G. The X-ray properties emerging from our analysis strenghten the hypothesis that SAX J1324.4−6200 belongs to the small group of persistent Be/XRBs. We also performed radio observations with the Parkes Murriyang telescope, to search for radio pulsations. However, no radio pulsations compatible with the rotational ephemeris of SAX J1324.4−6200 were detected. We rule out the hypothesis that SAX J1324.4−6200 is a γ-ray binary where the emission is produced by interactions between the pulsar and the companion winds. Other models commonly used to account for the production of γ-rays in accreting pulsars cannot reproduce the bright emission from SAX J1324.4−6200. We examined other possible mechanisms behind the γ-ray emission and note that there is a ∼0.5% chance probability that an unknown extragalactic active galactic nucleus (AGN) observed through the Galactic plane may coincidentally fall within the Fermi/LAT error circle of the source and be responsible for the γ-ray emission. [ABSTRACT FROM AUTHOR]

Fast radio bursts (FRBs) can present a variety of polarization properties and some of them are characterized by narrow spectra. In this work, we study spectral properties from the perspective of intrinsic radiation mechanisms and absorption through the waves propagating in the magnetosphere. The intrinsic radiation mechanisms are considered by invoking quasi-periodic bunch distribution and perturbations on charged bunches moving on curved trajectories. The narrowband emission is likely to reflect some quasi-periodic structure on the bulk of bunches, which may be due to quasi-periodically sparking in a "gap" or quasi-monochromatic Langmuir waves. A sharp spike would appear in the spectrum if the perturbations were to induce a monochromatic oscillation of bunches; however, it is difficult to create a narrow spectrum because the Lorentz factor has large fluctuations, so the spike disappears. Both the bunching mechanism and perturbations scenarios share the same polarization properties, with a uniformly distributed bulk of bunches. We investigated the absorption effects, including Landau damping and curvature self-absorption in the magnetosphere, which are significant at low frequencies. Subluminous O-mode photons cannot escape from the magnetosphere due to the Landau damping, leading to a height-dependent lower frequency cut-off. The spectra can be narrow when the frequency cut-off is close to the characteristic frequency of curvature radiation, however, such conditions cannot always be met. The spectral index is 5/3 at low-frequency bands due to the curvature self-absorption is not as steep as what is seen in observations. The intrinsic radiation mechanisms are more likely to generate the observed narrow spectra of FRBs, rather than the absorption effects. [ABSTRACT FROM AUTHOR]

After a successful supernova, a proto-neutron star (PNS) cools by emitting neutrinos on ∼1–100 s time-scales. Provided that there are neutrino emission 'hotspots' or 'cold-spots' on the surface of the rotating PNS, we can expect a periodic modulation in the number of neutrinos observable by detectors. We show that Fourier transform techniques can be used to determine the PNS rotation rate from the neutrino arrival times. Provided there is no spin-down, a 1-parameter Discrete Fourier Transform (DFT) is sufficient to determine the spin period of the PNS. If the PNS is born as a magnetar with polar magnetic field strength B 0 ≳ 1015 G and is 'slowly' rotating with an initial spin period ≳100 ms, then it can spin-down to periods of the order of seconds during the cooling phase. We propose a modified DFT technique with three frequency parameters to detect spin-down. Due to lack of neutrino data from a nearby supernova except the ∼20 neutrinos detected from SN1987A, we use toy models and one physically motivated modulating function to generate neutrino arrival times. We use the false alarm rate (FAR) to quantify the significance of the Fourier power spectrum peaks. We show that PNS rotation and spin-down are detected with |$\rm FAR\,\lt\, 2~{{\ \rm per\ cent}}$| (2σ) for periodic signal content |$\rm M\gtrsim 13-15~{{\ \rm per\ cent}}$| if 5 × 103 neutrinos are detected in ∼3 s and with |$\rm FAR\,\lt\, 1{{\ \rm per\ cent}}$| for |$\rm M\,\ge 5{{\ \rm per\ cent}}$| if 5 × 104 neutrinos are detected in ∼3 s. Since we can expect ∼104−105 neutrino detections from a supernova at 10 kpc, detection of PNS rotation and spin-down is possible using the neutrinos from the next Galactic supernova. [ABSTRACT FROM AUTHOR]

Neutron star magnetic field evolution is mediated through the Hall effect and Ohmic dissipation in the crust while ambipolar diffusion is taking place in the core. These effects have been studied in detail in either part of the star, however, their combined, simultaneous evolution and interplay has not been explored in detail yet. Here, we present simulation results of the simultaneous evolution of the magnetic field in the core due to ambipolar diffusion and the crust due to Hall effect and Ohmic decay, under the assumption of axial symmetry. We find that a purely poloidal field generates a toroidal field in the crust, due to the Hall effect, that sinks into the core. A purely toroidal field remains toroidal and spreads into the core and the crust. Finally, for a mixed poloidal–toroidal field, the north–south symmetry is broken due to the Hall effect in the crust, however, ambipolar diffusion, tends to restore it. We examine the role of ambipolar diffusion to the magnetic field decay and we compare the rate of the conversion of magnetic field energy into heat, finding that it enhances the magnetic field decay in neutron stars. [ABSTRACT FROM AUTHOR]

The shape and polarization properties of the radio pulse profiles of radio-loud magnetars provide a unique opportunity to investigate their magnetospheric properties. Gaussian process regression analysis was used to investigate the variation in the total intensity shape of the radio pulse profiles of the magnetar Swift J1818.0–1607. The observed profile shape was found to evolve through three modes between MJDs 59104 and 59365. The times at which these transitions occurred coincided with changes in the amplitude of modulations in the spin-down rate. The amount of linear and circular polarization was also found to vary significantly with time. Lomb–Scargle periodogram analysis of the spin-down rate revealed three possibly harmonically related frequencies. This could point to the magnetar experiencing seismic activity. However, no profile features exhibited significant periodicity, suggesting no simple correlations between the profile variability and fluctuations of the spin-down on shorter time-scales within the modes. Overall, this implies that the mode changes seen are a result of local magnetospheric changes, with other theories, such as precession, less able to explain these observations. [ABSTRACT FROM AUTHOR]

At the high temperatures inside recently formed neutron stars (⁠|$T\gtrsim 5\times 10^{8}\, \text{K}$|⁠), the particles in their cores are in the 'strong-coupling' regime, in which collisional forces make them behave as a single, stably stratified, and thus non-barotropic fluid. In this regime, axially symmetric hydromagnetic quasi-equilibrium states are possible, which are only constrained to have a vanishing azimuthal Lorentz force. In these states, the particle species deviate from chemical (β) equilibrium, which tends to be restored by β decays (Urca reactions), inducing fluid motions that change the magnetic field configuration. If the stars remained hot for a sufficiently long time, this evolution would eventually lead to a chemical equilibrium state, in which the fluid is barotropic and the magnetic field, if axially symmetric, satisfies the non-linear Grad–Shafranov equation. Here, we present a numerical scheme that decouples the magnetic and thermal evolution, enabling to efficiently perform, for the first time, long-term magnetothermal simulations in this regime for different magnetic field strengths and geometries. Our results demonstrate that, even for magnetar-strength fields |$\gtrsim 10^{16} \, \mathrm{G}$|⁠ , the feedback from the magnetic evolution on the thermal evolution is negligible. Thus, as the core passively cools, the Urca reactions quickly become inefficient at restoring chemical equilibrium, so the magnetic field evolves very little, and the Grad–Shafranov state is not attained. Therefore, any substantial evolution of the core magnetic field must occur later, in the 'weak-coupling' regime (⁠|$T\lesssim 5\times 10^8 \, \mathrm{K}$|⁠), when Urca reactions are frozen and ambipolar diffusion becomes relevant. [ABSTRACT FROM AUTHOR]

In some cases, the merger of two neutron stars can produce a rapidly rotating and highly magnetized millisecond magnetar. A significant proportion of the rotational energy deposited to the emerging ejecta can produce a late-time radio brightening from interacting with the ambient medium. Detection of this late-time radio emission from short GRBs can have profound implications for understanding the physics of the progenitor. We report the radio observations of five short GRBs – 050709, 061210, 100625A, 140903A, and 160821B using the legacy Giant Metrewave Radio Telescope (GMRT) at 1250, 610, and 325 MHz frequencies and the upgraded-GMRT (uGMRT) at band 5 (1050–1450 MHz) and band 4 (550–900 MHz) after ∼2–11 yr from the time of the burst. The GMRT observations at low frequencies are particularly important to detect the signature of merger ejecta emission at the peak. These observations are the most delayed searches associated with some GRBs for any late-time low-frequency emission. We find no evidence for such an emission. We find that none of these GRBs is consistent with maximally rotating magnetar with a rotational energy of |$\sim 10^{53}\, {\rm erg}$|⁠. However, magnetars with lower rotational energies cannot be completely ruled out. Despite the non-detection, our study underscores the power of radio observations in the search for magnetar signatures associated with short GRBs. However, only future radio observatories may be able to detect these signatures or put more stringent constraints on the model. [ABSTRACT FROM AUTHOR]

Kilonovae are optical transients following the merger of neutron star binaries, which are powered by the r -process heating of merger ejecta. However, if a merger remnant is a long-lived supramassive neutron star supported by its uniform rotation, it will inject energy into the ejecta through spin-down power. The energy injection can boost the peak luminosity of a kilonova by many orders of magnitudes, thus significantly increasing the detectable volume. Therefore, even if such events are only a small fraction of the kilonova population, they could dominate the detection rates. However, after many years of optical sky surveys, no such event has been confirmed. In this work, we build a boosted kilonova model with rich physical details, including the description of the evolution and stability of a proto neutron star, and the energy absorption through X-ray photoionization. We simulate the observation prospects and find the only way to match the absence of detection is to limit the energy injection by the newly born magnetar to only a small fraction of the neutron star rotational energy, thus they should collapse soon after the merger. Our result indicates that most supramassive neutron stars resulting from binary neutron star mergers are short lived and they are likely to be rare in the Universe. [ABSTRACT FROM AUTHOR]

This study was focused on the investigation of a magnetic field penetrating from the core of a neutron star to its surface. The range of possible field configurations in the intermediate solid crust is less limited owing to the elastic force acting on the force balance. When the Lorentz force is excessively strong, the magnetoelastic equilibrium does not hold, and thus, the magnetic field becomes constrained. By numerically solving for the magnetoelastic equilibrium in a thin crust, the range of the magnetic field at the core–crust interface was determined, while assuming the exterior to be fixed as a dipole in vacuum. The results revealed that the toroidal component should be smaller than the poloidal component at the core–crust interface for the surface dipole, B 0 > 2.1 × 1014 G. Consequently, a strong toroidal field, for example, B ∼ 1016 G, as suggested by free precession of magnetars should be confined to a deep interior core and should be reduced to B ∼ 1014 G at the bottom of the crust. The findings of this study provide insights into the interior field structure of magnetars. Further investigations on more complicated geometries with higher multipoles and exterior magnetosphere are necessary. [ABSTRACT FROM AUTHOR]

Fast Radio Bursts (FRBs) are extra-galactic origin milli-second duration bright radio bursts. Theoretically, FRBs may produce optical counterparts with durations from milliseconds to hours. The FRB optical counterparts may be detectable in future large field telescopes, including the China Space Station Telescope (CSST), the 2.5-meter Wide Field Survey Telescope (WFST) lead by the University of Science and Technology of China (USTC) and the Purple Mountain Observatory (PMO), and the Earth 2.0 (ET). The fast radio burst optical counterparts are grouped into millisecond time-scale optical counterparts, hourly time-scale optical counterparts, and optical afterglow for our study. The first two can be generated by the high-energy extension of the radio radiation of fast radio bursts and the inverse Compton scattering of high-energy electrons. The event rates highly depend on the optical-to-radio flux ratio η ν. For millisecond duration optical counterparts, the detection rate of WFST, CSST, and ET can reach hundreds per year in an ideal case. If η ν ∼ 10 − 3 , the corresponding annual detection rates of WFST and CSST are in the order of 1, and the annual detection rate of ET is 19.5. For the hourly timescale optical counterparts, ideally, the age of the supernova remnant is 5 years, η ν is about 10 − 6 , and the annual detection rates are above 100. The X-ray counterpart of FRB 200428 indicates that FRBs may produce relativistic outflow, which will interact with the interstellar medium to produce optical afterglows. Combined with the standard afterglow model, the detectability of optical afterglow is explored with a simulation of fast radio bursts following the redshift and energy distribution from the literature. With a total energy-radio energy ratio similar to FRB 200428, (ζ = 10 5), the estimated annual detection rates of CSST, WFST, and ET are 1.3, 1.0, and 67, respectively. [ABSTRACT FROM AUTHOR]

We present here timing and spectral analyses of SGR J1830−0645 based on an AstroSat observation carried out on 2020 October 16, about a week after the onset of its first detected X-ray outburst. Using data taken with the Soft X-ray Telescope (SXT) and Large Area X-ray Proportional Counter (LAXPC), we have detected 0.9–10 keV coherent pulsations at a period of ∼10.4 s. The pulse profiles were single-peaked, asymmetric, and consisted of minor peaks attributable to hotspots on the neutron star surface. The pulsed fraction evolved significantly with energy, increasing to energies around 5 keV with a steep drop thereafter. The 0.9–25 keV SXT–LAXPC energy spectrum is best described with two thermal components having temperatures ∼0.46 and ∼1.1 keV (emission radii of ∼2.4 and ∼0.65 km, respectively, assuming a distance of 4 kpc) along with a power-law component having a photon index of ∼0.39. We report the detection of 67 X-ray bursts having an average duration of ∼33 ms. The brightest burst lasted for about 90 ms and had a 3–25 keV fluence of ∼5 × 10−9 erg cm−2. [ABSTRACT FROM AUTHOR]

We report the serendipitous discovery of an extremely intermittent radio pulsar, PSR J1710−3452, with a relatively long spin period of 10.4 s. The object was discovered through the detection of 97 bright radio pulses in only one out of 66 epochs of observations spanning almost three years. The bright pulses have allowed the source to be localized to a precision of 0.5 arcsec through radio imaging. We observed the source location with the Swift X-ray telescope but did not detect any significant X-ray emission. We did not identify any high-energy bursts or multifrequency counterparts for this object. The solitary epoch of detection hinders the calculation of the surface magnetic field strength, but the long period and the microstructure in the single-pulses resembles the emission of radio-loud magnetars. If this is indeed a magnetar, it is located at a relatively high Galactic latitude (2.9°), making it potentially one of the oldest and the most intermittent magnetars known in the Galaxy. The very short activity window of this object is unique and may point towards a yet undetected population of long period, highly transient radio emitting neutron stars. [ABSTRACT FROM AUTHOR]

We study dynamics of relativistic coronal mass ejections (CMEs), from launching by shearing of foot-points (either slowly – the 'Solar flare' paradigm, or suddenly – the 'star quake' paradigm), to propagation in the preceding magnetar wind. For slow shear, most of the energy injected into the CME is first spent on the work done on breaking through the overlaying magnetic field. At later stages, sufficiently powerful CMEs may lead to the 'detonation' of a CME and opening of the magnetosphere beyond some equipartition radius r eq, where the decreasing energy of the CME becomes larger than the decreasing external magnetospheric energy. Post-CME magnetosphere relaxes via the formation of a plasmoid-mediated current sheet, initially at ∼ r eq, and slowly reaching the light cylinder. Both the location of the foot-point shear and the global magnetospheric configuration affect the frequent/weak versus rare/powerful CME dichotomy – to produce powerful flares, the slow shear should be limited to field lines that close in near the star. After the creation of a topologically disconnected flux tube, the tube quickly (at ∼ the light cylinder) comes into force-balance with the preceding wind and is passively advected/frozen in the wind afterward. For fast shear (a local rotational glitch), the resulting large amplitude Alfvén waves lead to the opening of the magnetosphere (which later recovers similarly to the slow shear case). At distances much larger than the light cylinder, the resulting shear Alfvén waves propagate through the wind non-dissipatively. [ABSTRACT FROM AUTHOR]

During a systematic search for new X-ray pulsators in the XMM–Newton archive, we discovered a high amplitude (⁠|$PF\simeq 86~{{\ \rm per\ cent}}$|⁠) periodic (⁠|$P\simeq 7.25\, \mathrm{s}$|⁠) modulation in the X-ray flux of 4XMM J045626.3–694723 (J0456 hereafter), a previously unclassified source in the Large Magellanic Cloud (LMC). The period of the modulation is strongly suggestive of a spinning neutron star (NS). The source was detected only during one out of six observations in 2018–2022. Based on an absorbed power-law spectral model with photon slope of Γ ≃ 1.9, we derive a 0.3–10 keV luminosity of |$L_\mathrm{X}\simeq 2.7\times 10^{34}\, \mathrm{erg}\, \mathrm{s}^{-1}$| for a distance of 50 kpc. The X-ray properties of J0456 are at variance with those of variable LMC X-ray pulsars hosted in high-mass X-ray binary systems with a Be-star companion. Based on Southern African Large Telescope (SALT) spectroscopic observations of the only optical object that matches the X-ray uncertainty region, we cannot completely rule out that J0456 is an NS accreting from a late-type (G8-K3) star, an as-yet-unobserved binary evolutionary outcome in the Magellanic Clouds (MCs). We show that the source properties are in better agreement with those of magnetars. J0456 may thus be the second known magnetar in the LMC after SGR 0526–66. [ABSTRACT FROM AUTHOR]

Fast radio bursts (FRBs) are mysterious bursts in the millisecond time-scale at radio wavelengths. Currently, there is little understanding about the classification of repeating FRBs, based on difference in physics, which is of great importance in understanding their origin. Recent works from the literature focus on using specific parameters to classify FRBs to draw inferences on the possible physical mechanisms or properties of these FRB subtypes. In this study, we use publicly available 1652 repeating FRBs from FRB 121102 detected with the Five-hundred-metre Aperture Spherical Telescope (FAST), and studied them with an unsupervised machine learning model. By fine-tuning the hyperparameters of the model, we found that there is an indication for four clusters from the bursts of FRB 121102 instead of the two clusters ('Classical' and 'Atypical') suggested in the literature. Wherein, the 'Atypical' cluster can be further classified into three sub-clusters with distinct characteristics. Our findings show that the clustering result we obtained is more comprehensive not only because our study produced results which are consistent with those in the literature but also because our work uses more physical parameters to create these clusters. Overall, our methods and analyses produced a more holistic approach in clustering the repeating FRBs of FRB 121102. [ABSTRACT FROM AUTHOR]

Using physics-informed neural networks (PINNs) to solve a specific boundary value problem is becoming more popular as an alternative to traditional methods. However, depending on the specific problem, they could be computationally expensive and potentially less accurate. The functionality of PINNs for real-world physical problems can significantly improve if they become more flexible and adaptable. To address this, our work explores the idea of training a PINN for general boundary conditions and source terms expressed through a limited number of coefficients, introduced as additional inputs in the network. Although this process increases the dimensionality and is computationally costly, using the trained network to evaluate new general solutions is much faster. Our results indicate that PINN solutions are relatively accurate, reliable, and well behaved. We applied this idea to the astrophysical scenario of the magnetic field evolution in the interior of a neutron star connected to a force-free magnetosphere. Solving this problem through a global simulation in the entire domain is expensive due to the elliptic solver's needs for the exterior solution. The computational cost with a PINN was more than an order of magnitude lower than the similar case solved with a finite difference scheme, arguably at the cost of accuracy. These results pave the way for the future extension to three-dimensional of this (or a similar) problem, where generalized boundary conditions are very costly to implement. [ABSTRACT FROM AUTHOR]

Swift J1818.0−1607 is a radio-emitting magnetar that was discovered in X-ray outburst in 2020 March. Starting 4 d after this outburst, we began a nearly 5-month multifrequency observing campaign at 2.2, 8.4, and 32 GHz using telescopes in the NASA Deep Space Network. Using a dual-frequency observing mode, we were able to observe Swift J1818.0−1607 simultaneously at either 2.2 and 8.4 GHz or 8.4 and 32 GHz. Over the course of the campaign, we find that the flux density increases substantially and the spectrum changes from uncharacteristically steep (α < −2.2) to the essentially flat (α ≈ 0) spectrum typical of radio-emitting magnetars. In addition to the expected profile evolution on time-scales of days to months, we find that Swift J1818.0−1607 also exhibits mode switching where the pulse profile changes between two distinct shapes on time-scales of seconds to minutes. For two of the radio observations, we also had accompanying X-ray observations using the Neutron Star Interior Composition Explorer telescope that occurred on the same day. We find a near anti-alignment (0.40 phase cycles) between the peaks of the radio and X-ray pulse profiles, which is most likely explained by an intrinsic misalignment between the X-ray- and radio-emitting regions. [ABSTRACT FROM AUTHOR]

Aims. Assuming fast radio bursts (FRBs) are produced by matter travelling ultra-relativistically in a localised region of a smooth bundle of streamlines, we study the constraints applied by geometry to the morphology and polarisation of the burst in time and frequency independently of the intrinsic radiative process. Methods. We express the problem only in terms of the local properties of direction and curvature of a streamline. This allowed us to cast the general results to any desired geometry. By applying this framework to two geometries inspired by pulsar and magnetar magnetospheres, we namely illustrate the dipolar polar-cap region and a magnetic dipole with an additional toroidal component. Results. Geometry constrains bursts to occur within an envelope in the frequency versus time plane (dynamic spectrum). This envelope notably characterises spectral occupancy and frequency drifts (both burst-to-burst and within an individual burst). We illustrate how one can simulate bursts by specifying some basic properties of an intrinsic emission process. In particular we show that the typical properties of one-off bursts can be produced in polar-cap geometry by a star with a spin period > 1 s, while bursts from repeating sources are better accounted for with an additional strong toroidal component and a sub-second spin period. Conclusions. We propose that a relationship between burst morphologies and the properties of the source, such as its spin period and magnetospheric properties, can be established at least qualitatively based on geometrical considerations. Our results favour models where repeaters are younger and faster magnetars with highly twisted magnetospheres. [ABSTRACT FROM AUTHOR]

While the dipole magnetic field axis of neutron stars is usually postulated to cross the star's centre, it may be displaced from this location, as it has been recently indicated in the millisecond pulsar J0030+0451. Under these conditions, the electromagnetic rocket effect may be activated, where the magnetic field exerts a net force, accelerating the star. This post-natal kick mechanism relies on asymmetric electromagnetic radiation from an off-centre dipole and may be relevant to the high spatial velocities of pulsars |$\sim 10^{3}$|  km s−1. Here, we explore its impact in young pulsars associated with supernova remnants, and we compare the observational data on characteristic quantities, such as the braking index and proper motion, with results obtained from the rocket effect. Using a Markov Chain Monte Carlo analysis, we explore the required conditions, for the initial spin periods and the distance between the magnetic axis and the star's centre, so that the velocity kick due to the rocket effect approaches the present velocity. We find that the electromagnetic rocket effect can account for typical pulsar transverse velocities assuming an initial spin period of 3.8  |$\rm {ms}$| and a dipole field whose distance from the centre of the star is approximately 7  |$\rm {km}$|. We also explore the influence of the rocket effect on the braking index of a neutron star, and we find that for the sample studied this impact is minimal. Finally, we apply the rocket effect model on the pulsars J0030+0451 and J0538+2817, which are likely candidates for this mechanism. [ABSTRACT FROM AUTHOR]

We study the stability of Hall magnetohydrodynamic with strong magnetic fields in which Landau quantization of electrons is important. We find that the strong-field Hall modes can be destabilized by the dependence of the differential magnetic susceptibility on magnetic field strength. This hydrodynamic instability, thermodynamic in origin and stabilized by magnetic domain formation, is studied using linear perturbation theory. It is found to have typical growth time of order ≲103 yr, with the growth time decreasing as a function of wavelength of the perturbation. The instability is self-limiting, turning off following a period of local field growth by a few per cent of the initial value. Finite temperature is also shown to limit the instability, with sufficiently high temperatures eliminating it altogether. Alfvén waves can show similar unstable behaviour on shorter time-scales. We find that Ohmic heating due to the large fields developed via the instability and magnetic domain formation is not large enough to account for observed magnetar surface temperatures. However, Ohmic heating is enhanced by the oscillatory differential magnetic susceptibility of Landau-quantized electrons, which could be important to magnetothermal simulations of neutron star crusts. [ABSTRACT FROM AUTHOR]

Neutron stars cool down during their lifetime through the combination of neutrino emission from the interior and photon cooling from the surface. Strongly magnetized neutron stars, called magnetars, are no exception, but the effect of their strong fields adds further complexities to the cooling theory. Besides other factors, modelling the outermost hundred meters (the envelope) plays a crucial role in predicting their surface temperatures. In this letter, we revisit the influence of envelopes on the cooling properties of neutron stars, with special focus on the critical effects of the magnetic field. We explore how our understanding of the relation between the internal and surface temperatures has evolved over the past two decades, and how different assumptions about the neutron star envelope and field topology lead to radically different conclusions on the surface temperature and its cooling with age. In particular, we find that relatively old magnetars with core-threading magnetic fields are actually much cooler than a rotation-powered pulsar of the same age. This is at variance with what is typically observed in crustal-confined models. Our results have important implications for the estimates of the X-ray luminosities of aged magnetars, and the subsequent population study of the different neutron star classes. [ABSTRACT FROM AUTHOR]

Fast radio bursts (FRBs) are millisecond-time-scale radio transients, the origins of which are predominantly extragalactic and likely involve highly magnetized compact objects. FRBs undergo multipath propagation, or scattering, from electron density fluctuations on sub-parsec scales in ionized gas along the line of sight. Scattering observations have located plasma structures within FRB host galaxies, probed Galactic and extragalactic turbulence, and constrained FRB redshifts. Scattering also inhibits FRB detection and biases the observed FRB population. We report the detection of scattering times from the repeating FRB 20190520B that vary by up to a factor of 2 or more on minutes to days-long time-scales. In one notable case, the scattering time varied from 7.9 ± 0.4 ms to less than 3.1 ms (⁠|$95{{\ \rm per\ cent}}$|  confidence) over 2.9 min at 1.45 GHz. The scattering times appear to be uncorrelated between bursts or with dispersion and rotation measure variations. Scattering variations are attributable to dynamic, inhomogeneous plasma in the circumsource medium, and analogous variations have been observed from the Crab pulsar. Under such circumstances, the frequency dependence of scattering can deviate from the typical power law used to measure scattering. Similar variations may therefore be detectable from other FRBs, even those with inconspicuous scattering, providing a unique probe of small-scale processes within FRB environments. [ABSTRACT FROM AUTHOR]

Magnetars exhibit a variety of transient high-energy phenomena in the form of bursts, outbursts, and giant flares. It is a common belief that these events originate in the sudden release of magnetic energy due to the rearrangement of a twisted magnetic field. We present global models of a 3D force-free (FF) non-linear twisted magnetar magnetosphere. We solve the FF equations following the Grad–Rubin approach in a compactified spherical domain. Appropriate boundary conditions are imposed at the surface of the star for the current distribution and the magnetic field. Our implementation is tested by reproducing various known analytical as well as axisymmetric numerical results. We then proceed to study general 3D models with non-axisymmetric current distributions, such as fields with localized twists that resemble hotspots at the surface of the star, and we examine characteristic quantities such as energy, helicity, and twist. Finally, we discuss implications on the available energy budget, the surface temperature, and the diffusion time-scale, which can be associated with observations. [ABSTRACT FROM AUTHOR]

The elasticity of neutron star crust is important for adequate interpretation of observations. To describe elastic properties one should rely on theoretical models. The most widely used is Coulomb crystal model (system of point-like charges on neutralizing uniform background), in some works it is corrected for electron screening. These models neglect finite size of nuclei. This approximation is well justified except for the innermost crustal layers, where nuclei size becomes comparable with the inter-nuclear spacing. Still, even in those dense layers it seems reasonable to apply the Coulomb crystal result, if one assumes that nuclei are spherically symmetric: Coulomb interaction between them should be the same as interaction between point-like charges. This argument is indeed correct; however, as we point here, shear of crustal lattice generates (microscopic) quadrupole electrostatic potential in a vicinity of lattice cites, which induces deformation on the nuclei. We analyse this problem analytically within compressible liquid drop model. In particular, for ground state crust composition the effective shear modulus is reduced for a factor of |$1-u^{5/3}/(2+3\, u-4\, u^{1/3})$|⁠ , where u is the ratio of the nuclei volume to the volume of the cell. This result is universal, i.e. it does not depend on the applied nucleon interaction model within applied approach. For the innermost layers of inner crust u ∼ 0.2 leading to reduction of the shear modulus by |$\sim 25{{\ \rm per\ cent}}$|⁠ , which can be important for correct interpretation of quasi-periodic oscillations in the tails of magnetar flares. [ABSTRACT FROM AUTHOR]

The long-term evolution of the internal, strong magnetic fields of neutron stars needs a specific numerical modelling. The diversity of the observed phenomenology of neutron stars indicates that their magnetic topology is rather complex and 3D simulations are required, for example, to explain the observed bursting mechanisms and the creation of surface hotspots. We present MATINS , a new 3D numerical code for magnetothermal evolution in neutron stars, based on a finite-volume scheme that employs the cubed-sphere system of coordinates. In this first work, we focus on the crustal magnetic evolution, with the inclusion of realistic calculations for the neutron star structure, composition, and electrical conductivity assuming a simple temperature evolution profile. MATINS follows the evolution of strong fields (1014 − 1015 Gauss) with complex non-axisymmetric topologies and dominant Hall-drift terms, and it is suitable for handling sharp current sheets. After introducing the technical description of our approach and some tests, we present long-term simulations of the non-linear field evolution in realistic neutron star crusts. The results show how the non-axisymmetric Hall cascade redistributes the energy over different spatial scales. Following the exploration of different initial topologies, we conclude that during a few tens of kyr, an equipartition of energy between the poloidal and toroidal components happens at small-scales. However, the magnetic field keeps a strong memory of the initial large scales, which are much harder to be restructured or created. This indicates that large-scale configuration attained during the neutron star formation is crucial to determine the field topology at any evolution stage. [ABSTRACT FROM AUTHOR]

We present LOFAR imaging observations from the April/May 2020 active episode of magnetar SGR 1935 + 2154. We place the earliest radio limits on persistent emission following the low-luminosity fast radio burst FRB 200428 from the magnetar. We also perform an image-plane search for transient emission and find no radio flares during our observations. We examine post-FRB radio upper limits in the literature and find that all are consistent with the multiwavelength afterglow predicted by the synchrotron maser shock model interpretation of FRB 200428. However, early optical observations appear to rule out the simple versions of the afterglow model with constant-density circumburst media. We show that these constraints may be mitigated by adapting the model for a wind-like environment, but only for a limited parameter range. In addition, we suggest that late-time non-thermal particle acceleration occurs within the afterglow model when the shock is no longer relativistic, which may prove vital for detecting afterglows from other Galactic FRBs. We also discuss future observing strategies for verifying either magnetospheric or maser shock FRB models via rapid radio observations of Galactic magnetars and nearby FRBs. [ABSTRACT FROM AUTHOR]

The short gamma ray bursts (GRBs) are the aftermath of the merger of binary compact objects (neutron star–neutron star or neutron star–black hole systems). With the simultaneous detection of gravitational wave (GW) signal from GW 170817 and GRB 170817A, the much-hypothesized connection between GWs and short GRBs has been proved beyond doubt. The resultant product of the merger could be a millisecond magnetar or a black hole depending on the binary masses and their equation of state. In the case of a magnetar central engine, fraction of the rotational energy deposited to the emerging ejecta produces late-time synchrotron radio emission from the interaction with the ambient medium. In this paper, we present an analysis of a sample of short GRBs located at a redshift of z ≤ 0.16 , which were observed at the late-time to search for the emission from merger ejecta. Our sample consists of seven short GRBs, which have radio upper limits available from very large array and Australian telescope compact array observations. We generate the model light curves using the standard magnetar model incorporating the relativistic correction. Using the model light curves and upper limits we constrain the number density of the ambient medium to be 10 - 5 – 10 - 3 cm - 3 for rotational energy of the magnetar E rot ∼ 5 × 10 51 erg. Variation in ejecta mass does not play a significant role in constraining the number density. [ABSTRACT FROM AUTHOR]

In the seconds following their formation in core-collapse supernovae, 'proto'-magnetars drive neutrino-heated magnetocentrifugal winds. Using a suite of two-dimensional axisymmetric magnetohydrodynamic simulations, we show that relatively slowly rotating magnetars with initial spin periods of P ⋆0 = 50–500 ms spin down rapidly during the neutrino Kelvin–Helmholtz cooling epoch. These initial spin periods are representative of those inferred for normal Galactic pulsars, and much slower than those invoked for gamma-ray bursts and superluminous supernovae. Since the flow is non-relativistic at early times, and because the Alfvén radius is much larger than the proto-magnetar radius, spin-down is millions of times more efficient than the typically used dipole formula. Quasi-periodic plasmoid ejections from the closed zone enhance spin-down. For polar magnetic field strengths B 0 ≳ 5 × 1014 G, the spin-down time-scale can be shorter than the Kelvin–Helmholtz time-scale. For B 0 ≳ 1015 G, it is of the order of seconds in early phases. We compute the spin evolution for cooling proto-magnetars as a function of B 0, P ⋆0, and mass (M). Proto-magnetars born with B 0 greater than |$\simeq 1.3\times 10^{15}\, {\rm \, G}\, (P_{\star 0}/{400\, \rm \, ms})^{-1.4}(M/1.4\, {\rm M}_\odot)^{2.2}$| spin down to periods >1 s in just the first few seconds of evolution, well before the end of the cooling epoch and the onset of classic dipole spin-down. Spin-down is more efficient for lower M and for larger P ⋆0. We discuss the implications for observed magnetars, including the discrepancy between their characteristic ages and supernova remnant ages. Finally, we speculate on the origin of 1E 161348−5055 in the remnant RCW 103, and the potential for other ultra-slowly rotating magnetars. [ABSTRACT FROM AUTHOR]

It is well established that magnetars are neutron stars with extreme magnetic fields and young ages, but the evolutionary pathways to their creation are still uncertain. Since most massive stars are in binaries, if magnetars are a frequent result of core-collapse supernovae, some fractions are expected to have a bound companion at the time of observation. In this paper, we utilize literature constraints, including deep Hubble Space Telescope imaging, to search for bound stellar companions to magnetars. The magnitude and colour measurements are interpreted in the context of binary population synthesis predictions. We find two candidates for stellar companions associated with CXOU J171405.7–381031 and SGR 0755–2933, based on their J–H colours and H -band absolute magnitudes. Overall, the proportion of the Galactic magnetar population with a plausibly stellar near-infrared (NIR) counterpart candidate, based on their magnitudes and colours, is between 5 and 10 per cent. This is consistent with a population synthesis prediction of 5 per cent, for the fraction of core-collapse neutron stars arising from primaries that remain bound to their companion after the supernova. These results are therefore consistent with magnetars being drawn in an unbiased way from the natal core-collapse neutron star population, but some contribution from alternative progenitor channels cannot be ruled out. [ABSTRACT FROM AUTHOR]

The recently discovered radio pulsar GLEAM-X J162759.5–523504.3 with an extremely long spin period was reported to have a radio luminosity that exceeds by orders of magnitude the spin-down power of the pulsar. In this Letter, we rigorously calculate the radio luminosity of the source taking into account the dependence of the opening angle of the pulsar-emission cone, first on the spin period alone and then on both the spin parameters and the observing frequency. We also revise the value of the spin-down power reported previously. Our analysis is based on the description of the spectral data in terms of two power-law indices as well as a single power-law index. Even if the pulsar's opening angle is treated as a frequency-independent parameter in line with the usual assumption, the period dependence of this parameter implies relatively small opening angles, and therefore, radio luminosities well below the spin-down power. Although we estimate higher radio luminosities in the physically more plausible case of a frequency-dependent opening angle, the spin-down power is again not exceeded by the highest possible radio luminosity. The radio efficiency of GLEAM-X J162759.5–523504.3 can therefore not be used in favour of a magnetar hypothesis. [ABSTRACT FROM AUTHOR]

We report on the the long term monitoring campaign of the seemingly youngest magnetar Swift J1818.0–1607 at radio and X-ray wavelengths over a span of one year. We obtained a coherent timing solution for the magnetar over the same time span. The frequency derivative of the magnetar shows systematic variation with the values oscillating about a mean value of −1.37 × 10−11 Hz s−1. The magnitude of the variation in the frequency derivative reduces with time before converging on the mean value. This corresponds to a characteristic age of ∼ 860 years, 2–4 times more than previously estimated. We were able to identify four states in the spin-frequency derivative that were quantified by the amount of modulation about the mean value and the transition between these states seem to be correlated with the change in the radio emission of the magnetar while no correlation is seen in the average radio profile variability on a shorter time-scale (days). The 0.5–12 keV X-ray flux shows a monotonic decrease that can be attributed to thermal emission from a hot spot on the surface of the neutron star that is reducing in size. Such decrease is consistent with what is seen in other magnetars. The potential correlation between the radio emission mode and the behaviour of the spin-down rate hints to a global change in the magnetopshere of the magnetar akin to the correlation seen in a subset of mode-changing radio pulsars and suggests a physical link between the two sub-populations. [ABSTRACT FROM AUTHOR]

We model neutron stars as magnetized hybrid stars with an abrupt hadron–quark phase transition in their cores, taking into account current constraints from nuclear experiments and multimessenger observations. We include magnetic field effects considering the Landau level quantization of charged particles and the anomalous magnetic moment of neutral particles. We construct the magnetized hybrid equation of state, and we compute the particle population, the matter magnetization and the transverse and parallel pressure components. We integrate the stable stellar models, considering the dynamical stability for rapid or slow hadron–quark phase conversion. Finally, we calculate the frequencies and damping times of the fundamental and g non-radial oscillation modes. The latter, a key mode to learn about phase transitions in compact objects, is only obtained for stars with slow conversions. For low magnetic fields, we find that one of the objects of the GW170817 binary system might be a hybrid star belonging to the slow extended stability branch. For magnetars, we find that a stronger magnetic field always softens the hadronic equation of state. Besides, only for some parameter combinations a stronger magnetic field implies a higher hybrid star maximum mass. Contrary to previous results, the incorporation of anomalous magnetic moment does not affect the studied astrophysical quantities. We discuss possible imprints of the microphysics of the equation of state that could be tested observationally in the future, and that might help infer the nature of dense matter and hybrid stars. [ABSTRACT FROM AUTHOR]

Extremely bright coherent radio bursts with millisecond duration, reminiscent of cosmological fast radio bursts, were codetected with anomalously-hard X-ray bursts from a Galactic magnetar SGR 1935 + 2154. We investigate the possibility that the event was triggered by the magnetic energy injection inside the magnetosphere, thereby producing magnetically-trapped fireball (FB) and relativistic outflows simultaneously. The thermal component of the X-ray burst is consistent with a trapped FB with an average temperature of ∼200–300 keV and size of ∼105 cm. Meanwhile, the non-thermal component of the X-ray burst and the coherent radio burst may arise from relativistic outflows. We calculate the dynamical evolution of the outflow, launched with an energy budget of 1039–1040 erg comparable to that for the trapped FB, for different initial baryon load η and magnetization σ0. If hard X-ray and radio bursts are both produced by the energy dissipation of the outflow, the outflow properties are constrained by combining the conditions for photon escape and the intrinsic timing offset ≲ 10 ms among radio and X-ray burst spikes. We show that the hard X-ray burst must be generated at r X ≳ 108 cm from the magnetar, irrespective of the emission mechanism. Moreover, we find that the outflow quickly accelerates up to a Lorentz factor of 102 ≲ Γ ≲ 103 by the time it reaches the edge of the magnetosphere and the dissipation occurs at 1012 cm ≲ r radio, X ≲ 1014 cm. Our results imply either extremely-clean (η ≳ 104) or highly-magnetized (σ0 ≳ 103) outflows, which might be consistent with the rarity of the phenomenon. [ABSTRACT FROM AUTHOR]

We investigate the magneto–elastic equilibrium of a neutron star crust and magnetic energy stored by the elastic force. The solenoidal motion driven by the Lorentz force can be controlled by the magnetic elastic force, so that conditions for the magnetic field strength and geometry are less restrictive. For equilibrium models, the minor solenoidal part of the magnetic force is balanced by a weak elastic force because the irrotational part is balanced by the dominant gravity and pressure forces. Therefore, a strong magnetic field may be confined in the interior, regardless of poloidal or toroidal components. We numerically calculated axially symmetric models with the maximum shear–strain and found that a magnetic energy >1046 erg can be stored in the crust, even for a normal surface dipole-field-strength (<1013 G). The magnetic energy much exceeds the elastic energy (1044−1045 erg). The shear–stress spatial distribution revealed that the elastic structure is likely to break down near the surface. In particular, the critical position is highly localized at a depth less than 100 m from the surface. [ABSTRACT FROM AUTHOR]

The origin and fate of magnetars [young, extremely magnetized neutron stars (NSs)] remains unsolved. Probing their evolution is therefore crucial for investigating possible links to other species of isolated NSs, such as the X-ray dim NSs (XDINSs) and rotating radio transients (RRATs). Here, we investigate the spin evolution of magnetars. Two avenues of evolution are considered: one with exponentially decaying B-fields, the other with sub- and superexponential decay. Using Monte Carlo methods, we synthesize magnetar populations using different input distributions and physical parameters, such as for the initial spin period, its time derivative, and the B-field decay time-scale. Additionally, we introduce a fade-away procedure that can account for the fading of old magnetars, and we briefly discuss the effect of alignment of the B-field and spin axes. Imposing the Galactic core-collapse supernova rate of ∼20 kyr−1 as a strict upper limit on the magnetar birthrate and comparing the synthetic populations to the observed one using both manual and automatic optimization algorithms for our input parameter study, we find that the B-field must decay exponentially or superexponentially with a characteristic decay time-scale of 0.5−10 kyr (with a best value of ∼4 kyr). In addition, the initial spin period must be less than 2 s. If these constraints are kept, we conclude that there are multiple choices of input physics that can reproduce the observed magnetar population reasonably well. We also conclude that magnetars may well be evolutionary linked to the population of XDINSs, whereas they are in general unlikely to evolve into RRATs. [ABSTRACT FROM AUTHOR]

A key tool astronomers have to investigate the nature of extragalactic transients is their position on their host galaxies. Galactocentric offsets, enclosed fluxes, and the fraction of light statistic are widely used at different wavelengths to help infer the nature of transient progenitors. Motivated by the proposed link between magnetars and fast radio bursts (FRBs), we create a face-on image of the Milky Way using best estimates of its size, structure, and colour. We place Galactic magnetars, pulsars, low-mass, and high-mass X-ray binaries on this image, using the available distance information. Galactocentric offsets, enclosed fluxes, and fraction of light distributions for these systems are compared to extragalactic transient samples. We find that FRBs follow the distributions for Galactic neutron stars closest, with 24 (75 per cent) of the Anderson–Darling tests we perform having a p -value greater than 0.05. This suggests that FRBs are located on their hosts in a manner consistent with Galactic neutron stars on the Milky Way's light, although we cannot determine which specific neutron star population is the best match. The Galactic distributions are consistent with other extragalactic transients much less often across the range of comparisons made, with type Ia SNe in second place, at only 33 per cent of tests exceeding 0.05. Overall, our results provide further support for FRB models invoking isolated young neutron stars, or binaries containing a neutron star. [ABSTRACT FROM AUTHOR]

In some short gamma-ray bursts, precursor flares occurring ∼ seconds prior to the main episode have been observed. These flares may then be associated with the last few cycles of the inspiral when the orbital frequency is a few hundred Hz. During these final cycles, tidal forces can resonantly excite quasi-normal modes in the inspiralling stars, leading to a rapid increase in their amplitude. It has been shown that these modes can exert sufficiently strong strains on to the neutron star crust to instigate yieldings. Due to the typical frequencies of g - modes being ∼100 Hz, their resonances with the orbital frequency match the precursor timings and warrant further investigation. Adopting realistic equations of state and solving the general-relativistic pulsation equations, we study g -mode resonances in coalescing quasi-circular binaries, where we consider various stellar rotation rates, degrees of stratification, and magnetic field structures. We show that for some combination of stellar parameters, the resonantly excited g 1 and g 2 modes may lead to crustal failure and trigger precursor flares. [ABSTRACT FROM AUTHOR]

The mechanism by which the radiation received from obliquely rotating neutron stars is generated remains an open question half a century after the discovery of pulsars. In contrast, considerable progress has recently been made in determining the structure of the magnetosphere that surrounds these objects: numerical computations based on the force-free, magnetohydrodynamic, and particle-in-cell formalisms have now firmly established that the magnetosphere of an oblique rotator entails a current sheet outside its light cylinder whose rotating distribution pattern moves with linear speeds exceeding the speed of light in vacuum. However, the role played by the superluminal motion of this current sheet in generating the multiwavelength, focused pulses of radiation that we receive from neutron stars is unknown. Here, we insert the description of the current sheet provided by the numerical simulations in the classical expression for the retarded potential and thereby calculate the radiation field generated by this source in the time domain. We find a radiation consisting of highly focused pulses whose salient features (brightness temperature, polarization, spectrum and profile with microstructure and with a phase lag between the radio and gamma-ray peaks) are strikingly similar to those of the emission received from pulsars. In addition, the flux density of this radiation diminishes with the distance D from the star as D −3/2 (rather than D −2) in certain latitudinal directions: a result that suggests that the high energetic requirements normally attributed to magnetars and the sources of fast radio bursts and gamma-ray bursts could be artefacts of the assumption that the radiation fields of all sources necessarily decay as predicted by the inverse-square law. [ABSTRACT FROM AUTHOR]

We examine the equilibrium of a magnetized neutron star crust. We calculate axially symmetric models in which an elastic force balances solenoidal motion driven by a Lorentz force. A large variety of equilibrium models are allowed by incorporating the elastic shear deformation; in addition, toroidal-magnetic-field-dominated models are available. These results remarkably differ from those in barotropic fluid stars. We demonstrate some models wherein the magnetic energy exceeds the elastic energy. The excess comes from the fact that a large amount of magnetic energy is associated with the irrotational part of the magnetic force, which is balanced with gravity and pressure. It is sufficient for equilibrium models that the minor solenoidal part is balanced by a weak elastic force. We find that the elasticity in the crust plays an important role on the magnetic field confinement. Further, we present the spatial distribution of the shear stress at the elastic limit, by which the crust-fracture location can be identified. The result has useful implications for realistic crust-quake models. [ABSTRACT FROM AUTHOR]

Magnetic field evolution in neutron-star crusts is driven by the Hall effect and Ohmic dissipation, for as long as the crust is sufficiently strong to absorb Maxwell stresses exerted by the field and thus makes the momentum equation redundant. For the strongest neutron-star fields, however, stresses build to the point of crustal failure, at which point the standard evolution equations are no longer valid. Here, we study the evolution of the magnetic field of the crust up to and beyond crustal failure, whence the crust begins to flow plastically. We perform global axisymmetric evolutions, exploring different types of failures affecting a limited region of the crust. We find that a plastic flow does not simply suppress the Hall effect even in the regime of a low plastic viscosity, but it rather leads to non-trivial evolution – in some cases even overreacting and enhancing the impact of the Hall effect. Its impact is more pronounced in the toroidal field, with the differences on the poloidal field being less substantial. We argue that both the nature of magnetar bursts and their spin-down evolution will be affected by plastic flow, so that observations of these phenomena may help us to constrain the way the crust fails. [ABSTRACT FROM AUTHOR]

We study magnetohydrodynamic stability of neutron star core matter composed of neutrons, protons, and leptons threaded by a magnetar-strength magnetic field 1014–1017 G, where quantum electrodynamical effects and Landau quantization of fermions are important. Stability is determined using the Friedman–Schutz formalism for the canonical energy of fluid perturbations, which we calculate for a magnetizable fluid with H ≠ B. Using this and the Euler–Heisenberg–Fermi–Dirac Lagrangian for a strongly magnetized fluid of Landau-quantized charged fermions, we calculate the local stability criteria for a neutron star core with a spherical axisymmetric geometry threaded by a toroidal field, accounting for magnetic and composition gradient buoyancy. We find that, for sufficiently strong fields B ≳ 1015 G, the magnetized fluid is unstable to a magnetosonic-type instability with growth times of the order of 10−3 s. The instability is triggered by sharp changes in the second-order field derivative of the Euler–Heisenberg–Fermi–Dirac Lagrangian that occur where additional Landau levels start being populated. These sharp changes are divergent at zero temperature, but are finite for non-zero temperature, so realistic neutron star core temperatures 5 × 107 K < T < 5 × 108 K are used. We conjecture that this mechanism could promote the formation of magnetic domains as predicted by Blandford and Hernquist and Suh and Mathews. [ABSTRACT FROM AUTHOR]