Open Access
Issue
A&A
Volume 692, December 2024
Article Number A67
Number of page(s) 12
Section Planets, planetary systems, and small bodies
DOI https://doi.org/10.1051/0004-6361/202452193
Published online 03 December 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

With the rapid advance in exoplanet science, it is now possible to characterize not only the fundamental parameters (mass and radius) of planets, but also their interior structure and bulk composition (e.g. Dorn et al. 2018; Nettelmann & Valencia 2021; Adibekyan et al. 2021a,b; Armstrong et al. 2020; Helled et al. 2022; Boley et al. 2023). These parameters are important for determining the habitability conditions of terrestrial planets (e.g. Shahar et al. 2019; Unterborn et al. 2022), and understanding the bulk composition is essential for unraveling the origins of planetary diversity (e.g. Owen & Wu 2017; Jin & Mordasini 2018; Bitsch et al. 2019; Venturini et al. 2020).

Information about the mass and radius of exoplanets alone is insufficient to accurately characterize the interior composition of rocky exoplanets (e.g. Valencia et al. 2007; Dorn et al. 2015; Unterborn et al. 2023). To break this degeneracy, the abundances of the main rocky-forming elements (Mg, Si, and Fe) in the host star are often used. Stellar atmospheres offer the only observational window into the remnants of the planet-forming disk’s composition, as these disks, with lifetimes typically lasting only a few million years, have long dissipated around most discovered planet-hosting stars.

When modeling the interiors of rocky planets, a one-to-one correlation between the composition of the planet and that of its host star is typically assumed (e.g. Thiabaud et al. 2015; Dorn et al. 2015; Unterborn et al. 2016). This assumption, to a first-order, appears to be valid for Earth (McDonough 2003), Venus (Aitta 2012), Mars (Khan et al. 2018), and the Sun (see Dorn et al. 2015; Unterborn et al. 2016). It may also be valid for some exoplanets (Santos et al. 2015; Schulze et al. 2021), but not necessarily for others (Bonomo et al. 2019; Plotnykov & Valencia 2020; Brinkman et al. 2023; Mah & Bitsch 2023).

Assuming a one-to-one compositional correlation is a reasonable starting point, but it breaks down even within our own Solar System, in the case of Mercury (Morgan & Anders 1980). Adibekyan et al. (2021a)1 conducted a statistical comparison of the iron-to-silicate mass fraction in low-mass exoplanets and their host stars, concluding that while there is a correlation between the composition of rocky planets and their host stars, the relationship is not one-to-one. The analysis focused on 22 low-mass (M < 10 M)) and small-sized (R < 2 R)) planets with mass and radius measurements precise to within 30%, all orbiting F, G, and K-type stars. However, the study by Adibekyan et al. (2021a) has potential limitations: (i) The planets’ interiors were modeled under the assumption that they consist only of a core and mantle, with no volatiles included, and (ii) the present-day stellar abundances were used as a proxy for the primordial protoplanetary disk composition. Some astrophysical processes, such as atomic diffusion, can alter a star’s surface composition as it evolves (Deal et al. 2018), meaning that the present-day stellar abundances could differ significantly from their primordial values.

In this work, we study the star-planet compositional link by addressing the aforementioned limitations. Additionally, our sample is approximately 50% larger than that of Adibekyan et al. (2021a) and includes more precise measurements of planetary mass and radius. The paper is structured as follows: in Sect. 2, we introduce the sample. Section 3 details the methods used to determine both the primordial composition of the planet-forming disk and the composition of the exoplanets. The star-planet compositional link is discussed in Sect. 4, and we provide a summary of the work in Sect. 5.

2 The sample

We started our sample selection using the PlanetS Catalog (Otegi et al. 2020)2, which provides a list of transiting exoplanets with relative uncertainties in mass and radius better than 25% and 8%, respectively. Given our focus on potentially rocky planets, we selected only those with radii below 2 R to exclude likely gas-rich sub-Neptunes (Fulton et al. 2017; Owen & Wu 2017; Venturini et al. 2020), orbiting F, G, and K-type stars with effective temperatures (Teff) ranging between 4600 and 6500 K. This temperature range allows for the precise determination of stellar atmospheric abundances of relevant elements. We excluded Kepler-36 b and Kepler-36 c which have masses determinations via transit timing variation (TTV) method as there might be a discrepancy in the density distributions of planets with mass determinations via TTV and Radial Velocity (RV) method (e.g. Mills & Mazeh 2017; Leleu et al. 2023; Adibekyan et al. 2024). Additionally, we excluded planets orbiting Kepler-65, TOI-1807, TOI-1416, and TOI-1444 due to the unavailability of high-resolution spectra for their host stars. This resulted in a sample of 32 exoplanets orbiting 30 stars, for which we could obtain stellar abundances determined from high-resolution spectra. The distribution of these planets on the Mass-Radius diagram is shown in Fig. 1, color-coded by their equilibrium temperatures (Teq), as reported in the PlanetS Catalog. The properties of the planets are presented in Table A.1.

thumbnail Fig. 1

Distribution of planets in the mass–radius diagram color-coded by their Teq. The curves display the expected mass–radius relationship for different compositions taken from Luo et al. (2024).

3 Present-day and primordial compositions

When conducting stellar spectroscopy, the present-day properties of stars are determined. However, certain properties, including the atmospheric chemical abundances, can change as stars age, even while they remain in the main sequence. In this section, we describe how we determined the primordial compositions of stars hosting planets, and how we used this information to compute the composition of the building blocks of those planets.

