skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Monte Carlo Radiation Transport for Astrophysical Transients Powered by Circumstellar Interaction
Abstract In this paper, we introduceSuperLite, an open-source Monte Carlo radiation transport code designed to produce synthetic spectra for astrophysical transient phenomena affected by circumstellar interaction.SuperLiteutilizes Monte Carlo methods for semi-implicit, semirelativistic radiation transport in high-velocity shocked outflows, employing multigroup structured opacity calculations. The code enables rapid post-processing of hydrodynamic profiles to generate high-quality spectra that can be compared with observations of transient events, including superluminous supernovae, pulsational pair-instability supernovae, and other peculiar transients. We present the methods employed inSuperLiteand compare the code’s performance to that of other radiative transport codes, such asSuperNuand CMFGEN. We show thatSuperLitehas successfully passed standard Monte Carlo radiation transport tests and can reproduce spectra of typical supernovae of Type Ia, Type IIP, and Type IIn.  more » « less
Award ID(s):
1907617
PAR ID:
10440282
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
953
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 132
Size(s):
Article No. 132
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present the implementation of a two-moment-based general-relativistic multigroup radiation transport module in theGeneral-relativisticmultigridnumerical (Gmunu) code. On top of solving the general-relativistic magnetohydrodynamics and the Einstein equations with conformally flat approximations, the code solves the evolution equations of the zeroth- and first-order moments of the radiations in the Eulerian-frame. An analytic closure relation is used to obtain the higher order moments and close the system. The finite-volume discretization has been adopted for the radiation moments. The advection in spatial space and frequency-space are handled explicitly. In addition, the radiation–matter interaction terms, which are very stiff in the optically thick region, are solved implicitly. The implicit–explicit Runge–Kutta schemes are adopted for time integration. We test the implementation with a number of numerical benchmarks from frequency-integrated to frequency-dependent cases. Furthermore, we also illustrate the astrophysical applications in hot neutron star and core-collapse supernovae modelings, and compare with other neutrino transport codes. 
    more » « less
  2. Abstract We present five far- and near-ultraviolet spectra of the Type II plateau supernova, SN 2022acko, obtained 5, 6, 7, 19, and 21 days after explosion, all observed with the Hubble Space Telescope/Space Telescope Imaging Spectrograph. The first three epochs are earlier than any Type II plateau supernova has been observed in the far-ultraviolet revealing unprecedented characteristics. These three spectra are dominated by strong lines, primarily from metals, which contrasts with the featureless early optical spectra. The flux decreases over the initial time series as the ejecta cool and line blanketing takes effect. We model this unique data set with the non–local thermodynamic equilibrium radiation transport codeCMFGEN, finding a good match to the explosion of a low-mass red supergiant with energyEkin= 6 × 1050erg. With these models we identify, for the first time, the ions that dominate the early ultraviolet spectra. We present optical photometry and spectroscopy, showing that SN 2022acko has a peak absolute magnitude ofV= − 15.4 mag and plateau length of ∼115 days. The spectra closely resemble those of SN 2005cs and SN 2012A. Using the combined optical and ultraviolet spectra, we report the fraction of flux as a function of bluest wavelength on days 5, 7, and 19. We create a spectral time-series of Type II supernovae in the ultraviolet, demonstrating the rapid decline of flux over the first few weeks of evolution. Future observations of Type II supernovae are required to map out the landscape of exploding red supergiants, with and without circumstellar material, which is best revealed in high-quality ultraviolet spectra. 
    more » « less
  3. Abstract We introduceWeakhub, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino–matter, neutrino–neutrino interactions and plasma process, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations ofβ-processes, incorporating various weak corrections and medium modifications due to the nuclear equation of state. Calculations of plasma processes, electron neutrino–antineutrino annihilation, and nuclear de-excitation are also included. We also present the detailed derivations of weak interactions and the coupling to the two-moment based general-relativistic multigroup radiation transport in the general-relativisticmultigridnumerical (Gmunu) code. We compare the neutrino opacity spectra for all interactions and estimate their contributions at hydrodynamical points in core-collapse supernovae and binary neutron star (BNS) postmerger remnants, and predict the effects of improved opacities in comparison to conventional ones for a BNS postmerger at a specific hydrodynamical point. We test the implementation of the conventional set of interactions by comparing it to an open-source neutrino libraryNuLibin a core-collapse supernova simulation. We demonstrate good agreement with discrepancies of less than ∼10% in luminosity for all neutrino species, while also highlighting the reasons contributing to the differences. To compare the advanced interactions to the conventional set in core-collapse supernova modeling, we perform simulations to analyze their impacts on neutrino signatures, hydrodynamical behaviors, and shock dynamics, showing significant deviations. 
    more » « less
  4. Abstract We present the Keck Infrared Transient Survey, a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roman SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. Here, we describe the first release of data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients <17 mag in a red optical band (usually ZTFror ATLASobands); a volume-limited sample including all transients within redshiftz< 0.01 (D≈ 50 Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete toz= 0.005. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events, luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction usingPypeIt, which requires minimal human interaction to ensure reproducibility. 
    more » « less
  5. Abstract It has become increasingly useful to answer questions in gravitational-wave astronomy usingtransdimensionalmodels, where the number of free parameters can be varied depending on the complexity required to fit the data. Given the growing interest in transdimensional inference, we introduce a new package for the Bayesian inference Library (Bilby), calledtBilby. ThetBilbypackage allows users to set up transdimensional inference calculations using the existingBilbyarchitecture with off-the-shelf nested samplers and/or Markov Chain Monte Carlo algorithms. Transdimensional models are particularly helpful when seeking to test theoretically uncertain predictions described by phenomenological models. For example, bursts of gravitational waves can be modeled using a superposition ofNwavelets, whereNis itself a free parameter. Short pulses are modeled with small values ofN, whereas longer, more complicated signals are represented with a large number of wavelets stitched together. Other transdimensional models have been used to describe instrumental noise and the population properties of gravitational-wave sources. We provide a few demonstrations oftBilby, including fitting the gravitational-wave signal GW150914 with a superposition ofNsine-Gaussian wavelets. We outline our plans to further develop thetBilbycode suite for a broader range of transdimensional problems. 
    more » « less