Researcher profile

Ilya Mandel

Ilya Mandel contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2025arXiv

Converged simulations of the nozzle shock in tidal disruption events

When debris from a star that experienced a tidal disruption events (TDE) after passing too close to a massive black hole returns to pericenter on the second passage, it is compressed, leading to the formation of nozzle shocks (in the orbital plane) and pancake shocks (perpendicular to the orbital plane). Resolving these shocks is a long-standing problem in the hydrodynamic simulations of parabolic TDEs. Excessive numerical energy dissipation or heating unrealistically expands the stream. In this Letter, we apply adaptive particle refinement to our 3D general relativistic smoothed particle simulations to locally increase the resolution near the pericenter. We achieve resolutions equivalent to $6.55\times10^{11}$ particles, allowing us to converge on the true energy dissipation. We conclude that only $4\times10^{-5}$ of the orbital energy is dissipated in nozzle shocks for a Sun-like star tidally disrupted by a $10^6$ solar-mass black hole, therefore the nozzle shocks are unlikely to be important in the evolution of TDEs.

preprint2022arXiv

Building better spin models for merging binary black holes: Evidence for non-spinning and rapidly spinning nearly aligned sub-populations

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 $χ\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 $χ\sim0.45$. We suggest that the inferred spin distribution is consistent with the hypothesis that all merging binaries form via the field formation scenario.

preprint2022arXiv

Common envelopes in massive stars II: The distinct roles of hydrogen and helium recombination

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.

preprint2022arXiv

Common envelopes in massive stars: Towards the role of radiation pressure and recombination energy in ejecting red supergiant envelopes

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.

preprint2022arXiv

