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.

We present multi-wavelength observations of two gap transients followed by the Carnegie Supernova Project-II and supplemented with data obtained by a number of different programs. Here in the first of two papers, we focus on the intermediate luminosity red transient (ILRT) designated SNhunt120, while in a companion paper we examine the luminous red novae AT 2014ej. Our data set for SNhunt120 consists of an early optical discovery, estimated to be within 3 days after outburst, the subsequent optical and near-infrared broadband followup extending over a $\sim$2 month period, two visual- and two near-infrared wavelength spectra, and Spitzer Space Telescope observations extending from early ($+$28 d) to late ($+$1155 d) phases. SNhunt120 resembles other ILRTs such as NGC 300-2008-OT and SN 2008S, and like these other ILRTs, SNhunt120 exhibits prevalent mid-infrared emission at both early and late phases. From the comparison of SNhunt120 and other ILRTs to electron-capture supernova simulations, we find that the current models underestimate the explosion kinetic energy and thereby produce synthetic light curves that over-estimate the luminosity. Finally, examination of pre-outburst Hubble Space Telescope images yields no progenitor detection.

We present optical and near-infrared broadband photometry and optical spectra of AT 2014ej from the the Carnegie Supernova Project-II. These observations are complemented with data from the CHilean Automatic Supernova sEarch, the Public ESO Spectroscopic Survey of Transient Objects, and from the Backyard Observatory Supernova Search. Observational signatures of AT 2014ej reveal that it is similar to other members of the gap-transient subclass known as luminous red novae (LRNe), including the ubiquitous double hump light curve and spectral properties akin to the LRN SN 2017jfs. A medium-dispersion, visual-wavelength spectrum of AT 2014ej taken the Magellan Clay telescope, exhibits a P Cygni H$\alpha$ feature characterized by a blue velocity at zero intensity of $\approx 110$ km s$^{-1}$ and a P Cygni minimum velocity of $\approx70$ km s$^{-1}$, and which we attribute to emission from a circumstellar wind. Inspection of pre-outbust Hubble Space Telescope images yields no conclusive progenitor detection. In comparison with a sample of LRNe from the literature, AT 2014ej lies among the brighter end of the luminosity distribution. Comparison of the ultra-violet, optical, infrared (UVOIR) light curves of well-observed LRNe to common-envelope evolution models from the literature, indicates the models under predict the luminosity of the comparison sample at all phases and also produce inconsistent time-scales of the secondary peak. Future efforts to model LRNe should expand upon the current parameter space explored and therefore may consider more massive systems and a wider range of dynamical timescales.

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.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

Ultra-luminous X-ray sources (ULXs) are X-ray sources located off-galactic centre and have luminosities exceeding the Eddington limit of a stellar-mass black hole ($L_X>10^{39}\;{\rm erg\,s}^{-1}$). With the discovery of X-ray pulsations in some of these objects (e.g. M82 X-2), we know a part of the ULX population has a neutron star (NS) accretor. We present systematic modelling of intermediate-mass X-ray binaries (IMXBs; donor-star mass range $2.0$ to $8.0\;$M$_{\odot}$) to explain the formation of this sub-population of ULXs. Using MESA, we explore the allowed initial parameter space of binary systems consisting of a neutron star and an intermediate-mass donor star, that could explain the observed properties of ULXs. 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. 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 NS mass ($1.3\;$M$_{\odot}$ and $2.0\;$M$_{\odot}$) we find that IMXBs can produce NS-ULXs with typical time-averaged isotropic-equivalent X-ray luminosities between $10^{39}$--$10^{41}\;{\rm erg\,s}^{-1}$ on a timescale of up to $\sim\!1.0\;$Myr for the lower luminosities, assuming geometrical beaming. Finally, we estimate their likelihood of detection, the types of white-dwarf remnants left behind by the donors, as well as the total amount of mass accreted by the neutron stars. 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 IMXBs in a super-Eddington mass-transfer phase.

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.

Radio pulsars in short-period eccentric binary orbits can be used to study both gravitational dynamics and binary evolution. The binary system containing PSR J1141$-$6545 includes a massive white dwarf (WD) companion that formed before the gravitationally bound young radio pulsar. We observe a temporal evolution of the orbital inclination of this pulsar that we infer is caused by a combination of a Newtonian quadrupole moment and Lense-Thirring precession of the orbit resulting from rapid rotation of the WD. Lense-Thirring precession, an effect of relativistic frame-dragging, is a prediction of general relativity. This detection is consistent with the evolutionary scenario in which the WD accreted matter from the pulsar progenitor, spinning up the WD to a period $< 200$ seconds.

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.

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.

Next-generation observations will revolutionize our understanding of binary black holes and will detect new sources, such as intermediate-mass black holes. Primary science goals include: Discover binary black holes throughout the observable Universe; Reveal the fundamental properties of black holes; Uncover the seeds of supermassive black holes.

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.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

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.

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.

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.

The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. By performing a series of three-dimensional hydrodynamical simulations, we investigate how CCSN explosions affect their binary companion. We find that the amount of removed stellar mass, the resulting impact velocity, and the chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Also, it is foud that the impact effects of CCSN ejecta on the structure of main-sequence (MS) companions, and thus their long term post-explosion evolution, is in general not be dramatic.

We present a suite of seven 3D supernova simulations of non-rotating low-mass progenitors using multi-group neutrino transport. Our simulations cover single star progenitors with zero-age main sequence masses between $9.6 M_\odot$ and $12.5 M_\odot$ and (ultra)stripped-envelope progenitors with initial helium core masses between $2.8 M_\odot$ and $3.5 M_\odot$. We find explosion energies between $0.1\,\mathrm{Bethe}$ and $0.4\,\mathrm{Bethe}$, which are still rising by the end of the simulations. Although less energetic than typical events, our models are compatible with observations of less energetic explosions of low-mass progenitors. In six of our models, the mass outflow rate already exceeds the accretion rate onto the proto-neutron star, and the mass and angular momentum of the compact remnant have closely approached their final value, barring the possibility of later fallback. While the proto-neutron star is still accelerated by the gravitational tug of the asymmetric ejecta, the acceleration can be extrapolated to obtain estimates for the final kick velocity. We obtain gravitational neutron star masses between $1.22 M_\odot$ and $1.44 M_\odot$, kick velocities between $11\, \mathrm{km}\, \mathrm{s}^{-1}$ and $695\, \mathrm{km}\, \mathrm{s}^{-1}$, and spin periods from $20\, \mathrm{ms}$ to $2.7\,\mathrm{s}$, which suggests that typical neutron star birth properties can be naturally obtained in the neutrino-driven paradigm. We find a loose correlation between the explosion energy and the kick velocity. There is no indication of spin-kick alignment, but a correlation between the kick velocity and the neutron star angular momentum, which needs to be investigated further as a potential point of tension between models and observations.

