Incomplete understanding of key kilonova microphysics (opacity, heating, thermalization)
Characterize the wavelength-dependent ejecta opacities, the radioactive nuclear heating rates, and the thermalization efficiencies governing kilonova emission in order to reduce modeling uncertainties and enable reliable inference of ejecta masses and velocities from observed light curves.
References
In addition, several pieces of kilonova physics are still not well understood even in state-of-the-art simulations, such as the uncertainty in wavelength-dependent opacities, nuclear heating rate, and thermalization efficiencies (Barnes et al. 2021; Bulla 2023; Brethauer et al. 2024; Sarin & Rosswog 2024).
                — Uniform Modeling of Observed Kilonovae: Implications for Diversity and the Progenitors of Merger-Driven Long Gamma-Ray Bursts
                
                (2409.02158 - Rastinejad et al., 3 Sep 2024) in Section 4.5 (Caveats to Kilonova Model)