Search

Scallop-shell stars, a recently discovered class of young M dwarfs, show complex optical light curves that are characterized by periodic dips as well as other features that are stable over tens to hundreds of rotation cycles. The origin of these features is not well-understood. 2MASS J05082729$-$2101444 is a $\sim$25 Myr old scallop-shell star that was identified using TESS data; it has a photometric period of 6.73h that has been attributed to rotation. Of the $\sim$50 recently confirmed scallop-shell stars, it is one of the few detected at radio frequencies between 1 and 8 GHz. We observed this rare system with the upgraded Giant Meterwave Radio Telescope at 575--720 MHz, covering 88% of the photometric period in each of the two observations scheduled almost a month apart in 2023. We detected $\sim$millijansky emission from the target in both epochs, with a significant circular polarization fraction: $|V/I|\sim$20--50%. The 3.5-min phase-folded light curves reveal unique variability in circular polarization, showing an $\sim$hour-long helicity reversal in both epochs, similar in amplitude, length, and (possibly) phase. These results suggest two emission components: The first is a persistent, moderately polarized component possibly ascribable to gyro-synchrotron emission driven by centrifugal breakout events. The second is a highly polarized, short burst-like component, likely due to an electron cyclotron maser (ECM), indicative of auroral emission and potentially responsible for the helicity reversal. To explain this, we discuss the different origins of the plasma responsible for the radio emission, including the possibility that the occulting material is acting as a plasma source. Future coordinated multifrequency radio and optical observations can further constrain the underlying scenario, as well as the magnetic geometry of the system, if we assume an ECM-like auroral emission., Comment: 9 pages, 5 figures, 1 table, Accepted for publication in Astronomy & Astrophysics Letters

The young stellar object [BHB2007]-1 has been extensively studied in the past at radio, millimeter, and infrared wavelengths. It shows a gap in the disk and previous observations claimed the possible emission from a forming sub-stellar object, in correspondence to the disk gap. Here, we analyze a set of 8 Karl Jansky Very Large Array (VLA) observations at 15 GHz and spread over a month. We infer a slowly variable emission from the star, with a $\sim 15 \text{-} 20\%$ circular polarization detected in two of the eight observations. The latter can be related to the magnetic fields in the system, while the unpolarized and moderately varying component can be indicative of free-free emission associated with jet induced shocks or interaction of the stellar wind with dense surrounding material. We discard any relevant short flaring activities when sampling the radio light curves down to 10 seconds and find no clear evidence of emission from the sub-stellar object inferred from past observations, although deeper observations could shed further light on this., Comment: 9 pages, 6 figures, 2 tables, Accepted for publication in A&A

We analyze a long-lived hyper-massive neutron star merger remnant (post-merger lifetime $>250$ ms) that has been obtained via large eddy simulations with a gradient subgrid-scale model. We find a clear helicoidal magnetic field structure that is governed by the toroidal component of the magnetic field. Although no jet emerges during the simulation time, we observe at late times a significant increase of the poloidal component of the magnetic field at all scales. We also compare with the results of several binary neutron star simulations with moderate resolution of $120$~m, that are evolved up to $50$~ms after the merger, which differ in terms of the initial topology and strength of the magnetic field. We find that the best choice is an isotropic small-scale magnetic field distribution that mimics the turbulent state that generically develops during the merger. This initial configuration reaches a closer agreement with our high-resolution simulation results than the purely dipolar large-scale fields that are commonly employed in these type of simulations. This provides a recipe to perform such simulations avoiding the computationally expensive grids required to faithfully capture the amplification of the magnetic field by Kelvin-Helmholtz instabilities., Comment: 13 pages, 11 figures

The amplification of astrophysical magnetic fields takes place via dynamo instability in turbulent environments. The presence of vorticity is crucial for the dynamo to happen. However, the role of vorticity is not yet fully understood. This work is an extension of previous research on the effect of an irrotational subsonic forcing on a magnetized medium in the presence of rotation or a differential velocity profile, aimed at exploring a wider parameter space in terms of Reynolds numbers, magnetic Prandtl number, forcing scale, cooling timescale in a Newtonian cooling. We study the effect of imposing either the acceleration or the velocity forcing function to be curl-free and evaluate the terms responsible for the evolution vorticity. We use Direct Numerical Simulations (DNS) to solve the fully compressible, resistive magnetohydrodynamical (MHD) equations with the Pencil Code. We study both isothermal and non-isothermal regimes and address the relative importance of different vorticity source terms.We report no small-scale dynamo for the models that do not include shear. We find a hydro instability, followed by a magnetic one, when a shearing velocity profile is applied. The vorticity production is found to be numerical in the purely irrotational case. Non-isothermality, rotation, shear or forcing in the form of a velocity curl-free, when included, contribute to increasing vorticity. Consistently with our previous study, we find that turbulence driven by subsonic expansion waves can amplify vorticity and magnetic field only in the presence of a background shearing profile. The presence of a cooling function make the instability happens on a shorter timescale. We estimate critical Reynolds and Magnetic Reynolds Numbes of 40 and 20, respectively., Comment: Submitted for publication. Comments are welcome

