Juan A. Miralles

Evolution of Proto-Neutron Stars

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 proto-neutron 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 proto-neutron star and on the expected signal in terrestrial detectors. We find that the differences in predicted luminosities and emitted neutrino energies do not depend much upon the details of the initial models or the underlying high-density EOS for early times (t<10 s), provided that opacities are calculated consistently with the EOS. The same holds true for models that allow for the presence of hyperons, except when the initial mass is significantly larger than the maximum mass for cold, catalyzed matter. For times larger than about 10 s, and prior to the occurrence of neutrino transparency, the neutrino luminosities decay exponentially with a time constant that is sensitive to the high-density properties of matter. We also find the average emitted neutrino energy increases during the first 5 s of evolution and then decreases nearly linearly with time. In general, increasing the proto-neutron star mass increases the average energy and the luminosity of neutrinos, as well as the overall evolutionary timescale. The influence of hyperons or variations in the dense matter EOS is increasingly important at later times. Metastable stars, those with hyperons that are unstable to collapse upon deleptonization, have relatively long evolution times, which increase the nearer the mass is to the maximum mass supportable by a cold, deleptonized star.

Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radiopulsars, 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.

Numerical 3+1 General Relativistic Magnetohydrodynamics: A Local Characteristic Approach

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 MHD 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 MHD 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 disks accreting onto a black hole and subject to the magnetorotational instability.

2D Cooling of magnetized neutron stars

Context. Many thermally emitting, isolated neutron stars have magnetic fields that are larger than 10 13 G. A realistic cooling model that includes the presence of high magnetic fields should be r econsidered. Aims. We investigate the effects of an 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 emissi on processes is solved for realistic models of the neutron star interior and crust. Results. The presence of the magnetic field a ffects 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 Joule heating on the thermal evolution of strongly magnetized neutron stars. Both magnetic fields and Joule heating play an important role in keeping mag netars 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 ra dio‐pulsars.

The Impact of Magnetic Field on the Thermal Evolution of Neutron Stars

The impact of strong magnetic fields B > 1013 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 magnetothermal evolution of magnetized neutron stars is needed to explain the observations.

Evidence for heating of neutron stars by magnetic-field decay.

We show the existence of a strong trend between neutron star (NS) 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 approximately 10(6) yr. We estimate the minimum temperature that a NS with a given magnetic field can reach in this interpretation.

A new spherically symmetric general relativistic hydrodynamical code

In this paper we present a full general relativistic one-dimensional hydro-code which incorporates a modern high-resolution shock-capturing algorithm, with an approximate Riemann solver, for the correct modelling of formation and propagation of strong shocks. The efficiency of this code in treating strong shocks is demonstrated by some numerical experiments. The interest of this technique in several astrophysical scenarios is discussed.

A new code for the Hall-driven magnetic evolution of neutron stars

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

RELATIVISTIC MAGNETOHYDRODYNAMICS: RENORMALIZED EIGENVECTORS AND FULL WAVE DECOMPOSITION RIEMANN SOLVER

We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wave front in the fluid rest frame. Proper expressions of the renormalized eigenvectors in conserved variables are obtained through the corresponding matrix transformations. Our work completes previous analysis that present different sets of right eigenvectors for non-degenerate and degenerate states, and can be seen as a relativistic generalization of earlier work performed in classical MHD. Based on the full wave decomposition (FWD) provided by the renormalized set of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) Riemann solver. Extensive testing against one- and two-dimensional standard numerical problems allows us to conclude that our solver is very robust. When compared with a family of simpler solvers that avoid the knowledge of the full characteristic structure of the equations in the computation of the numerical fluxes, our solver turns out to be less diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The amount of operations needed by the FWD solver makes it less efficient computationally than those of the HLL family in one-dimensional problems. However, its relative efficiency increases in multidimensional simulations.

Convective instability in proto-neutron stars

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.