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
Impact of Systematic Modeling Uncertainties on Kilonova Property Estimation
Abstract The precise atomic structure and therefore the wavelength-dependent opacities of lanthanides are highly uncertain. This uncertainty introduces systematic errors in modeling transients like kilonovae and estimating key properties such as mass, characteristic velocity, and heavy metal content. Here, we quantify how atomic data from across the literature as well as choices of thermalization efficiency ofr-process radioactive decay heating impact the light curve and spectra of kilonovae. Specifically, we analyze the spectra of a grid of models produced by the radiative transfer codeSedonathat span the expected range of kilonova properties to identify regions with the highest systematic uncertainty. Our findings indicate that differences in atomic data have a substantial impact on estimates of lanthanide mass fraction, spanning approximately 1 order of magnitude for lanthanide-rich ejecta, and demonstrate the difficulty in precisely measuring the lanthanide fraction in lanthanide-poor ejecta. Mass estimates vary typically by 25%–40% for differing atomic data. Similarly, the choice of thermalization efficiency can affect mass estimates by 20%–50%. Observational properties such as color and decay rate are highly model dependent. Velocity estimation, when fitting solely based on the light curve, can have a typical error of ∼100%. Atomic data of lightr-process elements can strongly affect blue emission. Even for well-observed events like GW170817, the total lanthanide production estimated using different atomic data sets can vary by a factor of ∼6.
more »
« less
- PAR ID:
- 10563400
- Publisher / Repository:
- ApJ
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 975
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 213
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract We study the production of very light elements (Z< 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10−6M⊙of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10−6M⊙of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavyr-process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817.more » « less
-
Abstract We present uniform modeling of eight kilonovae, five following short gamma-ray bursts (GRBs; including GRB 170817A) and three following long GRBs. We model their broadband afterglows to determine the relative contributions of afterglow and kilonova emission. We fit the kilonovae using a three-component model inMOSFiT, and report population median ejecta masses for the total, blue (κB = 0.5 cm2g−1), purple (κP = 3 cm2g−1), and red (κR = 10 cm2g−1) components. The kilonova of GW170817 is near the sample median in most derived properties. We investigate trends between the ejecta masses and the isotropic-equivalent and beaming-correctedγ-ray energies (Eγ,iso,Eγ), as well as rest-frame durations (T90,rest). We find long GRB kilonovae have higher median red ejecta masses (Mej,R ≳ 0.05M⊙) compared to on-axis short GRB kilonovae (Mej,R ≲ 0.02M⊙). We also observe a weak scaling between the total and red ejecta masses withEγ,isoandEγ, though a larger sample is needed to establish a significant correlation. These findings imply a connection between merger-driven long GRBs and larger tidal dynamical ejecta masses, which may indicate that their progenitors are asymmetric compact object binaries. We produce representative kilonova light curves, and find that the planned depths and cadences of the Rubin and Roman Observatory surveys will be sufficient for order-of-magnitude constraints onMej,B(and, for Roman,Mej,PandMej,R) of future kilonovae atz ≲ 0.1.more » « less
-
Context.Type Ia supernovae (SNe Ia) are a key probe in modern cosmology, as they can be used to measure luminosity distances at gigaparsec scales. Models of their light curves are used to project heterogeneous observed data onto a common basis for analysis. Aims.The SALT model currently used for SN Ia cosmology describes SNe as having two sources of variability, accounted for by a color parameterc, and a “stretch” parameterx1. We extend the model to include an additional parameter we labelx2, to investigate the cosmological impact of currently unaddressed light-curve variability. Methods.We constructed a new SALT model, that we dub “SALT3+”. This model was trained by an improved version of theSALTshakercode, using training data combining a selection of the second data release of cosmological SNe Ia from the Zwicky Transient Facility and the existing SALT3 training compilation. Results.We find additional, coherent variability in supernova light curves beyond SALT3. Most of this variation can be described as phase-dependent variation ing − randr − icolor curves, correlated with a boost in the height of the secondary maximum ini-band. These behaviors correlate with spectral differences, particularly in line velocity. We find that fits with the existing SALT3 model tend to address this excess variation with the color parameter, leading to less informative measurements of supernova color. We find that neglecting the new parameter in light-curve fits leads to a trend in Hubble residuals withx2of 0.039 ± 0.005 mag, representing a potential systematic uncertainty. However, we find no evidence of a bias in current cosmological measurements. Conclusions.We conclude that extended SN Ia light-curve models promise mild improvement in the accuracy of color measurements, and corresponding cosmological precision. However, models with more parameters are unlikely to substantially affect current cosmological results.more » « less
An official website of the United States government