This work presents a methodological approach to generate realistic $\gamma$-ray light curves of pulsars, resembling reasonably well the observational ones observed by the Fermi-Large Area Telescope instrument, fitting at the same time their high-energy spectra. The theoretical light curves are obtained from a spectral and geometrical model of the synchro-curvature emission. Despite our model relies on a few effective physical parameters, the synthetic light curves present the same main features observed in the observational $\gamma$-ray light curve zoo, such as the different shapes, variety in the number of peaks, and a diversity of peak widths. The morphological features of the light curves allows us to statistically compare the observed properties. In particular, we find that the proportion on the number of peaks found in our synthetic light curves is in agreement with the observational one provided by the third Fermi-LAT pulsar catalog. We also found that the detection probability due to beaming is much higher for orthogonal rotators (approaching 100%) than for small inclination angles (less than 20%).The small variation on the synthetic skymaps generated for different pulsars indicates that the geometry dominates over timing and spectral properties in shaping the gamma-ray light curves. This means that geometrical parameters like the inclination angle can be in principle constrained by gamma-ray data alone independently on the specific properties of a pulsar. At the same time, we find that $\gamma$-ray spectra seen by different observers can slightly differ, opening the door to constraining the viewing angle of a particular pulsar., Comment: 14 pages, 8 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The inflated radii observed in hundreds of Hot Jupiters represent a long-standing open issue. The observed correlation between radii and irradiation strength, and the occasional extreme cases, nearly double the size of Jupiter, remain without a comprehensive quantitative explanation. In this investigation, we delve into this issue within the framework of Ohmic dissipation, one of the most promising mechanisms for explaining the radius anomaly. Using the evolutionary code MESA, we simulate the evolution of irradiated giant planets, spanning the range 1 to 8 Jupiter masses, incorporating an internal source of Ohmic dissipation located beneath the radiative-convective boundary. Our modeling is based on physical parameters, and accounts for the approximated conductivity and the evolution of the magnetic fields, utilizing widely-used scaling laws. We compute the radius evolution across a spectrum of masses and equilibrium temperatures, considering varying amounts of Ohmic dissipation, calculated with the internal conductivity profile and an effective parametrization of the currents, based on the typical radius of curvature of the field lines. Our analysis reveals that this internal Ohmic dissipation can broadly reproduce the range of observed radii using values of radius of curvature up to about one order of magnitude lower than what we estimate from the Juno measurements of the Jovian magnetosphere and from MHD dynamo simulations presented herein. The observed trend with equilibrium temperature can be explained if the highly-irradiated planets have more intense and more small-scale magnetic fields. This suggests the possibility of an interplay between atmospherically induced currents and the interior, via turbulence, in agreement with recent box simulations of turbulent MHD in atmospheric columns., Comment: 13 pages, 14 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 perform a long and accurate Large-Eddy Simulation of a binary neutron star merger, following the newly formed remnant up to 110 milliseconds. The combination of high-order schemes, high-resolution and the gradient subgrid-scale model allow us to have among the highest effective resolutions ever achieved. Our results show that, although the magnetic fields are strongly amplified by the Kelvin-Helmholtz instability, they are coherent only over very short spatial scales until t \gtrsim 30 ms. Around that time, magnetic winding becomes more efficient leading to a linear growth of the toroidal component and slowly ordering the field to more axisymmetric, large scales. The poloidal component only starts to grow at small scales at much later times t \gtrsim 90 ms, in a way compatible with the magneto-rotational instability. No strong large-scale poloidal field or jet is produced in the timescales spanned by our simulation, although there is an helicoidal structure gradually developing at late times. We highlight that soon after the merger the topology is always strongly dominated by toroidal structures, with a complex distribution in the meridional plane and highly turbulent perturbations. Thus, starting with strong purely dipolar fields before the merger is largely inconsistent with the outcomes of a realistic evolution. Finally, we confirm the universality of the evolved topology, even when starting with very different magnetic fields confined to the outermost layers of each neutron star., Comment: 18 pages, 14 figures. arXiv admin note: text overlap with arXiv:2112.08413

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

This work concentrates on the effect of an irrotational forcing on a magnetized flow in the presence of rotation, baroclinicity, shear, or a combination of them. By including magnetic field in the model we can evaluate the occurrence of dynamo on both small and large scales. We aim at finding what are the minimum ingredients needed to trigger a dynamo instability and what is the relation between dynamo and the growth of vorticity. We use the Pencil code to run resistive MHD direct numerical simulations. We report no dynamo in all cases where only rotation is included, regardless on the equation of state. Conversely, the inclusion of a background sinusoidal shearing profile leads to an hydrodynamic instability that produces an exponential growth of the vorticity at all scales, starting from small ones. This is know as vorticity dynamo. The onset of this instability occurs after a rather long temporal evolution of several thousand turbulent turnover times. The vorticity dynamo in turn drives an exponential growth of the magnetic field, first at small scales, then also at large one. The instability then saturates and the magnetic field approximately reaches equipartition with the turbulent kinetic energy. During the saturation phase we can observe a winding of the magnetic field in the direction of the shearing flow. By varying the intensity of the shear we see that the growth rates of this instability change. The inclusion of the baroclinic term delays the onset of the vorticity dynamo but leads to a more rapid growth. We demonstrate how in the presence of shear, even a purely irrotational forcing amplifies the field to equipartition. At the same time, we confirm how this forcing alone does not lead to vorticity nor magnetic field growth, and this picture does not change in the presence of rotation or baroclinicity up to $256^3$ meshpoints., Comment: 14 pages, 8 figures, submitted for publication, comments welcome

