ABSTRACT
Millisecond pulsars are mainly characterized by their spin periods, B-fields and masses - quantities that are largely affected by previous interactions with a companion star in a binary ...system. In this paper, we investigate the formation mechanism of millisecond pulsars by considering the pulsar recycling process in both intermediate-mass X-ray binaries (IMXBs) and low-mass X-ray binaries (LMXBs). The IMXBs mainly lead to the formation of binary millisecond pulsars with a massive carbon-oxygen (CO) or an oxygen-neon-magnesium white dwarf (ONeMg WD) companion, whereas the LMXBs form recycled pulsars with a helium white dwarf (He WD) companion. We discuss the accretion physics leading to the spin-up line in the -diagram and demonstrate that such a line cannot be uniquely defined. We derive a simple expression for the amount of accreted mass needed for any given pulsar to achieve its equilibrium spin and apply this to explain the observed differences of the spin distributions of recycled pulsars with different types of companions. From numerical calculations we present further evidence for significant loss of rotational energy in accreting X-ray millisecond pulsars in LMXBs during the Roche-lobe decoupling phase (Tauris 2012) and demonstrate that the same effect is negligible in IMXBs. We examine the recycling of pulsars with CO WD companions via Case BB Roche-lobe overflow (RLO) of naked helium stars in post-common envelope binaries. We find that such pulsars typically accrete of the order of 0.002-0.007 M⊙ which is just about sufficient to explain their observed spin periods. We introduce isochrones of radio millisecond pulsars in the -diagram to follow their spin evolution and discuss their true ages from comparison with observations. Finally, we apply our results of the spin-up process to complement our investigation of the massive pulsar PSR J1614−2230 from Paper I, confirming that this system formed via stable Case A RLO in an IMXB and enabling us to put new constraints on the birth masses of a number of recycled pulsars.
Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs) that have been spun up to high rotation rates via accretion of matter from a companion star in a low-mass X-ray binary ...(LMXB). Here we investigate binary evolution in close LMXBs to study the formation of radio MSPs with low-mass helium white dwarf companions (He WDs) in tight binaries with orbital periods Psuborbapprox =2-9 h. We find that severe fine-tuning is necessary to reproduce the observed MSPs in tight binaries with He WD companions of mass <0.20 M, which suggests that something needs to be modified or is missing in the standard input physics of LMXB modelling. The results of the thermal evolution of the (proto) He WDs are reported in a follow-up paper (Paper II).
A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either ...through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study for the first time the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities of Z = 0.02, 0.01, 0.001 and 0.0002, producing WDs with masses between ~ 0.16−0.45 M⊙. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. In addition, we confirm that the threshold mass at which this dichotomy occurs depends on metallicity. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances, i.e. the observed presence of metals in their atmospheres. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime. This is of great importance as helium is a crucial ingredient in the driving of the κ-mechanism suggested for the newly observed ELM proto-WD pulsators. However, we find that the number of hydrogen shell flashes and, as a result, the hydrogen envelope mass at the beginning of the cooling track, are not influenced significantly by rotational mixing. In addition to being dependent on proto-WD mass and metallicity, the hydrogen envelope mass of the newly formed proto-WDs depends on whether or not the donor star experiences a temporary contraction when the H-burning shell crosses the hydrogen discontinuity left behind by the convective envelope. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency.
Formation of Double Neutron Star Systems Tauris, T. M.; Kramer, M.; Freire, P. C. C. ...
Astrophysical journal/The Astrophysical journal,
09/2017, Letnik:
846, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers ...are anticipated to be discovered using the Square Kilometre Array searching for radio pulsars, and the high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase, and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most . We investigate DNS masses, spins, and velocities, and in particular correlations between spin period, orbital period, and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes that may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short γ-ray burst.
Context. Ultra-luminous X-ray sources (ULXs) are those X-ray sources located away from the centre of their host galaxy with luminosities exceeding the Eddington limit of a stellar-mass black hole ( L ...X > 10 39 erg s −1 ). Observed X-ray variability suggests that ULXs are X-ray binary systems. The discovery of X-ray pulsations in some of these objects (e.g. M82 X-2) suggests that a certain fraction of the ULX population may have a neutron star as the accretor. Aims. We present systematic modelling of low- and intermediate-mass X-ray binaries (LMXBs and IMXBs; donor-star mass range 0.92–8.0 M ⊙ and neutron-star accretors) to explain the formation of this sub-population of ULXs. Methods. Using MESA, we explored the allowed initial parameter space of binary systems consisting of a neutron star and a low- or intermediate-mass donor star that could explain the observed properties of ULXs. These donors are transferring mass at super-Eddington rates while the accretion is limited locally in the accretion disc by the Eddington limit. Thus, our simulations take into account beaming effects and also include stellar rotation, tides, general angular momentum losses, and a detailed and self-consistent calculation of the mass-transfer rate. Results. Exploring the initial parameters that lead to the formation of neutron-star ULXs, we study the conditions that lead to dynamical stability of these systems, which depends strongly on the response of the donor star to mass loss. Using two values for the initial neutron star mass (1.3 M ⊙ and 2.0 M ⊙ ), we present two sets of mass-transfer calculation grids for comparison with observations of NS ULXs. We find that LMXBs/IMXBs can produce NS-ULXs with typical time-averaged isotropic-equivalent X-ray luminosities of between 10 39 and 10 41 erg s −1 on a timescale of up to ∼1.0 Myr for the lower luminosities. Finally, we estimate their likelihood of detection, the types of white-dwarf remnants left behind by the donors, and the total amount of mass accreted by the neutron stars. Conclusions. We show that observed super-Eddington luminosities can be achieved in LMXBs/IMXBs undergoing non-conservative mass transfer while assuming geometrical beaming. We also compare our results to the observed pulsating ULXs and infer their initial parameters. Our results suggest that a large subset of the observed pulsating ULX population can be explained by LMXBs/IMXBs in a super-Eddington mass-transfer phase.
Context. Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs), formed via type Ib/c core-collapse supernovae (SNe), which have been spun up to high rotation rates via ...accretion from a companion star in a low-mass X-ray binary (LMXB). In an alternative formation channel, NSs are produced via the accretion-induced collapse (AIC) of a massive white dwarf (WD) in a close binary. Aims. Here we investigate binary evolution leading to AIC and examine if NSs formed in this way can subsequently be recycled to form MSPs and, if so, how they can observationally be distinguished from pulsars formed via the standard core-collapse SN channel in terms of their masses, spins, orbital periods and space velocities. Methods. Numerical calculations with a detailed stellar evolution code were used for the first time to study the combined pre- and post-AIC evolution of close binaries. We investigated the mass transfer onto a massive WD (treated as a point mass) in 240 systems with three different types of non-degenerate donor stars: main-sequence stars, red giants, and helium stars. When the WD is able to accrete sufficient mass (depending on the mass-transfer rate and the duration of the accretion phase) we assumed it collapses to form a NS and we studied the dynamical effects of this implosion on the binary orbit. Subsequently, we followed the mass-transfer epoch which resumes once the donor star refills its Roche lobe and calculated the continued LMXB evolution until the end. Results. We show that recycled pulsars may form via AIC from all three types of progenitor systems investigated and find that the final properties of the resulting MSPs are, in general, remarkably similar to those of MSPs formed via the standard core-collapse SN channel. However, as a consequence of the fine-tuned mass-transfer rate necessary to make the WD grow in mass, the resultant MSPs created via the AIC channel preferentially form in certain orbital period intervals. In addition, their predicted small space velocities can also be used to identify them observationally. The production time of NSs formed via AIC can exceed 10 Gyr which can therefore explain the existence of relatively young NSs in globular clusters. Our calculations are also applicable to progenitor binaries of SNe Ia under certain conditions.
Recent discoveries of weak and fast optical transients raise the question of their origin. We investigate the minimum ejecta mass associated with core-collapse supernovae (SNe) of Type Ic. We show ...that mass transfer from a helium star to a compact companion can produce an ultra-stripped core which undergoes iron core collapse and leads to an extremely fast and faint SN Ic. In this Letter, a detailed example is presented in which the pre-SN stellar mass is barely above the Chandrasekhar limit, resulting in the ejection of only ~0.05-0.20 M sub(middot in circle) of material and the formation of a low-mass neutron star (NS). We compute synthetic light curves of this case and demonstrate that SN 2005ek could be explained by our model. We estimate that the fraction of such ultra-stripped to all SNe could be as high as 10 super(?3)-10 super(?2). Finally, we argue that the second explosion in some double NS systems (for example, the double pulsar PSR J0737?3039B) was likely associated with an ultra-stripped SN Ic.
Context. The recently detected gravitational wave signals (GW150914 and GW151226) of the merger event of a pair of relatively massive stellar-mass black holes (BHs) calls for an investigation of the ...formation of such progenitor systems in general. Aims. We analyse the common-envelope (CE) stage of the traditional formation channel in binaries where the first-formed compact object undergoes an in-spiral inside the envelope of its evolved companion star and ejects the envelope in this process. Methods. We calculated envelope binding energies of donor stars with initial masses between 4 and 115M sub(middot in circle) for metallicities of Z= Z sub(Milky Way) = Z sub(middot in circle)/ 2 and Z= Z sub(middot in circle)/ 50, and derived minimum masses of in-spiralling objects needed to eject these envelopes. Results. In addition to producing double white dwarf and double neutron star binaries, CE evolution may also produce massive BH-BH systems with individual BH component masses of up to ~50 - 60M sub(middot in circle), in particular for donor stars evolved to giants beyond the Hertzsprung gap. However, the physics of envelope ejection of massive stars remains uncertain. We discuss the applicability of the energy-budget formalism, the location of the bifurcation point, the recombination energy, and the accretion energy during in-spiral as possible energy sources, and also comment on the effect of inflated helium cores. Conclusions. Massive stars in a wide range of metallicities and with initial masses of up to at least 115M sub(middot in circle) may shed their envelopes and survive CE evolution, depending on their initial orbital parameters, similarly to the situation for intermediate- and low-mass stars with degenerate cores. In addition to being dependent on stellar radius, the envelope binding energies and lambda-values also depend on the applied convective core-overshooting parameter, whereas these structure parameters are basically independent of metallicity for stars with initial masses below 60M sub(middot in circle). Metal-rich stars > or = 60M sub(middot in circle) become luminous blue variables and do not evolve to reach the red giant stage. We conclude that based on stellar structure calculations, and in the view of the usual simple energy budget analysis, events like GW150914 and GW151226 might be produced by the CE channel. Calculations of post-CE orbital separations, however, and thus the estimated LIGO detection rates, remain highly uncertain.
Abstract
We report the discovery of PSR J1757−1854, a 21.5-ms pulsar in a highly-eccentric, 4.4-h orbit with a neutron star (NS) companion. PSR J1757−1854 exhibits some of the most extreme ...relativistic parameters of any known pulsar, including the strongest relativistic effects due to gravitational-wave damping, with a merger time of 76 Myr. Following a 1.6-yr timing campaign, we have measured five post-Keplerian parameters, yielding the two component masses (mp = 1.3384(9) M⊙ and mc = 1.3946(9) M⊙) plus three tests of general relativity, which the theory passes. The larger mass of the NS companion provides important clues regarding the binary formation of PSR J1757−1854. With simulations suggesting 3-σ measurements of both the contribution of Lense–Thirring precession to the rate of change of the semimajor axis and the relativistic deformation of the orbit within ∼7–9 yr, PSR J1757−1854 stands out as a unique laboratory for new tests of gravitational theories.
Context. The origin and number of the Galactic supergiant X-ray binaries is currently not well understood. They consist of an evolved massive star and a neutron star or black-hole companion. X-rays ...are thought to be generated from the accretion of wind material donated by the supergiant, while mass transfer due to Roche-lobe overflow is mostly disregarded because the high mass ratios of these systems are thought to render this process unstable. Aims. We investigate how the proximity of supergiant donor stars to the Eddington limit, and their advanced evolutionary stage, may influence the evolution of massive and ultra-luminous X-ray binaries with supergiant donor stars (SGXBs and ULXs). Methods. We constructed models of massive stars with different internal hydrogen and helium gradients (H/He gradients) and different hydrogen-rich envelope masses, and exposed them to slow mass-loss to probe the response of the stellar radius. In addition, we computed the corresponding Roche-lobe overflow mass-transfer evolution with our detailed binary stellar evolution code, approximating the compact objects as point masses. Results. We find that a H/He gradient in the layers beneath the surface, as it is likely present in the well-studied donor stars of observed SGBXs, can enable mass transfer in SGXBs on a nuclear timescale with a black-hole or a neutron star accretor, even for mass ratios in excess of 20. In our binary evolution models, the donor stars rapidly decrease their thermal equilibrium radius and can therefore cope with the inevitably strong orbital contraction imposed by the high mass ratio. We find that the orbital period derivatives of our models agree well with empirical values. We argue that the SGXB phase may be preceded by a common-envelope evolution. The envelope inflation near the Eddington limit means that this mechanism more likely occurs at high metallicity. Conclusion. Our results open a new perspective for understanding that SGBXs are numerous in our Galaxy and are almost completely absent in the Small Magellanic Cloud. Our results may also offer a way to find more ULX systems, to detect mass transfer on nuclear timescales in ULX systems even with neutron star accretors, and shed new light on the origin of the strong B-field in these neutron stars.
PDF
