The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has opened interesting debates about the number, age and nature of pulsars expected in the Galactic center region. In this work, we present extensive X-ray monitoring of the outburst of SGR J1745-2900 using the Chandra X-ray Observatory, the only instrument with the spatial resolution to distinguish the magnetar from the supermassive black hole (2.4" angular distance). It was monitored from its outburst onset in April 2013 until August 2019, collecting more than fifty Chandra observations for a total of more than 2.3 Ms of data. Soon after the outburst onset, the magnetar emission settled onto a purely thermal emission state that cooled from a temperature of about 0.9 to 0.6 keV over 6 years. The pulsar timing properties showed at least two changes in the period derivative, increasing by a factor of about 4 during the outburst decay. We find that the long-term properties of this outburst challenge current models for the magnetar outbursts.

It is generally accepted that the non-linear, dynamical evolution of magnetic fields in the interior of neutron stars plays a key role in the explanation of the observed phenomenology. Understanding the transfer of energy between toroidal and poloidal components, or between different scales, is of particular relevance. In this letter, we present the first 3D simulations of the Hall instability in a neutron star crust, confirming its existence for typical magnetar conditions. We confront our results to estimates obtained by a linear perturbation analysis, which discards any interpretation as numerical instabilities and confirms its physical origin. Interestingly, the Hall instability creates locally strong magnetic structures that occasionally can make the crust yield to the magnetic stresses and generates coronal loops, similarly as solar coronal loops find their way out through the photosphere. This supports the viability of the mechanism, which has been proposed to explain magnetar outbursts.

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.

We present three-dimensional force-free electrodynamics simulations of magnetar magnetospheres that demonstrate the instability of certain degenerate, high energy equilibrium solutions of the Grad-Shafranov equation. This result indicates the existence of an unstable branch of twisted magnetospheric solutions and allows to formulate an instability criterion. The rearrangement of magnetic field lines as a consequence of this instability triggers the dissipation of up to 30% of the magnetospheric energy on a thin layer above the magnetar surface. During this process, we predict an increase of the mechanical stresses onto the stellar crust, which can potentially result in a global mechanical failure of a significant fraction of it. We find that the estimated energy release and the emission properties are compatible with the observed giant flare events. The newly identified instability is a candidate for recurrent energy dissipation, which could explain part of the phenomenology observed in magnetars.

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

After two major outbursts in 2006 and 2011, on 2017 May 16 the magnetar CXOU J164710.2-455216, hosted within the massive star cluster Westerlund I, emitted a short (20 ms) burst, which marked the onset of a new active phase. We started a long-term monitoring campaign with Swift (45 observations), Chandra (5 observations) and NuSTAR (4 observations) from the activation until 2018 April. During the campaign, Swift BAT registered the occurrence of multiple bursts, accompanied by two other enhancements of the X-ray persistent flux. The long time span covered by our observations allowed us to study the spectral as well as the timing evolution of the source. After 11 months since the 2017 May outburst onset, the observed flux was about 15 times higher than its historical minimum level and a factor of 3 higher than the level reached after the 2006 outburst. This suggests that the crust has not fully relaxed to the quiescent level, or that the source quiescent level has changed following the multiple outburst activities in the past 10 years or so. This is another case of multiple outbursts from the same source on a yearly time scale, a somehow recently discovered behaviour in magnetars.

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.

We study the long-term quasi-steady evolution of the force-free magnetosphere of a magnetar coupled to its internal magnetic field. We find that magnetospheric currents can be maintained on long timescales of the order of thousands of years. Meanwhile, the energy, helicity and twist stored in the magnetosphere all gradually increase over the course of this evolution, until a critical point is reached, beyond which a force-free magnetosphere cannot be constructed. At this point, some large-scale magnetospheric rearrangement, possibly resulting in an outburst or a flare, must occur, releasing a large fraction of the stored energy, helicity and twist. After that, the quasi-steady evolution should continue in a similar manner from the new initial conditions. The timescale for reaching this critical point depends on the overall magnetic field strength and on the relative fraction of the toroidal field. The energy stored in the force-free magnetosphere is found to be up to $\sim 30\%$ larger than the corresponding vacuum energy. This implies that for a $10^{14}$ G field at the pole, the energy budget available for fast magnetospheric events is of the order of a few $10^{44}$ erg. The spindown rate is estimated to increase by up to $\sim 60\%$, since the dipole content in the magnetosphere is enhanced by the currents present there. A rough estimate of the braking index $n$ reveals that it is systematically $n < 3$ for the most part of the evolution, consistent with actual measurements for pulsars and early estimates for several magnetars.

We present a detailed analysis of the properties of twisted, force-free magnetospheres of non-rotating neutron stars, which are of interest in the modelling of magnetar properties and evolution. In our models the magnetic field smoothly matches to a current-free (vacuum) solution at some large external radius, and they are specifically built to avoid pathological surface currents at any of the interfaces. By exploring a large range of parameters, we find a few remarkable general trends. We find that the total dipolar moment can be increased by up to $40\%$ with respect to a vacuum model with the same surface magnetic field, due to the contribution of magnetospheric currents to the global magnetic field. Thus, estimates of the surface magnetic field based on the large-scale dipolar braking torque are slightly overestimating the surface value by the same amount. Consistently, there is a moderate increase in the total energy of the model with respect to the vacuum solution of up to $25\%$, which would be the available energy budget in the event of a fast, global magnetospheric reorganization commonly associated with magnetar flares. We have also found the interesting result of the existence of a critical twist ($\varphi_{\rm max} \lesssim 1.5$ rad), beyond which we cannot find any more numerical solutions. Combining the models considered in this paper with the evolution of the interior of neutron stars will allow us to study the influence of the magnetosphere on the long-term magnetic, thermal, and rotational evolution.

We extend our previous study of equilibrium solutions of non-rotating force-free magnetospheres of neutron stars. We show that multiple solutions exist for the same sets of parameters, implying that the solutions are degenerate. We are able to obtain configurations with disconnected field lines, however, in nearly all cases these correspond to degenerate higher energy solutions. We carry out a wide parametric search in order to understand the properties of the solutions. We confirm our previous results that the lower energy solutions have up to $\sim 25\%$ more energy than the vacuum case, helicity of the order of $\sim 5$ (in some defined units), maximum twist of $\sim 1.5$ rad, and a dipole strength that is up to $\sim 40\%$ larger than the vacuum dipole. Including the degenerate higher energy solutions allows for larger theoretical limits of up to $\sim 80\%$ more energy with respect to the vacuum case, helicity of the order of $\sim 8$, and a dipole strength that can now be up to four times that of the vacuum dipole, while the twist can be significantly larger and even diverge for configurations with disconnected domains. The higher energy solutions are probably unstable, therefore, it is unlikely that such magnetospheres exist under normal conditions in magnetars and high magnetic field pulsars.

We study the coupling of the force-free magnetosphere to the long-term internal evolution of a magnetar. We allow the relation between the poloidal and toroidal stream functions - that characterizes the magnetosphere - to evolve freely without constraining its particular form. We find that, on time-scales of the order of kyr, the energy stored in the magnetosphere gradually increases, as the toroidal region grows and the field lines expand outwards. This continues until a critical point is reached beyond which force-free solutions for the magnetosphere can no longer be constructed, likely leading to some large-scale magnetospheric reorganization. The energy budget available for such events can be as high as several $10^{45}\,$erg for fields of $10^{14}\,$G. Subsequently, starting from the new initial conditions, the evolution proceeds in a similar manner. The time-scale to reach the critical point scales inversely with the magnetic field amplitude. Allowing currents to pass through the last few meters below the surface, where the magnetic diffusivity is orders of magnitude larger than in the crust, should give rise to a considerable amount of energy deposition through Joule heating. We estimate that the effective surface temperature could increase locally from $\sim 0.1\,$keV to $\sim 0.3 - 0.6\,$keV, in good agreement with observations. Similarly, the power input from the interior into the magnetosphere could be as high as $10^{35} - 10^{36}\,$erg/s, which is consistent with peak luminosities observed during magnetar outbursts. Therefore, a detailed treatment of currents flowing through the envelope may be needed to explain the thermal properties of magnetars.

We report on a new NuSTAR observation and on the ongoing Swift XRT monitoring campaign of the peculiar source 1E 161348-5055, located at the centre of the supernova remnant RCW 103, which is recovering from its last outburst in June 2016. The X-ray spectrum at the epoch of the NuSTAR observation can be described by either two absorbed blackbodies ($kT_{BB_1}$ ~ 0.5 keV, $kT_{BB_2}$ ~ 1.2 keV) or an absorbed blackbody plus a power law ($kT_{BB_1}$ ~ 0.6 keV, $\Gamma$ ~ 3.9). The observed flux was ~ 9 $\times$ 10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, ~ 3 times lower than what observed at the outburst onset, but about one order of magnitude higher than the historical quiescent level. A periodic modulation was detected at the known 6.67 hr periodicity. The spectral decomposition and evolution along the outburst decay are consistent with 1E 161348-5055 being a magnetar, the slowest ever detected.

The mass of the stellar precursor of supernova (SN) 1987A and the burst of neutrinos observed at the moment of the explosion are consistent with the core-collapse formation of a neutron star. However, no compelling evidence for the presence in SN 1987A of a compact object of any kind has been found yet in any band of the electromagnetic spectrum, prompting questions on whether the neutron star survived and, if it did, on its properties. Starting from the analysis of recent Chandra observations, here we appraise the current observational situation. We derived limits on the X-ray luminosity of a compact object with a nonthermal, Crab-pulsar-like spectrum of the order of $\approx$$(1$-$5)\times10^{35}$ erg s$^{-1}$, corresponding to limits on the rotational energy loss of a possible X-ray pulsar in SN 1987A of $\approx$$(0.5$-$1.5)\times10^{38}$ erg s$^{-1}$. However, a much brighter X-ray source cannot be excluded if, as is likely, it is enshrouded in a cloud of absorbing matter with metallicity similar to that expected in the outer layers of a massive star towards the end of its life. We found that other limits obtained from various arguments and observations in other energy ranges either are unbinding or allow a similar maximum luminosity of the order of $\approx$$10^{35}$ erg s$^{-1}$. We conclude that while a pulsar alike the one in the Crab Nebula in both luminosity and spectrum is hardly compatible with the observations, there is ample space for an `ordinary' X-ray-emitting young neutron star, born with normal initial spin period, temperature and magnetic field, to be hiding inside the evolving remnant of SN 1987A.

We present the results of the systematic study of all magnetar outbursts observed to date, through a reanalysis of data acquired in about 1100 X-ray observations. We track the temporal evolution of the outbursts soft X-ray spectral properties and the luminosities of the single spectral components as well as of the total emission. We model empirically all outburst light curves, and estimate the characteristic decay time-scales as well as the energetics involved. We investigate the link between different parameters (e.g. the luminosity at the peak of the outburst and in quiescence, the maximum luminosity increase, the decay time-scale and energy of the outburst, the neutron star surface dipolar magnetic field and characteristic age, etc.), and unveil several correlations among these quantities. We discuss our results in the context of the internal crustal heating and twisted bundle models for magnetar outbursts. This study is complemented by the Magnetar Outburst Online Catalogue (http://magnetars.ice.csic.es), an interactive data base where the user can plot any combination of the parameters derived in this work, and download all data.

Observations of thermal radiation from neutron stars can potentially provide information about the states of supranuclear matter in the interiors of these stars with the aid of the theory of neutron-star thermal evolution. We review the basics of this theory for isolated neutron stars with strong magnetic fields, including most relevant thermodynamic and kinetic properties in the stellar core, crust, and blanketing envelopes.

In a core-collapse supernova, a huge amount of energy is released in the Kelvin-Helmholtz phase subsequent to the explosion, when the proto-neutron star cools and deleptonizes as it loses neutrinos. Most of this energy is emitted through neutrinos, but a fraction of it can be released through gravitational waves. We model the evolution of a proto-neutron star in the Kelvin-Helmholtz phase using a general relativistic numerical code, and a recently proposed finite temperature, many-body equation of state; from this we consistently compute the diffusion coefficients driving the evolution. To include the many-body equation of state, we develop a new fitting formula for the high density baryon free energy at finite temperature and intermediate proton fraction. We estimate the emitted neutrino signal, assessing its detectability by present terrestrial detectors, and we determine the frequencies and damping times of the quasi-normal modes which would characterize the gravitational wave signal emitted in this stage.

We report on 3.5 years of Chandra monitoring of the Galactic Centre magnetar SGR J1745-2900 since its outburst onset in April 2013. The magnetar spin-down has shown at least two episodes of period derivative increases so far, and it has slowed down regularly in the past year or so. We observed a slightly increasing trend in the time evolution of the pulsed fraction, up to about 55 per cent in the most recent observations. SGR J1745-2900 has not reached the quiescent level yet, and so far the overall outburst evolution can be interpreted in terms of a cooling hot region on the star surface. We discuss possible scenarios, showing in particular how the presence of a shrinking hot spot in this source is hardly reconcilable with internal crustal cooling and favors the untwisting bundle model for this outburst. Moreover, we also show how the emission from a single uniform hot spot is incompatible with the observed pulsed fraction evolution for any pair of viewing angles, suggesting an anisotropic emission pattern.

We revisit the various approximations employed to study the long-term evolution of the magnetic field in neutron star cores and discuss their limitations and possible improvements. A recent controversy on the correct form of the induction equation and the relevant evolution timescale in superconducting neutron star cores is addressed and clarified. We show that this ambiguity in the estimation of timescales arises as a consequence of nominally large terms that appear in the induction equation, but which are, in fact, mostly irrotational. This subtlety leads to a discrepancy by many orders of magnitude when velocity fields are absent or ignored. Even when internal velocity fields are accounted for, only the solenoidal part of the electric field contributes to the induction equation, which can be substantially smaller than the irrotational part. We also argue that stationary velocity fields must be incorporated in the slow evolution of the magnetic field as the next level of approximation.

We study ambipolar diffusion in strongly magnetised neutron stars, with special focus on the effects of neutrino reaction rates and the impact of a superfluid/superconducting transition in the neutron star core. For axisymmetric magnetic field configurations, we determine the deviation from $\beta-$equilibrium induced by the magnetic force and calculate the velocity of the slow, quasi-stationary, ambipolar drift. We study the temperature dependence of the velocity pattern and clearly identify the transition to a predominantly solenoidal flow. For stars without superconducting/superfluid constituents and with a mixed poloidal-toroidal magnetic field of typical magnetar strength, we find that ambipolar diffusion proceeds fast enough to have a significant impact on the magnetic field evolution only at low core temperatures, $T \lesssim 1-2\times10^8$ K. The ambipolar diffusion timescale becomes appreciably shorter when fast neutrino reactions are present, because the possibility to balance part of the magnetic force with pressure gradients is reduced. We also find short ambipolar diffusion timescales in the case of superconducting cores for $T \lesssim 10^9$ K, due to the reduced interaction between protons and neutrons. In the most favourable scenario, with fast neutrino reactions and superconducting cores, ambipolar diffusion results in advection velocities of several km/kyr. This velocity can substantially reorganize magnetic fields in magnetar cores, in a way that can only be confirmed by dynamical simulations.

We study the torsional magneto-elastic oscillations of relativistic superfluid magnetars and explore the effects of a phase transition in the crust-core interface (nuclear pasta) which results in a weaker elastic response. Exploring various models with different extension of nuclear pasta phases, we find that the differences in the oscillation spectrum present in purely elastic modes (weak magnetic field), are smeared out with increasing strength of the magnetic field. For magnetar conditions, the main characteristic and features of models without nuclear pasta are preserved. We find in general two classes of magneto-elastic oscillations which exhibit a different oscillation pattern. For B_p < 4 x 10^{14} G, the spectrum is characterised by the turning points and edges of the continuum which are confined into the star's core, and have no constant phase. Increasing the magnetic field, we find, in addition, several magneto-elastic oscillations which reach the surface and have an angular structure similar to crustal modes. These global magneto-elastic oscillations show a constant phase and become dominant when B_p > 5 x 10^{14} G. We do not find any evidence of fundamental pure crustal modes in the low frequency range (below 200 Hz) for B_p >= 10^{14} G.

We study the evolution of the rotation rate of a proto-neutron star, born in a core-collapse supernova, in the first seconds of its life. During this phase, the star evolution can be described as a sequence of stationary configurations, which we determine by solving the neutrino transport and the stellar structure equations in general relativity. We include in our model the angular momentum loss due to neutrino emission. We find that the requirement of a rotation rate not exceeding the mass-shedding limit at the beginning of the evolution implies a strict bound on the rotation rate at later times. Moreover, assuming that the proto-neutron star is born with a finite ellipticity, we determine the emitted gravitational wave signal, and estimate its detectability by present and future ground-based interferometric detectors.

Confronting theoretical models with observations of thermal radiation emitted by neutron stars is one of the most important ways to understand the properties of both, superdense matter in the interiors of the neutron stars and dense magnetized plasmas in their outer layers. Here we review the theory of thermal emission from the surface layers of strongly magnetized neutron stars, and the main properties of the observational data. In particular, we focus on the nearby sources for which a clear thermal component has been detected, without being contaminated by other emission processes (magnetosphere, accretion, nebulae). We also discuss the applications of the modern theoretical models of the formation of spectra of strongly magnetized neutron stars to the observed thermally emitting objects.

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.

The observation of several neutron stars in the center of supernova remnants and with significantly lower values of the dipolar magnetic field than the average radio-pulsar population has motivated a lively debate about their formation and origin, with controversial interpretations. A possible explanation requires the slow rotation of the proto-neutron star at birth, which is unable to amplify its magnetic field to typical pulsar levels. An alternative possibility, the hidden magnetic field scenario, considers the accretion of the fallback of the supernova debris onto the neutron star as responsible for the submergence (or screening) of the field and its apparently low value. In this paper we study under which conditions the magnetic field of a neutron star can be buried into the crust due to an accreting, conducting fluid. For this purpose, we consider a spherically symmetric calculation in general relativity to estimate the balance between the incoming accretion flow and the magnetosphere. Our study analyses several models with different specific entropy, composition, and neutron star masses. The main conclusion of our work is that typical magnetic fields of a few times 1e12 G can be buried by accreting only 1e-3 - 1e-2 solar mass, a relatively modest amount of mass. In view of this result, the Central Compact Object scenario should not be considered unusual, and we predict that anomalously weak magnetic fields should be common in very young (< few kyr) neutron stars.

We study the timing and spectral properties of the low-magnetic field, transient magnetar SWIFT J1822.3-1606 as it approached quiescence. We coherently phase-connect the observations over a time-span of ~500 days since the discovery of SWIFT J1822.3-1606 following the Swift-BAT trigger on 2011 July 14, and carried out a detailed pulse phase spectroscopy along the outburst decay. We follow the spectral evolution of different pulse phase intervals and find a phase and energy-variable spectral feature, which we interpret as proton cyclotron resonant scattering of soft photon from currents circulating in a strong (10^14 G) small-scale component of the magnetic field near the neutron star surface, superimposed to the much weaker (~3 x 10^13 G) magnetic dipole field. We discuss also the implications of the pulse-resolved spectral analysis for the emission regions on the surface of the cooling magnetar.

A large fraction of Gamma Ray Bursts (GRBs) displays an X-ray plateau phase within <10^{5} s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, between January 2005 and August 2014. From the derived initial magnetic field distribution for the possible magnetars left behind by the GRBs, we study the evolution and properties of a simulated GRB-magnetar population using numerical simulations of magnetic field evolution, coupled with Monte Carlo simulations of Pulsar Population Synthesis in our Galaxy. We find that if the GRB X-ray plateaus are powered by the rotational energy of a newly formed magnetar, the current observational properties of the Galactic magnetar population are not compatible with being formed within the GRB scenario (regardless of the GRB type or rate at z=0). Direct consequences would be that we should allow the existence of magnetars and "super-magnetars" having different progenitors, and that Type Ib/c SNe related to Long GRBs form systematically neutron stars with higher initial magnetic fields. We put an upper limit of <16 "super-magnetars" formed by a GRB in our Galaxy in the past Myr (at 99% c.l.). This limit is somewhat smaller than what roughly expected from Long GRB rates, although the very large uncertainties do not allow us to draw strong conclusion in this respect.

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.

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.

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.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

Magnetars are neutron stars in which a strong magnetic field is the main energy source. About two dozens of magnetars, plus several candidates, are currently known in our Galaxy and in the Magellanic Clouds. They appear as highly variable X-ray sources and, in some cases, also as radio and/or optical pulsars. Their spin periods (2-12 s) and spin-down rates (~10^{-13}-10^{-10} s/s) indicate external dipole fields of ~10^{13-15} G, and there is evidence that even stronger magnetic fields are present inside the star and in non-dipolar magnetospheric components. Here we review the observed properties of the persistent emission from magnetars, discuss the main models proposed to explain the origin of their magnetic field and present recent developments in the study of their evolution and connection with other classes of neutron stars.

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.

We monitored the quiescent thermal emission from neutron stars in low-mass X-ray binaries after active periods of intense activity in x-rays (outbursts). The theoretical modeling of the thermal relaxation of the neutron star crust may be used to establish constraints on the crust composition and transport properties, depending on the astrophysical scenarios assumed. We numerically simulated the thermal evolution of the neutron star crust and compared them with inferred surface temperatures for five sources: MXB 1659-29, KS 1731-260, EXO 0748-676, XTE J1701-462 and IGR J17480-2446. We find that the evolution of MXB 1659-29, KS 1731-260 and EXO 0748-676 can be well described within a deep crustal cooling scenario. Conversely, we find that the other two sources can only be explained with models beyond crustal cooling. For the peculiar emission of XTE J1701-462 we propose alternative scenarios such as residual accretion during quiescence, additional heat sources in the outer crust, and/or thermal isolation of the inner crust due to a buried magnetic field. We also explain the very recent reported temperature of IGR J17480-2446 with an additional heat deposition in the outer crust from shallow sources.

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.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

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.

We report on the quiescent state of the Soft Gamma Repeater SGR 0501+4516 observed by XMM-Newton on 2009 August 30. The source exhibits an absorbed flux ~75 times lower than that measured at the peak of the 2008 outburst, and a rather soft spectrum, with the same value of the blackbody temperature observed with ROSAT back in 1992. This new observation is put into the context of all existing X-ray data since its discovery in August 2008, allowing us to complete the study of the timing and spectral evolution of the source from outburst until its quiescent state. The set of deep XMM-Newton observations performed during the few-years timescale of its outburst allows us to monitor the spectral characteristics of this magnetar as a function of its rotational period, and their evolution along these years. After the first ~10 days, the initially hot and bright surface spot progressively cooled down during the decay. We discuss the behaviour of this magnetar in the context of its simulated secular evolution, inferring a plausible dipolar field at birth of 3x10^14 G, and a current (magneto-thermal) age of ~10 kyr.

We study the outburst of the newly discovered X-ray transient 3XMM J185246.6+003317, re-analysing all available XMM-Newton, observations of the source to perform a phase-coherent timing analysis, and derive updated values of the period and period derivative. We find the source rotating at P=11.55871346(6) s (90% confidence level; at epoch MJD 54728.7) but no evidence for a period derivative in the 7 months of outburst decay spanned by the observations. This translates in a 3sigma upper limit for the period derivative of Pdot<1.4x10^{-13} s/s, which, assuming the classical magneto-dipolar braking model, gives a limit on the dipolar magnetic field of B_dip<4.1x10^{13} G . The X-ray outburst and spectral characteristics of 3XMM J185246.6+003317 confirms the identification as a magnetar, but the magnetic field upper limit we derive defines it as the third "low-B" magnetar discovered in the past three years, after SGR 0418+5729 and Swift J1822.3-1606. We have also obtained an upper limit to the quiescent luminosity (< 4x10^{33} erg/s), in line with the expectations for an old magnetar. The discovery of this new low field magnetar reaffirms the prediction of about one outburst per year from the hidden population of aged magnetars.

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.

The center of our Galaxy hosts a supermassive black hole, Sagittarius (Sgr) A*. Young, massive stars within 0.5 pc of SgrA* are evidence of an episode of intense star formation near the black hole a few Myr ago, which might have left behind a young neutron star traveling deep into SgrA*'s gravitational potential. On 2013 April 25, a short X-ray burst was observed from the direction of the Galactic center. Thanks to a series of observations with the Chandra and the Swift satellites, we pinpoint the associated magnetar at an angular distance of 2.4+/-0.3 arcsec from SgrA*, and refine the source spin period and its derivative (P=3.7635537(2) s and \dot{P} = 6.61(4)x10^{-12} s/s), confirmed by quasi simultaneous radio observations performed with the Green Bank (GBT) and Parkes antennas, which also constrain a Dispersion Measure of DM=1750+/-50 pc cm^{-3}, the highest ever observed for a radio pulsar. We have found that this X-ray source is a young magnetar at ~0.07-2 pc from SgrA*. Simulations of its possible motion around SgrA* show that it is likely (~90% probability) in a bound orbit around the black hole. The radiation front produced by the past activity from the magnetar passing through the molecular clouds surrounding the Galactic center region, might be responsible for a large fraction of the light echoes observed in the Fe fluorescence features.

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.

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.

We consider static spherically symmetric Lovelock black holes and generalize the dimensionally continued black holes in such a way that they asymptotically for large r go over to the d-dimensional Schwarzschild black hole in dS/AdS spacetime. This means that the master algebraic polynomial is not degenerate but instead its derivative is degenerate. This family of solutions contains an interesting class of pure Lovelock black holes which are the Nth order Lovelock {\Lambda}-vacuum solu- tions having the remarkable property that their thermodynamical parameters have the universal character in terms of the event horizon radius. This is in fact a characterizing property of pure Lovelock theories. We also demonstrate the universality of the asymptotic Einstein limit for the Lovelock black holes in general.

We report on the long term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729, using all the available X-ray data obtained with RXTE, SWIFT, Chandra, and XMM-Newton observations, from the discovery of the source in June 2009, up to August 2012. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: \dot{P}=4(1)x10^{-15} s/s, significant at ~3.5 sigma confidence level. This leads to a surface dipolar magnetic field of B_dip ~6x 10^{12} G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ~1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ~0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about 3 orders of magnitude: from about 1.4x10^{-11} to 1.2x10^{-14} erg/s/cm^2 . Deep radio, millimeter, optical and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ~550 kyr, and a dipolar magnetic field at birth of ~10^{14} G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such result in several other astrophysical contexts.

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.

We report on the long term X-ray monitoring with Swift, RXTE, Suzaku, Chandra, and XMM-Newton of the outburst of the newly discovered magnetar Swift J1822.3-1606 (SGR 182-1606), from the first observations soon after the detection of the short X-ray bursts which led to its discovery (July 2011), through the first stages of its outburst decay (April 2012). Our X-ray timing analysis finds the source rotating with a period of P = 8.43772016(2) s and a period derivative Pdot = 8.3(2) x 10e-14 s s-1, which entails an inferred dipolar surface magnetic field of 2.7 x 10e13 G at the equator. This measurement makes Swift J1822.3-1606 the second lowest magnetic field magnetar (after SGR 0418+5729; Rea et al. 2010). Following the flux and spectral evolution from the beginning of the outburst, we find that the flux decreased by about an order of magnitude, with a subtle softening of the spectrum, both typical of the outburst decay of magnetars. By modeling the secular thermal evolution of Swift J1822.3+1606, we find that the observed timing properties of the source, as well as its quiescent X-ray luminosity, can be reproduced if it was born with a poloidal and crustal toroidal fields of Bp ~ 1.5 x 10e14 G and Btor ~ 7 x 10e14 G, respectively, and if its current age is ~550 kyr (more details in Rea et al. 2012).

(Abridged) The rotational evolution of isolated neutron stars is dominated by the magnetic field anchored to the solid crust of the star. Assuming that the core field evolves on much longer timescales, the crustal field evolves mainly though Ohmic dissipation and the Hall drift, and it may be subject to relatively rapid changes with remarkable effects on the observed timing properties. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable consequences in the braking index, which is the most sensitive quantity to reflect small variations of the timing properties that are caused by magnetic field rearrangements. By performing axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art microphysical inputs, we find that the effect of the magnetic field evolution on the braking index can be divided into three qualitatively different stages depending on the age and the internal temperature: a first stage that may be different for standard pulsars (with n~3) or low field neutron stars that accreted fallback matter during the supernova explosion (systematically n<3); in a second stage, the evolution is governed by almost pure Ohmic field decay, and a braking index n>3 is expected; in the third stage, at late times, when the interior temperature has dropped to very low values, Hall oscillatory modes in the neutron star crust result in braking indices of high absolute value and both positive and negative signs. Models with strong (1e14 G) multipolar or toroidal components, even with a weak (~1e12 G) dipolar field are consistent with the observed trend of the timing properties.

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.

The radio pulsar models based on the existence of an inner accelerating gap located above the polar cap rely on the existence of a small scale, strong surface magnetic field $B_s$. This field exceeds the dipolar field $B_d$, responsible for the braking of the pulsar rotation, by at least one order of magnitude. Neither magnetospheric currents nor small scale field components generated during neutron star's birth can provide such field structures in old pulsars. While the former are too weak to create $B_s \gtrsim 5\times 10^{13}$G$\;\gg B_d$, the ohmic decay time of the latter is much shorter than $10^6$ years. We suggest that a large amount of magnetic energy is stored in a toroidal field component that is confined in deeper layers of the crust, where the ohmic decay time exceeds $10^7$ years. This toroidal field may be created by various processes acting early in a neutron star's life. The Hall drift is a non-linear mechanism that, due to the coupling between different components and scales, may be able to create the demanded strong, small scale, magnetic spots. Taking into account both realistic crustal microphysics and a minimal cooling scenario, we show that, in axial symmetry, these field structures are created on a Hall time scale of $10^3$-$10^4$ years. These magnetic spots can be long-lived, thereby fulfilling the pre-conditions for the appearance of the radio pulsar activity. Such magnetic structures created by the Hall drift are not static, and dynamical variations on the Hall time scale are expected in the polar cap region.

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.

We report on the long term X-ray monitoring with Swift, RXTE, Suzaku, Chandra and XMM-Newton of the outburst of the newly discovered magnetar Swift J1822.3-1606 (SGR 1822-1606), from the first observations soon after the detection of the short X-ray bursts which led to its discovery, through the first stages of its outburst decay (covering the time-span from July 2011, until end of April 2012). We also report on archival ROSAT observations which witnessed the source during its likely quiescent state, and on upper limits on Swift J1822.3-1606's radio-pulsed and optical emission during outburst, with the Green Bank Telescope (GBT) and the Gran Telescopio Canarias (GTC), respectively. Our X-ray timing analysis finds the source rotating with a period of P=8.43772016(2) s and a period derivative \dot{P}=8.3(2)x10^{-14} s s^{-1} , which entails an inferred dipolar surface magnetic field of B~2.7x10^{13} G at the equator. This measurement makes Swift J1822.3-1606 the second lowest magnetic field magnetar (after SGR 0418+5729; Rea et al. 2010). Following the flux and spectral evolution from the beginning of the outburst, we find that the flux decreased by about an order of magnitude, with a subtle softening of the spectrum, both typical of the outburst decay of magnetars. By modeling the secular thermal evolution of Swift J1822.3-1606, we find that the observed timing properties of the source, as well as its quiescent X-ray luminosity, can be reproduced if it was born with a poloidal and crustal toroidal fields of B_{p}~1.5x10^{14} G and B_{tor}~7x10^{14} G, respectively, and if its current age is ~550 kyr.

The availability of a large amount of observational data recently collected from magnetar outbursts is now calling for a complete theoretical study of outburst characteristics. In this letter (the first of a series dedicated to model magnetar outbursts), we tackle the long-standing open issue of whether or not short bursts and glitches are always connected to long-term radiative outbursts. We show that the recent detection of short bursts and glitches seemingly unconnected to outbursts is only misleading our understanding of these events. We show that, in the framework of the starquake model, neutrino emission processes in the magnetar crust limit the temperature, and therefore the luminosity. This natural limit to the maximum luminosity makes outbursts associated with bright persistent magnetars barely detectable. These events are simply seen as a small luminosity increase over the already bright quiescent state, followed by a fast return to quiescence. In particular, this is the case for 1RXS J1708-4009, 1E 1841-045, SGR 1806-20, and other bright persistent magnetars. On the other hand, a similar event (with the same energetics) in a fainter source will drive a more extreme luminosity variation and longer cooling time, as for sources such as XTE J1810-197, 1E 1547-5408 and SGR 1627-41. We conclude that the non-detection of large radiative outbursts in connection with glitches and bursts from bright persistent magnetars is not surprising per se, nor it needs of any revision on the glitches and burst mechanisms as explained by current theoretical models.

This volume contains most of the links to the presentations delivered at this international workshop. This meeting was the second in the series following the previous I Encuentro Ib\'erico de Compstar, held at the University of Coimbra, Portugal in 2010. The main purpose of this meeting was to strengthen the scientific collaboration between the participants of the Iberian and the rest of the southern European branches of the European Nuclear Astrophysics network, formerly, COMPSTAR. This ESF (European Science Foundation) supported network has been crucial in helping to make a broader audience for the the most interesting and relevant research lines being developed currently in Nuclear Astrophysics, especially related to the physics of neutron stars. The program of the meeting was tailored to theoretical descriptions of the physics of neutron stars although some input from experimental observers and other condensed matter and optics areas of interest was also included.

High magnetic fields are a distinguishing feature of neutron stars and the existence of sources (the soft gamma repeaters and the anomalous X-ray pulsars) hosting an ultra-magnetized neutron star (or magnetar) has been recognized in the past few decades. Magnetars are believed to be powered by magnetic energy and not by rotation, as with normal radio pulsars. Until recently, the radio quietness and magnetic fields typically above the quantum critical value (Bq~4.4x10^{13} G), were among the characterizing properties of magnetars. The recent discovery of radio pulsed emission from a few of them, and of a low dipolar magnetic field soft gamma repeater, weakened further the idea of a clean separation between normal pulsars and magnetars. In this Letter we show that radio emission from magnetars might be powered by rotational energy, similarly to what occurs in normal radio pulsars. The peculiar characteristics of magnetars radio emission should be traced in the complex magnetic geometry of these sources. Furthermore, we propose that magnetar radio activity or inactivity can be predicted from the knowledge of the star's rotational period, its time derivative and the quiescent X-ray luminosity.

We model neutrino emission from a newly born neutron star subsequent to a supernova explosion to study its sensitivity to the equation of state, neutrino opacities, and convective instabilities at high baryon density. We find the time period and spatial extent over which convection operates is sensitive to the behavior of the nuclear symmetry energy at and above nuclear density. When convection ends within the proto-neutron star, there is a break in the predicted neutrino emission that may be clearly observable.

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 necessary and sufficient condition for the thermodynamical universality of the static spherically symmetric Lovelock black hole is that it is the pure Lovelock {\Lambda}-vacuum solution. By universality we mean the thermodynamical parameters: temperature and entropy always bear the same relationship to the horizon radius irrespective of the Lovelock order and the spacetime dimension. For instance, the entropy always goes in terms of the horizon radius as rh and r^2 respectively for h odd and even dimensions. This universality uniquely identifies the pure Lovelock black hole with {\Lambda}.

Highly magnetized neutron stars (NSs) are characterized by a bewildering range of astrophysical manifestations. Here, building on our simulations of the evolution of magnetic stresses in the NS crust and its ensuing fractures (Perna & Pons 2011), we explore in detail, for the middle-age and old NSs, the dependence of starquake frequency and energetics on the relative strength of the poloidal (B_p) and toroidal (B_tor) components. We find that, for B_p >~10^{14}G, since a strong crustal toroidal field B_tor B_p is quickly formed on a Hall timescale, the initial toroidal field needs to be B_tor >> B_p to have a clear influence on the outbursting behaviour. For initial fields B_p <~ 10^{14}G, it is very unlikely that a middle-age (t~10^5 years) NS shows any bursting activity. This study allows us to solve the apparent puzzle of how NSs with similar dipolar magnetic fields can behave in a remarkably different way: an outbursting 'magnetar' with a high X-ray luminosity, or a quiet, low-luminosity, "high-$B$" radio pulsar. As an example, we consider the specific cases of the magnetar 1E2259+586 and the radio pulsar PSRJ1814-1744, which at present have a similar dipolar field ~6x10^{13}G. We determine for each object an initial magnetic field configuration that reproduces the observed timing parameters at their current age. The same two configurations also account for the differences in quiescent X-ray luminosity and for the 'magnetar/outbursting' behaviour of 1E2259+586 but not of PSRJ1814-1744. We further use the theoretically predicted surface temperature distribution to compute the light-curve for these objects. In the case of 1E2259+586, for which data are available, our predicted temperature distribution gives rise to a pulse profile whose double-peaked nature and modulation level is consistent with the observations.

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.

SGR 0418+5729 is a transient Soft Gamma-ray Repeater which underwent a major outburst in June 2009, during which the emission of short bursts was observed. Its properties appeared quite typical of other sources of the same class until long-term X-ray monitoring failed to detect any period derivative. The present upper limit on $\dot P$ implies that the surface dipole field is $B_p\lesssim 7.5\times 10^{12}\ {\rm G}$ (Rea et al 2010), well below those measured in other Soft Gamma-ray Repeaters (SGRs) and in the Anomalous X-ray Pulsars (AXPs), a group of similar sources. Both SGRs and AXPs are currently believed to be powered by ultra-magnetized neutron stars (magnetars, $B_p\approx 10^{14}$--$10^{15}\ {\rm G}$). SGR 0418+5729 hardly seems to fit in such a picture. We show that the magneto-rotational properties of SGR 0418+5729 can be reproduced if this is an aged magnetar, $\approx 1\ {\rm Myr}$ old, which experienced substantial field decay. The large initial toroidal component of the internal field required to match the observed properties of SGR 0418+5729 ensures that crustal fractures, and hence bursting activity, can still occur at present time. The thermal spectrum observed during the outburst decay is compatible with the predictions of a resonant compton scattering model (as in other SGRs/AXPs) if the field is low and the magnetospheric twist moderate.

We report on a 63ks Chandra observation of the X-ray transient Swift J195509.6+261406 discovered as the afterglow of what was first believed to be a long duration Gamma-Ray Burst (GRB 070610). The outburst of this source was characterized by unique optical flares on timescales of second or less, morphologically similar to the short X-ray bursts usually observed from magnetars. Our Chandra observation was performed ~2 years after the discovery of the optical and X-ray flaring activity of this source, catching it in its quiescent state. We derive stringent upper limits on the quiescent emission of Swif J195509.6+261406 which argues against the possibility of this object being a typical magnetar. Our limits show that the most viable interpretation on the nature of this peculiar bursting source, is a binary system hosting a black hole or a neutron star with a low mass companion star (< 0.12 M_{\odot}), and with an orbital period smaller than a few hours.

Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) are young neutron stars (NSs) characterized by high X-ray quiescent luminosities, outbursts, and, in the case of SGRs, sporadic giant flares. They are believed to be powered by ultra-strong magnetic fields (hence dubbed magnetars). The diversity of their observed behaviours is however not understood, and made even more puzzling by the discovery of magnetar-like bursts from "low-field" pulsars. Here we perform long-term 2D simulations that follow the evolution of magnetic stresses in the crust; these, together with recent calculations of the breaking stress of the neutron star crust, allow us to establish when starquakes occur. For the first time, we provide a quantitative estimate of the burst energetics, event rate, and location on the neutron star surface, which bear a direct relevance for the interpretation of the overall magnetar phenomenology. Typically, an "SGR-like" object tends to be more active than an "AXP-like" object or a "high-$B$ radio pulsar", but there is no fundamental separation among what constitutes the apparent different classes. Among the key elements that create the variety of observed phenomena, age is more important than a small variation in magnetic field strength. We find that outbursts can also be produced in old, lower-field pulsars (B ~ a few x 10^{12} G), but those events are much less frequent than in young, high-field magnetars.

We find general relativistic solutions of equilibrium magnetic field configurations in magnetars, extending previous results of Colaiuda et al. (2008). Our method is based on the solution of the relativistic Grad-Shafranov equation, to which Maxwell's equations can be reduced in some limit. We obtain equilibrium solutions with the toroidal magnetic field component confined into a finite region inside the star, and the poloidal component extending to the exterior. These so-called twisted-torus configurations have been found to be the final outcome of dynamical simulations in the framework of Newtonian gravity, and appear to be more stable than other configurations. The solutions include higher order multipoles, which are coupled to the dominant dipolar field. We use arguments of minimal energy to constrain the ratio of the toroidal to the poloidal field.

We perform population synthesis studies of different types of neutron stars taking into account the magnetic field decay. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G]) > \sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible. We apply the obtained field distribution and the model of decay to study long-term evolution of neuton stars till the stage of accretion from the interstellar medium. It is shown that though the subsonic propeller stage can be relatively long, initially highly magnetized neutron stars ($B_0 > \sim 10^{13}$ G) reach the accretion regime within the Galactic lifetime if their kick velocities are not too large. The fact that in previous studies made $>$10 years ago, such objects were not considered results in a slight increase of the Accretor fraction in comparison with earlier conclusions. Most of the neutron stars similar to the Magnificent seven are expected to become accreting from the interstellar medium after few billion years of their evolution. They are the main predecestors of accreting isolated neutron stars.

The non-linear Hall term present in the induction equation in the electron-magneto-hydrodynamics limit is responsible for the Hall drift of the magnetic field and, in some cases, for the so-called Hall instability. We investigate whether or not the growth rates and eigenfunctions found in the linear analysis are consistent with the results of non-linear numerical simulations. Following the linear analysis of Rheinhardt & Geppert, we study the same cases for which the Hall instability was predicted by solving the non-linear Hall induction equation using a two-dimensional conservative and divergence-free finite difference scheme that overcomes intrinsic difficulties of pseudo-spectral methods and can describe situations with arbitrarily high magnetic Reynolds numbers. We show that unstable modes can grow to the level of the background field without being overwhelmed by the Hall cascade, and cause a complete rearrangement of the field geometry. We confirm both the growth rates and eigenfunctions found in the linearized analysis and hence the instability. In the non-linear regime, after the unstable modes grow to the background level, the naturally selected modes become stable and oscillatory. Later on, the evolution tends to select the modes with the longest possible wavelengths, but this process occurs on the magnetic diffusion timescale. We confirm the existence of the Hall instability. We argue against using the misleading terminology that associates the non-linear Hall term with a turbulent Hall cascade, since small-scale structures are not created everywhere. The field evolves instead in a Burgers-like manner, forms local structures with strong gradients which become shocks in the zero resistivity limit, and Hall waves are launched and propagated through the entire domain.

We perform population synthesis studies of different types of neutron stars (thermally emitting isolated neutron stars, normal radio pulsars, magnetars) taking into account the magnetic field decay and using results from the most recent advances in neutron star cooling theory. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. For this purpose, we fix a baseline neutron star model (all microphysics input), and other relevant parameters to standard values (velocity distribution, mass spectrum, birth rates ...), allowing to vary the initial magnetic field strength. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G])>\sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life (only about a factor of 2 for low field neutron stars but more than an order of magnitude for magnetars), and a mass distribution function dominated by low mass objects. This model explains satisfactorily all known populations. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible.

We study the mutual influence of thermal and magnetic evolution in a neutron star's crust in axial symmetry. Taking into account realistic microphysical inputs, we find the heat released by Joule effect consistent with the circulation of currents in the crust, and we incorporate its effects in 2D cooling calculations. We solve the induction equation numerically using a hybrid method (spectral in angles, but a finite--differences scheme in the radial direction), coupled to the thermal diffusion equation. We present the first long term 2D simulations of the coupled magneto-thermal evolution of neutron stars. This substantially improves previous works in which a very crude approximation in at least one of the parts (thermal or magnetic diffusion) has been adopted. Our results show that the feedback between Joule heating and magnetic diffusion is strong, resulting in a faster dissipation of the stronger fields during the first million years of a NS's life. As a consequence, all neutron stars born with fields larger than a critical value (about 5 10^13 G) reach similar field strengths (approximately 2-3 10^{13} G) at late times. Irrespectively of the initial magnetic field strength, after $10^6$ years the temperature becomes so low that the magnetic diffusion timescale becomes longer than the typical ages of radio--pulsars, thus resulting in apparently no dissipation of the field in old NS. We also confirm the strong correlation between the magnetic field and the surface temperature of relatively young NSs discussed in preliminary works. The effective temperature of models with strong internal toroidal components are systematically higher than those of models with purely poloidal fields, due to the additional energy reservoir stored in the toroidal field that is gradually released as the field dissipates.

We report on a new mechanism for heat conduction in the neutron star crust. We find that collective modes of superfluid neutron matter, called superfluid phonons (sPhs), can influence heat conduction in magnetized neutron stars. They can dominate the heat conduction transverse to magnetic field when the magnetic field $B \gsim 10^{13}$ G. At density $\rho \simeq 10^{12}-10^{14} $ g/cm$^3$ the conductivity due to sPhs is significantly larger than that due to lattice phonons and is comparable to electron conductivity when temperature $\simeq 10^8$ K. This new mode of heat conduction can limit the surface anisotropy in highly magnetized neutron stars. Cooling curves of magnetized neutron stars with and without superfluid heat conduction could show observationally discernible differences.

Context: Many thermally emitting isolated neutron stars have magnetic fields larger than 10^13 G. A realistic cooling model that includes the presence of high magnetic fields should be reconsidered. Aims: We investigate the effects of anisotropic temperature distribution and Joule heating on the cooling of magnetized neutron stars. Methods: The 2D heat transfer equation with anisotropic thermal conductivity tensor and including all relevant neutrino emission processes is solved for realistic models of the neutron star interior and crust. Results: The presence of the magnetic field affects significantly the thermal surface distribution and the cooling history during both, the early neutrino cooling era and the late photon cooling era. Conclusions: There is a large effect of the Joule heating on the thermal evolution of strongly magnetized neutron stars. Both magnetic fields and Joule heating play a key role in keeping magnetars warm for a long time. Moreover, this effect is important for intermediate field neutron stars and should be considered in radio-quiet isolated neutron stars or high magnetic field radio-pulsars.

We present two-dimensional simulations for the cooling of neutron stars with strong magnetic fields (B > 1e13 Gauss). We study how the cooling curves are influenced by magnetic field decay. We show that the Joule heating effects are very large and in some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field and predict the surface temperature distribution for the early and late stages of the evolution of isolated neutron stars, comparing our results with available observational data of isolated neutron stars.

We find numerical solutions of the coupled system of Einstein-Maxwell's equations with a linear approach, in which the magnetic field acts as a perturbation of a spherical neutron star. In our study, magnetic fields having both poloidal and toroidal components are considered, and higher order multipoles are also included. We evaluate the deformations induced by different field configurations, paying special attention to those for which the star has a prolate shape. We also explore the dependence of the stellar deformation on the particular choice of the equation of state and on the mass of the star. Our results show that, for neutron stars with mass M = 1.4 Msun and surface magnetic fields of the order of 10^15 G, a quadrupole ellipticity of the order of 10^(-6) - 10^(-5) should be expected. Low mass neutron stars are in principle subject to larger deformations (quadrupole ellipticities up to 10^(-3) in the most extreme case). The effect of quadrupolar magnetic fields is comparable to that of dipolar components. A magnetic field permeating the whole star is normally needed to obtain negative quadrupole ellipticities, while fields confined to the crust typically produce positive quadrupole ellipticities.

The impact of strong magnetic fields B>10e13 G on the thermal evolution of neutron stars is investigated, including crustal heating by magnetic field decay. For this purpose, we perform 2D cooling simulations with anisotropic thermal conductivity considering all relevant neutrino emission processes for realistic neutron stars. The standard cooling models of neutron stars are called into question by showing that the magnetic field has relevant (and in many cases dominant) effects on the thermal evolution. The presence of the magnetic field significantly affects the thermal surface distribution and the cooling history of these objects during both, the early neutrino cooling era and the late photon cooling era. The minimal cooling scenario is thus more complex than generally assumed. A consistent magneto-thermal evolution of magnetized neutron stars is needed to explain the observations.

We present 2D simulations of the cooling of neutron stars with strong magnetic fields (B \geq 10^{13} G). We solve the diffusion equation in axial symmetry including the state of the art microphysics that controls the cooling such as slow/fast neutrino processes, superfluidity, as well as possible heating mechanisms. We study how the cooling curves depend on the the magnetic field strength and geometry. Special attention is given to discuss the influence of magnetic field decay. We show that Joule heating effects are very large and in some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field for the early and late stages of the evolution of isolated neutron stars.

In this article we study the possibility that, due to non-linear couplings, unstable g-modes associated to convective motions excite stable oscillating g-modes. This problem is of particular interest, since gravitational waves emitted by a newly born proto-neutron star pulsating in its stable g-modes would be in the bandwidth of VIRGO and LIGO. Our results indicate that nonlinear saturation of unstable modes occurs at relatively low amplitudes, and therefore, even if there exists a coupling between stable and unstable modes, it does not seem to be sufficiently effective to explain, alone, the excitation of the oscillating g-modes found in hydrodynamical simulations.

We study the non--linear evolution of magnetic fields in neutron star crusts with special attention to the influence of the Hall drift. Our goal is to understand the conditions for fast dissipation due to the Hall term in the induction equation. We study the interplay of Ohmic dissipation and Hall drift in order to find a timescale for the overall crustal field decay. We solve numerically the Hall induction equation by means of a hybrid method (spectral in angles but finite differences in the radial coordinate). The microphysical input consists of the most modern available crustal equation of state, composition and electrical conductivities. We present the first long term simulations of the non--linear magnetic field evolution in realistic neutron star crusts with a stratified electron number density and temperature dependent conductivity. We show that Hall drift influenced Ohmic dissipation takes place in neutron star crusts on a timescale of 1 Myr. When the initial magnetic field has magnetar strength, the fast Hall drift results in an initial rapid dissipation stage that lasts 10-50 kyr. The interplay of the Hall drift with the temporal variation and spatial gradient of conductivity tends to favor the displacement of toroidal fields toward the inner crust, where stable configurations can last for 1 Myr. We show that the thermally emitting isolated neutron stars, as the Magnificent Seven, are very likely descendants of neutron stars born as magnetars.

We derive the relativistic equations for stellar perturbations, including in a consistent way shear viscosity in the stress-energy tensor, and we numerically integrate our equations in the case of large viscosity. We consider the slow rotation approximation, and we neglect the coupling between polar and axial perturbations. In our approach, the frequency and damping time of the emitted gravitational radiation are directly obtained. We find that, approaching the inviscid limit from the finite viscosity case, the continuous spectrum is regularized. Constant density stars, polytropic stars, and stars with realistic equations of state are considered. In the case of constant density stars and polytropic stars, our results for the viscous damping times agree, within a factor two, with the usual estimates obtained by using the eigenfunctions of the inviscid limit. For realistic neutron stars, our numerical results give viscous damping times with the same dependence on mass and radius as previously estimated, but systematically larger of about 60%.

We show the existence of a strong trend between neutron star surface temperature and the dipolar component of the magnetic field extending through three orders of field magnitude, a range that includes magnetars, radio-quiet isolated neutron stars, and many ordinary radio pulsars. We suggest that this trend can be explained by the decay of currents in the crust over a time scale of few Myr. We estimate the minimum temperature that a NS with a given magnetic field can reach in this interpretation.

We present a unified explanation for the observed properties of the isolated neutron star RX J0720.4-3125 by obtaining, for the first time, a self-consistent model that accounts simultaneously for the observed X-ray spectrum and optical excess, the pulsed fraction, the observed spectral feature around 0.3 keV, and the long--term spectral evolution. By fitting the parameters of our realistic self--consistent models to all archival XMM--Newton observations and available optical data, we show that all observed properties are consistent with a normal neutron star with a proper radius of about 12 km, a temperature at the magnetic pole of about 100 eV and a magnetic field strength of 2-3 x 10^{13} G, with no need to invoke additional emission and absorption components nor exotic internal composition. The observed variability of the effective temperature, strength of the spectral feature, and pulsed fraction are in good agreement with the predictions of our model in which the star is subject to free precession, producing changes in the angle between the magnetic field and the rotation axis of tens of degrees with a periodicity of 7 years, as pointed out by other authors on the basis of phenomenological models. In addition to the evidence of internal toroidal components, we also find strong evidence of non-dipolar magnetic fields, since all spectral properties are better reproduced with models with strong quadrupolar components.

The thermal emission from isolated neutron stars is not well understood. The X-ray spectrum is very close to a blackbody but there is a systematic optical excess flux with respect to the extrapolation to low energy of the best blackbody fit. This fact, in combination with the observed pulsations in the X-ray flux, can be explained by anisotropies in the surface temperature distribution.We study the thermal emission from neutron stars with strong magnetic fields in order to explain the origin of the anisotropy. We find (numerically) stationary solutions in axial symmetry of the heat transportequations in the neutron star crust and the condensed envelope. The anisotropy in the conductivity tensor is included consistently. The presence of magnetic fields of the expected strength leads to anisotropy in the surface temperature. Models with toroidal components similar to or larger than the poloidal field reproduce qualitatively the observed spectral properties and variability of isolated neutron stars. Our models also predict spectral features at energies between 0.2 and 0.6 keV.

We have extended the procedure to find the exact solution of the Riemann problem in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics in which the magnetic field of the initial states is tangential to the discontinuity and orthogonal to the flow velocity. The wave pattern produced after the break up of the initial discontinuity is analogous to the non--magnetic case and we show that the problem can be understood as a purely relativistic hydrodynamical problem with a modified equation of state. The new degree of freedom introduced by the non-zero component of the magnetic field results in interesting effects consisting in the change of the wave patterns for given initial thermodynamical states, in a similar way to the effects arising from the introduction of tangential velocities. Secondly, when the magnetic field dominates the thermodynamical pressure and energy, the wave speeds approach the speed of light leading to fast shocks and fast and arbitrarily thin rarefaction waves. Our approach is the first non-trivial exact solution of a Riemann problem in relativistic magnetohydrodynamics and it can also be of great interest to test numerical codes against known analytical or exact solutions.

We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e. Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic magnetohydrodynamics equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic magnetohydrodynamics in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick accretion disks around a black hole prone to the magnetorotational instability.

Within a fully relativistic framework, we derive and solve numerically the perturbation equations of relativistic stars, including the stresses produced by a non-vanishing shear viscosity in the stress-energy tensor. With this approach, the real and imaginary parts of the frequency of the modes are consistently obtained. We find that, approaching the inviscid limit from the finite viscosity case, the continuous spectrum is regularized and we can calculate the quasi-normal modes for stellar models that do not admit solutions at first order in perturbation theory when the coupling between the polar and axial perturbations is neglected. The viscous damping time is found to agree within factor 2 with the usual estimate obtained by using the eigenfunctions of the inviscid limit and some approximation for the energy dissipation integrals. We find that the frequencies and viscous damping times for relativistic $r-$modes lie between the Newtonian and Cowling results. We compare the results obtained with homogeneous, polytropic and realistic equations of state and find that the frequencies depend only on the rotation rate and on the compactness parameter (M/R), being almost independent of the equation of state. Our numerical results for realistic neutron stars give viscous damping times with the same dependence on mass and radius as previously estimated, but systematically larger of about 60%.

We investigate the thermal emission from magnetic neutron star surfaces in which the cohesive effects of the magnetic field have produced the condensation of the atmosphere and the external layers. This may happen for sufficiently cool atmospheres with moderately intense magnetic fields. The thermal emission from an isothermal bare surface of a neutron star shows no remarkable spectral features, but it is significantly depressed at energies below some threshold energy. However, since the thermal conductivity is very different in the normal and parallel directions to the magnetic field lines, the presence of the magnetic field is expected to produce a highly anisotropic temperature distribution, depending on the magnetic field geometry. In this case, the observed flux of such an object looks very similar to a BB spectrum, but depressed in a nearly constant factor at all energies. This results in a systematic underestimation of the area of the emitter (and therefore its size) by a factor 5-10 (2-3).

We have developed a formalism to study non-adiabatic, non-radial oscillations of non-rotating compact stars in the frequency domain, including the effects of thermal diffusion in the framework of general relativistic perturbation theory. When a general equation of state depending on temperature is used, the perturbations of the fluid result in heat flux which is coupled with the spacetime geometry through the Einstein field equations. Our results show that the frequency of the first pressure (p) and gravity (g) oscillation modes is significantly affected by thermal diffusion, while that of the fundamental (f) mode is basically unaltered due to the global nature of that oscillation. The damping time of the oscillations is generally much smaller than in the adiabatic case (more than two orders of magnitude for the p- and g-modes) reflecting the effect of thermal dissipation. Both the isothermal and adiabatic limits are recovered in our treatment and we study in more detail the intermediate regime. Our formalism finds its natural astrophysical application in the study of the oscillation properties of newly born neutron stars, neutron stars with a deconfined quark core phase, or strange stars which are all promising sources of gravitational waves with frequencies in the band of the first generation and advanced ground-based interferometric detectors.

We study the effects of rotation on the quasi normal modes (QNMs) of a newly born proto neutron star (PNS) at different evolutionary stages, until it becomes a cold neutron star (NS). We use the Cowling approximation, neglecting spacetime perturbations, and consider different models of evolving PNS. The frequencies of the modes of a PNS are considerably lower than those of a cold NS, and are further lowered by rotation; consequently, if QNMs were excited in a sufficiently energetic process, they would radiate waves that could be more easily detectable by resonant-mass and interferometric detectors than those emitted by a cold NS. We find that for high rotation rates, some of the g-modes become unstable via the CFS instability; however, this instability is likely to be suppressed by competing mechanisms before emitting a significant amount of gravitational waves.

We study the conditions for convective instability in rotating, non-magnetic proto--neutron stars. The criteria that determine stability of nascent neutron stars are analogous to the Solberg--Hoiland conditions but including the presence of lepton gradients. Our results show that, for standard angular velocity profiles, convectively unstable modes with wave-vectors parallel to the rotation axis are suppressed by a stable angular momentum profile, while unstable modes with wave-vectors perpendicular to the axis remain unaltered. Since the wave-vector is perpendicular to the velocity perturbation, the directional selection of the unstable modes may result in fluid motions along the direction of the rotation axis. This occurs in rigidly rotating stars as well as in the inner core of differentially rotating stars. Our results provide a natural source of asymmetry for proto--neutron stars with the only requirement that angular velocities be of the order of the convective characteristic frequency.

We investigate the evolution of rigidly and differentially rotating protoneutron stars (PNSs) during the first twenty seconds of their life. We solve the equations describing stationary axisymmetric configurations in general relativity coupled to a finite temperature, relativistic equation of state, to obtain a sequence of quasi-equilibrium configurations describing the evolution of newly born neutron stars. Our estimates show that the scale of variation of the angular velocity in a PNSs is of the order of 7-10 km. We obtain the maximum rotation frequency that can be reached as the protoneutron stars deleptonizes and cools down, as well as other relevant parameters such as total angular momentum or |T/W|. Our study shows that imposing physical constraints (conservation of baryonic mass and angular momentum) and choosing reasonable thermodynamical profiles as the star evolves gives results consistent with the energetics of more complex simulations of non-rotating PNSs. It appears to be unlikely that PNSs formed in nearly axisymmetric core collapse reach the critical angular velocity to undergo the bar mode instability. They could, however, undergo secular or low |T/W| rotational instabilities a few seconds after birth, resulting in a strong emission of gravitational waves retarded with respect to the neutrino luminosity peak. We also found that the geometry of strongly differentially rotating protoneutron stars can become toroidal-like for large values of the angular velocity, before reaching the mass shedding limit.

In this paper we study the effects of rotation on the frequencies of the quasi-normal modes of a proto-neutron star (PNS) born in a gravitational collapse during the first minute of life. Our analysis starts a few tenths of seconds after the PNS formation, when the stellar evolution can be described by a sequence of equilibrium configurations. We use the evolutionary models developed by Pons et al. (1999; 2001) that describe how a non rotating star cools down and contracts while neutrino diffusion and thermalization processes dominate the stellar dynamics. For assigned values of the evolution time, we set the star into slow rotation and integrate the equations of stellar perturbations in the Cowling approximation, both in the time domain and in the frequency domain, to find the quasi-normal mode frequencies. We study the secular instability of the g-modes, that are present in the oscillation spectrum due to the intense entropy and composition gradients that develop in the stellar interior, and we provide an estimate of the growth time of the unstable modes based on a post-Newtonian formula.

We study the gravitational radiation associated to the non--radial oscillations of newly born, hot neutron stars. The frequencies and damping times of the relevant quasi--normal modes are computed for two different models of proto--neutron stars, at different times of evolution, from its birth until it settles down as a cold neutron star. We find that the oscillation properties of proto--neutron stars are remarkably different from those of their cold, old descendants and that this affects the characteristic features of the gravitational signal emitted during the post-collapse evolution. The consequences on the observability of these signals by resonant--mass and interferometric detectors are analyzed. We find that gravitational waves from the pulsations of a newborn proto--neutron star in the galaxy could be detected with a signal to noise ratio of 5 by the first generation interferometers, if the energy stored in the modes is greater than $\sim 10^{-8} M_\odot c^2$, or by a resonant antenna if it is greater than $\sim 10^{-4} M_\odot c^2$. In addition since at early times the frequency of the spacetime modes is much lower than that of a cold neutron star, they would be also detectable with the same signal to noise ratio if a comparable amount of energy is radiated into these modes.

We compute the overlap function between Post-Newtonian (PN) templates and gravitational signals emitted by binary systems composed of one neutron star and one point mass, obtained by a perturbative approach. The calculations are performed for different stellar models and for different detectors, to estimate how effectual and faithful the PN templates are, and to establish whether effects related to the internal structure of neutron stars may possibly be extracted by the matched filtering technique.

A phase transition occurring in the inner core of a neutron star could be associated to a density discontinuity that would affect the frequency spectrum of the non-radial oscillation modes in two ways. Firstly, it would produce a softening of the equation of state, leading to more compact equilibrium configurations and changing the frequency of the fundamental and pressure modes of the neutron star. Secondly, a new non-zero frequency g-- mode would appear, associated to each discontinuity. These discontinuity g--modes have typical frequencies larger than those of g--modes previously studied in the literature (thermal, core g-- modes, or g--modes due to chemical inhomogeneities in the outer layers), and smaller than that of the fundamental mode; therefore they should be distinguishable from the other modes of non radial oscillation. In this paper we investigate how high density discontinuities change the frequency spectrum of the non-radial oscillations, in the framework of the general relativistic theory of stellar perturbations. Our purpose is to understand whether a gravitational signal, emitted at the frequencies of the quasi normal modes, may give some clear information on the equation of state of the neutron star and, in particular, on the parameters that characterize the density discontinuity. We discuss some astrophysical processes that may be associated to the excitation of these modes, and estimate how much gravitational energy should the modes convey to produce a signal detectable by high frequency gravitational detectors.

We extend our approach for the exact solution of the Riemann problem in relativistic hydrodynamics to the case in which the fluid velocity has components tangential to the initial discontinuity. As in one-dimensional flows, we here show that the wave-pattern produced in a multidimensional relativistic Riemann problem can be predicted entirely by examining the initial conditions. Our method is logically very simple and allows for a numerical implementation of an exact Riemann solver which is both straightforward and computationally efficient. The simplicity of the approach is also important for revealing special relativistic effects responsible for a smooth transition from one wave-pattern to another when the tangential velocities in the initial states are suitably varied. While the content of this paper is focussed on a flat spacetime, the local Lorentz invariance allows its use also in fully general relativistic calculations.

The stability properties of newly born neutron stars, or proto--neutron stars, are considered. We take into account dissipative processes, such as neutrino transport and viscosity, in the presence of a magnetic field. In order to find the regions of the star subject to different sorts of instability, we derive the general instability criteria and apply it to evolutionary models of PNSs. The influence of the magnetic field on instabilities is analyzed and the critical magnetic field stabilizing the star is obtained. In the light of our results, we estimate of the maximum poloidal magnetic field that might be present in young pulsars or magnetars.

The effects that the structure of a neutron star would have on the gravitational emission of a binary system are studied in a perturbative regime, and in the frequency domain. Assuming that a neutron star is perturbed by a point mass moving on a close, circular orbit, we solve the equations of stellar perturbations in general relativity to evaluate the energy lost by the system in gravitational waves. We compare the energy output obtained for different stellar models with that found by assuming that the perturbed object is a black hole with the same mass, and we discuss the role played by the excitation of the stellar modes. Ouresults indicate that the stellar structure begins to affect the emitted power when the orbital velocity is v >0.2c (about 185 Hz for a binary system composed of two canonical neutron stars). We show that the differences between different stellar models and a black hole are due mainly to the excitation of the quasinormal modes of the star. Finally, we discuss to what extent and up to which distance the perturbative approach can be used to describe the interaction of a star and a pointlike massive body.

We discuss efforts to determine the mass, radius, and surface composition of the nearby compact object RX J185635-3754 from its multi-wavelength spectral energy distribution. We compute non-magnetized model atmospheres and emergent spectra for selected compositions and gravities, and discuss efforts to fit existing and new observational data from ROSAT, EUVE and the HST. The spectral energy distribution matches that expected from a heavy-element dominated atmosphere, but not from a uniform temperature blackbody. Non-magnetic light element atmospheres cannot be simultaneously reconciled with the optical and X-ray data. We extend previous studies, which were limited to one fixed neutron star mass and radius. For uniform temperature models dominated by heavy elements, the redshift z is constrained to be 0.3 < z < 0.4 and the best-fit mass and radius are M approx 0.9 solar masses and R approx 6 km (for a 61 pc distance). These values for M and R together are not permitted for any plausible equation of state, including that of a self-bound strange quark star. A simplified two-temperature model allows masses and radii up to about 50% larger, or a factor of 2 in the case of a black body. The observed luminosity is consistent with the thermal emission of an isolated neutron star no older than about 1 million years, the age inferred from available proper motion and parallax information

We compute the energy spectra of the gravitational signals emitted when a pointlike mass moves on a closed orbit around a non rotating neutron star, inducing a perturbation of its gravitational field and its internal structure. The Einstein equations and the hydrodynamical equations are perturbed and numerically integrated in the frequency domain. The results are compared with the energy spectra computed by the quadrupole formalism which assumes that both masses are pointlike, and accounts only for the radiation emitted because the orbital motion produces a time dependent quadrupole moment. The results of our perturbative approach show that, in general, the quadrupole formalism overestimates the amount of emitted radiation, especially when the two masses are close. However, if the pointlike mass is allowed to move on an orbit so tight that the keplerian orbital frequency resonates with the frequency of the fundamental quasi-normal mode of the star (2w_K=w_f), this mode can be excited and the emitted radiation can be considerably larger than that computed by the quadrupole approach.

Neutrino fluxes from proto-neutron stars with and without quarks are studied. Observable differences become apparent after 10--20 s of evolution. Sufficiently massive stars containing negatively-charged, strongly interacting, particles collapse to black holes during the first minute of evolution. Since the neutrino flux vanishes when a black hole forms, this is the most obvious signal that quarks (or other types of strange matter) have appeared. The metastability timescales for stars with quarks are intermediate between those containing hyperons and kaon condensates.

The existence of compact dimensions which are accessible only to gravity represents an intriguing possible solution to the hierarchy problem. At present the strongest constraint on the existence of such compact Gravity-Only Dimensions (GODs) comes from SN1987a. Here we report on the first self-consistent simulations of the early, neutrino-emitting phase of a proto-neutron star which include energy losses due to the coupling of the Kaluza-Klein modes of the graviton which arise in a world with GODs. We compare the neutrino signals from these simulations to that from SN1987a and use a rigorous probabilistic analysis to derive improved bounds for the radii of such GODs. We find that the possibility that there are two compact extra dimensions with radii larger than 0.66 $\mu$m is excluded at the 95% confidence level---as is the possibility that there are three compact extra dimensions larger than 0.8 nm.

The main stages in the evolution of a neutron star, from its birth as a proto-neutron star, to its old age as a cold, catalyzed configuration, are described. A proto-neutron star is formed in the aftermath of a successful supernova explosion and its evolution is dominated by neutrino diffusion. Its neutrino signal is a valuable diagnostic of its internal structure and composition. During its transformation from a hot, lepton-rich to a cold, catalyzed remnant, the possibility exists that it can collapse into a black hole, which abruptly terminates neutrino emissions. The essential microphysics, reviewed herein, that controls its evolution are the equation of state of dense matter and its associated neutrino opacities. Several simulations of the proto-neutron star evolution, involving different assumptions about the composition of dense matter, are described. After its evolution into a nearly isothermal neutron star a hundred or so years after its birth, it may be observable through its thermal emission in X-rays during its life in the next million years. Its surface temperature will depend upon the rapidity of neutrino emission processes in its core, which depends on the composition of dense matter and whether or not its constituents exhibit superfluidity and superconductivity. Observations of thermal emission offer the best hope of a determination of the radius of a neutron star. The implications for the underlying dense matter equation of state of an accurate radius determination are explored.

We present simulations of the evolution of a proto-neutron star in which kaon-condensed matter might exist, including the effects of finite temperature and trapped neutrinos. The phase transition from pure nucleonic matter to the kaon condensate phase is described using Gibbs' rules for phase equilibrium, which permit the existence of a mixed phase. A general property of neutron stars containing kaon condensates, as well as other forms of strangeness, is that the maximum mass for cold, neutrino-free matter can be less than the maximum mass for matter containing trapped neutrinos or which has a finite entropy. A proto-neutron star formed with a baryon mass exceeding that of the maximum mass of cold, neutrino-free matter is therefore metastable, that is, it will collapse to a black hole at some time during the Kelvin-Helmholtz cooling stage. The effects of kaon condensation on metastable stars are dramatic. In these cases, the neutrino signal from a hypothetical galactic supernova (distance $\sim8.5$ kpc) will stop suddenly, generally at a level above the background in the SuperK and SNO detectors, which have low energy thresholds and backgrounds. This is in contrast to the case of a stable star, for which the signal exponentially decays, eventually disappearing into the background. We find the lifetimes of kaon-condensed metastable stars to be restricted to the range 40--70 s and weakly dependent on the proto-neutron star mass, in sharp contrast to the significantly larger mass dependence and range (1--100 s) of hyperon-rich metastable stars.

We study the equation of state (EOS) of kaon-condensed matter including the effects of temperature and trapped neutrinos. It is found that the order of the phase transition to a kaon-condensed phase, and whether or not Gibbs' rules for phase equilibrium can be satisfied in the case of a first order transition, depend sensitively on the choice of the kaon-nucleon interaction. The main effect of finite temperature, for any value of the lepton fraction, is to mute the effects of a first order transition, so that the thermodynamics becomes similar to that of a second order transition. Above a critical temperature, found to be at least 30--60 MeV depending upon the interaction, the first order transition disappears. The phase boundaries in baryon density versus lepton number and baryon density versus temperature planes are delineated. We find that the thermal effects on the maximum gravitational mass of neutron stars are as important as the effects of trapped neutrinos, in contrast to previously studied cases in which the matter contained only nucleons or in which hyperons and/or quark matter were considered. Kaon-condensed EOSs permit the existence of metastable neutron stars, because the maximum mass of an initially hot, lepton-rich protoneutron star is greater than that of a cold, deleptonized neutron star. The large thermal effects imply that a metastable protoneutron star's collapse to a black hole could occur much later than in previously studied cases that allow metastable configurations.

The linear hydrodynamic stability of proto-neutron stars (PNSs) is considered taking into account dissipative processes such as neutrino transport and viscosity. We obtain the general instability criteria which differ essentially from the well-known Ledoux criterion used in previous studies. We apply the criteria to evolutive models of PNSs that, in general, can be subject to the various known regimes such as neutron fingers and convective instabilities. Our results indicate that the fingers instability arises in a more extended region of the stellar volume and lasts a longer time than expected.

We study the mathematical character of the angular moment equations of radiative transfer in spherical symmetry and conclude that the system is hyperbolic for general forms of the closure relation found in the literature. Hyperbolicity and causality preservation lead to mathematical conditions allowing to establish a useful characterization of the closure relations. We apply numerical methods specifically designed to solve hyperbolic systems of conservation laws (the so-called Godunov-type methods), to calculate numerical solutions of the radiation transport equations in a static background. The feasibility of the method in any kind of regime, from diffusion to free-streaming, is demonstrated by a number of numerical tests and the effect of the choice of the closure relation on the results is discussed.

We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on the tangential velocities is analysed, and the impact of this result on the development of multidimensional relativistic hydrodynamic codes (of Godunov type) is discussed.

Our contribution concerns with the numerical solution of the 3D general relativistic hydrodynamical system of equations within the framework of the 3+1 formalism. We summarize the theoretical ingredients which are necessary in order to build up a numerical scheme based on the solution of local Riemann problems. Hence, the full spectral decomposition of the Jacobian matrices of the system, i.e., the eigenvalues and the right and left eigenvectors, is explicitly shown. An alternative approach consists in using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits an accurate description of numerical general relativistic hydrodynamics.

We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on the tangential velocities is analysed. This solution has been used to build up an exact Riemann solver implemented in a multidimensional relativistic (Godunov-type) hydro-code.

We derive and analyze a simplified formulation of the numerical viscosity terms appearing in the expression of the numerical fluxes associated to several High-Resolution Shock-Capturing schemes. After some algebraic pre-processing, we give explicit expressions for the numerical viscosity terms of two of the most widely used flux formulae, which implementation saves computational time in multidimensional simulations of relativistic flows. Additionally, such treatment explicitely cancells and factorizes a number of terms helping to amortiguate the growing of round-off errors. We have checked the performance of our formulation running a 3D relativistic hydrodynamical code to solve a standard test-bed problem and found that the improvement in efficiency is of high practical interest in numerical simulations of relativistic flows in Astrophysics.

An extensive study of the effects of correlations on both charged and neutral current weak interaction rates in dense matter is performed. Both strong and electromagnetic correlations are considered.The propagation of particle-hole interactions in the medium plays an important role in determining the neutrino mean free paths. The effects due to Pauli-Blocking and density, spin, and isospin correlations in the medium significantly reduce the neutrino cross sections. Due to the lack of experimental information at high density, these correlations are necessarily model dependent. For example, spin correlations in nonrelativistic models are found to lead to larger suppressions of neutrino cross sections compared to those of relativistic models. This is due to the tendency of the nonrelativistic models to develop spin instabilities. Notwithstanding the above caveats, and the differences between nonrelativistic and relativistic approaches such as the spin- and isospin-dependent interactions and the nucleon effective masses, suppressions of order 2--3, relative to the case in which correlations are ignored, are obtained. Neutrino interactions in dense matter are especially important for supernova and early neutron star evolution calculations. The effects of correlations for protoneutron star evolution are calculated. Large effects on the internal thermodynamic properties of protoneutron stars, such as the temperature, are found. These translate into significant early enhancements in the emitted neutrino energies and fluxes, especially after a few seconds. At late times, beyond about 10 seconds, the emitted neutrino fluxes decrease more rapidly compared to simulations without the effects of correlations, due to the more rapid onset of neutrino transparency in the protoneutron star

We present a general and practical procedure to solve the general relativistic hydrodynamic equations by using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits accurate numerical calculations of general relativistic hydrodynamics problems using the numerical tools developed for the special relativistic case with negligible computational cost. The feasibility of the method has been confirmed by a number of numerical experiments.

We study the thermal and chemical evolution during the Kelvin-Helmholtz phase of the birth of a neutron star, employing neutrino opacities that are consistently calculated with the underlying equation of state (EOS). Expressions for the diffusion coefficients appropriate for general relativistic neutrino transport in the equilibrium diffusion approximation are derived. The diffusion coefficients are evaluated using a field-theoretical finite temperature EOS that includes the possible presence of hyperons. The variation of the diffusion coefficients is studied as a function of EOS and compositional parameters. We present results from numerical simulations of protoneutron star cooling for internal stellar properties as well as emitted neutrino energies and luminosities. We discuss the influence of the initial stellar model, the total mass, the underlying EOS, and the addition of hyperons on the evolution of the protoneutron star and upon the expected signal in terrestrial detectors.

We calculate the Legendre expansion of the rate of the process $\nu + \bar{\nu} \leftrightarrow e^+ + e^-$ up to 3rd order extending previous results of other authors which only consider the 0th and 1st order terms. Using different closure relations for the moment equations of the radiative transfer equation we discuss the physical implications of taking into account quadratic and cubic terms on the energy deposition outside the neutrinosphere in a simplified model. The main conclusion is that 2nd order is necessary in the semi-transparent region and gives good results if an appropriate closure relation is used.

The early evolution of a protoneutron star depends both on the equation of state and neutrino interactions at high density. We identify the important sources of neutrino opacity and through model calculations show that in-medium effects play an important role in determining the neutrino mean free paths. The effects due to Pauli-blocking and many-body correlations due to strong interactions reduce the neutrino cross sections by large factors compared to the case in which these effects are ignored. We discuss these results in the context of neutrino transport in a protoneutron star.

We study the diffusive transport of neutrinos in a newly born neutron star to explore its sensitivity to dense matter properties. Energy and lepton number which are trapped during the catastrophic implosion diffuse out on the time scale of a few tens of seconds. Results for different dense matter models are presented.