While magnetism in exoplanets remains largely unknown, Hot Jupiters have been considered as natural candidates to harbour intense magnetic fields, both due to their large masses and their high energy budgets coming from irradiation as a consequence of their vicinity to their host stars. In this work we perform MHD simulations of a narrow day-side atmospheric column of ultra-hot Jupiters, suitable for very high local temperatures (T > 3000 K). Since the conductivity in this regime is very high, the dominant effect is winding due to the intense zonal winds. By including a forcing that mimics the wind profiles obtained in global circulation models, the shear layer induces a strong toroidal magnetic field (locally reaching hundreds of gauss), supported by meridional currents. Such fields and the sustaining currents don$'$t depend on the internally generated field, but are all confined in the thin (less than a scale-height) shear layer around 1 bar. Additionally, we add random perturbations that induce turbulent motions, which lead to further (but much smaller) magnetic field generation to a broader range of depths. These results allow an evaluation of the currents induced by the atmospheric dynamo. Although here we use ideal MHD and the only resistivity comes from the numerical scheme, we estimate a-posteriori the amount of Ohmic heat deposited in the outer layers, which could be employed in evolutionary models for Hot Jupiters' inflated radii., Comment: Accepted 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

Several important properties of rotation-powered millisecond pulsars (MSPs), such as their mass-radius ratio, equation of state and magnetic field topology, can be inferred from precise observations and modelling of their X-ray light curves. In the present study, we model the thermal X-ray signals originated in MSPs, all the way from numerically solving the surrounding magnetospheres up to the ray tracing of the emitted photons and the final computation of their light curves and spectra. Our modelled X-ray signals are then compared against a set of very accurate NICER observations of four target pulsars: PSR J0437-4715, PSR J1231-1411, PSR J2124-3358 and PSR J0030+0451. We find very good simultaneous fits for the light curve and spectral distribution in all these pulsars. The magnetosphere is solved by performing general relativistic force-free simulations of a rotating neutron star (NS) endowed with a simple centered dipolar magnetic field, for many different stellar compactness and pulsar misalignments. From these solutions, we derive an emissivity map over the surface of the star, which is based on the electric currents in the magnetosphere. In particular, the emission regions (ERs) are determined in this model by spacelike four-currents that reach the NS. We show that this assumption, together with the inclusion of the gravitational curvature on the force-free simulations, lead to non-standard ERs facing the closed-zone of the pulsar, in addition to other ERs within the polar caps. The combined X-ray signals from these two kinds of ERs (both antipodal) allow to approximate the non-trivial interpulses found in all the target MSPs light curves., Comment: 16 pages, 9 figures

The detected high-energy pulsars' population is growing in number, and thus, having agile and physically relevant codes to analyze it consistently is important. Here, we update our existing synchro-curvature radiation model by including a better treatment of the particle injection, particularly where the large pitch angle particles dominate the spectra, and by implementing a fast and accurate minimization technique. The latter allows a large improvement in computational cost, needed to test model enhancements and to apply the model to a larger pulsar population. We successfully fit the sample of pulsars with X-ray and $\gamma$-ray data. Our results indicate that, for every emitting particle, the spatial extent of their trajectory where the pitch angle is large and most of the detected X-ray radiation is produced is a small fraction of the light cylinder. We also confirm with this new approach that synchrotron radiation is not negligible for most of the gamma-ray pulsars detected. In addition, with the results obtained, we argue that J0357+3205 and J2055+2539 are MeV-pulsar candidates and are suggested for exhaustive observations in this energy band., Comment: 12 pages, 7 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The period and the period derivative of a pulsar are critical magnitudes for defining the properties of the magnetospheric size and plasma dynamics. The pulsar light cylinder, the magnetic field intensity nearby it, and the curvature radius all depend on these timing properties, and shape the observed high-energy synchro-curvature emission. Therefore, the radiative properties of pulsars are inextricably linked to them. This fact poses the question of how well does a given pulsar's spectral energy distribution embeds information of the timing parameters, and if so, whether we can deduce them if they have not been measured directly. This is relevant to possibly constrain the timing properties of potential pulsar candidates among unidentified $\gamma$-ray sources. We consider well-measured pulsar spectra blinding us from the knowledge of their timing properties, and address this question by using our radiative synchro-curvature model that was proven able to fit the observed spectra of the pulsar population. We find that in the majority of the cases studied (8/13), the spin period is constrained within a range of about one order of magnitude, within which the real period lies. In the other cases, there is degeneracy and no period range can be constrained. This can be used to facilitate the blind search of pulsed signals., Comment: 8 pages, 4 figures. Accepted for publication in The Astrophysical Journal

