Search

A generic characteristic of self-graviting systems is that they have a negative heat capacity. An important example of this behaviour is given by the Schwarzschild black hole. The case of charged and rotating black holes is even more interesting, since a change of sign of the specific heat takes place through an infinite discontinuity. This has been usually associated with a black hole thermodynamic phase transition appearing at the points where the heat capacity diverges, the so-called Davies points. This aspect of black hole thermodynamics has been addressed from different perspectives, motivating different interpretations since its discovery in the 1970s. In this paper, a physical reinterpretation of the heat capacity is provided for spherically symmetric and static black holes. Our analysis is partially based on a reformulation of the black hole heat capacity using the Newman--Penrose formalism. The application to the Reissner--Nordstr\"om--de Sitter black hole case reveals a clear physical interpretation of the Newman--Penrose scalars evaluated at the event horizon. This allows us to write the heat capacity as a balance of pressures defined at the horizon. In particular, a matter pressure (coming from the energy-momentum tensor) and a thermal pressure (coming from the holographic energy equipartition of the horizon). The Davies point is identified with the point where the Komar thermal energy density matches the matter pressure at the horizon. We also compare the black hole case with the case of self-graviting objects and their corresponding thermal evolutions. We conclude that the heat capacity of black holes and self-graviting systems can be understood qualitatively in similar terms., Comment: 9 pages, 2 figures

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., Comment: 9 pages, 9 figures; updated version accepted for publication by MNRAS

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., Comment: 12 pages, 9 figures; added missing reference; accepted for publication in MNRAS

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., Comment: 10 pages, 8 figures; corrected typos; accepted for publication in MNRAS

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., Comment: 6 pages, version accepted by MNRAS

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., Comment: 14 pages, 11 figures, version accepted for publication in MNRAS

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., Comment: 17 pages, 9 figures; corrected typos; accepted for publication in MNRAS

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., Comment: 8 pages (emulated), 7 figures; ApJ in press

We analyze the influence of the magnetic field in the convexity properties of the relativistic magnetohydrodynamics system of equations. To this purpose we use the approach of Lax, based on the analysis of the linearly degenerate/genuinely non-linear nature of the characteristic fields. Degenerate and non-degenerate states are discussed separately and the non-relativistic, unmagnetized limits are properly recovered. The characteristic fields corresponding to the material and Alfv\'en waves are linearly degenerate and, then, not affected by the convexity issue. The analysis of the characteristic fields associated with the magnetosonic waves reveals, however, a dependence of the convexity condition on the magnetic field. The result is expressed in the form of a generalized fundamental derivative written as the sum of two terms. The first one is the generalized fundamental derivative in the case of purely hydrodynamical (relativistic) flow. The second one contains the effects of the magnetic field. The analysis of this term shows that it is always positive leading to the remarkable result that the presence of a magnetic field in the fluid reduces the domain of thermodynamical states for which the EOS is non-convex., Comment: 14 pages. Submitted to Classical and Quantum Gravity

Isolated neutron stars show a diversity in timing and spectral properties, which has historically led to a classification in different sub-classes. The magnetic field plays a key role in many aspects of the neutron star phenomenology: it regulates the braking torque responsible for their timing properties and, for magnetars, it provides the energy budget for the outburst activity and high quiescent luminosities (usually well above the rotational energy budget). We aim at unifying this observational variety by linking the results of the state-of-the-art 2D magneto-thermal simulations with observational data. The comparison between theory and observations allows to place two strong constraints on the physical properties of the inner crust. First, strong electrical currents must circulate in the crust, rather than in the star core. Second, the innermost part of the crust must be highly resistive, which is in principle in agreement with the presence of a novel phase of matter so-called nuclear pasta phase., Comment: 12 pages, 3 figures. Proceedings of the invited talk at SEA2014, as best Spanish PhD thesis in Astronomy and Astrophysics 2013

We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in parameter space that are more consistent with the observational data. In particular, we find that any broad distribution of birth spin periods with $P_0\lesssim 0.5$ s can fit the data, and that if the alignment angle is allowed to vary consistently with the torque model, realistic magnetospheric models are favoured compared to models with classical magneto-dipolar radiation losses. Assuming that the initial magnetic field is given by a lognormal distribution, our optimal model has mean strength $\langle\log B_0{\rm [G]}\rangle \approx 13.0-13.2$ with width $\sigma (\log B_0) = 0.6-0.7$. However, there are strong correlations between parameters. This degeneracy in the parameter space can be broken by an independent estimate of the pulsar birth rate or by future studies correlating this information with the population in other observational bands (X-rays and $\gamma$-rays)., Comment: 10 pages, 9 figures, submitted and accepted to MNRAS, comments welcome

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

