Search for: All records

Formed in the aftermath of a core-collapse supernova or neutron star merger, a hot proto–neutron star (PNS) launches an outflow driven by neutrino heating lasting for up to tens of seconds. Though such winds are considered potential sites for the nucleosynthesis of heavy elements via the rapid neutron capture process (r-process), previous work has shown that unmagnetized PNS winds fail to achieve the necessary combination of high entropy and/or short dynamical timescale in the seed nucleus formation region. We present three-dimensional general-relativistic magnetohydrodynamical simulations of PNS winds which include the effects of a dynamically strong (B≳ 10¹⁵G) dipole magnetic field. After initializing the magnetic field, the wind quickly develops a helmet-streamer configuration, characterized by outflows along open polar magnetic field lines and a “closed” zone of trapped plasma at lower latitudes. Neutrino heating within the closed zone causes the thermal pressure of the trapped material to rise in time compared to the polar outflow regions, ultimately leading to the expulsion of this matter from the closed zone on a timescale of ∼60 ms, consistent with the predictions of Thompson. The high entropies of these transient ejecta are still growing at the end of our simulations and are sufficient tomore »enable a successful second-peakr-process in at least a modest ≳1% of the equatorial wind ejecta.

AT 2022cmc is a luminous optical transient (νLν ≳ 1045 erg s−1) accompanied by decaying non-thermal X-rays (peak duration tX ≲ days and isotropic energy EX,iso ≳ 1053 erg) and a long-lived radio/mm synchrotron afterglow, which has been interpreted as a jetted tidal disruption event (TDE). Both an equipartition analysis and a detailed afterglow model reveal the radio/mm emitting plasma to be expanding mildly relativistically (Lorentz factor $\Gamma \gtrsim \, \mathrm{ few}$ ) with an opening angle θj ≃ 0.1 and roughly fixed energy Ej,iso ≳ few × 1053 erg into an external medium of density profile n ∝ R−k with k ≃ 1.5–2, broadly similar to that of the first jetted TDE candidate Swift J1644+57 and consistent with Bondi accretion at a rate of ∼$10^{-3}\,\dot{M}_{\rm Edd}$ on to a 106 M⊙ black hole before the outburst. The rapidly decaying optical emission over the first days is consistent with fast-cooling synchrotron radiation from the same forward shock as the radio/mm emission, while the bluer slowly decaying phase to follow likely represents a separate thermal emission component. Emission from the reverse shock may have peaked during the first days, but its non-detection in the optical band places an upper bound Γj ≲ 100 on themore »Lorentz factor of the unshocked jet. Although a TDE origin for AT 2022cmc is indeed supported by some observations, the vast difference between the short-lived jet activity phase tX ≲ days and the months-long thermal optical emission also challenges this scenario. A stellar core-collapse event giving birth to a magnetar or black hole engine of peak duration ∼1 d offers an alternative model also consistent with the circumburst environment, if interpreted as a massive star wind.

Gamma-ray bursts (GRBs) have historically been divided into two classes. Short-duration GRBs are associated with binary neutron star mergers (NSMs), while long-duration bursts are connected to a subset of core-collapse supernovae (SNe). GRB 211211A recently made headlines as the first long-duration burst purportedly generated by an NSM. The evidence for an NSM origin was excess optical and near-infrared emission consistent with the kilonova observed after the gravitational-wave-detected NSM GW170817. Kilonovae derive their unique electromagnetic signatures from the properties of the heavy elements synthesized by rapid neutron capture (ther-process) following the merger. Recent simulations suggest that the “collapsar” SNe that trigger long GRBs may also producer-process elements. While observations of GRB 211211A and its afterglow rule out an SN typical of those that follow long GRBs, an unusual collapsar could explain both the duration of GRB 211211A and ther-process-powered excess in its afterglow. We use semianalytic radiation transport modeling to evaluate low-mass collapsars as the progenitors of GRB 211211A–like events. We compare a suite of collapsar models to the afterglow-subtracted emission that followed GRB 211211A, and find the best agreement for models with high kinetic energies and an unexpected pattern of⁵⁶Ni enrichment. We discuss how core-collapse explosions could produce suchmore »ejecta, and how distinct our predictions are from those generated by more straightforward kilonova models. We also show that radio observations can distinguish between kilonovae and the more massive collapsar ejecta we consider here.

