Double Neutron Star Systems Research Articles

MeerKAT observations of pair-plasma induced birefringence in the double pulsar eclipses

ABSTRACT PSR J0737−3039A/B is unique among double neutron star systems. Its near-perfect edge-on orbit causes the fast spinning pulsar A to be eclipsed by the magnetic field of the slow spinning pulsar B. Using high-sensitivity MeerKAT radio observations combined with updated constraints on the system geometry, we studied the impact of these eclipses on the incident polarization properties of pulsar A. Averaging light curves together after correcting for the rotation of pulsar B revealed enormous amounts of circular polarization and rapid changes in the linear polarization position angle, which occur at phases where emission from pulsar A is partially transmitted through the magnetosphere of pulsar B. These behaviours confirm that the eclipse mechanism is the result of synchrotron absorption in a relativistic pair-plasma confined to the closed-field region of pulsar B’s truncated dipolar magnetic field. We demonstrate that changes in circular polarization handedness throughout the eclipses are directly tied to the average line of sight magnetic field direction of pulsar B, from which we unambiguously determine the complete magnetic and viewing geometry of the pulsar.

VLBA Astrometry of the Galactic Double Neutron Stars PSR J0509+3801 and PSR J1930–1852: A Preliminary Transverse Velocity Distribution of Double Neutron Stars and its Implications

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

SN 2023zaw: An Ultrastripped, Nickel-poor Supernova from a Low-mass Progenitor

We present SN 2023zaw—a subluminous (M r = −16.7 mag) and rapidly evolving supernova (t 1/2,r = 4.9 days), with the lowest nickel mass (≈0.002 M ⊙) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He i and Ca near-infrared emission lines with velocities of ∼10,000−12,000 km s−1. The late-time spectra show prominent narrow He i emission lines at ∼1000 km s−1, indicative of interaction with He-rich circumstellar material. SN 2023zaw is located in the spiral arm of a star-forming galaxy. We perform radiation-hydrodynamical and analytical modeling of the lightcurve by fitting with a combination of shock-cooling emission and nickel decay. The progenitor has a best-fit envelope mass of ≈0.2 M ☉ and an envelope radius of ≈50 R ⊙. The extremely low nickel mass and low ejecta mass (≈0.5 M ⊙) suggest an ultrastripped SN, which originates from a mass-losing low-mass He-star (zero-age main-sequence mass < 10 M ⊙) in a close binary system. This is a channel to form double neutron star systems, whose merger is detectable with LIGO. SN 2023zaw underscores the existence of a previously undiscovered population of extremely low nickel mass (<0.005 M ☉) stripped-envelope supernovae, which can be explored with deep and high-cadence transient surveys.

Discovery and timing of ten new millisecond pulsars in the globular cluster Terzan 5

We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856–1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected timing solutions for all ten discoveries, covering nearly two decades of archival observations from the Green Bank Telescope for all but one. Highlights include PSR J1748−2446ao which is an eccentric (e = 0.32) wide-orbit (orbital period Pb = 57.55 d) system. We were able to measure the rate of advance of periastron (ω̇) for this system allowing us to determine a total mass of 3.17 ± 0.02 M⊙. With a minimum companion mass (Mc) of ∼0.8 M⊙, PSR J1748−2446ao is a candidate double neutron star (DNS) system. If confirmed to be a DNS, it would be the fastest spinning pulsar (P = 2.27 ms) and the longest orbital period measured for any known DNS system. PSR J1748−2446ap has the second highest eccentricity for any recycled pulsar (e ∼ 0.905) and for this system we can measure the total mass (1.997 ± 0.006 M⊙) and estimate the pulsar and companion masses, (1.700−0.045+0.015 M⊙ and 0.294−0.014+0.046 M⊙, respectively). PSR J1748−2446ar is an eclipsing redback (minimum Mc ∼ 0.34 M⊙) system whose properties confirm it to be the counterpart to a previously published source identified in radio and X-ray imaging. We were also able to detect ω̇ for PSR J1748−2446au leading to a total mass estimate of 1.82 ± 0.07 M⊙ and indicating that the system is likely the result of Case A Roche lobe overflow. With these discoveries, the total number of confirmed pulsars in Terzan 5 is 49, the highest for any globular cluster so far. These discoveries further enhance the rich set of pulsars known in Terzan 5 and provide scope for a deeper understanding of binary stellar evolution, cluster dynamics and ensemble population studies.

