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 magnetar Swift ,J1818.0-1607 was discovered in March 2020 when Swift detected a 9 ms hard X-ray burst and a long-lived outburst. Prompt X-ray observations revealed a spin period of 1.36 s, soon confirmed by the discovery of radio pulsations. We report here on the analysis of the Swift burst and follow-up X-ray and radio observations. The burst average luminosity was $L_{\rm burst} \sim2\times 10^{39}$ erg/s (at 4.8 kpc). Simultaneous observations with XMM-Newton and NuSTAR three days after the burst provided a source spectrum well fit by an absorbed blackbody ($N_{\rm H} = (1.13\pm0.03) \times 10^{23}$ cm$^{-2}$ and $kT = 1.16\pm0.03$ keV) plus a power-law ($\Gamma=0.0\pm1.3$) in the 1-20 keV band, with a luminosity of $\sim$$8\times10^{34}$ erg/s, dominated by the blackbody emission. From our timing analysis, we derive a dipolar magnetic field $B \sim 7\times10^{14}$ G, spin-down luminosity $\dot{E}_{\rm rot} \sim 1.4\times10^{36}$ erg/s and characteristic age of 240 yr, the shortest currently known. Archival observations led to an upper limit on the quiescent luminosity $<$$5.5\times10^{33}$ erg/s, lower than the value expected from magnetar cooling models at the source characteristic age. A 1 hr radio observation with the Sardinia Radio Telescope taken about 1 week after the X-ray burst detected a number of strong and short radio pulses at 1.5 GHz, in addition to regular pulsed emission; they were emitted at an average rate 0.9 min$^{-1}$ and accounted for $\sim$50% of the total pulsed radio fluence. We conclude that Swift ,J1818.0-1607 is a peculiar magnetar belonging to the small, diverse group of young neutron stars with properties straddling those of rotationally and magnetically powered pulsars. Future observations will make a better estimation of the age possible by measuring the spin-down rate in quiescence.

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.

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.

The Q & U Bolometric Interferometer for Cosmology, QUBIC, is an innovative experiment designed to measure the polarization of the Cosmic Microwave Background and in particular the signature left therein by the inflationary expansion of the Universe. The expected signal is extremely faint, thus extreme sensitivity and systematic control are necessary in order to attempt this measurement. QUBIC addresses these requirements using an innovative approach combining the sensitivity of Transition Edge Sensor cryogenic bolometers, with the deep control of systematics characteristic of interferometers. This makes QUBIC unique with respect to others classical imagers experiments devoted to the CMB polarization. In this contribution we report a description of the QUBIC instrument including recent achievements and the demonstration of the bolometric interferometry performed in lab. QUBIC will be deployed at the observation site in Alto Chorrillos, in Argentina at the end of 2019.

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.

QUBIC (Q \& U Bolometric Interferometer for Cosmology) is an international ground-based experiment dedicated in the measurement of the polarized fluctuations of the Cosmic Microwave Background (CMB). It is based on bolometric interferometry, an original detection technique which combine the immunity to systematic effects of an interferometer with the sensitivity of low temperature incoherent detectors. QUBIC will be deployed in Argentina, at the Alto Chorrillos mountain site near San Antonio de los Cobres, in the Salta province. The QUBIC detection chain consists in 2048 NbSi Transition Edge Sensors (TESs) cooled to 350mK.The voltage-biased TESs are read out with Time Domain Multiplexing based on Superconducting QUantum Interference Devices (SQUIDs) at 1 K and a novel SiGe Application-Specific Integrated Circuit (ASIC) at 60 K allowing to reach an unprecedented multiplexing (MUX) factor equal to 128. The QUBIC experiment is currently being characterized in the lab with a reduced number of detectors before upgrading to the full instrument. I will present the last results of this characterization phase with a focus on the detectors and readout system.

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.

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.

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.

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.

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.

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.

In this paper we describe QUBIC, an experiment that takes up the challenge posed by the detection of primordial gravitational waves with a novel approach, that combines the sensitivity of state-of-the art bolometric detectors with the systematic effects control typical of interferometers. The so-called "self-calibration" is a technique deeply rooted in the interferometric nature of the instrument and allows us to clean the measured data from instrumental effects. The first module of QUBIC is a dual band instrument (150 GHz and 220 GHz) that will be deployed in Argentina during the Fall 2018.

QUBIC is an instrument aiming at measuring the B mode polarisation anisotropies at medium scales angular scales (30-200 multipoles). The search for the primordial CMB B-mode polarization signal is challenging, because of many difficulties: smallness of the expected signal, instrumental systematics that could possibly induce polarization leakage from the large E signal into B, brighter than anticipated polarized foregrounds (dust) reducing to zero the initial hope of finding sky regions clean enough to have a direct primordial B-modes observation. The QUBIC instrument is designed to address all aspects of this challenge with a novel kind of instrument, a Bolometric Interferometer, combining the background-limited sensitivity of Transition-Edge-Sensors and the control of systematics allowed by the observation of interference fringe patterns, while operating at two frequencies to disentangle polarized foregrounds from primordial B mode polarization. Its characteristics are described in details in this Technological Design Report.

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.

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.

$\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.

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.

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.

$\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.

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.

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.

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.

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.

