[Paper Review] Chemical Distribution of the Dynamical Ejecta in the Neutron Star Merger GW170817

Shreya Anand, P. T. H. Pang+12 more|arXiv (Cornell University)|Jul 20, 2023

Gamma-ray bursts and supernovaePhysics and Astronomy3 citations

TL;DR

This study uses precise Hubble Space Telescope and VLBI inclination constraints (19–25°) to refine kilonova modeling of GW170817’s AT2017gfo, revealing that the observer's line of sight intersects a lanthanide-poor, polar dynamical ejecta component with an average electron fraction $\bar{Y}_{\rm e,dyn} = 0.22$. The results confirm strong angular dependence in ejecta composition, with lanthanide-rich material confined to a 35° half-opening angle around the equatorial plane, resolving prior tensions with numerical-relativity simulations and supporting robust $r$-process nucleosynthesis up to the third peak.

ABSTRACT

GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{ m e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{ m e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{ m e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{ m e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{ m e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{ m rd}$ peak in the tidal dynamical ejecta.

Motivation & Objective

• To improve inference of ejecta properties in the neutron star merger GW170817 by incorporating precise inclination angle constraints from Hubble and VLBI observations.
• To resolve discrepancies between prior one-dimensional kilonova models and numerical-relativity simulations by employing multi-dimensional, inclination-dependent ejecta modeling.
• To determine whether the observed early blue emission in AT2017gfo can be explained by a non-spherical, compositionally stratified dynamical ejecta without additional powering mechanisms.
• To assess the robustness of $r$-process nucleosynthesis up to the third peak based on inferred electron fraction distributions in the dynamical ejecta.
• To demonstrate that early-time photometry combined with inclination constraints enables reliable inference of chemical composition and ejecta geometry in kilonovae.

Proposed method

• Employed a grid-based, multi-dimensional kilonova modeling framework that incorporates state-of-the-art heating rates, thermalization efficiencies, and frequency-dependent opacities.
• Parameterized the dynamical ejecta using $\bar{Y}_{\rm e,dyn}$, the average electron fraction, with angular dependence assumed based on numerical-relativity simulations: $Y_{\rm e,dyn}$ varies from ~0.17 at the equator to ~0.41 at the pole.
• Used a feed-forward neural network trained on synthetic light curves to enable fast Bayesian parameter estimation across the full parameter space.
• Applied a tight prior on the viewing angle ($19^\circ < \theta_{\rm obs} < 25^\circ$) derived from proper motion measurements of the relativistic jet in AT2017gfo.
• Fitted the model to multi-band photometry of AT2017gfo to infer ejecta masses, velocities, and electron fractions, with constraints on both dynamical and disk-wind components.
• Validated results against numerical-relativity predictions and compared with recent claims of spherical ejecta based on Sr II spectral features.

Experimental results

Research questions

• RQ1 Does the inclusion of precise inclination constraints significantly alter the inferred properties of the dynamical ejecta in GW170817’s kilonova AT2017gfo?
• RQ2 Can the early blue emission in AT2017gfo be explained by a non-spherical, compositionally stratified dynamical ejecta without invoking additional powering mechanisms?
• RQ3 What is the angular dependence of the electron fraction in the dynamical ejecta, and how does it affect the production of $r$-process elements?
• RQ4 How do the inferred ejecta parameters compare with predictions from numerical-relativity simulations, and does the new analysis resolve prior tensions?
• RQ5 To what extent can early-time photometry alone, combined with inclination constraints, be used to infer the chemical composition of kilonova ejecta?

Key findings

• The best-fit model with the tight inclination prior yields $\bar{Y}_{\rm e,dyn} = 0.22$, a significant reduction from the previous estimate of $0.24$ without inclination constraints.
• The dynamical ejecta exhibit strong angular dependence: $Y_{\rm e,dyn} \sim 0.17$ in the equatorial plane (lanthanide-rich, $Y_{\rm e} \lesssim 0.25$) and $Y_{\rm e,dyn} \sim 0.41$ at the poles (lanthanide-poor, $Y_{\rm e} \gtrsim 0.25$), with a half-opening angle of $35^\circ$ for the lanthanide-rich component.
• The observer’s line of sight is aligned with the lanthanide-poor, polar dynamical ejecta, which explains the early blue emission in AT2017gfo without requiring additional powering sources.
• The inferred ejecta masses and velocities—$M_{\rm ej,dyn} = 1.3 \times 10^{-3}\,M_\odot$, $\bar{v}_{\rm ej,dyn} = 0.14\,c$—are in good agreement with numerical-relativity predictions, resolving prior tensions.
• The range of electron fractions, extending down to $Y_{\rm e,dyn} \sim 0.17$, enables robust production of $r$-process nuclei up to the third peak, consistent with nucleosynthetic models.
• The results contradict a recent claim of spherical ejecta with minimal $Y_{\rm e}$ variation, instead supporting a highly structured, anisotropic composition in the dynamical ejecta.