Search

We present a comprehensive, configurable open-source framework for estimating the rate of electromagnetic detection of kilonovae (KNe) associated with gravitational wave detections of binary neutron star (BNS) mergers. We simulate the current LIGO-Virgo-KAGRA (LVK) observing run (O4) using up-to-date sensitivity and up-time values as well as the next observing run (O5) using predicted sensitivities. We find the number of discoverable kilonovae during LVK O4 to be ${ 1}_{- 1}^{+ 4}$ or ${ 2 }_{- 2 }^{+ 3 }$, (at 90% confidence) depending on the distribution of NS masses in coalescing binaries, with the number increasing by an order of magnitude during O5 to ${ 19 }_{- 11 }^{+ 24 }$. Regardless of mass model, we predict at most five detectable KNe (at 95% confidence) in O4. We also produce optical and near-infrared light curves that correspond to the physical properties of each merging system. We have collated important information for allocating observing resources and directing search and follow-up observations including distributions of peak magnitudes in several broad bands and timescales for which specific facilities can detect each KN. The framework is easily adaptable, and new simulations can quickly be produced as input information such as merger rates and NS mass distributions are refined. Finally, we compare our suite of simulations to the thus-far completed portion of O4 (as of October 14, 2023), finding a median number of discoverable KNe of 0 and a 95-percentile upper limit of 2, consistent with no detection so far in O4., Comment: 16 pages, 13 figures, submitted to MNRAS

We motivate the cosmological science case of measuring Type Ia supernovae with the Nancy Grace Roman Space Telescope as part of the High Latitude Time Domain Survey. We discuss previously stated requirements for the science, and a baseline survey strategy. We discuss the various areas that must still be optimized and point to the other white papers that consider these topics in detail. Overall, the baseline case should enable an exquisite measurement of dark energy using SNe Ia from z=0.1 to z>2, and further optimization should only strengthen this once-in-a-generation experiment.

We present results from simulating slitless spectroscopic observations with the Nancy Grace Roman Space Telescope's (Roman) Wide-Field Instrument (WFI) P127 prism spanning 0.75 $\mu m$ to 1.8 $\mu m$. We quantify the efficiency of recovered Type Ia supernovae (SNe Ia) redshifts, as a function of P127 prism exposure time, to guide planning for future observing programs with the Roman prism. Generating the two-dimensional dispersed images and extracting one-dimensional spectra is done with the slitless spectroscopy package pyLINEAR along with custom-written software. From the analysis of 1698 simulated SN Ia P127 prism spectra, we show the efficiency of recovering SN redshifts to $z\lesssim3.0$, highlighting the exceptional sensitivity of the Roman P127 prism. Redshift recovery is assessed by setting a requirement of $\sigma_z = (\left|z - z_\mathrm{true} \right|)/(1+z) \leq 0.01$. We find that 3 hr exposures are sufficient for meeting this requirement, for $\gtrsim 50\%$ of the sample of mock SNe Ia at $z\approx2$ and within $\pm5$ days of rest-frame maximum light in the optical. We also show that a 1 hr integration of Roman can achieve the same precision in completeness to a depth of $24.4 \pm 0.06$ AB mag (or $z\lesssim 1$). Implications for cosmological studies with Roman P127 prism spectra of SNe Ia are also discussed., Comment: ApJ accepted. Updated to match accepted version