3.1 Present-day abundances of the host stars

Given the importance of a consistent characterization of the host stars, we compiled the stellar parameters and chemical abundances of the sample stars from five published articles where these properties were derived using the same methodology and specific tools, making them directly comparable and usable as a single set. This methodology was also applied to K2-131, EPIC 249893012, and HD 20329 in this study. We re-derived the parameters and abundances of EPIC 249893012 and HD 20329 using higher quality spectra obtained with the ESPRESSO spectrograph (Pepe et al. 2021). These parameters and their corresponding references are listed in Table A.2.

In all cases, the stellar atmospheric parameters (Teff, surface gravity log ɡ, microturbulent velocity Vtur, and metallicity [Fe/H]) and the abundances were determined following the methods outlined in our previous works (e.g. Santos et al. 2013; Sousa 2014; Adibekyan et al. 2012, 2015). In our local thermody-namic equilibrium (LTE) analysis, we measured the equivalent widths (EW) of iron lines using the ARES v2 code3 (Sousa et al. 2015), ensuring excitation and ionization equilibrium. We employed a grid of Kurucz model atmospheres (Kurucz 1993) and the 2014 version of the radiative transfer code MOOG (Sneden 1973). The abundances of Mg and Si were derived using the classical curve-of-growth analysis method under the LTE assumption. Although the EWs of the spectral lines were automatically measured with ARES, we visually inspected the EW measurements to ensure accuracy.

Since determining the abundances of C and O is particularly challenging for stars cooler than about 5200 K, we decided to maintain homogeneity by predicting the abundances of these elements using a machine-learning model developed in Adibekyan et al. (2024). This model predicts the abundances of C and O based on the input abundances of Fe, Si, and Mg wth a precision of about 0.07 dex.

3.2 Primordial abundances of the host stars

We used a grid-based modeling approach (similar to Deal et al. 2021) to determine the ages and primordial abundances of the planet hosting stars. The grid of stellar models is computed with the public Code d’Evolution Stellaire Adaptatif et Modulaire4 (Cesam2k20, Morel & Lebreton 2008; Marques et al. 2013; Deal et al. 2018). The parameter space of the grid varies in mass, primordial iron abundance [Fe/H]prim and helium mass fraction Yprim, and the core overshoot parameter αov. In terms of chemical transport, the grid includes atomic diffusion with radiative accelerations (Deal et al. 2018), the transport induced by rotation (Marques et al. 2013), and we added an extra turbulent mixing calibrated on the Sun to reproduce its lithium surface abundance. The transport of chemicals is followed for H, He, Li, Be, B, C, N, O, Ne, Na, Mg, Si, Al, S, Ca, Fe, and some of their isotopes. The input physics is the same as the models presented in Deal et al. (2020) for atomic diffusion, rotation and the microphysics, except that the primordial metal mixture follows Asplund et al. (2021), the T(τ) relation in the atmosphere follows Vernazza et al. (1981), and the convection formalism follows Böhm-Vitense (1958) with a solar-calibrated mixing-length parameter αMLT = 1.998.

The stellar property inferences are performed with the Bayesian Asteroseismic Inference on a Massive Scale public code5 (AIMS, Rendle et al. 2019). The stellar properties are inferred using log g, Teff and [Fe/H] as constraints, in a configuration of the code where no seismic constraints are required. Note that the observed [Fe/H] is directly compared to the one predicted by the models and not to the total metallicity [M/H] for a better accuracy (Deal et al. 2018; Moedas et al. 2022). From the posterior distributions obtained for each star, we derive the fundamental stellar properties (mass, radius, and age) and the primordial chemical composition from each elements followed by the code. The provided error bars are the 16th and 84th per-centiles of the distributions. The ages and primordial abundances of the sample stars are presented in Table A.3.

In Fig. 2, we display the amplitude of the iron abundance change as a function of the age and temperature of the stars. The figure illustrates that the oldest stars experience the most significant reduction in iron abundance and that, for a given age, atomic diffusion is more pronounced in hotter stars due to their narrower convective envelopes. Furthermore, the figure highlights that the temporal change in metallicity can be several times larger than the uncertainties in the present-day [Fe/H] measurements.

The overall trend depicted in Fig. 2 is consistent when considering other elements used in this work. This consistency is important as it suggests that the changes in abundance ratios remain negligible, as illustrated in Fig. 3.

It is important to note that the impact of atomic diffusion generally depends on the competing transport processes included in the models (Moedas et al. 2022). A different treatment of rotation or additional mixing could also influence the results (see e.g. Deal et al. 2020). Furthermore, for hotter stars above ~6500 K, atomic diffusion may cause more pronounced changes in abundance ratios. However, since we have added extra mixing calibrated to reproduce the Sun’s lithium abundance, any variations in the efficiency of rotation or diffusion would be compensated for by this calibration, resulting in a negligible impact on our results.

thumbnail Fig. 2

Difference between primordial and present-day abundances of iron as a function of age. Crosses show the data points, color-coded by Teff. Red error bars represent the uncertainties of present-day [Fe/H].

thumbnail Fig. 3

Difference between primordial and present-day abundance ratios of [Fe/Si] and [Fe/Mg] as a function of age. Crosses show the data points. Red error bars represent the uncertainties of present-day [Fe/Si] and [Fe/Mg].

thumbnail Fig. 4