The detection of a binary neutron star merger in 2017 through both gravitational waves and electromagnetic emission opened a new era of multimessenger astronomy. The understanding of the magnetic field amplification triggered by the Kelvin-Helmholtz instability during the merger is still a numerically unresolved problem because of the relevant small scales involved. One of the uncertainties comes from the simplifications usually assumed in the initial magnetic topology of merging neutron stars. We perform high-resolution, convergent large-eddy simulations of binary neutron star mergers, following the newly formed remnant for up to $30$ milliseconds. Here we specifically focus on the comparison between simulations with different initial magnetic configurations, going beyond the widespread-used aligned dipole confined within each star. The results obtained show that the initial topology is quickly forgotten, in a timescale of few miliseconds after the merger. Moreover, at the end of the simulations, the average intensity ($B\sim 10^{16}$ G) and the spectral distribution of magnetic energy over spatial scales barely depend on the initial configuration. This is expected due to the small-scale efficient dynamo involved, and thus it holds as long as: (i) the initial large-scale magnetic field is not unrealistically high (as often imposed in mergers studies); (ii) the turbulent instability is numerically (at least partially) resolved, so that the amplified magnetic energy is distributed across a wide range of scales and becomes orders-of-magnitude larger than the initial one., Comment: 8 pages, 5 figures

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

The detection of binary neutron star mergers represents one of the most important and complex astrophysical discoveries of the recent years. One of the unclear aspects of the problem is the turbulent magnetic field amplification, initially triggered by the Kelvin-Helmholtz instability at much smaller scales than any reachable numerical resolution nowadays. Here we present numerical simulations of the first ten milliseconds of a binary neutron star merger. First, we confirm in detail how the simulated amplification depends on the numerical resolution and is distributed on a broad range of scales, as expected from turbulent MHD theory. We find that an initial large-scale magnetic field of $10^{11}\,$G inside each star is amplified in the remnant to root-mean-square values above $10^{16}\,$G within the first $5$ milliseconds for our highest-resolution run. Then, we run large eddy simulations, exploring the performance of the subgrid-scale gradient model, already tested successfully in previous turbulent box simulations. We show that the addition of this model is especially important in the induction equation, since it leads to an amplification of the magnetic field comparable to a higher-resolution run, but with a greatly reduced computational cost. In the first 10 milliseconds, there is no clear hint for an ordered, large-scale magnetic field, which should indeed occur in longer timescales through magnetic winding and the magneto-rotational instability., Comment: 14 pages, 7 figures

We present the design, manufacturing, and testing of a 37-element array of corrugated feedhorns for Cosmic Microwave Background (CMB) measurements between $140$ and $170$ GHz. The array was designed to be coupled to Kinetic Inductance Detector arrays, either directly (for total power measurements) or through an orthomode transducer (for polarization measurements). We manufactured the array in platelets by chemically etching aluminum plates of $0.3$ mm and $0.4$ mm thickness. The process is fast, low-cost, scalable, and yields high-performance antennas compared to other techniques in the same frequency range. Room temperature electromagnetic measurements show excellent repeatability with an average cross polarization level about $-20$ dB, return loss about $-25$ dB, first sidelobes below $-25$ dB and far sidelobes below $-35$ dB. Our results qualify this process as a valid candidate for state-of-the-art CMB experiments, where large detector arrays with high sensitivity and polarization purity are of paramount importance in the quest for the discovery of CMB polarization $B$-modes.

The detection of binary neutron star mergers represents one of the most important astrophysical discoveries of the recent years. Due to the extreme matter and gravity conditions and the rich dynamics developed, it becomes a tremendous challenge to accurately simulate numerically all the scales present during the collision. Here we present how to study such systems by using large eddy simulations with a self-consistent subgrid-scale gradient model, that we generalized to the special relativistic case in a previous work and now extend to the general relativistic case. Adapted from nonrelativistic scenarios, the so-called gradient model allows to capture part of the effects of the hidden dynamics on the resolved scales, by means of a physically-agnostic, mathematically-based Taylor expansion of the nonlinear terms in the conservative evolution equations' fluxes. We assess the validity of this approach in bounding-box simulations of the magnetic Kelvin-Helmholtz instability. Several resolutions and a broad range of scenarios are considered in order to carefully test the performance of the model under three crucial aspects: (i) highly curved backgrounds, (ii) jumps on the fluid density profiles and (iii) strong shocks. The results suggest our extension of the gradient subgrid-scale model to general relativistic magnetohydrodynamics is a promising approach for studying binary neutron stars mergers, and potentially to other relevant astrophysical scenarios., Comment: 15 pages, 7 figures Submitted to PRD

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., Comment: 10 pages, 7 figures, accepted for publication on A&A