NASA will launch the Nancy Grace Roman Space Telescope (Roman) in the second half of this decade, which will allow for a generation-defining measurement of dark energy through multiple probes, including Type Ia supernovae (SNe Ia). To improve decisions on survey strategy, we have created the first simulations of realistic Roman images that include artificial SNe Ia injected as point sources in the images. Our analysis combines work done on Roman simulations for weak gravitational lensing studies as well as catalog-level simulations of SN samples. We have created a time series of images over two years containing $\sim$ 1,050 SNe Ia, covering a 1 square degree subarea of a planned 5 square degree deep survey. We have released these images publicly for community use along with input catalogs of all injected sources. We create secondary products from these images by generating coadded images and demonstrating recovery of transient sources using image subtraction. We perform first-use analyses on these images in order to measure galaxy-detection efficiency, point source-detection efficiency, and host-galaxy association biases. The simulated images can be found here: https://roman.ipac.caltech.edu/sims/SN_Survey_Image_sim.html., Comment: 12 pages, 12 figures. Submitted to MNRAS. For simulated images see https://roman.ipac.caltech.edu/sims/SN_Survey_Image_sim.html

Common envelope evolution, the key orbital tightening phase of the traditional formation channel for close binaries, is a multistage process that presents many challenges to the establishment of a fully descriptive, predictive theoretical framework. In an approach complementary to global 3D hydrodynamical modeling, we explore the range of applicability for a simplified drag formalism that incorporates the results of local hydrodynamic "wind tunnel" simulations into a semi-analytical framework in the treatment of the common envelope dynamical inspiral phase using a library of realistic giant branch stellar models across the low, intermediate, and high mass regimes. In terms of a small number of key dimensionless parameters, we characterize a wide range of common envelope events, revealing the broad range of applicability of the drag formalism as well its self-similar nature across mass regimes and ages. Limitations arising from global binary properties and local structural quantities are discussed together with the opportunity for a general prescriptive application for this formalism., Comment: 16 pages, 12 figures

Supernova (SN) feedback plays a vital role in the evolution of galaxies. While modern cosmological simulations capture the leading structures within galaxies, they struggle to provide sufficient resolution to study small-scale stellar feedback, such as the detailed evolution of SN remnants. It is thus common practice to assume subgrid models that are rarely extended to low metallicities, and which routinely use the standard solar abundance pattern. With the aid of 1-d hydrodynamical simulations, we extend these models to consider low metallicities and non-solar abundance patterns as derived from spectra of Milky Way stars. For that purpose, a simple, yet effective framework has been developed to generate non-solar abundance pattern cooling functions. We find that previous treatments markedly over-predict SN feedback at low metallicities and show that non-negligible changes in the evolution of SN remnants of up to $\approx 50\%$ in $cooling\; mass$ and $\approx 27\%$ in $momentum\; injection\; from\; SN\; remnants$ arise from non-solar abundance patterns. We use our simulations to quantify these results as a function of metallicity and abundance pattern variations and present analytic formulae to accurately describe the trends. These formulae have been designed to serve as subgrid models for SN feedback in cosmological hydrodynamical simulations., Comment: You can find the code to generate the cooling functions for an arbitrary composition at https://github.com/pikarpov/CoolingCurve

High-resolution spectroscopy of stars on the red giant branch (RGB) of the globular cluster M15 has revealed a large (~1 dex) dispersion in the abundances of r-process elements, like Ba and Eu. Neutron star mergers (NSMs) have been proposed as a major source of the r-process. However, most NSM models predict a delay time longer than the timescale for cluster formation. One possibility is that a NSM polluted the surfaces of stars in M15 long after the cluster finished forming. In this case, the abundances of the polluting elements would decrease in the first dredge-up as stars turn on to the RGB. We present Keck/DEIMOS abundances of Ba in 66 stars along the entire RGB and the top of the main sequence. The Ba abundances have no trend with stellar luminosity (evolutionary phase). Therefore, the stars were born with the Ba they have today, and Ba did not originate in a source with a delay time longer than the timescale for cluster formation. In particular, if the source of Ba was a neutron star merger, it would have had a very short delay time. Alternatively, if Ba enrichment took place before the formation of the cluster, an inhomogeneity of a factor of 30 in Ba abundance needs to be able to persist over the length scale of the gas cloud that formed M15, which is unlikely., Comment: 10 pages, 2 figures, accepted to ApJL

