Search

The interstellar medium (ISM) of disk galaxies is turbulent, and yet the fundamental nature of ISM turbulence, the energy cascade, is not understood in detail. In this study, we use high-resolution simulations of a hydrodynamical, gravitationally stratified, supernova (SNe)-driven, multiphase ISM to probe the nature of a galactic turbulence cascade. Through the use of kinetic energy flux transfer functions split into interactions between compressible $\mathbf{u}_c$ and incompressible $\mathbf{u}_s$ modes, we show that there exists a large-to-small-scale cascade in both $\mathbf{u}_c$ and $\mathbf{u}_s$ when mediated by an additional $\mathbf{u}_s$ mode. But the $\mathbf{u}_s$ cascade is highly non-local. Moreover, there is a $\mathbf{u}_c$ mediated component of the $\mathbf{u}_s$ cascade that proceeds in the opposite direction -- an inverse cascade from small-to-large scales. The cascade feeds flux into scales well beyond the scale height, energizing the winds and fueling the direct cascades. Both the strongly non-local and the inverse $\mathbf{u}_s$ cascades happen on scales that have a power law $\mathbf{u}_s$ energy spectrum, highlighting how degenerate the spectrum is to the true underlying physical processes. We directly show that the inverse cascade comes from $\mathbf{u}_s$ modes interacting with expanding SNe remnants (SNRs) and that $\mathbf{u}_s$ modes are generated to leading order via baroclinic, highly corrugated cooling layers between warm $(T\lesssim 10^4\,\rm{K})$ and hot $(T\gg10^4\,\rm{K})$ gas in these SNRs. Finally, we outline a complete phenomenology for SNe-driven turbulence in a galactic disk, estimate a $10^{-16}\,\rm{G}$ Biermann field generated from SNR cooling layers, and highlight the strong deviations that SNe-driven turbulence has from the conventional Kolmogorov model., Comment: Main text 24 pages, 13 figures. Appendix, 6 pages, 7 figures. Submitted to ApJ. Comments welcome

We present a novel method that enables us to estimate the acceleration of individual millisecond pulsars (MSPs) using only their spin period and its time derivative. For our binary MSP sample, we show that one can obtain an empirical calibration of the magnetic braking term that relies only on observed quantities. We find that such a model for magnetic braking is only valid for MSPs with small surface magnetic field strengths ($<3\times10^8$ G) and large characteristic ages ($>$ 5 Gyr). With this method we are able to effectively double the number of pulsars with line-of-sight acceleration measurements, from 28 to 54 sources. This expanded dataset leads to an updated measurement of the total density in the midplane, which we find to be $\rho_0$ = 0.086 $\pm$ 0.001 stat. $\pm$ 0.006 sys. M$_\odot$/pc$^3$, and an updated measurement of the local dark matter density, which we calculate to be $\rho_{0,\mathrm{DM}}$ = 0.014 $\pm$ 0.005 stat. $\pm$ 0.006 sys. M$_\odot$/pc$^3$ (0.53 $\pm$ 0.30 GeV/cm$^3$). This updated value for $\rho_{0,\mathrm{DM}}$ is in good agreement with literature values derived from kinematic estimates. We show that each new acceleration measurement improves the precision of $\rho_{0,\mathrm{DM}}$ by the same amount as roughly 10$^5$ stars. The pulsar accelerations are very asymmetric above and below the disk; we show that the shape and size of this asymmetry can be largely explained by the north-south asymmetry of disk star counts and the offset in the Milky Way disk and halo centers of mass due to the Large Magellanic Cloud., Comment: Comments welcome

We present a detailed analysis of AT 2020nov, a tidal disruption event (TDE) in the center of its host galaxy, located at a redshift of $z = 0.083$. AT 2020nov exhibits unique features, including double-peaked Balmer emission lines, a broad UV/optical flare, and a peak log luminosity in the extreme ultra-violet (EUV) estimated at $\sim$$45.66^{+0.10}_{-0.33} \; \mathrm{erg} \, \mathrm{s^{-1}}$. A late-time X-ray flare was also observed, reaching an absorbed luminosity of $1.67 \times 10^{43} \; \mathrm{erg} \, \mathrm{s^{-1}}$ approximately 300 days after the UV/optical peak. Multi-wavelength coverage, spanning optical, UV, X-ray, and mid-infrared (MIR) bands, reveals a complex spectral energy distribution (SED) that includes MIR flaring indicative of dust echoes, suggesting a dust covering fraction consistent with typical TDEs. Spectral modeling indicates the presence of an extended, quiescent disk around the central supermassive black hole (SMBH) with a radius of $\sim$$5.06^{+0.59}_{-0.77} \times 10^4 \; \mathrm{R_g}$. The multi-component SED model, which includes a significant EUV component, suggests that the primary emission from the TDE is reprocessed by this extended disk, producing the observed optical and MIR features. The lack of strong AGN signatures in the host galaxy, combined with the quiescent disk structure, highlights AT 2020nov as a rare example of a TDE occurring in a galaxy with a dormant but extended pre-existing accretion structure., Comment: 39 pages, 20 figures, submitted to ApJ

We expand the theoretical framework by Gottlieb el al. (2023), which connects binary merger populations with long and short binary gamma-ray bursts (lbGRBs and sbGRBs), incorporating kilonovae as a key diagnostic tool. We show that lbGRBs, powered by massive accretion disks around black holes (BHs), should be accompanied by bright, red kilonovae. In contrast, sbGRBs - if also powered by BHs - would produce fainter, red kilonovae, potentially biasing against their detection. However, magnetized hypermassive neutron star (HMNS) remnants that precede BH formation can produce jets with power ($P_{\rm NS} \approx 10^{51}\,{\rm erg\,s^{-1}}$) and Lorentz factor ($\Gamma>10$), likely compatible with sbGRB observations, and would result in distinctly bluer kilonovae, offering a pathway to identifying the sbGRB central engine. Recent modeling by Rastinejad et al. (2024) found luminous red kilonovae consistently accompany lbGRBs, supporting lbGRB originating from BH-massive disk systems, likely following a short-lived HMNS phase. The preferential association of sbGRBs with comparably luminous kilonovae argues against the BH engine hypothesis for sbGRBs, while the bluer hue of these KNe provides additional support for an HMNS-driven mechanism. Within this framework, BH-NS mergers likely contribute exclusively to the lbGRB population with red kilonovae. Our findings suggest that GW170817 may, in fact, have been an lbGRB to on-axis observers. Finally, we discuss major challenges faced by alternative lbGRB progenitor models, such as white dwarf-NS or white dwarf-BH mergers and accretion-induced collapse forming magnetars, which fail to align with observed GRB timescales, energies, and kilonova properties.