Difference in iron-to-silicate mass (top) and water mass (bottom) fractions of planet building blocks determined from the present-day and primordial stellar abundances. The differences are plotted as a function of firon,prim star $f_{{\rm{iron,prim }}}^{{\rm{star }}}$ and fwater,prim star $f_{{\rm{water,prim }}}^{{\rm{star }}}$, respectively.

3.3 Composition of planet building blocks

We converted the derived primordial atmospheric abundances of C, O, Mg, Si, and Fe into absolute abundances (Adibekyan 2019) using solar reference values from Asplund et al. (2021) and assuming a hydrogen atom number ϵH of 1012. These absolute stellar abundances were then used to estimate the iron-to-silicate mass fraction (firon star )$\left( {f_{{\rm{iron }}}^{{\rm{star }}}} \right)$ and the water mass fraction (fwater star )$\left( {f_{{\rm{water }}}^{{\rm{star }}}} \right)$ of the protoplanetary disks using the stoichiometric model described in Santos et al. (2015, 2017). Note that firon star $f_{{\rm{iron }}}^{{\rm{star }}}$ does not depend on water availability. For the protosolar disk, the model predicts an firon $f_{{\rm{iron }}}^ \oplus $ of 32.3±2.1% and an fwater$f_{{\rm{water}}}^ \oplus $ of 58.2±3.2% using solar abundances from Asplund et al. (2021). The present-day and primordial compositions of planet building blocks are displayed in Table A.4.

In Fig. 4, we illustrate the impact of using primordial stellar abundances instead of present-day values on firon star $f_{{\rm{iron }}}^{{\rm{star }}}$ and fwaterstar $f_{{\rm{water}}}^{{\rm{star }}}$. The largest difference in the iron-to-silicate mass fraction is approximately 0.4%, with the lowest uncertainty for this parameter around 1% and a mean uncertainty of about 2.2%. A similar trend is observed for the water mass fraction, with the largest difference of about 0.8%, which is approximately 6 times and 8 times smaller than the minimum and mean uncertainty of this parameter, respectively. This behavior is consistent with the results from Sect. 3.2, where we demonstrated that while individual element abundances experience significant temporal changes, their ratios remain relatively stable over time. Consequently, when using Fe/Mg and Fe/Si mineralogical ratios to determine the composition of planet-building blocks, it is generally acceptable to use present-day stellar atmospheric abundances for these low temperature stars.

In Fig. 5, we show how the scaled densities of planets (ρ/ρEarth−like) depend on the primordial iron-to-silicate mass fraction (firon,primstar )$\left( {f_{{\rm{iron,prim}}}^{{\rm{star }}}} \right)$ and the primordial water mass fraction (firon,primstar )$\left( {f_{{\rm{iron,prim}}}^{{\rm{star }}}} \right)$. The scaled density is calculated by dividing the bulk density (ρ) of the planet by the density of a planet with Earth-like composition (ρEarth−like, Luo et al. 2024) for the same planetary mass. Throughout this paper, unless stated otherwise, we performed Ordinary Least Squares (OLS) regression using the statsmodels package to quantify relationships. The slopes and associated uncertainties were then used in t-statistics (via SCIPY) to test the null hypothesis that the data can be modeled with a zero slope. The null hypothesis is rejected if the resulting two-tailed p-values (P(t-stat)) are below the typical threshold of 0.05. We opted for OLS over Orthogonal Distance Regression (ODR), which was used in Adibekyan et al. (2021a), based on the discussion in Brinkman et al. (2024), which concluded that there is a non-negligible probability of obtaining a statistically significant relation when using ODR on a sample similar to ours.

Compared to Adibekyan et al. (2021a), the top panel of Fig. 5 shows two low-density planets (TOI-733b and HD 23472 c) orbiting stars with solar firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$. In Fig. 1, these two planets are the farthest from the Earth-like composition curve, and our interior analysis suggests that they contain a significant amount of water (above 10%, see Table A.5), potentially in the form of steam. We excluded these two planets from the regression analysis shown in Fig. 5. It is important to note that firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$ refers to the iron-to-silicate mass ratio and does not account for water in the calculation.

The top panel of Fig. 5 depicts that low-density planets are likely formed in protoplanetary disks with a low iron-to-silicate mass fraction, with an OLS regression p-value of 0.002 suggesting a statistically significant correlation between scaled density and firon,prim star $f_{{\rm{iron,prim }}}^{{\rm{star }}}$. This correlation was previously noted in Adibekyan et al. (2021a). It is important to distinguish that firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$ and [Fe/H] of the host stars are different parameters. Furthermore, there is no statistically significant correlation (p-value ~ 0.08) between scaled density and [Fe/H].

The bottom panel of Fig. 5 suggests a possible negative correlation between scaled density and fwater,primstar $f_{{\rm{water,prim}}}^{{\rm{star }}}$. However, this trend is not statistically significant, with a p-value of 0.07. Overall, both low iron content and high water content are likely factors contributing to the low density of rocky planets.

Interestingly, the density gap between super-Mercuries and super-Earths appears less pronounced compared to the findings of Adibekyan et al. (2021a), which is in line with the results of some resent studies (Brinkman et al. 2024; Schulze et al. 2024). This discrepancy can be attributed to updated mass and radius determinations and stronger precision constraints on these parameters for the sample planets. Notably, the four highest-density planets in our sample orbit stars with above-average firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$.

