Search

The braking torque that dictates the timing properties of magnetars is closely tied to the large-scale dipolar magnetic field on their surface. The formation of this field has been a topic of ongoing debate. One proposed mechanism, based on macroscopic principles, involves an inverse cascade within the neutron star's crust. However, this phenomenon has not been observed in realistic simulations. In this study, we provide compelling evidence supporting the feasibility of the inverse cascading process in the presence of an initial helical magnetic field within realistic neutron star crusts and discuss its contribution to the amplification of the large-scale magnetic field. Our findings, derived from a systematic investigation that considers various coordinate systems, peak wavenumber positions, crustal thicknesses, magnetic boundary conditions, and magnetic Lundquist numbers, reveal that the specific geometry of the crustal domain - with its extreme aspect ratio - requires an initial peak wavenumber from small-scale structures for the inverse cascade to occur. However, this extreme aspect ratio limits the inverse cascade to magnetic field structures on scales comparable to the neutron star's crust, making the formation of a large-scale dipolar surface field unlikely. Despite this limitation, the inverse cascade can significantly impact the magnetic field evolution in the interior of the crust, potentially explaining the observed characteristics of highly magnetized objects with weak surface dipolar fields, such as low-field magnetars or central compact objects., Comment: Submitted to A&A, 12 pages, 8 Figures. Comments are welcome

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., Comment: Submitted to PRL, 6 pages, 3 figures. Comments are welcome

This doctoral thesis investigates the long-term evolution of the strong magnetic fields within isolated neutron stars (NSs), the most potent magnetic objects in the universe. Their magnetic influence extends beyond their surface to encompass the magnetised plasma in their vicinity. The overarching magnetic configuration significantly impacts the observable characteristics of the highly magnetised NSs, i.e., magnetars. Conversely, the internal magnetic field undergoes prolonged evolution spanning thousands to millions of years, intricately linked to thermal evolution. The diverse observable phenomena associated with NSs underscore the complex 3D nature of their magnetic structure, thereby requiring sophisticated numerical simulations. A central focus of this thesis involves a thorough exploration of state-of-the-art 3D coupled magneto-thermal evolution models. This marks a pioneering achievement as we conduct, for the first time, the most realistic 3D simulations to date, spanning the first million years of a NS's life using the newly developed code MATINS, which adeptly accounts for both Ohmic dissipation and Hall drift within the NS's crust. Our simulations incorporate highly accurate temperature-dependent microphysical calculations and adopt the star's structure based on a realistic equation of state. To address axial singularities in 3D simulations, we employ the cubed-sphere coordinates. We also account for corresponding relativistic factors in the evolution equations and use the latest envelope model from existing literature, in addition to an initial magnetic field structure derived from proton-NS dynamo simulations. Within this framework, we quantitatively simulate the thermal luminosity, timing properties, and magnetic field evolution, pushing the boundaries of numerical modeling capabilities and enabling the performance of several astrophysical studies within this thesis., Comment: PhD Thesis, defended on November 8th, 2023, at the Universitat Aut\`onoma de Barcelona, Spain

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., Comment: 60 pages, 10 figures, 7 tables; Accepted for publication; first four authors are corresponding authors

We conducted a study on the thermal properties of stellar matter with the nuclear energy density functional BCPM. This functional is based on microscopic Brueckner-Hartree-Fock calculations and has demonstrated success in describing cold neutron stars. To enhance its applicability in astrophysics, in this study we extend the BCPM equation of state to finite temperature for $\beta$-stable neutrino-free matter, taking into consideration the hot inner crust. Such an equation of state holds significant importance for hot compact objects, particularly those resulting from a binary neutron star merger event. Our exploration has shown that with increasing temperature there is a fast decrease of the crust-core transition density, suggesting that for hot stars it is not realistic to assume a fixed value of this density. The microscopic calculations also reveal that the presence of nuclear clusters persists up to $T=7.21$ MeV, identified as the limiting temperature of the crust. Above this threshold, the manifestation of clusters is not anticipated. Below this temperature, clusters within the inner crust are surrounded by uniform matter with varying densities, allowing for the distinction between the upper and lower transition density branches. Moreover, we computed mass--radius relations of neutron stars, assuming an isothermal profile for $\beta$-stable neutron star matter at various temperature values. Our findings highlight the significant influence of the hot inner crust on the mass--radius relationship, leading to the formation of larger and more inflated neutron stars. Consequently, under our prescription, the final outcome is a unified equation of state at finite temperature., Comment: Submitted to A&A, 9 pages, 6 figures

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., Comment: 26 pages, 16 figures; published in MNRAS

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., Comment: 8 pages, 6 figures, submitted for publication in MNRAS, comments are welcome

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., Comment: 11 pages, 10 figures, submitted for publication in MNRAS

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., Comment: 5 pages, 4 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Comment: 21 pages, 17 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Comment: 27 pages, 18 figures, 3 tables Accepted in Computer Physics Communications

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., Comment: 9 pages, 4 figures, ApJ Letter in press