The energy provided in the radioactive decay of thorium (Th) and uranium (U) isotopes, embedded in planetary mantles, sustains geodynamics important for surface habitability such as the generation of a planetary magnetic dynamo. In order to better understand the thermal evolution of nearby exoplanets, stellar photospheric abundances can be used to infer the material composition of orbiting planets. Here we constrain the intrinsic dispersion of the r-process element europium (Eu) (measured in relative abundance [Eu/H]) as a proxy for Th and U in local F, G, and K type dwarf stars. Adopting stellar-chemical data from two high quality spectroscopic surveys, we have determined a small intrinsic scatter of 0.025 dex in [Eu/H] within the disk. We further investigate the stellar anti-correlation in [Eu/$\alpha$] vs [$\alpha$/H] at late metallicities to probe in what regimes planetary radiogenic heating may lead to periods of extended dynamo collapse. We find that only near-solar metallicity stars in the disk have Eu inventories supportive of a persistent dynamo in attendant planets, supporting the notion of a ``metallicity Goldilocks zone'' in the galactic disk. The observed anti-correlation further provides novel evidence regarding the nature of r-processes injection by substantiating $\alpha$ element production is decoupled from Eu injection. This suggests either a metallicity-dependent r-process in massive core-collapse supernovae, or that neutron-star merger events dominate r-process production in the recent universe., Comment: 11 pages, 4 figures

We propose a new predictive theory for the analysis of common envelope (CE) events which incorporates the effects of relevant hydrodynamical processes into a simple analytical framework. We introduce the ejection and dynamical parameters $\xi$ and $\beta$, which define if envelope ejection is energetically or hydrodynamically favorable, respectively, during CE inspiral. When combined, these parameters offer a detailed narrative of how inspiral begins, proceeds, and ends that is consistent with preliminary comparisons to 3D hydrodynamical models. This physically-motivated framework impacts predictions for CE outcomes, especially for systems that have energy excess, and offers promise as a potential alternative for the treatment of CE in binary population synthesis., Comment: 12 pages, 3 figures. Accepted for publication in ApJ Letters

The Earth sits inside a 300pc-wide void that was carved by a series of supernova explosions that went off tens of millions of years ago, pushing away interstellar gas and creating a bubble-like structure. The $^{60}$Fe peak deposits found in the deep-sea crust have been interpreted by the imprints left by the ejecta of supernova explosions occurring about 2-3 and 5-6 Myr ago. It is likely that the $^{60}$Fe peak at about 2-3 Myr originated from a supernova occurring in the Upper Centaurus Lupus association in Scorpius Centaurus ($\approx$140 pc) or the Tucana Horologium association ($\approx$70 pc). Whereas, the $\approx$ 5-6 Myr peak is likely attributed to the solar system's entrance into the bubble. In this {\it Letter}, we show that the supernova source responsible for synthesizing the $^{60}$Fe peak deposits $\approx$ 2-3 Myr ago can consistently explain the cosmic-ray spectrum and the large-scale anisotropy between 100 TeV and 100 PeV. The cosmic-ray knee could then potentially be attributed entirely to a single nearby "Pevatron" source. Matching the intensity and shape of the cosmic-ray spectrum allows us to place stringent constraints on the cosmic-ray energy content from the supernova as well as on the cosmic-ray diffusion coefficient. Making use of such constraints we provide a robust estimate of the temporal variation of terrestrial ionizing cosmic radiation levels and discuss their implications in the development of early life on Earth by plausibly influencing the mutation rate and, as such, conceivably assisting in the evolution of complex organisms., Comment: 10 pages, 5 figures, submitted to the Astrophysical Journal Letters

Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4\sigma$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in \textit{Gaia} DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6\%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 80 kyr, a companion mass of 0.3 $M_\odot$ (indicative of a white dwarf or low-mass star), and a semi-major axis of 2000 AU. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios.

Observations of tidal disruption events (TDEs) show signs of Nitrogen enrichment reminiscent of other astrophysical sources such as active galactic nuclei (AGN) and star-forming galaxies. Given that TDEs probe the gas from a single star, it is possible to test if the observed enrichment is consistent with expectations from the CNO cycle by looking at the observed Nitrogen/Carbon (N/C) abundance ratios. Given that $\approx 20\%$ of solar mass stars (and an even larger fraction of more massive stars) live in close binaries, it is worthwhile to also consider what TDEs from stars influenced by binary evolution would look like. We show here that TDEs from stars stripped of their Hydrogen-rich (and Nitrogen-poor) envelopes through previous binary-induced mass loss can produce much higher observable N/C enhancements than even TDEs from massive stars. Additionally, we predict that the time-dependence of the N/C abundance ratio in the mass fallback rate of stripped stars will follow the inverse behavior of main-sequence stars, enabling a more accurate characterization of the disrupted star., Comment: 11 pages, 5 figures, submitted to ApJL

Stars grazing supermassive black holes on bound orbits may produce periodic flares over many passages, known as repeating partial tidal disruption events (TDEs). Here, we present 3D hydrodynamic simulations of sun-like stars over multiple tidal encounters. The star is significantly restructured and becomes less concentrated as a result of mass loss and tidal heating. The vulnerability to mass loss depends sensitively on the stellar density structure, and the strong correlation between the fractional mass loss $\Delta M/M_*$ and the ratio of the central and average density $\rho_{\mathrm{c}}/\bar\rho$, which was initially derived in disruption simulations of main-sequence stars, also applies for stars strongly reshaped by tides. Over multiple orbits, the star loses progressively more mass in each encounter and is doomed to a complete disruption. Throughout its lifetime, the star may produce numerous weak flares (depending on the initial impact parameter), followed by a couple of luminous flares whose brightness increases exponentially. Flux-limited surveys are heavily biased toward the brightest flares, which may appear similar to the flare produced by the same star undergoing a full disruption on its first tidal encounter. This places new challenges on constraining the intrinsic TDE rates, which need to take repeating TDEs into account. Other types of stars with different initial density structures (e.g., evolved stars with massive cores) follow distinct evolution tracks, which might explain the diversity of the long-term luminosity evolution seen in recently uncovered repeaters., Comment: 15 pages, 5 figures, 1 table. Updated to accepted version (ApJ)

