Exoplanet characterization is one of the main foci of current exoplanetary science. For super-Earths and sub-Neptunes, we mostly rely on mass and radius measurements, which allow to derive the body's mean density and give a rough estimate of the planet's bulk composition. However, the determination of planetary interiors is a very challenging task. In addition to the uncertainty in the observed fundamental parameters, theoretical models are limited due to the degeneracy in determining the planetary composition. We aim to study several aspects that affect internal characterization of super-Earths and sub-Neptunes: observational uncertainties, location on the M-R diagram, impact of additional constraints as bulk abundances or irradiation, and model assumptions. We use a full probabilistic Bayesian inference analysis that accounts for observational and model uncertainties. We employ a Nested Sampling scheme to efficiently produce the posterior probability distributions for all the planetary structural parameter of interest. We include a structural model based on self-consistent thermodynamics of core, mantle, high-pressure ice, liquid water, and H-He envelope. Regarding the effect of mass and radius uncertainties on the determination of the internal structure, we find three different regimes: below the Earth-like composition line and above the pure-water composition line smaller observational uncertainties lead to better determination of the core and atmosphere mass respectively, and between them structure characterization only weakly depends on the observational uncertainties. We show that small variations in the temperature or entropy profiles lead to radius variations that are comparable to the observational uncertainty, suggesting that uncertainties linked to model assumptions can become more relevant to determine the internal structure than observational uncertainties.

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.