Search

PSR J1227$-$6208 is a 34.53-ms recycled pulsar with a massive companion. This system has long been suspected to belong to the emerging class of massive recycled pulsar-ONeMg white dwarf systems such as PSR J2222$-$0137, PSR J1528$-$3146 and J1439$-$5501. Here we present an updated emission and timing analysis with more than 11 years of combined Parkes and MeerKAT data, including 19 hours of high-frequency data from the newly installed MeerKAT S-band receivers. We measure a scattering timescale of 1.22 ms at 1 GHz with a flat scattering index 3.33<$\beta$<3.62, and a mean flux density of 0.53-0.62 mJy at 1 GHz with a steep spectral index 2.06<$\alpha$<2.35. Around 15% of the emission is linearly and circularly polarised, but the polarisation angle does not follow the rotating vector model. Thanks to the sensitivity of MeerKAT, we successfully measure a rate of periastron advance of 0.0171(11) deg/yr, and a Shapiro delay with an orthometric amplitude of 3.6$\pm$0.5 $\mu$s and an orthometric shape of 0.85$\pm$0.05. The main source of uncertainty in our timing analysis is chromatic correlated dispersion measure noise, which we model as a power law in the Fourier space thanks to the large frequency coverage provided by the Parkes UWL receiver. Assuming general relativity and accounting for the measurements across all the implemented timing noise models, the total mass, companion mass, pulsar mass and inclination angle are constrained at 2.3

The mergers of double neutron stars (DNSs) systems are believed to drive the majority of short $\gamma$-ray bursts (SGRBs), while also serving as production sites of heavy r-process elements. Despite being key to i) confirming the nature of the extragalactic SGRBs, ii) addressing the poorly-understood r-process enrichment in the ultra-faint dwarf galaxies (UFDGs), and iii) probing the formation process of DNS systems, the space velocity distribution of DNSs is still poorly constrained due to the small number of DNSs with well-determined astrometry. In this work, we determine new proper motions and parallaxes of two Galactic DNSs -- PSR J0509+3801 and PSR J1930-1852, using the Very Long Baseline Array, and estimate the transverse velocities $v_\perp$ of all the 11 isolated Galactic DNSs having proper motion measurements in a consistent manner. Our correlation analysis reveals that the DNS $v_\perp$ is tentatively correlated with three parameters: spin period, orbital eccentricity, and companion mass. With the preliminary $v_\perp$ distribution, we obtain the following findings. Firstly, the refined $v_\perp$ distribution is confirmed to agree with the observed displacements of the localized SGRBs from their host galaxy birth sites. Secondly, we estimate that around 11% and 25% of DNSs remain gravitationally bound to UFDGs with escape velocities of 15$\mathrm{~km~s^{-1}}$ and 25$\mathrm{~km~s^{-1}}$, respectively. Hence, the retained DNSs might indeed be responsible for the r-process enrichment confirmed so far in a few UFDGs. Finally, we discuss how a future ensemble of astrometrically determined DNSs may probe the multimodality of the $v_\perp$ distribution., Comment: 17 pages, 3 figures, 4 tables, accepted for publication in ApJ

We report discovery and characterization of a main-sequence G star orbiting a dark object with mass $1.90\pm 0.04 M_{\odot}$. The system was discovered via Gaia astrometry and has an orbital period of 731 days. We obtained multi-epoch RV follow-up over a period of 639 days, allowing us to refine the Gaia orbital solution and precisely constrain the masses of both components. The luminous star is a $\gtrsim 12$ Gyr-old, low-metallicity halo star near the main-sequence turnoff ($T_{\rm eff}\approx 6000$ K; $\log(g/\left[{\rm cm\,s^{-2}}\right])\approx 4.0$; $\rm [Fe/H]\approx-1.25$; $M\approx0.79 M_{\odot}$) with a highly enhanced lithium abundance. The RV mass function sets a minimum companion mass for an edge-on orbit of $M_2 > 1.67 M_{\odot}$, well above the Chandrasekhar limit. The Gaia inclination constraint, $i=68.7\pm 1.4$ deg, then implies a companion mass of $M_2=1.90\pm0.04 M_{\odot}$. The companion is most likely a massive neutron star: the only viable alternative is two massive white dwarfs in a close binary, but this scenario is disfavored on evolutionary grounds. The system's low eccentricity ($e=0.122\pm 0.002$) disfavors dynamical formation channels and implies that the neutron star likely formed with little mass loss ($\lesssim1\,M_{\odot}$) and with a weak natal kick ($v_{\rm kick}\lesssim 20\,\rm km\,s^{-1}$). The current orbit is too small to have accommodated the neutron star progenitor as a red supergiant or super-AGB star. The simplest formation scenario -- isolated binary evolution -- requires the system to have survived unstable mass transfer and common envelope evolution with a donor-to-accretor mass ratio $>10$. The system, which we call Gaia NS1, is likely a progenitor of symbiotic X-ray binaries and long-period millisecond pulsars. Its discovery challenges binary evolution models and bodes well for Gaia's census of compact objects in wide binaries., Comment: 19 pages, 12 figures, accepted to OJAp

Knowledge of the formation of neutron stars (NSs) in supernova (SN) explosions is of fundamental importance in wide areas of contemporary astrophysics: X-ray binaries, magnetars, radio pulsars, and, not least, double NS systems which merge and become gravitational wave sources. A recent study by Richardson et al. reported that the NS in the Be-star/X-ray binary SGR 0755-2933 (CPD -29 2176) descended from an ultra-stripped SN. Using the same observational data as Richardson et al., however, we find that the majority of progenitor solutions for SGR 0755-2933 are of normal Type Ib/c SNe, which allows for up to several solar masses of material to be ejected in the SN event. To correctly probe the SN explosion physics and inferring pre-SN conditions in a binary system, a full kinematic analysis based on post-SN data is always needed., Comment: Nature, accepted for publication (10 pages, 5 figures, incl. Extended Data), Matters Arising, Authors' version