We determine the properties of the binary star V106 in the old open cluster NGC6791. We identify the system to be a blue straggler cluster member by using a combination of ground-based and Kepler photometry and multi-epoch spectroscopy. The properties of the primary component are found to be $M_{\rm p}\sim1.67 \rm M_{\odot}$, more massive than the cluster turn-off, with $R_{\rm p}\sim1.91 \rm R_{\odot}$ and $T_{\rm eff}=7110\pm100$ K. The secondary component is highly oversized and overluminous for its low mass with $M_{\rm s}\sim0.182 \rm M_{\odot}$, $R_{\rm s}\sim0.864 \rm R_{\odot}$ and $T_{\rm eff}=6875\pm200$ K. We identify this secondary star as a bloated (proto) extremely low-mass helium white dwarf. These properties of V106 suggest that it represents a typical Algol-paradox system and that it evolved through a mass-transfer phase which provides insight into its past evolution. We present a detailed binary stellar evolution model for the formation of V106 using the MESA code and find that the mass-transfer phase only ceased about 40 Myr ago. Due to the short orbital period (P=1.4463 d) another mass-transfer phase is unavoidable once the current primary star evolves towards the red giant phase. We argue that V106 will evolve through a common-envelope phase within the next 100 Myr and merge to become a single over-massive giant. The high mass will make it appear young for its true age, which is revealed by the cluster properties. Therefore, V106 is potentially a prototype progenitor of old field giants masquerading as young.

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

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.

The recent discoveries of many double neutron star systems and their detection as LIGO-Virgo merger events call for a detailed understanding of their origin. Explosions of ultra-stripped stars in binary systems have been shown to play a key role in this context and have also generated interest as a potential explanation for rapidly evolving hydrogen-free transients. Here we present the first attempt to model such explosions based on binary evolution calculations that follow the mass transfer to the companion to obtain a consistent core-envelope structure as needed for reliable predictions of the supernova transient. We simulate the explosion in 2D and 3D, and confirm the modest explosion energies ~10^50erg and small kick velocities reported earlier in 2D models based on bare carbon-oxygen cores. The spin-up of the neutron star by asymmetric accretion is small in 3D with no indication of spin-kick alignment. Simulations up to shock breakout show the mixing of sizeable amounts of iron group material into the helium envelope. In view of recent ideas for a mixing-length treatment (MLT) of Rayleigh-Taylor instabilities in supernovae, we perform a detailed analysis of the mixing, which reveals evidence for buoyancy-drag balance, but otherwise does not support the MLT approximation. The mixing may have implications for the spectroscopic signatures of ultra-stripped supernovae that need to be investigated in the future. Our stellar evolution calculation also predicts presupernova mass loss due to an off-centre silicon deflagration flash, which suggests that supernovae from extremely stripped cores may show signs of interactions with circumstellar material.

Radio pulsars found in binary systems with short orbital periods are usually fast spinning as a consequence of recycling via mass transfer from their companion stars; this process is also thought to decrease the magnetic field of the neutron star being recycled. Here, we report on timing observations of the recently discovered binary PSR J1755$-$2550 and find that this pulsar is an exception: with a characteristic age of 2.1 Myr, it is relatively young; furthermore, with a spin period of 315 ms and a surface magnetic field strength at its poles of 0.88$\times$10$^{12}$ G the pulsar shows no sign of having been recycled. Based on its timing and orbital characteristics, the pulsar either has a massive white dwarf (WD) or a neutron star (NS) companion. To distinguish between these two cases, we searched radio observations for a potential recycled pulsar companion and analysed archival optical data for a potential WD companion. Neither work returned conclusive detections. We apply population synthesis modelling and find that both solutions are roughly equally probable. Our population synthesis also predicts a minimum mass of 0.90 M$_{\odot}$ for the companion star to PSR J1755$-$2550 and we simulate the systemic runaway velocities for the resulting WDNS systems which may merge and possibly produce Ca-rich supernovae. Whether PSR J1755$-$2550 hosts a WD or a NS companion star, it is certainly a member of a rare subpopulation of binary radio pulsars.

Pulsars in relativistic binary systems have emerged as fantastic natural laboratories for testing theories of gravity, the most prominent example being the double pulsar, PSR J0737$-$3039. The HTRU-South Low Latitude pulsar survey represents one of the most sensitive blind pulsar surveys taken of the southern Galactic plane to date, and its primary aim has been the discovery of new relativistic binary pulsars. Here we present our binary pulsar searching strategy and report on the survey's flagship discovery, PSR J1757$-$1854. A 21.5-ms pulsar in a relativistic binary with an orbital period of 4.4 hours and an eccentricity of 0.61, this double neutron star (DNS) system is the most accelerated pulsar binary known, and probes a relativistic parameter space not yet explored by previous pulsar binaries.

To confirm the nature of the donor star in the ultra-compact X-ray binary candidate 47 Tuc X9, we obtained optical spectra (3,000$-$10,000 {\AA}) with the Hubble Space Telescope / Space Telescope Imaging Spectrograph. We find no strong emission or absorption features in the spectrum of X9. In particular, we place $3\sigma$ upper limits on the H$\alpha$ and HeII $\lambda 4686$ emission line equivalent widths $-$EW$_{\mathrm{H\alpha}} \lesssim 14$ {\AA} and $-$EW$_{\mathrm{HeII}} \lesssim 9$ {\AA}, respectively. This is much lower than seen for typical X-ray binaries at a similar X-ray luminosity (which, for $L_{\mathrm{2-10 keV}} \approx 10^{33}-10^{34}$ erg s$^{-1}$ is typically $-$EW$_{\mathrm{H\alpha}} \sim 50$ {\AA}). This supports our previous suggestion (by Bahramian et al.) of an H-poor donor in X9. We perform timing analysis on archival far-ultraviolet, $V$ and $I$-band data to search for periodicities. In the optical bands we recover the seven-day superorbital period initially discovered in X-rays, but we do not recover the orbital period. In the far-ultraviolet we find evidence for a 27.2 min period (shorter than the 28.2 min period seen in X-rays). We find that either a neutron star or black hole could explain the observed properties of X9. We also perform binary evolution calculations, showing that the formation of an initial black hole / He-star binary early in the life of a globular cluster could evolve into a present-day system such as X9 (should the compact object in this system indeed be a black hole) via mass-transfer driven by gravitational wave radiation.

PSR~J1618$-$3921 is an $11.99$-ms pulsar in a $22.7$-d orbit around a likely low-mass He white dwarf companion, discovered in a survey of the intermediate Galactic latitudes at 1400 MHz conducted with the Parkes radio telescope in the late 1990s. Although PSR~J1618$-$3921 was discovered more than 15 years ago, only limited information has been published about this pulsar which has a surprisingly large orbital eccentricity ($e \simeq 0.027$), considering its high spin frequency and likely small companion mass. The focus of this work is a precise measurement of the spin, astrometric and orbital characteristics of PSR J1618$-$3921. This was done with timing observations made at the Nan\c{c}ay Radio Telescope, from 2009 to 2017. We analyzed the timing data recorded at the Nan\c{c}ay Radio Telescope over several years to characterize the properties of PSR~J1618$-$3921. A rotation ephemeris for this pulsar was obtained by analyzing the times of arrival of the radio pulses at the telescope. We confirm the unusual eccentricity of PSR~J1618$-$3921 and discuss several hypotheses regarding its formation, in the context of other discoveries of recycled pulsars in eccentric orbits.