We present detailed radio observations of the tidal disruption event (TDE) ASASSN-19bt/AT2019ahk, obtained with the Australia Telescope Compact Array (ATCA), the Atacama Large Millimeter/submillimeter Array (ALMA), and the MeerKAT radio telescopes, spanning 40 to 1464 days after the onset of the optical flare. We find that ASASSN-19bt displays unusual radio evolution compared to other TDEs, as the peak brightness of its radio emission increases rapidly until 457 days post-optical discovery and then plateaus. Using a generalized approach to standard equipartition techniques, we estimate the energy and corresponding physical parameters for two possible emission geometries: a non-relativistic spherical outflow and a relativistic outflow observed from an arbitrary viewing angle. We find that the non-relativistic solution implies a continuous energy rise in the outflow from $E\sim10^{46}$ erg to $E\sim10^{49}$ erg with $\beta \approx 0.05$, while the off-axis relativistic jet solution instead suggests $E\approx10^{52}$ erg with $\Gamma\sim10$ erg at late times in the maximally off-axis case. We find that neither model provides a holistic explanation for the origin and evolution of the radio emission, emphasizing the need for more complex models. ASASSN-19bt joins the population of TDEs that display unusual radio emission at late times. Conducting long-term radio observations of these TDEs, especially during the later phases, will be crucial for understanding how these types of radio emission in TDEs are produced., Comment: 25 pages. Submitted to ApJ

Triple stellar systems allow us to study stellar processes that cannot be attained in binary stars. The evolutionary phases in which the stellar members undergo mass exchanges can alter the hierarchical layout of these systems. Yet, the lack of a self-consistent treatment of common-envelope (CE) in triple star-systems hinders the comprehensive understanding of their long-term fate. This letter examines the conditions predicted around binaries embedded within CEs using local 3D hydrodynamical simulations. We explore varying the initial binary separation, the flow Mach number, and the background stellar density gradients as informed by a wide array of CE conditions, including those invoked to explain the formation of the triple system hosting PSR J0337+1715. We find that the stellar density gradient governs the gaseous drag force, which determines the final configuration of the embedded binary. We observe a comparable net drag force on the center of mass but an overall reduction in the accretion rate of the binary compared to the single object case. We find that for most CE conditions, and in contrast to the uniform background density case, the binary orbital separation increases with time, softening the binary and preventing it from subsequently merging. We conclude that binaries spiraling within CEs become more vulnerable to be disrupted by tidal interactions. This can have profound implications on the final outcomes of triple star-systems.

The vision of 2030STEM is to address systemic barriers in institutional structures and funding mechanisms required to achieve full inclusion in Science, Technology, Engineering, and Mathematics (STEM) and accelerate leadership pathways for individuals from underrepresented populations across STEM sectors. 2030STEM takes a systems-level approach to create a community of practice that can test, learn and promote programs and policies that affirm and value cultural identities in STEM. To achieve parity and full representation in the STEM workforce, a variety of changes are needed across academia and STEM professional industries (e.g., business, finance, biotech, government) to accelerate underrepresented groups into positions of leadership throughout the STEM ecosystem. Through a series of subject matter interviews, roundtables, and curated analysis four major themes have surfaced, which, if implemented, could exponentially accelerate the creation of critical pathways to leadership, break down pre-existing barriers and biases, intentionally elevate the voices, value, and research of underrepresented groups in STEM, and implement new structural strategies at scale. This white paper provides a summary of innovative new practices designed to accelerate inclusion, including expanding on known global toolkits, new funding strategies, and the structural changes required throughout various STEM professions to propel pathways to leadership for underrepresented groups., Comment: 13 pages

We present a start-to-end simulation aimed at studying the long-term fate of high-mass X-ray binaries and whether a Thorne-\.Zytkow object (T\.ZO) might ultimately be assembled. We analyze results from a 3D hydrodynamical simulation that models the eventual fate of LMC X-4, a compact high-mass X-ray binary system, after the primary fills its Roche lobe and engulfs the neutron star companion. We discuss the outcome of this engulfment within the standard paradigm of T\.ZO formation. The post-merger angular momentum content of the stellar core is a key ingredient, as even a small amount of rotation can break spherical symmetry and produce a centrifugally supported accretion disk. Our findings suggest the inspiraling neutron star, upon merging with the core, can accrete efficiently via a disk at high rates ($\approx 10^{-2}M_\odot/{\rm s}$), subsequently collapsing into a black hole and triggering a bright transient with a luminosity and duration typical of an ultra-long gamma-ray burst. We propose that the canonical framework for T\.ZO formation via common envelope needs to be revised, as the significant post-merger accretion feedback will unavoidably unbind the vast majority of the surrounding envelope.

Thorne-\.{Z}ytkow objects (T\.{Z}Os), hypothetical merger products in which a neutron star is embedded in a stellar core, are traditionally considered steady-state configurations. Their assembly, especially through dynamical channels, is not well-understood. The predominant focus in the literature has been on the observational signatures related to the evolution and long-term fate of T\.{Z}Os, with their initial formation often treated as a given. However, the foundational calculations supporting the existence of T\.{Z}Os assume non-rotating spherically-symmetric initial conditions that we find to be inconsistent with a binary merger scenario. In this work, we explore the implications of post-merger dynamics in T\.{Z}O formation scenarios with field binary progenitors, specifically the role that angular momentum transport during the common envelope phase plays in constraining possible merger products, using the tools of stellar evolution and three-dimensional hydrodynamics. We also propose an alternative steady-state outcome for these mergers: the thin-envelope T\.{Z}O, an equilibrium solution consisting of a low-mass spherical envelope supported by the accretion disk luminosity of a central stellar-mass black hole. These configurations may be of interest to upcoming time-domain surveys as potential X-ray sources that may be preceded by a series of bright transient events., Comment: 14 pages, 8 figures, Accepted to ApJ

We present analytical and numerical models of the bright long GRB 210822A at $z=1.736$. The intrinsic extreme brightness exhibited in the optical, which is very similar to other bright GRBs (e.g., GRBs 080319B, 130427A, 160625A 190114C, and 221009A), makes GRB 210822A an ideal case for studying the evolution of this particular kind of GRB. We use optical data from the RATIR instrument starting at $T+315.9$ s, with publicly available optical data from other ground-based observatories, as well as Swift/UVOT, and X-ray data from the Swift/XRT instrument. The temporal profiles and spectral properties during the late stages align consistently with the conventional forward shock model, complemented by a reverse shock element that dominates optical emissions during the initial phases ($T<300$ s). Furthermore, we observe a break at $T=80000$s that we interpreted as evidence of a jet break, which constrains the opening angle to be about $\theta_\mathrm{j}=(3-5)$ degrees. Finally, we apply a machine-learning technique to model the multi-wavelength light curve of GRB 210822A using the AFTERGLOWPY library. We estimate the angle of sight $\theta_{obs}=(6.4 \pm 0.1) \times 10^{-1}$ degrees, the energy $E_0=(7.9 \pm 1.6)\times 10^{53}$ ergs, the electron index $p=2.54 \pm 0.10$, the thermal energy fraction in electrons $\epsilon_\mathrm{e}=(4.63 \pm 0.91) \times 10^{-5}$ and in the magnetic field $\epsilon_\mathrm{B}= (8.66 \pm 1.01) \times 10^{-6}$, the efficiency $\chi = 0.89 \pm 0.01$, and the density of the surrounding medium $n_\mathrm{0} = 0.85 \pm 0.01 cm^{-3}$., Comment: Resubmitted to MNRAS after moderate revision, 12 pages, 6 figures