The recently reported observation of VFTS 243 is the first example of a massive black-hole binary system with negligible binary interaction following black-hole formation. The black-hole mass ($\approx 10\ M_{\odot}$) and near-circular orbit ($e\approx 0.02$) of VFTS 243 suggest that the progenitor star experienced complete collapse, with energy-momentum being lost predominantly through neutrinos. VFTS 243 enables us to constrain the natal kick and neutrino-emission asymmetry during black-hole formation. At 68% C.L., the natal kick velocity (mass decrement) is $\lesssim 10$ km/s ($\lesssim 1.0\ M_{\odot}$), with a full probability distribution that peaks when $\approx 0.3\ M_{\odot}$ were ejected, presumably in neutrinos, and the black hole experienced a natal kick of $4$ km/s. The neutrino-emission asymmetry is $\lesssim 4$%, with best fit values of $\sim$0-0.2%. Such a small neutrino natal kick accompanying black-hole formation is in agreement with theoretical predictions., Comment: 15 pages (including supplemental material), 2 main figures and 5 supplemental figures. Accepted for publication in PRL

The majority of all stars are members of a binary system. The evolution of such binary stars and their subsequent production of pairs of compact objects in tight orbits, such as double neutron stars and double black holes, play a central role in modern astrophysics, Binary evolution leads to the formation of different types of violent cosmic events such as novae, supernova explosions, gamma-ray bursts, mass transfer and accretion processes in X-ray binaries, and the formation of exotic radio millisecond pulsars. In some cases, the binary systems terminate as spectacular collisions between neutron stars and/or black holes. These collisions lead to powerful emission of gravitational waves, as detected by LIGO since 2015. The coming decade is expected to reveal a large number of discoveries of binary compact systems, as well as their progenitors and merger remnants, from major instruments such as the radio Square-Kilometre Array; the gravitational wave observatories LIGO-Virgo-KAGRA-IndIGO and LISA; the astrometric space observatory Gaia; the James Webb Space Telescope; and the X-ray space observatories eXTP, STROBE-X, and Athena. In this light, it is important to have a modern textbook on the physics of binary stars evolution, from ordinary stars to X-ray binaries and gravitational wave sources. The scope of this book is that the reader (student or educated expert) will learn the physics of binary interactions, from stellar birth to compact objects, and relate this knowledge to the latest observations. The reader will learn about stellar structure and evolution, and detailed binary interactions covering a broad range of phenomena, including mass transfer and orbital evolution, formation and accretion onto compact objects (white dwarfs, neutron stars and black holes), and their observational properties. Exercises are provided throughout the book., Comment: Summary of the 864-pages textbook "Physics of Binary Star Evolution - From Stars to X-ray Binaries and Gravitational Wave Sources", by Thomas Tauris and Ed van den Heuvel, to be published by Princeton University Press in June 2023

The formation path to ultra-compact X-ray binaries (UCXBs) with black hole (BH) accretors is still unclear. In the classical formation scenario, it is difficult to eject the massive envelope of the progenitor star of the BH via common envelope process. Given that some neutron stars (NSs) in binary systems evidently have birth masses close to $\sim 2.0\;M_\odot$, we explore here the possibility that BH-UCXBs may form via accretion-induced collapse (AIC) of accreting NSs, assuming that these previously evolved in LMXBs to masses all the way up to the maximum limit of a NS. We demonstrate this formation path by modelling a few cases of NS-UCXBs with initial NS masses close to the maximum mass of a NS that evolve into BH-UCXBs after the NS accretes material from its He~WD companion. We follow the evolution of the post-AIC BH-UCXB and, based on simple arguments, we anticipate that there is about one BH-UCXB with an AIC origin and a He~WD donor within the current sample of known UCXBs and that 2--5 such BH-UCXBs may be detected in gravitational waves by LISA. In addition, we find that the X-ray luminosity of NS-UCXBs near their orbital period minimum exceeds $\sim 10^{39}\;{\rm erg\;s^{-1}}$ and thus such systems may appear as ultraluminous X-ray sources., Comment: ApJ accepted, 12 pages, 3 figures

Recent discoveries of gravitational wave (GW) events most likely originating from black hole (BH) + neutron star (NS) mergers reveal the existence of BH+NS binaries. The formation of BH+NS binaries and their merger rates through isolated binary evolution have been investigated extensively with population synthesis simulations. A detailed stellar evolution modelings of the formation of this population, however, is missing in the literature. In this work, we perform the first complete 1D model of more than 30 BH+NS progenitor systems which are calculated self-consistently until the iron core collapse with infall velocity exceeds 1000 km s^-1. Focusing on the progenitors of BH- NS GW sources, we apply the MESA code starting from a post-common envelope binary with short orbital period (< 1 day) consisting of a BH and a zero-age main-sequence helium star that experiences stable mass transfer. These NS masses could be significantly larger depending on the exact mass cut during the supernova explosion. These BH+NS systems are likely to merge and produce GW events within a Hubble time. System C is a potential progenitor of a GW200115-like event, while Systems A and B are possible candidates for a GW200105-like event and may represent the final destiny of the X-ray binary SS433., Comment: 14pages, 9 figures, Accepted to be published on ApJ