We use observations of heavy elements in very metal-poor stars ([Fe/H] < -2.5) in order to place constraints on the viability of collapsar models as a significant source of the r-process. We combine bipolar explosion nucleosynthesis calculations with recent disk calculations to make predictions of the observational imprints these explosions would leave on very metal-poor stars. We find that a source of low (~ 0.1-0.5 $M_\odot$) Fe mass which also yields a relatively high (> 0.08 $M_\odot$) r-process mass would, after subsequently mixing and forming new stars, result in [r/Fe] abundances up to three orders of magnitude higher than those seen in stars. In order to match inferred abundances, 10-10$^3 M_\odot$ of Fe would need be efficiently incorporated into the r-process ejecta. We show that Fe enhancement and hence [r/Fe] dilution from other nearby supernovae is not able to explain the observations unless significant inflow of pristine gas occurs before the ejecta are able to form new stars. Finally, we show that the inferred [Eu/Fe] abundances require levels of gas mixing which are in conflict with other properties of r-process enhanced metal-poor stars. Our results suggest that early r-process production is likely to be spatially uncorrelated with Fe production, a condition which can be satisfied by neutron star mergers due to their large kick velocities and purely r-process yields., Comment: 6 pages, 2 figures, accepted for publication in ApJL

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network., Comment: White paper submitted to the Astro2020 Decadal Survey

The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow up kilonova observations have shown that NSMs produce copious amount of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of Double Neutron Stars (DNSs), with the levels of r-process enrichment seen in ultra faint dwarf (UFD) galaxies such as Reticulum II and Tucana III. Assuming that such dwarf systems have lost a significant fraction of their stellar mass through tidal stripping, we conclude that contrary to most current models, it is the DNSs with rather large natal kicks but very short merging timescales that can enrich UFD-type galaxies. These binaries are either on highly eccentric orbits, or form with very short separations due to an additional mass-transfer between the first-born neutron star and a naked helium star, progenitor of the second-born neutron star. These DNSs are born with a frequency that agrees with the statistics of the r-process UFDs, and merge well within the virial radius of their host halos, therefore contributing significantly to their r-process enrichment., Comment: Accepted for publication in ApJ