We present the photometric and spectroscopic evolution of SN 2022oqm, a nearby multi-peaked hydrogen- and helium-weak calcium-rich transient (CaRT). SN 2022oqm was detected 13.1 kpc from its host galaxy, the face-on spiral galaxy NGC 5875. Extensive spectroscopic coverage reveals an early hot (T >= 40,000 K) continuum and carbon features observed $\sim$1~day after discovery, SN Ic-like photospheric-phase spectra, and strong forbidden calcium emission starting 38 days after discovery. SN 2022oqm has a relatively high peak luminosity (MB = -17 mag) for (CaRTs), making it an outlier in the population. We determine that three power sources are necessary to explain the light curve (LC), with each corresponding to a distinct peak. The first peak is powered by an expanding blackbody with a power law luminosity, suggesting shock cooling by circumstellar material (CSM). Subsequent LC evolution is powered by a double radioactive decay model, consistent with two sources of photons diffusing through optically thick ejecta. From the LC, we derive an ejecta mass and 56Ni mass of ~0.6 solar masses and ~0.09 solar masses. Spectroscopic modeling suggests 0.6 solar masses of ejecta, and with well-mixed Fe-peak elements throughout. We discuss several physical origins for SN 2022oqm and find either a surprisingly massive white dwarf progenitor or a peculiar stripped envelope model could explain SN 2022oqm. A stripped envelope explosion inside a dense, hydrogen- and helium-poor CSM, akin to SNe Icn, but with a large 56Ni mass and small CSM mass could explain SN 2022oqm. Alternatively, helium detonation on an unexpectedly massive white dwarf could also explain SN 2022oqm., Comment: 35 pages, 17 figures, 7 tables, Accepted for Publication in ApJ

Collapsing massive stars lead to a broad range of astrophysical transients, whose multi-wavelength emission is powered by a variety of processes including radioactive decay, activity of the central engine, and interaction of the outflows with a dense circumstellar medium. These transients are also candidate factories of neutrinos with energy up to hundreds of PeV. We review the energy released by such astrophysical objects across the electromagnetic wavebands as well as neutrinos, in order to outline a strategy to optimize multi-messenger follow-up programs. We find that, while a significant fraction of the explosion energy can be emitted in the infrared-optical-ultraviolet (UVOIR) band, the optical signal alone is not optimal for neutrino searches. Rather, the neutrino emission is strongly correlated with the one in the radio band, if a dense circumstellar medium surrounds the transient, and with X-rays tracking the activity of the central engine. Joint observations of transients in radio, X-rays, and neutrinos will crucially complement those in the UVOIR band, breaking degeneracies in the transient parameter space. Our findings call for heightened surveys in the radio and X-ray bands to warrant multi-messenger detections., Comment: 26 pages, including 7 figures and 2 appendices. Minor revisions, conclusions unchanged. Matches version accepted for publication in PRD

We study the properties of galaxies hosting mid-infrared outbursts in the context of a catalog of five hundred thousand galaxies from the Sloan Digital Sky Survey. We find that nuclear obscuration, as inferred by the surrounding dust mass, does not correlate with host galaxy type, stellar properties (e.g. total mass and mean age), or with the extinction of the host galaxy as estimated by the Balmer decrement. This implies that nuclear obscuration may not be able to explain any over-representation of tidal disruption events in particular host galaxies. We identify a region in the galaxy catalog parameter space that contains all unobscured tidal disruption events but only harbors $\lesssim $ 11\% of the mid-infrared outburst hosts. We find that mid-infrared outburst hosts appear more centrally concentrated and have higher galaxy S\'ersic indices than galaxies hosting active galactic nuclei (AGN) selected using the BPT classification. We thus conclude that the majority of mid-infrared outbursts are not hidden tidal disruption events but are instead consistent with being obscured AGN that are highly variable, such as changing-look AGN., Comment: 17 pages, 9 figures, accepted in ApJ Letters

We present optical, ultraviolet, and infrared data of the type II supernova (SN II) 2020jfo at 14.5 Mpc. This wealth of multiwavelength data allows to compare different metrics commonly used to estimate progenitor masses of SN II for the same object. Using its early light curve, we infer SN 2020jfo had a progenitor radius of $\approx$700 $R_{\odot}$, consistent with red supergiants of initial mass $M_{\rm ZAMS}=$11-13 $M_{\odot}$. The decline in its late-time light curve is best fit by a ${}^{56}$Ni mass of 0.018$\pm$0.007 $M_{\odot}$ consistent with that ejected from SN II-P with $\approx$13 $M_{\odot}$ initial mass stars. Early spectra and photometry do not exhibit signs of interaction with circumstellar matter, implying that SN 2020jfo experienced weak mass loss within the final years prior to explosion. Our spectra at $>$250 days are best fit by models from 12 $M_{\odot}$ initial mass stars. We analyzed integral field unit spectroscopy of the stellar population near SN 2020jfo, finding its massive star population had a zero age main sequence mass of 9.7$\substack{+2.5\\-1.3} M_{\odot}$. We identify a single counterpart in pre-explosion imaging and find it has an initial mass of at most $7.2\substack{+1.2\\-0.6} M_{\odot}$. We conclude that the inconsistency between this mass and indirect mass indicators from SN 2020jfo itself is most likely caused by extinction with $A_{V}=2$-3 mag due to matter around the progenitor star, which lowered its observed optical luminosity. As SN 2020jfo did not exhibit extinction at this level or evidence for interaction with circumstellar matter between 1.6-450 days from explosion, we conclude that this material was likely confined within $\approx$3000 $R_{\odot}$ from the progenitor star., Comment: 31 pages, 12 figures, accepted to MNRAS