Additionally, HD 20329 b, which orbits a star with a very low firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$, has a scaled density similar to that of Earth. This suggests that HD 20329 b might have a significantly higher fironplanet $f_{{\rm{iron}}}^{{\rm{planet }}}$ (see. Fig. 6) than what would be expected from its host star’s protoplanetary disk composition, potentially classifying it as a super-Mercury (Schulze et al. 2021). Adibekyan et al. (2021a) previously proposed that the formation of high-density super-Mercuries might be linked to the composition of the host star, a hypothesis that has been further explored theoretically in Johansen & Dorn (2022); Mah & Bitsch (2023).

3.4 Interior composition of exoplanets

To compare stellar abundances with planetary abundances, we use inference methods to estimate the confidence regions for the planetary iron-to-silicate mass fraction (firon planet )$\left( {f_{{\rm{iron }}}^{{\rm{planet }}}} \right)$, the water mass fraction (fwater,primstar )$\left( {f_{{\rm{water,prim}}}^{{\rm{star }}}} \right)$, and the iron-to-magnesium ratio ((Fe/)planet).

We use an interior model with 3 components: core, mantle, and water (mainly steam). The model is based on Dorn et al. (2015) with updates from Luo et al. (2024). The water can be present in the core, the mantle, and the surface, depending on the specific thermal state of the planets. In brief, the interior model uses different sets of equation of states for the iron in the core (Hakim et al. 2018; Miozzi et al. 2020; Luo et al. 2024), for silicates in the mantle (Fischer et al. 2011; Faik et al. 2018; Hemley et al. 1992; Musella et al. 2019; Connolly 2009), and surface water (mainly in steam phase for the planets of interest Haldemann et al. 2020). We assume adiabatic temperature profiles in the core, mantle and steam layer, except for pressures below 0.1 bar in the steam layer where we assume an isothermal profile. Additionally, we assumed a pure iron core and fixed the Mg/Si number ratio to 1.04, which is representative of Earth-like composition. The planet radius is set at 1 mbar.

Water can be added to the mantle melts, while the solid mantle is assumed to be dry. The addition of water reduces the density, for which we follow Bajgain et al. (2015). For small water mass fractions, this reduction is nearly independent of pressure and temperature. Further, water can be present in the core for which we use the equation of states of Luo et al. (2024), that accounts for the reduction in density. Further details on the model can be found in Luo et al. (2024). The composition of the planets are presented in Table A.5.

thumbnail Fig. 5

Scaled density of planets as a function of firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$ (top) and fwater,primstar $f_{{\rm{water,prim}}}^{{\rm{star }}}$ (bottom) color-coded by the equilibrium temperature of the planets. The positions of Earth is indicated with its symbols in black. The black solid lines represent the results of OLS regression without considering TOI-733 b and HD 23472 c. The red dashed-dotted line represent the visual separation between super-Earths and super-Mercuries. The black dashed line represent the result of the linear regression without considering TOI-733 b, HD 23472 c and the super-Mercuries located above the red line (p-value of about 10−5).

thumbnail Fig. 6

Iron-to-silicate and water mass fractions of planets and their host stars color-coded by the equilibrium temperature of the planets. The black solid line represents the results of OLS regression.

4 Star-planet compositional link

Having independently determined the iron-to-silicate mass fractions from both the stellar composition, and the planetary mass and radius, we can now examine the compositional link between solar-type stars and their low-mass planets. The top panel of Fig. 6 shows the relationship between planetary and stellar firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$. The OLS regression analysis, with a slope of 2.6±0.9 and p-value of 0.009, indicates a significant correlation between these two parameters. Given the ongoing discussion in the literature about this relationship, and to ensure the robustness of this result, we applied different regression methods (Appendix B), accounting for the uncertainties in the parameters. While the exact form of the relation varies with the methodology used, the presence of a steeper-than-one-to-one correlation remains evident. The figure also reveals that the range of planetary iron-to-silicate mass fractions is significantly larger than that of the protoplanetary disk, confirming the previous findings (Plotnykov & Valencia 2020; Adibekyan et al. 2021a).

In the bottom panel of Fig. 6, we display the relationship between planetary and protoplanetary water mass fractions. The absence of a clear correlation is somewhat expected, given that fwater,primstar $f_{{\rm{water,prim}}}^{{\rm{star }}}$ is calculated with the assumption that the planet-building blocks are located outside the water snowline. The low water content observed in most of the studied planets may actually suggest that they formed inside the water snowline. Several recent studies have also suggested a low water content in close-in super-Earths (Müller et al. 2024; Rogers et al. 2024).

We further explore the planet-star compositional link by examining the ratio of fironplanet $f_{{\rm{iron}}}^{{\rm{planet }}}$ to firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$ as a function of planetary radius (top panel of Fig. 7) and planetary mass (bottom panel of Fig. 7). There is a clear correlation between this ratio and planetary radius (p-value < 0.001), and no correlation with planetary mass (p-value of about 0.9). This is most likely due to the fact that planetary density depends on radius with a cubic relation, while the dependence on mass is only linear. However, the top panel also shows that for a given radius, denser planets have a higher ratio of fironplanet $f_{{\rm{iron}}}^{{\rm{planet }}}$ to firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$.

As planetary interiors are often modeled under the assumption that the Fe/Mg ratio in planets is similar to that of their host stars, Fig. 8 illustrates how the Fe/Mg planet-to-host ratio varies with planetary radius. The figure reveals that the Fe/Mg ratio in planets can be up to three times lower or up to ten times higher than in their host stars, indicating a large variation between the formed planets and the protoplanetary disk. This substantial variation suggests that processes occurring within the protoplanetary disk during planet formation play a significant role in determining the final composition of rocky planets. Additionally, Fig. 8 shows a clear dependence of the (Fe/)planet/(Fe/)star ratio on planetary radius.