The detection of double black hole (BH+BH) mergers provides a unique possibility to understand their physical properties and origin. To date, the LIGO-Virgo-KAGRA network of high-frequency gravitational wave observatories have announced the detection of more than 85 BH+BH merger events (Abbott et al. 2022a). An important diagnostic feature that can be extracted from the data is the distribution of effective inspiral spins of the BHs. This distribution is in clear tension with theoretical expectations from both an isolated binary star origin, which traditionally predicts close-to aligned BH component spins (Kalogera 2000; Farr et al. 2017), and formation via dynamical interactions in dense stellar environments that predicts a symmetric distribution of effective inspiral spins (Mandel & O'Shaughnessy 2010; Rodriguez et al. 2016b). Here it is demonstrated that isolated binary evolution can convincingly explain the observed data if BHs have their spin axis tossed during their formation process in the core collapse of a massive star, similarly to the process evidently acting in newborn neutron stars. BH formation without spin-axis tossing, however, has difficulties reproducing the observed data - even if alignment of spins prior to the second core collapse is disregarded. Based on simulations with only a minimum of assumptions, constrains from empirical data can be made on the spin magnitudes of the first- and second-born BHs, thereby serving to better understand massive binary star evolution prior to the formation of BHs., Comment: 18 pages, 11 figures, 1 table, ApJ accepted version (minor edits)

Neutron star-helium white dwarf (NS+He WD) binaries are important evolutionary products of close-orbit binary star systems. They are often observed as millisecond pulsars and may continue evolving into ultra-compact X-ray binaries (UCXBs) and continuous gravitational wave (GW) sources that will be detected by space-borne GW observatories, such as LISA, TianQin and Taiji. Nevertheless, the stability of NS+He WD binaries undergoing mass transfer is not well studied and still under debate. In this paper, we model the evolution of NS+He WD binaries with WD masses ranging from 0.17-0.45 $M_{\odot}$, applying the detailed stellar evolution code mesa. Contrary to previous studies based on hydrodynamics, we find that apparently all NS+He WD binaries undergo stable mass transfer. We find for such UCXBs that the larger the WD mass, the larger the maximum mass-transfer rate and the smaller the minimum orbital period during their evolution. Finally, we demonstrate numerically and analytically that there is a tight correlation between WD mass and GW frequency for UCXBs, independent of NS mass., Comment: Accepted by ApJ, 9 pages, 7 figures

Double white dwarf (DWD) binaries are important for studies of common-envelope (CE) evolution, Type Ia supernova progenitors and Galactic sources of low-frequency gravitational waves (GWs). PTF J0533+0209 is a DWD system with a short orbital period of Porb ~ 20 min and thus a so-called LISA verification source. The formation of this system and other DWDs is still under debate. In this paper, we discuss the possible formation scenarios of this binary and argue that it is not likely to have formed through CE evolution. Applying a new magnetic braking prescription, we use the MESA code to model the formation of this system through stable mass transfer. We find a model which can well reproduce the observed WD masses and orbital period but not the effective temperature and hydrogen abundance of the low-mass He WD component. We discuss the possibility of using H flashes to mitigate this discrepancy. Finally, we discuss the future evolution of this system into a AM CVn binary such as those that will be detected by space-borne GW observatories like LISA, TianQin and Taiji., Comment: Accepted by ApJ, 11 pages, 4 figures

Recent discoveries of double neutron star (DNS) mergers and ultra-stripped supernovae (SNe) raise the questions of their origin and connection. We present the first 1D~model of a DNS progenitor system which is calculated self-consistently until an ultra-stripped iron core collapse. We apply the \texttt{MESA} code starting from a post-common envelope binary consisting of a $1.35\;M_\odot$ NS and a $3.20\;M_\odot$ zero-age main-sequence helium star and continue the modelling via Case~BB Roche-lobe overflow until the infall velocity of the collapsing iron core exceeds $1000\;{\rm km\,s}^{-1}$. The exploding star has a total mass of $\sim 1.90\;M_\odot$, consisting of a $\sim 0.29\;M_\odot$ He-rich envelope embedding a CO core of $\sim 1.61\;M_\odot$ and an iron-rich core of $\sim 1.50\;M_\odot$. The resulting second-born NS has an estimated mass of $\sim 1.44\;M_\odot$ and we discuss the fate of the post-SN system, as well as the mild recycling of the first-born NS. Depending on the initial conditions, this family of systems is anticipated to reproduce the DNS mergers detected by the LIGO-network., Comment: 8 pages, 5 figures, ApJL in press

The first three observing runs with Advanced LIGO and Virgo have resulted in the detection of binary black hole mergers (BBH) with highly unequal mass components, which are difficult to reconcile with standard formation paradigms. The most representative of these is GW190814, a highly asymmetric merger between a 23 M$_{\odot}$ black hole and a 2.6 M$_{\odot}$ compact object. Here, we explore recent results suggesting that a sizeable fraction of stars with pre-collapse carbon-oxygen core masses above 10 M$_{\odot}$, and extending up to at least 30 M$_{\odot}$, may produce objects inside the so-called lower mass gap that bridges the division between massive pulsars and BHs in Galactic X-ray binaries. We demonstrate that such an explosion landscape would naturally cause a fraction of massive binaries to produce GW190814-like systems instead of symmetric-mass BBHs. We present examples of specific evolutionary channels leading to the formation of GW190814 and GW200210, a 24+2.8 M$_{\odot}$ merger discovered during the O3b observing run. We estimate the merger-rate density of these events in our scenario to be $\mathcal{O}$(5%) of the total BBH merger rate. Finally, we discuss the broader implications of this formation channel for compact object populations, and its possible relevance to less asymmetric merger events such as GW200105 and GW200115, Comment: A&A Letters, in press

The origin and fate of magnetars (young, extremely magnetized neutron stars, NSs) remain unsolved. Probing their evolution is therefore crucial for investigating possible links to other species of isolated NSs, such as the X-ray dim NSs (XDINSs) and rotating radio transients (RRATs). Here we investigate the spin evolution of magnetars. Two avenues of evolution are considered: one with exponentially decaying B-fields, the other with sub- and super-exponential decay. Using Monte Carlo methods, we synthesize magnetar populations using different input distributions and physical parameters, such as for the initial spin period, its time derivative and the B-field decay timescale. Additionally, we introduce a fade-away procedure that can account for the fading of old magnetars, and we briefly discuss the effect of alignment of the B-field and spin axes. Imposing the Galactic core-collapse supernova rate of ~20 kyr^-1 as a strict upper limit on the magnetar birthrate and comparing the synthetic populations to the observed one using both manual and automatic optimization algorithms for our input parameter study, we find that the B-field must decay exponentially or super-exponentially with a characteristic decay timescale of 0.5-10 kyr (with a best value of ~4 kyr). In addition, the initial spin period must be less than 2 sec. If these constraints are kept, we conclude that there are multiple choices of input physics that can reproduce the observed magnetar population reasonably well. We also conclude that magnetars may well be evolutionary linked to the population of XDINSs, whereas they are in general unlikely to evolve into RRATs., Comment: 18 pages, 11 figures, +appendices, accepted for publication in MNRAS

Close-orbit low-mass X-ray binaries (LMXBs), radio binary millisecond pulsars (BMSPs) with extremely low-mass helium WDs (ELM He~WDs) and ultra-compact X-ray binaries (UCXBs) are all part of the same evolutionary sequence. It is therefore of uttermost importance to understand how these populations evolve from one specie to another. Moreover, UCXBs are important gravitational wave (GW) sources and can be detected by future space-borne GW observatories. However, the formation and evolutionary link between these three different populations of neutron star (NS) binaries are not fully understood. In particular, a peculiar fine-tuning problem has previously been demonstrated for the formation of these systems. In this investigation, we test a newly suggested magnetic braking prescription and model the formation and evolution of LMXBs. We compute a grid of binary evolution models and present the initial parameter space of the progenitor binaries which successfully evolve all the way to produce UCXBs. We find that the initial orbital period range of LMXBs, which evolve into detached NS+ELM~He~WD binaries and later UCXBs, becomes significantly wider compared to evolution with a standard magnetic braking prescription, and thus helps to relieve the fine-tuning problem. However, we also find that formation of wide-orbit BMSPs is prohibited for strong versions of this new magnetic braking prescription, which therefore calls for a revision of the prescription. Finally, we present examples of the properties of UCXBs as Galactic GW sources and discuss their detection by the LISA, TianQin and Taiji observatories., Comment: 12 pages, 13 figures, accepted for publication in MNRAS

Thompson et al. (Reports, 1 November 2019, p. 637, Science) interpreted the unseen companion of the red giant star 2MASS J05215658+4359220 as most likely a black hole. We argue that if the red giant is about one solar mass, its companion can be a close binary consisting of two main-sequence stars. This would explain why no X-ray emission is detected from the system., Comment: 3 pages, Author version of Technical Comment published in Science on 8 May, 2020

The application of standard accretion theory to observations of X-ray binaries provides valuable insights into neutron star properties, such as their spin period and magnetic field. However, most studies concentrate on relatively old systems, where the neutron star is in its late propeller, accretor, or nearly spin equilibrium phase. Here we use an analytic model from standard accretion theory to illustrate the evolution of high-mass X-ray binaries early in their life. We show that a young neutron star is unlikely to be an accretor because of the long duration of ejector and propeller phases. We apply the model to the recently discovered ~4000 yr old high-mass X-ray binary XMMU J051342.6-672412 and find that the system's neutron star, with a tentative spin period of 4.4 s, cannot be in the accretor phase and has a magnetic field B > (a few)x10^13 G, which is comparable to the magnetic field of many older high-mass X-ray binaries and is much higher than the spin equilibrium inferred value of (a few)x10^11 G. The observed X-ray luminosity could be the result of thermal emission from a young cooling magnetic neutron star or a small amount of accretion that can occur in the propeller phase., Comment: 7 pages, 4 figures; accepted for publication in MNRAS

The compact binary radio pulsar system J0453+1559 (Martinez et al. 2015) consists of a recycled pulsar as primary component of $1.559(5)\,M_\odot$ and an unseen companion star of $1.174(4)\,M_\odot$. Because of the relatively large orbital eccentricity of $e = 0.1125$, it was argued that the companion is a neutron star, making it the neutron star with the lowest accurately determined mass to date. However, a direct observational determination of the nature of the companion is currently not feasible. Moreover, state-of-the-art stellar evolution and supernova modeling are contradictive concerning the possibility to produce such a low-mass neutron star remnant. Here we challenge the neutron star interpretation by reasoning that the lower-mass component could instead be a white dwarf born in a thermonuclear electron-capture supernova (tECSN) event, in which oxygen-neon deflagration in the degenerate stellar core of an ultra-stripped progenitor ejects several $0.1\,M_\odot$ of matter and leaves a bound ONeFe white dwarf as the second-formed compact remnant. We determine the ejecta mass and remnant kick needed in this scenario to explain the properties of PSR~J0453+1559 by a neutron star--white dwarf system. More work on tECSNe is needed to assess the viability of this scenario., Comment: ApJL, accepted for publication (6 pages, 4 figures)

Millisecond pulsars are old, fast spinning neutron stars thought to have evolved from classical pulsars in binary systems, where the rapid rotation is caused by the accretion of matter and angular momentum from their companion. During this transition between classical and millisecond pulsars, there is a magnetic field reduction of around 4 orders of magnitude, which is not well understood. According to the standard scenario, the magnetic field is reduced as a consequence of accretion, either through ohmic dissipation or through screening by the accreted matter. We explored an alternative hypothesis in which the magnetic field is reduced through ambipolar diffusion before the accretion. This is particularly effective during the long epoch in which the pulsar has cooled, but has not yet started accreting. This makes the final magnetic field dependent on the evolution time of the companion star and thus its initial mass. We use observed binary systems to constrain the time available for the magnetic field decay based on the current pulsar companion: a helium white dwarf, a carbon-oxygen white dwarf, or another neutron star. Based on a simplified model without baryon pairing, we show that the proposed process agrees with the general distribution of observed magnetic field strengths in binaries, but is not able to explain some mildly recycled pulsars where no significant decay appears to have occurred. We discuss the possibility of other formation channels for these systems and the conditions under which the magnetic field evolution would be set by the neutron star crust rather than the core., Comment: 10 pages, 4 figures, 2 tables

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). We construct models of massive stars with different internal hydrogen/helium gradients and different hydrogen-rich envelope masses, and expose them to slow mass loss to probe the response of the stellar radius. In addition, we compute the corresponding Roche-lobe overflow mass-transfer evolution with our detailed binary stellar evolution code, approximating the compact objects as point masses. We find that a hydrogen/helium gradient in the layers beneath the surface, as it is likely present in the well-studied donor stars of observed SGBXs, can enable nuclear timescale mass-transfer in SGXBs with a BH or a NS 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. Our results open a new perspective for understanding the large number of Galactic SGXBs, and their almost complete absence in the SMC. They may also offer a way to obtain more ULX systems, to find nuclear timescale mass-transfer in ULX systems even with neutron star accretors, and shed new light on the origin of the strong B-field in these neutron stars., Comment: 23 pages, 21 figures, we are thankful for any comments an this draft

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