Dynamical perturbations from supermassive black hole (SMBH) binaries can increase the rates of tidal disruption events (TDEs). However, most previous work focuses on TDEs from the heavy black hole in the SMBH binary (SMBHB) system. In this work, we focus on the lighter black holes in SMBHB systems and show that they can experience a similarly dramatic increase in their TDE rate due to perturbations from a more massive companion. While the increase in TDEs around the more massive black hole is mostly due to chaotic orbital perturbations, we find that, around the smaller black hole, the eccentric Kozai-Lidov (EKL) mechanism is dominant and capable of producing a comparably large number of TDEs. In this instance, the mass derived from the light curve and spectra of TDEs triggered by the lighter SMBH companion are expected to be significant smaller than the SMBH mass estimated from galaxy scaling relations, which are dominated by the more massive companion. This apparent inconsistency can help find SMBHB candidates that are not currently accreting as active galactic nuclei (AGN) and that are at separations too small to be resolved as two distinct sources. TDEs thus provide us with an exciting opportunity to study SMBHB, in particular when they might be detected from galaxies hosting SMBHs that are too massive to disrupt sun-like stars., Comment: 14 pages, 7 figures, 1 table, submitted to ApJ, corrected ref. links

Tidal disruption events (TDEs) take place when a star ventures too close to a supermassive black hole (SMBH) and becomes ruptured. One of the leading proposed physical mechanisms often invoked in the literature involves weak two-body interactions experienced by the population of stars within the host SMBH's sphere of influence, commonly referred to as two-body relaxation. This process can alter the angular momentum of stars at large distances and place them into nearly radial orbits, thus driving them to disruption. On the other hand, gravitational perturbations from an SMBH companion via the eccentric Kozai-Lidov (EKL) mechanism have also been proposed as a promising stellar disruption channel. Here we demonstrate that the combination of EKL and two-body relaxation in SMBH binaries is imperative for building a comprehensive picture of the rates of TDEs. Here we explore how the density profile of the surrounding stellar distribution and the binary orbital parameters of an SMBH companion influence the rate of TDEs. We show that this combined channel naturally produces disruptions at a rate that is consistent with observations and also naturally forms repeated TDEs, where a bound star is partially disrupted on multiple orbits. Recent observations show stars being disrupted in short-period orbits, which is challenging to explain when these mechanisms are considered independently. However, the diffusive effect of two-body relaxation, combined with the secular nature of the eccentricity excitations from EKL, is found to drive stars on short eccentric orbits at a much higher rate., Comment: Fixed affiliation

Collapsing stars constitute the main black hole (BH) formation channel, and are occasionally associated with the launch of relativistic jets that power $ \gamma $-ray bursts (GRBs). Thus, collapsars offer an opportunity to infer the natal (before spin-up/down by accretion) BH spin directly from observations. We show that once the BH saturates with large-scale magnetic flux, the jet power is dictated by the BH spin and mass accretion rate. Core-collapse simulations by Halevi et al. 2023 and GRB observations favor stellar density profiles that yield an accretion rate $ \dot{m} \approx 10^{-2} M_\odot~{\rm s^{-1}} $, weakly dependent on time. This leaves the spin as the main factor that governs the jet power. By comparing the jet power to characteristic GRB luminosities, we find that the majority of BHs associated with jets are likely born slowly spinning with a dimensionless spin $ a \simeq 0.2 $, or $ a \lesssim 0.5 $ for wobbling jets, with the main uncertainty originating in the unknown $ \gamma $-ray radiative efficiency. This result could be applied to the entire core-collapse BH population, unless an anti-correlation between the stellar magnetic field and angular momentum is present. In a companion paper, Jacquemin-Ide et al. 2023, we show that regardless of the natal spin, the extraction of BH rotational energy leads to spin-down to $ a \lesssim 0.2 $, consistent with gravitational-wave observations. We verify our results by performing the first 3D general relativistic magnetohydrodynamic simulations of collapsar jets with characteristic GRB energies, powered by slowly spinning BHs. We find that jets of typical GRB power struggle to escape from the star, providing the first numerical indication that many jets fail to generate a GRB.

WR 20a is the most massive close-in binary known in our Galaxy. It is composed of two $\approx$80 M$_\odot$ Wolf-Rayet stars with a short period of $\approx$3.7 days in the open cluster Westerlund 2. As such, WR 20a presents us with a unique laboratory for studying the currently uncertain physics of binary evolution and compact object formation as well as for studying the wind collision region in an massive eclipsing binary system. We use deep Chandra observations of WR 20a to study the time variability of the wind collision region between the two Wolf-Rayet stars and are able to produce an X-ray light curve covering $\approx$2/3 of its orbital period. We find that the X-ray light curve is asymmetric because the flux of one peak is 2.5$\sigma$ larger than the flux of the other peak. This asymmetry could be caused by asymmetric mass-loss from the two stars or by the lopsidedness of the wind collision region due to the unusually fast rotation of the system. The X-ray light curve is also shifted in phase space when compared to the optical light curves measured by TESS and ASAS-SN. Additionally, we explore the ultimate fate of this system by modeling the resultant binary black hole merger expected at the end of the two stars' lives. We conclude that this system will evolve to be a representative of the sub-population of LIGO progenitors of fast-spinning binary black hole merger events., Comment: 13 pages, 5 figures, submitted to ApJ

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 $\approx7$ 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., Comment: We're currently addressing concerns raised by the community about the analysis presented in this paper. Once we address them, we will replace this submission accordingly. 9 pages, 4 figures

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 non-inspiral GW sources is also anticipated. We report the discovery of a new class of non-inspiral GW sources - the end states of massive stars - that can produce the brightest simulated stochastic GW burst signal in LVK bands known to date, and could be detectable in the 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, $\sim 10-100$ Hz, over a characteristic jet activity timescale, $\sim 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 estimated rate of $ \sim 10 $ events/year. The detection rate drops to $ \sim 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 multi-messenger events., Comment: For movies of the simulation, see https://oregottlieb.com/gw.html

Momentum feedback from isolated supernova remnants (SNRs) have been increasingly recognized by modern cosmological simulations as a resolution-independent means to implement the effects of feedback in galaxies, such as turbulence and winds. However, the integrated momentum yield from SNRs is uncertain due to the effects of SN clustering and interstellar medium (ISM) inhomogeneities. In this paper, we use spatially-resolved observations of the prominent 10-kpc star-forming ring of M31 to test models of mass-weighted ISM turbulence driven by momentum feedback from isolated, non-overlapping SNRs. We use a detailed stellar-age distribution (SAD) map from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, observationally-constrained SN delay-time distributions, and maps of the atomic and molecular hydrogen to estimate the mass-weighted velocity dispersion using the Martizzi et al. ISM turbulence model. Our estimates are within a factor of 2 of the observed mass-weighted velocity dispersion in most of the ring, but exceed observations at densities $\lesssim 0.2$ cm$^{-3}$ and SN rates $>2.1\times 10^{-4}$ SN yr$^{-1}$ kpc$^{-2}$, even after accounting for plausible variations in stellar-age distribution models and ISM scale height assumptions. We conclude that at high SN rates the momentum deposited is most likely suppressed by the non-linear effects of SN clustering, while at low densities, SNRs reach pressure equilibrium before the cooling phase. These corrections should be introduced in models of momentum-driven feedback and ISM turbulence., Comment: 14 pages, 7 figures, published in ApJ