The Relativistic Spin Precession in the Compact Double Neutron Star System PSR J1946+2052

Abstract We observe systematic profile changes in the visible pulsar of the compact double neutron star system PSR J1946+2052 using observations with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The interpulse of PSR J1946+2052 changed from a single-peak to a double-peak shape from 2018–2021. We attribute this evolution as the result of the relativistic spin precession of the pulsar. With the high sensitivity of FAST, we also measure significant polarization for the first time, allowing us to model this with the precessional rotating vector model. Assuming, to the first order, a circular hollow-cone-like emission beam pattern and taking the validity of general relativity (GR), we derive the binary’s orbital inclination angle ( 63 ∘ − 3 ∘ + 5 ∘ ) and pulsar’s spin geometry. The pulsar’s spin vector and the orbital angular momentum vector are found to be only slightly misaligned ( 0. ∘ 21 − 0. ∘ 10 + 0. ∘ 28 ). The quoted uncertainties do not reflect the systematic uncertainties introduced by our model assumptions. By simulating future observations of profile and polarization evolution, we estimate that we could constrain the precession rate within a 43% uncertainty in 9 yr. Hence, we suggest that the system’s profile evolution could be combined with precise pulsar timing to test GR in the future.

Electron-capture supernovae in NS + He star systems and the double neutron star systems

ABSTRACT Electron-capture-supernovae (EC-SNe) provide an alternative channel for producing neutron stars (NSs). They play an important role in the formation of double NS (DNS) systems and the chemical evolution of galaxies, and contribute to the NS mass distribution in observations. It is generally believed that EC-SNe originate from e-captures on $\rm ^{24}Mg$ and $\rm ^{20}Ne$ in the massive degenerate oxygen–neon (ONe) cores with masses close to the Chandrasekhar limit (MCh). However, the origin of EC-SNe is still uncertain. In this paper, we systematically studied the EC-SNe in NS + He star systems by considering the explosive oxygen burning that may occur in the near-MCh ONe core. We provided the initial parameter spaces for producing EC-SNe in the initial orbital period − initial He star mass (log$P_{\rm orb}^{\rm i}-M_{\rm He}^{\rm i}$) diagram, and found that both $M_{\rm He}^{\rm i}$ and minimum $P_{\rm orb}^{\rm i}$ for EC-SNe increase with metallicity. Then, by considering NS kicks added to the newborn NS, we investigated the properties of the formed DNS systems after the He star companions collapse into NSs, such as the orbital periods, eccentricities, and spin periods of recycle pulsars (Pspin), etc. The results show that most of the observed DNS systems can be produced by NS kicks of $\lesssim$50 km s−1. In addition, we found that NSs could accrete more material if the residual H envelope on the He star companions is considered, which can form the mildly recycled pulsars (Pspin ∼ 20 ms) in DNS systems.

The FAST Galactic Plane Pulsar Snapshot Survey – V. PSR J1901+0658 in a double neutron star system

ABSTRACT Double neutron star (DNS) systems offer excellent opportunities to test gravity theories. We report the timing results of PSR J1901+0658, the first pulsar discovered in the FAST Galactic Plane Pulsar Snapshot (GPPS) Survey. Based on timing observations by FAST over 5 yr, we obtain the phase-coherent timing solution and derive the precise measurements of its position, spin parameters, orbital parameters, and dispersion measure. It has a period of 75.7 ms, a period derivative of 2.169(6) × 10−19 s s−1, and a characteristic age of 5.5 Gyr. This pulsar is in an orbit with a period of 14.45 d and an eccentricity of 0.366. One post-Keplerian parameter, periastron advance, has been well-measured as being 0.00531(9) deg yr−1, from which the total mass of this system is derived to be 2.79(7) M⊙. The pulsar has the mass upper limit of 1.68 M⊙, so the lower limit for the companion mass is 1.11 M⊙. Because PSR J1901+0658 is a partially recycled pulsar in an eccentric binary orbit with such a large companion mass, it should be in a DNS system according to the evolution history of the binary system.

A Relativistic Double Neutron Star Binary PSR J1846-0513

