Are pulsars born with a hidden magnetic field?
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 …
Crust-magnetosphere coupling during magnetar evolution and implications for the surface temperature
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 …
The exact solution of the Riemann problem with non-zero tangential velocities in 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, and the impact of this result on the development of multidimensional relativistic hydrodynamic codes (of Godunov type) is discussed.
Gravitational Waves from Rotating Proto-Neutron Stars
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 bec…
"Gravitational waves from newly born, hot neutron stars"
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 ar…
Gravitational waves from neutron stars at different evolutionary stages
We study how the internal structure of a neutron star and the physical processes that may occur during its evolution affect the quasi-normal mode spectrum, and consequently the gravitational radiation it emits. We discuss whether these modes can be excited and how much energy they should carry for the gravitational signal to be detectable by the first generation of interferometric antennas or by the new generation of high-frequency gravitational detectors, interferometric or resonant, that are under investigation.
Effects of Strong and Electromagnetic Correlations on Neutrino Interactions in Dense Matter
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 …
Kaon condensation in proto-neutron star matter
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 t…
Properties of the isolated neutron star RX J185635-3754
Abstract Despite improved observational capabilities, the surface properties of neutron stars remain uncertain. In principle, multiwavelength spectra can reveal their compositions and angular diameters and constrain the interior equation of state. The Chandra LETG spectrum of the brightest isolated neutron star is a featureless continuum. The X-ray blackbody fit underpredicts the observed optical fluxes, but the spectral energy distribution can be fit with a two blackbody model. No pulsations are seen. We discuss model atmosphere fits. We conclude there is no need to invoke exotic physics: RX J185635-3754 is likely a 0.5 × 10 6 years old pulsar.
The force-free twisted magnetosphere of a neutron star – II. Degeneracies of the Grad–Shafranov equation
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 dipo…
Evolution of Proto-Neutron stars with kaon condensates
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 …
The Large Observatory For x-ray Timing
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 WideFi…
Rotational effects on the oscillation frequencies of newly born proto-neutron stars
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 th…
Instability of twisted magnetar magnetospheres
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 f…
The force-free twisted magnetosphere of a neutron star
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 ma…
Long-term evolution of the force-free twisted magnetosphere of a magnetar
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 sho…
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 Jac…
An Exact Riemann Solver for Multidimensional Special Relativistic Hydrodynamics
We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics (Marti and Muller, 1994) 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. This solution has been used to build up an exact Riemann solver implemented in a multidimensional relativistic (Godunov-type) hydro-code.
Hyperbolic character of the angular moment equations of radiative transfer and numerical methods
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 demonstrat…
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.
Evolutionary sequences of rotating protoneutron stars
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 mom…
Neutron star formation with presence of hyperons
We study the influence of hyperons during the early stages of the birth of a neutron star (Kelvin-Helmholtz phase), employing neutrino opacities calculated consistently with the equation of state by considering all possible neutrino-hyperon reactions. Our results from numerical simulations of newly born neutron stars, or proto-neutron stars, show an increasingly important influence of hyperons at later times. It is remarkable the existence of metastable stars, which are stable at birth but become unstable during the evolution as the deleptonization proceeds and the hyperon concentration increases. We also present results from hydrodynamical simulations of the collapse to a black hole of met…
Accretion-driven gravitational radiation from nonrotating compact objects. Infalling quadrupolar shells
This paper reports results from numerical simulations of the gravitational radiation emitted from non--rotating compact objects(both neutron stars and Schwarzschild black holes) as a result of the accretion of matter. A hybrid procedure is adopted: we evolve, in axisymmetry, the linearized equations describing metric and fluid perturbations, coupled with a nonlinear hydrodynamics code that calculates the motion of the accreting matter. The initial matter distribution is shaped in the form of extended quadrupolar shells of dust or perfect fluid. Self--gravity and radiation reaction effects of the accreting fluid are neglected. This idealized setup is used to understand the qualitative featur…