Planet formation theory suggests that planet bulk compositions are likely to reflect the chemical abundance ratios of their host star's photosphere. Variations in the abundance of particular chemical species in stellar photospheres between different galactic stellar populations demonstrate that there are differences among the expected solid planet bulk compositions. We aim to present planetary mass-radius relations of solid planets for kinematically differentiated stellar populations, namely, the thin disc, thick disc, and halo. Using two separate internal structure models, we generated synthetic planets using bulk composition inputs derived from stellar abundances. We explored two scenarios, specifically iron-silicate planets at 0.1 AU and silicate-iron-water planets at 4 AU. We show that there is a persistent statistical difference in the expected mass-radius relations of solid planets among the different galactic stellar populations. At 0.1 AU for silicate-iron planets, there is a 1.51 to 2.04\% mean planetary radius difference between the thick and thin disc stellar populations, whilst for silicate-iron-water planets past the ice line at 4 AU, we calculate a 2.93 to 3.26\% difference depending on the models. Between the halo and thick disc, we retrieve at 0.1 AU a 0.53 to 0.69\% mean planetary radius difference, and at 4 AU we find a 1.24 to 1.49\% difference depending on the model. Future telescopes (such as PLATO) will be able to precisely characterize solid exoplanets and demonstrate the possible existence of planetary mass-radius relationship variability between galactic stellar populations.

We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.