We report the timing analysis of PSR J1846−0513, a pulsar discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in Commensal Radio Astronomy FAST Survey. The pulsar possesses a spin period of 23.36 ms and a spin-down rate ( Ṗ ) of 1.0106(3) × 10−18 s s−1, and it is located in an eccentric orbit (e ∼0.208) with an orbital period of 0.61 days. The characteristic age and surface magnetic field of the pulsar are found to be 366.62 Myr and 4.9178 × 109 G, respectively, indicating that it is a recycled pulsar. Using over two years of timing data, we measure the periastron advance ω̇ = 0.8956(8) deg yr−1. By assuming that this effect is purely relativistic, we have estimated the total mass M = 2.6287(35)M ⊙ and obtained an upper limit for the pulsar mass and a lower limit for the companion’s mass. Our results indicate that this is a double neutron star system.

Total r-process Yields of Milky Way Neutron Star Mergers

While it is now known that the mergers of double neutron star binary systems (NSMs) are copious producers of heavy elements, there remains much speculation about whether they are the sole or even principal site of rapid neutron-capture (r-process) nucleosynthesis, one of the primary ways in which heavy elements are produced. The occurrence rates, delay times, and galactic environments of NSMs hold sway over estimating their total contribution to the elemental abundances in the solar system and the Galaxy. Furthermore, the expected elemental yields of NSMs may depend on the merger parameters themselves—such as their stellar masses and radii—which are not currently considered in many galactic chemical evolution models. Using the characteristics of the observed sample of double neutron star (DNS) systems in the Milky Way as a guide, we predict the expected nucleosynthetic yields that a population of DNSs would produce upon merger, and we compare that nucleosynthetic signature to the heavy-element abundance pattern of solar system elements. We find that with our current models, the present DNS population favors the production of lighter r-process elements, while underproducing the heaviest elements relative to the solar system. This inconsistency could imply an additional site for the heaviest elements or a population of DNSs much different from that observed today.

PSR J2150+3427: A Possible Double Neutron Star System

PSR J2150+3427 is a 0.654 s pulsar discovered by the Commensal Radio Astronomy FAST Survey. From the follow-up observations, we find that the pulsar is in a highly eccentric orbit (e = 0.601) with an orbital period of 10.592 days and a projected semimajor axis of 25.488 lt-s. Using 2.7 yr of timing data, we also measured the rate of periastron advance = 0.0115(4) deg yr−1. An estimate for the total mass of the system using the gives M tot = 2.59(13)M ⊙, which is consistent with most of the known double neutron star (DNS) systems and one neutron star (NS)–white dwarf (WD) system named B2303+46. Combining with the mass function of the system gives the masses of M p < 1.67 and M c > 0.98 M ⊙ for the pulsar and the companion star, respectively. This constraint, along with the spin period and orbital parameters, suggests that it is possibly a DNS system, and we cannot entirely rule out the possibility of an NS–WD system. Future timing observations will vastly improve the uncertainty in , and are likely to allow the detection of additional relativistic effects, which can be used to modify the values of M p and M c . With a spin-down luminosity of = 5.07(6) × 1029 erg s−1, PSR J2150+3427 is a very low-luminosity pulsar, with only the binary pulsar J2208+4610 having a smaller .

Accretion spin-up of the massive component in the neutron star stripping model for short gamma-ray bursts

ABSTRACT In this paper, we use analytical methods to study the last stages of the double neutron star (NS) system evolution. Depending on the initial masses of the components, this evolution can occur either in the framework of the merging scenario or in the NS stripping model. The main new ingredient of this work, compared with previous calculations, is accounting for accretion spin-up of the massive component. This effect leads to a significant decrease in the duration of the stable mass transfer of matter in the stripping mechanism. Within the framework of the Newtonian approximation, we determine the boundary between the merging and stripping scenarios. It is shown that this boundary weakly depends on the total mass of the system and the specific form of the NS equation of state, and is determined mainly by the initial mass ratio of the components. The stripping scenario is realized at M2/M1 ≲ 0.8, so it should make a large contribution to the population of gravitational wave events from NS–NS coalescing binaries that are close to us, and their accompanying short gamma-ray bursts. Nevertheless, the value obtained requires further clarification, taking into account relativistic effects, possible non-conservative mass transfer, etc.

The MPIfR–MeerKAT Galactic Plane Survey – I. System set-up and early results