The prime candidate sources for the upcoming space-borne gravitational wave (GW) observatory LISA are the numerous Galactic tight binaries of white dwarfs (WDs) and neutron stars (NSs), many of which will coalesce and undergo mass transfer, leading to simultaneous emission of X-rays and GWs. Here, detailed and coherent numerical stellar models are explored for the formation and evolution of these systems, including finite-temperature effects and complete calculations of mass transfer from a WD to a NS accretor. Evolutionary tracks of characteristic strain amplitude are computed, and the unique pattern of their evolution in the GW frequency-dynamical chirp mass parameter space enables a firm identification of the nature of the systems. Furthermore, it is demonstrated that a precise detection of the chirp allows determination of the NS mass to an accuracy of a few per cent, with applications to constraining its equation-of-state, in particular for dual-line GW sources observed simultaneously at high and low frequencies., Comment: Published in PRL. This version includes corrected Eq.(7) published in Erratum

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

The first gravitational wave detections of mergers between black holes and neutron stars represent a remarkable new regime of high-energy transient astrophysics. The signals observed with LIGO-Virgo detectors come from mergers of extreme physical objects which are the end products of stellar evolution in close binary systems. To better understand their origin and merger rates, we have performed binary population syntheses at different metallicities using the new grid-based binary population synthesis code ComBinE. Starting from newborn pairs of stars, we follow their evolution including mass loss, mass transfer and accretion, common envelopes and supernova explosions. We apply the binding energies of common envelopes based on dense grids of detailed stellar structure models, make use of improved investigations of the subsequent Case BB Roche-lobe overflow and scale supernova kicks according to the stripping of the exploding stars. We demonstrate that all the double black hole mergers, GW150914, LVT151012, GW151226, GW170104, GW170608 and GW170814, as well as the double neutron star merger GW170817, are accounted for in our models in the appropriate metallicity regime. Our binary interaction parameters are calibrated to match the accurately determined properties of Galactic double neutron star systems, and we discuss their masses and types of supernova origin. Using our default values for the input physics parameters, we find a double neutron star merger rate of about 3.0 Myr^-1 for Milky-Way equivalent galaxies. Our upper limit to the merger-rate density of double neutron stars is R=400 yr^-1 Gpc^-3 in the local Universe (z=0)., Comment: 36 pages, 26 figures, 8 tables, plus 9 pages appendix. Accepted 2018 August 6 by MNRAS after revision according to referee report (in particular, including further discussions on the progenitor binary of GW170817)

