-
A Gap in the Densities of Small Planets Orbiting M Dwarfs: Rigorous Statistical Confirmation Using the Open-source Code RhoPop
Authors:
J. G. Schulze,
Ji Wang,
J. A. Johnson,
B. S. Gaudi,
R. Rodriguez Martinez,
C. T. Unterborn,
W. R. Panero
Abstract:
Using mass-radius-composition models, small planets ($\mathrm{R}\lesssim 2 \mathrm{R_\oplus}$) are typically classified into three types: iron-rich, nominally Earth-like, and those with solid/liquid water and/or atmosphere. These classes are generally expected to be variations within a compositional continuum. Recently, however, Luque & Pallé observed that potentially Earth-like planets around M d…
▽ More
Using mass-radius-composition models, small planets ($\mathrm{R}\lesssim 2 \mathrm{R_\oplus}$) are typically classified into three types: iron-rich, nominally Earth-like, and those with solid/liquid water and/or atmosphere. These classes are generally expected to be variations within a compositional continuum. Recently, however, Luque & Pallé observed that potentially Earth-like planets around M dwarfs are separated from a lower-density population by a density gap. Meanwhile, the results of Adibekyan et al. hint that iron-rich planets around FGK stars are also a distinct population. It therefore remains unclear whether small planets represent a continuum or multiple distinct populations. Differentiating the nature of these populations will help constrain potential formation mechanisms. We present the RhoPop software for identifying small-planet populations. RhoPop employs mixture models in a hierarchical framework and a nested sampler for parameter and evidence estimates. Using RhoPop, we confirm the two populations of Luque & Pallé with $>4σ$ significance. The intrinsic scatter in the Earth-like subpopulation is roughly half that expected based on stellar abundance variations in local FGK stars, perhaps implying M dwarfs have a smaller spread in the major rock-building elements (Fe, Mg, Si) than FGK stars. We apply RhoPop to the Adibekyan et al. sample and find no evidence of more than one population. We estimate the sample size required to resolve a population of planets with Mercury-like compositions from those with Earth-like compositions for various mass-radius precisions. Only 16 planets are needed when $σ_{M_p} = 5\%$ and $σ_{R_p} = 1\%$. At $σ_{M_p} = 10\%$ and $σ_{R_p} = 2.5\%$, however, over 154 planets are needed, an order of magnitude increase.
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
Fizzy Super-Earths: Impacts of Magma Composition on the Bulk Density and Structure of Lava Worlds
Authors:
Kiersten M. Boley,
Wendy R. Panero,
Cayman T. Unterborn,
Joseph G. Schulze,
Romy Rodrıguez Martınez,
Ji Wang
Abstract:
Lava worlds are a potential emerging population of Super-Earths that are on close-in orbits around their host stars with likely partially molten mantles. To date, few studies address the impact of magma on the observed properties of a planet. At ambient conditions magma is less dense than solid rock; however, it is also more compressible with increasing pressure. Therefore, it is unclear how large…
▽ More
Lava worlds are a potential emerging population of Super-Earths that are on close-in orbits around their host stars with likely partially molten mantles. To date, few studies address the impact of magma on the observed properties of a planet. At ambient conditions magma is less dense than solid rock; however, it is also more compressible with increasing pressure. Therefore, it is unclear how large-scale magma oceans affect planet observables, such as bulk density. We update ExoPlex, a thermodynamically self-consistent planet interior software, to include anhydrous, hydrous (2.2 wt \% H_2O), and carbonated magmas (5.2 wt\% CO_2). We find that Earth-like planets with magma oceans larger than \sim 1.5 R_{\oplus} and \sim 3.2 M_{\oplus} are modestly denser than an equivalent mass solid planet. From our model, three classes of mantle structures emerge for magma ocean planets: (1) mantle magma ocean, (2) surface magma ocean, and (3) one consisting of a surface magma ocean, solid rock layer, and a basal magma ocean. The class of planets in which a basal magma ocean is present may sequester dissolved volatiles on billion-year timescales, in which a 4 M_{\oplus} mass planet can trap more than 130 times the mass of water than in Earth's present-day oceans and 1000 times the carbon in the Earth's surface and crust.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
A Comparison of the Composition of Planets in Single- and Multi-Planet Systems Orbiting M dwarfs
Authors:
Romy Rodríguez Martínez,
David V. Martin,
B. Scott Gaudi,
Joseph G. Schulze,
Anusha Pai Asnodkar,
Kiersten M. Boley,
Sarah Ballard
Abstract:
We investigate and compare the composition of M-dwarf planets in systems with only one known planet (``singles") to those residing in multi-planet systems (``multis") and the fundamental properties of their host stars. We restrict our analysis to planets with directly measured masses and radii, which comprise a total of 70 planets: 30 singles and 40 multis in 19 systems. We compare the bulk densit…
▽ More
We investigate and compare the composition of M-dwarf planets in systems with only one known planet (``singles") to those residing in multi-planet systems (``multis") and the fundamental properties of their host stars. We restrict our analysis to planets with directly measured masses and radii, which comprise a total of 70 planets: 30 singles and 40 multis in 19 systems. We compare the bulk densities for the full sample, which includes planets ranging in size from $0.52 R_{\oplus}$ to $12.8R_\oplus$, and find that single planets have significantly lower densities on average than multis, which we cannot attribute to selection biases. We compare the bulk densities normalized by an Earth model for planets with $R_{p} < 6R_{\oplus}$, and find that multis are also denser with 99\% confidence. We calculate and compare the core/water mass fractions (CMF/WMF) of low-mass planets ($M_p <10 M_{\oplus}$), and find that the likely rocky multis (with $R_p <1.6 R_{\oplus}$) have lower CMFs than singles. We also compare the [Fe/H] metallicity and rotation period of all single versus multi-planet host stars with such measurements in the literature and find that multi-planet hosts are significantly more metal-poor than those hosting a single planet. Moreover, we find that host star metallicity decreases with increasing planet multiplicity. In contrast, we find only a modest difference in the rotation period. The significant differences in planetary composition and metallicity of the host stars point to different physical processes governing the formation of single- and multi-planet systems in M dwarfs.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
The Nominal Range of Rocky Planet Masses, Radii, Surface Gravities and Bulk Densities
Authors:
Cayman T. Unterborn,
Steven J. Desch,
Jonas Haldemann,
Alejandro Lorenzo,
Joseph G. Schulze,
Natalie R. Hinkel,
Wendy R. Panero
Abstract:
The two primary observable quantities of an exoplanet--its mass and radius--alone are not sufficient to probe a rocky exoplanet's interior composition and mineralogy. To overcome this, host-star abundances of the primary planet-building elements (Mg, Si, Fe) are typically used as a proxy for the planet's bulk composition. The majority of small exoplanet hosts, however, do not have available abunda…
▽ More
The two primary observable quantities of an exoplanet--its mass and radius--alone are not sufficient to probe a rocky exoplanet's interior composition and mineralogy. To overcome this, host-star abundances of the primary planet-building elements (Mg, Si, Fe) are typically used as a proxy for the planet's bulk composition. The majority of small exoplanet hosts, however, do not have available abundance data. Here we present the open-source ExoPlex mass-radius-composition solver. Unlike previous open-source mass-radius solvers, ExoPlex calculates the core chemistry and equilibrium mantle mineralogy for a bulk composition, including effects of mantle FeO content, core light elements and surface water/ice. We utilize ExoPlex to calculate the planetary radii, surface gravities and bulk densities for 10$^6$ model planets up to 2 R$_\oplus$ across these geochemistries, adopting the distribution of FGK stellar abundances to estimate of the range of bulk exoplanet compositions. We outline the $99.7\%$ distribution of radii, surface gravity and bulk densities that define planets as "nominally rocky." Planets outside this range require compositions outside those expected from stellar abundance data, likely making them either Fe-enriched super-Mercuries, or volatile-enriched mini-Neptunes. We apply our classification scheme to a sample of 85 well-resolved exoplanets without available host-star abundances. We estimate only 9 planets are within the "nominally rocky planet zone" at $>70\%$ confidence, while $\sim20\%$ and $\sim30\%$ of this sample can be reasonably classified as super-Mercuries or volatile-rich, respectively. Our results provide observers with a self-consistent way to broadly classify a planet as likely rocky, Mercury-like or volatile-enriched, using mass and radius measurements alone.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
A Reanalysis of the Composition of K2-106b: an Ultra-short Period Super-Mercury Candidate
Authors:
Romy Rodríguez Martínez,
B. Scott Gaudi,
Joseph G. Schulze,
Lorena Acuña,
Jared Kolecki,
Jennifer A. Johnson,
Anusha Pai Asnodkar,
Kiersten M. Boley,
Magali Deleuil,
Olivier Mousis,
Wendy R. Panero,
Ji Wang
Abstract:
We present a reanalysis of the K2-106 transiting planetary system, with a focus on the composition of K2-106b, an ultra-short period, super-Mercury candidate. We globally model existing photometric and radial velocity data and derive a planetary mass and radius for K2-106b of $M_{p} = 8.53\pm1.02~M_{\oplus}$ and $R_{p} = 1.71^{+0.069}_{-0.057}~R_{\oplus}$, which leads to a density of…
▽ More
We present a reanalysis of the K2-106 transiting planetary system, with a focus on the composition of K2-106b, an ultra-short period, super-Mercury candidate. We globally model existing photometric and radial velocity data and derive a planetary mass and radius for K2-106b of $M_{p} = 8.53\pm1.02~M_{\oplus}$ and $R_{p} = 1.71^{+0.069}_{-0.057}~R_{\oplus}$, which leads to a density of $ρ_{p} = 9.4^{+1.6}_{-1.5}$ $\rm g~cm^{-3}$, a significantly lower value than previously reported in the literature. We use planet interior models that assume a two-layer planet comprised of a liquid, pure Fe core and iron-free, $\rm MgSiO_{3}$ mantle, and we determine the range of core mass fractions that are consistent with the observed mass and radius. We use existing high-resolution spectra of the host star to derive Fe/Mg/Si abundances ([Fe/H]$=-0.03 \pm 0.01$, [Mg/H]$= 0.04 \pm 0.02$, [Si/H]$=0.03 \pm 0.06$) to infer the composition of K2-106b. We find that although K2-106b has a high density and core mass fraction ($44^{+12}_{-15}\%$) compared to the Earth ($33\%$), its composition is consistent with what is expected assuming that it reflects the relative refractory abundances of its host star. K2-106b is therefore unlikely to be a super-Mercury, as has been suggested in previous literature.
△ Less
Submitted 16 August, 2022;
originally announced August 2022.
-
Analytic Estimates of the Achievable Precision on the Physical Properties of Transiting Planets Using Purely Empirical Measurements
Authors:
Romy Rodriguez Martinez,
Daniel J. Stevens,
B. Scott Gaudi,
Joseph G. Schulze,
Wendy R. Panero,
Jennifer A. Johnson,
Ji Wang
Abstract:
We present analytic estimates of the fractional uncertainties on the mass, radius, surface gravity, and density of a transiting planet, using only empirical or semi-empirical measurements. We first express these parameters in terms of transit photometry and radial velocity (RV) observables, as well as the stellar radius $R_{\star}$, if required. In agreement with previous results, we find that, as…
▽ More
We present analytic estimates of the fractional uncertainties on the mass, radius, surface gravity, and density of a transiting planet, using only empirical or semi-empirical measurements. We first express these parameters in terms of transit photometry and radial velocity (RV) observables, as well as the stellar radius $R_{\star}$, if required. In agreement with previous results, we find that, assuming a circular orbit, the surface gravity of the planet ($g_p$) depends only on empirical transit and RV parameters; namely, the planet period $P$, the transit depth $δ$, the RV semi-amplitude $K_{\star}$, the transit duration $T$, and the ingress/egress duration $τ$. However, the planet mass and density depend on all these quantities, plus $R_{\star}$. Thus, an inference about the planet mass, radius, and density must rely upon an external constraint such as the stellar radius. For bright stars, stellar radii can now be measured nearly empirically by using measurements of the stellar bolometric flux, the effective temperature, and the distance to the star via its parallax, with the extinction $A_V$ being the only free parameter. For any given system, there is a hierarchy of achievable precisions on the planetary parameters, such that the planetary surface gravity is more accurately measured than the density, which in turn is more accurately measured than the mass. We find that surface gravity provides a strong constraint on the core mass fraction of terrestrial planets. This is useful, given that the surface gravity may be one of the best measured properties of a terrestrial planet.
△ Less
Submitted 22 January, 2021;
originally announced January 2021.
-
On the Probability that a Rocky Planet's Composition Reflects its Host Star
Authors:
J. G. Schulze,
Ji Wang,
J. A. Johnson,
B. S. Gaudi,
C. T. Unterborn,
W. R. Panero
Abstract:
The bulk density of a planet, as measured by mass and radius, is a result of planet structure and composition. Relative proportions of iron core, rocky mantle, and gaseous envelopes are degenerate for a given density. This degeneracy is reduced for rocky planets without significant gaseous envelopes when the structure is assumed to be a differentiated iron core and rocky mantle, in which the core…
▽ More
The bulk density of a planet, as measured by mass and radius, is a result of planet structure and composition. Relative proportions of iron core, rocky mantle, and gaseous envelopes are degenerate for a given density. This degeneracy is reduced for rocky planets without significant gaseous envelopes when the structure is assumed to be a differentiated iron core and rocky mantle, in which the core mass fraction (CMF) is a first-order description of a planet's bulk composition. A rocky planet's CMF may be derived both from bulk density and by assuming the planet reflects the host star's major rock-building elemental abundances (Fe, Mg, and Si). Contrasting CMF measures, therefore, shed light on the outcome diversity of planet formation from processes including mantle stripping, out-gassing, and/or late-stage volatile delivery. We present a statistically rigorous analysis of the consistency of these two CMF measures accounting for observational uncertainties of planet mass and radius and host-star chemical abundances. We find that these two measures are unlikely to be resolvable as statistically different unless the bulk density CMF is at least 40% greater than or 50% less than the CMF as inferred from the host star. Applied to 11 probable rocky exoplanets, Kepler-107c has a CMF as inferred from bulk density that is significantly greater than the inferred CMF from its host star (2$σ$) and is therefore likely an iron-enriched super-Mercury. K2-229b, previously described as a super-Mercury, however, does not meet the threshold for a super-Mercury at a 1- or 2- $σ$ level.
△ Less
Submitted 28 May, 2021; v1 submitted 17 November, 2020;
originally announced November 2020.