We report on a detailed spectral characterization of the non-thermal X-ray emission for a large sample of gamma-ray pulsars in the second Fermi-LAT catalogue. We outline the criteria adopted for the selection of our sample, its completeness, and critically describe different approaches to estimate the spectral shape and flux of pulsars. We perform a systematic modelling of the pulsars' X-ray spectra using archival observations with XMM-Newton, Chandra, and NuSTAR and extract the corresponding non-thermal X-ray spectral distributions. This set of data is made available online and is useful to confront with predictions of theoretical models., Comment: MNRAS in press

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., Comment: Accepted for Living Reviews in Computational Astrophysics. 64 pages, 20 figures

In this work, we introduce the use of the differential geometry Frenet-Serret equations to describe a magnetic line in a pulsar magnetosphere. These equations, which need to be solved numerically, fix the magnetic line in terms of their tangent, normal, and binormal vectors at each position, given assumptions on the radius of curvature and torsion. Once the representation of the magnetic line is defined, we provide the relevant set of transformations between reference frames; the ultimate aim is to express the map of the emission directions in the star co-rotating frame. In this frame, an emission map can be directly read as a light curve seen by observers located at a certain fixed angle with respect to the rotational axis. We provide a detailed step-by-step numerical recipe to obtain the emission map for a given emission process, and give a set of simplified benchmark tests. Key to our approach is that it offers a setting to achieve an effective description of the system's geometry {\it together} with the radiation spectrum. This allows to compute multi-frequency light curves produced by a specific radiation process (and not just geometry) in the pulsar magnetosphere, and intimately relates with averaged observables such as the spectral energy distribution., Comment: Accepted for publication in MNRAS

We apply a synchro-curvature spectral emission model based on characterizing the dynamics of magnetospheric particles to fit the phase-average spectra of the most extended database for the non-thermal spectra of pulsars. We consider 36 pulsars with well-determined non-thermal spectra from X-rays to gamma-rays. The sample includes Crab and Crab twin, for which the spectra extends even to the optical/ultraviolet and infrared energies. We find that the model --with just three physical parameters and a global scaling-- can fit the observations well across eight orders of magnitude for 18 of the 36 pulsars studied. Additionally, we find a set of 8 pulsars for which the model still provides arguably good fits and another set of 10 pulsars for which the model fails in reproducing the spectra. We discuss why, propose and provide physical interpretations for a simple model extension (related to the geometry of the accelerating system with regards to the observer) that allows dealing with all such cases, ultimately providing very good fits for all pulsars. The extended model is still austere, adding only two additional parameters to the former set, of the same kind of the ones previously used. We use these fits to discuss issues going from the observed spectral origin, to the extent of the dominance of synchrotron or curvature regimes, the use of a model as predictor for searching new non-thermal pulsars starting from gamma-ray surveys, and how the model offers a setting where phase shifts between X-ray and gamma-ray light curves would naturally arise., Comment: 20 pages, 15 figures, 2 tables. Accepted for publication in MNRAS on August 14

MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at play, which is instead effectively modeled by a set of additional terms in the evolution equations, dubbed as sub-grid-scale model. Here we extend such approach, commonly used in non-relativistic/non-magnetic/incompressible fluid dynamics, applying the so-called gradient model to a general set of balance-law equations, that includes the relevant case in which a non-trivial inversion of conserved to primitive fields is needed. In particular, we focus on the relativistic compressible ideal MHD scenario, providing for the first time (and for any equation of state) all the additional sub-grid-scale terms. As an application, we consider box simulations of the relativistic Kelvin-Helmholtz instability, which is also the first mechanism responsible for the magnetic field amplification in binary neutron star mergers and cannot yet be fully captured by the finest-grid and longest simulations available. The performance of our model is numerically assessed by comparing it to the residuals arising from the filtering of high-resolution simulations. We find that the model can fit very well those residuals from resolutions a few times higher. Although the application shown here explicitly considers the Minkowski metric, it can be directly extended to general relativity, thus settling the basis to implement the gradient sub-grid model in a GRMHD binary merger. Our results suggest that this approach will be potentially able to unveil much better the small-scale dynamics achievable in full GRMHD simulations., Comment: 15 pages, 6 figures. Final version. arXiv admin note: text overlap with arXiv:1904.04099

Recent detections by the gravitational wave facilities LIGO/Virgo have opened a window to study the internal structure of neutron stars through the gravitational waves emitted during their coalescence. In this work we explore, through numerical simulations, the gravitational radiation produced by the merger of binary neutron stars with dark matter particles trapped on their interior, focusing on distinguishable imprints produced by these dark matter cores. Our results reveal the presence of a strong m = 1 mode in the waveforms during the post-merger stage, together with other relevant features. Comparison of our results with observations might allow us to constraint the amount of dark matter in the interior of neutron star., Comment: 15 pages, 9 figures

Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales. Here we extend the sub-grid-scale gradient model, so far used in the momentum and induction equations, to account also for the unresolved scales in the energy evolution equation of a compressible ideal MHD fluid with a generic equation of state. We assess the model by considering box simulations of the turbulence triggered across a shear layer by the Kelvin-Helmholtz instability, testing cases where the small-scale dynamics cannot be fully captured by the resolution considered, such that the efficiency of the simulated dynamo effect depends on the resolution employed. This lack of numerical convergence is actually a currently common issue in several astrophysical problems, where the integral and fastest-growing-instability scales are too far apart to be fully covered numerically. We perform a-priori and a-posteriori tests of the extended gradient model. In the former, we find that, for many different initial conditions and resolutions, the gradient model outperforms other commonly used models in terms of correlation with the residuals coming from the filtering of a high-resolution run. In the second test, we show how a low-resolution run with the gradient model is able to quantitatively reproduce the evolution of the magnetic energy (the integrated value and the spectral distribution) coming from higher-resolution runs. This extension is the first step towards the implementation in relativistic magnetohydrodynamics., Comment: 25 pages, 23 figures, 2 tables. Submitted as a first revision

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., Comment: 6 pages, 5 figures

