Search for: All records

Abstract We present nucleosynthesis and light-curve predictions for a new site of the rapid neutron capture process (r-process) from magnetar giant flares (GFs). Motivated by observations indicating baryon ejecta from GFs, J. Cehula et al. proposed that mass ejection occurs after a shock is driven into the magnetar crust during the GF. We confirm using nuclear reaction network calculations that these ejecta synthesize moderate yields of third-peakr-process nuclei and more substantial yields of lighterr-nuclei, while leaving a sizable abundance of free neutrons in the outermost fastest expanding ejecta layers. The finalr-process mass fraction and distribution are sensitive to the relative efficiencies ofα-capture andn-capture freeze-outs. We use our nucleosynthesis output in a semianalytic model to predict the light curves of novae breves, the transients following GFs powered by radioactive decay. For a baryonic ejecta mass similar to that inferred of the 2004 Galactic GF from SGR 1806-20, we predict a peak UV/optical luminosity of ∼10³⁹–10⁴⁰erg s⁻¹at ∼10–15 minutes, rendering such events potentially detectable to several Mpc following a gamma-ray trigger by wide-field transient monitors such as ULTRASAT/UVEX. The peak luminosity and timescale of the transient increase with the GF strength due to the larger ejecta mass. Although GFs likely contribute 1%–10% of the total Galacticr-process budget, their short delay-times relative to star formation make them an attractive source to enrich the earliest generations of stars. 

Abstract Relativistic jets from a Kerr black hole (BH) following the core collapse of a massive star (“collapsar”) is a leading model for gamma-ray bursts (GRBs). However, the two key ingredients for a Blandford–Znajek-powered jet—rapid rotation and a strong magnetic field—seem mutually exclusive. Strong fields in the progenitor star’s core transport angular momentum outward more quickly, slowing down the core before collapse. Through innovative multidisciplinary modeling, we first use MESA stellar evolution models followed to core collapse to explicitly show that the small length scale of the instabilities—likely responsible for angular momentum transport in the core (e.g., Tayler–Spruit)—results in a lownetmagnetic flux fed to the BH horizon, far too small to power GRB jets. Instead, we propose a novel scenario in which collapsar BHs acquire their magnetic “hair” from their progenitor proto–neutron star (PNS), which is likely highly magnetized from an internal dynamo. We evaluate the conditions for the BH accretion disk to pin the PNS magnetosphere to its horizon immediately after the collapse. Our results show that the PNS spin-down energy released before collapse matches the kinetic energy of Type Ic-BL supernovae, while the nascent BH’s spin and magnetic flux produce jets consistent with observed GRB characteristics. We map our MESA models to 3D general-relativistic magnetohydrodynamic simulations and confirm that accretion disks confine the strong magnetic flux initiated near a rotating BH, enabling the launch of successful GRB jets, whereas a slower-spinning BH or one without a disk fails to do so. 

Abstract Eruptive mass loss in massive stars is known to occur, but the mechanism(s) are not yet well understood. One proposed physical explanation appeals to opacity-driven super-Eddington luminosities in stellar envelopes. Here, we present a 1D model for eruptive mass loss and implement this model in theMESAstellar evolution code. The model identifies regions in the star where the energy associated with a local super-Eddington luminosity exceeds the binding energy of the overlaying envelope. The material above such regions is ejected from the star. Stars with initial masses of 10−100M_⊙at solar and SMC metallicities are evolved through core helium burning, with and without this new eruptive mass-loss scheme. We find that eruptive mass loss of up to ∼10⁻²M_⊙yr⁻¹can be driven by this mechanism, and occurs in a vertical band on the H-R diagram between $3.5\lesssim \log \left(T_{\mathrm{eff}}/\mathrm{K}\right)\lesssim 4.0$. This predicted eruptive mass loss prevents stars of initial masses ≳20M_⊙from evolving to become red supergiants (RSGs), with the stars instead ending their lives as blue supergiants, and offers a possible explanation for the observed lack of RSGs in that mass regime. 

