A young cooling neutron star in the remnant of Supernova 1987A

We describe the cooling theory for isolated neutron stars that are several tens of years old. Their cooling differs greatly from the cooling of older stars that has been well studied in the literature. It is sensitive to the physics of the inner stellar crust and even to the thermal conductivity of the stellar core, which is never important at later cooling stages. The absence of observational evidence for the formation of a neutron star during the explosion of Supernova 1987A is consistent with the fact that the star was actually born there. It may still be hidden in the dense center of the supernova remnant. If, however, the star is not hidden, then it should have a low thermal luminosity (below ～10³⁴ erg s^?1) and a short internal thermal relaxation time (shorter than 13 yr). This requires that the star undergo intense neutrino cooling (e.g., via the direct Urca process) and have a thin crust with strong superfluidity of free neutrons and/or an anomalously high thermal conductivity.     

Simulations of glitches in isolated pulsars

Many radio pulsars exhibit glitches wherein the star's spin rate increases fractionally by ∼10⁻¹⁰–10⁻⁶. Glitches are ascribed to variable coupling between the neutron star crust and its superfluid interior. With the aim of distinguishing among different theoretical explanations for the glitch phenomenon, we study the response of a neutron star to two types of perturbations to the vortex array that exists in the superfluid interior: (1) thermal motion of vortices pinned to inner crust nuclei, initiated by sudden heating of the crust, (e.g., a starquake), and (2) mechanical motion of vortices (e.g., from crust cracking by superfluid stresses). Both mechanisms produce acceptable fits to glitch observations in four pulsars, with the exception of the 1989 glitch in the Crab pulsar, which is best fitted by the thermal excitation model. The two models make different predictions for the generation of internal heat and subsequent enhancement of surface emission. The mechanical glitch model predicts a negligible temperature increase. For a pure and highly conductive crust, the thermal glitch model predicts a surface temperature increase of as much as ∼2 per cent, occurring several weeks after the glitch. If the thermal conductivity of the crust is lowered by a high concentration of impurities, however, the surface temperature increases by ∼10 per cent about a decade after a thermal glitch. A thermal glitch in an impure crust is consistent with the surface emission limits following the 2000 January glitch in the Vela pulsar. Future surface emission measurements coordinated with radio observations will constrain glitch mechanisms and the conductivity of the crust.     

Minimal models of cooling neutron stars with accreted envelopes

We study the 'minimal' cooling scenario of superfluid neutron stars with nucleon cores, where the direct Urca process is forbidden and enhanced cooling is produced by neutrino emission due to the Cooper pairing of neutrons. Extending our recent previous work, we include the effects of surface accreted envelopes of light elements. We employ the phenomenological density-dependent critical temperatures   T _cp(ρ)  and   T _cnt(ρ)  of singlet-state proton and triplet-state neutron pairing in a stellar core, as well as the critical temperature   T _cns(ρ)  of singlet-state neutron pairing in a stellar crust. We show that the presence of accreted envelopes simplifies the interpretation of observations of thermal radiation from isolated neutron stars in the scenario of our recent previous work and widens the class of models for nucleon superfluidity in neutron star interiors consistent with the observations.     

Heating and cooling of magnetars with accreted envelopes

We study the thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in an internal layer. We focus on the effect of magnetized (   B ≳ 10¹⁴  G) non-accreted and accreted outermost envelopes composed of different elements, from iron to hydrogen or helium. We discuss a combined effect of thermal conduction and neutrino emission in the outer neutron star crust and calculate the cooling of magnetars with a dipole magnetic field for various locations of the heat layer, heat rates and magnetic field strengths. Combined effects of strong magnetic fields and light-element composition simplify the interpretation of magnetars in our model: these effects allow one to interpret observations assuming less extreme (therefore, more realistic) heating. Massive magnetars, with fast neutrino cooling in their cores, can have higher thermal surface luminosity.     

A numerical study of the r-mode instability of rapidly rotating nascent neutron stars

