skip to main content


Title: Axisymmetric Radiative Transfer Models of Kilonovae
Abstract

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g., in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide “curtaining,” which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

 
more » « less
NSF-PAR ID:
10361448
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
910
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 116
Size(s):
["Article No. 116"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The electromagnetic emission from the nonrelativistic 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 ofr-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 ofr-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 a factor of 2–4 uncertainties in the inferred ejecta mass based on peak infrared 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.

     
    more » « less
  2. Abstract

    As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first “kilonova” associated with a binary neutron star merger (NSM) suggests that these sites are hosts of the rapid neutron capture (“r”) process. However, it is unknown just how robust elemental production can be in mergers. Identifying signposts of the production of particular nuclei is critical for fully understanding merger-driven heavy-element synthesis. In this study, we investigate the properties of very neutron-rich nuclei for which superheavy elements (Z≥ 104) can be produced in NSMs and whether they can similarly imprint a unique signature on kilonova light-curve evolution. A superheavy-element signature in kilonovae represents a route to establishing a lower limit on heavy-element production in NSMs as well as possibly being the first evidence of superheavy-element synthesis in nature. Favorable NSM conditions yield a mass fraction of superheavy elementsXZ≥104≈ 3 × 10−2at 7.5 hr post-merger. With this mass fraction of superheavy elements, we find that the component of kilonova light curves possibly containing superheavy elements may appear similar to those arising from lanthanide-poor ejecta. Therefore, photometric characterizations of superheavy-element rich kilonova may possibly misidentify them as lanthanide-poor events.

     
    more » « less
  3. null (Ed.)
    ABSTRACT We present and discuss the optical spectrophotometric observations of the nearby (z = 0.087) Type I superluminous supernova (SLSN I) SN 2017gci, whose peak K-corrected absolute magnitude reaches Mg = −21.5 mag. Its photometric and spectroscopic evolution includes features of both slow- and of fast-evolving SLSN I, thus favoring a continuum distribution between the two SLSN-I subclasses. In particular, similarly to other SLSNe I, the multiband light curves (LCs) of SN 2017gci show two re-brightenings at about 103 and 142 d after the maximum light. Interestingly, this broadly agrees with a broad emission feature emerging around 6520 Å after ∼51 d from the maximum light, which is followed by a sharp knee in the LC. If we interpret this feature as Hα, this could support the fact that the bumps are the signature of late interactions of the ejecta with a (hydrogen-rich) circumstellar material. Then we fitted magnetar- and CSM-interaction-powered synthetic LCs on to the bolometric one of SN 2017gci. In the magnetar case, the fit suggests a polar magnetic field Bp ≃ 6 × 1014 G, an initial period of the magnetar Pinitial ≃ 2.8 ms, an ejecta mass $M_{\rm ejecta}\simeq 9\, \mathrm{M}_\odot $ and an ejecta opacity $\kappa \simeq 0.08\, \mathrm{cm}^{2}\, \rm{g}^{-1}$. A CSM-interaction scenario would imply a CSM mass $\simeq 5\, \mathrm{M}_\odot $ and an ejecta mass $\simeq 12\, \mathrm{M}_\odot $. Finally, the nebular spectrum of phase  + 187 d was modeled, deriving a mass of $\sim 10\, {\rm M}_\odot$ for the ejecta. Our models suggest that either a magnetar or CSM interaction might be the power sources for SN 2017gci and that its progenitor was a massive ($40\, {\rm M}_\odot$) star. 
    more » « less
  4. ABSTRACT

    Neutron star merger accretion discs can launch neutron-rich winds of >10−2M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ∼0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $\sim 1\, \mathrm{d}$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($\lesssim 1\, \mathrm{d}$) light curves are largely unaffected.

     
    more » « less
  5. Kilonovae, one source of electromagnetic emission associated with neutron star mergers, are powered by the decay of radioactive isotopes in the neutron-rich merger ejecta. Models for kilonova emission consistent with the electromagnetic counterpart to GW170817 predict characteristic abundance patterns, determined by the relative balance of different types of material in the outflow. Assuming that the observed source is prototypical, this inferred abundance pattern in turn must matchr-process abundances deduced by other means, such as what is observed in the solar system. We report on analysis comparing the input mass-weighted elemental compositions adopted in our radiative transfer simulations to the mass fractions of elements in the Sun, as a practical prototype for the potentially universal abundance signature from neutron star mergers. We characterize the extent to which our parameter inference results depend on our assumed composition for the dynamical and wind ejecta and examine how the new results compare to previous work. We find that a dynamical ejecta composition calculated using the FRDM2012 nuclear mass and FRLDM fission models with extremely neutron-rich ejecta (Ye= 0.035) along with moderately neutron-rich (Ye= 0.27) wind ejecta composition yields a wind-to-dynamical mass ratio ofMw/Md= 0.47, which best matches the observed AT2017gfo kilonova light curves while also producing the best-matching abundance of neutron capture elements in the solar system, though, allowing for systematics, the ratio may be as high as of order unity.

     
    more » « less