Extreme eccentricities of triple systems: Analytic results

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 $η$) and the double-averaging approximation (parametrised by the period ratio, $ε_{\rm SA}$) are relaxed, for circular outer orbits. We recover known results in both limiting cases (either $η$ or $ε_{\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.

preprint2022arXiv

Merging stellar-mass binary black holes

The LIGO and Virgo detectors have directly observed gravitational waves from mergers of pairs of stellar-mass black holes, along with a smaller number of mergers involving neutron stars. These observations raise the hope that compact object mergers could be used as a probe of stellar and binary evolution, and perhaps of stellar dynamics. This colloquium-style article summarises the existing observations, describes theoretical predictions for formation channels of merging stellar-mass black-hole binaries along with their rates and observable properties, and presents some prospects for gravitational-wave astronomy.

preprint2022arXiv

Rates of Compact Object Coalescences

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.

preprint2022arXiv

Stellar response after stripping as a model for common-envelope outcomes

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.

preprint2021arXiv

Cygnus X-1 contains a 21-solar mass black hole -- implications for massive star winds

The evolution of massive stars is influenced by the mass lost to stellar winds over their lifetimes. These winds limit the masses of the stellar remnants (such as black holes) that the stars ultimately produce. We use radio astrometry to refine the distance to the black hole X-ray binary Cygnus X-1, which we find to be $2.22^{+0.18}_{-0.17}$ kiloparsecs. When combined with previous optical data, this implies a black hole mass of $21.2\pm2.2$ solar masses, higher than previous measurements. The formation of such a high-mass black hole in a high-metallicity system constrains wind mass loss from massive stars.

preprint2021arXiv

Impact of Massive Binary Star and Cosmic Evolution on Gravitational Wave Observations II: Double Compact Object Rates and Properties

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.

preprint2021arXiv

Rapid stellar and binary population synthesis with COMPAS

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.

preprint2021arXiv

Wind mass-loss rates of stripped stars inferred from Cygnus X-1

Recent observations of the high-mass X-ray binary Cygnus X-1 have shown that both the companion star (41 solar masses) and the black hole (21 solar masses) are more massive than previously estimated. Furthermore, the black hole appears to be nearly maximally spinning. Here we present a possible formation channel for the Cygnus X-1 system that matches the observed system properties. In this formation channel, we find that the orbital parameters of Cygnus X-1, combined with the observed metallicity of the companion, imply a significant reduction in mass loss through winds relative to commonly used prescriptions for stripped stars.

preprint2020arXiv

Age dissection of the Milky Way discs: red giants in the Kepler field

[Abridged] Ensemble studies of red-giant stars with exquisite asteroseismic, spectroscopic, and astrometric constraints offer a novel opportunity to recast and address long-standing questions concerning the evolution of stars and of the Galaxy. Here, we infer masses and ages for nearly 5400 giants with available Kepler light curves and APOGEE spectra, and discuss some of the systematics that may affect the accuracy of the inferred stellar properties. First, we look at age-chemical-abundances relations. We find a dearth of young, metal-rich stars, and the existence of a significant population of old (8-9 Gyr), low-[$α$/Fe], super-solar metallicity stars, reminiscent of the age and metallicity of the well-studied open cluster NGC6791. The age-chemo-kinematic properties of these stars indicate that efficient radial migration happens in the thin disk. We find that ages and masses of the nearly 400 $α$-element-rich red-giant-branch (RGB) stars in our sample are compatible with those of an old (~11 Gyr), nearly coeval, chemical-thick disk population. Using a statistical model, we show that 95% of the population was born within ~1.5 Gyr. Moreover, we find a difference in the vertical velocity dispersion between low- and high-[$α$/Fe] populations, confirming their different chemo-dynamical histories. We then exploit the almost coeval $α$-rich population to gain insight into processes that may have altered the mass of a star along its evolution, which are key to improve the mapping of the observed stellar mass to age. We find evidence for a mean integrated RGB mass loss <$Δ$M>= 0.10 $\pm$ 0.02 Msun and that the occurrence of massive (M $\gtrsim$ 1.1 Msun) $α$-rich stars is of the order of 5% on the RGB, and significantly higher in the RC, supporting the scenario in which most of these stars had undergone interaction with a companion.

preprint2020arXiv

An alternative interpretation of GW190412 as a binary black hole merger with a rapidly spinning secondary

The LIGO-Virgo collaboration recently reported the properties of GW190412, a binary black hole merger with unequal component masses (mass ratio $0.25^{+0.06}_{-0.04}$ when using the EOBNR PHM approximant) and a non-vanishing effective spin aligned with the orbital angular momentum. They used uninformative priors to infer that the more massive black hole had a dimensionless spin magnitude between 0.17 and 0.59 at 90\% confidence. We argue that, within the context of isolated binary evolution, it is more natural to assume a priori that the first-born, more massive black hole has a negligible spin, while the spin of the less massive black hole is preferentially aligned with the orbital angular momentum if it is spun up by tides. Under this astrophysically motivated prior, we conclude that the lower mass black hole had a dimensionless spin component between 0.64 and 0.99 along the orbital angular momentum.

preprint2020arXiv

Be X-ray binaries in the SMC as indicators of mass transfer efficiency

Be X-ray binaries (BeXRBs) consist of rapidly rotating Be stars with neutron star companions accreting from the circumstellar emission disk. We compare the observed population of BeXRBs in the Small Magellanic Cloud with simulated populations of BeXRB-like systems produced with the COMPAS population synthesis code. We focus on the apparently higher minimal mass of Be stars in BeXRBs than in the Be population at large. Assuming that BeXRBs experienced only dynamically stable mass transfer, their mass distribution suggests that at least 30% of the mass donated by the progenitor of the neutron star is typically accreted by the B-star companion. We expect these results to affect predictions for the population of double compact object mergers. A convolution of the simulated BeXRB population with the star formation history of the Small Magellanic Cloud shows that the excess of BeXRBs is most likely explained by this galaxy&#39;s burst of star formation around 20--40 Myr ago.

preprint2020arXiv

Constraining the masses of microlensing black holes and the mass gap with Gaia DR2

Context: Gravitational microlensing is sensitive to compact-object lenses in the Milky Way, including white dwarfs, neutron stars or black holes, and could potentially probe a wide range of stellar remnant masses. However, the mass of the lens can be determined only in very limited cases, due to missing information on both source and lens distances and their proper motions. Aims: We aim at improving the mass estimates in the annual parallax microlensing events found in the 8 years of OGLE-III observations towards the Galactic Bulge (Wyrzykowski et al. 2016) with the use of Gaia Data Release 2 (DR2). Methods: We use Gaia DR2 data on distances and proper motions of non-blended sources and recompute the masses of lenses in parallax events. We also identify new events in that sample which are likely to have dark lens; the total number of such events is now 18. Results: The derived distribution of masses of dark lenses is consistent with a continuous distribution of stellar remnant masses. A mass gap between neutron-star and black-hole masses in the range between 2 and 5 solar masses is not favoured by our data, unless black holes receive natal-kicks above 20-80 km/s. We present 8 candidates for objects with masses within the putative mass gap, including a spectacular multi-peak parallax event with mass of $2.4^{+1.9}_{-1.3}\ M_\odot$ located just at 600 pc. The absence of an observational mass gap between neutron stars and black holes, or, conversely, the evidence for black hole natal kicks if a mass gap is assumed, can inform future supernova modelling efforts.

preprint2020arXiv

Detecting Double Neutron Stars with LISA

We estimate the properties of the double neutron star (DNS) population that will be observable by the planned space-based interferometer LISA. By following the gravitational radiation driven evolution of DNSs generated from rapid population synthesis of massive binary stars, we estimate that around 35 DNSs will accumulate a signal-to-noise ratio above 8 over a four-year LISA mission. The observed population mainly comprises Galactic DNSs (94 per cent), but detections in the LMC (5 per cent) and SMC (1 per cent) may also be expected. The median orbital frequency of detected DNSs is expected to be 0.8 mHz, and many of them will be eccentric (median eccentricity of $0.11$). The orbital properties will provide insights into DNS progenitors and formation channels. LISA is expected to localise these DNSs to a typical angular resolution of $2^\circ$, with best-constrained sources localised to a few arcminutes. These localisations may allow neutron star natal kick magnitudes to be constrained through the Galactic distribution of DNSs, and make it possible to follow up the sources with radio pulsar searches. However, LISA is also expected to resolve $\sim 10^4$ Galactic double white dwarfs, many of which may have binary parameters that resemble DNSs; we discuss how the combined measurement of binary eccentricity, chirp mass, and sky location may aid the identification of a DNS. We expect the best-constrained DNSs to have eccentricities known to a few parts in a thousand, chirp masses measured to better than 1 per cent fractional uncertainty, and sky localisation at the level of a few arcminutes.

preprint2020arXiv

Luminous Red Novae: population models and future prospects

A class of optical transients known as Luminous Red Novae (LRNe) have recently been associated with mass ejections from binary stars undergoing common-envelope evolution. We use the population synthesis code COMPAS to explore the impact of a range of assumptions about the physics of common-envelope evolution on the properties of LRNe. In particular, we investigate the influence of various models for the energetics of LRNe on the expected event rate and light curve characteristics, and compare with the existing sample. We find that the Galactic rate of LRNe is $\sim 0.2$ yr$^{-1}$, in agreement with the observed rate. In our models, the luminosity function of Galactic LRNe covers multiple decades in luminosity and is dominated by signals from stellar mergers, consistent with observational constraints from iPTF and the Galactic sample of LRNe. We discuss how observations of the brightest LRNe may provide indirect evidence for the existence of massive ($> 40$ M$_\odot$) red supergiants. Such LRNe could be markers along the evolutionary pathway leading to the formation of double compact objects. We make predictions for the population of LRNe observable in future transient surveys with the Large Synoptic Survey Telescope and the Zwicky Transient Facility. In all plausible circumstances, we predict a selection-limited observable population dominated by bright, long-duration events caused by common envelope ejections. We show that the Large Synoptic Survey Telescope will observe $20$--$750$ LRNe per year, quickly constraining the luminosity function of LRNe and probing the physics of common-envelope events.

preprint2020arXiv

Merger rates in primordial black hole clusters without initial binaries

Primordial black holes formed through the collapse of cosmological density fluctuations have been hypothesised as contributors to the dark matter content of the Universe. At the same time, their mergers could contribute to the recently observed population of gravitational-wave sources. We investigate the scenario in which primordial black holes form binaries at late times in the Universe. Specifically, we re-examine the mergers of primordial black holes in small clusters of ~30 objects in the absence of initial binaries. Binaries form dynamically through Newtonian gravitational interactions. These binaries act as heat sources for the cluster, increasing the cluster&#39;s velocity dispersion, which inhibits direct mergers through gravitational-wave two-body captures. Meanwhile, three-body encounters of tight binaries are too rare to tighten binaries sufficiently to allow them to merge through gravitational-wave emission. We conclude that in the absence of initial binaries, merger rates of primordial black holes in the Bird et al. (2016) initial cluster configuration are at least an order of magnitude lower than previously suggested, which makes gravitational-wave detections of such sources improbable.

preprint2020arXiv

The origin of spin in binary black holes: Predicting the distributions of the main observables of Advanced LIGO

We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of the main observables of Advanced LIGO. We used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. By assuming efficient angular momentum transport, we are able to present a model that is capable of simultaneously predicting the three main gravitational-wave observables: the effective inspiral spin parameter $χ_{eff}$, the chirp mass $M_{chirp}$, and the cosmological redshift of merger $z_{merger}$. We find an excellent agreement between our model and the ten events from the first two advanced detector observing runs. We make predictions for the third observing run O3 and for Advanced LIGO design sensitivity. We expect around 80% of events with $χ_{eff} < 0.1$, while the remaining 20% of events with $χ_{eff} \ge 0.1$ are split into ~10% with $M_{chirp} < 15$ M$_\odot$ and ~10% with $M_{chirp} \ge 15$ M$_\odot$. In conclusion, the favorable comparison of the existing LIGO/Virgo observations with our model predictions gives support to the idea that the majority, if not all of the observed mergers, originate from the evolution of isolated binaries. The first-born black hole has negligible spin because it lost its envelope after it expanded to become a giant star, while the spin of the second-born black hole is determined by the tidal spin up of its naked helium star progenitor by the first-born black hole companion after the binary finished the common-envelope phase.

preprint2020arXiv

The spectral evolution of AT 2018dyb and the presence of metal lines in tidal disruption events

We present light curves and spectra of the tidal disruption event (TDE) ASASSN-18pg / AT 2018dyb spanning a period of one year. The event shows a plethora of strong emission lines, including the Balmer series, He II, He I and metal lines of O III $λ$3760 and N III $λλ$ 4100, 4640 (blended with He II). The latter lines are consistent with originating from the Bowen fluorescence mechanism. By analyzing literature spectra of past events, we conclude that these lines are common in TDEs. The spectral diversity of optical TDEs is thus larger than previously thought and includes N-rich events besides H- and He-rich events. We study how the spectral lines evolve with time, by means of their width, relative strength, and velocity offsets. The velocity width of the lines starts at $\sim$ 13000 km s$^{-1}$ and decreases with time. The ratio of He II to N III increases with time. The same is true for ASASSN-14li, which has a very similar spectrum to AT 2018dyb but its lines are narrower by a factor of $>$2. We estimate a black hole mass of $M_{\rm BH}$ = $3.3^{+5.0}_{-2.0}\times 10^6$ $M_{\odot}$ by using the $M$-$σ$ relation. This is consistent with the black hole mass derived using the MOSFiT transient fitting code. The detection of strong Bowen lines in the optical spectrum is an indirect proof for extreme ultraviolet and (reprocessed) X-ray radiation and favors an accretion origin for the TDE optical luminosity. A model where photons escape after multiple scatterings through a super-Eddington thick disk and its optically thick wind, viewed at an angle close to the disk plane, is consistent with the observations.

preprint2019arXiv

Astrophysical science metrics for next-generation gravitational-wave detectors

The second generation of gravitational-wave detectors are being built and tuned all over the world. The detection of signals from binary black holes is beginning to fulfill the promise of gravitational-wave astronomy. In this work, we examine several possible configurations for third-generation laser interferometers in existing km-scale facilities. We propose a set of astrophysically motivated metrics to evaluate detector performance. We measure the impact of detector design choices against these metrics, providing a quantitative cost-benefit analyses of the resulting scientific payoffs.

preprint2019arXiv

Common Envelope Wind Tunnel: The Effects of Binary Mass Ratio and Implications for the Accretion-Driven Growth of LIGO Binary Black Holes

We present three-dimensional local hydrodynamic simulations of flows around objects embedded within stellar envelopes using a &#34;wind tunnel&#34; formalism. Our simulations model the common envelope dynamical inspiral phase in binary star systems in terms of dimensionless flow characteristics. We present suites of simulations that study the effects of varying the binary mass ratio, stellar structure, equation of state, relative Mach number of the object&#39;s motion through the gas, and density gradients across the gravitational focusing scale. For each model, we measure coefficients of accretion and drag experienced by the embedded object. These coefficients regulate the coupled evolution of the object&#39;s masses and orbital tightening during the dynamical inspiral phase of the common envelope. We extrapolate our simulation results to accreting black holes with masses comparable to that of the population of LIGO black holes. We demonstrate that the mass and spin accrued by these black holes per unit orbital tightening are directly related to the ratio of accretion to drag coefficients. We thus infer that the mass and dimensionless spin of initially non-rotating black holes change by of order $1\%$ and 0.05, respectively, in a typical example scenario. Our prediction that the masses and spins of black holes remain largely unmodified by a common envelope phase aids in the interpretation of the properties of the growing observed population of merging binary black holes. Even if these black holes passed through a common envelope phase during their assembly, features of mass and spin imparted by previous evolutionary epochs should be preserved.

preprint2019arXiv

Eclipses of continuous gravitational waves as a probe of stellar structure

Although gravitational waves only interact weakly with matter, their propagation is affected by a gravitational potential. If a gravitational wave source is eclipsed by a star, measuring these perturbations provides a way to directly measure the distribution of mass throughout the stellar interior. We compute the expected Shapiro time delay, amplification and deflection during an eclipse, and show how this can be used to infer the mass distribution of the eclipsing body. We identify continuous gravitational waves from neutron stars as the best candidates to detect this effect. When the Sun eclipses a far-away source, depending on the depth of the eclipse the time-delay can change by up to $\sim 0.034$ ms, the gravitational-wave strain amplitude can increase by $\sim 4$%, and the apparent position of the source in the sky can vary by $4&#39;&#39;$. Accreting neutron stars with Roche-lobe filling companion stars have a high probability of exhibiting eclipses, producing similar time delays but undetectable changes in amplitude and sky location. Even for the most rapidly rotating neutron stars, this time delay only corresponds to a few percent of the phase of the gravitational wave, making it an extremely challenging measurement. However, if sources of continuous gravitational waves exist just below the limit of detection of current observatories, next-generation instruments will be able to observe them with enough precision to measure the signal of an eclipsing star. Detecting this effect would provide a new direct probe to the interior of stars, complementing asteroseismology and the detection of solar neutrinos.

preprint2019arXiv

Explosions Driven by the Coalescence of a Compact Object with the Core of a Massive-Star Companion Inside a Common Envelope: Circumstellar Properties, Light Curves, and Population Statistics

We model explosions driven by the coalescence of a black hole or neutron star with the core of its massive-star companion. Upon entering a common envelope phase, a compact object may spiral all the way to the core. The concurrent release of energy is likely to be deposited into the surrounding common envelope, powering a merger-driven explosion. We use hydrodynamic models of binary coalescence to model the common envelope density distribution at the time of coalescence. We find toroidal profiles of material, concentrated in the binary&#39;s equatorial plane and extending to many times the massive star&#39;s original radius. We use the spherically-averaged properties of this circumstellar material (CSM) to estimate the emergent light curves that result from the interaction between the blast wave and the CSM. We find that typical merger-driven explosions are brightened by up to three magnitudes by CSM interaction. From population synthesis models we discover that the brightest merger-driven explosions, $M_V \sim -18$ to $-19$, are those involving black holes because they have the most massive and extended CSM. Black hole coalescence events are also common; they represent about 50% of all merger-driven explosions and approximately 0.3% of the core-collapse rate. Merger-driven explosions offer a window into the highly-uncertain physics of common envelope interactions in binary systems by probing the properties of systems that merge rather than eject their envelopes.

preprint2018arXiv

Comment on &#34;An excess of massive stars in the local 30 Doradus starburst&#34;

Schneider et al. (Science, 2018) used an ad hoc statistical method in their calculation of the stellar initial mass function. Adopting an improved approach, we reanalyse their data and determine a power law exponent of $2.05_{-0.14}^{+0.13}$. Alternative assumptions regarding data set completeness and the star formation history model can shift the inferred exponent to $2.11_{-0.19}^{+0.17}$ and $2.15_{-0.13}^{+0.13}$, respectively.