A star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events (TDEs) have now been identified in the optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three TDEs. The continuum polarization is independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with optical photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code we generate predictions for the degree of polarization as a function of disk mass and viewing angle, and we show that the predicted levels are compatible with the observations, for extended reprocessing envelopes of $\sim$1000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for TDE models, and opens an important new line of exploration in the study of TDEs., Comment: Author's version of paper to appear in Nature Astronomy. In the journal version the detailed discussion on the ISP determination will be moved from the Methods section to a Supplementary Information section. 58 pages in double spacing format, including 5 Figures, 10 Extended Data Figures and 2 Tables

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 hydrodynamics 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\gtrsim 1 \mathrm{M_\odot}$), extended ($\approx$ a few $\mathrm{R_\odot}$), evolved star on a grazing, bound orbit around the SMBH. Future hydrodynamics simulations of multiple tidal interactions will enable realistic models on the individual flare structure and the evolution over multiple flares., Comment: 19 pages, 9 figures. Updated to accepted version (ApJ)

About a hundred tidal disruption events (TDEs) have been observed and they exhibit a wide range of emission properties both at peak and over their lifetimes. Some TDEs peak 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. At late time, some optical TDEs rebrighten in X-rays, while others maintain strong UV/optical emission components. In this Letter, we conduct three-dimensional general relativistic radiation magnetohydrodynamics simulations of TDE accretion disks at varying accretion rates ranging from a few to a few tens of the Eddington accretion rate. 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. (2018), which predicts that the observed emission largely depends on the viewing angle of the observer with respect to the disk orientation (X-ray strong when viewed face-on and UV/optically strong when viewed disk-on). What is more, 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 the funnel to effectively open up and allow more X-rays to 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., Comment: 1 table and 4 figures, submitted to ApJL

TIC 470710327, a massive compact hierarchical triple-star system, was recently identified by NASA's Transiting Exoplanet Survey Satellite (TESS). 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 (ZLK) resonances 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 (ZAMS), when the stars begin to expand. Any evidence that the tertiary is a highly-magnetised ($\sim 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 co-planar 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., Comment: 5 pages, 3 figures. Accepted for publication at MNRAS Letters

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 simulate the flow in the vicinity of an engulfed SB in three-dimensional hydrodynamics. 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_\odot$ 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 $\approx10R_\odot$ and up to a few decades at the tip of the red giant branch. No SBs can eject the envelope of a $1M_\odot$ star before it evolves to $\approx10R_\odot$, if the orbit of the SB is the only energy source contributing to the ejection. In contrast, SBs as small as $\approx10M_\text{Jup}$ can eject the envelope at the tip of the red giant branch. 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., Comment: 20 pages, 15 figures. Accepted for publication in ApJ. Comments welcome. Reproducibility repository at DOI 10.5281/zenodo.6368226

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 mass for this event. Sahu et al. (2022) inferred a lens mass of 7.1 $\pm$ 1.3 Msun, consistent with a BH, while Lam et al. (2022) inferred 1.6-4.2 Msun, 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 sub-populations: 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 ($\gtrsim$ 70%) isolated BHs in the Milky Way are of binary origin. However, non-interactions lead to most massive BHs ($\gtrsim$ 15-20 Msun) 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., Comment: 13 pages, 8 figures, 1 table. Accepted for publication in ApJL

Gravitational wave (GW) emissions from extreme-mass-ratio inspirals (EMRIs) are promising sources for low-frequency GW-detectors. They result from a compact object, such as a stellar-mass black-hole (BH), captured by a supermassive black hole (SMBH). Several physical processes have been proposed to form EMRIs. In particular, weak two-body interactions over a long time scale (i.e., relaxation processes) have been proposed as a likely mechanism to drive the BH orbit to high eccentricity. Consequently, it is captured by the SMBH and becomes an EMRI. Here we demonstrate that EMRIs are naturally formed in SMBH binaries. Gravitational perturbations from an SMBH companion, known as the eccentric Kozai-Lidov (EKL) mechanism, combined with relaxation processes, yield a significantly more enhanced rate than any of these processes operating alone. Since EKL is sensitive to the orbital configuration, two-body relaxation can alter the orbital parameters, rendering the system in a more EKL-favorable regime. As SMBH binaries are expected to be prevalent in the Universe, this process predicts a substantially high EMRI rate., Comment: 11 pages, 6 figures, accepted to ApJ-Letters

The X-ray source 3XMM~J150052.0+015452 was discovered as a spectacular tidal disruption event candidate during a prolonged ($>11$ yrs) outburst (Lin et al. 2017). It exhibited unique quasi-soft X-ray spectra of characteristic temperature $kT\sim0.3$ keV for several years at the peak, but in a recent Chandra observation (10 yrs into the outburst) a super-soft X-ray spectrum of $kT\sim0.15$ keV was detected. Such dramatic spectral softening could signal the transition from the super-Eddington to thermal state or the temporary presence of a warm absorber. Here we report on our study of four new XMM-Newton follow-up observations of the source. We found that they all showed super-soft spectra, suggesting that the source had remained super-soft for $>5$ yrs. Then its spectral change is best explained as due to the super-Eddington to thermal spectral state transition. The fits to the thermal state spectra suggested a smaller absorption toward the source than that obtained in Lin et al. (2017). This led us to update the modeling of the event as due to the disruption of a 0.75 msun star by a massive black hole of a few$\times10^5$ msun. We also obtained two HST images in the F606W and F814W filters and found that the dwarf star-forming host galaxy can be resolved into a dominant disk and a smaller bulge. No central point source was clearly seen in either filter, ruling out strong optical emission associated with the X-ray activity., Comment: 9 pages, 2 figures, Accepted for publication in ApJL

