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 imprint of the crustal magnetic field on the thermal spectra and pulse profiles of isolated neutron stars

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the-art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double-peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50–60 per cent level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated ‘observed’ spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ∼1–2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements.

A VLT/FLAMES survey for massive binaries in Westerlund 1. III. The WC9d binary W239 and implications for massive stellar evolution

Context. There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars – frequency, mass ratio & orbital separation – are still poorly constrained. Aims. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-IR photometric observations in order to search for binary signatures. Comparison of these data to synthetic non-LTE model atmosphere spectra were used to derive the fundamental properties of the WC9 primary. Results. We found W239 to have an orbital period of only ~5.05 days, making it one of the most compact WC binaries yet identified. Analysis of the long term near-IR lightcurve reveals a significant flare between 2004-6. We interpret this as evidence for a third massive stellar component in the system in a long period (>6 yr), eccentric orbit, with dust production occuring at periastron leading to the flare. The presence of a near-IR excess characteristic of hot (~1300 K) dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short (<1 yr) period binaries, although confirmation will require further observations. Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent with other Galactic examples. Conclusions. The simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss – likely via case A mass transfer or a contact configuration – or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants.

Central compact objects and the hidden magnetic field scenario

Central compact objects (CCOs) are X-ray sources lying close to the centre of supernova remnants, with inferred values of the surface magnetic fields significantly lower (≲1011 G) than those of standard pulsars. In this paper, we revise the hidden magnetic field scenario, presenting the first 2D simulations of the submergence and re-emergence of the magnetic field in the crust of a neutron star. A post-supernova accretion stage of about 10−4–10−3 M⊙ over a vast region of the surface is required to bury the magnetic field into the inner crust. When accretion stops, the field re-emerges on a typical time-scale of 1–100 kyr, depending on the submergence conditions. After this stage, the surface magnetic field is restored close to its birth values. A possible observable consequence of the hidden magnetic field is the anisotropy of the surface temperature distribution, in agreement with observations of several of these sources. We conclude that the hidden magnetic field model is viable as an alternative to the antimagnetar scenario, and it could provide the missing link between CCOs and the other classes of isolated neutron stars.

IRAS 18357-0604 – an analogue of the galactic yellow hypergiant IRC +10420?

Context. Yellow hypergiants represent a short-lived evolutionary episode experienced by massive stars as they transit to and from a red supergiant phase. As such, their properties provide a critical test of stellar evolutionary theory, while recent observations unexpectedly suggest that a subset may explode as Type II supernovae. Aims. The galactic yellow hypergiant IRC +10420 is a cornerstone system for understanding this phase since it is the strongest post-RSG candidate known, has demonstrated real-time evolution across the Hertzsprung-Russell diagram and been subject to extensive mass loss. In this paper we report on the discovery of a twin of IRC +10420 - IRAS 18357-0604. Methods. Optical and near-IR spectroscopy are used to investigate the physical properties of IRAS 18357-0604 and also provide an estimate of its systemic velocity, while near- to mid-IR photometry probes the nature of its circumstellar environment. Results. These observations reveal pronounced spectral similarities between IRAS 18357-0604 and IRC +10420, suggesting comparable temperatures and wind geometries. IR photometric data reveals a similarly dusty circumstellar environment, although historical mass loss appears to have been heavier in IRC +10420. The systemic velocity implies a distance compatible with the red supergiant-dominated complex at the base of the Scutum Crux arm; the resultant luminosity determination is consistent with a physical association but suggests a lower initial mass than inferred for IRC +10420 (≲20 M⊙ versus ~40 M⊙). Evolutionary predictions for the physical properties of supernova progenitors derived from ~18–20 M⊙ stars – or ~12–15 M⊙ stars that have experienced enhanced mass loss as red supergiants – compare favourably with those of IRAS 18357-0604, which in turn appears to be similar to the the progenitor of SN2011dh; it may therefore provide an important insight into the nature of the apparently H-depleted yellow hypergiant progenitors of some Type IIb SNe.

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.

Population synthesis of isolated neutron stars with magneto-rotational evolution

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 P0 ≲ 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 〈log B0[G]〉 ≈ 13.0–13.2 with width σ(log B0) = 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 γ-rays).

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.

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.

The Gaia-ESO Survey: The analysis of high-resolution UVES spectra of FGK-type stars

Context. The ongoing Gaia-ESO Public Spectroscopic Survey is using FLAMES at the VLT to obtain high-quality medium-resolution Giraffe spectra for about 105 stars and high-resolution UVES spectra for about 5000 stars. With UVES, the Survey has already observed 1447 FGK-type stars. Aims. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO second internal release and will be part of its first public release of advanced data products. Methods. The final parameter scale is tied to the scale defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. In addition, a set of open and globular clusters is used to evaluate the physical soundness of the results. Each of the implemented methodologies is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted medians of those from the individual methods. Results. The recommended results successfully reproduce the atmospheric parameters of the benchmark stars and the expected Teff-log  g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log  g and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50−100 K for Teff, 0.10−0.25 dex for log  g and 0.05−0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex. Conclusions. The Gaia-ESO sample of high-resolution spectra of FGK-type stars will be among the largest of its kind analyzed in a homogeneous way. The extensive list of elemental abundances derived in these stars will enable significant advances in the areas of stellar evolution and Milky Way formation and evolution.

The Gaia-ESO Survey: radial metallicity gradients and age-metallicity relation of stars in the Milky Way disk

We study the relationship between age, metallicity, and α-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpcto 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages >9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more α-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

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.

Population synthesis of isolated neutron stars with magneto-rotational evolution – II. From radio-pulsars to magnetars

Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include X-ray thermal emitting pulsars within the same evolutionary model as radio-pulsars. We find that the cumulative distribution of the number of X-ray pulsars can be well reproduced by several models that, simultaneously, reproduce the characteristics of the radio-pulsar distribution. However, even considering the most favourable magneto-thermal evolution models with fast field decay, lognormal distributions of the initial magnetic field overpredict the number of visible sources with periods longer than 12 s. We then show that the problem can be solved with different distributions of magnetic field, such as a truncated lognormal distribution, or a binormal distribution with two distinct populations. We use the observational lack of isolated neutron stars (NSs) with spin periods P > 12 s to establish an upper limit to the fraction of magnetars born with B > 1015 G (less than 1 per cent). As future detections keep increasing the magnetar and high-B pulsar statistics, our approach can be used to establish a severe constraint on the maximum magnetic field at birth of NSs.

Are pulsars born with a hidden magnetic field?

The observation of several neutron stars in the centre 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 protoneutron 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 on to the neutron star as responsible for the submergence (or screening) of the field and its apparently low value. In this paper, we study under which conditions the magnetic field of a neutron star can be buried into the crust due to an accreting, conducting fluid. For this purpose, we consider a spherically symmetric calculation in general relativity to estimate the balance between the incoming accretion flow and the magnetosphere. Our study analyses several models with different specific entropy, composition, and neutron star masses. The main conclusion of our work is that typical magnetic fields of a few times 1012 G can be buried by accreting only 10−3–10−2 M⊙, a relatively modest amount of mass. In view of this result, the central compact object scenario should not be considered unusual, and we predict that anomalously weak magnetic fields should be common in very young (< few kyr) neutron stars.