-
Constraining light QCD axions with isolated neutron star cooling
Authors:
Antonio Gómez-Bañón,
Kai Bartnick,
Konstantin Springmann,
José A. Pons
Abstract:
The existence of lighter-than-expected QCD axions can lead to a new ground state of matter, where the sourced axion field reduces the nucleon effective mass. The presence of the axion field has structural consequences, in particular, it results in a thinner (or even prevents its existence) heat-blanketing envelope, significantly altering the cooling patterns of neutron stars. We exploit the anomal…
▽ More
The existence of lighter-than-expected QCD axions can lead to a new ground state of matter, where the sourced axion field reduces the nucleon effective mass. The presence of the axion field has structural consequences, in particular, it results in a thinner (or even prevents its existence) heat-blanketing envelope, significantly altering the cooling patterns of neutron stars. We exploit the anomalous cooling behavior to constrain previously uncharted regions of the axion parameter space by comparing model predictions with existing data from isolated neutron stars. Notably, this analysis does not require the lighter-than-expected QCD axion to be the dark matter candidate.
△ Less
Submitted 14 August, 2024;
originally announced August 2024.
-
On the Origin of Magnetar Fields: Chiral Magnetic Instability in Neutron Star Crusts
Authors:
Clara Dehman,
José A. Pons
Abstract:
We investigate the chiral magnetic instability in the crust of a neutron star as a potential mechanism for amplifying magnetic fields. This instability may become active when small deviations from chemical equilibrium are sustained over decades, driven by the star's gradual spin-down or residual heat loss. Our findings suggest that this mechanism can produce strong, large-scale magnetic fields con…
▽ More
We investigate the chiral magnetic instability in the crust of a neutron star as a potential mechanism for amplifying magnetic fields. This instability may become active when small deviations from chemical equilibrium are sustained over decades, driven by the star's gradual spin-down or residual heat loss. Our findings suggest that this mechanism can produce strong, large-scale magnetic fields consistent with models that align with observational data. Additionally, this instability naturally generates magnetic helicity in the star's crust, which is crucial for forming and maintaining strong dipolar toroidal fields, often invoked to explain magnetar observational phenomena. Our results offer a microphysically-based alternative to classical hydrodynamical dynamos for the origin of magnetar magnetic fields, addressing a long-standing debate in the field.
△ Less
Submitted 9 August, 2024;
originally announced August 2024.
-
Unveiling the optimization process of Physics Informed Neural Networks: How accurate and competitive can PINNs be?
Authors:
Jorge F. Urbán,
Petros Stefanou,
José A. Pons
Abstract:
This study investigates the potential accuracy boundaries of physics-informed neural networks, contrasting their approach with previous similar works and traditional numerical methods. We find that selecting improved optimization algorithms significantly enhances the accuracy of the results. Simple modifications to the loss function may also improve precision, offering an additional avenue for enh…
▽ More
This study investigates the potential accuracy boundaries of physics-informed neural networks, contrasting their approach with previous similar works and traditional numerical methods. We find that selecting improved optimization algorithms significantly enhances the accuracy of the results. Simple modifications to the loss function may also improve precision, offering an additional avenue for enhancement. Despite optimization algorithms having a greater impact on convergence than adjustments to the loss function, practical considerations often favor tweaking the latter due to ease of implementation. On a global scale, the integration of an enhanced optimizer and a marginally adjusted loss function enables a reduction in the loss function by several orders of magnitude across diverse physical problems. Consequently, our results obtained using compact networks (typically comprising 2 or 3 layers of 20-30 neurons) achieve accuracies comparable to finite difference schemes employing thousands of grid points. This study encourages the continued advancement of PINNs and associated optimization techniques for broader applications across various fields.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Constraints on the dense matter equation of state from young and cold isolated neutron stars
Authors:
Alessio Marino,
Clara Dehman,
Konstantinos Kovlakas,
Nanda Rea,
Jose A. Pons,
D. Viganò
Abstract:
Neutron stars are the dense and highly magnetic relics of supernova explosions of massive stars. The quest to constrain the Equation of State (EoS) of ultra-dense matter and thereby probe the behavior of matter inside neutron stars, is one of the core goals of modern physics and astrophysics. A promising method involves investigating the long-term cooling of neutron stars, and comparing theoretica…
▽ More
Neutron stars are the dense and highly magnetic relics of supernova explosions of massive stars. The quest to constrain the Equation of State (EoS) of ultra-dense matter and thereby probe the behavior of matter inside neutron stars, is one of the core goals of modern physics and astrophysics. A promising method involves investigating the long-term cooling of neutron stars, and comparing theoretical predictions with various sources at different ages. However, limited observational data, and uncertainties in source ages and distances, have hindered this approach. In this work, re-analyzing XMM-Newton and Chandra data from dozens of thermally emitting isolated neutron stars, we have identified three sources with unexpectedly cold surface temperatures for their young ages. To investigate these anomalies, we conducted magneto-thermal simulations across diverse mass and magnetic fields, considering three different EoS. We found that the "minimal" cooling model, failed to explain the observations, regardless the mass and the magnetic field, as validated by a machine learning classification method. The existence of these young cold neutron stars suggests that any dense matter EoS must be compatible with a fast cooling process at least in certain mass ranges, eliminating a significant portion of current EoS options according to recent meta-modelling analysis.
△ Less
Submitted 8 April, 2024;
originally announced April 2024.
-
3D code for MAgneto-Thermal evolution in Isolated Neutron Stars, MATINS: thermal evolution and lightcurves
Authors:
Stefano Ascenzi,
Daniele Viganò,
Clara Dehman,
José A. Pons,
Nanda Rea,
Rosalba Perna
Abstract:
The thermal evolution of isolated neutron stars is a key element in unraveling 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 magn…
▽ More
The thermal evolution of isolated neutron stars is a key element in unraveling 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 which 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.
△ Less
Submitted 2 September, 2024; v1 submitted 28 January, 2024;
originally announced January 2024.
-
Axion sourcing in dense stellar matter via CP-violating couplings
Authors:
Filippo Anzuini,
Antonio Gómez-Bañón,
José A. Pons,
Andrew Melatos,
Paul D. Lasky
Abstract:
Compact objects such as neutron stars and white dwarfs can source axion-like particles and QCD axions due to CP-violating axion-fermion couplings. The magnitude of the axion field depends on the stellar density and on the strength of the axion-fermion couplings. We show that even CP-violating couplings one order of magnitude smaller than existing constraints source extended axion field configurati…
▽ More
Compact objects such as neutron stars and white dwarfs can source axion-like particles and QCD axions due to CP-violating axion-fermion couplings. The magnitude of the axion field depends on the stellar density and on the strength of the axion-fermion couplings. We show that even CP-violating couplings one order of magnitude smaller than existing constraints source extended axion field configurations. For axion-like particles, the axion energy is comparable to the magnetic energy in neutron stars with inferred magnetic fields of the order of $10^{13}$ G, and exceeds by more than one order of magnitude the magnetic energy content of white dwarfs with inferred fields of the order of $10^{4}$ G. On the other hand, the energy stored in the QCD axion field is orders of magnitude lower due to the smallness of the predicted CP-violating couplings. It is shown that the sourced axion field can polarize the photons emitted from the stellar surface, and stimulate the production of photons with energies in the radio band.
△ Less
Submitted 22 November, 2023;
originally announced November 2023.
-
One- and two-argument equation of state parametrizations with continuous sound speed for neutron star simulations
Authors:
Gaël Servignat,
Philip J. Davis,
Jérôme Novak,
Micaela Oertel,
José A. Pons
Abstract:
We describe two fitting schemes that aim to represent the high-density part of realistic equations of state for numerical simulations such as neutron star oscillations. The low-density part of the equation of state is represented by an arbitrary polytropic crust, and we propose a generic procedure to stitch any desired crust to the high-density fit, which is performed on the internal energy, press…
▽ More
We describe two fitting schemes that aim to represent the high-density part of realistic equations of state for numerical simulations such as neutron star oscillations. The low-density part of the equation of state is represented by an arbitrary polytropic crust, and we propose a generic procedure to stitch any desired crust to the high-density fit, which is performed on the internal energy, pressure and sound speed of barotropic equations of state that describe cold neutron stars in $β$-equilibrium. An extension of the fitting schemes to equations of state with an additional compositional argument is proposed. In particular we develop a formalism that ensures the existence of a $β$-equilibrium at low densities. An additional feature of this low-density model is that it can be, in principle, applied to any parametrization. The performance of the fits is checked on mass, radius and tidal deformability as well as on the dynamical radial oscillation frequencies. To that end, we use a pseudospectral isolated neutron star evolution code based on a non-conservative form of the hydrodynamical equations. A comparison to existing parametrizations is proposed, as far as possible, and to published radial frequency values in the literature. The static and dynamic quantities are well reproduced by the fitting schemes. Our results suggest that, even though the radius is very sensitive to the choice of the crust, this choice has little influence on the oscillation frequencies of a neutron star.
△ Less
Submitted 17 May, 2024; v1 submitted 5 November, 2023;
originally announced November 2023.
-
Solving the Pulsar Equation using Physics-Informed Neural Networks
Authors:
Petros Stefanou,
Jorge F. Urbán,
José A. Pons
Abstract:
In this study, Physics-Informed Neural Networks (PINNs) are skilfully applied to explore a diverse range of pulsar magneto-spheric models, specifically focusing on axisymmetric cases. The study successfully reproduced various axisymmetric models found in the literature, including those with non-dipolar configurations, while effectively characterizing current sheet features. Energy losses in all st…
▽ More
In this study, Physics-Informed Neural Networks (PINNs) are skilfully applied to explore a diverse range of pulsar magneto-spheric models, specifically focusing on axisymmetric cases. The study successfully reproduced various axisymmetric models found in the literature, including those with non-dipolar configurations, while effectively characterizing current sheet features. Energy losses in all studied models were found to exhibit reasonable similarity, differing by no more than a factor of three from the classical dipole case. This research lays the groundwork for a reliable elliptic Partial Differential Equation solver tailored for astrophysical problems. Based on these findings, we foresee that the utilization of PINNs will become the most efficient approach in modelling three-dimensional magnetospheres. This methodology shows significant potential and facilitates an effortless generalization, contributing to the advancement of our understanding of pulsar magnetospheres.
△ Less
Submitted 10 October, 2023; v1 submitted 12 September, 2023;
originally announced September 2023.
-
3D evolution of neutron star magnetic-fields from a realistic core-collapse turbulent topology
Authors:
Clara Dehman,
Daniele Viganò,
Stefano Ascenzi,
Jose A. Pons,
Nanda Rea
Abstract:
We present the first 3D fully coupled magneto-thermal simulations of neutron stars (including the most realistic background structure and microphysical ingredients so far) applied to a very complex initial magnetic field topology in the crust, similar to what recently obtained by proto-neutron star dynamo simulations. In such configurations, most of the energy is stored in the toroidal field, whil…
▽ More
We present the first 3D fully coupled magneto-thermal simulations of neutron stars (including the most realistic background structure and microphysical ingredients so far) applied to a very complex initial magnetic field topology in the crust, similar to what recently obtained by proto-neutron star dynamo simulations. In such configurations, most of the energy is stored in the toroidal field, while the dipolar component is a few percent of the mean magnetic field. This initial feature is maintained during the long-term evolution (1e6 yr), since the Hall term favours a direct cascade (compensating for Ohmic dissipation) rather than a strong inverse cascade, for such an initial field topology. The surface dipolar component, responsible for the dominant electromagnetic spin-down torque, does not show any increase in time, when starting from this complex initial topology. This is at contrast with the timing properties of young pulsars and magnetars which point to higher values of the surface dipolar fields. A possibility is that the deep-seated magnetic field (currents in the core) is able to self-organize in large scales (during the collapse or in the early life of a neutron star). Alternatively, the dipolar field might be lower than is usually thought, with magnetosphere substantially contributing to the observed high spin-down, via e.g., strong winds or strong coronal magnetic loops, which can also provide a natural explanation to the tiny surface hotspots inferred from X-ray data.
△ Less
Submitted 10 May, 2023;
originally announced May 2023.
-
Modelling Force-Free Neutron Star Magnetospheres using Physics-Informed Neural Networks
Authors:
Jorge F. Urbán,
Petros Stefanou,
Clara Dehman,
José A. Pons
Abstract:
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…
▽ More
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 classical methods. These results pave the way for the future extension to 3D of this (or a similar) problem, where generalised boundary conditions are very costly to implement.
△ Less
Submitted 21 March, 2023;
originally announced March 2023.
-
How bright can old magnetars be? Assessing the impact of magnetized envelopes and field topology on neutron star cooling
Authors:
Clara Dehman,
José A. Pons,
Daniele Viganò,
Nanda Rea
Abstract:
Neutron stars cool down during their lifetime through the combination of neutrino emission from the interior and photon cooling from the surface. Strongly magnetised 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 rol…
▽ More
Neutron stars cool down during their lifetime through the combination of neutrino emission from the interior and photon cooling from the surface. Strongly magnetised 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.
△ Less
Submitted 5 January, 2023;
originally announced January 2023.
-
Magnetic dynamo caused by axions in neutron stars
Authors:
Filippo Anzuini,
José A. Pons,
Antonio Gómez-Bañón,
Paul D. Lasky,
Federico Bianchini,
Andrew Melatos
Abstract:
The coupling between axions and photons modifies Maxwell's equations, introducing a dynamo term in the magnetic induction equation. In neutron stars, for critical values of the axion decay constant and axion mass, the magnetic dynamo mechanism increases the total magnetic energy of the star. We show that this generates substantial internal heating due to enhanced dissipation of crustal electric cu…
▽ More
The coupling between axions and photons modifies Maxwell's equations, introducing a dynamo term in the magnetic induction equation. In neutron stars, for critical values of the axion decay constant and axion mass, the magnetic dynamo mechanism increases the total magnetic energy of the star. We show that this generates substantial internal heating due to enhanced dissipation of crustal electric currents. These mechanisms would lead magnetized neutron stars to increase their magnetic energy and thermal luminosity by several orders of magnitude, in contrast to observations of thermally-emitting neutron stars. To prevent the activation of the dynamo, bounds on the allowed axion parameter space can be derived.
△ Less
Submitted 19 November, 2022;
originally announced November 2022.
-
Modelling 3D Force-Free Neutron Star Magnetospheres
Authors:
Petros Stefanou,
Jose A. Pons,
Pablo Cerdá-Durán
Abstract:
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-Rub…
▽ More
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 localised 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 timescale, which can be associated with observations.
△ Less
Submitted 21 March, 2023; v1 submitted 16 November, 2022;
originally announced November 2022.
-
Constraining the nature of the 18-min periodic radio transient GLEAM-X J162759.5-523504.3 via multi-wavelength observations and magneto-thermal simulations
Authors:
N. Rea,
F. Coti Zelati,
C. Dehman,
N. Hurley-Walker,
D. De Martino,
A. Bahramian,
D. A. H. Buckley,
J. Brink,
A. Kawka,
J. A. Pons,
D. Vigano,
V. Graber,
M. Ronchi,
C. Pardo,
A. Borghese,
E. Parent
Abstract:
We observed the periodic radio transient GLEAM-X J162759.5-523504.3 (GLEAM-X J1627) using the Chandra X-ray Observatory for about 30-ks on January 22-23, 2022, simultaneously with radio observations from MWA, MeerKAT and ATCA. Its radio emission and 18-min periodicity led the source to be tentatively interpreted as an extreme magnetar or a peculiar highly magnetic white dwarf. The source was not d…
▽ More
We observed the periodic radio transient GLEAM-X J162759.5-523504.3 (GLEAM-X J1627) using the Chandra X-ray Observatory for about 30-ks on January 22-23, 2022, simultaneously with radio observations from MWA, MeerKAT and ATCA. Its radio emission and 18-min periodicity led the source to be tentatively interpreted as an extreme magnetar or a peculiar highly magnetic white dwarf. The source was not detected in the 0.3-8 keV energy range with a 3-sigma upper-limit on the count rate of 3x10^{-4} counts/s. No radio emission was detected during our X-ray observations either. Furthermore, we studied the field around GLEAM-X J1627 using archival ESO and DECam data, as well as recent SALT observations. Many sources are present close to the position of GLEAM-X J1627, but only two within the 2" radio position uncertainty. Depending on the assumed spectral distribution, the upper limits converted to an X-ray luminosity of L_{X}<6.5x10^{29} erg/s for a blackbody with temperature kT=0.3 keV, or L_{X}<9x10^{29} erg/s for a power-law with photon index Gamma = 2 (assuming a 1.3 kpc distance). Furthermore, we performed magneto-thermal simulations for neutron stars considering crust- and core-dominated field configurations. Based on our multi-band limits, we conclude that: i) in the magnetar scenario, the X-ray upper limits suggest that GLEAM-X J1627 should be older than ~1 Myr, unless it has a core-dominated magnetic field or has experienced fast-cooling; ii) in the white dwarf scenario, we can rule out most binary systems, a hot sub-dwarf and a hot magnetic isolated white dwarf (T>10.000 K), while a cold isolated white dwarf is still compatible with our limits.
△ Less
Submitted 4 October, 2022;
originally announced October 2022.
-
3D code for MAgneto-Thermal evolution in Isolated Neutron Stars, MATINS: The Magnetic Field Formalism
Authors:
Clara Dehman,
Daniele Viganò,
José A. Pons,
Nanda Rea
Abstract:
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 three-dimensional simulations are required, for example, to explain the observed bursting mechanisms and the creation of surface hotspots. We present MATINS, a…
▽ More
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 three-dimensional simulations are required, for example, to explain the observed bursting mechanisms and the creation of surface hotspots. We present MATINS, a new three dimensions numerical code for magneto-thermal 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 (1e14-1e15 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.
△ Less
Submitted 26 September, 2022;
originally announced September 2022.
-
Thermal luminosity degeneracy of magnetized neutron stars with and without hyperon cores
Authors:
F. Anzuini,
A. Melatos,
C. Dehman,
D. Viganò,
J. A. Pons
Abstract:
The dissipation of intense crustal electric currents produces high Joule heating rates in cooling neutron stars. Here it is shown that Joule heating can counterbalance fast cooling, making it difficult to infer the presence of hyperons (which accelerate cooling) from measurements of the observed thermal luminosity $L_γ$. Models with and without hyperon cores match $L_γ$ of young magnetars (with po…
▽ More
The dissipation of intense crustal electric currents produces high Joule heating rates in cooling neutron stars. Here it is shown that Joule heating can counterbalance fast cooling, making it difficult to infer the presence of hyperons (which accelerate cooling) from measurements of the observed thermal luminosity $L_γ$. Models with and without hyperon cores match $L_γ$ of young magnetars (with poloidal-dipolar field $B_{\textrm{dip}} \gtrsim 10^{14}$ G at the polar surface and $L_γ \gtrsim 10^{34}$ erg s$^{-1}$ at $t \lesssim 10^5$ yr) as well as mature, moderately magnetized stars (with $B_{\textrm{dip}} \lesssim 10^{14}$ G and $10^{31} \ \textrm{erg s}^{-1} \lesssim L_γ \lesssim 10^{32}$ erg s$^{-1}$ at $t \gtrsim 10^5$ yr). In magnetars, the crustal temperature is almost independent of hyperon direct Urca cooling in the core, regardless of whether the latter is suppressed or not by hyperon superfluidity. The thermal luminosities of light magnetars without hyperons and heavy magnetars with hyperons have $L_γ$ in the same range and are almost indistinguishable. Likewise, $L_γ$ data of neutron stars with $B_{\textrm{dip}} \lesssim 10^{14}$ G but with strong internal fields are not suitable to extract information about the equation of state as long as hyperons are superfluid, with maximum amplitude of the energy gaps of the order $\approx 1$ MeV.
△ Less
Submitted 29 May, 2022;
originally announced May 2022.
-
Fast cooling and internal heating in hyperon stars
Authors:
F. Anzuini,
A. Melatos,
C. Dehman,
D. Viganò,
J. A. Pons
Abstract:
Neutron star models with maximum mass close to $2 \ M_{\odot}$ reach high central densities, which may activate nucleonic and hyperon direct Urca neutrino emission. To alleviate the tension between fast theoretical cooling rates and thermal luminosity observations of moderately magnetized, isolated thermally-emitting stars (with $L_γ \gtrsim 10^{31}$ erg s$^{-1}$ at $t \gtrsim 10^{5.3}$ yr), some…
▽ More
Neutron star models with maximum mass close to $2 \ M_{\odot}$ reach high central densities, which may activate nucleonic and hyperon direct Urca neutrino emission. To alleviate the tension between fast theoretical cooling rates and thermal luminosity observations of moderately magnetized, isolated thermally-emitting stars (with $L_γ \gtrsim 10^{31}$ erg s$^{-1}$ at $t \gtrsim 10^{5.3}$ yr), some internal heating source is required. The power supplied by the internal heater is estimated for both a phenomenological source in the inner crust and Joule heating due to magnetic field decay, assuming different superfluidity models and compositions of the outer stellar envelope. It is found that a thermal power of $W(t) \approx 10^{34}$ erg s$^{-1}$ allows neutron star models to match observations of moderately magnetized, isolated stars with ages $t \gtrsim 10^{5.3}$ yr. The requisite $W(t)$ can be supplied by Joule heating due to crust-confined initial magnetic configurations with (i) mixed poloidal-toroidal fields, with surface strength $B_{\textrm{dip}} = 10^{13}$ G at the pole of the dipolar poloidal component and $\sim 90$ per cent of the magnetic energy stored in the toroidal component; and (ii) poloidal-only configurations with $B_{\textrm{dip}} = 10^{14}$ G.
△ Less
Submitted 26 October, 2021;
originally announced October 2021.
-
Magneto-thermal evolution of neutron stars with coupled Ohmic, Hall and ambipolar effects via accurate finite-volume simulations
Authors:
Daniele Viganò,
Alberto García-García,
José A. Pons,
Clara Dehman,
Vanessa Graber
Abstract:
Simulating the long-term evolution of temperature and magnetic fields in neutron stars is a major effort in astrophysics, having significant impact in several topics. A detailed evolutionary model requires, at the same time, the numerical solution of the heat diffusion equation, the use of appropriate numerical methods to control non-linear terms in the induction equation, and the local calculatio…
▽ More
Simulating the long-term evolution of temperature and magnetic fields in neutron stars is a major effort in astrophysics, having significant impact in several topics. A detailed evolutionary model requires, at the same time, the numerical solution of the heat diffusion equation, the use of appropriate numerical methods to control non-linear terms in the induction equation, and the local calculation of realistic microphysics coefficients. Here we present the latest extension of the magneto-thermal 2D code in which we have coupled the crustal evolution to the core evolution, including ambipolar diffusion. It has also gained in modularity, accuracy, and efficiency. We revise the most suitable numerical methods to accurately simulate magnetar-like magnetic fields, reproducing the Hall-driven magnetic discontinuities. From the point of view of computational performance, most of the load falls on the calculation of microphysics coefficients. To a lesser extent, the thermal evolution part is also computationally expensive because it requires large matrix inversions due to the use of an implicit method. We show two representative case studies: (i) a non-trivial multipolar configuration confined to the crust, displaying long-lived small-scale structures and discontinuities; and (ii) a preliminary study of ambipolar diffusion in normal matter. The latter acts on timescales that are too long to have relevant effects on the timescales of interest but sets the stage for future works where superfluid and superconductivity need to be included.
△ Less
Submitted 16 April, 2021;
originally announced April 2021.
-
Analyzing the Galactic pulsar distribution with machine learning
Authors:
Michele Ronchi,
Vanessa Graber,
Alberto Garcia-Garcia,
Jose A. Pons,
Nanda Rea
Abstract:
We explore the possibility of inferring the properties of the Galactic neutron star population through machine learning. In particular, in this paper we focus on their dynamical characteristics and show that an artificial neural network is able to estimate with high accuracy the parameters which control the current positions of a mock population of pulsars. For this purpose, we implement a simplif…
▽ More
We explore the possibility of inferring the properties of the Galactic neutron star population through machine learning. In particular, in this paper we focus on their dynamical characteristics and show that an artificial neural network is able to estimate with high accuracy the parameters which control the current positions of a mock population of pulsars. For this purpose, we implement a simplified population-synthesis framework (where selection biases are neglected at this stage) and concentrate on the natal kick-velocity distribution and the distribution of birth distances from the Galactic plane. By varying these and evolving the pulsar trajectories in time, we generate a series of simulations that are used to train and validate a suitably structured convolutional neural network. We demonstrate that our network is able to recover the parameters governing the kick-velocity and Galactic height distribution with a mean relative error of about $10^{-2}$. We discuss the limitations of our idealized approach and study a toy problem to introduce selection effects in a phenomenological way by incorporating the observed proper motions of 216 isolated pulsars. Our analysis highlights that increasing the sample of pulsars with accurate proper motion measurements by a factor of $\sim$10, one of the future breakthroughs of the Square Kilometer Array, we might succeed in constraining the birth spatial and kick-velocity distribution of the neutron stars in the Milky Way with high precision through machine learning.
△ Less
Submitted 23 June, 2021; v1 submitted 14 January, 2021;
originally announced January 2021.
-
On the rate of crustal failures in young magnetars
Authors:
Clara Dehman,
Daniele Viganò,
Nanda Rea,
Jose A. Pons,
Rosalba Perna,
Alberto Gracía-Gracía
Abstract:
The activity of magnetars is powered by their intense and dynamic magnetic fields and has been proposed as the trigger to extragalactic Fast Radio Bursts. Here we estimate the frequency of crustal failures in young magnetars, by computing the magnetic stresses in detailed magneto-thermal simulations including Hall drift and Ohmic dissipation. The initial internal topology at birth is poorly known…
▽ More
The activity of magnetars is powered by their intense and dynamic magnetic fields and has been proposed as the trigger to extragalactic Fast Radio Bursts. Here we estimate the frequency of crustal failures in young magnetars, by computing the magnetic stresses in detailed magneto-thermal simulations including Hall drift and Ohmic dissipation. The initial internal topology at birth is poorly known but is likely to be much more complex than a dipole. Thus, we explore a wide range of initial configurations, finding that the expected rate of crustal failures varies by orders of magnitude depending on the initial magnetic configuration. Our results show that this rate scales with the crustal magnetic energy, rather than with the often used surface value of the dipolar component related to the spin-down torque. The estimated frequency of crustal failures for a given dipolar component can vary by orders of magnitude for different initial conditions, depending on how much magnetic energy is distributed in the crustal non-dipolar components, likely dominant in newborn magnetars. The quantitative reliability of the expected event rate could be improved by a better treatment of the magnetic evolution in the core and the elastic/plastic crustal response, here not included. Regardless of that, our results are useful inputs in modelling the outburst rate of young Galactic magnetars, and their relation with the Fast Radio Bursts in our and other galaxies.
△ Less
Submitted 1 October, 2020;
originally announced October 2020.
-
The X-ray outburst of the Galactic Center magnetar over six years of Chandra observations
Authors:
N. Rea,
F. Coti Zelati,
D. Vigano',
A. Papitto,
F. Baganoff,
A. Borghese,
S. Campana,
P. Esposito,
D. Haggard,
G. L. Israel,
S. Mereghetti,
R. Mignani,
R. Perna,
J. A. Pons,
G. Ponti,
L. Stella,
D. F. Torres,
R. Turolla,
S. Zane
Abstract:
The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has…
▽ More
The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has opened interesting debates about the number, age and nature of pulsars expected in the Galactic center region. In this work, we present extensive X-ray monitoring of the outburst of SGR J1745-2900 using the Chandra X-ray Observatory, the only instrument with the spatial resolution to distinguish the magnetar from the supermassive black hole (2.4" angular distance). It was monitored from its outburst onset in April 2013 until August 2019, collecting more than fifty Chandra observations for a total of more than 2.3 Ms of data. Soon after the outburst onset, the magnetar emission settled onto a purely thermal emission state that cooled from a temperature of about 0.9 to 0.6 keV over 6 years. The pulsar timing properties showed at least two changes in the period derivative, increasing by a factor of about 4 during the outburst decay. We find that the long-term properties of this outburst challenge current models for the magnetar outbursts.
△ Less
Submitted 24 March, 2020; v1 submitted 16 March, 2020;
originally announced March 2020.
-
Nonaxisymmetric Hall instability: A key to understanding magnetars
Authors:
Konstantinos N. Gourgouliatos,
José A. Pons
Abstract:
It is generally accepted that the non-linear, dynamical evolution of magnetic fields in the interior of neutron stars plays a key role in the explanation of the observed phenomenology. Understanding the transfer of energy between toroidal and poloidal components, or between different scales, is of particular relevance. In this letter, we present the first 3D simulations of the Hall instability in…
▽ More
It is generally accepted that the non-linear, dynamical evolution of magnetic fields in the interior of neutron stars plays a key role in the explanation of the observed phenomenology. Understanding the transfer of energy between toroidal and poloidal components, or between different scales, is of particular relevance. In this letter, we present the first 3D simulations of the Hall instability in a neutron star crust, confirming its existence for typical magnetar conditions. We confront our results to estimates obtained by a linear perturbation analysis, which discards any interpretation as numerical instabilities and confirms its physical origin. Interestingly, the Hall instability creates locally strong magnetic structures that occasionally can make the crust yield to the magnetic stresses and generates coronal loops, similarly as solar coronal loops find their way out through the photosphere. This supports the viability of the mechanism, which has been proposed to explain magnetar outbursts.
△ Less
Submitted 10 January, 2020;
originally announced January 2020.
-
The long-term enhanced brightness of the magnetar 1E 1547.0-5408
Authors:
Francesco Coti Zelati,
Alice Borghese,
Nanda Rea,
Daniele Viganò,
Teruaki Enoto,
Paolo Esposito,
José A. Pons,
Sergio Campana,
Gian Luca Israel
Abstract:
We present the evolution of the X-ray emission properties of the magnetar 1E 1547.0-5408 since February 2004 over a time period covering three outbursts. We analyzed new and archival observations taken with the Swift, NuSTAR, Chandra and XMM-Newton X-ray satellites. The source has been observed at a relatively steady soft X-ray flux of $\approx$ 10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ (0.3-10 keV) over…
▽ More
We present the evolution of the X-ray emission properties of the magnetar 1E 1547.0-5408 since February 2004 over a time period covering three outbursts. We analyzed new and archival observations taken with the Swift, NuSTAR, Chandra and XMM-Newton X-ray satellites. The source has been observed at a relatively steady soft X-ray flux of $\approx$ 10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ (0.3-10 keV) over the last 9 years, which is about an order of magnitude fainter than the flux at the peak of the last outburst in 2009, but a factor of $\sim$ 30 larger than the level in 2006. The broad-band spectrum extracted from two recent NuSTAR observations in April 2016 and February 2019 showed a faint hard X-ray emission up to $\sim$ 70 keV. Its spectrum is adequately described by a flat power law component, and its flux is $\sim$ $7 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ (10-70 keV), that is a factor of $\sim$ 20 smaller than at the peak of the 2009 outburst. The hard X-ray spectral shape has flattened significantly in time, which is at variance with the overall cooling trend of the soft X-ray component. The pulse profile extracted from these NuSTAR pointings displays variability in shape and amplitude with energy (up to $\approx$ 25 keV). Our analysis shows that the flux of 1E 1547.0-5408 is not yet decaying to the 2006 level and that the source has been lingering in a stable, high-intensity state for several years. This might suggest that magnetars can hop among distinct persistent states that are probably connected to outburst episodes and that their persistent thermal emission can be almost entirely powered by the dissipation of currents in the corona.
△ Less
Submitted 19 December, 2019;
originally announced December 2019.
-
Magnetic, thermal and rotational evolution of isolated neutron stars
Authors:
José A. Pons,
Daniele Viganò
Abstract:
The strong magnetic field of neutron stars is intimately coupled to the observed temperature and spectral properties, as well as to the observed timing properties (distribution of spin periods and period derivatives). Thus, a proper theoretical and numerical study of the magnetic field evolution equations, supplemented with detailed calculations of microphysical properties (heat and electrical con…
▽ More
The strong magnetic field of neutron stars is intimately coupled to the observed temperature and spectral properties, as well as to the observed timing properties (distribution of spin periods and period derivatives). Thus, a proper theoretical and numerical study of the magnetic field evolution equations, supplemented with detailed calculations of microphysical properties (heat and electrical conductivity, neutrino emission rates) is crucial to understand how the strength and topology of the magnetic field vary as a function of age, which in turn is the key to decipher the physical processes behind the varied neutron star phenomenology. In this review, we go through the basic theory describing the magneto-thermal evolution models of neutron stars, focusing on numerical techniques, and providing a battery of benchmark tests to be used as a reference for present and future code developments. We summarize well-known results from axisymmetric cases, give a new look at the latest 3D advances, and present an overview of the expectations for the field in the coming years.
△ Less
Submitted 8 November, 2019;
originally announced November 2019.
-
Instability of twisted magnetar magnetospheres
Authors:
J. F. Mahlmann,
T. Akgün,
J. A. Pons,
M. A. Aloy,
P. Cerdá-Durán
Abstract:
We present three-dimensional force-free electrodynamics simulations of magnetar magnetospheres that demonstrate the instability of certain degenerate, high energy equilibrium solutions of the Grad-Shafranov equation. This result indicates the existence of an unstable branch of twisted magnetospheric solutions and allows to formulate an instability criterion. The rearrangement of magnetic field lin…
▽ More
We present three-dimensional force-free electrodynamics simulations of magnetar magnetospheres that demonstrate the instability of certain degenerate, high energy equilibrium solutions of the Grad-Shafranov equation. This result indicates the existence of an unstable branch of twisted magnetospheric solutions and allows to formulate an instability criterion. The rearrangement of magnetic field lines as a consequence of this instability triggers the dissipation of up to 30% of the magnetospheric energy on a thin layer above the magnetar surface. During this process, we predict an increase of the mechanical stresses onto the stellar crust, which can potentially result in a global mechanical failure of a significant fraction of it. We find that the estimated energy release and the emission properties are compatible with the observed giant flare events. The newly identified instability is a candidate for recurrent energy dissipation, which could explain part of the phenomenology observed in magnetars.
△ Less
Submitted 24 November, 2019; v1 submitted 31 July, 2019;
originally announced August 2019.
-
Triggering magnetar outbursts in 3D force-free simulations
Authors:
Federico Carrasco,
Daniele Viganò,
Carlos Palenzuela,
Jose A. Pons
Abstract:
In this letter, we present the first 3D force-free general relativity simulations of the magnetosphere dynamics related to the magnetar outburst/flare phenomenology. Starting from an initial dipole configuration, we adiabatically increase the helicity by twisting the footprints of a spot on the stellar surface and follow the succession of quasi-equilibrium states until a critical twist is reached.…
▽ More
In this letter, we present the first 3D force-free general relativity simulations of the magnetosphere dynamics related to the magnetar outburst/flare phenomenology. Starting from an initial dipole configuration, we adiabatically increase the helicity by twisting the footprints of a spot on the stellar surface and follow the succession of quasi-equilibrium states until a critical twist is reached. Twisting beyond that point triggers instabilities that results in the rapid expansion of magnetic field lines, followed by reconnection, as observed in previous axi-symmetric simulations. If the injection of magnetic helicity goes on, the process is recurrent, periodically releasing a similar amount of energy, of the order of a few % of the total magnetic energy. From our current distribution, we estimate the local temperature assuming that dissipation occurs mainly in the highly resistive outermost layer of the neutron star. We find that the temperature smoothly increases with injected twist, being larger for spots located in the tropical regions than in polar regions, and rather independent of their sizes. After the injection of helicity ceases, the magnetosphere relaxes to a new stable state, in which the persistent currents maintain the footprints area slightly hotter than before the onset of the instability.
△ Less
Submitted 25 January, 2019;
originally announced January 2019.
-
The multi-outburst activity of the magnetar in Westerlund I
Authors:
A. Borghese,
N. Rea,
R. Turolla,
J. A. Pons,
P. Esposito,
F. Coti Zelati,
V. Savchenko,
E. Bozzo,
R. Perna,
S. Zane,
S. Mereghetti,
S. Campana,
R. P. Mignani,
M. Bachetti,
G. Rodriguez,
F. Pintore,
A. Tiengo,
D. Gotz,
G. L. Israel,
L. Stella
Abstract:
After two major outbursts in 2006 and 2011, on 2017 May 16 the magnetar CXOU J164710.2-455216, hosted within the massive star cluster Westerlund I, emitted a short (20 ms) burst, which marked the onset of a new active phase. We started a long-term monitoring campaign with Swift (45 observations), Chandra (5 observations) and NuSTAR (4 observations) from the activation until 2018 April. During the…
▽ More
After two major outbursts in 2006 and 2011, on 2017 May 16 the magnetar CXOU J164710.2-455216, hosted within the massive star cluster Westerlund I, emitted a short (20 ms) burst, which marked the onset of a new active phase. We started a long-term monitoring campaign with Swift (45 observations), Chandra (5 observations) and NuSTAR (4 observations) from the activation until 2018 April. During the campaign, Swift BAT registered the occurrence of multiple bursts, accompanied by two other enhancements of the X-ray persistent flux. The long time span covered by our observations allowed us to study the spectral as well as the timing evolution of the source. After 11 months since the 2017 May outburst onset, the observed flux was about 15 times higher than its historical minimum level and a factor of 3 higher than the level reached after the 2006 outburst. This suggests that the crust has not fully relaxed to the quiescent level, or that the source quiescent level has changed following the multiple outburst activities in the past 10 years or so. This is another case of multiple outbursts from the same source on a yearly time scale, a somehow recently discovered behaviour in magnetars.
△ Less
Submitted 7 January, 2019;
originally announced January 2019.
-
A Simflowny-based high-performance 3D code for the generalized induction equation
Authors:
Daniele Viganò,
David Martínez-Gómez,
José A. Pons,
Carlos Palenzuela,
Federico Carrasco,
Borja Miñano,
Antoni Arbona,
Carles Bona,
Joan Massó
Abstract:
In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage…
▽ More
In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage and a consensus on the general dynamics at least for some simple initial data, only few attempts have been performed in 3D, due to the computational costs and the need for a proper numerical treatment of the intrinsic non-linearity of the equation. Here, we carefully analyze the general induction equation, studying its characteristic structure, and we present a new Cartesian 3D code, generated by the user-friendly, publicly available {\em Simflowny} platform. The code uses high-order numerical schemes for the time and spatial discretization, and relies on the highly-scalable {\em SAMRAI} architecture for the adaptive mesh refinement. We present the application of the code to several benchmark tests, showing the high order of convergence and accuracy achieved and the capabilities in terms of magnetic shock resolution and three-dimensionality. This paper paves the way for the applications to a realistic, 3D long-term evolution of neutron stars interior and, possibly, of other astrophysical sources.
△ Less
Submitted 11 December, 2018; v1 submitted 20 November, 2018;
originally announced November 2018.
-
Crust-magnetosphere coupling during magnetar evolution and implications for the surface temperature
Authors:
Taner Akgün,
Pablo Cerdá-Durán,
Juan Antonio Miralles,
José A. Pons
Abstract:
We study the coupling of the force-free magnetosphere to the long-term internal evolution of a magnetar. We allow the relation between the poloidal and toroidal stream functions - that characterizes the magnetosphere - to evolve freely without constraining its particular form. We find that, on time-scales of the order of kyr, the energy stored in the magnetosphere gradually increases, as the toroi…
▽ More
We study the coupling of the force-free magnetosphere to the long-term internal evolution of a magnetar. We allow the relation between the poloidal and toroidal stream functions - that characterizes the magnetosphere - to evolve freely without constraining its particular form. We find that, on time-scales of the order of kyr, the energy stored in the magnetosphere gradually increases, as the toroidal region grows and the field lines expand outwards. This continues until a critical point is reached beyond which force-free solutions for the magnetosphere can no longer be constructed, likely leading to some large-scale magnetospheric reorganization. The energy budget available for such events can be as high as several $10^{45}\,$erg for fields of $10^{14}\,$G. Subsequently, starting from the new initial conditions, the evolution proceeds in a similar manner. The time-scale to reach the critical point scales inversely with the magnetic field amplitude. Allowing currents to pass through the last few meters below the surface, where the magnetic diffusivity is orders of magnitude larger than in the crust, should give rise to a considerable amount of energy deposition through Joule heating. We estimate that the effective surface temperature could increase locally from $\sim 0.1\,$keV to $\sim 0.3 - 0.6\,$keV, in good agreement with observations. Similarly, the power input from the interior into the magnetosphere could be as high as $10^{35} - 10^{36}\,$erg/s, which is consistent with peak luminosities observed during magnetar outbursts. Therefore, a detailed treatment of currents flowing through the envelope may be needed to explain the thermal properties of magnetars.
△ Less
Submitted 9 October, 2018; v1 submitted 24 July, 2018;
originally announced July 2018.
-
Gazing at the ultra-slow magnetar in RCW 103 with NuSTAR and Swift
Authors:
A. Borghese,
F. Coti Zelati,
P. Esposito,
N. Rea,
A. De Luca,
M. Bachetti,
G. L. Israel,
R. Perna,
J. A. Pons
Abstract:
We report on a new NuSTAR observation and on the ongoing Swift XRT monitoring campaign of the peculiar source 1E 161348-5055, located at the centre of the supernova remnant RCW 103, which is recovering from its last outburst in June 2016. The X-ray spectrum at the epoch of the NuSTAR observation can be described by either two absorbed blackbodies ($kT_{BB_1}$ ~ 0.5 keV, $kT_{BB_2}$ ~ 1.2 keV) or a…
▽ More
We report on a new NuSTAR observation and on the ongoing Swift XRT monitoring campaign of the peculiar source 1E 161348-5055, located at the centre of the supernova remnant RCW 103, which is recovering from its last outburst in June 2016. The X-ray spectrum at the epoch of the NuSTAR observation can be described by either two absorbed blackbodies ($kT_{BB_1}$ ~ 0.5 keV, $kT_{BB_2}$ ~ 1.2 keV) or an absorbed blackbody plus a power law ($kT_{BB_1}$ ~ 0.6 keV, $Γ$ ~ 3.9). The observed flux was ~ 9 $\times$ 10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, ~ 3 times lower than what observed at the outburst onset, but about one order of magnitude higher than the historical quiescent level. A periodic modulation was detected at the known 6.67 hr periodicity. The spectral decomposition and evolution along the outburst decay are consistent with 1E 161348-5055 being a magnetar, the slowest ever detected.
△ Less
Submitted 1 May, 2018;
originally announced May 2018.
-
Can a bright and energetic X-ray pulsar be hiding amid the debris of SN 1987A?
Authors:
P. Esposito,
N. Rea,
D. Lazzati,
M. Matsuura,
R. Perna,
J. A. Pons
Abstract:
The mass of the stellar precursor of supernova (SN) 1987A and the burst of neutrinos observed at the moment of the explosion are consistent with the core-collapse formation of a neutron star. However, no compelling evidence for the presence in SN 1987A of a compact object of any kind has been found yet in any band of the electromagnetic spectrum, prompting questions on whether the neutron star sur…
▽ More
The mass of the stellar precursor of supernova (SN) 1987A and the burst of neutrinos observed at the moment of the explosion are consistent with the core-collapse formation of a neutron star. However, no compelling evidence for the presence in SN 1987A of a compact object of any kind has been found yet in any band of the electromagnetic spectrum, prompting questions on whether the neutron star survived and, if it did, on its properties. Starting from the analysis of recent Chandra observations, here we appraise the current observational situation. We derived limits on the X-ray luminosity of a compact object with a nonthermal, Crab-pulsar-like spectrum of the order of $\approx$$(1$-$5)\times10^{35}$ erg s$^{-1}$, corresponding to limits on the rotational energy loss of a possible X-ray pulsar in SN 1987A of $\approx$$(0.5$-$1.5)\times10^{38}$ erg s$^{-1}$. However, a much brighter X-ray source cannot be excluded if, as is likely, it is enshrouded in a cloud of absorbing matter with metallicity similar to that expected in the outer layers of a massive star towards the end of its life. We found that other limits obtained from various arguments and observations in other energy ranges either are unbinding or allow a similar maximum luminosity of the order of $\approx$$10^{35}$ erg s$^{-1}$. We conclude that while a pulsar alike the one in the Crab Nebula in both luminosity and spectrum is hardly compatible with the observations, there is ample space for an `ordinary' X-ray-emitting young neutron star, born with normal initial spin period, temperature and magnetic field, to be hiding inside the evolving remnant of SN 1987A.
△ Less
Submitted 13 March, 2018;
originally announced March 2018.
-
The force-free twisted magnetosphere of a neutron star - II: Degeneracies of the Grad-Shafranov equation
Authors:
Taner Akgün,
Pablo Cerdá-Durán,
Juan Antonio Miralles,
José A. Pons
Abstract:
We extend our previous study of equilibrium solutions of non-rotating force-free magnetospheres of neutron stars. We show that multiple solutions exist for the same sets of parameters, implying that the solutions are degenerate. We are able to obtain configurations with disconnected field lines, however, in nearly all cases these correspond to degenerate higher energy solutions. We carry out a wid…
▽ More
We extend our previous study of equilibrium solutions of non-rotating force-free magnetospheres of neutron stars. We show that multiple solutions exist for the same sets of parameters, implying that the solutions are degenerate. We are able to obtain configurations with disconnected field lines, however, in nearly all cases these correspond to degenerate higher energy solutions. We carry out a wide parametric search in order to understand the properties of the solutions. We confirm our previous results that the lower energy solutions have up to $\sim 25\%$ more energy than the vacuum case, helicity of the order of $\sim 5$ (in some defined units), maximum twist of $\sim 1.5$ rad, and a dipole strength that is up to $\sim 40\%$ larger than the vacuum dipole. Including the degenerate higher energy solutions allows for larger theoretical limits of up to $\sim 80\%$ more energy with respect to the vacuum case, helicity of the order of $\sim 8$, and a dipole strength that can now be up to four times that of the vacuum dipole, while the twist can be significantly larger and even diverge for configurations with disconnected domains. The higher energy solutions are probably unstable, therefore, it is unlikely that such magnetospheres exist under normal conditions in magnetars and high magnetic field pulsars.
△ Less
Submitted 10 October, 2018; v1 submitted 28 September, 2017;
originally announced September 2017.
-
Chandra monitoring of the Galactic Centre magnetar SGR J1745-2900 during the initial 3.5 years of outburst decay
Authors:
F. Coti Zelati,
N. Rea,
R. Turolla,
J. A. Pons,
A. Papitto,
P. Esposito,
G. L. Israel,
S. Campana,
S. Zane,
A. Tiengo,
R. P. Mignani,
S. Mereghetti,
F. K. Baganoff,
D. Haggard,
G. Ponti,
D. F. Torres,
A. Borghese,
J. Elfritz
Abstract:
We report on 3.5 years of Chandra monitoring of the Galactic Centre magnetar SGR J1745-2900 since its outburst onset in April 2013. The magnetar spin-down has shown at least two episodes of period derivative increases so far, and it has slowed down regularly in the past year or so. We observed a slightly increasing trend in the time evolution of the pulsed fraction, up to about 55 per cent in the…
▽ More
We report on 3.5 years of Chandra monitoring of the Galactic Centre magnetar SGR J1745-2900 since its outburst onset in April 2013. The magnetar spin-down has shown at least two episodes of period derivative increases so far, and it has slowed down regularly in the past year or so. We observed a slightly increasing trend in the time evolution of the pulsed fraction, up to about 55 per cent in the most recent observations. SGR J1745-2900 has not reached the quiescent level yet, and so far the overall outburst evolution can be interpreted in terms of a cooling hot region on the star surface. We discuss possible scenarios, showing in particular how the presence of a shrinking hot spot in this source is hardly reconcilable with internal crustal cooling and favors the untwisting bundle model for this outburst. Moreover, we also show how the emission from a single uniform hot spot is incompatible with the observed pulsed fraction evolution for any pair of viewing angles, suggesting an anisotropic emission pattern.
△ Less
Submitted 18 July, 2017; v1 submitted 5 July, 2017;
originally announced July 2017.
-
Long-term evolution of the force-free twisted magnetosphere of a magnetar
Authors:
Taner Akgün,
Pablo Cerdá-Durán,
Juan Antonio Miralles,
José A. Pons
Abstract:
We study the long-term quasi-steady evolution of the force-free magnetosphere of a magnetar coupled to its internal magnetic field. We find that magnetospheric currents can be maintained on long timescales of the order of thousands of years. Meanwhile, the energy, helicity and twist stored in the magnetosphere all gradually increase over the course of this evolution, until a critical point is reac…
▽ More
We study the long-term quasi-steady evolution of the force-free magnetosphere of a magnetar coupled to its internal magnetic field. We find that magnetospheric currents can be maintained on long timescales of the order of thousands of years. Meanwhile, the energy, helicity and twist stored in the magnetosphere all gradually increase over the course of this evolution, until a critical point is reached, beyond which a force-free magnetosphere cannot be constructed. At this point, some large-scale magnetospheric rearrangement, possibly resulting in an outburst or a flare, must occur, releasing a large fraction of the stored energy, helicity and twist. After that, the quasi-steady evolution should continue in a similar manner from the new initial conditions. The timescale for reaching this critical point depends on the overall magnetic field strength and on the relative fraction of the toroidal field. The energy stored in the force-free magnetosphere is found to be up to $\sim 30\%$ larger than the corresponding vacuum energy. This implies that for a $10^{14}$ G field at the pole, the energy budget available for fast magnetospheric events is of the order of a few $10^{44}$ erg. The spindown rate is estimated to increase by up to $\sim 60\%$, since the dipole content in the magnetosphere is enhanced by the currents present there. A rough estimate of the braking index $n$ reveals that it is systematically $n < 3$ for the most part of the evolution, consistent with actual measurements for pulsars and early estimates for several magnetars.
△ Less
Submitted 10 October, 2018; v1 submitted 24 June, 2017;
originally announced June 2017.
-
Evolution of a proto-neutron star with a nuclear many-body equation of state: Neutrino luminosity and gravitational wave frequencies
Authors:
Giovanni Camelio,
Alessandro Lovato,
Leonardo Gualtieri,
Omar Benhar,
José A. Pons,
Valeria Ferrari
Abstract:
In a core-collapse supernova, a huge amount of energy is released in the Kelvin-Helmholtz phase subsequent to the explosion, when the proto-neutron star cools and deleptonizes as it loses neutrinos. Most of this energy is emitted through neutrinos, but a fraction of it can be released through gravitational waves. We model the evolution of a proto-neutron star in the Kelvin-Helmholtz phase using a…
▽ More
In a core-collapse supernova, a huge amount of energy is released in the Kelvin-Helmholtz phase subsequent to the explosion, when the proto-neutron star cools and deleptonizes as it loses neutrinos. Most of this energy is emitted through neutrinos, but a fraction of it can be released through gravitational waves. We model the evolution of a proto-neutron star in the Kelvin-Helmholtz phase using a general relativistic numerical code, and a recently proposed finite temperature, many-body equation of state; from this we consistently compute the diffusion coefficients driving the evolution. To include the many-body equation of state, we develop a new fitting formula for the high density baryon free energy at finite temperature and intermediate proton fraction. We estimate the emitted neutrino signal, assessing its detectability by present terrestrial detectors, and we determine the frequencies and damping times of the quasi-normal modes which would characterize the gravitational wave signal emitted in this stage.
△ Less
Submitted 11 August, 2017; v1 submitted 6 April, 2017;
originally announced April 2017.
-
The relevance of ambipolar diffusion for neutron star evolution
Authors:
Andrea Passamonti,
Taner Akgün,
José A. Pons,
Juan A. Miralles
Abstract:
We study ambipolar diffusion in strongly magnetised neutron stars, with special focus on the effects of neutrino reaction rates and the impact of a superfluid/superconducting transition in the neutron star core. For axisymmetric magnetic field configurations, we determine the deviation from $β-$equilibrium induced by the magnetic force and calculate the velocity of the slow, quasi-stationary, ambi…
▽ More
We study ambipolar diffusion in strongly magnetised neutron stars, with special focus on the effects of neutrino reaction rates and the impact of a superfluid/superconducting transition in the neutron star core. For axisymmetric magnetic field configurations, we determine the deviation from $β-$equilibrium induced by the magnetic force and calculate the velocity of the slow, quasi-stationary, ambipolar drift. We study the temperature dependence of the velocity pattern and clearly identify the transition to a predominantly solenoidal flow. For stars without superconducting/superfluid constituents and with a mixed poloidal-toroidal magnetic field of typical magnetar strength, we find that ambipolar diffusion proceeds fast enough to have a significant impact on the magnetic field evolution only at low core temperatures, $T \lesssim 1-2\times10^8$ K. The ambipolar diffusion timescale becomes appreciably shorter when fast neutrino reactions are present, because the possibility to balance part of the magnetic force with pressure gradients is reduced. We also find short ambipolar diffusion timescales in the case of superconducting cores for $T \lesssim 10^9$ K, due to the reduced interaction between protons and neutrons. In the most favourable scenario, with fast neutrino reactions and superconducting cores, ambipolar diffusion results in advection velocities of several km/kyr. This velocity can substantially reorganize magnetic fields in magnetar cores, in a way that can only be confirmed by dynamical simulations.
△ Less
Submitted 15 November, 2016; v1 submitted 29 July, 2016;
originally announced August 2016.
-
Quasi-periodic oscillations in superfluid, relativistic magnetars with nuclear pasta phases
Authors:
Andrea Passamonti,
José A. Pons
Abstract:
We study the torsional magneto-elastic oscillations of relativistic superfluid magnetars and explore the effects of a phase transition in the crust-core interface (nuclear pasta) which results in a weaker elastic response. Exploring various models with different extension of nuclear pasta phases, we find that the differences in the oscillation spectrum present in purely elastic modes (weak magneti…
▽ More
We study the torsional magneto-elastic oscillations of relativistic superfluid magnetars and explore the effects of a phase transition in the crust-core interface (nuclear pasta) which results in a weaker elastic response. Exploring various models with different extension of nuclear pasta phases, we find that the differences in the oscillation spectrum present in purely elastic modes (weak magnetic field), are smeared out with increasing strength of the magnetic field. For magnetar conditions, the main characteristic and features of models without nuclear pasta are preserved. We find in general two classes of magneto-elastic oscillations which exhibit a different oscillation pattern. For B_p < 4 x 10^{14} G, the spectrum is characterised by the turning points and edges of the continuum which are confined into the star's core, and have no constant phase. Increasing the magnetic field, we find, in addition, several magneto-elastic oscillations which reach the surface and have an angular structure similar to crustal modes. These global magneto-elastic oscillations show a constant phase and become dominant when B_p > 5 x 10^{14} G. We do not find any evidence of fundamental pure crustal modes in the low frequency range (below 200 Hz) for B_p >= 10^{14} G.
△ Less
Submitted 29 July, 2016; v1 submitted 7 June, 2016;
originally announced June 2016.
-
The force-free twisted magnetosphere of a neutron star
Authors:
Taner Akgün,
Juan Antonio Miralles,
José A. Pons,
Pablo Cerdá-Durán
Abstract:
We present a detailed analysis of the properties of twisted, force-free magnetospheres of non-rotating neutron stars, which are of interest in the modelling of magnetar properties and evolution. In our models the magnetic field smoothly matches to a current-free (vacuum) solution at some large external radius, and they are specifically built to avoid pathological surface currents at any of the int…
▽ More
We present a detailed analysis of the properties of twisted, force-free magnetospheres of non-rotating neutron stars, which are of interest in the modelling of magnetar properties and evolution. In our models the magnetic field smoothly matches to a current-free (vacuum) solution at some large external radius, and they are specifically built to avoid pathological surface currents at any of the interfaces. By exploring a large range of parameters, we find a few remarkable general trends. We find that the total dipolar moment can be increased by up to $40\%$ with respect to a vacuum model with the same surface magnetic field, due to the contribution of magnetospheric currents to the global magnetic field. Thus, estimates of the surface magnetic field based on the large-scale dipolar braking torque are slightly overestimating the surface value by the same amount. Consistently, there is a moderate increase in the total energy of the model with respect to the vacuum solution of up to $25\%$, which would be the available energy budget in the event of a fast, global magnetospheric reorganization commonly associated with magnetar flares. We have also found the interesting result of the existence of a critical twist ($\varphi_{\rm max} \lesssim 1.5$ rad), beyond which we cannot find any more numerical solutions. Combining the models considered in this paper with the evolution of the interior of neutron stars will allow us to study the influence of the magnetosphere on the long-term magnetic, thermal, and rotational evolution.
△ Less
Submitted 10 October, 2018; v1 submitted 7 May, 2016;
originally announced May 2016.
-
Spin evolution of a proto-neutron star
Authors:
G. Camelio,
L. Gualtieri,
J. A. Pons,
V. Ferrari
Abstract:
We study the evolution of the rotation rate of a proto-neutron star, born in a core-collapse supernova, in the first seconds of its life. During this phase, the star evolution can be described as a sequence of stationary configurations, which we determine by solving the neutrino transport and the stellar structure equations in general relativity. We include in our model the angular momentum loss d…
▽ More
We study the evolution of the rotation rate of a proto-neutron star, born in a core-collapse supernova, in the first seconds of its life. During this phase, the star evolution can be described as a sequence of stationary configurations, which we determine by solving the neutrino transport and the stellar structure equations in general relativity. We include in our model the angular momentum loss due to neutrino emission. We find that the requirement of a rotation rate not exceeding the mass-shedding limit at the beginning of the evolution implies a strict bound on the rotation rate at later times. Moreover, assuming that the proto-neutron star is born with a finite ellipticity, we determine the emitted gravitational wave signal, and estimate its detectability by present and future ground-based interferometric detectors.
△ Less
Submitted 12 July, 2016; v1 submitted 10 January, 2016;
originally announced January 2016.
-
Simulated magnetic field expulsion in neutron star cores
Authors:
J. G. Elfritz,
J. A. Pons,
N. Rea,
K. Glampedakis,
D. Viganò
Abstract:
The study of long-term evolution of neutron star (NS) magnetic fields is key to understanding the rich diversity of NS observations, and to unifying their nature despite the different emission mechanisms and observed properties. Such studies in principle permit a deeper understanding of the most important parameters driving their apparent variety, e.g. radio pulsars, magnetars, x-ray dim isolated…
▽ More
The study of long-term evolution of neutron star (NS) magnetic fields is key to understanding the rich diversity of NS observations, and to unifying their nature despite the different emission mechanisms and observed properties. Such studies in principle permit a deeper understanding of the most important parameters driving their apparent variety, e.g. radio pulsars, magnetars, x-ray dim isolated neutron stars, gamma-ray pulsars. We describe, for the first time, the results from self-consistent magneto-thermal simulations considering not only the effects of the Hall-driven field dissipation in the crust, but adding a complete set of proposed driving forces in a superconducting core. We emphasize how each of these core-field processes drive magnetic evolution and affect observables, and show that when all forces are considered together in vectorial form, the net expulsion of core magnetic flux is negligible, and will have no observable effect in the crust (consequently in the observed surface emission) on megayear time-scales. Our new simulations suggest that strong magnetic fields in NS cores (and the signatures on the NS surface) will persist long after the crustal magnetic field has evolved and decayed, due to the weak combined effects of dissipation and expulsion in the stellar core.
△ Less
Submitted 22 December, 2015;
originally announced December 2015.
-
Are pulsars born with a hidden magnetic field?
Authors:
Alejandro Torres-Forné,
Pablo Cerdá-Durán,
José A. Pons,
José A. Font
Abstract:
The observation of several neutron stars in the center of supernova remnants and with significantly lower values of the dipolar magnetic field than the average radio-pulsar population has motivated a lively debate about their formation and origin, with controversial interpretations. A possible explanation requires the slow rotation of the proto-neutron star at birth, which is unable to amplify its…
▽ More
The observation of several neutron stars in the center of supernova remnants and with significantly lower values of the dipolar magnetic field than the average radio-pulsar population has motivated a lively debate about their formation and origin, with controversial interpretations. A possible explanation requires the slow rotation of the proto-neutron star at birth, which is unable to amplify its magnetic field to typical pulsar levels. An alternative possibility, the hidden magnetic field scenario, considers the accretion of the fallback of the supernova debris onto the neutron star as responsible for the submergence (or screening) of the field and its apparently low value. In this paper we study under which conditions the magnetic field of a neutron star can be buried into the crust due to an accreting, conducting fluid. For this purpose, we consider a spherically symmetric calculation in general relativity to estimate the balance between the incoming accretion flow and the magnetosphere. Our study analyses several models with different specific entropy, composition, and neutron star masses. The main conclusion of our work is that typical magnetic fields of a few times 1e12 G can be buried by accreting only 1e-3 - 1e-2 solar mass, a relatively modest amount of mass. In view of this result, the Central Compact Object scenario should not be considered unusual, and we predict that anomalously weak magnetic fields should be common in very young (< few kyr) neutron stars.
△ Less
Submitted 11 December, 2015; v1 submitted 12 November, 2015;
originally announced November 2015.
-
The Outburst Decay of the Low Magnetic Field Magnetar SWIFT J1822.3-1606: Phase-resolved Analysis and Evidence for a Variable Cyclotron Feature
Authors:
Guillermo A. Rodríguez Castillo,
Gian Luca Israel,
Andrea Tiengo,
David Salvetti,
Roberto Turolla,
Silvia Zane,
Nanda Rea,
Paolo Esposito,
Sandro Mereghetti,
Rosalba Perna,
Luigi Stella,
José A. Pons,
Sergio Campana,
Diego Götz,
Sara Motta
Abstract:
We study the timing and spectral properties of the low-magnetic field, transient magnetar SWIFT J1822.3-1606 as it approached quiescence. We coherently phase-connect the observations over a time-span of ~500 days since the discovery of SWIFT J1822.3-1606 following the Swift-BAT trigger on 2011 July 14, and carried out a detailed pulse phase spectroscopy along the outburst decay. We follow the spec…
▽ More
We study the timing and spectral properties of the low-magnetic field, transient magnetar SWIFT J1822.3-1606 as it approached quiescence. We coherently phase-connect the observations over a time-span of ~500 days since the discovery of SWIFT J1822.3-1606 following the Swift-BAT trigger on 2011 July 14, and carried out a detailed pulse phase spectroscopy along the outburst decay. We follow the spectral evolution of different pulse phase intervals and find a phase and energy-variable spectral feature, which we interpret as proton cyclotron resonant scattering of soft photon from currents circulating in a strong (10^14 G) small-scale component of the magnetic field near the neutron star surface, superimposed to the much weaker (~3 x 10^13 G) magnetic dipole field. We discuss also the implications of the pulse-resolved spectral analysis for the emission regions on the surface of the cooling magnetar.
△ Less
Submitted 30 October, 2015;
originally announced October 2015.
-
Constraining the GRB-magnetar model by means of the Galactic pulsar population
Authors:
Nanda Rea,
Miguel Gullon,
Jose' A. Pons,
Rosalba Perna,
Maria G. Dainotti,
Juan A. Miralles,
Diego F. Torres
Abstract:
A large fraction of Gamma Ray Bursts (GRBs) displays an X-ray plateau phase within <10^{5} s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, betw…
▽ More
A large fraction of Gamma Ray Bursts (GRBs) displays an X-ray plateau phase within <10^{5} s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, between January 2005 and August 2014. From the derived initial magnetic field distribution for the possible magnetars left behind by the GRBs, we study the evolution and properties of a simulated GRB-magnetar population using numerical simulations of magnetic field evolution, coupled with Monte Carlo simulations of Pulsar Population Synthesis in our Galaxy. We find that if the GRB X-ray plateaus are powered by the rotational energy of a newly formed magnetar, the current observational properties of the Galactic magnetar population are not compatible with being formed within the GRB scenario (regardless of the GRB type or rate at z=0). Direct consequences would be that we should allow the existence of magnetars and "super-magnetars" having different progenitors, and that Type Ib/c SNe related to Long GRBs form systematically neutron stars with higher initial magnetic fields. We put an upper limit of <16 "super-magnetars" formed by a GRB in our Galaxy in the past Myr (at 99% c.l.). This limit is somewhat smaller than what roughly expected from Long GRB rates, although the very large uncertainties do not allow us to draw strong conclusion in this respect.
△ Less
Submitted 6 October, 2015;
originally announced October 2015.
-
Neutron stars - cooling and transport
Authors:
A. Y. Potekhin,
J. A. Pons,
Dany Page
Abstract:
Observations of thermal radiation from neutron stars can potentially provide information about the states of supranuclear matter in the interiors of these stars with the aid of the theory of neutron-star thermal evolution. We review the basics of this theory for isolated neutron stars with strong magnetic fields, including most relevant thermodynamic and kinetic properties in the stellar core, cru…
▽ More
Observations of thermal radiation from neutron stars can potentially provide information about the states of supranuclear matter in the interiors of these stars with the aid of the theory of neutron-star thermal evolution. We review the basics of this theory for isolated neutron stars with strong magnetic fields, including most relevant thermodynamic and kinetic properties in the stellar core, crust, and blanketing envelopes.
△ Less
Submitted 18 October, 2017; v1 submitted 22 July, 2015;
originally announced July 2015.
-
Population Synthesis of Isolated Neutron Stars with magneto-rotational evolution II: from radio-pulsars to magnetars
Authors:
M. Gullón,
J. A. Pons,
J. A. Miralles,
D. Viganò,
N. Rea,
R. Perna
Abstract:
Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include $X$-ray…
▽ More
Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include $X$-ray thermal emitting pulsars within the same evolutionary model as radio-pulsars. We find that the cumulative distribution of the number of X-ray pulsars can be well reproduced by several models that, simultaneously, reproduce the characteristics of the radio-pulsar distribution. However, even considering the most favourable magneto-thermal evolution models with fast field decay, log-normal distributions of the initial magnetic field over-predict the number of visible sources with periods longer than 12 s. We then show that the problem can be solved with different distributions of magnetic field, such as a truncated log-normal distribution, or a binormal distribution with two distinct populations. We use the observational lack of isolated NSs with spin periods P>12 s to establish an upper limit to the fraction of magnetars born with B > 10^{15} G (less than 1\%). As future detections keep increasing the magnetar and high-B pulsar statistics, our approach can be used to establish a severe constraint on the maximum magnetic field at birth of NSs.
△ Less
Submitted 20 July, 2015;
originally announced July 2015.
-
The Spanish Square Kilometre Array White Book
Authors:
J. A. Acosta-Pulido,
I. Agudo,
A. Alberdi,
J. Alcolea,
E. J. Alfaro,
A. Alonso-Herrero,
G. Anglada,
P. Arnalte-Mur,
Y. Ascasibar,
B. Ascaso,
R. Azulay,
R. Bachiller,
A. Baez-Rubio,
E. Battaner,
J. Blasco,
C. B. Brook,
V. Bujarrabal,
G. Busquet,
M. D. Caballero-Garcia,
C. Carrasco-Gonzalez,
J. Casares,
A. J. Castro-Tirado,
L. Colina,
F. Colomer,
I. de Gregorio-Monsalvo
, et al. (94 additional authors not shown)
Abstract:
The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on…
▽ More
The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.
△ Less
Submitted 17 June, 2015; v1 submitted 10 June, 2015;
originally announced June 2015.
-
The X-ray outburst of the Galactic Centre magnetar SGR J1745-2900 during the first 1.5 year
Authors:
F. Coti Zelati,
N. Rea,
A. Papitto,
D. Viganò,
J. A. Pons,
R. Turolla,
P. Esposito,
D. Haggard,
F. K. Baganoff,
G. Ponti,
G. L. Israel,
S. Campana,
D. F. Torres,
A. Tiengo,
S. Mereghetti,
R. Perna,
S. Zane,
R. P. Mignani,
A. Possenti,
L. Stella
Abstract:
In 2013 April a new magnetar, SGR 1745-2900, was discovered as it entered an outburst, at only 2.4 arcsec angular distance from the supermassive black hole at the Centre of the Milky Way, Sagittarius A*. SGR 1745-2900 has a surface dipolar magnetic field of ~ 2x10^{14} G, and it is the neutron star closest to a black hole ever observed. The new source was detected both in the radio and X-ray bands…
▽ More
In 2013 April a new magnetar, SGR 1745-2900, was discovered as it entered an outburst, at only 2.4 arcsec angular distance from the supermassive black hole at the Centre of the Milky Way, Sagittarius A*. SGR 1745-2900 has a surface dipolar magnetic field of ~ 2x10^{14} G, and it is the neutron star closest to a black hole ever observed. The new source was detected both in the radio and X-ray bands, with a peak X-ray luminosity L_X ~ 5x10^{35} erg s^{-1}. Here we report on the long-term Chandra (25 observations) and XMM-Newton (8 observations) X-ray monitoring campaign of SGR 1745-2900, from the onset of the outburst in April 2013 until September 2014. This unprecedented dataset allows us to refine the timing properties of the source, as well as to study the outburst spectral evolution as a function of time and rotational phase. Our timing analysis confirms the increase in the spin period derivative by a factor of ~2 around June 2013, and reveals that a further increase occurred between 2013 Oct 30 and 2014 Feb 21. We find that the period derivative changed from 6.6x10^{-12} s s^{-1} to 3.3x10^{-11} s s^{-1} in 1.5 yr. On the other hand, this magnetar shows a slow flux decay compared to other magnetars and a rather inefficient surface cooling. In particular, starquake-induced crustal cooling models alone have difficulty in explaining the high luminosity of the source for the first ~200 days of its outburst, and additional heating of the star surface from currents flowing in a twisted magnetic bundle is probably playing an important role in the outburst evolution.
△ Less
Submitted 4 March, 2015;
originally announced March 2015.
-
Magnetic fields in Neutron Stars
Authors:
Daniele Viganò,
Jose A. Pons,
Juan A. Miralles,
Nanda Rea
Abstract:
Isolated neutron stars show a diversity in timing and spectral properties, which has historically led to a classification in different sub-classes. The magnetic field plays a key role in many aspects of the neutron star phenomenology: it regulates the braking torque responsible for their timing properties and, for magnetars, it provides the energy budget for the outburst activity and high quiescen…
▽ More
Isolated neutron stars show a diversity in timing and spectral properties, which has historically led to a classification in different sub-classes. The magnetic field plays a key role in many aspects of the neutron star phenomenology: it regulates the braking torque responsible for their timing properties and, for magnetars, it provides the energy budget for the outburst activity and high quiescent luminosities (usually well above the rotational energy budget). We aim at unifying this observational variety by linking the results of the state-of-the-art 2D magneto-thermal simulations with observational data. The comparison between theory and observations allows to place two strong constraints on the physical properties of the inner crust. First, strong electrical currents must circulate in the crust, rather than in the star core. Second, the innermost part of the crust must be highly resistive, which is in principle in agreement with the presence of a novel phase of matter so-called nuclear pasta phase.
△ Less
Submitted 27 January, 2015;
originally announced January 2015.
-
Neutron stars - thermal emitters
Authors:
A. Y. Potekhin,
A. De Luca,
J. A. Pons
Abstract:
Confronting theoretical models with observations of thermal radiation emitted by neutron stars is one of the most important ways to understand the properties of both, superdense matter in the interiors of the neutron stars and dense magnetized plasmas in their outer layers. Here we review the theory of thermal emission from the surface layers of strongly magnetized neutron stars, and the main prop…
▽ More
Confronting theoretical models with observations of thermal radiation emitted by neutron stars is one of the most important ways to understand the properties of both, superdense matter in the interiors of the neutron stars and dense magnetized plasmas in their outer layers. Here we review the theory of thermal emission from the surface layers of strongly magnetized neutron stars, and the main properties of the observational data. In particular, we focus on the nearby sources for which a clear thermal component has been detected, without being contaminated by other emission processes (magnetosphere, accretion, nebulae). We also discuss the applications of the modern theoretical models of the formation of spectra of strongly magnetized neutron stars to the observed thermally emitting objects.
△ Less
Submitted 24 December, 2015; v1 submitted 26 September, 2014;
originally announced September 2014.
-
Population synthesis of isolated Neutron Stars with magneto--rotational evolution
Authors:
Miguel Gullón,
Juan A. Miralles,
Daniele Viganò,
José A. Pons
Abstract:
We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in par…
▽ More
We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in parameter space that are more consistent with the observational data. In particular, we find that any broad distribution of birth spin periods with $P_0\lesssim 0.5$ s can fit the data, and that if the alignment angle is allowed to vary consistently with the torque model, realistic magnetospheric models are favoured compared to models with classical magneto-dipolar radiation losses. Assuming that the initial magnetic field is given by a lognormal distribution, our optimal model has mean strength $\langle\log B_0{\rm [G]}\rangle \approx 13.0-13.2$ with width $σ(\log B_0) = 0.6-0.7$. However, there are strong correlations between parameters. This degeneracy in the parameter space can be broken by an independent estimate of the pulsar birth rate or by future studies correlating this information with the population in other observational bands (X-rays and $γ$-rays).
△ Less
Submitted 26 June, 2014;
originally announced June 2014.