The presence of rocky exoplanets with a large refractory carbon inventory is predicted by chemical evolution models of protoplanetary disks of stars with photospheric C/O >0.65, and by models studying the radial transport of refractory carbon. High-pressure high-temperature laboratory experiments show that most of the carbon in these exoplanets differentiates into a graphite outer shell. Our aim is to evaluate the effects of a graphite outer shell on the thermal evolution of rocky exoplanets containing a metallic core and a silicate mantle. We implement a parameterized model of mantle convection to determine the thermal evolution of rocky exoplanets with graphite layer thicknesses up to 1000 km. We find that, due to the high thermal conductivity of graphite, conduction is the dominant heat transport mechanism in a graphite layer for the long-term evolution (>200 Myr). The conductive graphite shell essentially behaves like a stagnant lid with a fixed thickness. Models of Kepler-37b (Mercury-size) and a Mars-sized exoplanet show that a planet with a graphite lid cools faster than a planet with a silicate lid, and a planet without a stagnant lid cools the fastest. A graphite lid needs to be approximately ten times thicker than a corresponding silicate lid in order to produce similar thermal evolution.

Carbon-enriched rocky exoplanets have been proposed around dwarf stars as well as around binary stars, white dwarfs and pulsars. However, the mineralogical make up of such planets is poorly constrained. We performed high-pressure high-temperature laboratory experiments ($P$ = 1$-$2 GPa, $T$ = 1523$-$1823 K) on carbon-enriched chemical mixtures to investigate the deep interiors of Pluto- to Mars-size planets the upper mantles of larger planets. Our results show that these exoplanets, when fully-differentiated, comprise a metallic core, a silicate mantle and a graphite layer on top of the silicate mantle. The silicate mineralogy (olivine, orthopyroxene, clinopyroxene and spinel) is largely unaffected by the amount of carbon. Metals are either two immiscible iron-rich alloys (S-rich and S-poor) or a single iron-rich alloy in the Fe-C-S system with immiscibility depending on the S/Fe ratio and core pressure. Graphite is the dominant carbon-bearing phase at the conditions of our experiments with no traces of silicon carbide or carbonates. If the bulk carbon content is higher than needed to saturate the mantle and the core, graphite would be in the form of an additional layer on top of the silicate mantle assuming differentiation. For a thick enough graphite layer, diamonds would form at the bottom of this layer due to high pressures. We model the interior structure of Kepler-37b and show that a mere 10 wt% graphite layer would decrease its derived mass by 7%, suggesting future space missions that determine both radius and mass of rocky exoplanets with insignificant gaseous envelopes could provide quantitative limits on their carbon content. Future observations of rocky exoplanets with graphite-rich surfaces would show low albedos due to the low reflectance of graphite. The absence of life-bearing elements other than carbon on the surface likely makes them uninhabitable.

Theoretical models predict the condensation of silicon carbide around host stars with C/O ratios higher than 0.65 (cf. C/O$_{\mathrm{Sun}}$ = 0.54), in addition to its observations in meteorites, interstellar medium and protoplanetary disks. Consequently, the interiors of rocky exoplanets born from carbon-enriched refractory material are often assumed to contain large amounts of silicon carbide. Here we aim to investigate the stability of silicon carbide in the interior of carbon-enriched rocky exoplanets and to derive the reaction leading to its transformation. We performed a high-pressure high-temperature experiment to investigate the reaction between a silicon carbide layer and a layer representative of the bulk composition of a carbon-enriched rocky exoplanet. We report the reaction leading to oxidation of silicon carbide producing quartz, graphite, and molten iron silicide. Combined with previous studies, we show that in order to stabilize silicon carbide, carbon saturation is not sufficient, and a complete reduction of Fe$^{2+}$ to Fe$^{0}$ in a planetary mantle is required, suggesting that future spectroscopic detection of Fe$^{2+}$ or Fe$^{3+}$ on the surface of rocky exoplanets would imply the absence of silicon carbide in their interiors.

More than a third of all exoplanets can be classified as super-Earths based on radius (1-2 $R_{\bigoplus}$) and mass (< 10 $M_{\bigoplus}$). Here we model mass-radius relations based on silicate mantle and iron core equations of state to infer to first order the structure and composition range of rocky super-Earths assuming insignificant gas envelopes. We develop a new equation of state of hexagonal close packed (hcp) iron for super-Earth conditions (SEOS) based on density functional theory results for pressures up to 137~TPa. A comparison of SEOS and extrapolated equations of state for iron from the literature reveals differences in density of up to 4% at 1~TPa and up to 20% at 10~TPa. Such density differences change the derived mass by up to 10\% for Earth-like super-Earths (core radius fraction of 0.5) and 20% for Mercury-like super-Earths (core radius fraction of 0.8). We find that the effect of temperature on mass (< 5%) is smaller than that resulting from the extrapolation of the equations of state of iron and lower mantle temperatures are too low to allow for rock and iron miscibility for R<1.75 $R_{\bigoplus}$. We find that modeling uncertainties dominate over observational uncertainties for many observed super-Earths. We illustrate these uncertainties explicitly for Kepler-36b with well-constrained mass and radius. Assuming a core composition of 0.8$\rho$ Fe (equivalent to 50 mol% S) instead of pure Fe leads to an increase of the core radius fraction from 0.53 to 0.64. Using a mantle composition of Mg$_{0.5}$Fe$_{0.5}$SiO$_3$ instead of MgSiO$_3$ leads to a decrease of the core radius fraction to 0.33. Effects of thermal structure and the choice of equation of state for the core material on the core radius of Kepler-36b are small but non-negligible, reaching 2% and 5%, respectively.