thumbnail Fig. 7

The ratio of iron-to-silicate mass fraction of planets and the planet building blocks as a function of planetary radius (top panel) and mass (bottom panel). The black solid lines represent the results of the linear regression. The symbols in the top panel are color-coded by the scaled density of the planets.

thumbnail Fig. 8

The Fe/Mg abundance ratio of planets and their host stars as a function of planetary radius. The black solid line in the plot represent the results of the OLS regression (p-value ~ 0.006) performed on the linearly scaled data.

5 Summary

In this work, we analyzed a sample of 32 low-mass planets orbiting 30 FGK-type stars to investigate the compositional link between stars and their planets. The selected planets from the PlanetS Catalog have relative uncertainties in mass and radius better than 25% and 8%, respectively. Compared to our previous study (Adibekyan et al. 2021a), the current sample includes about 50% more planets, and our analysis was enhanced by determining the primordial compositions of the hosts and using a planetary interior model that accounts also for presence water (primarily in steam form).

We derived the primordial stellar abundances of C, O, Mg, Si, and Fe from their present-day atmospheric values, demonstrating that while individual elemental abundances can vary significantly over time, the ratios between them remain relatively stable. Note that this may not be true for hotter stars, and that the exact values are subject to the assumptions and limitations of the modeling.

We then used the primordial abundances as input for the stoi-chiometric model presented in Santos et al. (2015) to estimate the iron-mass fraction (firon,primstar )$\left( {f_{{\rm{iron,prim}}}^{{\rm{star }}}} \right)$ and water-mass fraction (fwater,primstar )$\left( {f_{{\rm{water,prim}}}^{{\rm{star }}}} \right)$ of the protoplanetary disks. When comparing the present-day and primordial compositions of the planet-building disks, we find negligible differences, as the abundance ratios of Fe/Mg and Fe/Si remain relatively constant over the evolution of stars in the main sequence. This consistency suggests that using present-day stellar abundances to infer the composition of protoplanetary disks is a valid approach, as the key elemental ratios critical for planet formation do not change significantly over time.

Our analysis reveal a strong correlation between the scaled density of rocky planets and the primordial iron-to-silicate mass fraction of their protoplanetary disks (firon,primstar )$\left( {f_{{\rm{iron,prim}}}^{{\rm{star }}}} \right)$. Although the density gap between super-Mercuries and super-Earths is less pronounced than previously observed in Adibekyan et al. (2021a), it remains visible.

We employed a three-component planetary interior model (Dorn et al. 2015; Luo et al. 2024) to derive the compositions of the planets based on their mass and radius. While we find a statistically significant correlation between the iron-to-silicate mass fraction of planets and their host stars, no such correlation is observed for the water mass fraction. This lack of correlation may be due to the assumption in the stoichiometric model that the building blocks of planets formed outside the water snowline, while the observed low water content suggests formation inside the snowline.

As already found in Adibekyan et al. (2021a), the range of fironplanet $f_{{\rm{iron}}}^{{\rm{planet }}}$ is significantly broader than that of their building blocks. This indicates that while the primordial composition of the protoplanetary disk sets a baseline for planetary composition, various processes during planet formation and evolution can lead to substantial deviations in the iron-to-silicate mass fraction of the resulting planets compared to their original building materials.

Additionally, we compared the Fe/Mg ratios in planets and their host stars, finding a variation from 0.3 to 10 times between the two. This ratio was also found to correlate with planetary radius, further indicating that processes within the protoplanetary disk significantly influence the final composition of rocky planets.

In conclusion, we confirm the previous results (Adibekyan et al. 2021a) suggesting that while a correlation exists between the composition of rocky planets and their host stars, the assumption of a strict one-to-one elemental abundance relationship is not supported by our observations.

Acknowledgements

This work was supported by Fundação para a Ciência e Tecnologia through national funds and by FEDER through COMPETE2020 -Programa Operacional Competitividade e Internacionalização by these grants: UIDB/04434/2020; UIDP/04434/2020; 2022.06962.PTDC. Funded/Co-funded by the European Union (ERC, FIERCE, 101052347). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. TLC is supported by Fundação para a Ciência e a Tecnologia (FCT) in the form of a work contract (CEECIND/00476/2018). We would like to thank the referee, Joe Schulze, for the constructive comments and suggestions, which have significantly improved both the quality and clarity of our work. This research has made use of the NASA Exoplanet Archive operated by the California Institute of Technology under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. In this work, we used the Python language and several scientific packages: Numpy (van der Walt et al. 2011), Scipy (Virtanen et al. 2020), Pandas (Wes McKinney 2010), statsmodels (Skipper Seabold & Josef Perktold 2010), and Matplotlib (Hunter 2007).

Appendix A Tables with compositions of stars and their planets

Table A.1

The main properties of the planets.

Table A.2

The stellar atmospheric parameters and present-day abundances of the planet host stars.

Table A.3

The ages and primordial abundances of the planet host stars.

Table A.4

The present-day and primordial compositions of planet building blocks.

Table A.5

The compositions of planets.

Appendix B The star-planet compositional relationship

