Search

Stellar-mass black holes (BHs) embedded in active galactic nuclei (AGN) may be major sources of astrophysical gravitational waves (GWs), contributing both to the observed LIGO-Virgo-KAGRA population of binary BH mergers and to future populations of LISA-band extreme mass ratio inspirals (EMRIs). The ability of these BHs to pair up into binaries, inspiral, and produce GWs will be shaped by the existence of migration traps, regions in the AGN where hydrodynamic torques vanish. Previous works have studied the existence and location of migration traps in AGN disks. Here, we investigate how individual BHs may escape such traps as an outcome of mergers, potentially suppressing hierarchical growth. We find that while GW recoil kicks are strong enough to kick merged BHs onto inclined orbits, gas drag quickly realigns them into the AGN disk. A more robust escape mechanism is gap opening: once a BH grows above a critical mass, its gravity disturbs the AGN gas sufficiently to eliminate the trap. In low-mass AGN relevant for LISA, gaps open easily and the resulting ``wet EMRI'' masses are unlikely to reflect protracted hierarchical mergers. In combination with our previous work, we find that migration traps only exist in a relatively narrow range of AGN luminosities between [10^{43.5},10^{45.5}] erg/s. We identify an even narrower AGN luminosity range for which stellar mass BHs can grow into the pair instability mass gap and beyond. This characteristic luminosity scale may assist in indirect tests of the ``AGN channel'' for binary BH mergers., Comment: "8 pages, 4 figures in main text. Comments welcome!