We report the discovery of PSR J1757$-$1854, a 21.5-ms pulsar in a highly-eccentric, 4.4-h orbit around 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 (GW) damping, with a merger time of 76 Myr. Following a 1.6-yr timing campaign, we have measured five post-Keplerian (PK) parameters, yielding the two component masses ($m_\text{p}=1.3384(9)\,\text{M}_\odot$ and $m_\text{c}=1.3946(9)\,\text{M}_\odot$) plus three tests of general relativity (GR), 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-$\sigma$ measurements of both the contribution of Lense-Thirring precession to the rate of change of the semi-major axis and the relativistic deformation of the orbit within $\sim7-9$ years, PSR J1757$-$1854 stands out as a unique laboratory for new tests of gravitational theories.

In this work, we report the discovery and characterization of PSR J1411+2551, a new binary pulsar discovered in the Arecibo 327 MHz Drift Pulsar Survey. Our timing observations of the radio pulsar in the system span a period of about 2.5 years. This timing campaign allowed a precise measurement of its spin period (62.4 ms) and its derivative (9.6 $\pm$ 0.7) $\times 10^{-20}\, \rm s\, s^{-1}$; from these, we derive a characteristic age of $\sim 10\,$Gyr and a surface magnetic field strength of 2.5 $\times 10^{9}$ G. These numbers indicate that this pulsar was mildly recycled by accretion of matter from the progenitor of the companion star. The system has an eccentric ($e\, = \, 0.17$) 2.61 day orbit. This eccentricity allows a highly significant measurement of the rate of advance of periastron, $\dot{\omega} = 0.07686 \pm 0.00046 ^{\circ}~{\rm yr}^{-1}$. Assuming general relativity accurately models the orbital motion, this implies a total system mass M = $2.538 \pm 0.022 M_{\odot}$. The minimum companion mass is $0.92\, M_{\odot}$ and the maximum pulsar mass is $1.62\, M_{\odot}$. The large companion mass and the orbital eccentricity suggest that PSR J1411+2551 is a double neutron star system; the lightest known to date including the DNS merger GW 170817. Furthermore, the relatively low orbital eccentricity and small proper motion limits suggest that the second supernova had a relatively small associated kick; this and the low system mass suggest that it was an ultra-stripped supernova.

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 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 $\sim 0.02\;M_{\odot}$. 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 which 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 $\gamma$-ray burst.

We present a high-precision timing analysis of PSR J1933$-$6211, a millisecond pulsar (MSP) with a 3.5-ms spin period and a white dwarf (WD) companion, using data from the Parkes radio telescope. Since we have accurately measured the polarization properties of this pulsar we have applied the matrix template matching approach in which the times of arrival are measured using full polarimetric information. We achieved a weighted root-mean-square timing residuals (rms) of the timing residuals of 1.23 $\rm \mu s$, 15.5$\%$ improvement compared to the total intensity timing analysis. After studying the scintillation properties of this pulsar we put constraints on the inclination angle of the system. Based on these measurements and on $\chi^2$ mapping we put a 2-$\sigma$ upper limit on the companion mass (0.44 M$_\odot$). Since this mass limit cannot reveal the nature of the companion we further investigate the possibility of the companion to be a He WD. Applying the orbital period-mass relation for such WDs, we conclude that the mass of a He WD companion would be about 0.26$\pm$0.01 M$_\odot$ which, combined with the measured mass function and orbital inclination limits, would lead to a light pulsar mass $\leqslant$ 1.0 M$_\odot$. This result seems unlikely based on current neutron star formation models and we therefore conclude that PSR J1933$-$6211 most likely has a CO WD companion, which allows for a solution with a more massive pulsar.

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.

In search for the counterpart to the Fermi-LAT source 3FGL J0838.8-2829, we performed a multi-wavelength campaign, in the X-ray band with Swift and XMM-Newton, performed infrared, optical (with OAGH, ESO-NTT and IAC80) and radio (ATCA) observations, as well as analysed archival hard X-ray data taken by INTEGRAL. We report on three X-ray sources consistent with the position of the Fermi-LAT source. We confirm the identification of the brightest object, RX J0838-2827, as a magnetic cataclysmic variable, that we recognize as an asynchronous system (not associated with the Fermi-LAT source). RX J0838-2827, is extremely variable in the X-ray and optical bands, and timing analysis reveals the presence of several periodicities modulating its X-ray and optical emission. The most evident modulations are interpreted as due to the binary system orbital period of ~1.64h and the white dwarf spin period of ~1.47h. A strong flux modulation at ~15h is observed at all energy bands, consistent with the beat frequency between spin and orbital periods. Optical spectra show prominent Hbeta, HeI and HeII emission lines Doppler-modulated at the orbital period and at the beat period. Therefore, RX J0838-2827, accretes through a diskless configuration and could be either a strongly asynchronous polar or a rare example of a pre-polar system in its way to reach synchronism. Among the other two X-ray sources, XMM J083850.4-282759 showed a variable X-ray emission, with a powerful flare lasting ~600s, similar to what is observed in transitional millisecond pulsars during the sub-luminous disc state: that would possibly associate this source with the Fermi-LAT source.

We report the discovery and the results of follow-up timing observations of PSR J2045+3633 and PSR J2053+4650, two binary pulsars found in the Northern High Time Resolution Universe pulsar survey being carried out with the Effelsberg radio telescope. Having spin periods of 31.7 ms and 12.6 ms respectively, and both with massive white dwarf companions, $M_{c}\, > \, 0.8\, M_{\odot}$, the pulsars can be classified as mildly recycled. PSR J2045+3633 is remarkable due to its orbital period (32.3 days) and eccentricity $e\, = \, 0.01721244(5)$ which is among the largest ever measured for this class. After almost two years of timing the large eccentricity has allowed the measurement of the rate of advance of periastron at the 5-$\sigma$ level, 0.0010(2)$^\circ~\rm yr^{-1}$. Combining this with a detection of the orthometric amplitude of the Shapiro delay, we obtained the following constraints on the component masses (within general relativity): $M_{p}\, = \, 1.33^{+0.30}_{-0.28}\, M_{\odot}$, and $M_{c}\, = \, 0.94^{+0.14}_{-0.13}\, M_{\odot}$. PSR J2053+4650 has a 2.45-day circular orbit inclined to the plane of the sky at an angle $i\, = \, 85.0^{+0.8}_{-0.9}\,{^\circ}$. In this nearly edge-on case the masses can be obtained from the Shapiro delay alone. Our timing observations resulted in a significant detection of this effect giving: $M_{p}\, = \, 1.40^{+0.21}_{-0.18}\, M_{\odot}$, and $M_{c}\, = \, 0.86^{+0.07}_{-0.06}\, M_{\odot}$.

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.

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.