Brinkman et al. (2024) recently suggested that both the slope and the significance of the relationship between the iron mass fraction of stars and their rocky planets depend on the fitting algorithm used. In this section, we apply several algorithms to test the robustness of this relationship.

The OLS regression used in the manuscript does not account for the uncertainties in and fironplanet $f_{{\rm{iron}}}^{{\rm{planet }}}$ and firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$. To address this, we performed a Monte Carlo (MC) test. Assuming Gaussian distributions for the uncertainties in these parameters, centered on their respective values, we randomly sampled values for fironplanet $f_{{\rm{iron}}}^{{\rm{planet }}}$ and firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$ neglecting the cases when the values are negative. After 10000 realizations, we obtained a median slope of 1.65±0.78, with a corresponding p-value (P(t-stat)) of 0.043.

As a second method, we applied the Weighted Least Squares (WLS) regression as implemented in statsmodels package, using the inverse square of the uncertainties in fironplanet $f_{{\rm{iron}}}^{{\rm{planet }}}$ as weights. This yielded a slope of 3.85±1.10 with a p-value of 0.002. We then repeated the previously described MC test for the WLS case. After 10 000 realizations, we obtained a median slope of 2.43±0.89, with a p-value of 0.01.

For the third method, we used Orthogonal Distance Regression (ODR), implemented in SCIPY, which accounts for uncertainties in both parameters during fitting. ODR resulted in a slope of 5.98±1.08, with a p-value of about 10−5. Repeating the MC test for ODR over 10000 realizations, we obtained a median slope of 5.85±1.07, leading to a p-value of ~ 10−5.

Our fourth method was Total Least Squares (TLS), a generalization of ODR that also considers observational errors in both variables. For the log-likelihood calculation, we used astroML package which follows the prescription of Hogg et al. (2010), and then used the scipy.optimize package to maximize it and find the best-fit parameters. After performing the MC test for the TLS method with 10 000 realizations, we obtained a median slope of 24.9±0.9, with a p-value below 10−22.

Finally, we applied PyStan, a Python interface to the Bayesian inference package Stan (Carpenter et al. 2017) which allowed us to account for uncertainties in both parameters within a Bayesian framework. We used wide priors for the slope, intercept, and sigma (the scale of the noise term), modeled as normal distributions centered at zero with a standard deviation of 100. From the posterior distributions of the slopes, we obtained a median value of 3.06, with a 95% credible interval of [0.8, 5.0]. Notably, less than 4% of the posterior slopes were smaller than 1.

Overall, these tests demonstrate a statistically significant relationship between fironplanet $f_{{\rm{iron}}}^{{\rm{planet }}}$ and firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$, with a slope greater than one. In Fig. B.1 we show the results of the aforementioned linear fits.

thumbnail Fig. B.1

Iron-to-silicate mass fractions of planets and their host stars, with different lines representing the results of the linear fits discussed in the manuscript.