ABSTRACT Galactic plane radio surveys play a key role in improving our understanding of a wide range of astrophysical phenomena. Performing such a survey using the latest interferometric telescopes produces large data rates necessitating a shift towards fully or quasi-real-time data analysis with data being stored for only the time required to process them. We present here the overview and set-up for the 3000-h Max-Planck-Institut für Radioastronomie (MPIfR)–MeerKAT Galactic Plane Survey (MMGPS). The survey is unique by operating in a commensal mode, addressing key science objectives of the survey including the discovery of new pulsars and transients and studies of Galactic magnetism, the interstellar medium and star formation rates. We explain the strategy coupled with the necessary hardware and software infrastructure needed for data reduction in the imaging, spectral, and time domains. We have so far discovered 78 new pulsars including 17 confirmed binary systems of which two are potential double neutron star systems. We have also developed an imaging pipeline sensitive to the order of a few tens of micro-Jansky ($\mu{\rm Jy}$) with a spatial resolution of a few arcseconds. Further science operations with an in-house built S-band receiver operating between 1.7 and 3.5 GHz are about to commence. Early spectral line commissioning observations conducted at S-band, targeting transitions of the key molecular gas tracer CH at 3.3 GHz already illustrate the spectroscopic capabilities of this instrument. These results lay a strong foundation for future surveys with telescopes like the Square Kilometre Array (SKA).

Double neutron star formation via consecutive type II supernova explosions

ABSTRACT Since the discovery of the first double neutron star (DNS) system, the number of these exotic binaries has reached 15. Here we investigate a channel of DNS formation in binary systems with components above the mass limit of Type II supernova explosion (SN II), that is, 8M⊙. We apply a spherically symmetric homologous envelope expansion model to account for mass-loss, and follow the dynamical evolution of the system numerically with a high-precision integrator. The first SN occurs in a binary system whose orbital parameters are pre-defined, then, the homologous expansion model is applied again in the newly formed system. Analysing 1 658 880 models we find that DNS formation via subsequent SN II explosions requires a fine-tuning of the initial parameters. Our model can explain DNS systems with a separation greater than 2.95 au. The eccentricity of the DNS systems spans a wide range thanks to the orbital circularization effect due to the second SN II explosion. The eccentricity of the DNS is sensitive to the initial eccentricity of the binary progenitor and the orbital position of the system preceding the second explosion. In agreement with the majority of the observations of DNS systems, we find the system centre-of-mass velocities to be less than 60 km s−1. Neutron stars that become unbound in either explosion gain a peculiar velocity in the range of 0.02 − 240 km s−1. In our model, the formation of tight DNS systems requires a post-explosion orbit-shrinking mechanism, possibly driven by the ejected envelopes.

New constraints on the kinematic, relativistic, and evolutionary properties of the PSR J1757−1854 double neutron star system

ABSTRACT PSR J1757−1854 is one of the most relativistic double neutron star binary systems known in our Galaxy, with an orbital period of $P_\text{b}=4.4\, \text{h}$ and an orbital eccentricity of e = 0.61. As such, it has promised to be an outstanding laboratory for conducting tests of relativistic gravity. We present the results of a 6-yr campaign with the 100-m Green Bank and 64-m Parkes radio telescopes, designed to capitalize on this potential. We identify secular changes in the profile morphology and polarization of PSR J1757−1854, confirming the presence of geodetic precession and allowing the constraint of viewing geometry solutions consistent with General Relativity. We also update PSR J1757−1854’s timing, including new constraints of the pulsar’s proper motion, post-Keplerian parameters, and component masses. We conclude that the radiative test of gravity provided by PSR J1757−1854 is fundamentally limited to a precision of 0.3 per cent due to the pulsar’s unknown distance. A search for pulsations from the companion neutron star is also described, with negative results. We provide an updated evaluation of the system’s evolutionary history, finding strong support for a large kick velocity of $w\ge 280\, \rm{km\,s}^{-1}$ following the second progenitor supernova. Finally, we reassess PSR J1757−1854’s potential to provide new relativistic tests of gravity. We conclude that a 3-σ constraint of the change in the projected semimajor axis ($\dot{x}$) associated with Lense–Thirring precession is expected no earlier than 2031. Meanwhile, we anticipate a 3-σ measurement of the relativistic orbital deformation parameter δθ as soon as 2026.

Investigating the distribution of double neutron stars and unconventional component mass

