Coalescence of neutron stars (NSs) gives rise to kilonova, thermal emission powered by radioactive decays of freshly synthesized r-process nuclei. Although observational properties are largely affected by bound–bound opacities of r-process elements, available atomic data have been limited. In this paper, we study element-to-element variation of the opacities in the ejecta of NS mergers by performing systematic atomic structure calculations of r-process elements for the first time. We show that the distributions of energy levels tend to be higher as electron occupation increases for each electron shell due to the larger energy spacing caused by larger effects of spin–orbit and electron–electron interactions. As a result, elements with a fewer number of electrons in the outermost shells tend to give larger contributions to the bound–bound opacities. This implies that Fe is not representative for the opacities of light r-process elements. The average opacities for the mixture of r-process elements are found to be κ ∼ 20–30 cm2 g−1 for the electron fraction of Ye ≤ 0.20, κ ∼ 3–5 cm2 g−1 for Ye = 0.25–0.35, and κ ∼ 1 cm2 g−1 for Ye = 0.40 at $T = 5000\!-\!10\, 000$ K, and they steeply decrease at lower temperature. We show that, even with the same abundance or Ye, the opacity in the ejecta changes with time by one order of magnitude from 1 to 10 d after the merger. Our radiative transfer simulations with the new opacity data confirm that ejecta with a high electron fraction (Ye ≳ 0.25, with no lanthanide) are needed to explain the early, blue emission in GW170817/AT2017gfo while lanthanide-rich ejecta (with a mass fraction of lanthanides ∼5 × 10−3) reproduce the long-lasting near-infrared emission.

No abstract

Ejected material from neutron star mergers give rise to electromagnetic emission powered by radioactive decays of r-process nuclei, which is so called kilonova or macronova. While properties of the emission are largely affected by opacities in the ejected material, available atomic data for r-process elements are still limited. We perform atomic structure calculations for r-process elements: Se (Z = 34), Ru (Z = 44), Te (Z = 52), Ba (Z = 56), Nd (Z = 60), and Er (Z = 68). We confirm that the opacities from bound-bound transitions of open f-shell, Lanthanide elements (Nd and Er) are higher than those of the other elements over a wide wavelength range. The opacities of open s-shell (Ba), p-shell (Se and Te), and d-shell (Ru) elements are lower than those of open f-shell elements and their transitions are concentrated in the ultraviolet wavelengths. We show that the optical brightness can be different by > 2 mag depending on the element abundances in the ejecta such that post-merger, Lanthanide-free ejecta produce brighter and bluer optical emission. Such blue emission from post-merger ejecta can be observed from the polar directions if the mass of the preceding dynamical ejecta in these regions is small. For the ejecta mass of 0.01 M ⊙ , observed magnitudes of the blue emission will reach 21.0 mag (100 Mpc) and 22.5 mag (200 Mpc) in g and r bands within a few days after the merger, which are detectable with 1m or 2m-class telescopes.

Abstract:We have measured extreme ultraviolet (EUV)

We present radiative transfer simulations for blue kilonovae hours after neutron star (NS) mergers by performing detailed opacity calculations for the first time. We calculate atomic structures and opacities of highly ionized elements (up to the 10th ionization) with atomic number Z = 20–56. We find that the bound–bound transitions of heavy elements are the dominant source of the opacities in the early phase (t < 1 day after the merger) and that the ions with a half-closed electron shell provide the highest contributions. The Planck mean opacity for lanthanide-free ejecta (with electron fraction of Y e = 0.30–0.40) can only reach around at t = 0.1 days, whereas that increases up to at t = 1 day. The spherical ejecta model with an ejecta mass of M ej = 0.05 M ⊙ gives the bolometric luminosity of at t ∼ 0.1 days. We confirm that the existing bolometric and multicolor data of GW170817 can be naturally explained by the purely radioactive model. The expected early UV signals reach 20.5 mag at t ∼ 4.3 hr for sources even at 200 Mpc, which is detectable by the facilities such as Swift and the Ultraviolet Transient Astronomy Satellite (ULTRASAT). The early-phase luminosity is sensitive to the structure of the outer ejecta, as also pointed out by Kasen et al. Therefore, the early UV observations give strong constraints on the structure of the outer ejecta and the presence of a heating source besides r-process nuclei.

Coalescence of binary neutron star give rise to electromagnetic emission, kilonova, powered by radioactive decays of r-process nuclei. Observations of kilonova associated with GW170817 provided unique opportunity to study the heavy element synthesis in the Universe. However, atomic data of r-process elements to decipher the light curves and spectral features of kilonova are not fully constructed yet. In this paper, we perform extended atomic calculations of neodymium (Nd, Z = 60) to study the impact of accuracies in atomic calculations to the astrophysical opacities. By employing multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods, we calculate energy levels and transition data of electric dipole transitions for Nd II, Nd III, and Nd IV ions. Compared with previous calculations, our new results provide better agreement with the experimental data. The accuracy of energy levels was achieved in the present work 10 %, 3 % and 11 % for Nd II, Nd III and Nd IV, respectively, comparing with the NIST database. We confirm that the overall properties of the opacity are not significantly affected by the accuracies of the atomic calculations. The impact to the Planck mean opacity is up to a factor of 1.5, which affects the timescale of kilonova at most 20%. However, we find that the wavelength dependent features in the opacity are affected by the accuracies of the calculations. We emphasize that accurate atomic calculations, in particular for low-lying energy levels, are important to provide predictions of kilonova light curves and spectra.

The observations of GW170817/AT2017gfo have provided us with evidence that binary neutron star mergers are sites of r-process nucleosynthesis. However, the observed signatures in the spectra of GW170817/AT2017gfo have not been fully decoded, especially in the near-infrared (NIR) wavelengths. In this paper, we investigate the kilonova spectra over the entire wavelength range with the aim of elemental identification. We systematically calculate the strength of bound–bound transitions by constructing a hybrid line list that is accurate for important strong transitions and complete for weak transitions. We find that the elements on the left side of the periodic table, such as Ca, Sr, Y, Zr, Ba, La, and Ce, tend to produce prominent absorption lines in the spectra. This is because such elements have a small number of valence electrons and low-lying energy levels, resulting in strong transitions. By performing self-consistent radiative transfer simulations for the entire ejecta, we find that La iii and Ce iii appear in the NIR spectra, which can explain the absorption features at λ ∼ 12000–14000 Å in the spectra of GW170817/AT2017gfo. The mass fractions of La and Ce are estimated to be >2 × 10−6 and ∼(1–100) × 10−5, respectively. An actinide element Th can also be a source of absorption as the atomic structure is analogous to that of Ce. However, we show that Th iii features are less prominent in the spectra because of the denser energy levels of actinides compared to those of lanthanides.

No abstract