Luminous fast blue optical transients (LFBOTs) such as AT2018cow form a rare class of engine-powered explosions of uncertain origin. A hallmark feature of these events is radio/millimeter synchrotron emission powered by the interaction of fast ≳0.1cejecta and dense circumstellar material (CSM) extending to large radii ≳10¹⁶cm surrounding the progenitor. Assuming this CSM to be an outflow from the progenitor, we show that dust grains up to ∼1μm in size can form in the outflow in the years before the explosion. This dusty CSM would attenuate the transient’s ultraviolet emission prior to peak light, before being destroyed by the rising luminosity, reddening the premaximum colors (consistent with the premaximum red-to-blue color evolution of the LFBOT candidate MUSSES2020J). Reradiation by the dust before being destroyed generates a near-infrared (NIR) “echo” of luminosity ∼10⁴¹–10⁴²erg s⁻¹lasting weeks, which is detectable over the transient’s rapidly fading blue continuum. We show that this dust echo is compatible with the previously unexplained NIR excess observed in AT2018cow. The gradual decay of the early NIR light curve can result from CSM, which is concentrated in a wide-angle equatorial outflow or torus, consistent with the highly aspherical geometry of AT2018cow’s ejecta. Premaximum optical/UV and NIR follow-up of LFBOTsmore »provide a new probe of their CSM environments and place additional constraints on their progenitors.

In 2021 August, the Fermi Large Area Telescope, H.E.S.S., and MAGIC detected GeV and TeVγ-ray emission from an outburst of recurrent nova RS Ophiuchi. This detection represents the first very high-energyγ-rays observed from a nova, and it opens a new window to study particle acceleration. Both H.E.S.S. and MAGIC described the observedγ-rays as arising from a single, external shock. In this paper, we perform detailed, multi-zone modeling of RS Ophiuchi’s 2021 outburst, including a self-consistent prescription for particle acceleration and magnetic field amplification. We demonstrate that, contrary to previous work, a single shock cannot simultaneously explain RS Ophiuchi’s GeV and TeV emission, in particular the spectral shape and distinct light-curve peaks. Instead, we put forward a model involving multiple shocks that reproduces the observedγ-ray spectrum and temporal evolution. The simultaneous appearance of multiple distinct velocity components in the nova optical spectrum over the first several days of the outburst supports the presence of distinct shocks, which may arise either from the strong latitudinal dependence of the density of the external circumbinary medium (e.g., in the binary equatorial plane versus the poles) or due to internal collisions within the white dwarf ejecta (which power theγ-ray emission in classical novae).

Despite recent progress, the astrophysical channels responsible for rapid neutron capture (r-process) nucleosynthesis remain an unsettled question. Observations of the kilonova following the gravitational-wave-detected neutron star merger GW170817 established mergers as one site of ther-process, but additional sources may be needed to fully explainr-process enrichment in the universe. One intriguing possibility is that rapidly rotating massive stars undergoing core collapse launchr-process-rich outflows off the accretion disks formed from their infalling matter. In this scenario,r-process winds are one component of the supernova (SN) ejecta produced by “collapsar” explosions. We present the first systematic study of the effects ofr-process enrichment on the emission from collapsar-generated SNe. We semianalytically modelr-process SN emission from explosion out to late times and determine its distinguishing features. The ease with whichr-process SNe can be identified depends on how effectively wind material mixes into the initiallyr-process-free outer layers of the ejecta. In many cases, enrichment produces a near-infrared (NIR) excess that can be detected within ∼75 days of explosion. We also discuss optimal targets and observing strategies for testing ther-process collapsar theory, and find that frequent monitoring of optical and NIR emission from high-velocity SNe in the first few months after explosion offers a reasonable chance ofmore »success while respecting finite observing resources. Such early identification ofr-process collapsar candidates also lays the foundation for nebular-phase spectroscopic follow-up in the NIR and mid-infrared, for example, with the James Webb Space Telescope.