During a common envelope episode in a binary system, the engulfed companion spirals to tighter orbital separations under the influence of drag from the surrounding envelope material. As this object sweeps through material with a steep radial gradient of density, net angular momentum is introduced into the flow, potentially leading to the formation of an accretion disk. The presence of a disk would have dramatic consequences for the outcome of the interaction because accretion might be accompanied by strong, polar outflows with enough energy to unbind the entire envelope. Without a detailed understanding of the necessary conditions for disk formation during common envelope, therefore, it is difficult to accurately predict the population of merging compact binaries. This paper examines the conditions for disk formation around objects embedded within common envelopes using the `wind tunnel' formalism developed by MacLeod et al. (2017). We find that the formation of disks is highly dependent on the compressibility of the envelope material. Disks form only in the most compressible of stellar envelope gas, found in envelopes' outer layers in zones of partial ionization. These zones are largest in low-mass stellar envelopes, but comprise small portions of the envelope mass and radius in all cases. We conclude that disk formation and associated accretion feedback in common envelope is rare, and if it occurs, transitory. The implication for LIGO black hole binary assembly is that by avoiding strong accretion feedback, common envelope interactions should still result in the substantial orbital tightening needed to produce merging binaries., Comment: 12 pages, 10 figures, submitted to ApJ

This paper examines the properties of flows around objects embedded within common envelopes in the simplified context of a "wind tunnel." We establish characteristic relationships between key common envelope flow parameters like the Mach number and density scale height. Our wind tunnel is a three-dimensional, cartesian geometry hydrodynamic simulation setup that includes the gravity of the primary and secondary stars and allows us to study the coefficients of drag and accretion experienced by the embedded object. Accretion and drag lead to a transformation of an embedded object and its orbit during a common envelope phase. We present two suites of simulations spanning a range of density gradients and Mach numbers -- relevant for flow near the limb of a stellar envelope to the deep interior. In one suite, we adopt an ideal gas adiabatic exponent of $\gamma=5/3$, in the other, $\gamma=4/3$. We find that coefficients of drag rise in flows with steeper density gradients and that coefficients of drag and accretion are consistently higher in the more compressible, $\gamma=4/3$ flow. We illustrate the impact of these newly derived coefficients by integrating the inspiral of a secondary object through the envelopes of $3M_\odot$ ($\gamma\approx5/3$) and $80M_\odot$ ($\gamma\approx4/3$) giants. In these examples, we find a relatively rapid initial inspiral because, near the stellar limb, dynamical friction drag is generated mainly from dense gas focussed from deeper within the primary-star's envelope. This rapid initial inspiral timescale carries potential implications for the timescale of transients from early common envelope interaction.

We construct a menu of objects that can give rise to bright flares when disrupted by massive black holes (BHs), ranging from planets to evolved stars. Through their tidal disruption, main sequence and evolved stars can effectively probe the existence of otherwise quiescent supermassive BHs and white dwarfs can probe intermediate mass BHs. Many low-mass white dwarfs possess extended hydrogen envelopes, which allow for the production of prompt flares in disruptive encounters with moderately massive BHs of $10^5$ to $10^7~M_\odot$--masses that may constitute the majority of massive BHs by number. These objects are a missing link in two ways: (1) for probing moderately massive BHs and (2) for understanding the hydrodynamics of the disruption of objects with tenuous envelopes. A flare arising from the tidal disruption of a $0.17~M_\odot$ white dwarf by a $10^5~M_\odot$ BH reaches a maximum between 0.6 and 11 days, with a peak fallback rate that is usually super-Eddington and results in a flare that is likely brighter than a typical tidal disruption event. Encounters stripping only the envelope can provide hydrogen-only fallback, while encounters disrupting the core evolve from H- to He-rich fallback. While most tidal disruption candidates observed thus far are consistent with the disruptions of main sequence stars, the rapid timescales of nuclear transients such as Dougie and PTF10iya are naturally explained by the disruption of low-mass white dwarfs. As the number of observed flares continues to increase, the menu presented here will be essential for characterizing nuclear BHs and their environments through tidal disruptions., Comment: 19 pages, 13 figures. Published in ApJ

We analyze data from several studies of metal-poor stars in the Milky Way, focusing on both strong (Eu) and weak (Sr) $r$-process elements. Because these elements were injected in an explosion, we calculate the mass swept up when the blast wave first becomes radiative, yielding a lower limit for the dilution of such elements and hence a lower limit on the ejecta mass which is incorporated into the next generation of stars. Our study demonstrates that in order to explain the largest enhancements in [Eu/Fe] observed in stars at low [Fe/H] metallicities, individual $r$-process production events must synthesize a minimum of $10^{-3.5} M_{\odot}$ of $r$-process material. We also show that if the site of Mg production is the same as that of Eu, individual injection events must synthesize up to $ \sim 10^{-3} M_{\odot}$ of $r$-process material. On the other hand, demanding that Sr traces Mg production results in $r$-process masses per event of $\sim 10^{-5} M_{\odot}$. This suggests that the astrophysical sites responsible for the genesis of the strong $r$-process elements need to operate at a drastically reduced rate when compared to core-collapse supernovae, while the synthesis of the weak $r$-process material is consistent with a supernova production site., Comment: Submitted to ApJL. 6 pages, 5 figures. Comments welcome!

This paper investigates the recent stellar merger transient M31LRN 2015 in the Andromeda galaxy. We analyze published optical photometry and spectroscopy along with a Hubble Space Telescope detection of the color and magnitude of the pre-outburst source. The transient outburst is consistent with dynamically driven ejecta at the onset of a common envelope episode, which eventually leads to the complete merger of a binary system. The light curve appears to contain two components: first $\sim10^{-2} M_\odot$ of fast ejecta driven by shocks at the onset of common envelope, and later, $\sim0.3 M_\odot$ of further ejecta as the secondary becomes more deeply engulfed within the primary. Just prior to merger, we find that the primary star is a $3-5.5 M_\odot$ sub-giant branch primary star with radius of $30-40R_\odot$. Its position in the color-magnitude diagram shows that it is growing in radius, consistent with a picture where it engulfs its companion. By matching the properties of the primary star to the transient outburst, we show that the optical transient lasts less than ten orbits of the original binary, which had a pre-merger period of $\sim 10$ days. We consider the possible orbital dynamics leading up to the merger, and argue that if the system merged due to the Darwin tidal instability it implies a lower mass main sequence companion of $0.1-0.6M_\odot$. This analysis represents a promising step toward a more detailed understanding of flows in common envelope episodes through direct observational constraints., Comment: Revised following report from referee. 21 pages, 11 figures

We present a method for modeling the evolution of detached double white dwarf (DWD) binaries hosting helium donors from the end of the common envelope (CE) phase to the onset of Roche Lobe overflow (RLOF). This is achieved by combining detailed stellar evolution calculations of extremely low mass (ELM) helium WDs possessing hydrogen envelopes with the the orbital shrinking of the binary driven by gravitational radiation. We show that the consideration of hydrogen fusion in these systems is crucial, as a significant fraction ($\approx$50%) of future donors are expected to still be burning when mass transfer commences. We apply our method to two detached eclipsing DWD systems, SDSS J0651+2844 and NLTT-11748, in order to demonstrate the effect that carbon-nitrogen-oxygen (CNO) flashes have on constraining the evolutionary history of such systems. We find that when CNO flashes are absent on the low mass WD ($M_{2}$ < $0.18 M_{\odot}$), such as in NLTT-11748, we are able to self consistently solve for the donor conditions at CE detachment given a reliable cooling age from the massive WD companion. When CNO flashes occur (0.18 $M_{\odot}$ < $M_{2}$ < 0.36 $M_{\odot}$), such as in SDSS J0651+2844, the evolutionary history is eradicated and we are unable to comment on the detachment conditions. We find that for any donor mass our models are able to predict the conditions at reattachment and comment on the stabilizing effects of hydrogen envelopes. This method can be applied to a population of detached DWDs with measured donor radii and masses., Comment: 5 pages, 5 figures, submitted to ApJL. Comments welcome

All evolved stars with masses M <2 solar masses undergo an initiating off-center helium core flash in their 0.48 solar mass He core as they ascend the red giant branch (RGB). This off-center flash is the first of a few successive helium shell subflashes that remove the core electron degeneracy over 2 Myrs, converting the object into a He burning star. Though characterized by Thomas over 40 years ago, this core flash phase has yet to be observationally probed. Using the Modules for Experiments in Stellar Astrophysics (MESA) code, we show that red giant asteroseismology enabled by space-based photometry (i.e. Kepler and CoRoT) can probe these stars during the flash. The rapid (< 100,000 years) contraction of the red giant envelope after the initiating flash dramatically improves the coupling of the p-modes to the core g-modes, making the detection of l=1 mixed modes possible for these 2 Myrs. This duration implies that 1 in 35 stars near the red clump in the HR diagram will be in their core flash phase. During this time, the star has a g-mode period spacing of Delta P_g = 70-100 seconds, lower than the Delta P_g=250 seconds of He burning stars in the red clump, but higher than the RGB stars at the same luminosity. This places them in an underpopulated part of the large frequency spacing (Delta nu) vs. Delta P_g diagram that should ease their identification amongst the thousands of observed red giants., Comment: 15 pages, 5 figures, To Appear in Astrophysical Journal Letters