Coalescence of neutron stars 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 neutron star 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 kappa ~ 20-30 cm^2 g^{-1} for the electron fraction of Ye < 0.20, kappa ~ 3-5 cm^2 g^{-1} for Ye = 0.25-0.35, and kappa ~ 1 cm^2 g^{-1} for Ye = 0.40 at T = 5,000-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 days 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 x 10^{-3}) reproduce the long-lasting near-infrared emission.

Using the multiconfiguration Dirac-Hartree-Fock and the relativistic configuration interaction methods, a consistent set of transition energies and radiative transition data for the lowest 546 (623, 701, 745) states of the $3p^4 3d$, $3s 3p^2 3d^2$, $3s 3p^3 4p$, $3s 3p^4$, $3s^2 3d^3$, $3s^2 3p^2 3d$, $3s^2 3p^2 4d$, $3s^2 3p^2 4s$, $3p^3 3d^2$, $3p^5$, $3s 3p 3d^3$, $3s 3p^3 3d$, $3s 3p^3 4s$, $3s^2 3p 3d^2$, %$3s^2 3p^2 4f$, $3s^2 3p^2 4p$, $3s^2 3p^3$ configurations in Mn~XI (Fe~XII, Co~XIII, Ni~XIV) is provided. The comparison between calculated excitation energies for the $n=4$ states and available experimental values for Fe XII indicate that the calculations are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Lines from these states are prominent in the soft X-rays. With the present calculations, several recent new identifications are confirmed. Other identifications involving $3p^2 4d$ levels in Fe~XII that were found questionable are discussed and a few new assignments are recommended. As some $n=4$ states of the other ions also show large discrepancies between experimental and calculated energies, we reassess their identification. The present study provides highly accurate atomic data for the $n=4$ states of P-like ions of astrophysical interest, for which experimental data are scarce.

Lanthanide elements play important roles as an opacity source in the ejected material from neutron star mergers. Accurate and complete atomic data are necessary to evaluate the opacities and to analyze the observed data. In this paper, we perform extended, ab-initio atomic calculations from Pr II (Z=59) to Gd II (Z=64). By using multiconfiguration Dirac-Hartree-Fock and relativistic configuration-interaction methods, implemented in the general-purpose relativistic atomic structure package GRASP2K, we calculate the energy levels and transition data of electric dipole transitions. These computations are based on strategies (with small variations) of Nd II published by Gaigalas et al. (2019). Accuracy of data is evaluated by comparing computed energy levels with the NIST database or other works. For the energy levels, we obtain the average relative accuracy of 8%, 12%, 6%, 8%, and 7% for Pr II, Pm II, Sm II, Eu II, and Gd II ions, respectively as compared with the NIST data. Accuracy of energy transfer to the wavelength as 3%, 14% and 11% for Pr II, Eu II and Gd II. Our computed E1 type transition probabilities are in good agreement with experimental values presented by other authors especially for strong transitions.

The main aim of this paper is to present accurate energy levels of the ground [Xe]$4f^{12}$ and first excited [Xe]$4f^{11}5d$ configurations of Er$^{2+}$. The energy level structure of the Er$^{2+}$ ion was computed using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction (RCI) methods, as implemented in the GRASP2018 program package. The Breit interaction, self-energy and vacuum polarization corrections were included in the RCI computations. The zero-first-order approach was used in the computations. Energy levels with the identification in $LS$ coupling for all (399) states belonging to the [Xe]$4f^{12}$ and [Xe]$4f^{11}5d$ configurations are presented. Electric dipole (E1) transition data between the levels of these two configurations are computed. The accuracy of the these data are evaluated by studying the behaviour of the transition rates as functions of the gauge parameter as well as by evaluating the cancellation factors. The core electron correlations were studied using different strategies. Root-mean-square deviations obtained in this study for states of the ground and excited configurations from the available experimental or semi-empirical data are 649 cm$^{-1}$, and 747 cm$^{-1}$, respectively.

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.

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 Msun, 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.