The first results of numerical analysis of classical r-modes of rapidly rotating compressible stellar models are reported. The full set of linear perturbation equations of rotating stars in Newtonian gravity is solved numerically without the slow rotation approximation. A critical curve of gravitational wave emission induced instability, which restricts the rotational frequencies of hot young neutron stars, is obtained. Taking the standard cooling mechanisms of neutron stars into account, we also show the 'evolutionary curves' along which neutron stars are supposed to evolve as cooling and spinning down proceed. Rotational frequencies of 1.4-M_⊙ stars suffering from this instability decrease to around 100 Hz when the standard cooling mechanism of neutron stars is employed. This result confirms the results of other authors, who adopted the slow rotation approximation.     

Afterburst thermal relaxation in neutron star crusts

We study thermal relaxation in a neutron star after internal heating events (outbursts) in the crust. We consider thin and thick spherically symmetric heaters, superfluid and non-superfluid crusts, stars with open and forbidden direct Urca processes in their cores. In particular, we analyze long-term thermal relaxation after deep crustal heating produced by nuclear transformations in fully or partly accreted crusts of transiently accreting neutron stars. This long-term relaxation has a typical relaxation time and an overall finite duration time for the crust to thermally equilibrate with the core. Neutron superfluidity in the inner crust greatly affects the relaxation if the heater is located in the inner crust. It shortens and unifies the time of emergence of thermal wave from the heater to the surface. This is important for the interpretation of observed outbursts of magnetars and transiently accreting neutron stars in quasi-persistent low-mass X-ray binaries.     

Self-similar pressure oscillations in neutron star envelopes as probes of neutron star structure

We study eigenmodes of acoustic oscillations of high multipolarity l ∼ 100–1000 and high frequency (∼100 kHz), localized in neutron star envelopes. We show that the oscillation problem is self-similar. Once the oscillation spectrum is calculated for a given equation of state (EOS) in the envelope and given stellar mass M and radius R , it can be rescaled to a star with any M and R (but the same EOS in the envelope). For l ≳ 300, the modes can be subdivided into the outer and inner ones. The outer modes are mainly localized in the outer envelope. The inner modes are mostly localized near the neutron drip point, being associated with the softening of the EOS after the neutron drip. We calculate oscillation spectra for the EOSs of cold-catalyzed and accreted matter and show that the spectra of the inner modes are essentially different. A detection and identification of high-frequency pressure modes would allow one to infer M and R and determine also the EOS in the envelope (accreted or ground state) providing a new, potentially powerful method to explore the main parameters and internal structure of neutron stars.     

Atmospheres and spectra of strongly magnetized neutron stars

We construct atmosphere models for strongly magnetized neutron stars with surface fields     and effective temperatures     . The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars, including radio pulsars, soft gamma-ray repeaters, and anomalous X-ray pulsars. In our models, the atmosphere is composed of pure hydrogen or helium and is assumed to be fully ionized. The radiative opacities include free–free absorption and scattering by both electrons and ions computed for the two photon polarization modes in the magnetized electron–ion plasma. Since the radiation emerges from deep layers in the atmosphere with     , plasma effects can significantly modify the photon opacities by changing the properties of the polarization modes. In the case where the magnetic field and the surface normal are parallel, we solve the full, angle-dependent, coupled radiative transfer equations for both polarization modes. We also construct atmosphere models for general field orientations based on the diffusion approximation of the transport equations and compare the results with models based on full radiative transport. In general, the emergent thermal radiation exhibits significant deviation from blackbody, with harder spectra at high energies. The spectra also show a broad feature     around the ion cyclotron resonance     , where Z and A are the atomic charge and atomic mass of the ion, respectively; this feature is particularly pronounced when     . Detection of the resonance feature would provide a direct measurement of the surface magnetic fields on magnetars.     

Time-variable emission from transiently accreting neutron stars in quiescence due to deep crustal heating