Diffuse radio emission was detected around the soft gamma-ray repeater SGR 1806-20, after its 2004 powerful giant flare. We study the possible extended X-ray emission at small scales around SGR 1806-20, in two observations by the High Resolution Camera Spectrometer (HRC-S) on board of the Chandra X-ray Observatory: in 2005, 115 days after the giant flare, and in 2013, during quiescence. We compare the radial profiles extracted from data images and PSF simulations, carefully considering various issues related with the uncertain calibration of the HRC PSF at sub-arcsecond scales. We do not see statistically significant excesses pointing to an extended emission on scales of arcseconds. As a consequence, SGR 1806-20 is compatible with being point-like in X-rays, months after the giant flare, as well as in quiescence.

According to the partially screened gap scenario, an efficient electron-positron pair creation, a general precondition of radio-pulsar activity, relies on the existence of magnetic spots, i.e., local concentrations of strong and small scale magnetic field structures at the surface of neutron stars. They have a strong impact on the surface temperature, which is potentially observable. Here we reinforce the idea that such magnetic spots can be formed by extracting magnetic energy from the toroidal field that resides in deep crustal layers, via Hall drift. We study and discuss the magneto-thermal evolution of qualitatively different neutron star models and initial magnetic field configurations that lead to the creation of magnetic spots. We find that magnetic spots can be created on a timescale of $10^4$ years with magnetic field strengths $\gtrsim 5\times 10^{13}$ G, provided almost the whole magnetic energy is stored in its toroidal component, and that the conductivity in the inner crust is not too large. The lifetime of the magnetic spots is at least $\sim$one million of years, being longer if the initial field permeates both core and crust.

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).

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.

The lack of X-ray pulsars with spin periods > 12 s raises the question of where the population of evolved high magnetic field neutron stars has gone. Unlike canonical radio pulsars, X-ray pulsars are not subject to physical limits to the emission mechanism nor observational biases against the detection of sources with longer periods. Here we show that a highly resistive layer in the innermost part of the crust of neutron stars naturally limits the spin period to a maximum value of about 10-20 s. This highly resistive layer is expected if the inner crust is amorphous and heterogeneous in nuclear charge, possibly due to the existence of a nuclear pasta phase. Our findings suggest that the maximum period of isolated X-ray pulsars can be the first observational evidence of an amorphous inner crust, which properties can be further constrained by future X-ray timing missions combined with more detailed models.

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.

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.

Central Compact Objects are peculiar young neutron stars, with very low external magnetic fields combined with high fluxes in the X-ray band and surface temperature anisotropies. However, in their crust the magnetic field can be strong, result of its burial during a short post-supernova hypercritical accretion episode. The implications of this latter scenario for the temperature anisotropy, pulsed fraction and luminosity are discussed.

CCOs are X-ray sources lying close the center of supernova remnants, with inferred values of the surface magnetic fields significantly lower (less than about 1e11 G) than those of standard pulsars. In this paper, we revise the hidden magnetic field scenario, presenting the first 2D simulations of the submergence and reemergence of the magnetic field in the crust of a neutron star. A post-supernova accretion stage of about 1e-4-1e-3 solar masses over a vast region of the surface is required to bury the magnetic field into the inner crust. When accretion stops, the field reemerges on a typical timescale of 1-100 kyr, depending on the submergence conditions. After this stage, the surface magnetic field is restored close to its birth values. A possible observable consequence of the hidden magnetic field is the anisotropy of the surface temperature distribution, in agreement with observations of several of these sources. We conclude that the hidden magnetic field model is viable as alternative to the anti-magnetar scenario, and it could provide the missing link between CCOs and the other classes of isolated neutron stars.

Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.

Nowadays, the analysis of the X-ray spectra of magnetically powered neutron stars or magnetars is one of the most valuable tools to gain insight into the physical processes occurring in their interiors and magnetospheres. In particular, the magnetospheric plasma leaves a strong imprint on the observed X-ray spectrum by means of Compton up-scattering of the thermal radiation coming from the star surface. Motivated by the increased quality of the observational data, much theoretical work has been devoted to develop Monte Carlo (MC) codes that incorporate the effects of resonant Compton scattering in the modeling of radiative transfer of photons through the magnetosphere. The two key ingredients in this simulations are the kinetic plasma properties and the magnetic field (MF) configuration. The MF geometry is expected to be complex, but up to now only mathematically simple solutions (self-similar solutions) have been employed. In this work, we discuss the effects of new, more realistic, MF geometries on synthetic spectra. We use new force-free solutions in a previously developed MC code to assess the influence of MF geometry on the emerging spectra. Our main result is that the shape of the final spectrum is mostly sensitive to uncertain parameters of the magnetospheric plasma, but the MF geometry plays an important role on the angle-dependence of the spectra.

The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.

The Anti Coincidence Shield (ACS) of the INTEGRAL SPI instrument provides an excellent sensitivity for the detection of Gamma Ray Bursts (GRBs) above ~ 75keV, but no directional and energy information is available. We studied the ACS response by using GRBs with known localizations and good spectral information derived by other satellites. We derived a count rate to flux conversion factor for different energy ranges and studied its dependence on the GRB direction and spectral hardness. For a typical GRB spectrum, we found that 1 ACS count corresponds on average to ~ 1E-10 erg/cm^2 in the 75keV-1MeV range, for directions orthogonal to the satellite pointing axis. This is broadly consistent with the ACS effective area derived from the Monte Carlo simulations, but there is some indication that the latter slightly overestimates the ACS sensitivity, especially for directions close to the instrument axis.