References

  1. Adibekyan, V. 2019, Geosciences, 9, 105 [Google Scholar]
  2. Adibekyan, V. Z., Sousa, S. G., Santos, N. C., et al. 2012, A&A, 545, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Adibekyan, V., Figueira, P., Santos, N. C., et al. 2015, A&A, 583, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Adibekyan, V., Dorn, C., Sousa, S. G., et al. 2021a, Science, 374, 330 [NASA ADS] [CrossRef] [Google Scholar]
  5. Adibekyan, V., Santos, N. C., Dorn, C., et al. 2021b, Commun. Byurakan Astrophys. Observ., 68, 447 [NASA ADS] [CrossRef] [Google Scholar]
  6. Adibekyan, V., Sousa, S. G., Delgado Mena, E., et al. 2024, A&A, 683, A159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Aitta, A. 2012, Icarus, 218, 967 [NASA ADS] [CrossRef] [Google Scholar]
  8. Armstrong, D. J., Lopez, T. A., Adibekyan, V., et al. 2020, Nature, 583, 39 [Google Scholar]
  9. Asplund, M., Amarsi, A. M., & Grevesse, N. 2021, A&A, 653, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bajgain, S., Ghosh, D. B., & Karki, B. B. 2015, Nat. Commun., 6, 8578 [NASA ADS] [CrossRef] [Google Scholar]
  11. Barros, S. C. C., Demangeon, O. D. S., Alibert, Y., et al. 2022, A&A, 665, A154 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bitsch, B., Raymond, S. N., & Izidoro, A. 2019, A&A, 624, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Böhm-Vitense, E. 1958, ZAp, 46, 108 [NASA ADS] [Google Scholar]
  14. Boley, K. M., Panero, W. R., Unterborn, C. T., et al. 2023, ApJ, 954, 202 [NASA ADS] [CrossRef] [Google Scholar]
  15. Bonomo, A. S., Zeng, L., Damasso, M., et al. 2019, Nat. Astron., 3, 416 [Google Scholar]
  16. Brinkman, C. L., Weiss, L. M., Dai, F., et al. 2023, AJ, 165, 88 [NASA ADS] [CrossRef] [Google Scholar]
  17. Brinkman, C. L., Polanski, A. S., Huber, D., et al. 2024, AJ, submitted, [arXiv:2409.08361] [Google Scholar]
  18. Carpenter, B., Gelman, A., Hoffman, M. D., et al. 2017, J. Statist. Softw., 76, 1 [NASA ADS] [CrossRef] [Google Scholar]
  19. Connolly, J. A. D. 2009, Geochem. Geophys. Geosyst., 10 [CrossRef] [Google Scholar]
  20. Damasso, M., Rodrigues, J., Castro-González, A., et al. 2023, A&A, 679, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Deal, M., Alecian, G., Lebreton, Y., et al. 2018, A&A, 618, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Deal, M., Goupil, M. J., Marques, J. P., Reese, D. R., & Lebreton, Y. 2020, A&A, 633, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Deal, M., Cunha, M. S., Keszthelyi, Z., Perraut, K., & Holdsworth, D. L. 2021, A&A, 650, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Dorn, C., Khan, A., Heng, K., et al. 2015, A&A, 577, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Dorn, C., Mosegaard, K., Grimm, S. L., & Alibert, Y. 2018, ApJ, 865, 20 [NASA ADS] [CrossRef] [Google Scholar]
  26. Faik, S., Tauschwitz, A., & Iosilevskiy, I. 2018, Comput. Phys. Commun., 227, 117 [NASA ADS] [CrossRef] [Google Scholar]
  27. Fischer, R. A., Campbell, A. J., Shofner, G. A., et al. 2011, Earth Planet. Sci. Lett., 304, 496 [CrossRef] [Google Scholar]
  28. Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, AJ, 154, 109 [Google Scholar]
  29. Hakim, K., Rivoldini, A., Van Hoolst, T., et al. 2018, Icarus, 313, 61 [Google Scholar]
  30. Haldemann, J., Alibert, Y., Mordasini, C., & Benz, W. 2020, A&A, 643, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Helled, R., Werner, S., Dorn, C., et al. 2022, Exp. Astron., 53, 323 [NASA ADS] [CrossRef] [Google Scholar]
  32. Hemley, R. J., Stixrude, L., Fei, Y., & Mao, H. K. 1992, in High-Pressure Research: Application to Earth and Planetary Sciences (American Geophysical Union (AGU)), 183 [Google Scholar]
  33. Hogg, D. W., Bovy, J., & Lang, D. 2010, arXiv e-prints [arXiv: 1008.4686] [Google Scholar]
  34. Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [NASA ADS] [CrossRef] [Google Scholar]
  35. Jin, S., & Mordasini, C. 2018, ApJ, 853, 163 [Google Scholar]
  36. Johansen, A., & Dorn, C. 2022, A&A, 662, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Khan, A., Liebske, C., Rozel, A., et al. 2018, J. Geophys. Res. (Planets), 123, 575 [NASA ADS] [CrossRef] [Google Scholar]
  38. Kurucz, R. L. 1993, SYNTHE spectrum synthesis programs and line data [Google Scholar]
  39. Leleu, A., Delisle, J. B., Udry, S., et al. 2023, A&A, 669, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Liu, Z., & Ni, D. 2023, A&A, 674, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Luo, H., Dorn, C., & Deng, J. 2024, Nat. Astron. submitted [arXiv:2401.16394] [Google Scholar]
  42. Mah, J., & Bitsch, B. 2023, A&A, 673, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. McDonough, W. F. 2003, Treatise Geochem., 2, 568 [NASA ADS] [Google Scholar]
  45. Mills, S. M., & Mazeh, T. 2017, ApJ, 839, L8 [Google Scholar]
  46. Miozzi, F., Matas, J., Guignot, N., et al. 2020, Minerals, 10, 100 [CrossRef] [Google Scholar]
  47. Moedas, N., Deal, M., Bossini, D., & Campilho, B. 2022, A&A, 666, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Morel, P., & Lebreton, Y. 2008, Ap&SS, 316, 61 [Google Scholar]
  49. Morgan, J. W., & Anders, E. 1980, PNAS, 77, 6973 [CrossRef] [Google Scholar]
  50. Müller, J., Bitsch, B., & Schneider, A. D. 2024, A&A, 688, A139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Musella, R., Mazevet, S., & Guyot, F. 2019, Phys. Rev. B, 99, 064110 [NASA ADS] [CrossRef] [Google Scholar]
  52. Nettelmann, N., & Valencia, D. 2021, in ExoFrontiers; Big Questions in Exoplanetary Science, ed. N. Madhusudhan, 16 [Google Scholar]
  53. Osborn, A., Armstrong, D. J., Cale, B., et al. 2021, MNRAS, 507, 2782 [NASA ADS] [CrossRef] [Google Scholar]
  54. Otegi, J. F., Bouchy, F., & Helled, R. 2020, A&A, 634, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Owen, J. E., & Wu, Y. 2017, ApJ, 847, 29 [Google Scholar]
  56. Pepe, F., Cristiani, S., Rebolo, R., et al. 2021, A&A, 645, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Plotnykov, M., & Valencia, D. 2020, MNRAS, 499, 932 [CrossRef] [Google Scholar]
  58. Rendle, B. M., Buldgen, G., Miglio, A., et al. 2019, MNRAS, 484, 771 [Google Scholar]
  59. Rogers, J. G., Schlichting, H. E., & Young, E. D. 2024, ApJ, 970, 47 [NASA ADS] [CrossRef] [Google Scholar]
  60. Santos, N. C., Sousa, S. G., Mortier, A., et al. 2013, A&A, 556, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Santos, N. C., Adibekyan, V., Mordasini, C., et al. 2015, A&A, 580, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Santos, N. C., Adibekyan, V., Dorn, C., et al. 2017, A&A, 608, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Schulze, J. G., Wang, J., Johnson, J. A., et al. 2021, Planet. Sci. J., 2, 113 [NASA ADS] [CrossRef] [Google Scholar]
  64. Schulze, J. G., Wang, J., Johnson, J. A., et al. 2024, Planet. Sci. J., 5, 71 [NASA ADS] [CrossRef] [Google Scholar]
  65. Shahar, A., Driscoll, P., Weinberger, A., & Cody, G. 2019, Science, 364, 434 [NASA ADS] [CrossRef] [Google Scholar]
  66. Skipper S., & Josef P. 2010, in Proceedings of the 9th Python in Science Conference, eds. S. Walt, & J. Millman, 92 [Google Scholar]
  67. Sneden, C. A. 1973, PhD thesis, The University OF Texas at Austin, USA [Google Scholar]
  68. Sousa, S. G. 2014, in Determination of Atmospheric Parameters of B, 297 [CrossRef] [Google Scholar]
  69. Sousa, S. G., Santos, N. C., Adibekyan, V., Delgado-Mena, E., & Israelian, G. 2015, A&A, 577, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Spaargaren, R. J., Wang, H. S., Mojzsis, S. J., Ballmer, M. D., & Tackley, P. J. 2023, ApJ, 948, 53 [CrossRef] [Google Scholar]
  71. Thiabaud, A., Marboeuf, U., Alibert, Y., Leya, I., & Mezger, K. 2015, A&A, 580, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Unterborn, C. T., Dismukes, E. E., & Panero, W. R. 2016, ApJ, 819, 32 [Google Scholar]
  73. Unterborn, C. T., Foley, B. J., Desch, S. J., et al. 2022, ApJ, 930, L6 [NASA ADS] [CrossRef] [Google Scholar]
  74. Unterborn, C. T., Desch, S. J., Haldemann, J., et al. 2023, ApJ, 944, 42 [Google Scholar]
  75. Valencia, D., Sasselov, D. D., & O’Connell, R. J. 2007, ApJ, 665, 1413 [Google Scholar]
  76. van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22 [Google Scholar]
  77. Venturini, J., Guilera, O. M., Haldemann, J., Ronco, M. P., & Mordasini, C. 2020, A&A, 643, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 [Google Scholar]
  79. Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261 [CrossRef] [Google Scholar]
  80. Wang, H. S., Quanz, S. P., Yong, D., et al. 2022, MNRAS, 513, 5829 [NASA ADS] [Google Scholar]
  81. Wes McKinney 2010, in Proceedings of the 9th Python in Science Conference, eds. Stéfan van der Walt & Jarrod Millman, 56 [Google Scholar]

