No Arabic abstract
In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.
Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on the order of 15 to 250 days. This is due to the anomalously long half-life of this isotope and the efficiency of fission thermalization compared to other nuclear channels. We estimate the fission fragment yield of this nucleus and outline the astrophysical conditions under which $^{254}$Cf has the greatest impact to the light curve. Future observations in the middle-IR which are bright during this regime could indicate the production of actinide nucleosynthesis.
During the second observing run of the Laser Interferometer gravitational- wave Observatory (LIGO) and Virgo Interferometer, a gravitational-wave signal consistent with a binary neutron star coalescence was detected on 2017 August 17th (GW170817), quickly followed by a coincident short gamma-ray burst trigger by the Fermi satellite. The Distance Less Than 40 (DLT40) Mpc supernova search performed pointed follow-up observations of a sample of galaxies regularly monitored by the survey which fell within the combined LIGO+Virgo localization region, and the larger Fermi gamma ray burst error box. Here we report the discovery of a new optical transient (DLT17ck, also known as SSS17a; it has also been registered as AT 2017gfo) spatially and temporally coincident with GW170817. The photometric and spectroscopic evolution of DLT17ck are unique, with an absolute peak magnitude of Mr = -15.8 pm 0.1 and an r-band decline rate of 1.1mag/d. This fast evolution is generically consistent with kilonova models, which have been predicted as the optical counterpart to binary neutron star coalescences. Analysis of archival DLT40 data do not show any sign of transient activity at the location of DLT17ck down to r~19 mag in the time period between 8 months and 21 days prior to GW170817. This discovery represents the beginning of a new era for multi-messenger astronomy opening a new path to study and understand binary neutron star coalescences, short gamma-ray bursts and their optical counterparts.
The neutron star (NS) merger GW170817 was followed over several days by optical-wavelength (blue) kilonova (KN) emission likely powered by the radioactive decay of light r-process nuclei synthesized by ejecta with a low neutron abundance (electron fraction Ye ~ 0.25-0.35). While the composition and high velocities of the blue KN ejecta are consistent with shock-heated dynamical material, the large quantity is in tension with the results of numerical simulations. We propose an alternative ejecta source: the neutrino-heated, magnetically-accelerated wind from the strongly-magnetized hypermassive NS (HMNS) remnant. A rapidly-spinning HMNS with an ordered surface magnetic field of strength B ~ 1-3e14 G and lifetime t_rem ~ 0.1-1 s can simultaneously explain the velocity, total mass, and electron fraction of the blue KN ejecta. The inferred HMNS lifetime is close to its Alfven crossing time, suggesting global magnetic torques could be responsible for bringing the HMNS into solid body rotation and instigating its gravitational collapse. Different origins for the KN ejecta may be distinguished by their predictions for the emission in the first hours after the merger, when the luminosity is enhanced by heating from internal shocks; the latter are likely generic to any temporally-extended ejecta source (e.g. magnetar or accretion disk wind) and are not unique to the emergence of a relativistic jet. The same shocks could mix and homogenizes the composition to a low but non-zero lanthanide mass fraction, X_La ~ 1e-3, as advocated by some authors, but only if the mixing occurs after neutrons are consumed in the r-process on a timescale >~ 1 s.
We present a rapid analytic framework for predicting kilonova light curves following neutron star (NS) mergers, where the main input parameters are binary-based properties measurable by gravitational wave detectors (chirp mass and mass ratio, orbital inclination) and properties dependent on the nuclear equation of state (tidal deformability, maximum NS mass). This enables synthesis of a kilonova sample for any NS source population, or determination of the observing depth needed to detect a live kilonova given gravitational wave source parameters in low latency. We validate this code, implemented in the public MOSFiT package, by fitting it to GW170817. A Bayes factor analysis overwhelmingly ($B>10^{10}$) favours the inclusion of an additional luminosity source in addition to lanthanide-poor dynamical ejecta during the first day. This is well fit by a shock-heated cocoon model, though differences in the ejecta structure, opacity or nuclear heating rate cannot be ruled out as alternatives. The emission thereafter is dominated by a lanthanide-rich viscous wind. We find the mass ratio of the binary is $q=0.92pm0.07$ (90% credible interval). We place tight constraints on the maximum stable NS mass, $M_{rm TOV}=2.17^{+0.08}_{-0.11}$ M$_odot$. For a uniform prior in tidal deformability, the radius of a 1.4 M$_odot$ NS is $R_{1.4}sim 10.7$ km. Re-weighting with a prior based on equations of state that support our credible range in $M_{rm TOV}$, we derive a final measurement $R_{1.4}=11.06^{+1.01}_{-0.98}$ km. Applying our code to the second gravitationally-detected neutron star merger, GW190425, we estimate that an associated kilonova would have been fainter (by $sim0.7$ mag at one day post-merger) and declined faster than GW170817, underlining the importance of tuning follow-up strategies individually for each GW-detected NS merger.
We present the first effort to aggregate, homogenize, and uniformly model the combined ultraviolet, optical, and near-infrared dataset for the electromagnetic counterpart of the binary neutron star merger GW170817. By assembling all of the available data from 18 different papers and 46 different instruments, we are able to identify and mitigate systematic offsets between individual datasets, and to identify clear outlying measurements, with the resulting pruned and adjusted dataset offering an opportunity to expand the study of the kilonova. The unified dataset includes 647 individual flux measurements, spanning 0.45 to 29.4 days post-merger, and thus has greater constraining power for physical models than any single dataset. We test a number of semi-analytical models and find that the data are well modeled with a three-component kilonova model: a blue lanthanide-poor component with Mej~0.020 Msol and vej~0.27c; an intermediate opacity purple component with Mej~0.047 Msol and vej~0.15c; and a red lanthanide-rich component with Mej~0.011 Msol and vej~0.14c. We further explore the possibility of ejecta asymmetry and its impact on the estimated parameters. From the inferred parameters we draw conclusions about the physical mechanisms responsible for the various ejecta components, the properties of the neutron stars, and, combined with an up-to-date merger rate, the implications for r-process enrichment via this channel. To facilitate future studies of this keystone event we make the unified dataset and our modeling code public.