Pulsar timing irregularities and the imprint of magnetic field evolution

Context. 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. Aims. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable consequences in the braking index n. This is the most sensitive quantity to reflect small variations of the timing properties that are caused by magnetic field rearrangements. Methods. We performed axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art microphysical inputs to calculate the evolution of the braking index. Relatively rapid magnetic field modifications can be expected only in the crust of neutron stars, where we focus our study. Results. 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 a high absolute value and both positive and negative signs. Conclusions. Current magneto-thermal evolution models predict a large contribution to the timing noise and, in particular, to the braking index, from temporal variations of the magnetic field. Models with strong (≳ 1014 G) multipolar or toroidal components, even with a weak (~1012 G) dipolar field are consistent with the observed trend of the timing properties.

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.

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

The accretion environment in Vela X-1 during a flaring period using XMM-Newton

We present analysis of 100 ks contiguous XMM-Newton data of the prototypical wind accretor Vela X-1. The observation covered eclipse egress between orbital phases 0.134 and 0.265, during which a giant flare took place, enabling us to study the spectral properties both outside and during the flare. This giant flare with a peak luminosity of 3.92+0.42-0.09 × 1037 erg s-1 allows estimates of the physical parameters of the accreted structure with a mass of ~1021 g. We have been able to model several contributions to the observed spectrum with a phenomenological model formed by three absorbed power laws plus three emission lines. After analysing the variations with orbital phase of the column density of each component, as well as those in the Fe and Ni fluorescence lines, we provide a physical interpretation for each spectral component. Meanwhile, the first two components are two aspects of the principal accretion component from the surface of the neutron star, the third component seems to be the X-ray light echo formed in the stellar wind of the companion.

Spectral features in isolated neutron stars induced by inhomogeneous surface temperatures

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.

The many lives of magnetized neutron stars

The magnetic field strength at birth is arguably one of the most important properties to determine the evolutionary path of a neutron star. Objects with very high fields, collectively known as magnetars, are characterized by high X-ray quiescent luminosities, occurrence of outbursts, and, for some of them, sporadic giant flares. While the magnetic field strength is believed to drive their collective behaviour, however, the diversity of their properties, and, especially, the observation of magnetar-like bursts from “low-field” pulsars, has been a theoretical puzzle. In this review, we discuss results of long-term simulations following the coupled evolution of the X-ray luminosity and the timing properties for a large, homogeneous sample of X-ray emitting isolated neutron stars, accounting for a range of initial magnetic field strengths, envelope compositions, and neutron star masses. In addition, by following the evolution of magnetic stresses within the neutron star crust, we can also relate the observed magnetar phenomenology to the physical properties of neutron stars, and in particular to their age and magnetic field strength and topology. The dichotomy of “high-B” field pulsars versus magnetars is naturally explained, and occasional outbursts from old, low B-field neutron stars are predicted. We conclude by speculating on the fate of old magnetars, and by presenting observational diagnostics of the neutron star crustal field topology.

Probing large-scale wind structures in Vela X–1 using off-states with INTEGRAL

Vela X–1 is the prototype of the class of wind-fed accreting pulsars in high-mass X-ray binaries hosting a supergiant donor. We have analysed in a systematic way 10 years of INTEGRAL data of Vela X–1 (22–50 keV) and we found that when outside the X-ray eclipse, the source undergoes several luminosity drops where the hard X-rays luminosity goes below ∼3 × 1035 erg s−1, becoming undetected by INTEGRAL. These drops in the X-ray flux are usually referred to as ‘off-states’ in the literature. We have investigated the distribution of these off-states along the Vela X–1 ∼ 8.9 d orbit, finding that their orbital occurrence displays an asymmetric distribution, with a higher probability to observe an off-state near the pre-eclipse than during the post-eclipse. This asymmetry can be explained by scattering of hard X-rays in a region of ionized wind, able to reduce the source hard X-ray brightness preferentially near eclipse ingress. We associate this ionized large-scale wind structure with the photoionization wake produced by the interaction of the supergiant wind with the X-ray emission from the neutron star. We emphasize that this observational result could be obtained thanks to the accumulation of a decade of INTEGRAL data, with observations covering the whole orbit several times, allowing us to detect an asymmetric pattern in the orbital distribution of off-states in Vela X–1.

Quiescent thermal emission from neutron stars in low-mass X-ray binaries

Context. 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). Aims. 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. Methods. 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, XTE J1701−462, EXO 0748−676  and IGR J17480−2446. Results. 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.

On the formation and evolution of the first Be star in a black hole binary MWC 656

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict ∼10–100 of B-BH binaries to currently reside in the Galactic disc, among which around 1/3 contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be-BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding black hole (BH) binaries, with BH X-ray novae formation issue being the prime example. Future evolution of MWC 656 with an ∼5 M⊙ BH and with an ∼13 M⊙ main-sequence companion on an ∼60 d orbit may lead to the formation of a coalescing BH–NS (neutron star) system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to ∼0.2 yr−1. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH–NS binaries.

On the Radial Onset of Clumping in the Wind of the B0I Massive Star QV Nor

We present an analysis of a 78 ks Chandra high-energy transmission gratings observation of the B0I star QV Nor, the massive donor of the wind-accreting pulsar 4U1538−52. The neutron star (NS) orbits its companion in a very close orbit (r < 1.4R*, in units of the stellar radii), thereby allowing probing of the innermost wind regions. The flux of the Fe Kα line during eclipse reduces to only ∼30% of the flux measured out of eclipse. This indicates that the majority of Fe fluorescence must be produced in regions close to the NS, at distances smaller than 1R* from its surface. The fact that the flux of the continuum decreases to only ∼3% during eclipse allows for a high contrast of the Fe Kα line fluorescence during eclipse. The line is not resolved and centered at 1.9368 0.0018 l = 0.0032 - + Å. From the inferred plasma speed limit of v < c l < 800 l D km s−1 and range of ionization parameters of log 1, 2 x = [- ], together with the stellar density profile, we constrain the location of the cold, dense material in the stellar wind of QV Nor using simple geometrical considerations. We then use the Fe Kα line fluorescence as a tracer of wind clumps and determine that these clumps in the stellar wind of QV Nor (B0I) must already be present at radii r < 1.25R*, close to the photosphere of the star.

Neutron Stars—Cooling and Transport

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.

The VLT-FLAMES Tarantula survey: XX. The nature of the X-ray bright emission-line star VFTS 399

Context: The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. 

Aims: One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity-in this study we aim to determine its physical nature and the cause of this behaviour. 

Methods: To accomplish this we performed a time-resolved analysis of optical, infrared and X-ray observations. 

Results: We found VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select He i lines-taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. 

Conclusions: VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor, sharing many observational characteristics with classical Be X-ray binaries. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass 25 M_⊙ for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating "single" O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.

A simple approach to the supernova progenitor–explosion connection

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parametrized approaches and find good agreement with semi-empirical measures for the ‘explodability’ of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

Progenitor-Dependent Explosion Dynamics in Self-Consistent, Axisymmetric Simulations of Neutrino-Driven Core-Collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M ⊙, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

The role of nuclear physics in supernovae and the evolution of neutron stars, Neutrino Opacities, Equation of State, Transport Coefficients, and Dark Matter Production

Thesis (Ph.D.)--University of Washington, 2016-08