Abstract Core-collapse supernovae (SNe) are candidate sites for rapid neutron capture process (r-process) nucleosynthesis. We explore the effects of enrichment fromr-process nuclei on the light curves of hydrogen-rich SNe and assess the detectability of these signatures. We modify the radiation hydrodynamics code, SuperNova Explosion Code, to include the approximate effects of opacity and radioactive heating fromr-process elements in the supernova (SN) ejecta. We present models spanning a range of totalr-process massesM_rand their assumed radial distribution within the ejecta, finding thatM_r≳ 10⁻²M_⊙is sufficient to induce appreciable differences in their light curves as compared to ordinary hydrogen-rich SNe (without anyr-process elements). The primary photometric signatures ofr-process enrichment include a shortening of the plateau phase, coinciding with the hydrogen-recombination photosphere retreating to ther-process-enriched layers, and a steeper post-plateau decline associated with a reddening of the SN colors. We compare ourr-process-enriched models to ordinary SNe models and observational data, showing that yields ofM_r≳ 10⁻²M_⊙are potentially detectable across several of the metrics used by transient observers, provided thatr-process-rich layers are mixed at least halfway to the ejecta surface. This detectability threshold can roughly be reproduced analytically using a two-zone (kilonova-within-an-SN) picture. Assuming that a small fraction of SNe produce a detectabler-process yield ofM_r≳ 10⁻²M_⊙, and respecting constraints on the total Galactic production rate, we estimate that ≳10³–10⁴SNe need be observed to find oner-enriched event, a feat that may become possible with the Vera Rubin Observatory. 

Abstract Kilonovae, the ultraviolet/optical/infrared counterparts to binary neutron star mergers, are an exceptionally rare class of transients. Optical follow-up campaigns are plagued by contaminating transients, which may mimic kilonovae but do not receive sufficient observations to measure the full photometric evolution. In this work, we present an analysis of the multiwavelength dataset of supernova (SN) 2025ulz, a proposed kilonova candidate following the low-significance detection of gravitational waves originating from the potential binary neutron star merger S250818k. Despite an early rapid decline in brightness, our multiwavelength observations of SN 2025ulz reveal that it is a type IIb SN. As part of this analysis, we demonstrate the capabilities of a novel quantitative scoring algorithm to determine the likelihood that a transient candidate is a kilonova, based primarily on its three-dimensional location and light-curve evolution. We also apply our scoring algorithm to other transient candidates in the localization volume of S250818k and find that, at all times after the discovery of SN 2025ulz, there are ≥4 candidates with a score comparable to SN 2025ulz, indicating that the kilonova search may have benefited from the additional follow-up of other candidates. During future kilonova searches, this type of scoring algorithm will be useful to rule out contaminating transients in real time, optimizing the use of valuable telescope resources. 

Abstract Gravitational-wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernovae (PISNe). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections of merging BHs. We useMESAto evolve single, nonrotating, massive helium cores with a metallicity ofZ= 10⁻⁵, until they either collapse to form a BH or explode as a PISN, without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the ±3σuncertainty in our high-resolution tabulated¹²C(α,γ)¹⁶O reaction rate probability distribution function. We extensively test temporal and spatial resolutions for resolving the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower-resolution models reach convergence. We establish a new lower edge of the upper mass gap asM_lower≃ ${60}_{-14}^{+32}$M_⊙from the ±3σuncertainty in the¹²C(α,γ)¹⁶O rate. We explore the effect of a larger 3αrate on the lower edge of the upper mass gap, findingM_lower≃ ${69}_{-18}^{+34}$M_⊙. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog. 

Abstract Periodic variables illuminate the physical processes of stars throughout their lifetime. Wide-field surveys continue to increase our discovery rates of periodic variable stars. Automated approaches are essential to identify interesting periodic variable stars for multiwavelength and spectroscopic follow-up. Here we present a novel unsupervised machine-learning approach to hunt for anomalous periodic variables using phase-folded light curves presented in the Zwicky Transient Facility Catalogue of Periodic Variable Stars by Chen et al. We use a convolutional variational autoencoder to learn a low-dimensional latent representation, and we search for anomalies within this latent dimension via an isolation forest. We identify anomalies with irregular variability. Most of the top anomalies are likely highly variable red giants or asymptotic giant branch stars concentrated in the Milky Way galactic disk; a fraction of the identified anomalies are more consistent with young stellar objects. Detailed spectroscopic follow-up observations are encouraged to reveal the nature of these anomalies.