When a massive star in a binary system explodes as a supernova its companion star may be polluted with heavy elements from the supernova ejecta. Such a pollution had been detected in a handful of post-supernova binaries (Gonzalez Hernandez et al. 2011), but none of them is associated with a supernova remnant. We report the discovery of a solar-type star in a close, eccentric binary system with a neutron star within the young Galactic supernova remnant RCW 86. Our discovery implies that the supernova progenitor was a moving star, which exploded near the edge of its wind bubble and lost most of its initial mass due to common-envelope evolution shortly before core collapse. We find that the solar-type star is strongly polluted with calcium and other elements, which places the explosion within the class of calcium-rich supernovae -- faint and fast transients (Filippenko et al. 2003; Kasliwal et al. 2012}, whose origin is strongly debated (Kawabata et al. 2010; Waldman et al. 2011), and provides the first observational evidence that supernovae of this type can arise from core-collapse explosions (Kawabata et al. 2010; Moriya et al. 2010).

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.

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. 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 that process. We calculate envelope binding energies of donor stars with initial masses between 4 and 115 Msun for metallicities of Z=Zsun/2 and Z=Zsun/50, and derive minimum masses of in-spiralling objects needed to eject these envelopes. We find that CE evolution, besides from producing WD-WD and NS-NS binaries, may, in principle, also produce massive BH-BH systems with individual BH component masses up to ~50-60 Msun, in particular for donor stars evolved to giants. 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. Massive stars in a wide range of metallicities and with initial masses up to at least 115 Msun may possibly shed their envelopes and survive CE evolution, depending on their initial orbital parameters, similarly to the situation for intermediate mass and low-mass stars with degenerate cores. We conclude that based on stellar structure calculations, and in the view of the usual simple energy budget analysis, events like GW150914 and GW151226 could possibly be produced from the CE channel. Calculations of post-CE orbital separations, however, and thus the estimated LIGO detection rates, remain highly uncertain. [Abridged]

We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting naturally for the presence of He, as well as traces of metals such as Mg and Ca, typically found in the atmospheres of ELM proto-WDs. Here we investigate whether rotational mixing can maintain a sufficient amount of He in the deeper driving region of the star, such that it can fuel, through HeII-HeIII ionization, the observed pulsations in this type of stars. Using state-of-the-art evolutionary models computed with MESA, we show that rotational mixing can indeed explain qualitatively the very existence and general properties of the known pulsating, mixed-atmosphere ELM proto-WDs. Moreover, such objects are very likely to pulsate again during their final WD cooling phase.

Timing results for the black-widow pulsar J2051-0827 are presented, using a 21-year dataset from four European Pulsar Timing Array telescopes and the Parkes radio telescope. This dataset, which is the longest published to date for a black-widow system, allows for an improved analysis that addresses previously unknown biases. While secular variations, as identified in previous analyses, are recovered, short-term variations are detected for the first time. Concurrently, a significant decrease of approx. 2.5x10-3 cm-3 pc in the dispersion measure associated with PSR J2051-0827 is measured for the first time and improvements are also made to estimates of the proper motion. Finally, PSR J2051-0827 is shown to have entered a relatively stable state suggesting the possibility of its eventual inclusion in pulsar timing arrays.

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.

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.

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.

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.

We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.

Recent advances in gravitational-wave astronomy make the direct detection of gravitational waves from the merger of two stellar-mass compact objects a realistic prospect. Evolutionary scenarios towards mergers of double compact objects generally invoke common-envelope evolution which is poorly understood, leading to large uncertainties in merger rates. We explore the alternative scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary which remain fully mixed due to their tidally induced high spin. We use the public stellar-evolution code MESA to systematically study this channel by means of detailed simulations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass ratios close to one, in two mass ranges, ~25...60msun and >~ 130msun, with pair instability supernovae (PISNe) being produced in-between. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for it. We map the initial parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, Kerr parameters and predict event rates. We typically find that for Z~<Z_sun/10, there is one BH+BH merger for ~1000 core-collapse supernovae. The advanced LIGO (aLIGO) detection rate is more uncertain and depends on the metallicity evolution. Deriving upper and lower limits from a local and a global approximation for the metallicity distribution of massive stars, we estimate aLIGO detection rates (at design limit) of ~19-550 yr^(-1) for BH+BH mergers below the PISN gap and of ~2.1-370 yr^(-1) above the PISN gap. Even with conservative assumptions, we find that aLIGO should soon detect BH+BH mergers from the MOB scenario and that these could be the dominant source for aLIGO detections.

The explosion of ultra-stripped stars in close binaries may explain new discoveries of weak and fast optical transients. We have demonstrated that helium star companions to neutron stars (NSs) may evolve into naked metal cores as low as ~1.5 Msun, barely above the Chandrasekhar mass limit, by the time they explode. Here we present a new systematic investigation of the progenitor evolution leading to such ultra-stripped supernovae (SNe), in some cases yielding pre-SN envelopes of less than 0.01 Msun. We discuss the nature of these SNe (electron-capture vs iron core-collapse) and their observational light-curve properties. Ultra-stripped SNe are highly relevant for binary pulsars, as well as gravitational wave detection of merging NSs by LIGO/VIRGO, since these events are expected to produce mainly low-kick NSs in the mass range 1.10-1.80 Msun.

The evolution of close-orbit progenitor binaries of double neutron star (DNS) systems leads to supernova (SN) explosions of ultra-stripped stars. The amount of SN ejecta mass is very limited from such, more or less, naked metal cores with envelope masses of only 0.01-0.2 Msun. The combination of little SN ejecta mass and the associated possibility of small NS kicks is quite important for the characteristics of the resulting DNS systems left behind. Here, we discuss theoretical predictions for DNS systems, based on Case BB Roche-lobe overflow prior to ultra-stripped SNe, and briefly compare with observations.

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.

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]

Gaia will identify several 1e5 white dwarfs, most of which will be in the solar neighborhood at distances of a few hundred parsecs. Ground-based optical follow-up spectroscopy of this sample of stellar remnants is essential to unlock the enormous scientific potential it holds for our understanding of stellar evolution, and the Galactic formation history of both stars and planets.

We present Hubble Space Telescope UV spectra of the 4.6 h period double white dwarf SDSS J125733.63+542850.5. Combined with Sloan Digital Sky Survey optical data, these reveal that the massive white dwarf (secondary) has an effective temperature T2 = 13030 +/- 70 +/- 150 K and a surface gravity log g2 = 8.73 +/- 0.05 +/- 0.05 (statistical and systematic uncertainties respectively), leading to a mass of M2 = 1.06 Msun. The temperature of the extremely low-mass white dwarf (primary) is substantially lower at T1 = 6400 +/- 37 +/- 50 K, while its surface gravity is poorly constrained by the data. The relative flux contribution of the two white dwarfs across the spectrum provides a radius ratio of R1/R2 = 4.2, which, together with evolutionary models, allows us to calculate the cooling ages. The secondary massive white dwarf has a cooling age of about 1 Gyr, while that of the primary low-mass white dwarf is likely to be much longer, possibly larger than 5 Gyrs, depending on its mass and the strength of chemical diffusion. These results unexpectedly suggest that the low-mass white dwarf formed long before the massive white dwarf, a puzzling discovery which poses a paradox for binary evolution.

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.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of very faint X-ray binaries, orbital period distribution of black hole X-ray binaries and neutron star spin up. For a summary, we refer to the paper.

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.

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]