The relativistic hydrodynamic system of equations for a perfect fluid obeying a causal equation of state is hyperbolic (Anile 1989 {\it Relativistic Fluids and Magneto-Fluids} (Cambridge: Cambridge University Press)). In this report, we derive the conditions for this system to be convex in terms of the fundamental derivative of the equation of state (Menikoff and Plohr 1989 {\it Rev. Mod. Phys.} {\bf 61} 75). The classical limit is recovered., Comment: 10 pages

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., Comment: 14 pages, 13 figures, Comput. Phys. Comm., in press

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., Comment: 16 pages, 17 figures, A&A accepted

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., Comment: 7 pages, 3 figures, oral contribution to "XXIII Trobades cientifiques de la Mediterrania - Supernovae: light in the darkness", Mao (Menorca-Spain) October 3 - 5, 2007

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., Comment: 4 pages, 3 figures, 1 table. Accepted for publication in APJ Letters

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., Comment: 4 pages, 3 figures, oral contribution to "EXOCT 2007, International Symposium on Exotic States of Nuclear Matter", Catania, Italy, 11-15 June 2007

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., Comment: 17 pages, 16 figures, version accepted for publication in A&A

Qualitatively new aspects of the (linear and non-linear) stability of sheared relativistic (slab) jets are analyzed. The linear problem has been solved for a wide range of jet models well inside the ultrarelativistic domain (flow Lorentz factors up to 20; specific internal energies $\approx 60c^2$). As a distinct feature of our work, we have combined the analytical linear approach with high-resolution relativistic hydrodynamical simulations, which has allowed us i) to identify, in the linear regime, resonant modes specific to the relativistic shear layer ii) to confirm the result of the linear analysis with numerical simulations and, iii) more interestingly, to follow the instability development through the non-linear regime. We find that very high-order reflection modes with dominant growth rates can modify the global, long-term stability of the relativistic flow. We discuss the dependence of these resonant modes on the jet flow Lorentz factor and specific internal energy, and on the shear layer thickness. The results could have potential applications in the field of extragalactic relativistic jets., Comment: Accepted for publication in Physical Review E. For better quality images, please check http://www.mpifr-bonn.mpg.de/staff/mperucho/Research.html

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., Comment: 17 pages, 5 figures, accepted for publication in the Journal of Fluid Mechanics

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., Comment: 18 pages, 4 figures, to appear in Astrophysical Journal

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., Comment: 25 pages, 14 figures. Submitted to Astrophysical Journal

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., Comment: 17 pages, 3 figures. Accepted for publication in Astrophysical Journal

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., Comment: 37 pags, 12 figures, accepted for publication in MNRAS

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., Comment: 12 pages, 2 figures, accepted for publication in Astronomy and Astrophysics

We consider Joule heating caused by dissipation of the magnetic field in the neutron star crust. This mechanism may be efficient in maintaining a relatively high surface temperature in very old neutron stars. Calculations of the thermal evolution show that, at the late evolutionary stage ($t \geq 10$ Myr), the luminosity of the neutron star is approximately equal to the energy released due to the field dissipation and is practically independent of the atmosphere models. At this stage, the surface temperature can be of the order of $3 \times 10^{4} - 10^{5}$K. Joule heating can maintain this high temperature during extremely long time ($\geq 100$ Myr), comparable with the decay time of the magnetic field., Comment: 13 pages (5 figures in the text). Accepted for publication in The Astrophysical Journal

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., Comment: 27 pages, 26 figures, accepted ApJ

The equation of state of an ideal Fermi gas is expressed in terms of Fermi-Dirac integrals. We give formulae for evaluation the Fermi-Dirac integrals of orders 1/2, 3/2, and 5/2 and their derivatives in various limits of non- and extreme degeneracy and relativity. We provide tables and a Fortran subroutine for numerical evaluation of the integrals and derivatives when a limit does not apply. The functions can be evaluated to better than 1% accuracy for any temperature and density using these methods., Comment: AAS LaTeX. No Figures