Transiently accreting neutron stars in quiescence ( L _X ≲10³⁴ erg s⁻¹) have been observed to vary in intensity by factors of few, over time-scales of days to years. If the quiescent luminosity is powered by a hot neutron star core, the core cooling time-scale is much longer than the recurrence time, and cannot explain the observed, more rapid variability. However, the non-equilibrium reactions which occur in the crust during outbursts deposit energy in isodensity shells, from which the thermal diffusion time-scale to the photosphere is days to years. The predicted magnitude of variability is too low to explain the observed variability unless – as is widely believed – the neutrons beyond the neutron-drip density are superfluid. Even then, the variability due to this mechanism in models with standard core neutrino cooling processes is less than 50 per cent – still too low to explain the reported variability. However, models with rapid core neutrino cooling can produce a variability by a factor as great as 20, on time-scales of days to years following an outburst. Thus, the factors of ∼ few intensity variability observed from transiently accreting neutron stars can be accounted for by this mechanism only if rapid core cooling processes are active.     

Crustal impurities and the internal temperature of a neutron star crust

When the upper part of a neutron star crystallizes to form the crust, the constituting ions are trapped in the lattice as a result of the low diffusion rates in the solid phase. As a consequence, the local composition of the crust corresponds to the equilibrium condition at the melting point and not at the present internal temperature. The inclusion of the small entropic contribution to the free energy does not lead to marked changes in our view of the microscopic structure of a neutron star crust, because the melting temperature is much smaller than the typical energies at play in the crystal cell. However, mixing between layers characterized by different nuclear species is found to be an important effect for the production of impurities. We show that one should expect an increase of the thermal diffusion time through the crust at small temperatures, because of the decrease of thermal conductivity in relatively thin impurity-rich layers acting as isolating shields.     

Resistive relaxation of a magnetically confined mountain on an accreting neutron star

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to  τ_I∼ 10⁵–10⁸ yr  (depending on the conductivity of the neutron star crust) for an accreted mass of   M _a= 1.2 × 10⁻⁴ M_⊙  . The magnetic dipole moment simultaneously re-emerges as the screening currents dissipate over  τ_I  . For non-axisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfvén time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease  τ_I  by one order of magnitude relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.     

Unveiling the thermal and magnetic map of neutron star surfaces though their X-ray emission: method and light-curve analysis

Recent Chandra and XMM–Newton observations of a number of X-ray 'dim' pulsating neutron stars have revealed quite unexpected features in the emission from these sources. Their soft thermal spectrum, believed to originate directly from the star surface, shows evidence for a phase-varying absorption line at some hundred eVs. The pulse modulation is relatively large (pulsed fractions in the range ∼12–35 per cent), the pulse shape is often non-sinusoidal, and the hard X-ray colour appears to be anticorrelated in phase with the total emission. Moreover, the prototype of this class, RX J0720.4−3125, has been found to undergo rather sensible changes in both its spectral and timing properties over a time-scale of a few years. All these new findings seem difficult to reconcile with the standard picture of a cooling neutron star endowed with a purely dipolar magnetic field, at least if surface emission is produced in an atmosphere on top of the crust. In this paper we explore how a dipolar+quadrupolar star-centred field influences the properties of the observed light curves. The phase-resolved spectrum has been evaluated accounting for both radiative transfer in a magnetized atmosphere and general relativistic ray-bending. We computed over 78 000 light curves, varying the quadrupolar components and the viewing geometry. A comparison of the data with our model indicates that higher-order multipoles are required to reproduce the observations.     

Almost analytic models of ultramagnetized neutron star envelopes