In this paper we describe QUBIC, an experiment that will observe the polarized microwave sky with a novel approach, which combines the sensitivity of state-of-the art bolometric detectors with the systematic effects control typical of interferometers. QUBIC unique features are the so-called "self-calibration", a technique that allows us to clean the measured data from instrumental effects, and its spectral imaging power, i.e. the ability to separate the signal in various sub-bands within each frequency band. QUBIC will observe the sky in two main frequency bands: 150 GHz and 220 GHz. A technological demonstrator is currently under testing and will be deployed in Argentina during 2019, while the final instrument is expected to be installed during 2020., Comment: Proceedings of the 2018 ICNFP conference, Crete. Published by Universe arXiv admin note: text overlap with arXiv:1801.03730

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., Comment: 23 pages, 13 figures. In press

The Cosmological Problem is considered in a five-dimensional (bulk) manifold with two time coordinates, obeying vacuum Einstein field equations. The evolution formalism is used there, in order to get a simple form of the resulting constraints. In the spatially flat case, this approach allows to find out the general solution, which happens to consist in a single metric. All the embedded Friedmann-Robertson-Walker (FRW) metrics can be obtained from this 'mother' metric ('M-metric') in the bulk, by projecting onto different four-dimensional hypersurfaces (branes). Having a time plane in the bulk allows to devise the specific curve which will be kept as the physical time coordinate in the brane. This method is applied for identifying FRW regular solutions, evolving from the infinite past (no Big Bang), even with an asymptotic initial state with non-zero radius (emergent universes). Explicit counter-examples are provided, showing that not every spatially-flat FRW metric can actually be embedded in a 3+2 bulk manifold. This implies that the extension of the Campbell theorem to the General Relativity case works only in its weaker form in this case, requiring as an extra assumption that the constraint equations hold at least in a single 4D hypersurface., Comment: Revised version, with a new section which provides some counter-examples showing that the basic assumption in Campbell's theorem is not always fulfilled when the extra dimension is time-like

The tremendous challenge of comparing our theoretical models with the gravitational-wave observations in the new era of multimessenger astronomy requires accurate and fast numerical simulations of complicated physical systems described by the Einstein and the matter equations. These requirements can only be satisfied if the simulations can be parallelized efficiently on a large number of processors and advanced computational strategies are adopted. To achieve this goal we have developed Simflowny, an open platform for scientific dynamical models which automatically generates parallel code for different simulation frameworks, allowing the use of HPC infrastructures to non-specialist scientists. One of these frameworks is SAMRAI, a mature patch-based structured adaptive mesh refinement infrastructure, capable of reaching exascale in some specific problems. Here we present the numerical techniques that we have implemented on this framework by using Simflowny in order to perform fast, efficient, accurate and highly-scalable simulations. These techniques involve high-order schemes for smooth and non-smooth solutions, Adaptive Mesh Refinement with arbitrary resolution ratios and an optimal strategy for the sub-cycling in time. We validate the automatically generated codes for the SAMRAI infrastructure with some simple test examples (i.e., wave equation and Newtonian MHD) and finally with the Einstein equations., Comment: 23 pages, 17 figures

$\gamma$-ray radiation from pulsars is usually thought to be mostly produced by the synchro-curvature losses of accelerated particles. Here we present a systematic study of all currently reported, good-quality Fermi-LAT pulsar spectral data. We do so by applying a model which follows the particle dynamics and consistently computes the emission of synchro-curvature radiation. By fitting observational data on a case by case basis, we are able to obtain constraints about the parallel electric field, the typical lengthscale over which particles emit the bulk of the detected radiation, and the number of involved particles. The model copes well with data of several dozens of millisecond and young pulsars. By correlating the inferred model parameters with the observed timing properties, some trends are discovered. First, a non-negligible part of the radiation comes from the loss of perpendicular momentum soon after pair creation. Second, the electric field strongly correlates with both the inverse of the emission lengthscale and the magnetic field at light cylinder, thus ruling out models with high-energy photon production close to the surface. These correlations unify young and millisecond pulsars under the same physical scenario, and predict that magnetars are intrinsically $\gamma$-ray quiet via syncrhro-curvature processes, since magnetospheric particles are not accelerated enough to emit a detectable $\gamma$-ray flux., Comment: Accepted for publication in MNRAS main journal. 25 pages, 4 tables, 16 figures

