Search for: All records

Abstract Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ∼10–100 Hz, over a characteristic jet activity timescale ∼10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ∼10 events yr −1 . The detection rate drops to ∼1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events. 

The engulfment of substellar bodies (SBs), such as brown dwarfs and planets, by giant stars is a possible explanation for rapidly rotating giants, lithium-rich giants, and the presence of SBs in close orbits around subdwarfs and white dwarfs. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an engulfed SB. We model the SB as a rigid body with a reflective surface because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram-pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory inside the star. We find that engulfment can increase the luminosity of a 1M_⊙star by up to a few orders of magnitude. The time for the star to return to its original luminosity is up to a few thousand years when the star has evolved to ≈10R_⊙and up to a few decades at the tip of the red giant branch (RGB). No SBs can eject the envelope of a 1M_⊙star before it evolves to ≈10R_⊙if the orbit of the SB is the only energy source contributing to the ejection. In contrast, SBs as small as ≈10M_Jupcan eject the envelope at the tip of the RGB. The numerical framework we introduce here can be used to study planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB.

 

The Milky Way is believed to host hundreds of millions of quiescent stellar-mass black holes (BHs). In the last decade, some of these objects have been potentially uncovered via gravitational microlensing events. All these detections resulted in a degeneracy between the velocity and the mass of the lens. This degeneracy has been lifted, for the first time, with the recent astrometric microlensing detection of OB110462. However, two independent studies reported very different lens masses for this event. Sahu et al. inferred a lens mass of 7.1 ± 1.3M_⊙, consistent with a BH, while Lam et al. inferred 1.6–4.2M_⊙, consistent with either a neutron star or a BH. Here, we study the landscape of isolated BHs formed in the field. In particular, we focus on the mass and center-of-mass speed of four subpopulations: isolated BHs from single-star origin, disrupted BHs of binary-star origin, main-sequence stars with a compact object companion, and double compact object mergers. Our model predicts that most (≳70%) isolated BHs in the Milky Way are of binary origin. However, noninteractions lead to most massive BHs (≳15–20M_⊙) being predominantly of single origin. Under the assumption that OB110462 is a free-floating compact object, we conclude that it is more likely to be a BH originally belonging to a binary system. Our results suggest that low-mass BH microlensing events can be useful to understand binary evolution of massive stars in the Milky Way, while high-mass BH lenses can be useful to probe single stellar evolution.

 

Metal-poor stars in the Milky Way (MW) halo display large star-to-star dispersion in theirr-process abundance relative to lighter elements. This suggests a chemically diverse and unmixed interstellar medium (ISM) in the early universe. This study aims to help shed light on the impact of turbulent mixing, driven by core-collapse supernovae (cc-SNe), on ther-process abundance dispersal in galactic disks. To this end, we conduct a series of simulations of small-scale galaxy patches which resolve metal-mixing mechanisms at parsec scales. Our setup includes cc-SNe feedback and enrichment fromr-process sources. We find that the relative rate of ther-process events to cc-SNe is directly imprinted on the shape of ther-process distribution in the ISM with more frequent events causing more centrally peaked distributions. We consider also the fraction of metals that is lost on galactic winds and find that cc-SNe are able to efficiently launch highly enriched winds, especially in smaller galaxy models. This result suggests that smaller systems, e.g., dwarf galaxies, may require higher levels of enrichment in order to achieve similar meanr-process abundances as MW-like progenitors systems. Finally, we are able to place novel constraints on the production rate ofr-process elements in the MW,$6\times {10}^{-7}{M}_{\odot }{\mathrm{yr}}^{-1}\lesssim {\dot{m}}_{\mathrm{rp}}\ll 4.7\times {10}^{-4}{M}_{\odot }{\mathrm{yr}}^{-1}$, imposed by accurately reproducing the mean and dispersion of [Eu/Fe] in metal-poor stars. Our results are consistent with independent estimates from alternate methods and constitute a significant reduction in the permitted parameter space.

 

Abstract Stars grazing supermassive black holes (SMBHs) on bound orbits may survive tidal disruption, causing periodic flares. Inspired by the recent discovery of the periodic nuclear transient ASASSN-14ko, a promising candidate for a repeating tidal disruption event (TDE), we study the tidal deformation of stars approaching SMBHs on eccentric orbits. With both analytical and hydrodynamic methods, we show the overall tidal deformation of a star is similar to that in a parabolic orbit provided that the eccentricity is above a critical value. This allows one to make use of existing simulation libraries from parabolic encounters to calculate the mass fallback rate in eccentric TDEs. We find the flare structures of eccentric TDEs show a complicated dependence on both the SMBH mass and the orbital period. For stars orbiting SMBHs with relatively short periods, we predict significantly shorter-lived duration flares than those in parabolic TDEs, which can be used to predict repeating events if the mass of the SMBH can be independently measured. Using an adiabatic mass-loss model, we study the flare evolution over multiple passages, and show the evolved stars can survive many more passages than main-sequence stars. We apply this theoretical framework to the repeating TDE candidate ASASSN-14ko and suggest that its recurrent flares originate from a moderately massive ( M ≳ 1 M ⊙ ), extended (likely ≈10 R ⊙ ), evolved star on a grazing, bound orbit around the SMBH. Future hydrodynamic simulations of multiple tidal interactions will enable realistic models on the individual flare structure and the evolution over multiple flares. 