Previous optical and radio observations of the binary millisecond pulsar PSR J1640+2224 have come to inconsistent conclusions about the identity of its companion, with some observations suggesting the companion is a low-mass helium-core (He-core) white dwarf (WD), while others indicate it is most likely a high-mass carbon-oxygen (CO) WD. Binary evolution models predict PSR J1640+2224 most likely formed in a low-mass X-ray binary (LMXB) based on the pulsar's short spin period and long-period, low-eccentricity orbit, in which case its companion should be a He-core WD with mass about $0.35 - 0.39 \, M_\odot$, depending on metallicity. If it is instead a CO WD, that would suggest the system has an unusual formation history. In this paper we present the first astrometric parallax measurement for this system from observations made with the Very Long Baseline Array (VLBA), from which we determine the distance to be $1520^{+170}_{-150}\,\mathrm{pc}$. We use this distance and a reanalysis of archival optical observations originally taken in 1995 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST) in order to measure the WD's mass. We also incorporate improvements in calibration, extinction model, and WD cooling models. We find that the existing observations are not sufficient to tightly constrain the companion mass, but we conclude the WD mass is $>0.4\,M_\odot$ with $>90\%$ confidence. The limiting factor in our analysis is the low signal-to-noise ratio of the original HST observations., Comment: 6 pages, 5 figures

We report on the results of a 4-year timing campaign of PSR~J2222$-0137$, a 2.44-day binary pulsar with a massive white dwarf (WD) companion, with the Nan\c{c}ay, Effelsberg and Lovell radio telescopes. Using the Shapiro delay for this system, we find a pulsar mass $m_{p}=1.76,\pm\,0.06,M_\odot$ and a WD mass $m_{c}\,=\,1.293\,\pm\,0.025\, M_\odot$. We also measure the rate of advance of periastron for this system, which is marginally consistent with the GR prediction for these masses. The short lifetime of the massive WD progenitor star led to a rapid X-ray binary phase with little ($< \, 10^{-2} \, M_\odot$) mass accretion onto the neutron star (NS); hence, the current pulsar mass is, within uncertainties, its birth mass; the largest measured to date. We discuss the discrepancy with previous mass measurements for this system; we conclude that the measurements presented here are likely to be more accurate. Finally, we highlight the usefulness of this system for testing alternative theories of gravity by tightly constraining the presence of dipolar radiation. This is of particular importance for certain aspects of strong-field gravity, like spontaneous scalarization, since the mass of PSR~J2222$-0137$ puts that system into a poorly tested parameter range., Comment: 17 pages in journal format, 3 figures, 2 tables, accepted for publication in the Astrophysical Journal. Some references fixed and updated

Gravitational waves from the binary black hole (BH) merger GW150914 may enlighten our understanding of ultra-luminous X-ray sources (ULXs), as BHs>30Msun can reach luminosities>4x10^39 erg s^-1 without exceeding their Eddington limit. It is then important to study variations of evolutionary channels for merging BHs, which might instead form accreting BHs and become ULXs. It was recently shown that massive binaries with mass ratios close to unity and tight orbits can undergo efficient rotational mixing and evolve chemically homogeneously, resulting in a compact BH binary. We study similar systems by computing ~120000 detailed binary models with the MESA code covering a wide range of initial parameters. For initial mass ratios M2/M1~0.1-0.4, primaries >40Msun can evolve chemically homogeneously, remaining compact and forming a BH without undergoing Roche-lobe overflow. The secondary then expands and transfers mass to the BH, initiating a ULX phase. We predict that ~1 out of 10^4 massive stars evolves this way, and that in the local universe 0.13 ULXs per Msun yr^-1 of star-formation rate are observable, with a strong preference for low-metallicities. At metallicities log Z>-3, BH masses in ULXs are limited to 60Msun due to the occurrence of pair-instability supernovae which leave no remnant, resulting in an X-ray luminosity cut-off. At lower metallicities, very massive stars can avoid exploding as pair-instability supernovae and instead form BHs with masses above 130Msun, producing a gap in the ULX luminosity distribution. After the ULX phase, neutron-star-BH binaries that merge in less than a Hubble time are produced with a formation rate <0.2 Gpc^-3 yr^-1. We expect that upcoming X-ray observatories will test these predictions, which together with additional gravitational wave detections will provide strict constraints on the origin of the most massive BHs that can be produced by stars., Comment: Accepted for publication in A&A. 19 Pages plus 16 pages of appendices. Abstract abridged