We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6x10^12 G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. Nonetheless, the modelling allows us to place constraints on the system geometry (i.e. the angles $\psi$ and $\xi$ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte-Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to $65\deg < \psi+\xi < 125\deg$ or $235\deg < \psi+\xi <295\deg$; (b) the temperature contrast between the poles and the equator must be at least a factor of $\sim6$, and (c) the size of the hottest region ranges between 0.2-0.7 km (including uncertainties on the source distance). Last, we interpret our findings within the context of internal and external heating models., Comment: 13 pages, 10 figures. Accepted for publication in MNRAS

$\gamma$-ray spectra of pulsars have been mostly studied in a phenomenological way, by fitting them to a cut-off power-law function. Here, we analyze a model where pulsed emission comes from synchro-curvature processes in a gap. We calculate the variation of kinetic energy of magnetospheric particles along the gap and the associated radiated spectra, considering an effective particle distribution. We fit the phase-averaged and phase-resolved Fermi-LAT spectra of the three brightest gamma-ray pulsars: Geminga, Crab, and Vela, and constrain the three free parameters we leave free in the model. Our best-fit models well reproduce the observed data, apart from residuals above a few GeV in some cases, range for which the inverse Compton scattering likely becomes the dominant mechanism. In any case, the flat slope at low-energy (< GeV) seen by Fermi-LAT both in the phase-averaged and phase-resolved spectra of most pulsars, including the ones we studied, requires that most of the detected radiation below ~ GeV is produced during the beginning of the particle trajectories, when radiation mostly come from the loss of perpendicular momentum., Comment: 11 pages, 9 figures, 2 tables. Published in MNRAS main journal

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., Comment: 12 pages, 3 figures. Proceedings of the invited talk at SEA2014, as best Spanish PhD thesis in Astronomy and Astrophysics 2013

The popular outer gap model of magnetospheric emission from pulsars has been widely applied to explain the properties observed in $\gamma$-rays. However, its quantitative predictions rely on a number of approximations and assumptions that are usually overlooked. Here we examine them, reviewing the main ingredients entering in the model, evaluating their range of uncertainties. Usually, in the quantitative applications of the model, key parameters like the radius of curvature and the energies of the interacting photons are taken to be a fixed, single value. Instead, here we explore their realistic ranges, and the impact of these on the consistency of the model itself. We conclude that the popular evaluation of the trans-field size of the gap as a function of period and period derivative, is unreliable and affected by a huge dispersion. Last, the exploration of the possible values for the radius of curvature, the local magnetic field and other quantities deserve more attention for quantitative applications of the outer gap model, like the calculation of $\gamma$-ray spectra, which is the subject of an accompanying paper., Comment: 19 pages, 5 figures, 3 tables. Accepted for pubblication in MNRAS main journal. References updated

One of the most important predictions of any gap model for pulsar magnetospheres is the predicted $\gamma$-ray spectra. In the outer gap model, the properties of the synchro-curvature radiation are sensitive to many parameters, whose realistic ranges have been studied in detail in an accompanying paper. There we demonstrated that the uncertainty in the radius of curvature, the magnetic field geometry, and the X-ray surface flux may affect by orders of magnitude the predicted flux and spectral peak in the $\gamma$-ray regime. Here, we present a systematic, numerical study of the impact of the different parameters on the particle dynamics along the gap and calculate the emitted synchro-curvature radiation along the trajectory. By integrating the emitted radiation along the gap and convolving it with a parametrized particle distribution, we discuss how the comparison with the wealth of {\em Fermi}-LAT data can be used to constrain the applicability of the model. The resulting spectra show very different energy peaks, fluxes and shapes, qualitatively matching the great variety of the observed {\em Fermi}-LAT pulsars. In particular, if we see a large fraction of photons emitted from the initial part of the trajectory, we show that the spectra will be flatter at the low-energy {\it Fermi}-LAT regime (100 MeV -- 1 GeV). This provides a solution for such observed flat spectra, while still maintain synchro-curvature radiation as the origin of these photons., Comment: 9 pages, 6 figures, 1 table. Accepted for pubblication in MNRAS main journal. References updated

We consider the fundamental problem of charged particles moving along and around a curved magnetic field line, revising the synchro-curvature radiation formulae introduced by Cheng and Zhang (1996). We provide more compact expressions to evaluate the spectrum emitted by a single particle, identifying the key parameter that controls the transition between the curvature-dominated and the synchrotron-dominated regime. This parameter depends on the local radius of curvature of the magnetic field line, the gyration radius, and the pitch angle. We numerically solve the equations of motion for the emitting particle by considering self-consistently the radiative losses, and provide the radiated spectrum produced by a particle when an electric acceleration is balanced by its radiative losses, as it is assumed to happen in the outer gaps of pulsar's magnetospheres. We compute the average spectrum radiated throughout the particle trajectory finding that the slope of the spectrum before the peak depends on the location and size of the emission region. We show how this effect could then lead to a variety of synchro-curvature spectra. Our results reinforce the idea that the purely synchrotron or curvature losses are, in general, inadequate to describe the radiative reaction on the particle motion, and the spectrum of emitted photons. Finally, we discuss the applicability of these calculations to different astrophysical scenarios., Comment: 9 pages, 5 figures, 2 tables. Accepted for publication in MNRAS main journal. References updated

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 $\sigma (\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 $\gamma$-rays)., Comment: 10 pages, 9 figures, submitted and accepted to MNRAS, comments welcome

The thermal X-ray spectra of several isolated neutron stars display deviations from a pure blackbody. The accurate physical interpretation of these spectral features bears profound implications for our understanding of the atmospheric composition, magnetic field strength and topology, and equation of state of dense matter. With specific details varying from source to source, common explanations for the features have ranged from atomic transitions in the magnetized atmospheres or condensed surface, to cyclotron lines generated in a hot ionized layer near the surface. Here we quantitatively evaluate the X-ray spectral distortions induced by inhomogeneous temperature distributions of the neutron star surface. To this aim, we explore several surface temperature distributions, we simulate their corresponding general relativistic X-ray spectra (assuming an isotropic, blackbody emission), and fit the latter with a single blackbody model. We find that, in some cases, the presence of a spurious 'spectral line' is required at a high significance level in order to obtain statistically acceptable fits, with central energy and equivalent width similar to the values typically observed. We also perform a fit to a specific object, RX J0806.4-4123, finding several surface temperature distributions able to model the observed spectrum. The explored effect is unlikely to work in all sources with detected lines, but in some cases it can indeed be responsible for the appearance of such lines. Our results enforce the idea that surface temperature anisotropy can be an important factor that should be considered and explored also in combination with more sophisticated emission models like atmospheres., Comment: 11 pages, 7 figures; accepted for publication in MNRAS

We present a long-term phase-coherent timing analysis and pulse-phase resolved spectroscopy for the two outbursts observed from the transient anomalous X-ray pulsar CXOU J164710.2-45521. For the first outburst we used 11 Chandra and XMM-Newton observations between September 2006 to August 2009, the longest baseline yet for this source. We obtain a coherent timing solution with $P=10.61065583(4)$ s, $\dot{P} = 9.72(1) \times 10^{-13}\;$s s$^{-1}$ and $\ddot{P} = -1.05(5)\times10^{-20}\; $s s$^{-2}$. Under the standard assumptions this implies a surface dipolar magnetic field of $\sim 10^{14}$ G, confirming this source as a standard-$B$ magnetar. We also study the evolution of the pulse profile (shape, intensity and pulsed fraction) as a function of time and energy. Using the phase-coherent timing solution we perform a phase-resolved spectroscopy analysis, following the spectral evolution of pulse-phase features, which hints at the physical processes taking place on the star. The results are discussed from the perspective of magneto-thermal evolution models and the untwisting magnetosphere model. Finally, we present similar analysis for the second, less intense, 2011 outburst. For the timing analysis we used Swift data together with 2 XMM-Newton and Chandra pointings. The results inferred for both outbursts are compared and briefly discussed in a more general framework., Comment: 13 pages, 14 figures

This work aims at studying how magnetic fields affect the observational properties and the long-term evolution of isolated neutron stars, which are the strongest magnets in the universe. The extreme physical conditions met inside these astronomical sources complicate their theoretical study, but, thanks to the increasing wealth of radio and X-ray data, great advances have been made over the last years. A neutron star is surrounded by magnetized plasma, the so-called magnetosphere. Modeling its global configuration is important to understand the observational properties of the most magnetized neutron stars, magnetars. On the other hand, magnetic fields in the interior are thought to evolve on long time-scales, from thousands to millions of years. The magnetic evolution is coupled to the thermal one, which has been the subject of study in the last decades. An important part of this thesis presents the state-of-the-art of the magneto-thermal evolution models of neutron stars during the first million of years, studied by means of detailed simulations. The numerical code here described is the first one to consistently consider the coupling of magnetic field and temperature, with the inclusion of both the Ohmic dissipation and the Hall drift in the crust., Comment: PhD thesis defended in University of Alicante, 20 September 2013

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50-60% level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated "observed" spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ~ 1-2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages, and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements., Comment: 12 pages, 7 figures, accepted to MNRAS

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radio-pulsars, thus confirming the idea that magnetic fields are involved in their X-ray emission. Here we present the results of 2D simulations of the fully-coupled evolution of temperature and magnetic field in neutron stars, including the state-of-the-art kinetic coefficients and, for the first time, the important effect of the Hall term. After gathering and thoroughly re-analysing in a consistent way all the best available data on isolated, thermally emitting neutron stars, we compare our theoretical models to a data sample of 40 sources. We find that our evolutionary models can explain the phenomenological diversity of magnetars, high-B radio-pulsars, and isolated nearby neutron stars by only varying their initial magnetic field, mass and envelope composition. Nearly all sources appear to follow the expectations of the standard theoretical models. Finally, we discuss the expected outburst rates and the evolutionary links between different classes. Our results constitute a major step towards the grand unification of the isolated neutron star zoo., Comment: 21 pages, 11 figures, MNRAS accepted