ABSTRACTThe mass distribution of double neutron stars (DNSs) can help us to understand the evolution of binary systems, which can also contain information about DNSs induced by gravitational wave radiation. Henceforth, for the 25 pairs of DNS systems so far discovered by radio observations, we have only statistically analysed 13 systems with precisely measured masses of two-component NSs, by employing the classical Gaussian mixture models based on the Akaike and Bayesian information criteria. The result infers that the mass distributions for recycled and non-recycled NSs follow a double-Gaussian and single-Gaussian distribution, respectively. In the conventional scenario of DNS formation, the mass of the first formed recycled NS is higher than that of the second-formed companion NS. However, there exist some unconventional cases with a reverse mass ratio. To explore the phenomenon of this ‘unconventional component mass’ of DNSs, we employed Monte Carlo sampling. We find that the unconventional DNS systems occupy about 27.7 per cent of all DNSs, which could be a useful reference for the future gravitational wave detection by LIGO/Virgo. Furthermore, we discuss a possible explanation for the unconventional DNS mass ratio, which involves the recycling process and its birth rate. This also has an implication for the mass distribution of the DNS progenitors.

An observationally derived kick distribution for neutron stars in binary systems

ABSTRACT Understanding the natal kicks received by neutron stars (NSs) during formation is a critical component of modelling the evolution of massive binaries. Natal kicks are an integral input parameter for population synthesis codes, and have implications for the formation of double NS systems and their subsequent merger rates. However, many of the standard observational kick distributions that are used are obtained from samples created only from isolated NSs. Kick distributions derived in this way overestimate the intrinsic NS kick distribution. For NSs in binaries, we can only directly estimate the effect of the natal kick on the binary system, instead of the natal kick received by the NS itself. Here, for the first time, we present a binary kick distribution for NSs with low-mass companions. We compile a catalogue of 145 NSs in low-mass binaries with the best available constraints on proper motion, distance, and systemic radial velocity. For each binary, we use a three-dimensional approach to estimate its binary kick. We discuss the implications of these kicks on system formation, and provide a parametric model for the overall binary kick distribution, for use in future theoretical modelling work. We compare our results with other work on isolated NSs and NSs in binaries, finding that the NS kick distributions fit using only isolated pulsars underestimate the fraction of NSs that receive low kicks. We discuss the implications of our results on modelling double NS systems, and provide suggestions on how to use our results in future theoretical works.

The Green Bank North Celestial Cap Survey. VII. 12 New Pulsar Timing Solutions

We present timing solutions for 12 pulsars discovered in the Green Bank North Celestial Cap 350 MHz pulsar survey, including six millisecond pulsars (MSPs), a double neutron star (DNS) system, and a pulsar orbiting a massive white dwarf companion. Timing solutions presented here include 350 and 820 MHz Green Bank Telescope data from initial confirmation and follow-up, as well as a dedicated timing campaign spanning 1 ryr PSR J1122−3546 is an isolated MSP, PSRs J1221−0633 and J1317−0157 are MSPs in black widow systems and regularly exhibit eclipses, and PSRs J2022+2534 and J2039−3616 are MSPs that can be timed with high precision and have been included in pulsar timing array experiments seeking to detect low-frequency gravitational waves. PSRs J1221−0633 and J2039−3616 have Fermi Large Area Telescope gamma-ray counterparts and also exhibit significant gamma-ray pulsations. We measure proper motions for three of the MSPs in this sample and estimate their space velocities, which are typical compared to those of other MSPs. We have detected the advance of periastron for PSR J1018−1523 and therefore measure the total mass of the DNS system, m tot = 2.3 ± 0.3 M ⊙. Long-term pulsar timing with data spanning more than 1 yr is critical for classifying recycled pulsars, carrying out detailed astrometry studies, and shedding light on the wealth of information in these systems post-discovery.

Missing for 20 yr: MeerKAT Redetects the Elusive Binary Pulsar M30B

