Comparative study of kilonova opacities for three elements of the sixth period (hafnium, osmium, and gold) from new atomic structure calculations in Hf I–IV, Os I–IV, and Au I–IV^★

¹ Physique Atomique et Astrophysique, Université de Mons, 7000 Mons, Belgium
e-mail: sirine.bennasr@umons.ac.be; sirinamarouan@gmail.com; pascal.quinet@umons.ac.be
² Institut d’Astronomie et d’Astrophysique, Université Libre de Bruxelles, CP 226, 1050 Brussels, Belgium
³ IPNAS, Université de Liège, Sart Tilman, 4000 Liège, Belgium

Received: 12 December 2023
Accepted: 15 February 2024

Abstract

Aims. It is now well established that a large amount of heavy (trans-iron) elements are produced during neutron star (NS) mergers. These elements can be detected in the spectra of the kilonova emitted from the post-merger ejected materials. Due to the high level densities that characterize the complex configurations belonging to heavy elements, thus giving rise to millions of absorption lines, the kilonova ejecta opacity is of significant importance. The elements that contribute the most to the latter are those with an unfilled nd subshell belonging to the fifth and the sixth rows of the periodic table, and those with an unfilled nf subshell belonging to the lanthanide and actinide groups. The aim of the present work is to make a new contribution to this field by performing large-scale atomic structure calculations in three specific sixth-row 5d elements, namely hafnium, osmium, and gold, in the first four charge stages (I–IV), and by computing the corresponding opacities, while focusing on the importance of the atomic models used.

Methods. The pseudo-relativistic Hartree–Fock (HFR) method, including extended sets of interacting configurations, was used for the atomic structure and radiative parameter calculations, while the expansion formalism was used to estimate the opacities.

Results. Theoretical energy levels, wavelengths, and oscillator strengths were computed for millions of spectral lines in Hf I–IV, Os I–IV, and Au I–IV ions, the reliability of these parameters being assessed through detailed comparisons with previously published experimental and theoretical results. The newly obtained atomic data were then used to calculate expansion opacities for typical kilonova conditions expected one day after the NS merger; these are a density of ρ = 10⁻¹³ g cm⁻³ and temperatures ranging from T = 5000 K to T = 15 000 K. Some agreements and differences were found when comparing our results with available data, highlighting the importance of using sufficiently complete atomic models for the determination of opacities.