Seeing pristine material from the donor star in a Type Ia supernova (SN Ia) explosion can reveal the nature of the binary system. In this paper, we present photometric and spectroscopic observations of SN 2020esm, one of the best-studied SNe of the class of "super-Chandrasekhar" SNe Ia (SC SNe Ia), with data obtained $-12$ to +360 days relative to peak brightness, obtained from a variety of ground- and space-based telescopes. Initially misclassified as a Type II supernova, SN 2020esm peaked at $M_{B} = -19.9$ mag, declined slowly ($\Delta m_{15}(B) = 0.92$ mag), and had particularly blue UV and optical colors at early times. Photometrically and spectroscopically, SN 2020esm evolved similarly to other SC SNe Ia, showing the usual low ejecta velocities, weak intermediate mass elements (IMEs), and the enhanced fading at late times, but its early spectra are unique. Our first few spectra (corresponding to a phase of $\gtrsim$10~days before peak) reveal a nearly-pure carbon/oxygen atmosphere during the first days after explosion. This composition can only be produced by pristine material, relatively unaffected by nuclear burning. The lack of H and He may further indicate that SN 2020esm is the outcome of the merger of two carbon/oxygen white dwarfs (WDs). Modeling its bolometric light curve, we find a $^{56}$Ni mass of $1.23^{+0.14}_{-0.14}$ M$_{\odot}$ and an ejecta mass of $1.75^{+0.32}_{-0.20}$ M$_{\odot}$, in excess of the Chandrasekhar mass. Finally, we discuss possible progenitor systems and explosion mechanisms of SN 2020esm and, in general, the SC SNe Ia class., Comment: 21 pages, 15 figures, accepted for publication in APJ

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 amongst metal poor stars, which suggests that the interstellar medium was not very well mixed at these early times. The purpose of this {\it Letter} is 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 pc scales, which enables us to measure the turbulent diffusion coefficient in various galaxy environments. By investigating the statistics of variations of $\alpha$ 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., Comment: 11 pages, 5 figures. Accepted for publication in ApJL

We present observations of SN 2020fqv, a Virgo-cluster Type II core-collapse supernova (CCSN) with a high temporal resolution light curve from the Transiting Exoplanet Survey Satellite (TESS) covering the time of explosion; ultraviolet (UV) spectroscopy from the Hubble Space Telescope (HST) starting 3.3 days post-explosion; ground-based spectroscopic observations starting 1.1~days post-explosion; along with extensive photometric observations. Massive stars have complicated mass-loss histories leading up to their death as CCSNe, creating circumstellar medium (CSM) with which the SNe interact. Observations during the first few days post-explosion can provide important information about the mass-loss rate during the late stages of stellar evolution. Model fits to the quasi-bolometric light curve of SN 2020fqv reveal ~0.23 $M_{\odot}$ of CSM confined within ~1450 $R_{\odot}$ ($10^{14}$ cm) from its progenitor star. Early spectra (<4 days post-explosion), both from HST and ground-based observatories, show emission features from high-ionization metal species from the outer, optically thin part of this CSM. We find that the CSM is consistent with an eruption caused by the injection of $\sim$$5\times 10^{46}$ erg into the stellar envelope $\sim$300 days pre-explosion, potentially from a nuclear burning instability at the onset of oxygen burning. Light-curve fitting, nebular spectroscopy, and pre-explosion \textit{HST} imaging consistently point to a red supergiant (RSG) progenitor with $M_{\rm ZAMS}$$\approx$$13.5$--$15 \, M_{\odot}$, typical for SN~II progenitor stars. This finding demonstrates that a typical RSG, like the progenitor of SN 2020fqv, has a complicated mass-loss history immediately before core collapse., Comment: Published in MNRAS

We present deep X-ray and radio observations of the Fast Blue Optical Transient (FBOT) AT2020xnd/ZTF20acigmel at $z=0.2433$ from $13$d to $269$d after explosion. AT2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT2018cow. AT2020xnd shows luminous radio emission reaching $L_{\nu}\approx8\times10^{29}$ergs$^{-1}$Hz$^{-1}$ at 20GHz and $75$d post explosion, accompanied by luminous and rapidly fading soft X-ray emission peaking at $L_{X}\approx6\times10^{42}$ergs$^{-1}$. Interpreting the radio emission in the context of synchrotron radiation from the explosion's shock interaction with the environment we find that AT2020xnd launched a high-velocity outflow ($v\sim$0.1-0.2$c$) propagating into a dense circumstellar medium (effective $\dot M\approx10^{-3}M_{\rm{sol}}$yr$^{-1}$ for an assumed wind velocity of $v_w=1000$kms$^{-1}$). Similar to AT2018cow, the detected X-ray emission is in excess compared to the extrapolated synchrotron spectrum and constitutes a different emission component, possibly powered by accretion onto a newly formed black hole or neutron star. These properties make AT2020xnd a high-redshift analog to AT2018cow, and establish AT2020xnd as the fourth member of the class of optically-luminous FBOTs with luminous multi-wavelength counterparts., Comment: 21 pages, 6 figures, 6 tables. Submitted to ApJ

Tidal disruption events (TDEs) provide a unique opportunity to probe the stellar populations around supermassive black holes (SMBHs). By combining light curve modeling with spectral line information and knowledge about the stellar populations in the host galaxies, we are able to constrain the properties of the disrupted star for three TDEs. The TDEs in our sample have UV spectra, and measurements of the UV N III to C III line ratios enabled estimates of the nitrogen-to-carbon abundance ratios for these events. We show that the measured nitrogen line widths are consistent with originating from the disrupted stellar material dispersed by the central SMBH. We find that these nitrogen-to-carbon abundance ratios necessitate the disruption of moderately massive stars ($\gtrsim 1 - 2 M_\odot$). We determine that these moderately massive disruptions are over-represented by a factor of $\gtrsim 10^2$ when compared to the overall stellar population of the post-starburst galaxy hosts. This implies that SMBHs are preferentially disrupting higher mass stars, possibly due to ongoing top-heavy star formation in nuclear star clusters or to dynamical mechanisms that preferentially transport higher mass stars to their tidal radii., Comment: 16 pages, 6 figures, 1 table, accepted to ApJ

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

Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed sub-population of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ~10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ~85% and ~60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction $f_\mathrm{B}\sim 0.05$ our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift $z \in [0,2.5]$, a tentative comparison with the metallicity distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances., Comment: Accepted for publication in A&A Letters, 12 pages, 6 figures, 1 table