A large number of low-mass (< 0.20 M_sun) helium white dwarfs (He WDs) have recently been discovered. The majority of these are orbiting another WD or a millisecond pulsar (MSP) in a close binary system; a few examples are found to show pulsations or to have a main-sequence star companion. There appear to be discrepancies between the current theoretical modelling of such low-mass He WDs and a number of key observed cases, indicating that their formation scenario remains to be fully understood. Here we investigate the formation of detached proto-He WDs in close-orbit low-mass X-ray binaries (LMXBs). Our prime focus is to examine the thermal evolution and the contraction phase towards the WD cooling track and investigate how this evolution depends on the WD mass. Our calculations are then compared to the most recent observational data. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in a large number of close-orbit LMXBs. Subsequently, we followed the evolution of the detached low-mass proto-He WDs, including stages with residual shell hydrogen burning and vigorous flashes caused by unstable CNO burning. We find that the time between Roche-lobe detachment until the low-mass proto-He WD reaches the WD cooling track is typically Delta_t_proto = 0.5 - 2 Gyr, depending systematically on the WD mass and therefore on its luminosity. The lowest WD mass for developing shell flashes is ~0.21 M_sun for progenitor stars of mass M2 <= 1.5 M_sun (and ~0.18 M_sun for M2 = 1.6 M_sun). The long timescale of low-mass proto-He WD evolution can explain a number of recent observations, including some MSP systems hosting He WD companions with very low surface gravities and high effective temperatures. We find no evidence for Delta_t_proto to depend on the occurrence of flashes and thus question the suggested dichotomy in thermal evolution of proto-WDs.

Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs) which have been spun up to high rotation rates via accretion of matter from a companion star in a low-mass X-ray binary (LMXB). However, many details of this recycling scenario remain to be understood. 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 P_orb = 2-9 hr. In particular, we examine: i) if such observed systems can be reproduced from theoretical modelling using standard prescriptions of orbital angular momentum losses (i.e. with respect to the nature and the strength of magnetic braking), ii) if our computations of the Roche-lobe detachments can match the observed orbital periods, and iii) if the correlation between WD mass and orbital period (M_WD, P_orb) is valid for systems with P_orb < 2 days. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in ~ 400 close LMXB systems with different initial values of donor star mass, NS mass, orbital period and the so-called gamma-index of magnetic braking. Subsequently, we followed the orbital and the interior evolution of the detached low-mass (proto) He WDs, including stages with residual shell hydrogen burning. We find that a severe fine-tuning is necessary to reproduce the observed MSPs in tight binaries with He WD companions of mass < 0.20 M_sun, which suggests that something needs to be modified or is missing in the standard input physics of LMXB modelling. We demonstrate that the theoretically calculated (M_WD, P_orb)-relation is in general also valid for systems with P_orb < 2 days, although with a large scatter in He WD masses between 0.15-0.20 M_sun. The results of the thermal evolution of the (proto) He WDs are reported in a follow-up paper (Paper II).

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

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.

Rotation-powered millisecond radio pulsars have been spun up to their present spin period by a $10^8$ - $10^9$ yr long X-ray-bright phase of accretion of matter and angular momentum in a low-to-intermediate mass binary system. Recently, the discovery of transitional pulsars that alternate cyclically between accretion and rotation-powered states on time scales of a few years or shorter, has demonstrated this evolutionary scenario. Here, we present a thorough statistical analysis of the spin distributions of the various classes of millisecond pulsars to assess the evolution of their spin period between the different stages. Accreting sources that showed oscillations exclusively during thermonuclear type I X-ray bursts (nuclear-powered millisecond pulsars) are found to be significantly faster than rotation-powered sources, while accreting sources that possess a magnetosphere and show coherent pulsations (accreting millisecond pulsars) are not. On the other hand, if accreting millisecond pulsars and eclipsing rotation-powered millisecond pulsars form a common class of transitional pulsars, these are shown to have a spin distribution intermediate between the faster nuclear-powered millisecond pulsars and the slower non-eclipsing rotation-powered millisecond pulsars. We interpret these findings in terms of a spin-down due to the decreasing mass-accretion rate during the latest stages of the accretion phase, and in terms of the different orbital evolutionary channels mapped by the various classes of pulsars. We summarize possible instrumental selection effects, showing that even if an unbiased sample of pulsars is still lacking, their influence on the results of the presented analysis is reduced by recent improvements in instrumentation and searching techniques.

Millisecond pulsars (MSPs) and normal non-recycled pulsars are both detected in $\gamma$-rays. However, it appears that a much larger fraction of known energetic and nearby MSPs are detected in $\gamma$-rays, in comparison with normal pulsars, thereby making undetected $\gamma$-ray MSPs exceptions. In this paper, we demonstrate that the viewing angles (i.e. between the pulsar spin axis and the line of sight) are well described by the orbital inclination angles which, for binary MSPs with helium white dwarf companions, can be determined using the relationship between the orbital period and the white dwarf mass. We use the predicted viewing angles, in complement with values obtained from other constraints when available, to identify the causes of non-detection of energetic and nearby MSPs from the point of view of beaming geometry and orientation. We find evidence for slightly different viewing angle distributions, and postulate that energetic and nearby MSPs are mainly undetected in $\gamma$-rays simply because they are seen under unfavourable (i.e. small) viewing angles. We finally discuss the magnetic fields of $\gamma$-ray detected pulsars and show that pulsars which are efficient at converting their rotational energy into $\gamma$-ray emission may have overestimated dipolar magnetic field strengths.