Tight binaries of helium white dwarfs (He WDs) orbiting millisecond pulsars (MSPs) will eventually "merge" due to gravitational damping of the orbit. The outcome has been predicted to be the production of long-lived ultra-compact X-ray binaries (UCXBs), in which the WD transfers material to the accreting neutron star (NS). Here we present complete numerical computations, for the first time, of such stable mass transfer from a He WD to a NS. We have calculated a number of complete binary stellar evolution tracks, starting from pre-LMXB systems, and evolved these to detached MSP+WD systems and further on to UCXBs. The minimum orbital period is found to be as short as 5.6 minutes. We followed the subsequent widening of the systems until the donor stars become planets with a mass of ~0.005 Msun after roughly a Hubble time. Our models are able to explain the properties of observed UCXBs with high helium abundances and we can identify these sources on the ascending or descending branch in a diagram displaying mass-transfer rate vs. orbital period., Comment: 6 pages, 4 figures, MNRAS Letters, in press

We investigate light-curve and spectral properties of ultra-stripped core-collapse supernovae. Ultra-stripped supernovae are the explosions of heavily stripped massive stars which lost their envelopes via binary interactions with a compact companion star. They eject only ~ 0.1 Msun and may be the main way to form double neutron-star systems which eventually merge emitting strong gravitational waves. We follow the evolution of an ultra-stripped supernova progenitor until iron core collapse and perform explosive nucleosynthesis calculations. We then synthesize light curves and spectra of ultra-stripped supernovae using the nucleosynthesis results and present their expected properties. Ultra-stripped supernovae synthesize ~ 0.01 Msun of radioactive 56Ni, and their typical peak luminosity is around 1e42 erg/s or -16 mag. Their typical rise time is 5 - 10 days. Comparing synthesized and observed spectra, we find that SN 2005ek, some of the so-called calcium-rich gap transients, and SN 2010X may be related to ultra-stripped supernovae. If these supernovae are actually ultra-stripped supernovae, their event rate is expected to be about 1 per cent of core-collapse supernovae. Comparing the double neutron-star merger rate obtained by future gravitational-wave observations and the ultra-stripped supernova rate obtained by optical transient surveys identified with our synthesized light-curve and spectral models, we will be able to judge whether ultra-stripped supernovae are actually a major contributor to the binary neutron star population and provide constraints on binary stellar evolution., Comment: 15 pages, 18 figures, 2 tables, accepted by Monthly Notices of the Royal Astronomical Society

Recycled pulsars are mainly characterized by their spin periods, B-fields and masses. All these quantities are affected by previous interactions with a companion star in a binary system. Therefore, we can use these quantities as fossil records and learn about binary evolution. Here, I briefly review the distribution of these observed quantities and summarize our current understanding of the pulsar recycling process., Comment: Brief summary of invited review talk @ MODEST-16. 4 pages, 3 figures. To appear in: "Cosmic-Lab: Star Clusters as Cosmic Laboratories for Astrophysics, Dynamics and Fundamental Physics", F.R Ferraro & B. Lanzoni eds, Mem. SAIt, Vol 87, n.4

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 between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. 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. 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. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. 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., Comment: 28 pages, 10 pages tables, 19 figures, accepted for publication in A&A

The mass function of neutron stars (NSs) contains information about the late evolution of massive stars, the supernova explosion mechanism, and the equation-of-state of cold, nuclear matter beyond the nuclear saturation density. A number of recent NS mass measurements in binary millisecond pulsar (MSP) systems increase the fraction of massive NSs (with $M > 1.8$ M$_{\odot}$) to $\sim 20\% $ of the observed population. In light of these results, we employ a Bayesian framework to revisit the MSP mass distribution. We find that a single Gaussian model does not sufficiently describe the observed population. We test alternative empirical models and infer that the MSP mass distribution is strongly asymmetric. The diversity in spin and orbital properties of high-mass NSs suggests that this is most likely not a result of the recycling process, but rather reflects differences in the NS birth masses. The asymmetry is best accounted for by a bimodal distribution with a low mass component centred at $1.393_{-0.029}^{+0.031}$ M$_{\odot}$ and dispersed by $0.064_{-0.025}^{+0.064}$ M$_{\odot}$, and a high-mass component with a mean of $1.807_{-0.132}^{+0.081}$ and a dispersion of $0.177_{-0.072}^{+0.115}$ M$_{\odot}$. We also establish a lower limit of $M_{max} \ge 2.018$ M$_{\odot}$ at 98% C.L. for the maximum NS mass, from the absence of a high-mass truncation in the observed masses. Using our inferred model, we find that the measurement of 350 MSP masses, expected after the conclusion of pulsar surveys with the Square-Kilometre Array, can result in a precise localization of a maximum mass up to 2.15 M$_{\odot}$, with a 5% accuracy. Finally, we identify possible massive NSs within the known pulsar population and discuss birth masses of MSPs., Comment: submitted to ApJ; 21 pages in aastex6 two-column format, 12 figures, 5 tables. Comments are welcome

Superluminous supernovae (SLSNe) have been suggested to be powered by strongly magnetized, rapidly rotating neutron stars which are often called magnetars. In this process, rotational energy of the magnetar is radiated via magnetic dipole radiation and heats the supernova ejecta. However, if magnetars are highly distorted in their geometric shape, rotational energy is mainly lost as gravitational wave radiation and thus such magnetars cannot power SLSNe. By simply comparing electromagnetic and gravitational wave emission timescales, we constrain upper limits to the ellipticity of magnetars by assuming that they power the observed SLSNe. We find that their ellipticity typically needs to be less than about a few 1e-3. This indicates that the toroidal magnetic field strengths in these magnetars are typically less than a few 1e16 G so that their distortions remain small. Because light-curve modelling of SLSNe shows that their dipole magnetic field strengths are of the order of 1e14 G, the ratio of poloidal to toroidal magnetic field strengths is found to be larger than ~ 0.01 in magnetars powering SLSNe., Comment: 5 pages, 1 figure, 1 table, accepted by Monthly Notices of the Royal Astronomical Society Letters