3

The latest version of the ARES code (ARES v2) can be downloaded at https://github.com/sousasag/ARES

All Tables

Table A.1

The main properties of the planets.

Table A.2

The stellar atmospheric parameters and present-day abundances of the planet host stars.

Table A.3

The ages and primordial abundances of the planet host stars.

Table A.4

The present-day and primordial compositions of planet building blocks.

Table A.5

The compositions of planets.

All Figures

thumbnail Fig. 1

Distribution of planets in the mass–radius diagram color-coded by their Teq. The curves display the expected mass–radius relationship for different compositions taken from Luo et al. (2024).

In the text
thumbnail Fig. 2

Difference between primordial and present-day abundances of iron as a function of age. Crosses show the data points, color-coded by Teff. Red error bars represent the uncertainties of present-day [Fe/H].

In the text
thumbnail Fig. 3

Difference between primordial and present-day abundance ratios of [Fe/Si] and [Fe/Mg] as a function of age. Crosses show the data points. Red error bars represent the uncertainties of present-day [Fe/Si] and [Fe/Mg].

In the text
thumbnail Fig. 4

Difference in iron-to-silicate mass (top) and water mass (bottom) fractions of planet building blocks determined from the present-day and primordial stellar abundances. The differences are plotted as a function of firon,prim star $f_{{\rm{iron,prim }}}^{{\rm{star }}}$ and fwater,prim star $f_{{\rm{water,prim }}}^{{\rm{star }}}$, respectively.

In the text
thumbnail Fig. 5

Scaled density of planets as a function of firon,primstar $f_{{\rm{iron,prim}}}^{{\rm{star }}}$ (top) and fwater,primstar $f_{{\rm{water,prim}}}^{{\rm{star }}}$ (bottom) color-coded by the equilibrium temperature of the planets. The positions of Earth is indicated with its symbols in black. The black solid lines represent the results of OLS regression without considering TOI-733 b and HD 23472 c. The red dashed-dotted line represent the visual separation between super-Earths and super-Mercuries. The black dashed line represent the result of the linear regression without considering TOI-733 b, HD 23472 c and the super-Mercuries located above the red line (p-value of about 10−5).

In the text
thumbnail Fig. 6

Iron-to-silicate and water mass fractions of planets and their host stars color-coded by the equilibrium temperature of the planets. The black solid line represents the results of OLS regression.

In the text
thumbnail Fig. 7

The ratio of iron-to-silicate mass fraction of planets and the planet building blocks as a function of planetary radius (top panel) and mass (bottom panel). The black solid lines represent the results of the linear regression. The symbols in the top panel are color-coded by the scaled density of the planets.

In the text
thumbnail Fig. 8

The Fe/Mg abundance ratio of planets and their host stars as a function of planetary radius. The black solid line in the plot represent the results of the OLS regression (p-value ~ 0.006) performed on the linearly scaled data.

In the text
thumbnail Fig. B.1

Iron-to-silicate mass fractions of planets and their host stars, with different lines representing the results of the linear fits discussed in the manuscript.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.