Recent ROSAT measurements show that the X-ray emission from isolated neutron stars is modulated at the stellar rotation period. To interpret these measurements, one needs precise calculations of the heat transfer through the thin insulating envelopes of neutron stars. We present nearly analytic models of the thermal structure of the envelopes of ultramagnetized neutron stars. Specifically, we examine the limit in which only the ground Landau level is filled. We use the models to estimate the amplitude of modulation expected from non-uniformities in the surface temperatures of strongly magnetized neutron stars. In addition, we estimate cooling rates for stars with fields B  ∼ 10¹⁵ − 10¹⁶ G, which are relevant to models that invoke 'magnetars' to account for soft γ-ray emission from some repeating sources.     

Magnetized relativistic jets and long-duration GRBs from magnetar spin-down during core-collapse supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spin-down of a newly formed millisecond,   B ∼ 10¹⁵ G  , magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spin-down powers  (∼10⁵¹–10⁵² erg s⁻¹)  , the magnetar wind is superfast at almost all latitudes, while for lower spin-down powers  (∼10⁵⁰ erg s⁻¹)  , the wind is subfast but still super-Alfvénic. In all cases, the rates at which the neutron star loses mass, angular momentum and energy are very similar to the corresponding free wind values (≲30 per cent differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated  (∼5–10°)  relativistic jet out along the rotation axis of the star. Nearly all of the spin-down power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.     

Crust–core coupling and r-mode damping in neutron stars: a toy model

r-modes in neutron stars with crusts are damped by viscous friction at the crust–core boundary. The magnitude of this damping, evaluated by Bildsten & Ushomirsky (BU) under the assumption of a perfectly rigid crust, sets the maximum spin frequency for neutron stars spun up by accretion in low-mass X-ray binaries (LMXBs). In this paper we explore the mechanical coupling between the core r-modes and the elastic crust, using a toy model of a constant-density neutron star having a crust with a constant shear modulus. We find that, at spin frequencies in excess of ≈50 Hz, the r-modes strongly penetrate the crust. This reduces the relative motion (slippage) between the crust and the core compared with the rigid-crust limit. We therefore revise down, by as much as a factor of 10²–10³ , the damping rate computed by BU, significantly reducing the maximal possible spin frequency of neutron stars with solid crusts. The dependence of the crust–core slippage on the spin frequency is complicated, and is very sensitive to the physical thickness of the crust. If the crust is sufficiently thick, the curve of the critical spin frequency for the onset of the r-mode instability becomes multivalued for some temperatures; this is related to avoided crossings between the r-mode and higher-order torsional modes in the crust. The critical frequencies are comparable to the observed spins of neutron stars in LMXBs and millisecond pulsars.     

What is the nature of RX J0720.4–3125?

RX J0720.4–3125 has recently been identified as a pulsating soft X-ray source in the ROSAT all-sky survey with a period of 8.391 s. Its spectrum is well characterized by a blackbody with a temperature of 8 × 10⁵ K. We propose that the radiation from this object is thermal emission from a cooling neutron star. For this blackbody temperature we can obtain a robust estimate of the object's age of ∼ 3 × 10⁵ yr, yielding a polar field ∼ 10¹⁴ G for magnetic dipole spin-down and a value of P compatible with current observations.     

On the dynamics of superfluid neutron star cores

We discuss the nature of the various modes of pulsation of superfluid neutron stars using comparatively simple Newtonian models and the Cowling approximation. The matter in these stars is described in terms of a two-fluid model, where one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars, and the other fluid represents a conglomerate of all other constituents (crust nuclei, protons, electrons, etc.). In our model, we incorporate the non-dissipative interaction known as the entrainment effect, whereby the momentum of one constituent (e.g. the neutrons) carries along part of the mass of the other constituent. We show that there is no independent set of pulsating g-modes in a non-rotating superfluid neutron star core, even though the linearized superfluid equations contain a well-defined (and real-valued) analogue to the so-called Brunt–Väisälä frequency. Instead, what we find are two sets of spheroidal perturbations whose nature is predominately acoustic. In addition, an analysis of the zero-frequency subspace (i.e. the space of time-independent perturbations) reveals two sets of degenerate spheroidal perturbations, which we interpret to be the missing g-modes, and two sets of toroidal perturbations. We anticipate that the degeneracy of all these zero-frequency modes will be broken by the Coriolis force in the case of rotating stars. To illustrate this we consider the toroidal pulsation modes of a slowly rotating superfluid star. This analysis shows that the superfluid equations support a new class of r-modes, in addition to those familiar from, for example, geophysical fluid dynamics. Finally, the role of the entrainment effect on the superfluid mode frequencies is shown explicitly via solutions to dispersion relations that follow from a 'local' analysis of the linearized superfluid equations.     

A Chandra X-ray observation of the globular cluster Terzan 1

We present a ∼19-ks Chandra Advanced CCD Imaging Spectrometer (ACIS)-S observation of the globular cluster Terzan 1. 14 sources are detected within 1.4 arcmin of the cluster centre with two of these sources predicted to be not associated with the cluster (background active galactic nuclei or foreground objects). The neutron star X-ray transient, X1732−304, has previously been observed in outburst within this globular cluster with the outburst seen to last for at least 12 yr. Here, we find four sources that are consistent with the ROSAT position for this transient, but none of the sources are fully consistent with the position of a radio source detected with the Very Large Array that is likely associated with the transient. The most likely candidate for the quiescent counterpart of the transient has a relatively soft spectrum and an unabsorbed 0.5–10 keV luminosity of  2.6 × 10³² erg s⁻¹  , quite typical of other quiescent neutron stars. Assuming standard core cooling, from the quiescent flux of this source we predict long (>400 yr) quiescent episodes to allow the neutron star to cool. Alternatively, enhanced core cooling processes are needed to cool down the core. However, if we do not detect the quiescent counterpart of the transient this gives an unabsorbed 0.5–10 keV luminosity upper limit of  8 × 10³¹ erg s⁻¹  . We also discuss other X-ray sources within the globular cluster. From the estimated stellar encounter rate of this cluster we find that the number of sources we detect is significantly higher than expected by the relationship of Pooley et al.     

Resonant cyclotron scattering and Comptonization in neutron star magnetospheres

Resonant cyclotron scattering of the surface radiation in the magnetospheres of neutron stars may considerably modify the emergent spectra and impede efforts to constrain neutron star properties. Resonant cyclotron scattering by a non-relativistic warm plasma in an inhomogeneous magnetic field has a number of unusual characteristics. (i) In the limit of high resonant optical depth, the cyclotron resonant layer is half opaque , in sharp contrast to the case of non-resonant scattering. (ii) The transmitted flux is on average Compton up-scattered by  ∼1 + 2β_T  , where  β_T  is the typical thermal velocity in units of the velocity of light; the reflected flux has on average the initial frequency. (iii) For both the transmitted and reflected fluxes, the dispersion of intensity decreases with increasing optical depth. (iv) The emergent spectrum is appreciably non-Planckian while narrow spectral features produced at the surface may be erased.
We derive semi-analytically modification of the surface Planckian emission due to multiple scattering between the resonant layers and apply the model to the anomalous X-ray pulsar 1E  1048.1 − 5937  . Our simple model fits just as well as the 'canonical' magnetar spectra model of a blackbody plus power law.     

Relativistic Correction of Neutrino Emission in Neutron Stars

By the relativistic mean field theory and relevant weak-interactional cooling theory, the relativistic cooling properties in the conventional and hyperonic neutron star matter are studied. Also a comparison between the relativistic and non-relativistic results after taking consideration of the gravity correction is performed. The results show that the relativistic effect of neutrino emission reduces the neutrino emissivity, luminosity, and the cooling rate of stellar objects, in comparison with the non-relativistic case. In the neutron star matter without hyperon, the amplitude of the cooling rate reduction caused by the relativistic effect is maximal after taking the gravity correction into consideration, it attains 56% for a 2 M_⊙ neutron star composed of conventional neutron star matter, and in the hyperonic matter the amplitude of reduction is minimal, about 38%.