PSR J2140−2311B is a 13 ms pulsar discovered in 2001 in a 7.8 hr Green Bank Telescope observation of the core-collapsed globular cluster M30 and predicted to be in a highly eccentric binary orbit. This pulsar has eluded detection since then; therefore, its precise orbital parameters have remained a mystery until now. In this work, we present the confirmation of this pulsar using observations taken with the UHF receivers of the MeerKAT telescope as part of the TRAPUM Large Survey Project. Taking advantage of the beamforming capability of our backends, we have localized it, placing it 1.′2(1) from the cluster center. Our observations have enabled the determination of its orbit: It is highly eccentric (e = 0.879) with an orbital period of 6.2 days. We also measured the rate of periastron advance, . Assuming that this effect is fully relativistic, general relativity provides an estimate of the total mass of the system, M TOT = 2.53 ± 0.08 M ⊙, consistent with the lightest double neutron star systems known. Combining this with the mass function of the system gives the pulsar and companion masses of m p < 1.43 M ⊙ and m c > 1.10 M ⊙, respectively. The massive, undetected companion could either be a massive white dwarf or a neutron star. M30B likely formed as a result of a secondary exchange encounter. Future timing observations will allow the determination of a phase-coherent timing solution, vastly improving our uncertainty in and likely enabling the detection of additional relativistic effects, which will determine m p and m c .

The MSPSRπ catalogue: VLBA astrometry of 18 millisecond pulsars

ABSTRACT With unparalleled rotational stability, millisecond pulsars (MSPs) serve as ideal laboratories for numerous astrophysical studies, many of which require precise knowledge of the distance and/or velocity of the MSP. Here, we present the astrometric results for 18 MSPs of the ‘MSPSR$\pi$’ project focusing exclusively on astrometry of MSPs, which includes the re-analysis of three previously published sources. On top of a standardized data reduction protocol, more complex strategies (i.e. normal and inverse-referenced 1D interpolation) were employed where possible to further improve astrometric precision. We derived astrometric parameters using sterne, a new Bayesian astrometry inference package that allows the incorporation of prior information based on pulsar timing where applicable. We measured significant (${&amp;gt;}3\, \sigma$) parallax-based distances for 15 MSPs, including 0.81 ± 0.02 kpc for PSR J1518+4904 – the most significant model-independent distance ever measured for a double neutron star system. For each MSP with a well-constrained distance, we estimated its transverse space velocity and radial acceleration. Among the estimated radial accelerations, the updated ones of PSR J1012+5307 and PSR J1738+0333 impose new constraints on dipole gravitational radiation and the time derivative of Newton’s gravitational constant. Additionally, significant angular broadening was detected for PSR J1643−1224, which offers an independent check of the postulated association between the HII region Sh 2-27 and the main scattering screen of PSR J1643−1224. Finally, the upper limit of the death line of γ-ray-emitting pulsars is refined with the new radial acceleration of the hitherto least energetic γ-ray pulsar PSR J1730−2304.

From 3D hydrodynamic simulations of common-envelope interaction to gravitational-wave mergers

Modelling the evolution of progenitors of gravitational-wave merger events in binary stars faces two major uncertainties: the common-envelope phase and supernova kicks. These two processes are critical for the final orbital configuration of double compact-object systems with neutron stars and black holes. Predictive one-dimensional models of common-envelope interaction are lacking and multidimensional simulations are challenged by the vast range of relevant spatial and temporal scales. Here, we present three-dimensional, magnetohydrodynamic simulations of the common-envelope interaction of an initially 10 M⊙ red supergiant primary star with a black-hole and a neutron-star companion. Employing the moving-mesh code AREPO and replacing the core of the primary star and the companion with point masses, we show that the high-mass regime is accessible to full ab initio simulations. About half of the common envelope is dynamically ejected at the end of our simulations and the ejecta mass fraction keeps growing. Almost complete envelope ejection seems possible if all ionised gas left over at the end of our simulation eventually recombines and the released energy continues to help unbind the envelope. We find that the dynamical plunge-in of both companions terminates at orbital separations that are too wide for gravitational waves to merge the systems in a Hubble time. However, the orbital separations at the end of our simulations are still decreasing such that the true final value at the end of the common-envelope phase remains uncertain. We discuss the further evolution of the system based on analytical estimates. A subsequent mass-transfer episode from the remaining 3 M⊙ core of the supergiant to the compact companion does not shrink the orbit sufficiently either. A neutron-star–neutron-star and neutron-star–black-hole merger is still expected for a fraction of the systems if the supernova kick aligns favourably with the orbital motion. For double neutron star (neutron-star–black-hole) systems we estimate mergers in about 9% (1%) of cases while about 77% (94%) of binaries are disrupted; that is, supernova kicks actually enable gravitational-wave mergers in the binary systems studied here. Assuming orbits smaller by one-third after the common-envelope phase enhances the merger rates by about a factor of two. However, the large post-common-envelope orbital separations found in our simulations mean that a reduction in predicted gravitational-wave merger events appears possible.