We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDEs). Motivated by recent hydrodynamical simulations, we assume that the debris streams promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radiusR_v∼ 10¹⁴cm (photosphere radius ∼10¹⁵cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin–Helmholtz contractionR_v∝t⁻¹, its temperature risingT_eff∝t^1/2while its total luminosity remains roughly constant; the optical luminosity decays as$\nu L_{\nu }\propto R_v^2{T}_{\mathrm{eff}}\propto t^{-3/2}$. Despite this similarity to the mass fallback rate${\dot{M}}_{\mathrm{fb}}\propto t^{-5/3}$, envelope heating from fallback accretion is subdominant compared to the envelope cooling luminosity except near optical peak (where they are comparable). Envelope contraction can be delayed by energy injection from accretion from the inner envelope onto the SMBH in a regulated manner, leading to a late-time flattening of the optical/X-ray light curves, similar to those observed in some TDEs. Eventually, as the envelope contracts to near the circularization radius, the SMBH accretion rate rises to its maximum, in tandem with the decreasing optical luminosity. This cooling-induced (rather than circularization-induced) delay of up to several hundred days may account for themore »delayed onset of thermal X-rays, late-time radio flares, and high-energy neutrino generation, observed in some TDEs. We compare the model predictions to recent TDE light-curve correlation studies, finding both agreement and points of tension.

The process of unstable mass transfer in a stellar binary can result in either a complete merger of the stars or successful removal of the donor envelope leaving a surviving more compact binary. Luminous red novae (LRNe) are the class of optical transients believed to accompany such merger/common envelope events. Past works typically model LRNe using analytic formulae for supernova light curves that make assumptions (e.g., radiation-dominated ejecta, neglect of hydrogen recombination energy) not justified in stellar mergers due to the lower velocities and specific thermal energy of the ejecta. We present a one-dimensional model of LRN light curves that accounts for these effects. Consistent with observations, we find that LRNe typically possess two light-curve peaks, an early phase powered by initial thermal energy of the hot, fastest ejecta layers and a later peak powered by hydrogen recombination from the bulk of the ejecta. We apply our model to a sample of LRNe to infer their ejecta properties (mass, velocity, and launching radius) and compare them to the progenitor donor star properties from pretransient imaging. We define the maximum luminosity achievable for a given donor star in the limit that the entire envelope is ejected, finding that several LRNemore »violate this limit. Shock interaction between the ejecta and predynamical mass loss may provide an additional luminosity source to alleviate this tension. Our model can also be applied to the merger of planets with stars or stars with compact objects.

Understanding the radii of massive stars throughout their evolution is important to answering numerous questions about stellar physics, from binary interactions on the main sequence to the pre-supernova radii. One important factor determining a star’s radius is the fraction of its mass in elements heavier than Helium (metallicity, Z). However, the metallicity enters stellar evolution through several distinct microphysical processes, and which dominates can change throughout stellar evolution and with the overall magnitude of Z. We perform a series of numerical experiments with 15 $\, \mathrm{M}_{\odot }$mesa models computed doubling separately the metallicity entering the radiative opacity, the equation of state, and the nuclear reaction network to isolate the impact of each on stellar radii. We explore separately models centred around two metallicity values: one near solar Z = 0.02 and another sub-solar Z ∼ 10−3, and consider several key epochs from the end of the main sequence to core carbon depletion. We find that the metallicity entering the opacity dominates at most epochs for the solar metallicity models, contributing to on average ∼60–90 per cent of the total change in stellar radius. Nuclear reactions have a larger impact (∼50–70 per cent) during most epochs in the subsolar Z models. The methodology introduced heremore »can be employed more generally to propagate known microphysics errors into uncertainties on macrophysical observables including stellar radii.

A growing number of core-collapse supernovae (SNe) that show evidence for interaction with dense circumstellar medium (CSM) are accompanied by “precursor” optical emission rising weeks to months prior to the explosion. The precursor luminosities greatly exceed the Eddington limit of the progenitor star, implying that they are accompanied by substantial mass loss. Here, we present a semi-analytic model for SN precursor light curves, which we apply to constrain the properties and mechanisms of the pre-explosion mass loss. We explore two limiting mass-loss scenarios: (1) an “eruption” arising from shock breakout following impulsive energy deposition below the stellar surface; and (2) a steady “wind,” due to sustained heating of the progenitor envelope. The eruption model, which resembles a scaled-down version of Type IIP SNe, can explain the luminosities and timescales of well-sampled precursors, for ejecta masses ∼ 0.1–1M_⊙and velocities ∼ 100–1000 km s⁻¹. By contrast, the steady wind scenario cannot explain the highest precursor luminosities ≳ 10⁴¹erg s⁻¹, under the constraint that the total ejecta mass does not exceed the entire progenitor mass (though the less luminous SN 2020tlf precursor can be explained by a mass-loss rate ∼ 1M_⊙yr⁻¹). However, shock interaction between the wind and pre-existing (earlier ejected) CSMmore »may boost its radiative efficiency and mitigate this constraint. In both the eruption and wind scenarios, the precursor ejecta forms compact (≲10¹⁵cm) optically thick CSM at the time of core collapse; though only directly observable via rapid post-explosion spectroscopy (≲ a few days before being overtaken by the SN ejecta), this material can boost the SN luminosity via shock interaction.