The radioactively powered transient following a binary neutron star (BNS) merger, known as a kilonova (KN), is expected to enter the steady-state nebular phase a few days after merger. Steady-state holds until thermal reprocessing time-scales become long, at which point the temperature and ionisation states need to be evolved time-dependently. We study the onset and significance of time-dependent effects using the non-local thermodynamic equilibrium (NLTE) spectral synthesis code SUMO. We employ a simple single-zone model with an elemental composition of Te, Ce, Pt and Th, scaled to their respective solar abundances. The atomic data are generated using the Flexible Atomic Code (FAC), and consist of energy levels and radiative transitions, including highly forbidden lines. We explore the KN evolution from 5 to 100 days after merger, varying ejecta mass and velocity. We also consider variations in the degree of electron magnetic field trapping, as well as radioactive power generation for alpha and beta decay (but omitting fission products). We find that the transition time, and magnitude of steady-state deviations are highly sensitive to these parameters. For typical KN ejecta, the deviations are minor within the time-frame studied. However, low density ejecta with low energy deposition show significant differences from ∼10 days. Important deviation of the ionisation structure solution impacts the temperature by altering the overall line cooling. Adiabatic cooling becomes important at t ≥ 60 days which, in addition to the temperature and ionisation effects, lead to the bolometric light curve deviating from the instantaneous radioactive power deposited.

We use magnetograms of 8 solar analogues of ages 30 Myr to 3.6 Gyr obtained from Zeeman Doppler Imaging (ZDI) and taken from the literature, together with two solar magnetograms, to drive magnetohydrodynamical (MHD) wind simulations and construct an evolutionary scenario of the solar wind environment and its angular momentum loss rate. With observed magnetograms of the radial field strength as the only variant in the wind model, we find that power law model fitted to the derived angular momentum loss rate against time, t, results in a spin down relation Ω ∝ t −0.51 , for angular speed Ω, which is remarkably consistent with the well-established Skumanich law Ω ∝ t −0.5 . We use the model wind conditions to estimate the magnetospheric standoff distances for an Earth-like test planet situated at 1 AU for each of the stellar cases, and to obtain trends of minimum and maximum wind ram pressure and average ram pressure in the solar system through time. The wind ram pressure declines with time as P ram ∝ t 2/3 , amounting to a factor of 50 or so over the present lifetime of the solar system.

A binary neutron star merger produces a rapidly evolving transient known as a kilonova (KN), which peaks a few days after merger. Modelling of KNe has often been approached assuming local thermodynamic equilibrium (LTE) conditions in the ejecta. We present the first analysis of non-local thermodynamic equilibrium (NLTE) level populations, using the spectral synthesis code sumo, and compare these to LTE values. We investigate the importance of the radiation field by conducting NLTE excitation calculations with and without radiative transfer. Level populations, in particular higher lying ones, start to show deviations from LTE several days after merger. Excitation is lower in NLTE for the majority of ions and states, and this tends to give lower expansion opacities. While the difference is small for the first few days, it grows to factors 2-10 after this. Our results are important both for demonstrating validity of LTE expansion opacities for an initial phase (less than a week), while highlighting the need for NLTE modelling during later phases. Considering also NLTE ionisation, our results indicate that NLTE can give both higher or lower opacities, depending on composition and wavelength, sometimes by orders of magnitudes.