Stellar-mass black holes can become embedded within the disks of active galactic nuclei (AGNs). Afterwards, their interactions are mediated by their gaseous surroundings. Here, we study the evolution of stellar-mass binary black holes (BBHs) embedded within AGN disks using three-dimensional hydrodynamic simulations and analytic methods, focusing on environments where the AGN disk scale heightHis ≳ the BBH sphere of influence. We model the local surroundings of the embedded BBHs using a wind tunnel formalism and characterize different accretion regimes based on the local properties of the disk. We develop prescriptions for accretion and drag for embedded BBHs. Using these prescriptions with AGN disk models that can represent the Toomre-unstable outer regions of AGN disks, we study the long-term evolution of BBHs as they migrate through the disk. We find that BBHs typically merge within ≲1–30 Myr, increasing their mass significantly in the process, allowing BBHs to enter (or cross) the pair-instability supernova mass gap. The BBH accretion rate often exceeds the Eddington limit, sometimes by several orders of magnitude. Many embedded BBHs will merge before migrating significantly in the disk. We also discuss possible electromagnetic signatures during and following the inspiral, finding that it is generally unlikely for the bolometric luminosity of the BBH to exceed the AGN luminosity.

 

Common-envelope evolution is a stage in binary system evolution in which a giant star engulfs a companion. The standard energy formalism is an analytical framework to estimate the amount of energy transferred from the companion's shrinking orbit into the envelope of the star that engulfed it. We show analytically that this energy transfer is larger than predicted by the standard formalism. As the orbit of the companion shrinks, the mass it encloses becomes smaller, and the companion is less bound than if the enclosed mass had remained constant. Therefore, more energy must be transferred to the envelope for the orbit to shrink further. We derive a revised energy formalism that accounts for this effect, and discuss its consequences in two contexts: the formation of neutron star binaries, and the engulfment of planets and brown dwarfs by their host stars. The companion mass required to eject the stellar envelope is smaller by up to 50% , leading to differences in common-envelope evolution outcomes. The energy deposition in the outer envelope of the star, which is related to the transient luminosity and duration, is up to a factor of ≈7 higher. Common-envelope efficiency values above unity, as defined in the literature, are thus not necessarily unphysical, and result at least partly from an incomplete description of the energy deposition. The revised energy formalism presented here can improve our understanding of stellar merger and common-envelope observations and simulations. 

ABSTRACT TIC 470710327, a massive compact hierarchical triple-star system, was recently identified by NASA’s Transiting Exoplanet Survey Satellite. TIC 470710327 is comprised of a compact (1.10 d) circular eclipsing binary, with total mass $\approx 10.9\!-\!13.2\, \rm {M_{\odot }}$, and a more massive $\approx 14\!-\!17\, \rm {M_{\odot }}$ eccentric non-eclipsing tertiary in a 52.04 d orbit. Here, we present a progenitor scenario for TIC 470710327 in which ‘2 + 2’ quadruple dynamics result in Zeipel–Lidov–Kozai oscillations that lead to a contact phase of the more massive binary. In this scenario, the two binary systems should form in a very similar manner, and dynamics trigger the merger of the more massive binary either during late phases of star formation or several Myr after the zero-age main sequence, when the stars begin to expand. Any evidence that the tertiary is a highly magnetized (∼1–10 kG), slowly rotating blue main-sequence star would hint towards a quadruple origin. Finally, our scenario suggests that the population of inclined compact multiple-stellar systems is reduced into coplanar systems, via mergers, late during star formation or early in the main sequence. The elucidation of the origin of TIC 470710327 is crucial in our understanding of multiple massive star formation and evolution. 

The ∼100 tidal disruption events (TDEs) observed so far exhibit a wide range of emission properties both at peak and over their lifetimes. Some TDEs radiate predominantly at X-ray energies, while others radiate chiefly at UV and optical wavelengths. While the peak luminosities across TDEs show distinct properties, the evolutionary behavior can also vary between TDEs with similar peak emission properties. In particular, for optical TDEs, while their UV and optical emissions decline somewhat following the fallback pattern, some events can greatly rebrighten in X-rays at late time. In this Letter, we conduct three-dimensional general relativistic radiation magnetohydrodynamics simulations of TDE accretion disks at varying accretion rates in the regime of super-Eddington accretion. We make use of Monte Carlo radiative transfer simulations to calculate the reprocessed spectra at various inclinations and at different evolutionary stages. We confirm the unified model proposed by Dai et al., which predicts that the observed emission largely depends on the viewing angle of the observer with respect to the disk orientation. Furthermore, we find that disks with higher accretion rates have elevated wind and disk densities, which increases the reprocessing of the high-energy radiation and thus generally augments the optical-to-X-ray flux ratio along a particular viewing angle. This implies that at later times, as the accretion level declines, we expect that more X-rays will leak out along intermediate viewing angles. Such dynamical model for TDEs can provide a natural explanation for the diversity in the emission properties observed in TDEs at peak and along their temporal evolution.

 

The extent to which turbulence mixes gas in the face of recurrent infusions of fresh metals by supernovae (SN) could help provide important constraints on the local star formation conditions. This includes predictions of the metallicity dispersion among metal-poor stars, which suggests that the interstellar medium was not very well mixed at these early times. The purpose of thisLetteris to help isolate, via a series of numerical experiments, some of the key processes that regulate turbulent mixing of SN elements in galactic disks. We study the gas interactions in small simulated patches of a galaxy disk with the goal of resolving the small-scale mixing effects of metals at parsec scales, which enables us to measure the turbulent diffusion coefficient in various galaxy environments. By investigating the statistics of variations ofαelements in these simulations, we are able to derive constraints not only on the allowed range of intrinsic yield variations in SN explosions but also on the star formation history of the Milky Way. We argue that the observed dispersion of [Mg/Fe] in metal-poor halo stars is compatible with the star-forming conditions expected in dwarf satellites or in an early low-star-forming Milky Way progenitor. In particular, metal variations in stars that have not been phase-mixed can be used to infer the star-forming conditions of disrupted dwarf satellites.