The progenitors of many CCSNe are expected to be in binary systems. After the SN explosion, the companion may suffer from mass stripping and be shock heated as a result of the impact of the SN ejecta. If the binary system is disrupted, the companion is ejected as a runaway and hypervelocity star. By performing a series of 3D hydrodynamical simulations of the collision of SN ejecta with the companion star, we investigate how CCSN explosions affect their companions. We use the BEC code to construct the detailed companion structure at the time of SN explosion. The impact of the SN blast wave on the companion is followed by means of 3D SPH simulations using the Stellar GADGET code. For main-sequence (MS) companions, we find that the amount of removed mass, impact velocity, and chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Their relationship can be approximately fitted by power laws, which is consistent with the results obtained from impact simulations of SNe~Ia. However, we find that the impact velocity is sensitive to the momentum profile of the outer SN ejecta and, in fact, may decrease with increasing ejecta mass, depending on the modeling of the ejecta. Because most companions to Ib/c CCSNe are in their MS phase at the moment of the explosion, combined with the strongly decaying impact effects with increasing binary separation, we argue that the majority of these SNe lead to inefficient mass stripping and shock heating of the companion star following the impact of the ejecta. Our simulations show that the impact effects of Ib/c SN ejecta on the structure of MS companions, and thus their long-term post-explosion evolution, is in general not dramatic. We find that at most 10% of their mass is lost, and their resulting impact velocities are less than 100 km/s., Comment: Accepted for publication in Astronomy and Astrophysics, some minor typographical errors are fixed, the affiliation of second author is corrected

The explosion of ultra-stripped stars in close binaries can lead to ejecta masses < 0.1 M_sun and may explain some of the recent discoveries of weak and fast optical transients. In Tauris et al. (2013), it was demonstrated that helium star companions to neutron stars (NSs) may experience mass transfer and evolve into naked ~1.5 M_sun metal cores, barely above the Chandrasekhar mass limit. Here we present a systematic investigation of the progenitor evolution leading to ultra-stripped supernovae (SNe). In particular, we examine the binary parameter space leading to electron-capture (EC SNe) and iron core-collapse SNe (Fe CCSNe), respectively, and determine the amount of helium ejected with applications to their observational classification as Type Ib or Type Ic. We mainly evolve systems where the SN progenitors are helium star donors of initial mass M_He = 2.5 - 3.5 M_sun in tight binaries with orbital periods of P_orb = 0.06 - 2.0 days, and hosting an accreting NS, but we also discuss the evolution of wider systems and of both more massive and lighter - as well as single - helium stars. In some cases we are able to follow the evolution until the onset of silicon burning, just a few days prior to the SN explosion. We find that ultra-stripped SNe are possible for both EC SNe and Fe CCSNe, and that the amount of helium ejected is correlated with P_orb - the tightest systems even having donors being stripped down to envelopes of less than 0.01 M_sun. We estimate the rise time of ultra-stripped SNe to be in the range 12 hr - 8 days, and light curve decay times between 1 and 50 days. Ultra-stripped SNe may produce NSs in the mass range 1.10 - 1.80 M_sun and are highly relevant for LIGO/VIRGO since most (possibly all) merging double NS systems have evolved through this phase. Finally, we discuss the low-velocity kicks which might be imparted on these resulting NSs at birth. [Abridged], Comment: 22 pages, 20 figures, 2 tables, MNRAS in press, refs. updated

In this Habilitationsschrift (Habilitation thesis) I present my research carried out over the last four years at the Argelander Institute for Astronomy (AIfA) and the Max Planck Institute for Radio Astronomy (MPIfR). The thesis summarizes my main findings and has been written to fulfill the requirements for the Habilitation qualification at the University of Bonn. Although my work is mainly focused on the topic of millisecond pulsars (MSPs), there is a fairly broad spread of research areas ranging from the formation of neutron stars (NSs) in various supernova (SN) events, to their evolution, for example, via accretion processes in binary and triple systems, and finally to their possible destruction in merger events. The thesis is organized in the following manner: A general introduction to neutron stars and millisecond pulsars is given in Chapter 1. A selection of key papers published in 2011-2014 are presented in Chapters 2-10, ordered within five main research areas (ultra-stripped SNe in close binaries, massive NSs in close binaries, spin-up of MSPs, formation of a triple MSP, accretion-induced collapse of white dwarfs). In Chapter 11, I give a summary and highlight ongoing projects and further outlook. Slight differences (typos or minor language editing) may in some cases appear in comparison with the actual journal versions of the published papers. Hence, for citing these papers please refer to the journals., Comment: 267 pages, 99 figures, Habilitation thesis, University of Bonn

Since their discovery in the late 1960's the population of known neutron stars (NSs) has grown to ~2500. The last five decades of observations have yielded many surprises and demonstrated that the observational properties of NSs are remarkably diverse. The surveys that will be performed with SKA (the Square Kilometre Array) will produce a further tenfold increase in the number of Galactic NSs known. Moreover, the SKA's broad spectral coverage, sub-arraying and multi-beaming capabilities will allow us to characterise these sources with unprecedented efficiency, in turn enabling a giant leap in the understanding of their properties. Here we review the NS population and outline our strategies for studying each of the growing number of diverse classes that are populating the "NS zoo". Some of the main scientific questions that will be addressed by the much larger statistical samples and vastly improved timing efficiency provided by SKA include: (i) the spin period and spin-down rate distributions (and thus magnetic fields) at birth, and the associated information about the SNe wherein they are formed; (ii) the radio pulsar-magnetar connection; (iii) the link between normal radio pulsars, intermittent pulsars and rotating radio transients; (iv) the slowest possible spin period for a radio pulsar (revealing the conditions at the pulsar death-line); (v) proper motions of pulsars (revealing SN kick physics); (vi) the mass distribution of NSs (vii) the fastest possible spin period for a recycled pulsar (constraining magnetosphere-accretion disc interactions, gravitational wave radiation and the equation-of-state); (viii) the origin of high eccentricity millisecond pulsars (MSPs); (ix) the formation channels for recently identified triple systems; and finally (x) how isolated MSPs are formed. We expect that the SKA will break new ground unveiling exotic systems that will challenge... [abridged], Comment: Comments: 28 pages, 4 figures, to be published in: "Advancing Astrophysics with the Square Kilometre Array", Proceedings of Science, PoS(AASKA14)039