The detection of the unusually heavy binary neutron star merger GW190425 marked a stark contrast to the mass distribution from known Galactic pulsars in double neutron star binaries and gravitational-wave source GW170817. We suggest here a formation channel for heavy binary neutron stars and light black hole - neutron star binaries in which massive helium stars, which had their hydrogen envelope removed during a common envelope phase, remain compact and avoid mass transfer onto the neutron star companion, possibly avoiding pulsar recycling. We present three-dimensional simulations of the supernova explosion of the massive stripped helium star and follow the mass fallback evolution and the subsequent accretion onto the neutron star companion. We find that fallback leads to significant mass growth in the newly formed neutron star. This can explain the formation of heavy binary neutron star systems such as GW190425, as well as predict the assembly of light black hole - neutron star systems such as GW200115. This formation avenue is consistent with the observed mass-eccentricity correlation of binary neutron stars in the Milky Way. Finally, avoiding mass transfer suggests an unusually long spin-period population of pulsar binaries in our Galaxy., Comment: 7 main pages, 3 main figures, plus appendices. Accepted for publication in ApJ Letters

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources., Comment: 86 pages, 9 figures

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $\gamma$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind., Comment: 19 pages, 14 figures, published in ApJ

Our understanding of the formation and evolution of binary black holes (BBHs) is significantly impacted by the recent discoveries made by the LIGO/Virgo collaboration. Of utmost importance is the detection of the most massive BBH system, GW190521. Here we investigate what it takes for field massive stellar binaries to account for the formation of such massive BBHs. Whether the high mass end of the BH mass function is populated by remnants of massive stars that either formed at extremely low metallicities and avoid the pair-instability mass gap or increase their birth mass beyond the pair-instability mass gap through the accretion of gas from the surrounding medium. We show that assuming that massive stars at very low metallicities can form massive BHs by avoiding pair-instability supernova, coupled with a correspondingly high formation efficiency for BBHs, can explain the observed BH mass function. To this end, one requires a relation between the initial and final mass of the progenitor stars at low metallicities that is shallower than what is expected from wind mass loss alone. On the other hand, assuming pair-instability operates at all metallicities, one can account for the observed BH mass function if at least about 10% of the BHs born at very low metallicities double their mass before they merge because of accretion of ambient gas. Such BBHs will have to spend about a Gyr within a parsec length-scale of their parent atomic cooling halos or a shorter timescale if they reside in the inner sub-parsecs of their host dark matter halos. Future stellar evolution calculations of massive stars at very low metallicity and hydrodynamical simulations of gas accretion onto BBHs born in atomic cooling halos can shed light on this debate., Comment: 10 pages, Comments are welcome

Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transient types, such as active galactic nuclei (AGN) and supernovae (SNe), which can make TDE classification ambiguous. A further complication is that "TDE signatures" have not been uniformly observed to similar sensitivities or even targeted across all candidates. This chapter reviews those events that are unusual relative to other TDEs, including the possibility of TDEs in pre-existing AGN, and summarises those characteristics thought to best distinguish TDEs from continuously accreting AGN, strongly flaring AGN, SNe, and Gamma-Ray Bursts (GRBs), as well as other potential impostors like stellar collisions, "micro-TDEs," and circumbinary accretion flows. We conclude that multiple observables should be used to classify any one event as a TDE. We also consider the TDE candidate population as a whole, which, for certain host galaxy or SMBH characteristics, is distinguishable statistically from non-TDEs, suggesting that at least some TDE candidates do in fact arise from SMBH-disrupted stars., Comment: 57 pages, 17 figures. Accepted for publication in Springer Space Science Reviews. Chapter in ISSI review "The Tidal Disruption of Stars by Massive Black Holes" vol. 79. Corrected several typos from previous submission

Stellar-mass black holes can become embedded within the gaseous disks of active galactic nuclei (AGNs). Afterwards, their interactions are mediated by their gaseous surroundings. In this work, we study the evolution of stellar-mass binary black holes (BBHs) embedded within AGN disks using a combination of three-dimensional hydrodynamic simulations and analytic methods, focusing on environments in which the AGN disk scale height $H$ is $\gtrsim$ 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, which range from wind-dominated to quasi-spherical. We use our simulations to develop prescriptions for mass accretion and drag for embedded BBHs. We use these prescriptions, along with AGN disk models that can represent the Toomre-unstable outer regions of AGN disks, to study the long-term evolution of the BBHs as they migrate through the disk. We find that BBHs typically merge within $\lesssim 5-30\,{\rm Myr}$, increasing their mass significantly in the process, allowing BBHs to enter (or cross) the pair-instability supernova mass gap. The rate at which gas is supplied to these BBHs often exceeds the Eddington limit, sometimes by several orders of magnitude. We conclude that most embedded BBHs will merge before migrating significantly in the disk. Depending on the conditions of the ambient gas and the distance to the system, LISA can detect the transition between the gas-dominated and gravitational wave dominated regime for inspiraling BBHs that are formed sufficiently close to the AGN ($\lesssim$ 0.1 pc). We also discuss possible electromagnetic signatures during and following the inspiral, finding that it is generally unlikely but not inconceivable for the bolometric luminosity of the BBH to exceed that of the host AGN., Comment: 19 pages, 8 figures, accepted to ApJ

We present Hubble Space Telescope imaging of a pre-explosion counterpart to SN 2019yvr obtained 2.6 years before its explosion as a type Ib supernova (SN Ib). Aligning to a post-explosion Gemini-S/GSAOI image, we demonstrate that there is a single source consistent with being the SN 2019yvr progenitor system, the second SN Ib progenitor candidate after iPTF13bvn. We also analyzed pre-explosion Spitzer/IRAC imaging, but we do not detect any counterparts at the SN location. SN 2019yvr was highly reddened, and comparing its spectra and photometry to those of other, less extinguished SNe Ib we derive $E(B-V)=0.51\substack{+0.27\\-0.16}$ mag for SN 2019yvr. Correcting photometry of the pre-explosion source for dust reddening, we determine that this source is consistent with a $\log(L/L_{\odot}) = 5.3 \pm 0.2$ and $T_{\mathrm{eff}} = 6800\substack{+400\\-200}$ K star. This relatively cool photospheric temperature implies a radius of 320$\substack{+30\\-50} R_{\odot}$, much larger than expectations for SN Ib progenitor stars with trace amounts of hydrogen but in agreement with previously identified SN IIb progenitor systems. The photometry of the system is also consistent with binary star models that undergo common envelope evolution, leading to a primary star hydrogen envelope mass that is mostly depleted but seemingly in conflict with the SN Ib classification of SN 2019yvr. SN 2019yvr had signatures of strong circumstellar interaction in late-time ($>$150 day) spectra and imaging, and so we consider eruptive mass loss and common envelope evolution scenarios that explain the SN Ib spectroscopic class, pre-explosion counterpart, and dense circumstellar material. We also hypothesize that the apparent inflation could be caused by a quasi-photosphere formed in an extended, low-density envelope or circumstellar matter around the primary star., Comment: 22 pages, 9 figures, submitted to MNRAS