In this invited talk at the 41st Li\`ege International Astrophysical Colloquium on "The eventful life of massive star multiples", I reviewed some aspects of our current understanding of neutron stars and black holes as end products of stellar evolution as well as the evolutionary paths leading to the formation of high-mass X-ray binaries.

We construct a semi-analytical model that describes the convective core mass evolution of massive stars experiencing mass loss during the main-sequence stage. We first conduct a suite of 1D stellar evolution calculations to build insight into how convective core masses behave under idealized mass loss. Based on these simulations, we find several universal relations between global properties of the star that hold regardless of the mass loss history. By combining these relations, we construct a semi-analytic framework that can predict the convective core mass evolution for arbitrary mass loss histories and hence the helium core mass at the end of the main sequence. Our formulae improve upon existing methods for predicting the core mass in rapid population synthesis codes., Comment: 12 pages, 10 figures

The latest binary black hole population estimates argue for a subpopulation of unequal component mass binaries with spins that are likely small but isotropically distributed. This implies a non-zero probability of detecting spin-induced orbital precession and higher order multipoles moments in the observed gravitational-wave signals. In this work we directly calculate the probability for precession and higher order multipoles in each significant gravitational-wave candidate observed by the LIGO--Virgo--KAGRA collaborations (LVK). We find that only one event shows substantial evidence for precession: GW200129_065458, and two events show substantial evidence for higher order multipoles: GW190412 and GW190814; any evidence for precession and higher order multipoles in other gravitational-wave signals is consistent with random fluctuations caused by noise. We then compare our observations with expectations from population models, and confirm that current population estimates from the LVK accurately predict the number of observed events with significant evidence for precession and higher order multipole moments. In particular, we find that this population model predicts that a binary with significant evidence for precession will occur once in every $\sim 50$ detections, and a binary with significant evidence for higher order multipole moments will occur once in every $\sim 70$ observations. However, we emphasise that since substantial evidence for precession and higher order multipole moments have only been observed in three events, any population model that includes a subpopulation of binaries yielding $\sim 2\%$ of events with detectable precession and higher order multipole moments will likely be consistent with the data., Comment: 13 pages, 5 figures

The three dynamically confirmed wind-fed black hole high-mass X-ray binaries (BH-HMXBs) are suggested to all contain a highly spinning black hole (BH). However, based on the theories of efficient angular momentum transport inside the stars, we expect that the first-born BHs in binary systems should have low spins, which is consistent with gravitational-wave observations. As a result, the origin of the high BH spins measured in wind-fed BH-HMXBs remains a mystery. In this paper, we conduct a binary population synthesis study on wind-fed BH-HMXBs at solar metallicity with the use of the newly developed code POSYDON, considering three scenarios for BH accretion: Eddington-limited, moderately super-Eddington, and fully conservative accretion. Taking into account the conditions for accretion-disk formation, we find that regardless of the accretion model, these systems are more likely to have already experienced a phase of Roche-lobe overflow after the BH formation. To account for the extreme BH spins, highly conservative accretion onto BHs is required, when assuming the accreted material carries the specific angular momentum at the innermost stable orbit. Besides, in our simulations we found that the systems with donor stars within the mass range of $10-20\,M_{\odot}$ are prevalent, posing a challenge in explaining simultaneously all observed properties of the BH-HMXB in our Galaxy, Cygnus X-1, and potentially hinting that the accretion efficiency onto non-degenerate stars, before the formation of the BH, is also more conservative than assumed in our simulations., Comment: Submitted to A&A. Comments are welcome!

Stellar models indicate that the core compactness of a star, which is a common proxy for its explodability in a supernova, does not increase monotonically with the star's mass. Rather, the core compactness dips sharply over a range of carbon-oxygen core masses; this range may be somewhat sensitive to the star's metallicity and evolutionary history. Stars in this compactness dip are expected to experience supernovae leaving behind neutron stars, whereas stars on either side of this range are expected to form black holes. This results in a hypothetical mass range in which black holes should seldom form. Quantitatively, when applied to binary stripped stars, these models predict a dearth of binary black holes with component masses $\approx 10 M_\odot - 15 M_\odot$. The population of gravitational-wave signals indicates potential evidence for a dip in the distribution of chirp masses of merging binary black holes near $\approx 10 M_\odot - 12 M_\odot$. This feature could be linked to the hypothetical component mass gap described above, but this interpretation depends on what assumptions are made of the binaries' mass ratios. Here, we directly probe the distribution of binary black hole component masses to look for evidence of a gap. We find no evidence for this feature using data from the third gravitational-wave transient catalogue (GWTC-3). If this gap does exist in nature, we find that it is unlikely to be resolvable by the end of the current (fourth) LIGO-Virgo-KAGRA (LVK) observing run., Comment: 16 pages, 11 figures

Wolf-Rayet stars in close binary systems can be tidally spun up by their companions, potentially leaving behind fast-spinning highly-magnetized neutron stars, known as ``magnetars", after core collapse. These newborn magnetars can transfer rotational energy into heating and accelerating the ejecta, producing hydrogen-poor superluminous supernovae (SLSNe). In this {\em{Letter}}, we propose that the magnetar wind of the newborn magnetar could significantly evaporate its companion star, typically a main-sequence or helium star, if the binary system is not disrupted by the {abrupt mass loss and} SN kick. The subsequent heating and acceleration of the evaporated star material along with the SN ejecta by the magnetar wind can produce a post-peak bump in the SLSN lightcurve. Our model can reproduce the primary peaks and post-peak bumps of four example observed multiband SLSN lightcurves, revealing that the mass of the evaporated material could be $\sim0.4-0.6\,M_\odot$ if the material is hydrogen-rich. {We propose that the magnetar could induce strongly enhanced evaporation from its companion star near the pericenter if the orbit of the post-SN binary is highly eccentric, ultimately generating multiple post-peak bumps in the SLSN lightcurves. This ``magnetar-star binary engine" model may offer a possible explanation for the evolution of polarization, along with the origin and velocity broadening of late-time hydrogen or helium broad spectral features observed in some bumpy SLSNe.} The diversity in the lightcurves and spectra of SLSNe may be attributed to the wide variety of companion stars and post-SN binary systems., Comment: Accepted by ApJL, 14 pages, 4 figures, 2 tables

Stars falling too close to massive black holes in the centres of galaxies can be torn apart by the strong tidal forces. Simulating the subsequent feeding of the black hole with disrupted material has proved challenging because of the range of timescales involved. Here we report a set of simulations that capture the relativistic disruption of the star, followed by one year of evolution of the returning debris stream. These reveal the formation of an expanding asymmetric bubble of material extending to hundreds of astronomical units -- an outflowing Eddington envelope with an optically thick inner region. Such envelopes have been hypothesised as the reprocessing layer needed to explain optical/UV emission in tidal disruption events, but never produced self-consistently in a simulation. Our model broadly matches the observed light curves with low temperatures, faint luminosities, and line widths of 10,000--20,000 km/s., Comment: 12 pages + appendices, 15 figures, accepted to ApJL. Data+movies on zenodo at https://doi.org/10.5281/zenodo.11438154 . Movies on YouTube... watch this one first: https://www.youtube.com/watch?v=KlffYCtYxFU then: https://www.youtube.com/watch?v=bPurbeQNgvI&t=19s and finally watch the expanding "bubble" with rho and T at photosphere: https://www.youtube.com/watch?v=D4dyZeUbfvM

We explore the evolution of massive stars (>8 solar masses) with 1-D models and present analytical fits to the masses and binding energies of the convective portions of their envelopes. These fits are given as functions of total mass, metallicity, and surface temperature (used as a proxy for evolutionary phase). They enable the application of the two-stage common envelope formalism (Hirai & Mandel 2022) in rapid binary population synthesis frameworks. We estimate that the degree of orbital hardening following common-envelope ejection spans 6 orders of magnitude and is a very strong function of the accretor mass, and, to a lesser extent, donor evolutionary phase., Comment: Minor changes (including to fits) to match version accepted for publication in ApJ

Accreting main-sequence stars expand significantly when the mass accretion timescale is much shorter than their thermal timescales. This occurs during mass transfer from an evolved giant star onto a main-sequence companion in a binary system, and is an important phase in the formation of compact binaries including X-ray binaries, cataclysmic variables, and gravitational-wave sources. In this study, we compute 1D stellar models of main-sequence accretors with different initial masses and accretion rates. The calculations are used to derive semi-analytical approximations to the maximum expansion radius. We assume that mass transfer remains fully conservative as long as the inflated accretor fits within its Roche lobe, leading stars to behave like hamsters, stuffing excess material behind their expanding cheeks. We suggest a physically motivated prescription for the mass growth of such "hamstars", which can be used to determine mass-transfer efficiency in rapid binary population synthesis models. With this prescription, we estimate that progenitors of high-mass X-ray binaries and gravitational-wave sources may have experienced highly non-conservative mass transfer. In contrast, for low-mass accretors, the accretion timescale can exceed the thermal timescale by a larger factor without causing significant radial expansion., Comment: 12 pages, 6 figures. Accepted for publication in ApJ Letters

Stars approaching supermassive black holes can be tidally disrupted. Despite being expected to emit X-rays, TDEs have been largely observed in optical bands, which is poorly understood. In this Letter, we simulate the tidal disruption of a $1~M_\odot$ main sequence star on an eccentric ($e=0.95$) orbit with a periapsis distance one or five times smaller than the tidal radius ($\beta = 1$ or $5$) using general relativistic smoothed particle hydrodynamics. We follow the simulation for up to a year post-disruption. We show that accretion disks in eccentric TDEs are masked by unbound material outflowing at $\sim10,000~$km/s. Assuming electron scattering opacity, this material would be visible as a $\sim100~$au photosphere at $\sim10^4~$K, in line with observations of candidate TDEs., Comment: 11 pages, 5 figures, submitted to ApJL

As gravitational wave detectors improve, observations of black hole (BH) mergers will provide the joint distribution of their masses and spins. This will be a critical benchmark to validate formation scenarios. Merging binary BHs formed through isolated binary evolution require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped stars are constrained by their surface critical rotation, restricting their angular momentum content at core-collapse. We use stripped star models at low metallicities ($Z_\odot/10$, $Z_\odot/50$ and $Z_\odot/250$) to determine the spins of BHs produced by critically rotating stellar progenitors. To study how such progenitors can arise, we consider their formation through chemically homogeneous evolution (CHE). We use a semianalytical model to study the final spins of CHE binaries, and compare our results against available detailed population synthesis models. We find that above BH masses of $\simeq 25M_\odot$, the dimensionless spin of critically rotating stripped stars ($a = Jc/(GM^2$)) is below unity. This results in an exclusion region at high chirp masses and effective spins that cannot be populated by binary evolution. CHE can produce binaries where both BHs hit this limit, producing a pile-up at the boundary of the excluded region. Highly spinning BHs arise from very low-metallicity CHE systems with short delay times, which merge at higher redshifts. On the other hand, the contribution of CHE to merging binary BHs in the third observing run of the LVK collaboration is expected to be dominated by systems with low spins ($\chi_\mathrm{eff}<0.5$) which merge near redshift zero. Owing to its higher projected sensitivity and runtime, the fourth observing run of the LVK collaboration can potentially place constraints on the high spin population and the existence of a limit set by critical rotation., Comment: Submitted to A&A. Abstract abridged

LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome., Comment: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

GW190521 is a remarkable gravitational-wave signal on multiple fronts: its source is the most massive black hole binary identified to date and could have spins misaligned with its orbit, leading to spin-induced precession -- an astrophysically consequential property linked to the binary's origin. However, due to its large mass, GW190521 was only observed during its final 3-4 cycles, making precession constraints puzzling and giving rise to alternative interpretations, such as eccentricity. Motivated by these complications, we trace the observational imprints of precession on GW190521 by dissecting the data with a novel time domain technique, allowing us to explore the morphology and interplay of the few observed cycles. We find that precession inference hinges on a quiet portion of the pre-merger data that is suppressed relative to the merger-ringdown. Neither pre-merger nor post-merger data alone are the sole driver of inference, but rather their combination: in the quasi-circular scenario, precession emerges as a mechanism to accommodate the lack of a stronger pre-merger signal in light of the observed post-merger. In terms of source dynamics, the pre-merger suppression arises from a tilting of the binary with respect to the observer. Establishing such a consistent picture between the source dynamics and the observed data is crucial for characterizing the growing number of massive binary observations and bolstering the robustness of ensuing astrophysical claims., Comment: 11 pages (excluding references), 9 figures

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

We use the rapid binary population synthesis code COMPAS to investigate commonly used prescriptions for the determination of mass transfer stability in close binaries and the orbital separations after stable mass transfer. The degree of orbital tightening during non-conservative mass transfer episodes is governed by the poorly-constrained angular momentum carried away by the ejected material. Increased orbital tightening drives systems towards unstable mass transfer leading to a common envelope. We find that the fraction of interacting binaries that will undergo only stable mass transfer throughout their lives fluctuates between a few and $\sim 20\%$ due to uncertainty in the angular momentum loss alone. If mass transfer is significantly non-conservative, stability prescriptions that rely on the assumption of conservative mass transfer under-predict the number of systems which experience unstable mass transfer and stellar mergers. This may substantially impact predictions about the rates of various transients, including luminous red novae, stripped-envelope supernovae, X-ray binaries, and the progenitors of coalescing compact binaries., Comment: 16 pages, 7 figures

Although there are estimated to be 100 million isolated black holes (BHs) in the Milky Way, only one has been found so far, resulting in significant uncertainty about their properties. The Galactic Bulge Time Domain Survey provides the only opportunity in the coming decades to grow this catalog by order(s) of magnitude. This can be achieved if 1) Roman's astrometric potential is fully realized in the observation strategy and software pipelines, 2) Roman's observational gaps of the Bulge are minimized, and 3) observations with ground-based facilities are taken of the Bulge to fill in gaps during non-Bulge seasons. A large sample of isolated BHs will enable a broad range of astrophysical questions to be answered, such as massive stellar evolution, origin of gravitational wave sources, supernova physics, and the growth of supermassive BHs, maximizing Roman's scientific return., Comment: 20 pages. Submitted in response to Nancy Grace Roman Space Telescope white paper call: https://roman.gsfc.nasa.gov/science/ccs_white_papers.html. v2 fixes a typo in Figure 5 axis label (days --> year)

The collapse of massive stars is one of the most-studied paths to black hole formation. In this chapter, we review black hole formation during the collapse of massive stars in the broader context of single and binary stellar evolution and the theory of supernova explosions. We provide a concise overview of the evolutionary channels that may lead to black hole formation -- the classical route of iron core collapse, collapse due to pair instability in very massive stars, and the hypothetical scenario of supermassive star collapse. We then review the current understanding of the parameter space for black hole formation and black hole birth properties that has emerged from theoretical and computational modelling of supernova explosions and transient observations. Finally, we discuss what the intricate interplay between stellar evolution, stellar explosions, and binary interactions implies for the formation of stellar-mass black holes., Comment: 53 pages, 9 figures. This chapter is the pre-print of the version currently in production. Please cite this chapter as the following: A.Heger, B. M\"uller, and I. Mandel. "Black holes as the end state of stellar evolution: Theory and simulations," in The Encyclopedia of Cosmology (Set 2): Black Holes, edited by Z. Haiman (World Scientific, New Jersey, 2023)

The majority of long duration ($>2$ s) gamma-ray bursts (GRBs) are believed to arise from the collapse of massive stars \cite{Hjorth+03}, with a small proportion created from the merger of compact objects. Most of these systems are likely formed via standard stellar evolution pathways. However, it has long been thought that a fraction of GRBs may instead be an outcome of dynamical interactions in dense environments, channels which could also contribute significantly to the samples of compact object mergers detected as gravitational wave sources. Here we report the case of GRB 191019A, a long GRB (T_90 = 64.4 +/- 4.5 s) which we pinpoint close (<100 pc projected) to the nucleus of an ancient (>1~Gyr old) host galaxy at z=0.248. The lack of evidence for star formation and deep limits on any supernova emission make a massive star origin difficult to reconcile with observations, while the timescales of the emission rule out a direct interaction with the supermassive black hole in the nucleus of the galaxy, We suggest that the most likely route for progenitor formation is via dynamical interactions in the dense nucleus of the host, consistent with the centres of such galaxies exhibiting interaction rates up to two orders of magnitude larger than typical field galaxies. The burst properties could naturally be explained via compact object mergers involving white dwarfs (WD), neutron stars (NS) or black holes (BH). These may form dynamically in dense stellar clusters, or originate in a gaseous disc around the supermassive black hole. Future electromagnetic and gravitational-wave observations in tandem thus offer a route to probe the dynamical fraction and the details of dynamical interactions in galactic nuclei and other high density stellar systems., Comment: Accepted to Nature Astronomy. This is the submitted version and will differ from the published version due to modifications in the refereeing process

We conduct binary population synthesis to investigate the formation of wind-fed high-mass X-ray binaries containing black holes (BH-HMXBs). We evolve multiple populations of high-mass binary stars and consider BH-HMXB formation rates, masses, spins and separations. We find that systems similar to Cygnus X-1 likely form after stable Case A mass transfer (MT) from the main sequence progenitors of black holes, provided such MT is characterised by low accretion efficiency, $\beta \lesssim 0.1$, with modest orbital angular momentum losses from the non-accreted material. Additionally, efficient BH-HMXB formation relies on a new simple treatment for Case A MT that allows donors to retain larger core masses compared to traditional rapid population-synthesis assumptions. At solar metallicity, our Preferred model yields $\mathcal{O}(1)$ observable BH-HMXBs in the Galaxy today, consistent with observations. In this simulation, $8\%$ of BH-HMXBs go on to merge as binary black holes or neutron star-black hole binaries within a Hubble time; however, none of the merging binaries have BH-HMXB progenitors with properties similar to Cygnus X-1. With our preferred settings for core mass growth, mass transfer efficiency and angular momentum loss, accounting for an evolving metallicity, and integrating over the metallicity-specific star formation history of the Universe, we find that BH-HMXBs may have contributed $\approx2$--$5$ BBH merger signals to detections reported in the third gravitational-wave transient catalogue of the LIGO-Virgo-KAGRA Collaboration. We also suggest MT efficiency should be higher during stable Case B MT than during Case A MT., Comment: 17 pages, 10 figures

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., Comment: Accepted for publication in MNRAS. 28 pages, 19 figures, 8 tables

Rapidly growing catalogs of compact binary mergers from advanced gravitational-wave detectors allow us to explore the astrophysics of massive stellar binaries. Merger observations can constrain the uncertain parameters that describe the underlying processes in the evolution of stars and binary systems in population models. In this paper, we demonstrate that binary black hole populations - namely, detection rates, chirp masses, and redshifts - can be used to measure cosmological parameters describing the redshift-dependent star formation rate and metallicity distribution. We present a method that uses artificial neural networks to emulate binary population synthesis computer models, and construct a fast, flexible, parallelisable surrogate model that we use for inference., Comment: 16 pages, 3 tables, 8 figures

Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black-hole spins, can play an important role in the early chemical enrichment of the Milky Way. We also show that none of the binary population synthesis models used in this paper, i.e., COMPAS, StarTrack, Brussels, ComBinE, and BPASS, can currently reproduce the elemental abundance observations. The predictions are problematic not only for neutron star mergers, but also for Type Ia supernovae, which may point to shortcomings in binary evolution models., Comment: Accepted for publication in the Astrophysical Journal Letters, 15 pages, 4 figures, 1 table

Transit and radial-velocity surveys over the past two decades have uncovered a significant population of short-period exoplanets. Among them are hot Jupiters, which are gas giant planets with orbital periods of a few days and found in 0.1-1% of Sun-like stars. Hot Jupiters are expected to be engulfed during their host star's radial expansion on the red giant branch. Planetary engulfment has been studied extensively as it may account for observed rapidly rotating and chemically enriched giant stars. We perform 3D hydrodynamical simulations of hot Jupiter engulfment by a 1 solar mass, 4 solar radii early red giant. Our "global" simulations simultaneously resolve the stellar envelope and planetary structure, modelling the hot Jupiter as a polytropic gas sphere. We find that approximately 90% of the hot Jupiter's mass is ablated in the convective part of the giant envelope, which would enhance the surface lithium abundance by 0.1 dex. The hot Jupiter is disrupted by a combination of ram pressure and tidal forces near the base of the convective envelope, with the deepest material penetrating to the radiative zone. The star experiences modest spin-up (~1 km/s), although engulfing a more massive companion could produce a rapidly rotating giant. Drag heating near the surface could exceed the unperturbed stellar luminosity and power an optical transient. For the amount of unbound ejecta recorded in the simulation, H-recombination could also power a transient that is around ten times the pre-engulfment luminosity, for several days., Comment: 14 pages, 11 figures

The speeds of young isolated pulsars are generally inferred from their observed 2-d velocities on the plane of the sky under the assumption that the unobserved radial velocity is not special, i.e., that the measured 2-d velocity is an isotropic projection of the full 3-d velocity. However, if pulsar spins are preferentially aligned with kicks, then the observer's viewing angle relative to the pulsar velocity vector is in fact special because the direction of the spin impacts the detectability of the pulsar. This means that the measured 2-d velocity of observable pulsars is not an isotropic projection, which affects inference on 3-d velocities. We estimate this effect and conclude that it could lead to a ~15% systematic over-estimate of neutron star natal kicks if young pulsars have high obliquity angles and narrow beams, but the exact correction factor depends on the distribution of beam-spin and spin-kick misalignment angles and beam widths., Comment: Updated to published version

Current prescriptions for supernova natal kicks in rapid binary population synthesis simulations are based on fits of simple functions to single pulsar velocity data. We explore a new parameterization of natal kicks received by neutron stars in isolated and binary systems developed by Mandel & M\"uller, which is based on 1D and 3D supernova simulations and accounts for the physical correlations between progenitor properties, remnant mass, and the kick velocity. We constrain two free parameters in this model using very long baseline interferometry velocity measurements of Galactic single pulsars. We find that the inferred values of natal kick parameters do not differ significantly between single and binary evolution scenarios. The best-fit values of these parameters are $v_{\rm ns} = 520$ km s$^{-1}$ for the scaling pre-factor for neutron star kicks, and $\sigma_{\rm ns}=0.3$ for the fractional stochastic scatter in the kick velocities., Comment: 8 pages, 8 figures, 1 table. Updated to match MNRAS version

We propose a new simple formalism to predict the orbital separations after common-envelope phases with massive star donors. We focus on the fact that massive red supergiants tend to have a sizeable radiative layer between the dense helium core and the convective envelope. Our formalism treats the common-envelope phase in two stages: dynamical in-spiral through the outer convective envelope and thermal timescale mass transfer from the radiative intershell. With fiducial choices of parameters, the new formalism typically predicts much wider separations compared to the classical energy formalism. Moreover, our formalism predicts that final separations strongly depend on the donor evolutionary stage and companion mass. Our formalism provides a physically-motivated alternative option for population synthesis studies to treat common-envelope evolution. This treatment will impact on predictions for massive-star binaries, including gravitational-wave sources, X-ray binaries and stripped-envelope supernovae., Comment: 10 pages, 5 figures, resubmitted to AAS journals

The three-body problem is famously chaotic, with no closed-form analytical solutions. However, hierarchical systems of three or more bodies can be stable over indefinite timescales. A system is considered hierarchical if the bodies can be divided into separate two-body orbits with distinct time- and length-scales, such that one orbit is only mildly affected by the gravitation of the other bodies. Previous work has mapped the stability of such systems at varying resolutions over a limited range of parameters, and attempts have been made to derive analytic and semi-analytic stability boundary fits to explain the observed phenomena. Certain regimes are understood relatively well. However, there are large regions of the parameter space which remain un-mapped, and for which the stability boundary is poorly understood. We present a comprehensive numerical study of the stability boundary of hierarchical triples over a range of initial parameters. Specifically, we consider the mass ratio of the inner binary to the outer third body ($q_{\rm out}$), mutual inclination ($i$), initial mean anomaly and eccentricity of both the inner and outer binaries ($e_{\rm in}$ and $e_{\rm out}$ respectively). We fit the dependence of the stability boundary on $q_{\rm out}$ as a threshold on the ratio of the inner binary's semi-major axis to the outer binary's pericentre separation $a_{\rm in}/R_{\rm p, out} \leq 10^{-0.6 + 0.04q_{\rm out}}q_{\rm out}^{0.32+0.1q_{\rm out}}$ for coplanar prograde systems. We develop an additional factor to account for mutual inclination. The resulting fit predicts the stability of $10^4$ orbits randomly initialised close to the stability boundary with $87.7\%$ accuracy., Comment: Accepted to PASA

The role of recombination during a common-envelope event has been long debated. Many studies have argued that much of hydrogen recombination energy, which is radiated in relatively cool and optically-thin layers, might not thermalise in the envelope. On the other hand, helium recombination contains 30% of the total recombination energy, and occurs much deeper in the stellar envelope. We investigate the distinct roles played by hydrogen and helium recombination in a common-envelope interaction experienced by a 12 solar mass red supergiant donor. We perform adiabatic, 3D hydrodynamical simulations that (i) include hydrogen, helium, and molecular hydrogen recombination, (ii) include hydrogen and helium recombination, (iii) include only helium recombination, and (iv) do not include recombination energy. By comparing these simulations, we find that the addition of helium recombination energy alone ejects 30% more envelope mass, and leads to a 16% larger post-plunge-in separation. Under the adiabatic assumption, adding hydrogen recombination energy increases the amount of ejected mass by a further 40%, possibly unbinding the entire envelope, but does not affect the post-plunge separation. Most of the ejecta becomes unbound at relatively high (>70%) degrees of hydrogen ionisation, where the hydrogen recombination energy is likely to expand the envelope instead of being radiated away., Comment: 10 pages and 6 figures. Accepted for publication in MNRAS

Triple stars and compact objects are ubiquitously observed in nature. Their long-term evolution is complex; in particular, the von-Zeipel-Lidov-Kozai (ZLK) mechanism can potentially lead to highly eccentric encounters of the inner binary. Such encounters can lead to a plethora of interacting binary phenomena, as well as stellar and compact-object mergers. Here we find explicit analytical formulae for the maximal eccentricity, $e_{\rm max}$, of the inner binary undergoing ZLK oscillations, where both the test particle limit (parametrised by the inner-to-outer angular momentum ratio $\eta$) and the double-averaging approximation (parametrised by the period ratio, $\epsilon_{\rm SA}$) are relaxed, for circular outer orbits. We recover known results in both limiting cases (either $\eta$ or $\epsilon_{\rm SA} \to 0$) and verify the validity of our model using numerical simulations. We test our results with two accurate numerical N-body codes, $\texttt{Rebound}$ for Newtonian dynamics and $\texttt{Tsunami}$ for general-relativistic (GR) dynamics, and find excellent correspondence. We discuss the implications of our results for stellar triples and both stellar and supermassive triple black hole mergers.

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, $Z$, and redshift $z, \mathcal{S}(Z,z)$. Considering these uncertainties we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal{S}(Z,z)$ variations can impact the predicted intrinsic and detectable merger rates by factors $10^2$-$10^4$. We find that BHBH rates are dominantly impacted by $\mathcal{S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal{S}(Z,z)$ for BHBH, BHNS and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95% of BHBH detections to contain a BH $\gtrsim 8\,\rm{M}_{\odot}$ and have mass ratios $\lesssim 4$. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties., Comment: 25 pages of which 12 are (rainbow) figures. Main results: Fig 2, 4, 5 and 6. Code/results publicly available at https://github.com/FloorBroekgaarden/Double-Compact-Object-Mergers. Comments welcome

The next generation of ground-based gravitational-wave detectors will observe coalescences of black holes and neutron stars throughout the cosmos, thousands of them with exceptional fidelity. The Science Book is the result of a 3-year effort to study the science capabilities of networks of next generation detectors. Such networks would make it possible to address unsolved problems in numerous areas of physics and astronomy, from Cosmology to Beyond the Standard Model of particle physics, and how they could provide insights into workings of strongly gravitating systems, astrophysics of compact objects and the nature of dense matter. It is inevitable that observatories of such depth and finesse will make new discoveries inaccessible to other windows of observation. In addition to laying out the rich science potential of the next generation of detectors, this report provides specific science targets in five different areas in physics and astronomy and the sensitivity requirements to accomplish those science goals. This report is the second in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book (this report), iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network., Comment: 69 pages, 18 figures

We perform 3D hydrodynamical simulations of a common-envelope event involving a 12 solar mass red supergiant donor. Massive stars are expected to be qualitatively different from low-mass stars as their envelopes have significant support from radiation pressure, which increases both the final separation and amount of mass ejected through the common-envelope interaction. We perform adiabatic simulations that include radiation energy through the equation of state, which results in ejecting 60 per cent more mass (up to two thirds of the total envelope mass becoming unbound, or more) and yield a 10 per cent larger final separation compared to simulations that assume an ideal gas. When also including recombination energy, we find that at least three quarters of the envelope, and possibly the entire envelope, may be unbound. The final separation further increases by almost 20 per cent. The additional amount of ejected material is mainly due to energy injected from helium recombination. Hydrogen recombination plays a comparatively small role, as it mainly occurs in gas that has already become unbound. We conclude that the internal energy of the envelope can be a significant energy source for ejecting the common envelope, but ultimately radiation transport and convection need to be included., Comment: Accepted for publication in MNRAS, 17 pages excluding appendices, 15 figures

GW200115 was the second merger of a black hole and a neutron star confidently detected through gravitational waves. Inference on the signal allows for a large black hole spin misaligned with the orbital angular momentum, but shows little support for aligned spin values. We show that this is a natural consequence of measuring the parameters of a black hole -- neutron star binary with non-spinning components while assuming the priors used in the LIGO-Virgo-KAGRA analysis. We suggest that, a priori, a non-spinning binary is more consistent with current astrophysical understanding., Comment: Minor updates to the version accepted to ApJ Letters

I present an analytical fit to the duration of the inspiral of two point masses driven by gravitational-wave emission. The fit is accurate to within 3% over the entire range of initial eccentricities from 0 to 0.99999.

Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It includes a number of tools for population processing in addition to the core binary evolution components. COMPAS is publicly available via the github repository https://github.com/TeamCOMPAS/COMPAS/, and is designed to allow for flexible modifications as evolutionary models improve. This paper describes the methodology and implementation of COMPAS. It is a living document which will be updated as new features are added to COMPAS; the current document describes COMPAS v02.21.00., Comment: Code publicly available via https://compas.science . Minor updates to match version accepted to ApJS

'Publish or perish' is an expression describing the pressure on academics to consistently publish research to ensure a successful career in academia. With a global pandemic that has changed the world, how has it changed academic productivity? Here we show that academics are posting just as many publications on the arXiv pre-print server as if there were no pandemic: 168,630 were posted in 2020, a +12.6% change from 2019 and $+1.4\sigma$ deviation above the predicted 162,577 $\pm$ 4,393. However, some immediate impacts are visible in individual research fields. Conference cancellations have led to sharp drops in pre-prints, but laboratory closures have had mixed effects. Only some experimental fields show mild declines in outputs, with most being consistent on previous years or even increasing above model expectations. The most significant change is a 50% increase ($+8\sigma$) in quantitative biology research, all related to the COVID-19 pandemic. Some of these publications are by biologists using arXiv for the first time, and some are written by researchers from other fields (e.g., physicists, mathematicians). While quantitative biology pre-prints have returned to pre-pandemic levels, 20% of the research in this field is now focussed on the COVID-19 pandemic, demonstrating a strong shift in research focus., Comment: Submitted to RSOS

Recent work paints a conflicting portrait of the distribution of black hole spins in merging binaries measured with gravitational waves. Some analyses find that a significant fraction of merging binaries contain at least one black hole with a spin tilt $>90^\circ$ with respect to the orbital angular momentum vector, which has been interpreted as a signature for dynamical assembly. Other analyses find the data are consistent with a bimodal population in which some binaries contain black holes with negligible spin while the rest contain black holes with spin vectors preferentially aligned with the orbital angular momentum vector. In this work, we scrutinize models for the distribution of black hole spins to pinpoint possible failure modes in which the model yields a faulty conclusion. We reanalyze data from the second LIGO--Virgo gravitational-wave transient catalog (GWTC-2) using a revised spin model, which allows for a sub-population of black holes with negligible spins. In agreement with recent results by Roulet et al., we show that the GWTC-2 detections are consistent with two distinct sub-populations. We estimate that $29-75\%$ (90\% credible interval) of merging binaries contain black holes with negligible spin $\chi \approx 0$. The remaining binaries are part of a second sub-population in which the spin vectors are preferentially (but not exactly) aligned to the orbital angular momentum. The black holes in this second sub-population are characterized by spins of $\chi\sim0.45$. We suggest that the inferred spin distribution is consistent with the hypothesis that all merging binaries form via the field formation scenario., Comment: 13 pages, 6 figures

We explore the effect of anisotropic wind driving on the properties of accretion onto black holes in close binaries. We specifically focus on line-driven winds, which are common in high-mass X-ray binaries. In close binary systems, the tidal force from the companion star can modify the wind structure in two different ways. One is the reduction of wind terminal velocity due to the weaker effective surface gravity. The other is the reduction in mass flux due to gravity darkening. We incorporate these effects into the so-called CAK theory in a simple way and investigate the wind flow around the accretor on the orbital scale. We find that a focused accretion stream is naturally formed when the Roche lobe filling factor is $\gtrsim0.8$-0.9, analogous to that of wind Roche lobe overflow, but only when the velocity reduction is taken into account. The formation of a stream is necessary to bring in sufficient angular momentum to form an accretion disc around the black hole. Gravity darkening effects reduce the amount of accreted angular momentum, but not enough to prevent the formation of a disc. Based on these results, we expect there to be a discrete step in the observability of high-mass X-ray binaries depending on whether the donor Roche lobe filling factor is below or above $\sim$0.8-0.9., Comment: accepted for publication in the Publications of the Astronomical Society of Australia

Binary neutron stars have been observed as millisecond pulsars, gravitational-wave sources, and as the progenitors of short gamma-ray bursts and kilonovae. Massive stellar binaries that evolve into merging double neutron stars are believed to experience a common-envelope episode. During this episode, the envelope of a giant star engulfs the whole binary. The energy transferred from the orbit to the envelope by drag forces or from other energy sources can eject the envelope from the binary system, leading to a stripped short-period binary. In this paper, we use one-dimensional single stellar evolution to explore the final stages of the common-envelope phase in progenitors of neutron star binaries. We consider an instantaneously stripped donor star as a proxy for the common-envelope phase and study the star's subsequent radial evolution. We determine a range of stripping boundaries which allow the star to avoid significant rapid re-expansion and which thus represent plausible boundaries for the termination of the common-envelope episode. We find that these boundaries lie above the maximum compression point, a commonly used location of the core/envelope boundary. We conclude that stars may retain fractions of a solar mass of hydrogen-rich material even after the common-envelope episode. If we consider orbital energy as the only energy source available, all of our models would overfill their Roche lobe after ejecting the envelope, whose binding energy includes gravitational, thermal, radiation, and recombination energy terms., Comment: 14 pages, 10 figures, 3 tables. Accepted for publication in MNRAS

Gravitational-wave detections are enabling measurements of the rate of coalescences of binaries composed of two compact objects -- neutron stars and/or black holes. The coalescence rate of binaries containing neutron stars is further constrained by electromagnetic observations, including Galactic radio binary pulsars and short gamma-ray bursts. Meanwhile, increasingly sophisticated models of compact objects merging through a variety of evolutionary channels produce a range of theoretically predicted rates. Rapid improvements in instrument sensitivity, along with plans for new and improved surveys, make this an opportune time to summarise the existing observational and theoretical knowledge of compact-binary coalescence rates., Comment: Invited review article for Living Reviews in Relativity. Updates relative to published version based on recent observations and models. The authors very much welcome further suggestions. All code and data are publicly available via https://zenodo.org/record/7017532 (version 7)

Interacting binaries are of general interest as laboratories for investigating the physics of accretion, which gives rise to the bulk of high-energy radiation in the Galaxy. They allow us to probe stellar evolution processes that cannot be studied in single stars. Understanding the orbital evolution of binaries is essential in order to model the formation of compact binaries. Here we focus our attention on studying orbital evolution driven by angular momentum loss through stellar winds in massive binaries. We run a suite of hydrodynamical simulations of binary stars hosting one mass losing star with varying wind velocity, mass ratio, wind velocity profile and adiabatic index, and compare our results to analytic estimates for drag and angular momentum loss. We find that, at leading order, orbital evolution is determined by the wind velocity and the binary mass ratio. Small ratios of wind to orbital velocities and large accreting companion masses result in high angular momentum loss and a shrinking of the orbit. For wider binaries and binaries hosting lighter mass-capturing companions, the wind mass-loss becomes more symmetric, which results in a widening of the orbit. We present a simple analytic formula that can accurately account for angular momentum losses and changes in the orbit, which depends on the wind velocity and mass ratio. As an example of our formalism, we compare the effects of tides and winds in driving the orbital evolution of high mass X-ray binaries, focusing on Vela X-1 and Cygnus X-1 as examples., Comment: 22 pages, 16 figures, submitted to ApJ, comments welcome

Observations of binary pulsars and pulsars in globular clusters suggest that at least some pulsars must receive weak natal kicks at birth. If all pulsars received strong natal kicks above \unit[50]{\kms}, those born in globular clusters would predominantly escape, while wide binaries would be disrupted. On the other hand, observations of transverse velocities of isolated radio pulsars indicate that only $5\pm2\%$ have velocities below \unit[50]{\kms}. We explore this apparent tension with rapid binary population synthesis modelling. We propose a model in which supernovae with characteristically low natal kicks (e.g., electron-capture supernovae) only occur if the progenitor star has been stripped via binary interaction with a companion. We show that this model naturally reproduces the observed pulsar speed distribution and without reducing the predicted merging double neutron star yield. We estimate that the zero-age main sequence mass range for non-interacting progenitors of electron-capture supernovae should be no wider than ${\approx}0.2 M_\odot$., Comment: 7 pages including figures, tables, and references. 3 figures including 5 pdfs, and 1 table

In this paper, a fast and parallelizable method based on Gaussian Processes (GPs) is introduced to emulate computer models that simulate the formation of binary black holes (BBHs) through the evolution of pairs of massive stars. Two obstacles that arise in this application are the a priori unknown conditions of BBH formation and the large scale of the simulation data. We address them by proposing a local emulator which combines a GP classifier and a GP regression model. The resulting emulator can also be utilized in planning future computer simulations through a proposed criterion for sequential design. By propagating uncertainties of simulation input through the emulator, we are able to obtain the distribution of BBH properties under the distribution of physical parameters., Comment: 24 pages, 11 figures

Active galactic nuclei (AGN) are prominent environments for stellar capture, growth and formation. These environments may catalyze stellar mergers and explosive transients, such as thermonuclear and core-collapse supernovae (SNe). SN explosions in AGN discs generate strong shocks, leading to unique observable signatures. We develop an analytical model which follows the evolution of the shock propagating in the disc until it eventually breaks out. We derive the peak luminosity, bolometric lightcurve, and breakout time. The peak luminosities may exceed $10^{45}$ erg s$^{-1}$ and last from hours to days. The brightest explosions occur in regions of reduced density; either off-plane, or in discs around low-mass central black holes ($\sim 10^6\ M_\odot$), or in starved subluminous AGNs. Explosions in the latter two sites are easier to observe due to a reduced AGN background luminosity. We perform suites of 1D Lagrangian radiative hydrodynamics SNEC code simulations to validate our results and obtain the luminosity in different bands, and 2D axisymmetric Eulerian hydrodynamics code HORMONE simulations to study the morphology of the ejecta and its deviation from spherical symmetry. The observed signature is expected to be a bright blue, UV, or X-ray flare on top of the AGN luminosity from the initial shock breakout, while the subsequent red part of the lightcurve will largely be unobservable. We estimate the upper limit for the total event rate to be $\mathcal{R}\lesssim 100\ \rm yr^{-1}\ Gpc^{-3}$ for optimal conditions and discuss the large uncertainties in this estimate. Future high-cadence transient searches may reveal these events. Some existing tidal disruption event candidates may originate from AGN supernovae., Comment: Accepted to MNRAS

We investigate the shape of the jet break in within-beam gamma-ray burst (GRB) optical afterglows for various lateral jet structure profiles. We consider cases with and without lateral spreading and a range of inclinations within the jet core half-opening angle, $\theta_c$. We fit model and observed afterglow lightcurves with a smoothly-broken power-law function with a free-parameter $\kappa$ that describes the sharpness of the break. We find that the jet break is sharper ($\kappa$ is greater) when lateral spreading is included than in the absence of lateral spreading. For profiles with a sharp-edged core, the sharpness parameter has a broad range of $0.1\lesssim\kappa\lesssim4.6$, whereas profiles with a smooth-edged core have a narrower range of $0.1\lesssim\kappa\lesssim2.2$ when models both with and without lateral spreading are included. For sharp-edged jets, the jet break sharpness depends strongly on the inclination of the system within $\theta_c$, whereas for smooth-edged jets, $\kappa$ is more strongly dependent on the size of $\theta_c$. Using a sample of 20 GRBs we find nine candidate smooth-edged jet structures and eight candidate sharp-edged jet structures, while the remaining three are consistent with either. The shape of the jet break, as measured by the sharpness parameter $\kappa$, can be used as an initial check for the presence of lateral structure in within-beam GRBs where the afterglow is well-sampled at and around the jet-break time., Comment: 13 pages, 5 figures, accepted for publication in MNRAS 15/07/21