The recent detection of hypervelocity stars (HVSs) as late-type B-stars and HVS candidate G/K-dwarfs raises the important question of their origin. In this Letter, we investigate the maximum possible velocities of such HVSs if they are produced from binaries which are disrupted via an asymmetric supernova explosion. We find that HVSs up to ~770 and ~1280 km/s are possible in the Galactic rest frame from this scenario for these two subclasses of HVSs, respectively. We conclude that whereas a binary origin cannot easily explain all of the observed velocities of B-type HVSs (in agreement with their proposed central massive black hole origin) it can indeed account for the far majority (if not all) of the recently detected G/K-dwarf HVS candidates., Comment: 6 pages, 6 figures, including appendix, in press, MNRAS Letters (Updated and a comment added on the spin axis of SN-induced HVSs)

Millisecond pulsars (MSPs) are believed to be old neutron stars, formed via Type Ib/c core-collapse supernovae, which have subsequently been spun up to high rotation rates via accretion from a companion star in a highly circularised low-mass X-ray binary. The recent discoveries of Galactic field binary MSPs in eccentric orbits, and mass functions compatible with that expected for helium white dwarf companions, PSR J2234+06 and PSR J1946+3417, therefore challenge this picture. Here we present a hypothesis for producing this new class of systems, where the MSPs are formed directly from a rotationally-delayed accretion-induced collapse of a super-Chandrasekhar mass white dwarf. We compute the orbital properties of the MSPs formed in such events and demonstrate that our hypothesis can reproduce the observed eccentricities, masses and orbital periods of the white dwarfs, as well as forecasting the pulsar masses and velocities. Finally, we compare this hypothesis to a triple star scenario., Comment: Published in MNRAS Letters (6 pages, 1 table and 3 figures), typo corrected in Fig. 1

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. 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. 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 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. We demonstrate that the final properties of these MSPs are, in general, remarkably similar to those of MSPs formed via the standard core-collapse SN channel. However, the resultant MSPs created via the AIC channel preferentially form in certain orbital period intervals. Finally, we discuss the link between AIC and young NSs in globular clusters. Our calculations are also applicable to progenitor binaries of SNe Ia under certain conditions. [Abridged], Comment: 26 pages, 20 figures, 2 tables. A few references added. A&A in press

Eclipsing binary millisecond pulsars (the so-called black widows and redbacks) can provide important information about accretion history, pulsar irradiation of their companion stars and the evolutionary link between accreting X-ray pulsars and isolated millisecond pulsars. However, the formation of such systems is not well understood, nor the difference in progenitor evolution between the two populations of black widows and redbacks. Whereas both populations have orbital periods between $0.1-1.0\;{\rm days}$ their companion masses differ by an order of magnitude. In this paper, we investigate the formation of these systems via evolution of converging low-mass X-ray binaries by employing the MESA stellar evolution code. Our results confirm that one can explain the formation of most of these eclipsing binary millisecond pulsars using this scenario. More notably, we find that the determining factor for producing either black widows or redbacks is the efficiency of the irradiation process, such that the redbacks absorb a larger fraction of the emitted spin-down energy of the radio pulsar (resulting in more efficient mass loss via evaporation) compared to that of the black widow systems. We argue that geometric effects (beaming) is responsible for the strong bimodality of these two populations. Finally, we conclude that redback systems do not evolve into black widow systems with time., Comment: 19 pages, 5 figures. Accepted for publication

Many physically motivated extensions to general relativity (GR) predict significant deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 +/- 0.04 solar mass pulsar in a 2.46-hr orbit with a 0.172 +/- 0.003 solar mass white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling., Comment: This is the authors' version of the work. The definite version is published in Science Online, 26 April 2013, Vol: 340, Issue: 6131 doi: 10.1126/science.1233232. 54 pages, 18 figures, 4 tables

Millisecond pulsars (MSPs) are mainly characterised by their spin periods, B-fields and masses - quantities which are largely affected by previous interactions with a companion star in a binary system. In this paper, we investigate the formation mechanism of MSPs 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 MSPs 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 PPdot-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 MSPs 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 0.002-0.007 M_sun which is just about sufficient to explain their observed spin periods. We introduce isochrones of radio MSPs in the PPdot-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 the massive pulsar J1614-2230 (Paper I) and put new constraints on the birth masses of a number of recycled pulsars. [Abridged], Comment: MNRAS in press, 32 pages, 14 figures, 4 tables, appendix. Version 2: minor typos corrected

Millisecond pulsars are old neutron stars that have been spun up to high rotational frequencies via accretion of mass from a binary companion star. An important issue for understanding the physics of the early spin evolution of millisecond pulsars is the impact of the expanding magnetosphere during the terminal stages of the mass-transfer process. Here I report binary stellar evolution calculations that show that the braking torque acting on a neutron star, when the companion star decouples from its Roche-lobe, is able to dissipate >50% of the rotational energy of the pulsar. This effect may explain the apparent difference in observed spin distributions between x-ray and radio millisecond pulsars and help account for the noticeable age discrepancy with their young white dwarf companions., Comment: Published in Science, 3 February, 2012 (30 pages, 4 figures, reference numbers corrected in SOM, Fig.2 rotated)

In this review I discuss the characteristics and the formation of all classes of millisecond pulsars (MSPs). The main focus is on the stellar astrophysics of X-ray binaries leading to the production of fully recycled MSPs with white dwarf (WD) or substellar semi-degenerate companions. Depending on the nature of the companion star MSPs are believed to form from either low-mass X-ray binaries (LMXBs) or intermediate-mass X-ray binaries (IMXBs). For each of these two classes of X-ray binaries the evolutionary status of the donor star -- or equivalently, the orbital period -- at the onset of the Roche-lobe overflow (RLO) is the determining factor for the outcome of the mass-transfer phase and thus the nature of the MSP formed. Furthermore, the formation of binary MSPs is discussed in context of the (P,P_dot)-diagram, as well as new interpretations of the Corbet diagram. Finally, I present new models of Case A RLO of IMXBs in order to reproduce the two solar mass pulsar PSR J1614-2230., Comment: 17 pages, 9 figures (color). Invited review at ESO workshop "Evolution of Compact Binaries", March 2011, Vina del Mar (Chile), ASP Conf. Series. Fig.7 is available in high resolution on request