Search

The $\gamma$-process in core-collapse supernovae (CCSNe) can produce a number of neutron-deficient stable isotopes heavier than iron (p-nuclei). However, current model predictions do to not fully reproduce the solar abundances. We investigate the impact of different explosion energies and parameters on the nucleosynthesis of p-nuclei, by studying stellar models with different initial masses and CCSN explosions. We find that the total p-nuclei yields are only marginally affected by the CCSN explosion prescriptions if the $\gamma$-process production is already efficient in the stellar progenitors due to a C-O shell merger. In most of CCSN explosions from progenitors without C-O shell merger, the $\gamma$-process yields increase with the explosion energy up to an order of magnitude, depending on the progenitor structure and the CCSN prescriptions. The trend of the p-nuclei production with the explosion energy is more complicated if we look at the production of single p-nuclei. The light p-nuclei tend to be the most enhanced with increasing the explosion energy. In particular, for the CCSN models where the $\alpha$-rich freeze-out component is ejected, the yields of the lightest p-nuclei increase by up to three orders of magnitude. We provide the first extensive study using different sets of massive stars of the impact of varying CCSN explosion prescriptions on the production of the p-nuclei. Unlike previous expectations and recent results in the literature, we find that the average production of p-nuclei tends to increase with the explosion energy. We also confirm that the pre-explosive production of p-nuclei in C-O shell mergers is a robust result, independently from the subsequent explosive nucleosynthesis. A realistic range of variations in the evolution of stellar progenitors and in the CCSN explosions might boost the CCSN contribution to the galactic chemical evolution of p-nuclei., Comment: 11 pages, 3 figures, 4 tables. Accepted for publication in A&A

Kilonovae are the electromagnetic transients created by the radioactive decay of freshly synthesized elements in the environment surrounding a neutron star merger. To study the fundamental physics in these complex environments, kilonova modeling requires, in part, the use of radiative transfer simulations. The microphysics involved in these simulations results in high computational cost, prompting the use of emulators for parameter inference applications. Utilizing a training set of 22248 high-fidelity simulations (composed of 412 unique ejecta parameter combinations evaluated at 54 viewing angles), we use a neural network to efficiently train on existing radiative transfer simulations and predict light curves for new parameters in a fast and computationally efficient manner. Our neural network can generate millions of new light curves in under a minute. We discuss our emulator's degree of off-sample reliability and parameter inference of the AT2017gfo observational data. Finally, we discuss tension introduced by multi-band inference in the parameter inference results, particularly with regard to the neural network's recovery of viewing angle., Comment: 16 pages, 11 figures

The interaction of $\beta$-particles with the weakly ionized plasma background is an important mechanism for powering the kilonova transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti (1971) following the seminal formulation of Bethe (1930), in a spectral solver of the Vlasov-Maxwell-Boltzmann equations. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon-matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric kilonova ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for $\beta$-particles in the kilonova are a bulk velocity parameter $\vec{u}=0$, a thermal velocity parameter $\vec{\alpha}=0.5c$, and a 9x9x9 mode velocity basis set (Hermite orders 0 to 8 in each dimension). For ejecta-interior sample zones, we estimate the ratio of thermalization from large-angle ($\gtrsim2.5^{\circ}$) bound excitation scattering to total thermalization is $\sim$0.002-0.003., Comment: Accepted to ApJ, 26 pages, 8 figures

The electromagnetic emission from the non-relativistic ejecta launched in neutron star mergers (either dynamically or through a disk wind) has the potential to probe both the total mass and composition of this ejecta. These observations are crucial in understanding the role of these mergers in the production of r-process elements in the universe. However, many properties of the ejecta can alter the light-curves and we must both identify which properties play a role in shaping this emission and understand the effects these properties have on the emission before we can use observations to place strong constraints on the amount of r-process elements produced in the merger. This paper focuses on understanding the effect of the velocity distribution (amount of mass moving at different velocities) for lanthanide-rich ejecta on the light-curves and spectra. The simulations use distributions guided by recent calculations of disk outflows and compare the velocity-distribution effects to those of ejecta mass, velocity and composition. Our comparisons show that uncertainties in the velocity distribution can lead to factor of 2-4 uncertainties in the inferred ejecta mass based on peak infra-red luminosities. We also show that early-time UV or optical observations may be able to constrain the velocity distribution, reducing the uncertainty in the ejecta mass., Comment: 21 pages, 16 figures, submitted to ApJ

HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<10"$ full width at half maximum) and broad spectral coverage (0.2--80 keV) with an effective area far superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. HEX-P is ideally suited to address important problems in the physics and astrophysics of supernova remnants (SNRs) and pulsar-wind nebulae (PWNe). For shell SNRs, HEX-P can greatly improve our understanding via more accurate spectral characterization and localization of non-thermal X-ray emission from both non-thermal-dominated SNRs and those containing both thermal and non-thermal components, and can discover previously unknown non-thermal components in SNRs. Multi-epoch HEX-P observations of several young SNRs (e.g., Cas A and Tycho) are expected to detect year-scale variabilities of X-ray filaments and knots, thus enabling us to determine fundamental parameters related to diffusive shock acceleration, such as local magnetic field strengths and maximum electron energies. For PWNe, HEX-P will provide spatially-resolved, broadband X-ray spectral data separately from their pulsar emission, allowing us to study how particle acceleration, cooling, and propagation operate in different evolution stages of PWNe. HEX-P is also poised to make unique and significant contributions to nuclear astrophysics of Galactic radioactive sources by improving detections of, or limits on, $^{44}$Ti in the youngest SNRs and by potentially discovering rare nuclear lines as evidence of double neutron star mergers. Throughout the paper, we present simulations of each class of objects, demonstrating the power of both the imaging and spectral capabilities of HEX-P to advance our knowledge of SNRs, PWNe, and nuclear astrophysics., Comment: 34 pages, 15 figures; part of a suite of papers describing the HEX-P hard X-ray mission concept

HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging (<10" FWHM) and broad spectral coverage (0.2-80 keV) with an effective area far superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. With the recent discoveries of over 40 ultra-high-energy gamma-ray sources (detected above 100 TeV) and neutrino emission in the Galactic Plane, we have entered a new era of multi-messenger astrophysics facing the exciting reality of Galactic PeVatrons. In the next decade, as more Galactic PeVatrons and TeV gamma-ray sources are expected to be discovered, the identification of their acceleration and emission mechanisms will be the most pressing issue in both particle and high-energy astrophysics. In this paper, along with its companion papers (Reynolds et al. 2023, Mori et al. 2023), we will present that HEX-P is uniquely suited to address important problems in various cosmic-ray accelerators, including Galactic PeVatrons, through investigating synchrotron X-ray emission of TeV-PeV electrons produced by both leptonic and hadronic processes., Comment: 31 pages, 12 figures, submitted to FrASS

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

Core-collapse explosions of massive stars leave behind neutron stars, with a known diversity that includes the "Central Compact Objects" (CCOs). Typified by the neutron star discovered near the centre of the Cas A supernova remnant (SNR), CCOs have been observed to shine only in X-rays. To address their supernova progenitors, we perform a systematic study of SNRs that contain a CCO and display X-ray emission from their shock-heated ejecta. We make use of X-ray data primarily using the Chandra X-ray observatory, complemented with XMM-Newton. This study uses a systematic approach to the analysis of each SNR aimed at addressing the supernova progenitor as well as the explosion properties (energy and ambient density). After fitting for the ejecta abundances estimated from a spatially resolved spectroscopic study, we compare the data to six nucleosynthesis models making predictions on supernova ejecta yields in core-collapse explosions. We find that the explosion models commonly used by the astrophysics community do not match the ejecta yields for any of the SNRs, suggesting additional physics, e.g. multi-dimensional explosion models or updated progenitor structures, are required. Overall we find low-mass ($\leq$25 solar masses) progenitors among the massive stars population and low-energy explosions ($<$10$^{51}$ ergs). We discuss degeneracies in our model fitting, particularly how altering the explosion energy affects the estimate of the progenitor mass. Our systematic study highlights the need for improving on the theoretical models for nucleosynthesis predictions as well as for sensitive, high-resolution spectroscopy observations to be acquired with next-generation X-ray missions., Comment: 29 pages, 64 figures, 17 Tables; to appear in MNRAS

The Compton Spectrometer and Imager (COSI) is a NASA Small Explorer (SMEX) satellite mission in development with a planned launch in 2027. COSI is a wide-field gamma-ray telescope designed to survey the entire sky at 0.2-5 MeV. It provides imaging, spectroscopy, and polarimetry of astrophysical sources, and its germanium detectors provide excellent energy resolution for emission line measurements. Science goals for COSI include studies of 0.511 MeV emission from antimatter annihilation in the Galaxy, mapping radioactive elements from nucleosynthesis, determining emission mechanisms and source geometries with polarization measurements, and detecting and localizing multimessenger sources. The instantaneous field of view for the germanium detectors is >25% of the sky, and they are surrounded on the sides and bottom by active shields, providing background rejection as well as allowing for detection of gamma-ray bursts and other gamma-ray flares over most of the sky. In the following, we provide an overview of the COSI mission, including the science, the technical design, and the project status., Comment: 8 pages

The Compton Spectrometer and Imager (COSI) is a selected Small Explorer (SMEX) mission launching in 2027. It consists of a large field-of-view Compton telescope that will probe with increased sensitivity the under-explored MeV gamma-ray sky (0.2-5 MeV). We will present the current status of cosipy, a Python library that will perform spectral and polarization fits, image deconvolution, and all high-level analysis tasks required by COSI's broad science goals: uncovering the origin of the Galactic positrons, mapping the sites of Galactic nucleosynthesis, improving our models of the jet and emission mechanism of gamma-ray bursts (GRBs) and active galactic nuclei (AGNs), and detecting and localizing gravitational wave and neutrino sources. The cosipy library builds on the experience gained during the COSI balloon campaigns and will bring the analysis of data in the Compton regime to a modern open-source likelihood-based code, capable of performing coherent joint fits with other instruments using the Multi-Mission Maximum Likelihood framework (3ML). In this contribution, we will also discuss our plans to receive feedback from the community by having yearly software releases accompanied by publicly-available data challenges.

The Interplanetary Network (IPN) is a detection, localization and alert system that utilizes the arrival time of transient signals in gamma-ray detectors on spacecraft separated by planetary baselines to geometrically locate the origin of these transients. Due to the changing astrophysical landscape and the new emphasis on time domain and multi-messenger astrophysics (TDAMM) from the Pathways to Discovery in Astronomy and Astrophysics for the 2020s, this Gamma-ray Transient Network Science Analysis Group was tasked to understand the role of the IPN and high-energy monitors in this new era. The charge includes describing the science made possible with these facilities, tracing the corresponding requirements and capabilities, and highlighting where improved operations of existing instruments and the IPN would enhance TDAMM science. While this study considers the full multiwavelength and multimessenger context, the findings are specific to space-based high-energy monitors. These facilities are important both for full characterization of these transients as well as facilitating follow-up observations through discovery and localization. The full document reports a brief history of this field, followed by our detailed analyses and findings in some 68 pages, providing a holistic overview of the role of the IPN and high-energy monitors in the coming decades., Comment: Terms of Reference and additional information on the Science Analysis Group are available at https://pcos.gsfc.nasa.gov/sags/gtn-sag.php

Kilonovae are a rare class of astrophysical transients powered by the radioactive decay of nuclei heavier than iron, synthesized in the merger of two compact objects. Over the first few days, the kilonova evolution is dominated by a large number of radioactive isotopes contributing to the heating rate. On timescales of weeks to months, its behavior is predicted to differ depending on the ejecta composition and merger remnant. However, late-time observations of known kilonovae are either missing or limited. Here we report observations of a luminous red transient with a quasi-thermal spectrum, following an unusual gamma-ray burst of long duration. We classify this thermal emission as a kilonova and track its evolution up to two months after the burst. At these late times, the recession of the photospheric radius and the rapidly-decaying bolometric luminosity ($L_{\rm bol}\propto t^{-2.7\pm 0.4}$) support the recombination of lanthanide-rich ejecta as they cool., Comment: 47 pages, 14 figures, 9 tables; submitted; a minor typo fixed

Kilonovae are a novel class of astrophysical transients, and the only observationally-confirmed site of rapid neutron capture nucleosynthesis (the r-process) in the Universe. To date, only a handful of kilonovae have been detected, with just a single spectroscopically-observed event (AT 2017gfo). Spectra of AT 2017gfo provided evidence for the formation of elements heavier than iron; however, these spectra were collected during the first ~ 10 days, when emission from light r-process elements dominates the observations. Heavier elements, if synthesised, are expected to shape the late-time evolution of the kilonova, beyond the phases for which we have spectral observations. Here we present spectroscopic observations of a rapidly-reddening thermal transient, following the gamma-ray burst, GRB 230307A. Early (2.4 day) optical spectroscopy identifies the presence of a hot (T ~ 6700 K) thermal continuum. By 29 days, this component has expanded and cooled significantly (T ~ 640 K), yet it remains optically thick, indicating the presence of high-opacity ejecta. We show that these properties can only be explained by the merger of compact objects, and further, leads us to infer the production of the heavy lanthanide elements. We identify several spectral features (in both absorption and emission), whose cause can be explained by newly-synthesised heavy elements. This event marks only the second recorded spectroscopic evidence for the synthesis of r-process elements, and the first to be observed at such late times., Comment: 41 pages, 11 figures, 4 tables. Submitted - comments welcome!

The LIGO HET Response (LIGHETR) project is an enterprise to follow up optical transients (OT) discovered as gravitational wave merger sources by the LIGO/Virgo collaboration (LVC). Early spectroscopy has the potential to constrain crucial parameters such as the aspect angle. The LIGHETR collaboration also includes the capacity to model the spectroscopic evolution of mergers to facilitate a real-time direct comparison of models with our data. The principal facility is the Hobby-Eberly Telescope. LIGHETR uses the massively-replicated VIRUS array of spectrographs to search for associated OTs and obtain early blue spectra and in a complementary role, the low-resolution LRS-2 spectrograph is used to obtain spectra of viable candidates as well as a densely-sampled series of spectra of true counterparts. Once an OT is identified, the anticipated cadence of spectra would match or considerably exceed anything achieved for GW170817 = AT2017gfo for which there were no spectra in the first 12 hours and thereafter only roughly once daily. We describe special HET-specific software written to facilitate the program and attempts to determine the flux limits to undetected sources. We also describe our campaign to follow up OT candidates during the third observational campaign of the LIGO and Virgo Scientific Collaborations. We obtained VIRUS spectroscopy of candidate galaxy hosts for 5 LVC gravitational wave events and LRS-2 spectra of one candidate for the OT associated with S190901ap. We identified that candidate, ZTF19abvionh = AT2019pip, as a possible Wolf-Rayet star in an otherwise unrecognized nearby dwarf galaxy., Comment: 26 pages, 15 figures

The electromagnetic emission from neutron star mergers is comprised of multiple components. Synchrotron emission from the disk-powered jet as well as thermal emission from the merger ejecta (powered by a variety of sources) are among the most studied sources. The low masses and high velocities of the merger ejecta quickly develop conditions where emission from collisionless shocks becomes critical and synchrotron emission from the merger ejecta constitutes a third component to the observed signal. The aim of this project is to examine shock development, magnetic field generation and particle acceleration in the case of mildly relativistic shocks, which are expected when the tidal ejecta of neutron star mergers drive a shock into the external medium. Using LANL's VPIC (vector particle-in-cell) code, we have run simulations of such mildly-relativistic, collisionless, weakly-magnetized plasmas and compute the resultant magnetic fields and particle energy spectra. We show the effects of varying plasma conditions, as well as explore the validity of using different proton to electron mass ratios in VPIC. Our results have implications for observing late-time electromagnetic counterparts to gravitational wave detections of neutron star mergers., Comment: Submitted To ApJ

In this work, we present a simple interpolation methodology for spectroscopic time series, based on conventional interpolation techniques (random forests) implemented in widely-available libraries. We demonstrate that our existing library of simulations is sufficient for training, producing interpolated spectra that respond sensitively to varied ejecta parameter, post-merger time, and viewing angle inputs. We compare our interpolated spectra to the AT2017gfo spectral data, and find parameters similar to our previous inferences using broadband light curves. However, the spectral observations have significant systematic short-wavelength residuals relative to our models, which we cannot explain within our existing framework. Similar to previous studies, we argue that an additional blue component is required. We consider a radioactive heating source as a third component characterized by light, slow-moving, lanthanide-free ejecta with $M_{\rm th} = 0.003~M_\odot$, $v_{\rm th} = 0.05$c, and $\kappa_{\rm th} = 1$ cm$^2$/g. When included as part of our radiative transfer simulations, our choice of third component reprocesses blue photons into lower energies, having the opposite effect and further accentuating the blue-underluminosity disparity in our simulations. As such, we are unable to overcome short-wavelength deficits at later times using an additional radioactive heating component, indicating the need for a more sophisticated modeling treatment., Comment: 11 pages, 7 figures, submitted to Phys. Rev. Research, presented at April APS session F13.00006

We systematically investigate the mergers of neutron star-white dwarf binaries from beginning to end, with focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a large accretion disk around the neutron star. The material evolves according to an advection-dominated accretion disk model with nuclear burning, neutrino-emissions, and disk-surface wind ejection. The extreme dynamics of the entire process has proven difficult for analytic analysis, and thus currently the properties are poorly understood. The outflows from the mergers are iron- and nickel-rich, giving rise to optical and infrared emissions powered from the decay of the radioactive iron-type isotopes, calculated via the SuperNu light-curve code. We find these systems capable of powering bright, yet short-lived, optical transients with the potential to power gamma-ray bursts., Comment: LA-UR-22-29879

The All-sky Medium Energy Gamma-ray Observatory eXplorer (AMEGO-X) is designed to identify and characterize gamma rays from extreme explosions and accelerators. The main science themes include: supermassive black holes and their connections to neutrinos and cosmic rays; binary neutron star mergers and the relativistic jets they produce; cosmic ray particle acceleration sources including Galactic supernovae; and continuous monitoring of other astrophysical events and sources over the full sky in this important energy range. AMEGO-X will probe the medium energy gamma-ray band using a single instrument with sensitivity up to an order of magnitude greater than previous telescopes in the energy range 100 keV to 1 GeV that can be only realized in space. During its three-year baseline mission, AMEGO-X will observe nearly the entire sky every two orbits, building up a sensitive all-sky map of gamma-ray sources and emission. AMEGO-X was submitted in the recent 2021 NASA MIDEX Announcement of Opportunity., Comment: 23 pages, 16 figures, Published Journal of Astronomical Telescopes, Instruments, and Systems

High Reynolds Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

We present a comprehensive optical and near-infrared census of the fields of 90 short gamma-ray bursts (GRBs) discovered in 2005-2021, constituting all short GRBs for which host galaxy associations are feasible ($\approx$ 60% of the total Swift short GRB population). We contribute 245 new multi-band imaging observations across 49 distinct GRBs and 25 spectra of their host galaxies. Supplemented by literature and archival survey data, the catalog contains 335 photometric and 40 spectroscopic data sets. The photometric catalog reaches $3\sigma$ depths of $\gtrsim 24-27$ mag and $\gtrsim 23-26$ mag for the optical and near-infrared bands, respectively. We identify host galaxies for 84 bursts, in which the most robust associations make up 54% (49/90) of events, while only a small fraction, 6.7%, have inconclusive host associations. Based on new spectroscopy, we determine 17 host spectroscopic redshifts with a range of $z\approx 0.15-1.6$ and find that $\approx$ 25-44% of Swift short GRBs originate from $z>1$. We also present the galactocentric offset catalog for 83 short GRBs. Taking into account the large range of individual measurement uncertainties, we find a median of projected offset of $\approx 7.9$ kpc, for which the bursts with the most robust associations have a smaller median of $\approx 4.9$ kpc. Our catalog captures more high-redshift and low-luminosity hosts, and more highly-offset bursts than previously found, thereby diversifying the population of known short GRB hosts and properties. In terms of locations and host luminosities, the populations of short GRBs with and without detectable extended emission are statistically indistinguishable. This suggests that they arise from the same progenitors, or from multiple progenitors which form and evolve in similar environments. All of the data products are available on the BRIGHT website., Comment: 53 pages, 9 figures, 6 tables, submitted

Turbulence plays an important role in astrophysical phenomena, including core-collapse supernovae (CCSN), but current simulations must rely on subgrid models since direct numerical simulation (DNS) is too expensive. Unfortunately, existing subgrid models are not sufficiently accurate. Recently, Machine Learning (ML) has shown an impressive predictive capability for calculating turbulence closure. We have developed a physics-informed convolutional neural network (CNN) to preserve the realizability condition of Reynolds stress that is necessary for accurate turbulent pressure prediction. The applicability of the ML subgrid model is tested here for magnetohydrodynamic (MHD) turbulence in both the stationary and dynamic regimes. Our future goal is to utilize this ML methodology (available on GitHub) in the CCSN framework to investigate the effects of accurately-modeled turbulence on the explosion of these stars., Comment: For our ML algorithm on GitHub, see https://github.com/pikarpov-LANL/Sapsan/wiki/Estimators\#physics-informed-cnn-for-turbulence-modeling

Our understanding of the convective-engine paradigm driving core-collapse supernovae has been used for 2 decades to predict the remnant mass distribution from stellar collapse. These predictions improve as our understanding of this engine increases. In this paper, we review our current understanding of convection (in particular, the growth rate of convection) in stellar collapse and study its effect on the remnant mass distribution. We show how the depth of the mass gap between neutron stars and black holes can help probe this convective growth. We include a study of the effects of stochasticity in both the stellar structure and the convective seeds caused by stellar burning. We study the role of rotation and its effect on the pair-instability mass gap. Under the paradigm limiting stellar rotation to those stars in tight binaries, we determine the effect of rotation on the remnant mass distribution., Comment: 17 pages, accepted by ApJ

Understanding astrophysical phenomena involving compact objects requires an insight about the engine behind core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of detected population of BH-BH, BH-NS and NS-NS mergers. To gain this understanding, we must tie our current knowledge of pre-SN stars properties and their potential explosions to the final NS or BH mass distribution. The timescale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. In this study we adopt the new formulas for the relation between the pre-SN star properties and its remnant from Fryer et al. 2022 into StarTrack population synthesis code and check how they impact double compact object (DCO) mergers formed via isolated binary evolution. The new formulas give one ability to test a wide spectrum of assumptions on the convection growth time. In particular, different variants allow for a smooth transition between having a deep lower mass gap and a remnant mass distribution filled by massive NSs and low mass BHs. In this paper we present distribution of masses, mass ratios and the local merger rate densities of DCO mergers for different variants of new remnant mass formulas. We test them together with different approaches to other highly uncertain processes. We find that mass distribution of DCO mergers up to m_1+m_2 < 35 Msun is sensitive to adopted assumption on SN convection growth timescale. Between the two extreme tested variants the probability of compact object formation within the lower mass gap may differ up to 2 orders of magnitude. The mass ratio distribution of DCO mergers is significantly influenced by SN model only for our standard mass transfer stability criteria., Comment: 22 pages, version accepted for publication in MNRAS

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe.

Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 to 100 Myr can be used to probe the origin of the Solar System. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: $^{26}$Al, $^{36}$Cl, $^{41}$Ca, $^{53}$Mn, $^{60}$Fe, $^{92}$Nb, $^{97}$Tc, $^{98}$Tc, $^{107}$Pd, $^{126}$Sn, $^{129}$I, $^{135}$Cs, $^{146}$Sm, $^{182}$Hf, and $^{205}$Pb. We probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25M$_{\odot}$, explosion energies between 3.4$\times$10$^{50}$ and 1.8$\times$10$^{52}$ ergs and compact remnant masses between 1.5M$_{\odot}$and 4.89M$_{\odot}$. We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as $^{92}$Nb and $^{97}$Tc, are the most affected by the remnant mass, $^{92}$Nb varying by five orders of magnitude. Isotopes synthesised primarily in explosive C-burning and explosive He-burning, such as $^{60}$Fe, are most affected by explosion energies. $^{60}$Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15M$_{\odot}$model. The final yield of each examined SLR is used to compare to literature models.

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

We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration (t_E~270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge. HST imaging, conducted at eight epochs over an interval of six years, reveals a clear relativistic astrometric deflection of the background star's apparent position. Ground-based photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of the Earth's motion on the microlensing light curve. Combining the HST astrometry with the ground-based light curve and the derived parallax, we obtain a lens mass of 7.1 +/- 1.3 Msun and a distance of 1.58 +/- 0.18 kpc. We show that the lens emits no detectable light, which, along with having a mass higher than is possible for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper motion for the BH. The proper motion is offset from the mean motion of Galactic-disk stars at similar distances by an amount corresponding to a transverse space velocity of ~45 km/s, suggesting that the BH received a 'natal kick' from its supernova explosion. Previous mass determinations for stellar-mass BHs have come from radial-velocity measurements of Galactic X-ray binaries, and from gravitational radiation emitted by merging BHs in binary systems in external galaxies. Our mass measurement is the first for an isolated stellar-mass BH using any technique., Comment: 37 pages, Published in ApJ

There are currently many large-field surveys operational and planned including the powerful Vera C. Rubin Observatory Legacy Survey of Space and Time. These surveys will increase the number and diversity of transients dramatically. However, for some transients, like supernovae (SNe), we can gain more understanding by directed observations (e.g. shock breakout, $\gamma$-ray detections) than by simply increasing the sample size. For example, the initial emission from these transients can be a powerful probe of these explosions. Upcoming ground-based detectors are not ideally suited to observe the initial emission (shock emergence) of these transients. These observations require a large field-of-view X-ray mission with a UV follow up within the first hour of shock breakout. The emission in the first one hour to even one day provides strong constraints on the stellar radius and asymmetries in the outer layers of stars, the properties of the circumstellar medium (e.g. inhomogeneities in the wind for core-collapse SNe, accreting companion in thermonuclear SNe), and the transition region between these two. This paper describes a simulation for the number of SNe that could be seen by a large field of view lobster eye X-ray and UV observatory., Comment: 13 pages, 7 figures, submitted to ApJ

The immense power of gamma-ray bursts (GRBs) make them ideal probes of the early universe. By using absorption lines in the afterglows of high-redshift GRBs, astronomers can study the evolution of metals in the early universe. With an understanding of the nature of GRB progenitors, the rate and properties of GRBs observed at high redshift can probe the star formation history and the initial mass function of stars at high redshift. This paper presents a detailed study of the metallicity- and mass-dependence of the properties of long-duration GRBs under the black-hole accretion disk paradigm to predict the evolution of these properties with redshift. These models are calibrated on the current GRB observations and then used to make predictions for new observations and new missions (e.g. the proposed Gamow mission) studying high-redshift GRBs., Comment: 21 pages, 20 figures, submitted to ApJ

Building efficient, accurate and generalizable reduced order models of developed turbulence remains a major challenge. This manuscript approaches this problem by developing a hierarchy of parameterized reduced Lagrangian models for turbulent flows, and investigates the effects of enforcing physical structure through Smoothed Particle Hydrodynamics (SPH) versus relying on neural networks (NN)s as universal function approximators. Starting from Neural Network (NN) parameterizations of a Lagrangian acceleration operator, this hierarchy of models gradually incorporates a weakly compressible and parameterized SPH framework, which enforces physical symmetries, such as Galilean, rotational and translational invariances. Within this hierarchy, two new parameterized smoothing kernels are developed in order to increase the flexibility of the learn-able SPH simulators. For each model we experiment with different loss functions which are minimized using gradient based optimization, where efficient computations of gradients are obtained by using Automatic Differentiation (AD) and Sensitivity Analysis (SA). Each model within the hierarchy is trained on two data sets associated with weekly compressible Homogeneous Isotropic Turbulence (HIT): (1) a validation set using weakly compressible SPH; and (2) a high fidelity set from Direct Numerical Simulations (DNS). Numerical evidence shows that encoding more SPH structure improves generalizability to different turbulent Mach numbers and time shifts, and that including the novel parameterized smoothing kernels improves the accuracy of SPH at the resolved scales., Comment: Phys. Rev. Fluids 8, 054602. Published 5 May 2023

The Compton Spectrometer and Imager (COSI) is a 0.2-5 MeV Compton telescope capable of imaging, spectroscopy, and polarimetry of astrophysical sources. Such capabilities are made possible by COSI's germanium cross-strip detectors, which provide high efficiency, high resolution spectroscopy and precise 3D positioning of photon interactions. Science goals for COSI include studies of 0.511 MeV emission from antimatter annihilation in the Galaxy, mapping radioactive elements from nucleosynthesis, determining emission mechanisms and source geometries with polarization, and detecting and localizing multimessenger sources. The instantaneous field of view (FOV) for the germanium detectors is >25% of the sky, and they are surrounded on the sides and bottom by active shields, providing background rejection as well as allowing for detection of gamma-ray bursts or other gamma-ray flares over >50% of the sky. We have completed a Phase A concept study to consider COSI as a Small Explorer (SMEX) satellite mission, and here we discuss the advances COSI-SMEX provides for astrophysics in the MeV bandpass., Comment: 9 pages, 4 figures, 1 table, author affiliations provided on final page. arXiv admin note: substantial text overlap with arXiv:1908.04334

Recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have non-uniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher then the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a 'flashlight' effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of $\zeta \approx 10^{-16} \rm s^{-1}$; however, within $r \approx 0.03$ pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved sub-grid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

Nuclear isomers are populated in the rapid neutron capture process (r process) of nucleosynthesis. The r process may cover a wide range of temperatures, potentially starting from several tens of GK (several MeV) and then cooling as material is ejected from the event. As the r-process environment cools, isomers can freeze out of thermal equilibrium or be directly populated as astrophysical isomers (astromers). Two key behaviors of astromers -- ground state<->isomer transition rates and thermalization temperatures -- are determined by direct transition rates between pairs of nuclear states. We perform a sensitivity study to constrain the effects of unknown transitions on astromer behavior. We also introduce a categorization of astromers that describes their potential effects in hot environments. We provide a table of neutron-rich isomers that includes the astromer type, thermalization temperature, and key unmeasured transition rates., Comment: 52 pages, 13 figures, 2 tables

There is a growing set of observational data demonstrating that cosmic rays exhibit small-scale anisotropies (5-30 deg) with amplitude deviations lying between 0.01-0.1 percent that of the average cosmic ray flux. A broad range of models have been proposed to explain these anisotropies ranging from finite-scale magnetic field structures to dark matter annihilation. The standard diffusion transport methods used in cosmic ray propagation do not capture the transport physics in a medium with finite-scale or coherent magnetic field structures. Here, we present a Monte Carlo transport method, applying it to a series of finite-scale magnetic field structures to determine the requirements of such fields in explaining the observed cosmic ray,small-scale anisotropies., Comment: 14 pages, 13 figures

The breakout of a supernova blast wave from its progenitor star provides strong constraints on the star and its immediate surroundings. These surroundings are shaped by mass loss from the star and can include a wide variety of inhomogeneities. Here we present results of multi-dimensional radiation-hydrodynamics calculations of the interactions of the supernova blast wave with inhomogeneities in the immediate surroundings of a massive Wolf-Rayet star, calculating the effect these interactions have on the shock breakout signal from supernovae., Comment: 14 pages, 12 figures, submitted to ApJ

Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. Post-bounce temperature and density evolution of ZAMS 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid (NuGrid) nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio-isotopes that are of particular interest to the next generation of gamma-ray astronomical observations; 43K, 47Ca, 44Sc, 47Sc, 48V, 48Cr, 51Cr, 52Mn, 59Fe, 56Co, 57Co, 57Ni. These nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43K, 47Ca, 44Sc, 47Sc, and 59Fe, those that are dependent on the progenitor structure are 48V, 51Cr, and 57Co, and those that probe neither are 48Cr, 52Mn, 57Ni, and 56Co. We discuss prospects of observing these radionuclides in supernova remnants., Comment: 17 pages, 19 figures, submitted to ApJ

One of the difficulties in developing accurate numerical models of radiation flow in a coupled radiation-hydrodynamics setting is accurately modeling the transmission across a boundary layer. The COAX experiment is a platform design to test this transmission including standard radiograph and flux diagnostics as well as a temperature diagnostic measuring the population of excitation levels and ionization states of a dopant embedded within the target material. Using a broad range of simulations, we study the experimental errors in this temperature diagnostic. We conclude with proposed physics experiments that show features that are much stronger than the experimental errors and provide the means to study transport models., Comment: 19 pages, 19 figures, accepted by HEDP

Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fundamental input for chemical evolution codes, are usually taken from theoretical nucleosynthesis calculations performed by the nuclear astrophysics community. Therefore, almost all chemical evolution predictions rely on the nuclear physics behind those calculations. In this proceedings, we highlight the impact of nuclear physics uncertainties on galactic chemical evolution predictions. We demonstrate that nuclear physics and galactic evolution uncertainties both have a significant impact on interpreting the origin of neutron-capture elements in our Solar System. Those results serve as a motivation to create and maintain collaborations between the fields of nuclear astrophysics and galaxy evolution., Comment: 8 pages, 2 figures, NPA-IX proceedings

We present a Chandra and XMM-Newton imaging and spectroscopic study of the supernova remnant (SNR) RCW 103 (G332.4-00.4) containing the Central Compact Object 1E 161348-5055. The high resolution Chandra X-ray images reveal enhanced emission in the south-eastern and north-western regions. Equivalent width line images of Fe L, Mg, Si, and S using XMM-Newton data were used to map the distribution of ejecta. The SNR was sectioned into 56 regions best characterized by two-component thermal models. The harder component (kT~0.6 keV) is adequately fitted by the VPSHOCK non-equilibrium ionization model with an ionization timescale ~ 1E11-1E12 cm^-3 s, and slightly enhanced abundances over solar values. The soft component (kT~0.2 keV), fitted by the APEC model, is well described by plasma in collisional ionization equilibrium with abundances consistent with solar values. Assuming a distance of 3.1 kpc and a Sedov phase of expansion into a uniform medium, we estimate an SNR age of 4.4 kyr, a swept-up mass M_sw=16 f_s^{-1/2} D_{3.1}^{5/2} solar masses, and a low explosion energy E=3.7E49 f_s^{-1/2} D_{3.1}^{5/2} erg. This energy could be an order of magnitude higher if we relax the Sedov assumption, the plasma has a low filling factor, the plasma temperature is under-estimated, or if the SNR is expanding into the progenitor's wind-blown bubble. Standard explosion models did not match the ejecta yields. By comparing the fitted abundances to the most recent core-collapse nucleosynthesis models, our best estimate yields a low-mass progenitor around 12-13 solar masses, lower than previously reported. We discuss degeneracies in the model fitting, particularly the effect of altering the explosion energy on the progenitor mass estimate., Comment: 21 pages, 11 figures, 6 tables; to appear in MNRAS

In this Astro2020 APC White Paper, we describe a Small Explorer (SMEX) mission concept called the Compton Spectrometer and Imager. COSI is a Compton telescope that covers the bandpass often referred to as the "MeV Gap" because it is the least explored region of the whole electromagnetic spectrum. COSI provides a significant improvement in sensitivity along with high-resolution spectroscopy, enabling studies of 511 keV electron-positron annihilation emission and measurements of several radioactive elements that trace the Galactic history of supernovae. COSI also measures polarization of gamma-ray bursts (GRBs), accreting black holes, and pulsars as well as detecting and localizing multimessenger sources. In the following, we describe the COSI science, the instrument, and its capabilities. We highlight many Astro2020 science WPs that describe the COSI science in depth., Comment: 11 pages, submitted to the Decadal Survey on Astronomy and Astrophysics (Astro2020) as an Activities, Projects, or State of the Profession Consideration (APC) White Paper

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band., Comment: Astro2020 APC White Paper Updated to make small change to author list in metadata

Black holes in binary star systems are vital for understanding the process of pr oducing gravitational wave sources, understanding how supernovae work, and for p roviding fossil evidence for the high mass stars from earlier in the Universe. At the present time, sample sizes of these objects, and especially of black hole s in binaries, are quite limited. Furthermore, more precise measurements of the binary parameters are needed, as well. With improvements primarily in X-ray an d radio astronomy capabilities, it should be possible to build much larger sampl es of much better measured black hole binaries., Comment: white paper submitted to Astro2020 Decadal Survey

Our understanding of the engines and progenitors of gamma-ray bursts has expanded through the ages as a broader set of diagnostics has allowed us to test our understanding of these objects. Here we review the history of the growth in our understanding, focusing on 3 leading engines and 9 potential progenitors. The first gravitational wave detection of a short burst is the latest in a series of breakthrough observations shaping this understanding and we study the importance of multi-diagnostic, multi-messenger observations on these engines and their progenitors. Our understanding based on a detailed study of nearby bursts can be applied to make predictions for the trends expected as we begin to observe high redshift bursts and we discuss these trends., Comment: 19 pages, 10 figures

We explore the impact of pulsar electromagnetic dipole and fallback accretion emission on the luminosity of a suite of kilonova models. The pulsar models are varied over pulsar magnetic field strength, pulsar lifetime, ejecta mass, and elemental abundances; the fallback models are varied over fallback accretion rate and ejecta mass. For the abundances, we use Fe and Nd as representatives of the wind and dynamical ejecta, respectively. We simulate radiative transfer in the ejecta in either 1D spherical or 2D cylindrical spatial geometry. For the grid of 1D simulations, the mass fraction of Nd is 0, $10^{-4}$, or $10^{-3}$ and the rest is Fe. Our models that fit the bolometric luminosity of AT 2017gfo (the kilonova associated with the first neutron star merger discovered in gravitational waves, GW170817) do not simultaneously fit the B, V, and I time evolution. However, we find that the trends of the evolution in B and V magnitudes are better matched by the fallback model relative to the pulsar model, implying the time dependence of the remnant source influences the color evolution. Further exploration of the parameter space and model deficiencies is needed before we can describe AT 2017gfo with a remnant source., Comment: 15 pages, 10 figures, accepted to ApJ. Simulation data available (under "Source term study") at: https://ccsweb.lanl.gov/astro/transient/transients_astro.html

Galactic binaries with orbital periods less than $\approx$1 hr are strong gravitational wave sources in the mHz regime, ideal for the Laser Interferometer Space Antenna (LISA). In fact, theory predicts that \emph{LISA} will resolve tens of thousands of Galactic binaries individually with a large fraction being bright enough for electromagnetic observations. This opens up a new window where we can study a statistical sample of compact Galactic binaries in both, the electromagnetic as well the gravitational wavebands. Using multi-messenger observations we can measure tidal effects in detached double WD systems, which strongly impact the outcome of WD mergers. For accreting WDs as well as NS binaries, multi-messenger observations give us the possibility to study the angular momentum transport due to mass transfer. In this white paper we present an overview of the opportunities for research on Galactic binaries using multi-messenger observations and summarize some recommendations for the 2020 time-frame., Comment: 5 pages, 1 figure, submitted to the Astro2020 call for science white papers