We report on the discovery of four millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) pulsar survey being conducted at the Parkes 64-m radio telescope. All four MSPs are in binary systems and are likely to have white dwarf companions. In addition, we present updated timing solutions for 12 previously published HTRU MSPs, revealing new observational parameters such as five proper motion measurements and significant temporal dispersion measure variations in PSR J1017-7156. We discuss the case of PSR J1801-3210, which shows no significant period derivative after four years of timing data. Our best-fit solution shows a period derivative of the order of $10^{-23}$, an extremely small number compared to that of a typical MSP. However, it is likely that the pulsar lies beyond the Galactic Centre, and an unremarkable intrinsic period derivative is reduced to close to zero by the Galactic potential acceleration. Furthermore, we highlight the potential to employ PSR J1801-3210 in the strong equivalence principle test due to its wide and circular orbit. In a broader comparison with the known MSP population, we suggest a correlation between higher mass functions and the presence of eclipses in `very low-mass binary pulsars', implying that eclipses are observed in systems with high orbital inclinations. We also suggest that the distribution of the total mass of binary systems is inversely-related to the Galactic height distribution. Finally, we report on the first detection of PSRs J1543-5149 and J1811-2404 as gamma-ray pulsars.

The millisecond pulsar in a triple system (PSR J0337+1715, recently discovered by Ransom et al.) is an unusual neutron star with two orbiting white dwarfs. The existence of such a system in the Galactic field poses new challenges to stellar astrophysics for understanding evolution, interactions and mass-transfer in close multiple stellar systems. In addition, this system provides the first precise confirmation for a very wide-orbit system of the white dwarf mass-orbital period relation. Here we present a self-consistent, semi-analytical solution to the formation of PSR J0337+1715. Our model constrains the peculiar velocity of the system to be less than 160 km/s and brings novel insight to, for example, common envelope evolution in a triple system, for which we find evidence for in-spiral of both outer stars. Finally, we briefly discuss our scenario in relation to alternative models.

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_sun of material and the formation of a low-mass neutron star. 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 0.001-0.01. Finally, we argue that the second explosion in some double neutron star systems (for example, the double pulsar PSR J0737-3039B) was likely associated with an ultra-stripped SN Ic.

We report on timing observations of the recently discovered binary pulsar PSR J1952+2630 using the Arecibo Observatory. The mildly recycled 20.7-ms pulsar is in a 9.4-hr orbit with a massive, M_WD > 0.93 M_sun, white dwarf (WD) companion. We present, for the first time, a phase-coherent timing solution, with precise spin, astrometric, and Keplerian orbital parameters. This shows that the characteristic age of PSR J1952+2630 is 77 Myr, younger by one order of magnitude than any other recycled pulsar-massive WD system. We derive an upper limit on the true age of the system of 50 Myr. We investigate the formation of PSR J1952+2630 using detailed modelling of the mass-transfer process from a naked helium star on to the neutron star following a common-envelope phase (Case BB Roche-lobe overflow). From our modelling of the progenitor system, we constrain the accretion efficiency of the neutron star, which suggests a value between 100 and 300% of the Eddington accretion limit. We present numerical models of the chemical structure of a possible oxygen-neon-magnesium WD companion. Furthermore, we calculate the past and the future spin evolution of PSR J1952+2630, until the system merges in about 3.4 Gyr due to gravitational wave emission. Although we detect no relativistic effects in our timing analysis we show that several such effects will become measurable with continued observations over the next 10 years; thus PSR J1952+2630 has potential as a testbed for gravitational theories.

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.

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]

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.

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagreement, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

For many years it has been recognized that the terminal stages of mass transfer in a low-mass X-ray binary (LMXB) should cause the magnetosphere of the accreting neutron star to expand, leading to a braking torque acting on the spinning pulsar. After the discovery of radio millisecond pulsars (MSPs) it was therefore somewhat a paradox (e.g. Ruderman et al. 1989) how these pulsars could retain their fast spins following the Roche-lobe decoupling phase, RLDP. Here I present a solution to this so-called "turn-off problem" which was recently found by combining binary stellar evolution models with torque computations (Tauris 2012). The solution is that during the RLDP the spin equilibrium of the pulsar is broken and therefore it remains a fast spinning object. I briefly discuss these findings in view of the two observed spin distributions in the populations of accreting X-ray millisecond pulsars (AXMSPs) and radio MSPs.

Although the first millisecond pulsars (MSPs) were discovered 30 years ago we still do not understand all details of their formation process. Here, we present new results from Tauris, Langer & Kramer (2012) on the recycling scenario leading to radio MSPs with helium or carbon-oxygen white dwarf companions via evolution of low- and intermediate mass X-ray binaries (LMXBs, IMXBs). We discuss the location of the spin-up line in the (P,Pdot)-diagram and estimate the amount of accreted mass needed to obtain a given spin period and compare with observations. Finally, we constrain the true ages of observed recycled pulsars via calculated isochrones in the (P,Pdot)-diagram.

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]

PSR J1738+0333 is one of the four millisecond pulsars known to be orbited by a white dwarf companion bright enough for optical spectroscopy. Of these, it has the shortest orbital period, making it especially interesting for a range of astrophysical and gravity related questions. We present a spectroscopic and photometric study of the white dwarf companion and infer its radial velocity curve, effective temperature, surface gravity and luminosity. We find that the white dwarf has properties consistent with those of low-mass white dwarfs with thick hydrogen envelopes, and use the corresponding mass-radius relation to infer its mass; M_WD = 0.181 +/- +0.007/-0.005 solar masses. Combined with the mass ratio q=8.1 +/- 0.2 inferred from the radial velocities and the precise pulsar timing ephemeris, the neutron star mass is constrained to M_PSR = 1.47 +/- +0.07/-0.06 solar masses. Contrary to expectations, the latter is only slightly above the Chandrasekhar limit. We find that, even if the birth mass of the neutron star was only 1.20 solar masses, more than 60% of the matter that left the surface of the white dwarf progenitor escaped the system. The accurate determination of the component masses transforms this system in a laboratory for fundamental physics by constraining the orbital decay predicted by general relativity. Currently, the agreement is within 1 sigma of the observed decay. Further radio timing observations will allow precise tests of white dwarf models, assuming the validity of general relativity.

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.

We report near-infrared radial velocity measurements of the recently identified donor star in the high mass X-ray binary system OAO 1657-415 obtained in the H band using ISAAC on the VLT. Cross-correlation methods were employed to construct a radial velocity curve with a semi-amplitude of $22.1 \pm 3.5$ km s$^{-1}$. Combined with other measured parameters of this system this provides a dynamically determined neutron star mass of $1.42 \pm 0.26 M$_{\odot}$ and a mass of $14.3 \pm 0.8 M$_{\odot}$ for the Ofpe/WN9 highly evolved donor star. OAO 1657-415 is an eclipsing High Mass X-ray binary pulsar with the largest eccentricity and orbital period of any within its class. Of the ten known eclipsing X-ray binary pulsars OAO 1657-415 becomes the ninth with a dynamically determined neutron star mass solution and only the second in an eccentric system. Furthermore, the donor star in OAO 1657-415 is much more highly evolved than the majority of the supergiant donors in other High Mass X-ray binaries (HMXBs), joining a small but growing list of HMXBs donors with extensive hydrogen depleted atmospheres. Considering the evolutionary development of OAO 1657-415, we have estimated the binding energy of the envelope of the mass donor and find that there is insufficient energy for the removal of the donor's envelope via spiral-in, ruling out a Common Envelope evolutionary scenario. With its non-zero eccentricity and relatively large orbital period the identification of a definitive evolutionary pathway for OAO1657-415 remains problematic, we conclude by proposing two scenarios which may account for OAO1657-415 current orbital configuration.

The recent discovery of a 2 M_sun binary millisecond pulsar (Demorest et al. 2010) has not only important consequences for the equation-of-state of nuclear matter at high densities but also raises the interesting question if the neutron star PSR J1614-2230 was born massive. The answer is vital for understanding neutron star formation in core collapse supernovae. Furthermore, this system raises interesting issues about the nature of the progenitor binary and how it evolved during its mass exchanging X-ray phase. In this paper we discuss the progenitor evolution of PSR J1614-2230. We have performed detailed stellar evolution modelling of intermediate-mass X-ray binaries undergoing Case A Roche-lobe overflow (RLO) and applied an analytic parameterization for calculating the outcome of either a common envelope evolution or the highly super-Eddington isotropic re-emission mode. We find two viable possibilities for the formation of the PSR J1614-2230 system: either it contained a 2.2-2.6 M_sun giant donor star and evolved through a common envelope and spiral-in phase or, more likely, it descended from a close binary system with a 4.0-5.0 M_sun main sequence donor star via Case A RLO. We conclude that the neutron star must have been born with a mass of ~1.95 M_sun or 1.7+-0.15 M_sun, respectively - which significantly exceeds neutron star birth masses in previously discovered radio pulsar systems. Based on the expected neutron star birth masses from considerations of stellar evolution and explosion models, we find that the progenitor star of PSR J1614-2230 is likely to have been more massive than 20 M_sun.

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.

We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, {\mu}_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.

In this chapter we review the formation and evolution of compact binaries with neutron star and/or black hole components (i.e. LMXBs, HMXBs, binary pulsars). After an introduction we discuss stellar evolution in binaries and the processes involved in the mass-transfer phases of close binaries (RLO, CE, ang.mom.loss) with radiative/convective hydrogen or helium donor stars. We also describe the effects of accretion, asymmetric SN and systems merging as a result of gravitational wave radiation.

We have performed detailed population synthesis on a large number (20 million) of binary systems in order to investigate the properties of massive double degenerate binaries. We have included new important results in our input physics in order to obtain more reliable estimates of the merging timescales and relative formation rates. These improvements include refined treatment of the binding energy in a common envelope, helium star evolution and reduced kicks imparted to newborn black holes. The discovery and observations of GRB afterglows and the identification of host galaxies have allowed comparisons of theoretical distributions of merger sites with the observed distribution of afterglow positions relative to host galaxies. To help investigate the physical nature of short- and long-duration GRBs, we compute the distances of merging neutron stars (NS) and/or black holes (BH) from the centers of their host galaxies, as predicted by their formation scenario combined with motion in galactic potentials. Furthermore, we estimate the formation rate and merging rate of these massive double degenerate binaries. The latter is very important for the prospects of detecting gravitational waves with LIGO/VIRGO. We find that the expected detection rate for LIGO II is ~850 per year for galactic field sources and that this rate is completely dominated by merging BHBH binaries. Even LIGO I may detect such an event (~0.25 per year). Our preferred model estimate the Galactic field NSNS merger rate to be ~1.5*10^(-6) per year. For BHBH systems this model predicts a merger rate of ~9.7*10^{-6} per year. Our studies also reveal an accumulating numerous population of very wide orbit BHBH systems which never merge (t >> t_Hubble).

We investigate the evolution of pulsars in the (P,P-dot) diagram. We first present analytical formulae to follow the evolution of a pulsar using simple exponential models for magnetic field decay and alignment. We then compare these evolutionary tracks with detailed model calculations using ohmic decay of crustal neutron star magnetic fields. We find that, after an initial phase with a small braking index, n, pulsars evolve with enhanced torque decay (n >> 3) for about 1 Myr. The long term evolution depends on the impurity parameter of the crust. If impurities are negligible in older isolated pulsars we expect their true age to be approximately equal to their observed characteristic age, tau = P/(2*P-dot). It is not possible from data to constrain model parameters of the neutron star crust.

According to the standard picture for binary interactions, the outcome of binaries surviving the evolution through a common envelope (CE) and spiral-in phase is determined by the internal structure of the donor star at the onset of the mass transfer, as well as the poorly-known efficiency parameter, eta_CE}, for the ejection of the H-envelope of the donor. In this Research Note we discuss the bifurcation point which separates the ejected, unprocessed H-rich material from the inner core region of the donor (the central part of the star which will later contract to form a compact object). We demonstrate that the exact location of this point is very important for evaluating the binding energy parameter, lambda, which is used to determine the post-CE orbital separation. Here we compare various methods to define the bifurcation point (core/envelope boundary) of evolved stars with masses 4, 7, 10 and 20 M_sun. We consider the specific nuclear energy production rate profile, the change in the mass-density gradient (Bisscheroux 1998), the inner region containing less than 10% hydrogen, the method suggested by Han et al. (1994) and the entropy profile. We also calculated effective polytropic index profiles. The entropy profile method measures the convective boundary (at the onset of flatness in the specific entropy) which is not equivalent to the core boundary for RGB stars. Hence, this method is not applicable for RGB stars, unless the actual bifurcation point of a CE is located at the bottom of the outer convection zone (resulting in larger values of lambda and larger post-CE orbital separations). On the AGB, where highly degenerate and condensed cores are formed, we find good agreement between the various methods, except for massive (20 M_sun) stars.

We investigate the effect of orbital decay caused by nuclear expansion of a (sub)giant star in synchronous binary system. We compare this effect with the presence of a magnetic stellar wind and show that the additional transfer of orbital angular momentum into spin angular momentum is relatively important -- especially since it has been shown that the effect of magnetic braking saturates at short orbital periods.

We have investigated the structure of evolved giant stars with masses 3-10 M_sun in order to evaluate the binding energy of the envelope to the core prior to mass transfer in close binary systems. This binding energy is expressed by a parameter lambda which is crucial for determining the outcome of binaries evolving through a common envelope (CE) and spiral-in phase. We discuss the lambda-parameter and the efficiency of envelope ejection in the CE-phase, and show that lambda depends strongly on the evolutionary stage (i.e. stellar radius) of the donor star at the onset of the mass transfer. The existence of this relation enables us to introduce a new approach for solving the energy equation. For a given observed binary system we can derive a unique solution for the original mass and age of the donor star, as well as the pre-CE orbital period. We find that the value of lambda is typically between 0.2 and 0.8. But in some cases, particularly on the asymptotic giant branch of lower-mass stars, it is possible that lambda > 5. A high value of lambda (rather than assuming a high efficiency parameter, eta_CE >1) is sufficient to explain the long final orbital periods observed among those binary millisecond pulsars which are believed to have evolved through a CE-phase. We also present a tabulation of lambda as a function of stellar radius and mass, which is useful for a quick estimation of the orbital decay during a common envelope and spiral-in phase.

We have investigated the formation of the binary radio pulsars PSR B2303+46 and PSR J1141-6545 via Monte Carlo simulations of a large number of interacting stars in binary systems. PSR B2303+46 has recently been shown (van Kerkwijk & Kulkarni 1999) to be the first neutron star - white dwarf binary system observed, in which the neutron star was born after the formation of the white dwarf. We discuss the formation process for such a system and are able to put constraints on the parameters of the initial ZAMS binary. We present statistical evidence in favor of a white dwarf companion to the binary pulsar PSR J1141-6545, just recently discovered in the Parkes Multibeam Survey. If this is confirmed by observations this system will be the second one known in which the neutron star was born after its white dwarf companion. We also predict a minimum space velocity of 150 km/s for PSR J1141-6545, and show it must have experienced an asymmetric SN in order to explain its low eccentricity. Finally, we estimate the birthrate of these systems relative to other binary pulsar systems and present the expected distribution of their orbital periods, eccentricities and velocities.

We have performed detailed numerical calculations of the non-conservative evolution of close X-ray binary systems with intermediate-mass (2.0-6.0 M_sun) donor stars and a 1.3 M_sun accreting neutron star. We calculated the thermal response of the donor star to mass loss, in order to determine its stability and follow the evolution of the mass transfer. Under the assumption of the "isotropic re-emission model" we demonstrate that in many cases it is possible for the binary to prevent a spiral-in and survive a highly super-Eddington mass-transfer phase (1 << M_dot/M_Edd < 10^5) on a sub-thermal timescale, if the convective envelope of the donor star is not too deep. These systems thus provide a new formation channel for binary millisecond pulsars with heavy CO white dwarfs and relatively short orbital periods (3-50 days). However, we conclude that to produce a binary pulsar with a O-Ne-Mg white dwarf or P_orb ~1 day (e.g. PSR B0655+64) the above scenario does not work, and a spiral-in phase is still considered the most plausible scenario for the formation of such a system.

We present an updated list of direct strong evidence in favour of kicks being imparted to newborn neutron stars. In particular we discuss the new cases of evidence resulting from recent observations of the X-ray binary Circinus X-1 and the newly discovered binary radio pulsar PSR J1141-6545. We conclude that the assumption that neutron stars receive a kick velocity at their formation is unavoidable (van den Heuvel & van Paradijs 1997).

We discuss the influence of tidal spin-orbit interactions on the orbital dynamics of close intermediate-mass X-ray binaries. In particular we consider here a process in which spin angular momentum of a contracting RLO donor star, in a synchronous orbit, is converted into orbital angular momentum and thus helps to stabilize the mass transfer by widening the orbit. Binaries which would otherwise suffer from dynamically unstable mass transfer (leading to the formation of a common envelope and spiral-in evolution) are thus shown to survive a phase of extreme mass transfer on a sub-thermal timescale. Furthermore, we discuss the orbital evolution prior to RLO in X-ray binaries with low-mass donors, caused by the competing effects of wind mass loss and tidal effects due to expansion of the (sub)giant.

Both the recently observed high runaway velocities of Cyg X-1 (~50 km/s) and X-ray Nova Sco 1994 (> ~100 km/s) and the relatively low radial velocities of the black hole X-ray binaries with low mass donor stars, can be explained by symmetric mass ejection in the supernovae (SNe) which formed the black holes in these systems. Assuming symmetric mass ejection in black hole formation, we estimate the amount of mass that must have been ejected in the stellar core collapse in order to explain the velocities of the above X-ray binaries. We find that at least 2.6 M_sun and 4.1 M_sun must have been ejected in the formation of Cyg X-1 and Nova Sco, respectively. A similar mass loss fraction (f = 0.35) for the black hole binaries with low mass donors, gives low velocities, in agreement with the observations. We conclude that the black holes in X-ray binaries are all consistent with being formed in a successful SN in which mass is ejected. A possible kick at the formation of the black hole is not needed to explain their space velocities.

We have performed detailed numerical calculations of the non-conservative evolution of close binary systems with low-mass (1.0-2.0 M_sun) donor stars and a 1.3 M_sun accreting neutron star. Rather than using analytical expressions for simple polytropes, we calculated the thermal response of the donor star to mass loss, in order to determine the stability and follow the evolution of the mass transfer. Tidal spin-orbit interactions and Reimers wind mass-loss were also taken into account. We have re-calculated the correlation between orbital period and white dwarf mass in wide binary radio pulsar systems. Furthermore, we find an anti-correlation between orbital period and neutron star mass under the assumption of the "isotropic re-emission" model and compare this result with observations. We conclude that the accretion efficiency of neutron stars is rather low and that they eject a substantial fraction of the transferred material even when accreting at a sub-Eddington level. The mass-transfer rate is a strongly increasing function of initial orbital period and donor star mass. For relatively close systems with light donors (P < 10 days and M_2 < 1.3 M_sun) the mass-transfer rate is sub-Eddington, whereas it can be highly super-Eddington by a factor of 10^4 for wide systems with relatively heavy donor stars (1.6 - 2.0 M_sun) as a result of their deep convective envelopes. We briefly discuss the evolution of X-ray binaries with donor stars in excess of 2 M_sun. Based on our calculations we present evidence that PSR J1603-7202 evolved through a phase with unstable mass transfer from a relatively heavy donor star and therefore is likely to host a CO white dwarf companion.

We have analyzed the kinematical parameters of Cir X-1 to constrain the nature of its companion star, the eccentricity of the binary and the pre-supernova parameter space. We argue that the companion is most likely to be a low-mass (< 2.0 M_sun) unevolved star and that the eccentricity of the orbit is 0.94 +/- 0.04. We have evaluated the dynamical effects of the supernova explosion and we find it must have been asymmetric. On average, we find that a kick of 740 km/s is needed to account for the recently measured radial velocity of +430 km/s (Johnston, Fender & Wu) for this extreme system. The corresponding minimum kick velocity is 500 km/s. This is the largest kick needed to explain the motion of any observed binary system. If Cir X-1 is associated with the supernova remnant G321.9-0.3 then we find a limiting minimum age of this remnant of 60000 yr. Furthermore, we predict that the companion star has lost 10% of its mass as a result of stripping and ablation from the impact of the supernova shell shortly after the explosion.

We propose a scenario to form low-mass, single, slow rotating white dwarfs from a solar-like star accompanied by a massive planet, or a brown dwarf, in a relatively close orbit (e.g. HD 89707). Such white dwarfs were recently found by Maxted & Marsh (1998). When the solar-like star ascends the giant branch it captures the planet and the subsequent spiral-in phase expels the envelope of the giant leaving a low-mass helium white dwarf remnant. In case the planet evaporizes, or fills its own Roche-lobe, the outcome is a single undermassive white dwarf. The observed distribution of planetary systems supports the applicability of this scenario.

If the theoretical relationship between white dwarf mass and orbital period for wide-orbit binary radio pulsars is assumed to be correct, then the neutron star mass of PSR J2019+2425 is shown to be ~1.20 M_sun. Hence the mass of the neutron star in this system prior to the mass transfer phase is expected to have been < 1.1 M_sun. Alternatively this system descends from the accretion induced collapse (AIC) of a massive white dwarf. We estimate the magnetic inclination angles of all the observed wide-orbit low-mass binary pulsars in the Galactic disk using the core-mass period relation and assuming that the spin axis of an accreting neutron star aligns with the orbital angular momentum vector in the recycling process of the pulsar. The large estimated magnetic inclination angle of PSR J2019+2425, in combination with its old age, gives for this system evidence against alignment of the magnetic field axis with the rotational spin axis. However, in the majority of the similar systems the distribution of magnetic inclination angles is concentrated toward low values (if the core-mass period relation is correct) and suggests that alignment has taken place.