# Planetary Refractory Composition and Volatile Accretion into Gas Giants in the Protoplanetary Disks of the Sun and WASP-12

G. S. Pekmezci, O. Mousis, J. I. Lunine

► To cite this version:

G. S. Pekmezci, O. Mousis, J. I. Lunine. Planetary Refractory Composition and Volatile Accretion into Gas Giants in the Protoplanetary Disks of the Sun and WASP-12. The Astrophysical Journal, 2020, 904, 10.3847/1538-4357/abbd90. insu-03667034

HAL Id: insu-03667034

https://insu.hal.science/insu-03667034v1

Submitted on 6 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Planetary Refractory Composition and Volatile Accretion into Gas Giants in the Protoplanetary Disks of the Sun and WASP-12

G. S. Pekmezci<fnref n="1" />, O. Mousis<fnref n="2" />, and J. I. Lunine<fnref n="3" />

Received 2020 March 30; revised 2020 August 24; accepted 2020 September 16; published 2020 November 26

## Abstract

We present a detailed theoretical exploration of the refractory compositions and volatile enrichments of planets forming in protoplanetary disks with solar-like conditions. The two cases of the Sun and WASP-12 are studied due to the availability of spectral measurements and their known planets. The distribution throughout their disks of solid compounds with a wide range of volatilities is computed by a comprehensive chemical thermodynamics code. After the calculation of refractory compounds down to the water snowline, the compositional distributions are documented for planets generated in certain locations of protoplanetary disks depending on thermodynamic conditions. These results are referred to proposed bulk compositions for solar terrestrial planets, and for the core of the hot Jupiter WASP-12b. The material left over after the formation of rocky components is collected and treated in calculations to determine the abundances of fundamental volatile molecules in the outer regions of the disks. The distributions of planetesimal volatile composition are then altered for four different cases of the carbon-to-oxygen ratios, and for oxidizing and reducing conditions, in order to adjust the best fit for the accretion zone of Jupiter and WASP-12b. We compare the Jovian results to in situ atmospheric measurements from Jupiter's atmosphere. Overall, this study proposes a holistic approach to estimate possible planetary interior and envelope compositions from hot toward cold disk zones, along with the mass of planetesimals accreted into the envelopes of gas giants.

*Unified Astronomy Thesaurus concepts: Exoplanets (498); Hot Jupiters (753); Solar system planets (1260); Atmospheric composition (2120); Inner planets (797); Cosmochemistry (331); Extrasolar rocky planets (511); Extrasolar gas giants (509); Protoplanetary disks (1300); Astrochemistry (75)*

## 1. Introduction

The two most important parameters governing the overall composition of a planet are the composition of its parent star and the chemical and physical conditions of the protoplanetary disk where it formed and accreted (Guillot 2008; Eistrup et al. 2016; Mordasini et al. 2016, 2012; Cridland et al. 2017; Santos et al. 2017; Li et al. 2020b). Chemical diversity of distant condensates is, therefore, key to adapting solar system concepts such as the Kuiper Belt, snowline positions, icy satellites, and heavy-element enrichment in giant planet atmospheres to extrasolar disks, especially when considering varying abundance ratios of carbon to oxygen (Madhusudhan et al. 2011b; Madhusudhan 2012; Eistrup et al. 2016; Bitsch & Battistini 2020).

The C/O ratio is an important factor in the volatile ice composition of planetesimals, mainly determining the availability of oxygen in water formation and condensation, independent of the location of planet formation (Johnson et al. 2012; Moses et al. 2013; Cridland et al. 2017), as well as controlling the carbide and silicate distribution (Thiabaud et al. 2015). On the other hand, the C/O ratio in the atmosphere of a gas giant influences the equilibrium composition of important trace species such as CO, which can be indicators of vertical transport (Moses et al. 2013).

According to the core accretion scenario, the atmospheric C/O ratio of a giant planet might be different from that of its parent star (Mousis et al. 2011; Moses et al. 2013; Cridland et al. 2019a, 2019b), providing constraints on the composition of the circumstellar disk, location of formation, migration, and planetesimal contribution (Eistrup et al. 2018; Cridland et al. 2017;

Thiabaud et al. 2015; Moses et al. 2013; Madhusudhan 2012). The expanded version of the core accretion model with the inclusion of disk evolution and migration explains the basic properties of the gas giants (Mousis et al. 2009a, 2009b, 2009c), although the best estimates for the Jovian water abundance by Li et al. (2020a) may require more radical ideas involving gaseous accretion of certain heavy elements (Bitsch & Battistini 2020; Mousis et al. 2019).

Another important parameter for the nature of chemical processes within a disk is the environmental redox state (Johnson et al. 2012), defined by what molecular forms carbon assumes in the gas phase through the circumstellar disks, resulting in a range of reducing and oxidizing conditions. Condensation of solar system disk material beyond the snowline indicates a strong correlation between these two parameters and the carbon contribution in refractory material (Johnson et al. 2012).

In this study, accordingly, we first focus on the rocky composition of objects as a function of the distance from the central star and the formation conditions. Refractory species that build up the basic constituents of terrestrial planets are examined under certain thermodynamic conditions with the assumption of chemical equilibrium.

Thermochemical equilibrium is an approach that has long been in use to calculate gas and grain chemistry inside regions of circumstellar disks with sufficiently high temperatures, showing a general agreement with solar, meteorite, and planetary observations (Bitsch & Battistini 2020; Mousis et al. 2012). Although our results for refractory compositions are consistent with our current understanding of the formation and composition of terrestrial planets, this approach has its

<footnote n="1">¹Research and Education Center for the Superiorly Intelligent, 16120, Nilüfer, Bursa, Turkey; gspekmezci@roma2.infn.it</footnote>

<footnote n="2">²Aix Marseille Université, CNRS, LAM, Marseille, France</footnote>

<footnote n="3">³Cornell Center for Astrophysics and Planetary Sciences, Space Sciences Building, Cornell University, Ithaca, NY 14853, USA</footnote>
limitation, demonstrably so for close-in gas giants because they
are believed not to have formed in their current location but to
have migrated inward from distant locations in their stellar
systems (Lin & Papaloizou 1986; Ward 1997; Alibert et al.
2005a). Afterwards, we study the atmospheric enrichment in
our own Jupiter and in the hot Jupiter detected in the WASP-12
system, connecting the compositions of the giant planets to the
disk metallicity. The latter system is chosen since its stellar
chemical data are measured almost fully by Fossati et al.
(2010), which is rare for a host star. And Madhusudhan et al.
(2011a) reported the C/O ratio in the atmosphere of WASP-
12b to be at least twice that of its parent star, making it an
interesting case. Since its atmospheric composition is yet to be
measured directly, this study aims to provide a prediction.

The volatile accretion approach was successfully applied to
the solar system (Mousis et al. 2009a, 2009b, 2009c, 2011,
2012), and it is developed further here with a greater availability
of chemical species. Results of the study show consistency with
the mentioned papers. The connection between the C/O ratio
and issues such as the core accretion model, nebular oxidation
state, and planetary enrichment are examined and discussed in
the following sections.

**2. Methodology**

Generic planetary compositions are assumed to have two
possible components: refractory, i.e., the metallic and rocky
phase, and atmospheric, i.e., the volatile and gaseous phase.
The treatment for the refractory part of this work holds some
similarities to the studies of Bitsch & Battistini (2020), Carter-
Bond et al. (2012), and Bond et al. (2010b), with an
equilibrium assumption for the chemical composition of the
disks, where planetesimals are formed and accrete into larger
planetary embryos. The assumption of equilibrium condensa-
tion is already supported by various analyses of early chondritic
meteorites (Kuchner & Seager 2005; Bond et al. 2010a, 2010b;
Li et al. 2020b).

To examine the changes in types and abundances of
dominating refractory species, we employ a comprehensive
chemical thermodynamics code coupled with a simplified
kinetics scheme to calculate gas and ice abundances depending
on temperature, pressure, and elemental abundances, in a
similar way to the study of Lodders & Fegley (2002).

In order to model the atmospheres of Jupiter and WASP-12b,
a method is used where the remnant gas content in the disk
after the formation of refractory species is taken as the initial
parameter. The leftover abundances of five basic atmospheric
elements, namely oxygen, carbon, sulfur, phosphorus, and
nitrogen, are taken as the foundation of major volatiles in the
atmospheres, which are enriched by the planetesimal-derived
materials they accreted. The atmospheric compositions are
modeled for four different cases: where the disk conditions are
either oxidizing or reducing, and where the C/O ratio is either
stellar or unity. The fractionation and estimated masses of
important ices, and ice-to-rock percentages calculated under
these four combinations, are documented for these two gas
giants. The dependence of ice and rock enrichment on the C/O
ratio and redox states is evaluated for the final expected
composition of the planets.

## 2.1. Modeling the Disks

This work is an exploration of possible planetary composi-
tions with specific assumptions. A quiescent disk profile is
drawn in order to focus on the direct relations between
chemical composition and basic disk properties, as practiced
also by Bitsch & Battistini (2020), Price et al. (2020), Eistrup
et al. (2018), Moriarty et al. (2014), Bond et al. (2010b), and
Dodson-Robinson et al. (2009). A static model shall work fine
to approximate the abundances of disk species including water
and carbon gases, which are not significantly affected by disk
dynamics (Price et al. 2020).

Circumstellar disks are considered to be geometrically thin,
symmetric around their midplanes, and to remain in hydrostatic
equilibrium. And the physical conditions are regarded as stable
in their local values, as practiced also by Price et al. (2020).

Below are the disk factors that are disregarded in our study:

1. Stellar irradiation is not taken as an energy source, considering its negligible effect on the temperature of the midplane.

2. Vertical, radial, and relative motions of the gas and dust are not considered.

3. Modifications in midplane conditions, such as variation in mass accretion, are not considered to affect the condensation of materials throughout disks.

4. Viscosity, opacity, photoionization, stellar pollution, gravitational settling, turbulence, and gas giant migration are also not considered.

Nevertheless, it is possible to find various disk models in the
literature that studied some of the factors dismissed in our
work, as indicated in the following:

Disk evolution by Cridland et al. (2019b, 2020), Li et al. (2020b), Price et al. (2020), Eistrup et al. (2016, 2018), Alessi et al. (2017), Ali-Dib (2017), Helling et al. (2014), Moriarty et al. (2014), Mordasini et al. (2012). Viscosity by Price et al. (2020), Booth & Clarke (2018), Cridland et al. (2016), Mordasini et al. (2016), Piso et al. (2016), Wehrstedt & Gail (2003). Type I and/or type II migration by Cridland et al. (2017, 2019a), Ali-Dib (2017), Mordasini et al. (2016). Radial drift by Booth & Ilee (2019), Booth & Clarke (2018), Piso et al. (2016). Stellar evolution by Artur de la Villarmois et al. (2019), Santos et al. (2017). Ionization by Cridland et al. (2016), Eistrup et al. (2016, 2018). Radiative transfer by Cridland et al. (2017). Photoevaporation by Ali-Dib (2017). Accretion heating by Piso et al. (2016). Gravitational instability by Ilee et al. (2011). And non-equilibrium chemistry by Cridland et al. (2016).

Some of those studies are not concerned with chemical
kinetics for the simplicity of physical dynamics, such as the
models by Ali-Dib (2017), Mordasini et al. (2016), Ali-Dib
et al. (2014), and Marboeuf et al. (2014).

## 2.2. Data

The Sun and WASP-12 are chosen as case studies due to the
availability of chemical composition data. As summarized in
Table 1, they are both G-type dwarfs and have similar physical
properties.

Table 2 shows the data for their elemental number densities
relative to H₂ with the references. The Sun's photospheric
chemical composition is fully measured and summarized by
<table><caption>**Table 1**
Basic Stellar Parameters for the Sun and WASP-12</caption><thead><tr><th>Parameter</th><th>Sun</th><th>WASP-12<sup>a</sup></th></tr></thead><tbody><tr><td>Spectral type</td><td>G2V</td><td>G0V</td></tr><tr><td>Mass (M<sub>☉</sub>)</td><td>1</td><td>1.35</td></tr><tr><td>Radius (R<sub>☉</sub>)</td><td>1</td><td>1.57</td></tr><tr><td>Effective temperature (K)</td><td>5778</td><td>6300</td></tr></tbody></table>

**Note.**

<table><caption>**Table 2**
Elemental Abundances for Stars Relative to H<sub>2</sub></caption><thead><tr><th></th><th>Sun<sup>a</sup></th><th>WASP-12<sup>b</sup></th></tr></thead><tbody><tr><td>C/O</td><td>0.55</td><td>0.45</td></tr><tr><td>H</td><td>2</td><td>2</td></tr><tr><td>He</td><td>0.19</td><td>0.19</td></tr><tr><td>O</td><td>1.07 × 10<sup>-3</sup></td><td>1.59 × 10<sup>-3</sup></td></tr><tr><td>C</td><td>5.90 × 10<sup>-4</sup></td><td>7.10 × 10<sup>-4</sup></td></tr><tr><td>N</td><td>1.48 × 10<sup>-4</sup></td><td>2.24 × 10<sup>-4</sup></td></tr><tr><td>Mg</td><td>8.73 × 10<sup>-5</sup></td><td>2.03 × 10<sup>-4</sup></td></tr><tr><td>Si</td><td>7.10 × 10<sup>-5</sup></td><td>1.61 × 10<sup>-4</sup></td></tr><tr><td>Fe</td><td>6.93 × 10<sup>-5</sup></td><td>2.18 × 10<sup>-4</sup></td></tr><tr><td>S</td><td>2.89 × 10<sup>-5</sup></td><td>3.32 × 10<sup>-5</sup></td></tr><tr><td>Al</td><td>6.18 × 10<sup>-6</sup></td><td>4.69 × 10<sup>-6</sup></td></tr><tr><td>Ca</td><td>4.80 × 10<sup>-6</sup></td><td>1.44 × 10<sup>-5</sup></td></tr><tr><td>Na</td><td>3.81 × 10<sup>-6</sup></td><td>4.69 × 10<sup>-6</sup></td></tr><tr><td>Ni</td><td>3.64 × 10<sup>-6</sup></td><td>4.58 × 10<sup>-6</sup></td></tr><tr><td>Cr</td><td>9.57 × 10<sup>-7</sup></td><td>3.65 × 10<sup>-6</sup></td></tr><tr><td>Mn</td><td>5.90 × 10<sup>-7</sup></td><td>7.78 × 10<sup>-7</sup></td></tr><tr><td>P</td><td>5.64 × 10<sup>-7</sup></td><td>6.13 × 10<sup>-7</sup></td></tr><tr><td>K</td><td>2.35 × 10<sup>-7</sup></td><td>1.52 × 10<sup>-6</sup></td></tr><tr><td>Co</td><td>2.14 × 10<sup>-7</sup></td><td>2.09 × 10<sup>-7</sup></td></tr><tr><td>Ti</td><td>1.95 × 10<sup>-7</sup></td><td>5.72 × 10<sup>-7</sup></td></tr><tr><td>Zn</td><td>7.96 × 10<sup>-8</sup></td><td>9.57 × 10<sup>-8</sup></td></tr><tr><td>Cu</td><td>3.40 × 10<sup>-8</sup></td><td>3.10 × 10<sup>-8</sup></td></tr><tr><td>V</td><td>1.87 × 10<sup>-8</sup></td><td>4.81 × 10<sup>-8</sup></td></tr><tr><td>Sc</td><td>3.10 × 10<sup>-9</sup></td><td>5.64 × 10<sup>-9</sup></td></tr><tr><td>Sr</td><td>1.63 × 10<sup>-9</sup></td><td>1.90 × 10<sup>-8</sup></td></tr></tbody></table>

**Notes.**

Asplund et al. (2009). The data for WASP-12 (Fossati et al. 2010) provide all elements necessary for the study except phosphorus. The elements Ca, Cr, Fe, Mg, Si, Sr, Ti, and V, each documented with their two lines, are assumed to have their average values. H and He are kept fixed to the solar values for WASP-12 as well. For the unknown phosphorus abundance of WASP-12, the value is scaled from the solar using the metallicity of the star (i.e., Fe abundance).

Since it is currently not possible to observe the chemical composition of disk midplanes (Eistrup et al. 2018), protoplanetary disks are assumed to directly reflect the chemistry of their parent stars while modeling the elemental composition of their planets, as also done by Bitsch & Battistini (2020) and Moriarty et al. (2014). Fortunately, it is possible to draw correlations between the refractory elements of the Sun's photosphere, meteorites, and rocky planets (Santos et al. 2017).

Moreover, the data for WASP-12 stand for its protostellar values as well, since detailed investigations about the

compositions of exoplanetary disks are not yet available. This assumption has been employed before, at least for the solar case (Johnson et al. 2012; Alibert et al. 2005a; Bond et al. 2010a).

## 2.3. The Software

The solid phases, which are silicate and metal compounds and formed in the disks at equilibrium composition, are calculated with HSC Chemistry 7.1. The program has a sophisticated set of algorithms automatically utilized within the framework of a thermochemical database containing enthalpy, entropy, and heat capacity data for 25,000+ chemical compounds (Roine 2006).

This program provides reliable computations for the effects of sundry variables on a chemical system at equilibrium (Roine 2006), and documents the equilibrium composition with a minimized Gibbs free energy for the defined system (Bond et al. 2010a). Once the user sets the quantity of raw materials and system conditions, the code will calculate the abundances of products to come out according to the initial conditions and chemical reactions to take place (Roine 2006). Studies concerned with solar nebula chemistry or supernova stellar outflows, whose results were correlated with mineralogy observations of interplanetary dust particles, have used the HSC software (e.g., Bond et al. 2010a; Moriarty et al. 2014).

## 2.4. Refractory Phase

After being confident with the limitations and performance capabilities of the program, the largest list of species possible is constituted using all the combinations found in the code's database for the universally most abundant and important refractory-forming elements: Al, C, Ca, Co, Cr, Cu, Fe, K, Mg, Mn, N, Na, Ni, O, P, S, Sc, Si, Sr, Ti, V, Zn, along with H and He. The input elements are entered only in the gas phase, and the system is restricted to these elements only. Inside the pool of possibilities, the alternatives listed below are included as output compounds:

1. Substances with the same composition, but different molecular or crystal structure (e.g., CuFeS<sub>2</sub>: both as ferrous sulfide and as chalcopyrite)

2. Various stoichiometric distributions for a certain species (e.g., FeSi<sub>2</sub>, FeSi<sub>2.33</sub>, FeSi<sub>2.43</sub>)

3. Compounds with different stoichiometric numbers (e.g., 3MgO.2SiO<sub>2</sub>.2H<sub>2</sub>O and 3MgO.4SiO<sub>2</sub>.H<sub>2</sub>O)

4. Organic species carrying up to three carbon atoms

5. Variations of compounds with an exact chemical formula (e.g., allotropic forms, organic isomers, excited states, duplets, and triplets).

The list below gives the species that are excluded due either to the numeric limitations of the code or to the conditions not considered for stellar disks, such as an ionized or aqueous environment. Liquids are excluded since they are generally found to occupy very small fields of stability in the disks, and they do not significantly affect condensation sequences (Bond et al. 2010a).

1. Both gaseous and aqueous ions

2. Radicals

3. Liquid compounds and species

4. Organic compounds containing three or more carbon atoms

5. Alloy combinations.

<footnote><sup>a</sup> Hebb et al. (2009).</footnote>

<footnote><sup>a</sup> Asplund et al. (2009).</footnote>

<footnote><sup>b</sup> Fossati et al. (2010).</footnote>
The complete list of input species, whose total number reached 1225, is given in Table A1 in the Appendix after alterations mentioned above are removed in order to avoid repetitions.

Temperature is taken within the range 100–2000 K, and pressure is kept fixed at $10^{-4}$ bar, a traditional value for solar nebula studies, because it is the characteristic total pressure at or around 1 au from the Sun (Cassen 1996; Fegley 2000; Lodders 2003; Bond et al. 2010b).

We use a solar temperature-pressure profile for WASP-12 since it is the same stellar type as the Sun. In this way, the same T-P conditions for both stars allow ready comparison of changes in chemical abundances between the two. Furthermore, as noted by Lodders & Fegley (2002), choosing a general T-P profile allows one to apply the resulting values of abundances to any other disks or bodies, and to determine the chemical compositions for any desired T-P range.

When a reaction is required to be calculated for a temperature value that is not available in the database, the code extrapolates the equation of the equilibrium constant to the required temperature. In order to avoid unexpected results, temperature values lower than 100 K are not taken into consideration.

We assume that the compounds appearing in temperature zones lower than 400 K in our results may be kinetically inhibited over the lifetime of the disk.

## 2.5. Volatile Phase

In the core accretion model, the rapid accretion of gases and gas-coupled solids takes place after the formation of solid cores by collisional accretion of planetesimals in disk zones beyond the water snowline (Mousis et al. 2009c). Accordingly, the material that does not get trapped in cores is free to enrich the gas in heavy elements, or to be trapped by clathration or adsorption into ice to enrich the envelopes of the giant planets at a later stage of the process.

Similar to the approach of Johnson et al. (2012), the volatile ice abundances in the envelopes of the planets are calculated taking the C, N, O, P, and S abundances remaining in the gas, following the formation of compounds consistent with mass conservation. This method seems to agree with the measured enrichments in the atmosphere of Jupiter (Alibert et al. 2005b).

Abundances of volatiles accreted into the gas giants are calculated by a FORTRAN code whose basic theory is explained briefly below, following the assumptions and formulations provided by Mousis et al. (2009c).

For the accretion disk, the turbulent model provided by Alibert et al. (2005a) is taken as a reference. Condensation leading to the formation of clathrates and pure condensates in the accretion zones of proto-giant planets is treated by assuming an initially homogenous gas phase of water consistent with the elemental stellar composition.

The actual clathration efficiency is currently unknown for the solar nebula. Therefore, among the set of trial runs from 0% to 100%, complete clathration is chosen for simplicity and assumed to be constant throughout the disk.

Hydrate, clathrate, and pure condensate stability curves, which were derived from the experimental work of Lunine & Stevenson (1985) and laboratory data of Lide (2002), are used to estimate the volatiles trapped in icy planetesimals for the temperature and pressure conditions at the distance of Jupiter.

Under these assumptions, the mass ratios of different ices can be estimated with respect to water in the planetesimals. The

mass ratio of the volatile relative to that of water is given by Mousis & Gautier (2004) as

$$m_i = \frac{X_i}{X_{H_2O}} \times \frac{\Sigma(r; T_i, P_i)}{\Sigma(r; T_{H_2O}, P_{H_2O})} \quad (1)$$

Here, X stands for the mass mixing ratio of the associated compound relative to H₂ in the disk. Σ is the disk surface density at the distance r, at the clathration epoch of the associated species, and at the epoch of water condensation.

$$M_i = m_i / \Sigma m_i \quad (2)$$

$$\sum M_i = 1. \quad (3)$$

The mass fraction $M_i$ is calculated from $m_i$, relative to the total volatile compounds involved in the formation of icy planetesimals.

Despite it being an oversimplification, it is assumed that all of the vapor ends up in the clathrate for every species, in order to follow the due implications.

The ranges of volatile abundance in planetary envelopes are estimated from the study of Mousis et al. (2011). The ices formed beyond the snowline carry an important fraction of volatiles affecting the O and C abundances by vaporization when they enter into the planetary envelopes (Madhusudhan et al. 2011b). Icy contributions are shaped by the initial elemental composition of the disk gas phase, rather than the location of their formation, or the thermodynamic conditions selected, rendering all planetesimals similar (Mousis et al. 2009b, 2012; Madhusudhan et al. 2011b). Therefore, the gas-phase composition is defined a priori, assuming that the protoplanetary disk has the stellar elemental abundances. The volatiles are then assumed to be entirely captured in the icy planetesimals (Johnson et al. 2012). While this is a simplifying assumption, it is not bad for giant planets formed beyond the amorphous-to-crystalline transition region. Mousis et al. (2020), on the other hand, provided a detailed study on the fractionation that occurs upon conversion of amorphous ice to crystalline ice.

Since the approach of the extended core scenario worked fine for the volatile enrichment in Jupiter's atmosphere (Mousis et al. 2009c, 2011), this work extrapolates that method to the hot Jupiter WASP-12b, seeking the most suitable matching of volatile abundances in its atmosphere. The abundances of species are evaluated at a fixed temperature of 300 K. This is a value higher than the condensation temperature of water in the protoplanetary disks, yet still cold enough for many volatile species to be trapped in meteorites. On the other hand, the equilibrium reactions involving refractory compounds can safely be regarded as quenched at a temperature of 300 K under stellar disk conditions. Hence, it is acceptable to assume that the leftover abundances of elements of concern would be more or less constant for the zones with temperatures lower than the quench temperature of their equilibrium reactions.

After collecting the results from the HSC code for refractory species, the ones with a molar fraction above $10^{-10}$ relative to H₂ are selected. The five primal biogenic elements O, C, N, S, and P, whose initial abundances are documented in Table 2, first incorporate into solids, and then their residual abundances are determined and given in Table 3. These elemental abundances are used to compute the main volatile molecules residing in the disks before condensing or being trapped into clathrates, as described in, e.g., Mousis et al. (2011, 2012).
<table><caption>**Table 3**
Remaining Molar Fractions of Five Elements Relative to H₂ Considered for Volatile Accretion after the Consumption by Refractory Phases in Stellar Disks</caption><thead><tr><th>Element/H<sub>2</sub></th><th>Sun</th><th>WASP-12</th></tr></thead><tbody><tr><td>O</td><td>6.88 × 10<sup>-4</sup></td><td>6.16 × 10<sup>-4</sup></td></tr><tr><td>C</td><td>5.92 × 10<sup>-4</sup></td><td>7.10 × 10<sup>-4</sup></td></tr><tr><td>N</td><td>1.49 × 10<sup>-4</sup></td><td>2.24 × 10<sup>-4</sup></td></tr><tr><td>P</td><td>5.64 × 10<sup>-7</sup></td><td>6.13 × 10<sup>-7</sup></td></tr><tr><td>S</td><td>1.45 × 10<sup>-5</sup></td><td>1.66 × 10<sup>-5</sup></td></tr></tbody></table>

To maintain the simplicity of the model, physical and chemical factors such as the angle and wavelength of stellar irradiation exposure, and photochemistry, which are known to influence the atmospheres (Visscher et al. 2006; Moses et al. 2011, 2013), are not taken into account.

## 2.6. Redox State

The alteration of the conditions, to be oxidizing or reducing, depends on the oxygen abundance. We consider the redox state of the gas phase only. Operationally, the circumstances will be oxidizing when the dominating gas in the disk is CO, and they will be reducing when the major gas is CH₄.

In the oxidizing case, O, C, N, S, and P are limited to the following molecular species: H₂O, CO, CO₂, CH₃OH, CH₄, N₂, NH₃, H₂S, and PH₃ (Madhusudhan et al. 2011a, 2011b; Mousis et al. 2009a, 2009b, 2012; Eistrup et al. 2016). The ratios between carbon-bearing species are fixed to the following values: CO/CO₂/CH₃OH/CH₄ = 70/10/2/1 (Johnson et al. 2012; Mousis et al. 2009a, 2009b, 2012; Madhusudhan et al. 2011a).

As for the nitrogen-bearing components, N₂ takes the majority with the ratio N₂/NH₃ = 10/1. Finally, sulfur-bearing components are treated with the ratio H₂S/S = 1/2 (Pasek et al. 2005; Mousis et al. 2009a, 2009b).

In the reducing case, carbon is considered to be consumed only by CH₄, which leaves no carbon for CO, CO₂, and CH₃OH. Nitrogen gas is now the minor nitride species with a ratio N₂/NH₃ = 1/10. The sulfur ratio remains the same as in the oxidizing case, H₂S/S = 1/2. Finally, in both cases, all phosphorus is found only in the form of PH₃.

The ratios given above are compatible with thermochemical calculations for the Sun's protoplanetary disk, consistent with the measurements of the interstellar medium including both solid and gas contributors, and close to the measurements made on comets (Frerking et al. 1982; Ohishi et al. 1992; Gibb et al. 2000; Ehrenfreund & Schutte 2000; Mousis et al. 2012).

## 2.7. C/O Ratio

Since the planetary models depend on stellar chemical compositions, the initial C/O ratio is important to a certain extent (Bitsch & Battistini 2020; Cridland et al. 2019b; Fortney 2012), as it is in the case of the equilibrium composition of a hot Jupiter, which is known to show sensitivity to the C/O ratio (Moses et al. 2013).

Observations show that the atmospheres of such planets may have C/O ratios higher than their host stars (Madhusudhan et al. 2011b). As one possible reason for this, Mousis et al. (2009c) proposed that the formation zone within the disk might have been deprived of oxygen, and carbon-rich planetesimals were accreted into their envelopes. Therefore, the C/O ratio for

both stars is converted to 1 by adjusting the oxygen abundance down to the amount of carbon. This method is preferred by Lodders (2010) as well because it does not add to the existing amount of carbon in the environment.

The leftover oxygen abundance after the formation of volatile species is calculated for both disks when the C/O ratio is unity. The initial amount of oxygen input to the software is adjusted to a value that will yield the targeted leftover oxygen abundance for the volatile compounds in each case. When the targeted oxygen abundance is reached, rocky phase compositions are calculated and documented for the Sun and WASP-12, consistent with the redox state.

For the case where C/O is 1, the O abundance is highly depleted with respect to the initial ratio defined by the solar elemental abundance, therefore oxygen is shared only among the carbonic species, leaving no space for water formation (Lodders 2010; Mousis et al. 2012). Since the elemental quantities are derived after the removal of refractory species in the disks, these C/O ratios refer to the gas-phase C/O ratios directly.

## 2.8. Comparison with Jovian Measurements

Following Mousis et al. (2012), once the composition of planetesimals is calculated for the oxidizing and reducing cases, the mass of heavy elements found in Jupiter's envelope is adjusted so that it can provide the best fit for volatile abundances measured by the Galileo probe, as reported by Wong et al. (2004) and Fletcher et al. (2009). When the abundances of these molecules are determined, the remaining oxygen content is left for water formation (Mousis et al. 2009c).

## 2.9. Enrichment Factors

Enrichment from planetesimals to the planets is considered separately for the ice and rocks. In order to estimate the mass of enrichment made in ices, the species H₂O, CH₄, CO, CO₂, and CH₃OH are taken into account. However, for the enrichment in rocks, all species that have molar fractions above 10⁻¹⁰ relative to H₂ are considered to contribute to the bulk content of the planet. The results are provided below.

# 3. Refractory Phase

The refractory composition throughout the disks is documented only for the species surpassing a fraction of 5% over the total mass, in order to make the plots interpretable. The distribution of species is presented as a function of temperature. The pressure is constant at 10⁻⁴ bar for each case, and components are given as percentage fractions over a total number of 840 refractory species, within the temperature range 100–2000 K.

Generally speaking, specific temperatures associated with refractory species are found here to be slightly lower by about 50–100 K than the results given by previous studies on solar nebula chemistry, such as those from Bond et al. (2010b), Fegley (2000), and Prinn (1993). This difference is mostly observed for the species appearing at temperatures lower than 700 K. The reason might depend on the preferences of the temperature-pressure profile of the solar disk, the data source of solar chemical composition, and/or the program code used for calculations.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The leading species of expected refractory composition in the disk of (a) the Sun and (b) WASP-12 on the basis of temperature. Dashed lines stand for alternative molecular combinations of species with a curve of the same color.</figcaption></figure>

The temperature values attributed to species are also different in this work from the others, especially when modeling the solar terrestrial planets in the following sections. Temperatures considered are those for which the aforementioned species are found at their maximum abundances, not the maximum temperatures at which the species form, which was how it was generally done in previous studies.

For example, Moriarty et al. (2014) noted that oxygen-bearing molecules appear below 1300 K, but our results show that oxygen can safely survive up to 1730 K by combining with aluminum, and even above 1900 K if incorporated in CaTiO₃ and K₂TiO₃.

Figure 1 depicts the dominant solid-phase species calculated with stellar elemental abundances. A zonal division between diverse refractory compounds formed out of certain elements is clearly observable throughout the disks, as noted also by Bond

et al. (2010b). The disks can be imagined as zonal divisions where some certain compounds dominate inward, while others dominate outward. Toward cooler regions, the formation of compounds shifts from refractory to volatile products (Lodders 2003; Bond et al. 2010a).

### 3.1. Magnesium-bearing Species

The dominant magnesium compound in the disks is Mg₄Si₆O₂₁H₁₂ (sepiolite). Being a magnesium silicate as well as a clay mineral, it contributes most to crustal structures formed at temperatures of up to 500 K. Below 400 K, it comprises 40%–50% of overall rocky compounds but is not available above 520 K.

Where sepiolite begins diminishing, magnesium coalesces with another silicate compound Mg₂SiO₄, which emerges above 400 K in two different molecular alignments of magnesium orthosilicate: forsterite and the general form, also
known as 2MgO·SiO$_2$. In the disk compositions, forsterite can reach up to 20% around 500 K, then gradually perishes around 1400 K.

MgSiO$_3$ (magnesium metasilicate) contributes in the disks around 520 K, exhibiting a relatively stable pattern between 550 K and 1300 K. When above 1400 K, magnesium metasilicate is lost quickly, along with Mg$_2$SiO$_4$ isomers. These two are often regarded as the major magnesium and silicon carriers in the disks (Lodders 2003), forming the majority of mantle material for rocky planets (Alessi et al. 2017).

Apart from MgO, magnesium tends to combine with silicon, and apart from disthene (Al$_2$SiO$_5$), silicon tends to combine with magnesium only.

### 3.2. Iron-bearing Species

The formation of magnetite (Fe$_3$O$_4$) begins around 320 K, whereas Fogley (1988) finds it to be 400 K for the solar nebula. With higher temperatures, Fe$_3$O$_4$ loses its oxygen, and free iron allotropes dominate in the exact trend. Iron remains roughly stable from 700 to 1300 K, making a contribution of about 30% for the solar case and 40% in WASP-12.

Metallic iron condensation coincides with magnesium silicates, as also noted by Lodders (2003). However, she mentions this temperature interval to be as narrow as 100 K, whereas it is about 1000 K here.

Pyrrhotite (FeS; iron 0.877 sulfide) appears in the solar case only, with a maximum 7% overall up to around 580 K, as suggested also by Lodders (2003). This temperature for FeS is given as 687 K by Prinn (1993), and 710 K by Fogley (2000). Fe$_3$O$_4$ and FeS contribute 70%-80% of the core materials in terrestrial planets (Alessi et al. 2017).

Finally, wüstite (FeO), also known as iron II oxide, is produced with a peak around 390 K for 5% in the solar case, and between 320-400 K with a peak of 9% in the WASP-12 case.

### 3.3. Calcium-bearing Species

The only calcium-carrying species in cold zones is calcium triiron pentaoxide (CaFe$_3$O$_5$), which appears in regions colder than 350 K with a contribution of 8-10%. In the interval 1750-1500 K, variations of calcium oxides combined with aluminum and silicates appear in the forms CaAl$_4$O$_7$, Ca$_2$SiO$_4$, Ca$_2$Al$_2$SiO$_7$, Ca$_3$SiO$_7$, and Ca$_3$SiO$_5$. Between 1600 and 1500 K, calcium silicates cover the entire refractory distribution in the disks along with K$_2$TiO$_3$ in WASP-12. For the solar case, calcium combines with titanium oxide as Ca$_2$TiO$_3$ above 1600 K, being the only rocky species in the range 1820-1740 K. Almost all of the calcium species mentioned here are confirmed in the study of Lodders (2003).

### 3.4. Aluminum-bearing Species

Aluminum oxide is observed in various forms of Al$_2$O$_3$. They can combine with calcium to form CaAl$_2$O$_7$, CaAl$_2$O$_4$, or Ca$_2$Al$_2$SiO$_7$ between 1750 and 1500 K. One of the major calcium-aluminum oxides that appeared in the study of Lodders (2003) was hibonite (CaAl$_{12}$O$_{19}$), which is not included in the calculations here.

Excluding CaFe$_3$O$_5$, calcium and aluminum appear only in hot regions for a narrow range of 200 K, as supported also by Miyazaki & Korenaga (2020). Above 1750 K, refractory components of disks are dominated only by CrC$_{24}$H$_{36}$ and K$_2$TiO$_3$ or CaTiO$_3$.

Miyazaki & Korenaga (2020) modeled a disk chemistry where Al, Ca, and Mg minerals drift toward inner regions from distant zones. Our study, however, shows that those minerals can naturally be produced at high temperatures without the necessity of being dragged from outer zones.

Although the order of abundances might change in two disks, and a few species might appear or disappear in one of them, the general behaviors of species are quite similar. Close abundances of elements result in similar disk refractory distributions, as also confirmed by Lodders (2010).

The majority of the silicates are already known to be the iron- and magnesium-bearing ones in the solar nebula (Léger et al. 2004). The species found to have a mass fraction of at least 5% in the disks of the Sun and WASP-12 but not mentioned in similar work by Bond et al. (2010a, 2010b) are listed in Table 4. Both studies worked with the same temperature and pressure values, therefore the only reason that Bond et al. (2010a) did not observe these species is that they did not include them in their calculations. Interestingly enough, the three species, namely Al$_2$O$_3$, CaTiO$_3$, and Mg$_2$SiO$_4$, are not taken into consideration by Bond et al. (2010a), but they were mentioned in an earlier study by Prinn (1993), as well as here being highly refractory species that are expected to survive in hot regions of the solar disk.

Table 5 lists the species that appeared in the results of Bond et al. (2010a) but do not reach a mass fraction of 5% in the result here, despite being included. Accordingly, solid component distributions in the disks documented by Bond et al. (2010a) and this work differ from each other fundamentally for hot and cold extremes of the temperature regime of interest. Still, these two studies are compatible in respect of the species dominating in a bigger part of the disks, which are pyroxene (MgSiO$_3$), olivine (Mg$_2$SiO$_4$), and elemental iron (Fe). This conclusion, as well as the iron peak around 1500 K, is also supported by Lodders (2003).

Another example is the study of Alessi et al. (2017), which documented planetary solid mass for five super-Earths, with the maximum mineral distributions of FeS, Fe$_3$O$_4$, FeAl$_2$O$_4$, CaMgSi$_2$O$_6$, MgSiO$_3$, and Mg$_2$SiO$_4$. However, only the relative abundances of the last two are comparable with our results for the Earth.

Since the redox state of the disks is adjusted with the oxygen abundance, there is only one set of refractory results for neutral, oxidizing, and reducing conditions when the C/O ratio is stellar. Oxygen-bearing species increase in abundance due to the increased availability of oxygen when the circumstances are oxidizing, and endure for a broader temperature region. The influence of reducing conditions is even less important, diminishing the abundances of a few oxidized species and increasing the fraction of elemental iron, only.

## 4. Planetary Models

In this section, further implications from the rocky species are presented for individual planets. Bulk compositions are predicted for solar terrestrial planets. For the case of WASP-12, we present a range of possible refractory compositions for planetesimals accreted in WASP-12b.

Since many types of planets, especially terrestrial ones, are thought to be common in various extrasolar planetary systems (Bond et al. 2010a, 2010b), and since planetary bulk compositions and mineral ratios of planetesimals, meteorites, and asteroids provide a great amount of information on their
<table><caption>**Table 4**
The Dominant Rocky Species Found in This Study but Not by Bond et al. (2010a)</caption><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>CaAl<sub>4</sub>O<sub>7</sub></td><td>Mg(OH)<sub>2</sub></td></tr><tr><td>Al<sub>2</sub>SiO<sub>5</sub></td><td>CaFe<sub>3</sub>O<sub>5</sub></td><td>Mg<sub>2</sub>SiO<sub>4</sub></td></tr><tr><td>Ca<sub>2</sub>Al<sub>2</sub>SiO<sub>7</sub></td><td>CaSiO<sub>3</sub></td><td>Mg<sub>4</sub>Si<sub>6</sub>O<sub>21</sub>H<sub>12</sub></td></tr><tr><td>Ca<sub>2</sub>SiO<sub>4</sub></td><td>CrC<sub>24</sub>H<sub>36</sub></td><td>MgO</td></tr><tr><td>Ca<sub>3</sub>Si<sub>2</sub>O<sub>7</sub></td><td>FeAl<sub>2</sub>O<sub>4</sub></td><td>MnC<sub>2</sub></td></tr><tr><td>Ca<sub>3</sub>SiO<sub>5</sub></td><td>FeO</td><td>NaH<sub>2</sub>SiO<sub>4.7</sub>H<sub>2</sub>O</td></tr><tr><td>CaAl<sub>2</sub>O<sub>4</sub></td><td>K<sub>2</sub>TiO<sub>3</sub></td><td>SiO<sub>2</sub></td></tr></table>

<table><caption>**Table 5**
The Rocky Species Found in the Results of Bond et al. (2010a) but Not in This Study</caption><tr><td>C</td><td>FeSiO<sub>3</sub></td><td>MgS</td></tr><tr><td>CaMgSi<sub>2</sub>O<sub>6</sub></td><td>Mg<sub>3</sub>Si<sub>2</sub>O<sub>5</sub>(OH)<sub>4</sub></td><td>NaAlSi<sub>3</sub>O<sub>8</sub></td></tr><tr><td rowspan="2">CaS</td><td>MgAl<sub>2</sub>O<sub>4</sub></td><td rowspan="2">SiC</td></tr><tr><td>Fe<sub>2</sub>SiO<sub>4</sub></td></tr></table>

places of origin, the results of this section are expected to be of use in interpreting exoplanet observations.

This approach may very well work also in the reverse direction since disk composition is correlated to the radial distance (Bond et al. 2010b). The mineral composition of a planetary body can provide a constraint on where it formed in the disk (Piso et al. 2016).

In the distributional charts of Figures 2–6 below, magnesium-bearing species are presented in green shades, iron-bearing ones in red shades, compounds with calcium in blue shades, and those with aluminum in brown shades. Other refractory species with titanium, chromium, or nickel are depicted in gray shades, and the rest of the materials with smaller fractions are collected under “others” and colored violet.

### 4.1. The Solar System

To present the refractory bulk composition for the solar planets, an approach is introduced on the assumption of formation in situ, because it is believed that planetary locations are important in determining their compositions (Lewis 1974; Kuchner & Seager 2005; Santos et al. 2017). This assumption is good for the terrestrial planets because the bulk of the accreted material from which they formed is probably local, and they probably formed in regions close to the Sun (Bitsch & Battistini 2020), whereas the giants are likely to have migrated before settling into their current orbits (Lin 2008). As such, the calculations are done down to 100 K only, but the giant planets fall beyond this temperature range.

First, estimated bulk compositions for the solar system's planets are collected from previous works, and the temperature value offering the closest ratios of species is taken as a reference to the birthplace temperature for the mentioned planet. Since refractory phases of the planetary bulks are not known in detail, the average value, or at least the existence of a characteristic species, is taken as a parameter to indicate the temperature of the location where it is formed and to trace the distribution of other species in its vicinity.

Among the species considered, those reaching a mass fraction down to 0.1% are documented by mass fraction and in a comparative manner between two different temperature profiles for each planet.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Mercury's refractory composition calculated for the temperature where elemental iron is found at the maximum abundance.</figcaption></figure>

A planetary formation model introduced by Bond et al. (2010a) tried to reproduce the Earth, but it failed to produce the terrestrial composition. The reason was blamed on the computer code they used, for not forming the planet at the distance where the most likely composition was found. Therefore, they suggested that the Earth might have migrated by at least 0.4 au. Migration is not necessarily predicted in the evolution of inner planets (Bitsch & Battistini 2020), and this scenario is not required for the Earth according to the approach offered in this work either.

Similarly, Elser et al. (2012) applied chemical equilibrium and dynamical simulations to circumstellar disk models; however, their models could not reproduce the bulk composition of terrestrial planets.

Ronco et al. (2015) also simulated solar rocky planets, but since they documented the compositions over elemental distributions only, their results are not comparable with ours.

#### 4.1.1. Mercury

Explanations for the inner structure of Mercury are not accurate. Given its high uncompressed density of $5.3 \text{ g cm}^{-3}$, iron-based species are assumed to be dominant (Goettel 1988). The traditional assumption for Mercury's bulk composition is that it consists of 65%–70% of metallic and 25%–30% of silicate materials (Kozlovskaja & Pechernikova 1992; Strom & Sprague 2003).

Figure 1(a) shows that the maximum elemental abundance of iron is found at the temperature zone of 1420 K, with a total fraction of 55%. The refractory composition for the region where elemental iron is found at its maximum is presented in Figure 2. The next major compounds are aluminum, calcium, and magnesium silicates such as Al<sub>2</sub>SiO<sub>5</sub>, CaSiO<sub>3</sub>, MgSiO<sub>3</sub>, Mg<sub>2</sub>SiO<sub>4</sub>, MgO, as well as a small amount of elemental nickel. The sum of these silicate phases covers 41% of the overall mass after the elemental iron, nickel, and chromium. This ratio of silicates is higher than the traditional estimations of 10%–15%. It is a remarkable difference; nevertheless, dynamic processes such as evaporation and collisional stripping very likely influenced Mercury's overall composition (Bitsch & Battistini 2020; Santos et al. 2017).

In this approach, Mercury most probably formed in a temperature zone of 1400–1440 K. An increase in the temperature does not result in much difference in iron fraction, because iron is not produced above 1440 K.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Distribution of dominant nickel forms in the solar nebula. Dark shade: Ni (FCC), mid-shade: Ni, light shade: Ni₃S₂.</figcaption></figure>

### 4.1.2. Venus

Venus is considered to be similar in chemical composition and internal structure to the Earth in the respects of mass, mean density, and solar distance (Stolper 1980; Treiman 2009). This assumption was partially confirmed thanks to the Venera and Vega missions (Treiman 2009).

Venusian lava is measured to contain an FeO content equal to the Earth's, which is 8%, indicating the presence of a similar metallic core (Taylor 1991). K/U ratios of both planets point out a similarity for their volatile to refractory elements (Taylor 1991). The similar abundance of potassium also suggests that Venus and the Earth have comparable abundances of Fe and Mn, which leads to comparable bulk mantles and core compositions, as well as sizes (Treiman 2009).

Venusian normative constituents and the sequence of minerals aligning with pressure are also the same as those of the Earth (Stolper 1980). Therefore, our results for the Earth can be applied also to Venus.

### 4.1.3. Earth

The leading elements with their fractions for the entire planet are iron (32.1%), oxygen (30.1%), silicon (15.1%), magnesium (13.9%), sulfur (2.9%), nickel (1.8%), calcium (1.5%), and aluminum (1.4%); the core is thought to consist of 88.8% iron, 5.8% nickel, and 4.5% sulfur (Morgan & Anders 1980).

While the other major elements such as iron, oxygen, and silicon are incorporated into complex molecules with varying fractions throughout the disk, nickel does not combine with other species in detectable fractions. Therefore, its abundance is taken as a proxy for determining temperature for when the proto-Earth appeared. There are 63 nickel-bearing compounds in the calculations, but it mostly remains in elemental form, apart from 1.5% as Ni₃S₂ at 150–300 K. As shown in Figure 3, nickel peaks at 1460 K with a fraction of 3.9%, together with its face-centered cubic (FCC) allotrope.

The compositions for the possibilities of the maximum nickel abundance in the disk, 3.9% at 1460 K, and its traditionally estimated abundance, 1.8% at 1318 K, can be compared depending on Figure 1(a) to decide which one suggests a more realistic bulk for the Earth. For both cases, metallic iron makes up one-third of the solid phase of the planet. The next most abundant material is magnesium silicates, MgSiO₃ and Mg₂SiO₄, where Ni is found in the desired abundance. For the maximum Ni contribution, silicates of aluminum, Al₂SiO₅, and calcium, Ca₂SiO₄, come forward.

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Earth's refractory composition calculated for the temperature where elemental nickel is found at the estimated abundance.</figcaption></figure>

Another important difference is that magnesium silicates comprise a fraction of only 0.2% in the case of maximum Ni. For the desired Ni fraction instead, metallic iron, nickel, chromium, and ZnCo₃ have mass fractions of 31.6%, 1.8%, 0.3%, and 0.1% respectively, which can safely account for the core, as demonstrated in Figure 4. 48.7% is distributed among silicates, with 21.6% of MgSiO₃, 18.8% of Mg₂SiO₄, 4.3% of CaSiO₃, 4.0% of Al₂SiO₅, and 0.2% of K₂TiO₃. The leading oxides are SiO₂ with 8.3% and MgO with 8.2%. The only carbon-bearing solid compound is MnC₂ with a fraction of 0.4%. Accordingly, the composition documented for the classical nickel fraction can account also for the distribution of other compounds in the bulk of the Earth.

### 4.1.4. Mars

Two substances are found in the literature with bulk abundance estimations of Mars. The first one is FeO, given in a range of 15%–19.5% by different methods (Taylor 2012). Although it is not close, the maximum FeO fraction found is 7.4% at 392 K throughout the solar disk, where the composition is dominated by Mg₄Si₆O₂₁H₁₂ with 45%, metallic iron with 14%, and MgO with 11%.

The second proxy is sulfur, which is claimed to have a core fraction of 16% (Rivoldini et al. 2011). In order to determine the disk region with the maximum sulfur content, combinations of sulfur with iron and nickel are included as well, since such combinations are likely under the physical conditions of the core.

The following compounds of Fe, Ni, and S included in the calculations: Fe, FeS, Fe₂S, Fe₂S₃, Fe₇S₈, Fe₉S₈, FeS₂, S, Ni, Ni₃S₄, Ni₆S₅, Ni₇S₆, Ni₉S₈, NiS, and NiS₂. The temperature where the sum of these species is at the maximum is 490 K with a bigger contribution from Fe–S molecules, and the related composition is shown in Figure 5. At that temperature, Mg₂SiO₄ is dominant with 33.5%, followed by metallic iron with 22.2%, then Mg₄Si₆O₂₁H₁₂ with 19.6%, and FeS with 8.6%. The rest consists of various silicates, metal oxides, and 1.5% of elemental nickel.

Since the leading magnesium silicate is Mg₂SiO₄ at 490 K while it is Mg₄Si₆O₁₂H₂₁ at 392 K, and since the abundance of elemental iron is bigger in the composition of 490 K, the option based on FeS content seems more probable.

Moriarty et al. (2014) mention 550 K as the temperature of the region where terrestrial planets form, which is only close to what we offer for the formation zone of Mars, but ~900 K lower than those of Mercury and the Earth.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Refractory composition of Mars calculated for the temperature where FeS is found at the maximum abundance.</figcaption></figure>

## 4.2. The System of WASP-12

Although WASP-12 is known to harbor at least one planet, in situ approaches used for solar planets cannot be utilized for its system. This is due to the fact that its detected planet is known to be a hot Jupiter, which is thought to have migrated to its current orbital distance after forming beyond snowlines according to the core accretion scenario. As the presence of a hot Jupiter may indicate a multiplanet system, the disk of WASP-12 is worth investigating to estimate the possible refractory compositions for any other planets that might be observed in the future. Indeed, both observations (Schneider 2020) and simulations (Raymond et al. 2005) confirm the existence and formation of terrestrial planets in extrasolar systems, despite the difficulty in revealing their chemical compositions in detail (Bond et al. 2010b).

To exemplify various planetary embryos in the disk of WASP-12, temperature zones with diverse compositions are chosen so that main compositional differences can be observed easily. The general view of refractory bulks is not radically different from that of solar ones, since the fundamental parameters such as mass, size, effective temperature, and elemental abundances are close to the solar values. The plots are restricted to the compounds that reach or exceed a mass fraction of 1%. The sample refractory compositions are considered down to 300 K, and the metallic, silicate, and oxide contributions differ in each selection. Hence, all representative cases fall into the category of rocky planets. It should be noted also that the classification according to the size, and more importantly the liquid and gas contributions, are not taken into account; therefore, the fractionations stand only for solid parts of the generated planets.

Six different refractory compositions are illustrated for the disk of WASP-12 in Figure 6. According to the results, the rocky part of a planet that formed around 1854 K should consist of 92% CrC<sub>24</sub>H<sub>36</sub> and 8% K<sub>2</sub>TiO<sub>3</sub>. This combination looks like an alternative planetary structure where CrC<sub>24</sub>H<sub>36</sub> comprises the mantle and nearly the entire planet, whereas K<sub>2</sub>TiO<sub>3</sub> covers the top as a thin crust. The stratification is not a mandatory expectation as both compounds can also be distributed arbitrarily.

For a planet formed at 1659 K, the entire rocky composition is shared by various oxides of calcium and/or aluminum, along with 28% K<sub>2</sub>TiO<sub>3</sub>. K<sub>2</sub>TiO<sub>3</sub> might precipitate as the planet's central core or make up a thin mantle. However, it seems more probable to find the compounds as mixed instead of stratified, since they all have the same chemical classification.

The composition at 1562 K is another example of the mixture of oxides with aluminum and/or calcium, including also some silicates. This combination is a candidate for terrestrial planets, probably with a large mantle and rich crust structure, but lacking a metallic core. Therefore, it can be categorized as a coreless planet (Seager & Elkins-Tanton 2008). K<sub>2</sub>TiO<sub>3</sub> is again the only oxide species bearing a metallic atom, with a ratio of 7%.

The model for 1026 K is a proper terrestrial planet with an elemental iron contribution of a mass fraction of 38% to make up the core, which may also lead to a high density (Bitsch & Battistini 2020). The rest of the planetary mass is divided among silicates of magnesium, aluminum, and calcium. This composition is compatible with the extrasolar planets of Kepler-10 and CoRoT-7b, which are known to resemble to Mercury by their iron cores making up to 65% of their bulks (Moriarty et al. 2014). Due to the presence of a metallic core, a magnetic field can also be expected if there is a motion to cause the dynamo effect.

The composition drawn for 441 K still has an iron core to cover 31% of the total mass; nevertheless, the major part consists of magnesium silicates. With a rich combination of FeS, FeO, Ca<sub>3</sub>Si<sub>2</sub>O<sub>7</sub>, and Cr<sub>2</sub>FeO<sub>4</sub>, this model is also placed in the category of terrestrial.

The iron core disappears for the composition found at 344 K. Iron is incorporated in oxide compounds such as Fe<sub>3</sub>O<sub>4</sub>, FeO, CaFe<sub>3</sub>O<sub>5</sub>, FeAl<sub>2</sub>O<sub>4</sub>, Cr<sub>2</sub>FeO<sub>4</sub>, or FeS. The dominant compound is Mg<sub>4</sub>Si<sub>6</sub>O<sub>21</sub>H<sub>12</sub> with 42%. This combination is reminiscent of the planet Mars, as far as what is known from analyses of its crust. A planet formed at that temperature could then be imagined as a coreless terrestrial planet resembling Mars.

Most of the sample refractory compositions from the disk of WASP-12 discussed above resemble the Earth, or at least remain within the definition of a terrestrial planet. This conclusion agrees with observations in optical and ultraviolet spectra of white dwarf stars, showing that extraterrestrial planetesimals are similar to the Earth's bulk composition since their masses are dominated by terrestrial planet-building elements such as oxygen, magnesium, silicon, and iron (Santos et al. 2017; Bond et al. 2010b).

As for the interior structure, all rocky planets can be differentiated into a silicate crust and a metallic core (Schaefer et al. 2012). The models presented above can also account for the core compositions of other planet types such as gas or ice giants.

# 5. Volatile Phase

The results for elemental enrichment, accreted volatile content, and accretion mass concerning the two disks are documented separately. The enrichments in nitrogen, sulfur, and phosphorus are not necessarily influenced by the changes in the C/O ratio when the conditions are reducing. Their reaction to the C/O ratio under oxidizing conditions is very small as well.

## 5.1. Jupiter

The oxygen fraction consumed by the refractory species is calculated to be 64% for the solar nebula, leaving 36% for volatiles.

The predictions for the composition of planetesimals accreted to the envelope of Jupiter are summarized in Table 6 (left half) corresponding to the four sets of the C/O ratio and redox state. The C/O ratio of the planetary envelope is calculated as 0.7 when the disk conditions are oxidizing, and 0.4 when reducing. Oxidizing circumstances lead to a C/O
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Diverse refractory compositions throughout the disk of WASP-12. Pie charts show the compositions for 1854 K, 1659 K, 1562 K, 1026 K, 441 K, and 344 K respectively, from left to right and top to bottom.</figcaption></figure>

ratio on Jupiter bigger than that of the Sun, and reducing circumstances to a C/O ratio lower than that of the Sun.

### 5.1.1. Enrichment of Basic Elements

Five fundamental elements, O, C, N, S, and P, are found in higher abundances in Jupiter's envelope than in the Sun. The highest enrichment is in oxygen, being 6 times, and observed when the conditions are reducing with the solar C/O ratio. Its lowest enrichment is 1.5 times under reducing conditions, with a C/O ratio of unity. The oxygen enrichment factor is bigger when C/O is 0.55 than when it is 1, since the disk becomes oxygen-deficient as C/O increases.

The highest carbon enrichment is 4.6 when the conditions are oxidizing and the C/O ratio is 1. This value is in agreement

with Mordasini et al. (2016), who mentioned it to be around 4 times. Its lowest level is 2.8 under the reducing conditions of the same C/O ratio.

### 5.1.2. Accretion Abundances

Water accretion is a small quantity of $0.8 M_\oplus$ under oxidizing conditions; in reducing conditions $CH_4$ leaves all available oxygen to $H_2O$ to comprise a mass of 12–16 $M_\oplus$. If C/O is 1, there is no water contribution for oxidizing conditions, and 4–5 $M_\oplus$ for reducing conditions.

Oxygen-bearing carbon species appear only under oxidizing conditions, comprising 6–7 $M_\oplus$ of $CO$, 3 $M_\oplus$ of $CO_2$, and 0.5 $M_\oplus$ of $CH_3OH$ for the solar C/O ratio. These values decrease slightly when the C/O ratio increases to 1.
<table><caption>**Table 6**
The Abundances of Accreted Volatile Species, and Enrichment Factors for Basic Atmospheric Elements for Generic Jupiter and WASP-12b for Minimum and Maximum Fits</caption><thead><tr><th rowspan="3"></th><th colspan="4">Sun</th><th colspan="4">WASP-12</th></tr><tr><th colspan="2">C/O: Solar</th><th colspan="2">C/O: 1</th><th colspan="2">C/O: Stellar</th><th colspan="2">C/O: 1</th></tr><tr><th>Oxidizing</th><th>Reducing</th><th>Oxidizing</th><th>Reducing</th><th>Oxidizing</th><th>Reducing</th><th>Oxidizing</th><th>Reducing</th></tr></thead><tbody><tr><td>C/O</td><td>0.74</td><td>0.39</td><td>1</td><td>1</td><td>1.1</td><td>0.50</td><td>1</td><td>1</td></tr><tr><td>Oxygen<sup>a</sup></td><td>2.7-3.2</td><td>4.4-6.0</td><td>2.0-2.5</td><td>1.5-2.1</td><td>1.5-1.9</td><td>2.7-3.6</td><td>1.6-2.0</td><td>1.3-1.7</td></tr><tr><td>Carbon<sup>a</sup></td><td>3.7-4.3</td><td>3.1-4.2</td><td>3.7-4.6</td><td>2.8-3.8</td><td>3.6-4.5</td><td>3.0-4.1</td><td>3.6-4.5</td><td>3.0-3.8</td></tr><tr><td>Nitrogen<sup>a</sup></td><td>2.9-3.3</td><td>4.5-6.2</td><td>2.8-3.5</td><td>4.5-6.2</td><td>2.8-3.5</td><td>4.6-6.2</td><td>2.8-3.5</td><td>4.9-6.2</td></tr><tr><td>Sulfur<sup>a</sup></td><td>3.3-3.8</td><td>2.4-3.3</td><td>3.1-3.8</td><td>2.4-3.3</td><td>3.0-3.8</td><td>2.4-3.3</td><td>3.0-3.8</td><td>2.6-3.3</td></tr><tr><td>Phosphorus<sup>b</sup></td><td>4.9-5.7</td><td>4.5-6.1</td><td>4.9-6.0</td><td>4.5-6.1</td><td>4.7-6.0</td><td>4.5-6.1</td><td>4.7-6.0</td><td>4.9-6.1</td></tr><tr><td>M(H<sub>2</sub>O) (M<sub>⊕</sub>)</td><td>0.8-0.9</td><td>11.7-15.9</td><td>0</td><td>4.1-5.5</td><td>0</td><td>10.5-14.2</td><td>0</td><td>5.3-6.6</td></tr><tr><td>M(CO) (M<sub>⊕</sub>)</td><td>6.1-7.1</td><td>0</td><td>5.0-6.2</td><td>0</td><td>5.5-7.0</td><td>0</td><td>5.9-7.4</td><td>0</td></tr><tr><td>M(CO<sub>2</sub>) (M<sub>⊕</sub>)</td><td>2.9-3.3</td><td>0</td><td>2.3-2.9</td><td>0</td><td>2.6-3.3</td><td>0</td><td>2.7-3.5</td><td>0</td></tr><tr><td>M(CH<sub>3</sub>OH) (M<sub>⊕</sub>)</td><td>0.5-0.6</td><td>0</td><td>0.4-0.5</td><td>0</td><td>0.4-0.6</td><td>0</td><td>0.5-0.6</td><td>0</td></tr><tr><td>M(O) (M<sub>⊕</sub>)</td><td>5.5-6.3</td><td>10.4-14.2</td><td>3.9-4.9</td><td>3.6-4.9</td><td>4.3-5.5</td><td>9.3-12.6</td><td>4.6-5.8</td><td>4.7-5.9</td></tr><tr><td>M(H<sub>2</sub>O)/M<sub>total</sub></td><td>0.1</td><td>0.6</td><td>0</td><td>0.4</td><td>0</td><td>0.6</td><td>0</td><td>0.4</td></tr><tr><td>M(CO)/M<sub>total</sub></td><td>0.5</td><td>0</td><td>0.5</td><td>0</td><td>0.5</td><td>0</td><td>0.5</td><td>0</td></tr><tr><td>M(CO<sub>2</sub>)/M<sub>total</sub></td><td>0.2</td><td>0</td><td>0.2</td><td>0</td><td>0.2</td><td>0</td><td>0.2</td><td>0</td></tr><tr><td>M(CH<sub>3</sub>OH)/M<sub>total</sub></td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>M<sub>Ice</sub> (M<sub>⊕</sub>)</td><td>11.6-13.5</td><td>17.7-24.1</td><td>9.9-12.4</td><td>9.7-13.1</td><td>11.6-14.7</td><td>17.9-24.4</td><td>11.9-15.0</td><td>13.0-16.4</td></tr><tr><td>M<sub>Rock</sub> (M<sub>⊕</sub>)</td><td>3.9-4.5</td><td>7.4-10.0</td><td>2.8-3.5</td><td>2.6-3.5</td><td>8.6-10.9</td><td>18.6-25.3</td><td>6.0-7.6</td><td>19.8-25.0</td></tr><tr><td>M<sub>Ice</sub>/M<sub>Rock</sub></td><td>3.0</td><td>2.4</td><td>3.6</td><td>3.8</td><td>1.3</td><td>1</td><td>2</td><td>0.7</td></tr></tbody></table>

**Notes.** Data acquired by the Galileo probe.

The total mass of ice accretion is naturally bigger with a lower C/O ratio (Pekmezci et al. 2019). The difference is more obvious under reducing conditions, where the ice mass is almost doubled as C/O falls to 0.55 from 1. Once the C/O ratio is 1, the redox state of the environment does not make an important difference to the ice accretion.

The dominant accretion under oxidizing circumstances is from CO and CO₂, and this is not affected if the C/O ratio is doubled. For the reducing cases, H₂O and CH₄ dominate. When C/O is doubled, H₂O loses two-thirds of its total fraction, while CH₄ nearly doubles.

Rock accretion is not affected by the increase in the C/O ratio under oxidizing conditions. However, it drops to one-third under reducing conditions, in agreement with the trend noted by Cridland et al. (2019b).

### 5.1.3. Accreted Planetesimal Composition

The mass fractions of the planetesimals accreted to Jupiter depending on the C/O ratio and redox state are depicted in Figure 7 in a comparable manner. In the case of solar C/O ratio and oxidizing conditions, the mass distribution among volatiles is 52% CO, 25% CO₂, 7% H₂O, 6% N₂, 4% CH₃OH, 4% H₂S, and others. For the case of solar C/O ratio and reducing

conditions, H₂O dominates with 65%, followed by CH₄ with 23%, NH₃ with 8%, and others. When the C/O ratio is 1 and the redox state is oxidizing, the distribution is given as follows: CO 50%, CO₂ 24%, CH₄ 9%, N₂ 7%, H₂S 4%, and others. Finally, when the C/O ratio is 1 and conditions are reducing, H₂O makes up 42%, CH₄ 38%, NH₃ 15%, and the others are less.

### 5.1.4. Comparison with Measurements

In the envelope of Jupiter, the dominant oxygen-carrying gas is H₂O, the major carbon-bearing gas is CH₄ followed by a trace amount of CO, the prime nitrogen-holder is NH₃, and the main sulfur gas is H₂S (Lodders & Fegley 2002; Fegley & Lodders 1994). Unlike in the solar disk, CH₄ overcomes CO, and NH₃ overcomes N₂ in its atmosphere. This is because the accreted CO and N₂ are processed into CH₄ and NH₃, respectively, due to the high temperature and pressure conditions provided by Jupiter's gravity (Prinn & Fegley 1981; Fegley & Lodders 1994; Bézard et al. 2002; Wong et al. 2008).

The in situ measurements collected from the Galileo probe revealed that the noble elements such as Ar, Kr, Xe, as well as C, N, and S are enriched in the atmosphere of Jupiter by about three times relative to their solar values (Madhusudhan et al.

<footnote>a Wong et al. (2004).</footnote>

<footnote>b Fletcher et al. (2009).</footnote>
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Distribution of planetesimal volatile composition in the accretion zone of Jupiter for oxidizing conditions (light shaded bars on the left) and reducing conditions (dark shaded bars on the right), and for C/O ratios of 0.55 and 1.</figcaption></figure>

2011a; Guillot 2008; Alibert et al. 2005b; Lodders & Fegley 1994, 2002; Hersant et al. 2001).

As noted by Mousis et al. (2012), phosphorus does not show an enrichment factor close to 3 in either of the sets considered. Varying in the range 4.5–6.1, phosphorus does not demonstrate dependence on either the C/O ratio or the redox state. Therefore, it is left out of the evaluation of the comparison between the measured and calculated values.

What is responsible for the supersolar abundances is thought to be a large quantity of planetesimal accretion into the planet from mixed clathrate hydrates, rocks, and ices that were composed of the solar disk material during its formation era, in accordance with the core accretion model (Hersant et al. 2001; Mousis et al. 2009c; Öberg et al. 2011). Nevertheless, an alternative enrichment model from the gas phase is given by Mousis et al. (2019).

The data from the Galileo probe for the 1 bar pressure zone in the atmosphere of Jupiter showed a low oxygen abundance, 30% of the solar value, contrary to expectations (Mousis et al. 2012). If low oxygen content in the Jovian atmosphere was not only a result of atmospheric dynamics and meteorological activities at the spot where the probe made the measurements, but also indicated a lack of oxygen abundance relative to carbon in the planetary bulk, then Jupiter may have formed under disk conditions where oxygen was depleted relative to carbon (Kuchner & Seager 2005; Mousis et al. 2012). Arguing against this are the results of the Juno microwave measurements, which suggest a water abundance at low northern latitudes in excess of solar (Li et al. 2020a), though this remains to be confirmed by additional data.

Since an oxygen enhancement factor bigger than 5 is ruled out by thermochemical models such as those designed by Visscher & Moses (2011) and Visscher et al. (2010), the case in which the C/O ratio is solar, the conditions are reducing, and

an enrichment by 4.4 to 6.0 times is required can safely be kept out of consideration. Among three other cases, the one requiring the biggest range of oxygen enrichment is where the C/O is solar and conditions are oxidizing, with a range of 2.7 to 3.2. The oxygen enrichment range for this case is calculated to be 6.8 to 7.2 by Mousis et al. (2012), nearly two times bigger than the value found here. This is a natural result since the initial oxygen abundance is smaller here, being the quantity remaining after its consumption by rocks.

The enrichment interval of oxygen calculated by Mousis et al. (2012) for C/O ~ 1 and the oxidizing case, 2.1–2.4 times the solar value, is close enough to our results, 2.0–2.5. Therefore, both studies suggest an average enrichment factor of 2.2 over the oxygen abundance in the solar nebula in order to maintain the compositional data measured by the Galileo probe in the atmosphere of Jupiter. This corresponds to the accretion of 3.9–4.9 $M_⊕$ of oxygen into the planet.

Unlike the study of Mousis et al. (2012), this study also documents reducing conditions. In the case where the C/O ratio is 1 and the conditions are reducing, the necessary oxygen enrichment factor is even lower than in the previous case, from 1.5 to 2.1. The oxygen mass required to be accreted into Jupiter is then 3.6–4.9 $M_⊕$. Therefore, consistently with the standard core accretion model, the minimum oxygen accretion required for explaining the enrichment of noble gases, carbon, nitrogen, and sulfur in the Jovian atmosphere by three times over the solar value is a C/O ratio of unity under reducing conditions, limiting oxygen accretion to only 4.3 $M_⊕$. However, if the enrichment factor of oxygen in the envelope of Jupiter is in the middle of the range of the values in Li et al. (2020a), 2 or 3 times solar, then oxidizing conditions are preferable for both C/O cases.

The bulk oxygen fraction and the C/O ratio of Jupiter remain to be better determined with microwave measurements by Juno at other latitudes. Should the water abundance turn out to be low, it could suggest that Jupiter is a carbon-rich planet (Mordasini et al. 2016). The challenge is to produce realistic models that obtain such a low oxygen abundance in the disk, so that little water is produced. A carbon-rich model for Jupiter would be consistent with the ideas of Lodders (2003) and would even relieve the problem from Juno gravity data that the amount of heavy elements in Jupiter’s envelope must be limited (Wahl et al. 2017).

## 5.2. WASP-12b

The model developed for Jupiter by Mousis et al. (2009c) is applied to the planet WASP-12 after the refractory formation and on the basis of the core accretion model. This method was utilized for Jupiter and Saturn before (e.g., Mousis et al. 2009c, 2012). It is then a good opportunity to test the validity of the estimation of volatile enrichment mass on extrasolar planets.

For the disk of WASP-12, the oxygen fraction consumed by the refractory species is calculated to be 39%, leaving 61% for volatiles.

The predictions for the composition of planetesimals accreted to the envelope of a generic hot Jupiter formed in the stellar disk of WASP-12, namely WASP-12b, are summarized in Table 6 (right half) corresponding to the four fits. The C/O ratio of the planetary envelope is calculated to be 1.1 when the disk conditions are oxidizing, and 0.5 when reducing. Oxidizing circumstances lead to a C/O ratio in WASP-12b bigger than in its parent star, and reducing
circumstances make the C/O ratio lower than that of its parent star.

### 5.2.1. Enrichment of Basic Elements

Five fundamental elements are found in higher abundances in the atmosphere of WASP-12b than in the parent star. The highest enrichment in oxygen is 3.6 when the calculations are done for the stellar C/O ratio under reducing conditions. The lowest one is 1.3, still for reducing conditions, but where C/O is 1. The range of oxygen enrichment for the oxidizing case shifts from 1.5–1.9 to 1.6–2.0 where C/O changes from stellar to 1.

### 5.2.2. Accretion Abundances

Accretion of carbon species for the stellar C/O ratio under oxidizing conditions is as follows: $5.5–7 M_⊕$ of CO, $2.5–3 M_⊕$ of CO$_2$, and $0.5 M_⊕$ of CH$_3$OH. These values slightly increase when the C/O ratio is 1.

Total ice contribution is found to be largest for the stellar C/O ratio under reducing conditions, with an average of $21 M_⊕$. It is not necessarily influenced by the C/O ratio for oxidizing circumstances in the disk of WASP-12.

Reducing conditions favor a larger contribution for the rocky part of a generic WASP-12b, which can increase up to $25.3 M_⊕$ for the stellar C/O ratio. The lowest rock accretion is found for a C/O ratio of unity but under oxidizing conditions, which can decrease down to $6 M_⊕$.

### 5.2.3. Composition of Accreted Planetesimals

The mass fractions of the accreted volatile species for the hot Jupiter WASP-12b depending on the C/O ratio and redox state are given in Figure 8 in a comparable manner. In the case of stellar C/O ratio and oxidizing conditions, the mass distribution among volatiles is as follows: CO 48%, CO$_2$ 22%, CH$_4$ 11%, N$_2$ 9%, CH$_3$OH 4%, H$_2$S 4%, and others, while there is no room left for H$_2$O.

For the case of stellar C/O ratio and reducing conditions, H$_2$O dominates with 58%, then comes CH$_4$ with 26%, NH$_3$ with 12%, and the others.

When the C/O ratio is 1 and the redox state is oxidizing, the distribution is given as follows: CO 49%, CO$_2$ 23%, N$_2$ 9%, CH$_4$ 8%, CH$_3$OH 4%, H$_2$S 4%, and others.

Finally, for unity C/O and reducing conditions, H$_2$O comprises 41%, CH$_4$ 36%, NH$_3$ 19%, and others.

CO and CO$_2$ are not affected if the C/O ratio is doubled under oxidizing circumstances. When C/O is doubled under reducing conditions, H$_2$O loses some one-third of its total fraction, while CH$_4$ increases nearly 1.5 times.

### 5.2.4. Comparison with Measurements

WASP-12b, a transiting hot Jupiter with an equilibrium temperature of 2600 K, has a carbon-rich dayside atmosphere detected by the thermal emission and photometric observations from Spitzer (Öberg et al. 2011; Madhusudhan et al. 2011b, 2011a; Madhusudhan 2012; Johnson et al. 2012). The C/O ratio of the atmosphere is greater than 1, whereas that of its parent star is 0.44 (Johnson et al. 2012; Madhusudhan et al. 2011a, 2011b). This difference points out an inconsistency between the planetary atmosphere and the stellar composition since the protoplanetary disk does not provide a high bulk C/O ratio, the star is not carbon-rich, and heavy-element enrichment could not be met only by

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Distribution of planetesimal volatile composition in the accretion zone of WASP-12b for oxidizing conditions (light shaded bars on the left) and reducing conditions (dark shaded bars on the right), and for C/O ratios of 0.45 and 1.</figcaption></figure>

accretion of icy planetesimals (Madhusudhan et al. 2011a; Johnson et al. 2012; Moses et al. 2013).

The C/O ratio in the atmosphere of WASP-12b is the result of the composition of a low abundance of H$_2$O and high abundances of CH$_4$ and CO relative to what could be expected by thermochemical equilibrium models with solar elemental values (Madhusudhan et al. 2011a, 2011b; Madhusudhan 2012; Moses et al. 2013). The composition of an object formed out of disk material at this equilibrium temperature should have principally been CO and H$_2$O, with a small contribution of CH$_4$ and CO$_2$ (Madhusudhan et al. 2011a).

According to the core accretion model, the planet is thought to have originated from an oxygen-depleted zone of the protostellar disk (Johnson et al. 2012), beyond the snowline, and to have accreted the gas on its way toward the star (Madhusudhan et al. 2011b; Öberg et al. 2011). When the C/O ratio is 1 in the atmosphere, the majority of the oxygen is tied up in CO, leaving less than 10% for H$_2$O formation, and CH$_4$ exceeds the solar contribution by 10 to 100 times (Madhusudhan et al. 2011b). The calculations of Madhusudhan et al. (2011b) with stellar abundances led to a C/O ratio of 0.27 in the planetary envelope. In order to produce 1 for the planetary C/O ratio, they had to decrease the stellar oxygen entry by 40%.

However, unlike the claim by Madhusudhan et al. (2011b), the results of this work show that it is not strictly necessary for the planet WASP-12b to have formed in an oxygen-deficient zone of the protoplanetary disk in order to maintain an atmospheric C/O ratio of 1. The calculations for the stellar C/O ratio and oxidizing case are actually compatible with the measurements and models in the literature suggesting the C/O ratio to be 1 for the planet, since this case results in the planetary C/O ratio being 1.06. It should be noted, however,
that the initial amount of oxygen for calculations of volatile
composition is smaller here than it is in the study of
Madhusudhan et al. (2011b).

As long as the conditions are oxidizing, the stellar elemental
abundances naturally yield a C/O ratio of 1.06 in the envelope of
the planet. The required mass of oxygen accreted is 4.3–5.4 M<sub>⊕</sub>.
The oxygen accretion for the oxidizing case where the C/O ratio
is 1 but not stellar is only slightly bigger, being 4.6–5.8 M<sub>⊕</sub>. In
actual fact, for this C/O ratio and reducing conditions, the
necessary oxygen accretion is also very close, being 4.7–5.9 M<sub>⊕</sub>.
Since a reducing environment is producing a planetary C/O ratio
of 0.5 with the stellar elemental abundances, this case can safely
be disregarded.

The ice accretion for the case of stellar C/O and oxidizing conditions is 0.8–9.7 M<sub>⊕</sub> according to Madhusudhan et al. (2011b) while it is calculated to be 11.6–14.7 M<sub>⊕</sub> in this study. For the case of unity C/O and oxidizing conditions, this accretion is found to be 0.3–15.9 M<sub>⊕</sub> by them, but it is calculated as 11.9–15.0 M<sub>⊕</sub> here.

Madhusudhan et al. (2011b) used the same methodology as
this work given by Mousis et al. (2009a, 2009b), but they
calculated the predefined disk composition for a low temper-
ature such as 30 K, whereas this study used a temperature of
300 K. Furthermore, they did not consider the oxygen sink into
the silicates, unlike the planetary elementary abundances
determined after the formation of rocky species here.

6. Discussion and Conclusion

Close-in giant planets are expected to have an atmospheric
spectrum with a combination of H<sub>2</sub>, CH<sub>4</sub>, H<sub>2</sub>O, CO, CO<sub>2</sub>, and
NH<sub>3</sub> (Lodders & Fegley 2002; Line et al. 2010; Öberg et al.
2011; Moses et al. 2013). If strong H<sub>2</sub>O bands are detected
from such a planet, in combination with weak CH<sub>4</sub> bands, this
is considered to be a sign of a C/O ratio smaller than 1
(Madhusudhan 2012).

In contrast, if CH₄ is detected with strong features, especially
accompanied by weak H₂O bands, then that is an indication of a
C/O ratio equal to or greater than 1 (Madhusudhan 2012). CO
bands, however, should always appear in the spectrum, regardless
of the C/O ratio (Madhusudhan 2012). The C/O ratios of a few
transiting exoplanets are known, but the James Webb Space
Telescope is expected to be able to more accurately determine
those of many more (Madhusudhan 2012).

In agreement with the anticipation of theoretical studies,
spectroscopic observations confirm the existence of extrasolar
close-in giant planets whose atmospheres have an abundance of
carbon exceeding that of oxygen (Madhusudhan et al. 2011b;
Mousis et al. 2012; Moses et al. 2013; Pekmezci et al. 2019).
A high C/O ratio leads to the formation of carbon-rich
and oxygen-poor planets with diverse interior compositions
(Lodders 2010; Öberg et al. 2011; Madhusudhan et al. 2011a).

For a C/O ratio lower than 1, H₂O and CO₂ are abundant, as
well as CO, which is the major constituent regardless of the
C/O ratio (Moses et al. 2013; Lodders 2010). For C/O higher
than 1, CH₄, HCN, and C₂H₂ are abundant (Moses et al. 2013),
indicating a formation history in either a water-poor or a
carbon-rich region (Lodders 2010; Pekmezci et al. 2019). The
species that carry neither oxygen nor carbon are not necessarily
affected by changes in the C/O ratio, as also supported by
Moses et al. (2013).

The dominant molecules as a function of the C/O ratio and
redox states calculated in this work are in agreement with the

expectations and observations mentioned by previous studies.
Encouraged by the results collected from every case studied
here, extrasolar giant planets can be expected to have different
C/O ratios from their parent stars. When the stellar elemental
compositions are used, the modeled hot Jupiters exhibit
atmospheric C/O ratios higher than their stars under oxidizing
conditions, and vice versa under reducing conditions. If the
model used here is to be proven worthy of consideration at least
roughly, then the C/O ratios of the atmospheres of hot Jupiters
can be generalized to be overstellar under oxidizing conditions,
and substellar under reducing ones, due to gas and ice accretion
into their envelopes during their formation and migration
epochs (Pekmezci et al. 2019).

There are various mechanisms suggested for the difference
in the C/O ratio of a gas giant and that of its protostellar disk.
For instance, accretion of solids rich in carbonic content but
poor in water ice, or accretion of gas rich in CO but poor in
H₂O, increases the carbon content in a planet (Lodders 2010;
Moses et al. 2013). Diffusion and turbulence, inward drift and
trapping of ices, heterogeneity throughout the disk depending
on place and time, differentiation of physical and chemical
characteristics, snowlines of various ices at different distances
from the parent star, and many other factors can cause planets
to be enriched in or deprived of carbon or oxygen at diverse
levels (Lodders 2010; Moses et al. 2013).

However, even if these factors play small roles in
determining the final atmospheric C/O ratio of a hot Jupiter,
divergence from the stellar C/O ratio might be a natural
consequence of the core accretion scenario. This might be
taken into consideration for the solar system as well, since bulk
oxygen contents of all of the solar giants still remain
undetermined (Wong et al. 2008). The results of this study
therefore may serve theoretical modeling and analyses of solar
and extrasolar terrestrial planets, especially for the determina-
tion of their chemical compositions and estimation of their
formation conditions. Future implications from direct measure-
ments, especially high-quality spectral analysis for host stars, as
well as from theoretical models and calculations will shed light
on the credibility of the assumptions made and the self-
consistent results reached in the study.

This work has been partly carried out thanks to the support of
the A\*MINDEX project (n° ANR-11-IDEX-0001-02) funded by the
"Investissements d'Avenir" French Government program, mana-
ged by the French National Research Agency (ANR).

The authors thank Dr. Torrence Johnson for important
advice and suggestions that enhanced the science reported here.
And Dr. Mohamad Ali-Dib was very helpful with technical
details.

G. S. Pekmezci is indebted to the Institut UTINAM and Observatoire des Sciences de l'Univers de Besançon for the accommodation and research facilities provided. She also expresses her gratitude to Prof. Roberto Buonanno for his valuable explanations and mentoring, and Prof. Atila Özgüç for his continuous support.

Appendix

The complete list of input species taken into account in
calculations are given in Table A1.
<table><caption>Table A1
The Species Taken into Account in Calculations</caption><thead><tr><th></th><th></th><th>1. Elemental and Uncategorized Species</th><th></th><th></th></tr></thead><tbody><tr><td>Al</td><td>CaZn</td><td>Fe<sub>3</sub>Si</td><td>MnSi</td><td>Sc</td></tr><tr><td>Al<sub>3</sub>Ni</td><td>CaZn<sub>2</sub></td><td>Fe<sub>5</sub>Si<sub>3</sub></td><td>Na</td><td>Sc<sub>5</sub>Si<sub>3</sub></td></tr><tr><td>Al<sub>3</sub>Ni<sub>2</sub></td><td>Co</td><td>FeSi</td><td>Na<sub>3</sub>AlH<sub>6</sub></td><td>Si</td></tr><tr><td>Al<sub>3</sub>Ti</td><td>Co<sub>2</sub>Al<sub>9</sub></td><td>FeSi<sub>2</sub></td><td>NaAlH<sub>4</sub></td><td>Si<sub>2</sub>H<sub>6</sub></td></tr><tr><td>Al<sub>5</sub>Co<sub>2</sub></td><td>Co<sub>2</sub>Si</td><td>FeTi</td><td>NaH</td><td>Sr</td></tr><tr><td>AlH<sub>3</sub></td><td>CoAl</td><td>K</td><td>Ni</td><td>SrH<sub>2</sub></td></tr><tr><td>AlH<sub>4</sub>K</td><td>CoAl<sub>3</sub></td><td>KAlH<sub>4</sub></td><td>Ni<sub>2</sub>H</td><td>Ti</td></tr><tr><td>AlNi<sub>3</sub></td><td>CoSi</td><td>KH</td><td>Ni<sub>2</sub>Si</td><td>Ti<sub>5</sub>Si<sub>3</sub></td></tr><tr><td>AlTi</td><td>CoSi<sub>2</sub></td><td></td><td>Ni<sub>3</sub>Ti</td><td>TiCr<sub>2</sub></td></tr><tr><td>C</td><td>Cr</td><td>Mg</td><td>Ni<sub>7</sub>Si<sub>13</sub></td><td>TiH</td></tr><tr><td>C<sub>60</sub></td><td>Cr<sub>3</sub>Si</td><td>Mg<sub>2</sub>Ni</td><td>NiAl</td><td>TiH<sub>2</sub></td></tr><tr><td>C<sub>70</sub></td><td>Cr<sub>5</sub>Si<sub>3</sub></td><td>Mg<sub>2</sub>Si</td><td>NiH</td><td>TiSi<sub>2</sub></td></tr><tr><td>Ca</td><td>CrH</td><td>Mg<sub>48</sub>Zn<sub>52</sub></td><td></td><td>V</td></tr><tr><td>Ca<sub>2</sub>Si</td><td>CrSi</td><td>MgCuAl</td><td>NiSi</td><td>V<sub>3</sub>Si</td></tr><tr><td>CaAl<sub>2</sub></td><td>CrSi<sub>2</sub></td><td>MgH<sub>2</sub></td><td>NiSi<sub>2</sub></td><td>V<sub>5</sub>Si<sub>3</sub></td></tr><tr><td>CaAl<sub>4</sub></td><td>Cu</td><td>MgNi<sub>2</sub></td><td>NiTi</td><td>VSi<sub>2</sub></td></tr><tr><td>CaH<sub>2</sub></td><td>Cu<sub>2</sub>Mg</td><td>MgZn<sub>2</sub></td><td>NiTi<sub>2</sub></td><td>Zn</td></tr><tr><td>CaMg<sub>2</sub></td><td>CuMg<sub>2</sub></td><td>Mn</td><td>P</td><td>Zn</td></tr><tr><td>CaSi</td><td>Fe</td><td>Mn<sub>3</sub>Si</td><td>P<sub>4</sub></td><td>ZnCo<sub>3</sub></td></tr><tr><td>CaSi<sub>2</sub></td><td>Fe<sub>2</sub>Ti</td><td>Mn<sub>5</sub>Si<sub>3</sub></td><td>S</td><td></td></tr><tr><td>(AlS)<sub>2</sub></td><td>Ca</td><td>FeO<sub>2</sub></td><td>N<sub>2</sub>O<sub>5</sub></td><td>S<sub>6</sub></td></tr><tr><td>Al<sub>2</sub></td><td>Ca<sub>2</sub></td><td>FeS</td><td>N<sub>3</sub></td><td>S<sub>7</sub></td></tr><tr><td>Al<sub>2</sub>C<sub>2</sub></td><td>CaCoO<sub>2</sub></td><td>H<sub>2</sub></td><td>Na<sub>2</sub>(CN)<sub>2</sub></td><td>S<sub>8</sub></td></tr><tr><td>Al<sub>2</sub>O</td><td>CaH</td><td>H<sub>2</sub>CN</td><td>Na<sub>2</sub></td><td>Sc<sub>2</sub></td></tr><tr><td>Al<sub>2</sub>O<sub>2</sub></td><td>CaO</td><td>H<sub>2</sub>CN<sub>2</sub></td><td>Na<sub>2</sub>O</td><td>Sc<sub>2</sub>O</td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>CaS</td><td>H<sub>2</sub>CNN</td><td>Na<sub>2</sub>O<sub>2</sub></td><td>Sc<sub>2</sub>O<sub>2</sub></td></tr><tr><td>Al<sub>2</sub>S</td><td>CCN</td><td>H<sub>2</sub>O</td><td>Na<sub>2</sub>O<sub>2</sub>H<sub>2</sub></td><td>ScC<sub>2</sub></td></tr><tr><td>AlC</td><td>CH</td><td>H<sub>2</sub>O<sub>2</sub></td><td>Na<sub>2</sub>SO<sub>4</sub></td><td>ScO</td></tr><tr><td>AlC<sub>2</sub></td><td>CH<sub>2</sub></td><td>H<sub>2</sub>S</td><td>NaCN</td><td>ScO<sub>2</sub></td></tr><tr><td>AlC<sub>3</sub>H<sub>9</sub></td><td>CH<sub>2</sub></td><td>H<sub>2</sub>S<sub>2</sub></td><td>NaCu</td><td></td></tr><tr><td>AlCuS</td><td>CH<sub>2</sub>N<sub>2</sub></td><td>H<sub>2</sub>SiO<sub>3</sub></td><td>NaH</td><td>ScS</td></tr><tr><td>AlCuS<sub>2</sub></td><td>CH<sub>2</sub>N<sub>4</sub></td><td>H<sub>2</sub>SO</td><td>NaNO<sub>2</sub></td><td>Si<sub>2</sub></td></tr><tr><td>AlH</td><td>CH<sub>3</sub></td><td>H<sub>2</sub>SO<sub>4</sub></td><td>NaNO<sub>3</sub></td><td>Si<sub>2</sub>C</td></tr><tr><td>AlH<sub>2</sub></td><td>CH<sub>3</sub>N</td><td>H<sub>3</sub>N<sub>2</sub></td><td>NaO</td><td>Si<sub>2</sub>H<sub>4</sub></td></tr><tr><td>AlH<sub>3</sub></td><td>CH<sub>3</sub>N<sub>2</sub></td><td>H<sub>3</sub>N<sub>2</sub></td><td>NaOH</td><td>Si<sub>2</sub>H<sub>6</sub></td></tr><tr><td>AlN</td><td>CH<sub>3</sub>N<sub>2</sub>CH<sub>3</sub></td><td>HCCN</td><td>NaOP</td><td>Si<sub>2</sub>N</td></tr><tr><td>AlO</td><td>CH<sub>3</sub>N<sub>3</sub></td><td>HCN</td><td>NaPO<sub>2</sub></td><td>Si<sub>2</sub>O<sub>2</sub></td></tr><tr><td>AlO<sub>2</sub></td><td>CH<sub>3</sub>NHNHCH<sub>3</sub></td><td>HCNH</td><td>NaPO<sub>3</sub></td><td></td></tr><tr><td>AlP</td><td>CH<sub>4</sub></td><td>HN<sub>3</sub></td><td>NCN</td><td>Si<sub>2</sub>P</td></tr><tr><td>AlP<sub>2</sub></td><td>CH<sub>4</sub>N</td><td>HNC</td><td>NH</td><td>Si<sub>2</sub>P<sub>2</sub></td></tr><tr><td>AlS</td><td>CH<sub>4</sub>N<sub>2</sub></td><td>HNO</td><td>NH<sub>2</sub></td><td>Si<sub>3</sub></td></tr><tr><td>AlS<sub>2</sub></td><td>CH<sub>5</sub>N</td><td>HNO<sub>2</sub></td><td>NH<sub>2</sub>NO<sub>2</sub></td><td>Si<sub>4</sub></td></tr><tr><td>C<sub>10</sub>N</td><td>CH<sub>5</sub>N<sub>2</sub></td><td>HNO<sub>3</sub></td><td>NH<sub>2</sub>OH</td><td>SiC</td></tr><tr><td>C<sub>10</sub>N<sub>2</sub></td><td>CH<sub>5</sub>N<sub>3</sub></td><td>HO<sub>2</sub></td><td>NH<sub>3</sub></td><td>SiC<sub>2</sub></td></tr><tr><td>C<sub>11</sub>N</td><td>CH<sub>6</sub>N<sub>2</sub></td><td>HO<sub>2</sub>S</td><td>Ni(C<sub>5</sub>H<sub>5</sub>)<sub>2</sub></td><td>SiH</td></tr><tr><td>C<sub>2</sub></td><td>CHN<sub>2</sub></td><td>HPO</td><td>Ni(OH)<sub>2</sub></td><td>SiH<sub>2</sub></td></tr><tr><td>C<sub>2</sub>H</td><td>CN</td><td>HS</td><td>Ni<sub>2</sub></td><td>SiH<sub>3</sub></td></tr><tr><td>C<sub>2</sub>H<sub>2</sub></td><td>CN<sub>2</sub></td><td>HS<sub>2</sub></td><td></td><td>SiH<sub>4</sub></td></tr><tr><td>C<sub>2</sub>H<sub>2</sub>N</td><td>CNN</td><td>HSO</td><td>NiH</td><td>SiN</td></tr><tr><td>C<sub>2</sub>H<sub>2</sub>N<sub>4</sub></td><td>CO</td><td>HSO<sub>3</sub></td><td>NiO</td><td>SiO</td></tr><tr><td>C<sub>2</sub>H<sub>3</sub></td><td>CO<sub>2</sub></td><td>HSOH</td><td>NiOH</td><td>SiO<sub>2</sub></td></tr><tr><td>C<sub>2</sub>H<sub>3</sub>N</td><td>Co<sub>2</sub></td><td>K<sub>2</sub>(CN)<sub>2</sub></td><td>NiS</td><td>SiOOH</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub></td><td>CoH</td><td>K<sub>2</sub></td><td>NO</td><td>SiP</td></tr><tr><td>C<sub>2</sub>H<sub>4</sub>N<sub>4</sub></td><td>CoO</td><td>K<sub>2</sub>(OH)<sub>2</sub></td><td>NO<sub>2</sub></td><td>SiP<sub>2</sub></td></tr><tr><td>C<sub>2</sub>H<sub>5</sub></td><td>CoPO<sub>2</sub></td><td>K<sub>2</sub>O</td><td>NO<sub>3</sub></td><td>SiS</td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>N</td><td>CoPO<sub>3</sub></td><td>K<sub>2</sub>O<sub>2</sub></td><td>O<sub>2</sub></td><td>SiS<sub>2</sub></td></tr><tr><td>C<sub>2</sub>H<sub>5</sub>N<sub>3</sub></td><td>COS</td><td>K<sub>2</sub>S</td><td>O<sub>3</sub></td><td>SN</td></tr><tr><td>C<sub>2</sub>H<sub>6</sub></td><td>CoS</td><td>K<sub>2</sub>SO<sub>4</sub></td><td>OH</td><td>SO</td></tr><tr><td>C<sub>2</sub>H<sub>6</sub>N</td><td>CP</td><td>KCN</td><td>P<sub>2</sub></td><td>SO<sub>2</sub></td></tr><tr><td>C<sub>2</sub>H<sub>7</sub>N</td><td>CP<sub>2</sub></td><td>KH</td><td>P<sub>2</sub>O<sub>3</sub></td><td>SO<sub>3</sub></td></tr><tr><td>C<sub>2</sub>H<sub>7</sub>N<sub>2</sub></td><td>Cr(CO)<sub>6</sub></td><td></td><td>P<sub>2</sub>O<sub>4</sub></td><td>SOH</td></tr><tr><td>C<sub>2</sub>H<sub>8</sub>N<sub>2</sub></td><td>Cr<sub>2</sub></td><td>KNa</td><td>P<sub>2</sub>O<sub>5</sub></td><td>Sr(OH)<sub>2</sub></td></tr><tr><td>C<sub>2</sub>N</td><td>Cr<sub>2</sub>O</td><td>KNO<sub>2</sub></td><td>P<sub>2</sub>S<sub>3</sub></td><td>Sr<sub>2</sub></td></tr><tr><td>C<sub>2</sub>N<sub>2</sub></td><td>Cr<sub>2</sub>O<sub>2</sub></td><td>KNO<sub>3</sub></td><td>P<sub>3</sub></td><td>SrCoO<sub>2</sub></td></tr></tbody></table>
<table><caption>Table A1
(Continued)</caption><thead><tr><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>C<sub>2</sub>O</td><td>Cr<sub>2</sub>O<sub>3</sub></td><td>KO</td><td>P<sub>3</sub>O<sub>6</sub></td><td>SrH</td></tr><tr><td>C<sub>2</sub>P</td><td>CrH</td><td>KOH</td><td>P<sub>4</sub></td><td>SrO</td></tr><tr><td>C<sub>2</sub>P<sub>2</sub></td><td>CrN</td><td>KS</td><td>P<sub>4</sub>O<sub>10</sub></td><td>SrOH</td></tr><tr><td>C<sub>2</sub>S<sub>2</sub></td><td>CrO</td><td>Mg(OH)<sub>2</sub></td><td>P<sub>4</sub>O<sub>6</sub></td><td>SrS</td></tr><tr><td>C<sub>3</sub></td><td>CrO<sub>2</sub></td><td>Mg<sub>2</sub></td><td>P<sub>4</sub>O<sub>7</sub></td><td>Ti<sub>2</sub></td></tr><tr><td>C<sub>3</sub>N</td><td>CrO<sub>3</sub></td><td>MgH</td><td>P<sub>4</sub>O<sub>8</sub></td><td>TiH</td></tr><tr><td>C<sub>3</sub>O<sub>2</sub></td><td>CrPO<sub>2</sub></td><td>MgN</td><td>P<sub>4</sub>O<sub>9</sub></td><td>TiN</td></tr><tr><td>C<sub>4</sub></td><td>CrPO<sub>3</sub></td><td>MgO</td><td>P<sub>4</sub>S<sub>3</sub></td><td>TiO</td></tr><tr><td>C<sub>4</sub>N</td><td>CrS</td><td>MgOH</td><td>P<sub>4</sub>S<sub>4</sub></td><td>TiO<sub>2</sub></td></tr><tr><td>C<sub>4</sub>N<sub>2</sub></td><td>CrS<sub>2</sub></td><td>MgS</td><td>P<sub>4</sub>S<sub>5</sub></td><td>TiS</td></tr><tr><td>C<sub>5</sub></td><td>CS</td><td>MnH</td><td>P<sub>4</sub>S<sub>7</sub></td><td>TiS<sub>2</sub></td></tr><tr><td>C<sub>5</sub>N</td><td>CS<sub>2</sub></td><td>MnO</td><td>P<sub>5</sub>S<sub>3</sub></td><td>V<sub>2</sub>O<sub>3</sub>(OH)<sub>4</sub></td></tr><tr><td>C<sub>5</sub>N<sub>4</sub></td><td>Cu<sub>2</sub></td><td>MnO<sub>2</sub></td><td>PH</td><td>V<sub>4</sub>O<sub>10</sub></td></tr><tr><td>C<sub>6</sub></td><td>Cu<sub>2</sub></td><td>MnOH</td><td>PH<sub>2</sub></td><td>VN</td></tr><tr><td>C<sub>6</sub>N</td><td>Cu<sub>2</sub>S</td><td>MnS</td><td>PH<sub>3</sub></td><td>VO</td></tr><tr><td>C<sub>6</sub>N</td><td>CuH</td><td>N<sub>2</sub></td><td>PN</td><td>VO<sub>2</sub></td></tr><tr><td>C<sub>6</sub>N<sub>2</sub></td><td>CuO</td><td>N<sub>2</sub>H</td><td>PO</td><td>VS</td></tr><tr><td>C<sub>7</sub></td><td>CuS</td><td>N<sub>2</sub>H<sub>2</sub></td><td>PO<sub>2</sub></td><td>Zn(OH)<sub>2</sub></td></tr><tr><td>C<sub>7</sub>O</td><td>Fe(C<sub>5</sub>H<sub>5</sub>)<sub>2</sub></td><td>N<sub>2</sub>H<sub>4</sub></td><td>PS</td><td>ZnH</td></tr><tr><td>C<sub>7</sub>N</td><td>Fe(CO)<sub>5</sub></td><td>N<sub>2</sub>H<sub>4</sub>.H<sub>2</sub>O</td><td>S<sub>2</sub></td><td>ZnO</td></tr><tr><td>C<sub>8</sub></td><td>Fe<sub>2</sub></td><td>N<sub>2</sub>O</td><td>S<sub>2</sub>O</td><td>ZnOH</td></tr><tr><td>C<sub>8</sub>N</td><td>FeC<sub>17</sub>H<sub>16</sub></td><td>N<sub>2</sub>O<sub>2</sub></td><td>S<sub>3</sub></td><td>ZnS</td></tr><tr><td>C<sub>8</sub>N<sub>2</sub></td><td>FeH</td><td>N<sub>2</sub>O<sub>3</sub></td><td>S<sub>4</sub></td><td></td></tr><tr><td>C<sub>9</sub>N</td><td>FeO</td><td>N<sub>2</sub>O<sub>4</sub></td><td>S<sub>5</sub></td><td></td></tr><tr><td></td><td></td><td>3. Carbides</td><td></td><td></td></tr><tr><td>Al<sub>4</sub>C<sub>3</sub></td><td>Co<sub>2</sub>C</td><td>CrC<sub>24</sub>H<sub>36</sub></td><td>Mn<sub>15</sub>C<sub>4</sub></td><td>Ni(C<sub>5</sub>H<sub>5</sub>)<sub>2</sub></td></tr><tr><td>Al<sub>4</sub>SiC<sub>4</sub></td><td>Cr(C<sub>5</sub>H<sub>5</sub>)<sub>2</sub></td><td>Fe(C<sub>5</sub>H<sub>5</sub>)<sub>2</sub></td><td>Mn<sub>23</sub>C<sub>6</sub></td><td>Ni<sub>3</sub>C</td></tr><tr><td>AlC<sub>3</sub>H<sub>9</sub></td><td>Cr<sub>23</sub>C<sub>6</sub></td><td>Fe<sub>2</sub>C</td><td>Mn<sub>3</sub>C</td><td>SiC</td></tr><tr><td>AlC<sub>8</sub>H<sub>19</sub></td><td>Cr<sub>3</sub>C<sub>2</sub></td><td>Fe<sub>3</sub>C</td><td>Mn<sub>5</sub>C<sub>2</sub></td><td>SrC<sub>2</sub></td></tr><tr><td>AlC<sub>9</sub>H<sub>21</sub></td><td>Cr<sub>4</sub>C</td><td>Mg<sub>2</sub>C<sub>3</sub></td><td>Mn<sub>7</sub>C<sub>3</sub></td><td>TiC</td></tr><tr><td>CaC<sub>2</sub></td><td>Cr<sub>7</sub>C<sub>3</sub></td><td>MgC<sub>2</sub></td><td>MnC<sub>2</sub></td><td>V<sub>2</sub>C</td></tr><tr><td>Co(C<sub>5</sub>H<sub>5</sub>)<sub>2</sub></td><td>CrC<sub>16</sub>H<sub>20</sub></td><td>Mn(C<sub>5</sub>H<sub>5</sub>)<sub>2</sub></td><td>Na<sub>2</sub>C<sub>2</sub></td><td>VC</td></tr><tr><td></td><td></td><td>4. Oxides</td><td></td><td></td></tr><tr><td>(NaPO<sub>3</sub>)<sub>3</sub></td><td>CoTiO<sub>3</sub></td><td>K<sub>3</sub>PO<sub>4.7</sub>H<sub>2</sub>O</td><td>Na<sub>2</sub>S<sub>2</sub>O<sub>3</sub></td><td>Sr<sub>2</sub>SiO<sub>4</sub></td></tr><tr><td>(NH<sub>4</sub>)<sub>2</sub>SO<sub>4</sub></td><td>CoZnTiO<sub>4</sub></td><td>K<sub>4</sub>P<sub>2</sub>O<sub>7</sub></td><td>Na<sub>2</sub>S<sub>2</sub>O<sub>3.5</sub>H<sub>2</sub>O</td><td>Sr<sub>2</sub>TiO<sub>4</sub></td></tr><tr><td>(NH<sub>4</sub>)<sub>2</sub>SO<sub>4.3</sub>NH<sub>3</sub></td><td>Cr(CO)<sub>6</sub></td><td>K<sub>4</sub>P<sub>2</sub>O<sub>7.3</sub>H<sub>2</sub>O</td><td>Na<sub>2</sub>S<sub>2</sub>O<sub>7</sub></td><td>Sr<sub>2</sub>V<sub>2</sub>O<sub>7</sub></td></tr><tr><td>12CaO.7Al<sub>2</sub>O<sub>3</sub></td><td>Cr<sub>2</sub>(SO<sub>4</sub>)<sub>3</sub></td><td>KH<sub>2</sub>PO<sub>4</sub></td><td>Na<sub>2</sub>SiO<sub>3</sub></td><td>Sr<sub>2</sub>VO<sub>4</sub></td></tr><tr><td>2Ca<sub>2</sub>SiO<sub>4</sub>.CaCO<sub>3</sub></td><td>Cr<sub>2</sub>FeO<sub>4</sub></td><td>KHSO<sub>4</sub></td><td>Na<sub>2</sub>SiO<sub>3.5</sub>H<sub>2</sub>O</td><td>Sr<sub>3</sub>(VO<sub>4</sub>)<sub>2</sub></td></tr><tr><td>2CaO.Al<sub>2</sub>O<sub>3</sub></td><td>Cr<sub>2</sub>O<sub>3</sub></td><td>KHSO<sub>5</sub></td><td>Na<sub>2</sub>SiO<sub>3.6</sub>H<sub>2</sub>O</td><td>Sr<sub>3</sub>Cr<sub>2</sub>O<sub>8</sub></td></tr><tr><td>2CaO.Al<sub>2</sub>O<sub>3</sub>.SiO<sub>2</sub></td><td>Cr<sub>3</sub>O<sub>4</sub></td><td>KMnO<sub>4</sub></td><td>Na<sub>2</sub>SiO<sub>3.8</sub>H<sub>2</sub>O</td><td>Sr<sub>3</sub>Fe<sub>2</sub>O<sub>6</sub></td></tr><tr><td>2CaO.Fe<sub>2</sub>O<sub>3</sub></td><td>Cr<sub>5</sub>O<sub>12</sub></td><td>KNaS<sub>2</sub>O<sub>7</sub></td><td>Na<sub>2</sub>SiO<sub>3.9</sub>H<sub>2</sub>O</td><td>Sr<sub>3</sub>SiO<sub>5</sub></td></tr><tr><td>2CaO.SiO<sub>2</sub></td><td>Cr<sub>8</sub>O<sub>21</sub></td><td>KNO<sub>2</sub></td><td>Na<sub>2</sub>SO<sub>3</sub></td><td>Sr<sub>3</sub>Ti<sub>2</sub>O<sub>7</sub></td></tr><tr><td>2CoO.SiO<sub>2</sub></td><td>CrO<sub>2</sub></td><td>KNO<sub>3</sub></td><td>Na<sub>2</sub>SO<sub>3.7</sub>H<sub>2</sub>O</td><td>Sr<sub>3</sub>V<sub>2</sub>O<sub>8</sub></td></tr><tr><td>2CoO.TiO<sub>2</sub></td><td>CrO<sub>3</sub></td><td>KO<sub>2</sub></td><td>Na<sub>2</sub>SO<sub>4</sub></td><td>Sr<sub>4</sub>Ti<sub>3</sub>O<sub>10</sub></td></tr><tr><td>2FeO.SiO<sub>2</sub></td><td>CrVO<sub>4</sub></td><td>KO<sub>3</sub></td><td>Na<sub>2</sub>SO<sub>4.10</sub>H<sub>2</sub>O</td><td>Sr<sub>7</sub>Fe<sub>10</sub>O<sub>22</sub></td></tr><tr><td>2MnO.TiO<sub>2</sub></td><td>Cu(VO<sub>3</sub>)<sub>2</sub></td><td>KOH</td><td>Na<sub>2</sub>SO<sub>4.7</sub>H<sub>2</sub>O</td><td>SrAl<sub>2</sub>Si<sub>2</sub>O<sub>8</sub></td></tr><tr><td>2Na<sub>2</sub>O.SiO<sub>2</sub></td><td>Cu<sub>2</sub>O</td><td>KOH.2H<sub>2</sub>O</td><td>Na<sub>2</sub>Ti<sub>6</sub>O<sub>13</sub></td><td>SrCO<sub>3</sub></td></tr><tr><td>2NiO.SiO<sub>2</sub></td><td>Cu<sub>2</sub>O.Al<sub>2</sub>O<sub>3</sub></td><td>KOH.H<sub>2</sub>O</td><td>Na<sub>2</sub>V<sub>2</sub>O<sub>6</sub></td><td>SrCoO<sub>3</sub></td></tr><tr><td>3Al<sub>2</sub>O<sub>3.2</sub>SiO<sub>2</sub></td><td>Cu<sub>2</sub>O.Fe<sub>2</sub>O<sub>3</sub></td><td>KPO<sub>3</sub></td><td>Na<sub>2</sub>V<sub>2</sub>O<sub>7</sub></td><td>SrCrO<sub>4</sub></td></tr><tr><td>3CaO.2SiO<sub>2</sub></td><td>Cu<sub>2</sub>P<sub>2</sub>O<sub>7</sub></td><td>KVO<sub>4</sub></td><td>Na<sub>2</sub>ZnO<sub>2</sub></td><td>SrFe<sub>12</sub>O<sub>19</sub></td></tr><tr><td>3CaO.Al<sub>2</sub>O<sub>3</sub></td><td>Cu<sub>2</sub>SO<sub>4</sub></td><td>Mg(NO<sub>3</sub>)<sub>2</sub></td><td>Na<sub>3</sub>HP<sub>2</sub>O<sub>7</sub></td><td>SrHPO<sub>4</sub></td></tr><tr><td>3CaO.Al<sub>2</sub>O<sub>3.3</sub>SiO<sub>2</sub></td><td>Cu<sub>3</sub>(PO<sub>4</sub>)<sub>2</sub></td><td>Mg(NO<sub>3</sub>)<sub>2.2</sub>H<sub>2</sub>O</td><td>Na<sub>3</sub>PO<sub>4</sub></td><td>SrO</td></tr><tr><td>3CaO.SiO<sub>2</sub></td><td>CuCO<sub>3</sub></td><td>Mg(NO<sub>3</sub>)<sub>2.6</sub>H<sub>2</sub>O</td><td>Na<sub>3</sub>PO<sub>4.12</sub>H<sub>2</sub>O</td><td>SrO.Al<sub>2</sub>O<sub>3</sub></td></tr><tr><td>3CaO.V<sub>2</sub>O<sub>5</sub></td><td>CuCrO<sub>2</sub></td><td>Mg(OH)<sub>2</sub></td><td>Na<sub>3</sub>VO<sub>4</sub></td><td>SrO<sub>2</sub></td></tr><tr><td>3MgO.2SiO<sub>2.2</sub>H<sub>2</sub>O</td><td>CuO</td><td>Mg(VO<sub>3</sub>)<sub>2</sub></td><td>Na<sub>4</sub>P<sub>2</sub>O<sub>7</sub></td><td>SrSiO<sub>3</sub></td></tr><tr><td>3MgO.4SiO<sub>2</sub>.H<sub>2</sub>O</td><td>CuO.Al<sub>2</sub>O<sub>3</sub></td><td>Mg<sub>2</sub>P<sub>2</sub>O<sub>7</sub></td><td>Na<sub>4</sub>P<sub>2</sub>O<sub>7.10</sub>H<sub>2</sub>O</td><td>SrSO<sub>4</sub></td></tr><tr><td>3Na<sub>2</sub>O.2SiO<sub>2</sub></td><td>CuO.Cr<sub>2</sub>O<sub>3</sub></td><td>Mg<sub>2</sub>SiO<sub>4</sub></td><td>Na<sub>4</sub>V<sub>2</sub>O<sub>7</sub></td><td>SrTi<sub>12</sub>O<sub>19</sub></td></tr><tr><td>3SrO.Al<sub>2</sub>O<sub>3</sub></td><td>CuO.CuSO<sub>4</sub></td><td>Mg<sub>2</sub>TiO<sub>4</sub></td><td>Na<sub>5</sub>P<sub>3</sub>O<sub>10</sub></td><td>SrTiO<sub>3</sub></td></tr><tr><td>3SrO.MgO.2SiO<sub>2</sub></td><td>CuO.Fe<sub>2</sub>O<sub>3</sub></td><td>Mg<sub>2</sub>V<sub>2</sub>O<sub>7</sub></td><td>Na<sub>5</sub>P<sub>3</sub>O<sub>10.6</sub>H<sub>2</sub>O</td><td>SrV<sub>2</sub>O<sub>6</sub></td></tr><tr><td>4CaO.3TiO<sub>2</sub></td><td>CuSO<sub>4</sub></td><td>Mg<sub>3</sub>(PO<sub>4</sub>)<sub>2</sub></td><td>Na<sub>6</sub>P<sub>2</sub>O<sub>8</sub></td><td>SrVO<sub>3</sub></td></tr><tr><td>4CaO.Al<sub>2</sub>O<sub>3</sub>.Fe<sub>2</sub>O<sub>3</sub></td><td>Fe(CO)<sub>5</sub></td><td>Mg<sub>3</sub>(VO<sub>4</sub>)<sub>2</sub></td><td>NaH<sub>2</sub>PO<sub>4</sub></td><td>Ti<sub>10</sub>O<sub>19</sub></td></tr><tr><td>4Na<sub>2</sub>O.5TiO<sub>2</sub></td><td>Fe(VO<sub>3</sub>)<sub>2</sub></td><td>Mg<sub>3</sub>Si<sub>2</sub>O<sub>5</sub>(OH)<sub>4</sub></td><td>NaH<sub>2</sub>SiO<sub>4.7</sub>H<sub>2</sub>O</td><td>Ti<sub>20</sub>O<sub>39</sub></td></tr><tr><td>4SrO.Al<sub>2</sub>O<sub>3</sub></td><td>FeO</td><td>Mg<sub>4</sub>Si<sub>6</sub>O<sub>21</sub>H<sub>12</sub></td><td>NaH<sub>2</sub>SiO<sub>4.8</sub>H<sub>2</sub>O</td><td>Ti<sub>2</sub>O<sub>3</sub></td></tr><tr><td>7MgO.8SiO<sub>2</sub>.H<sub>2</sub>O</td><td>Fe<sub>2</sub>(CO)<sub>9</sub></td><td>Mg<sub>4</sub>Si<sub>6</sub>O<sub>23</sub>H<sub>14</sub></td><td>NaH<sub>3</sub>P<sub>2</sub>O<sub>7</sub></td><td>Ti<sub>3</sub>O<sub>2</sub></td></tr><tr><td>Al<sub>2</sub>(SO<sub>4</sub>)<sub>3</sub></td><td>Fe<sub>2</sub>(SO<sub>4</sub>)<sub>3</sub></td><td>Mg<sub>6</sub>Si<sub>4</sub>O<sub>10</sub>(OH)<sub>8</sub></td><td>NaHSO<sub>4</sub></td><td>Ti<sub>3</sub>O<sub>5</sub></td></tr><tr><td>Al<sub>2</sub>CO</td><td>Fe<sub>2</sub>Al<sub>4</sub>Si<sub>5</sub>O<sub>18</sub></td><td></td><td>NaHSO<sub>4</sub>.H<sub>2</sub>O</td><td>Ti<sub>4</sub>O<sub>7</sub></td></tr></tbody></table>
<table><caption>Table A1
(Continued)</caption><thead><tr><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>Fe<sub>2</sub>O<sub>3</sub></td><td>Mg<sub>7</sub>Si<sub>8</sub>O<sub>22</sub>(OH)<sub>2</sub></td><td>NaMnO<sub>4</sub></td><td>Ti<sub>5</sub>O<sub>9</sub></td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub>.2SiO<sub>2</sub></td><td>Fe<sub>2</sub>SiO<sub>4</sub></td><td>MgMn<sub>2</sub>O<sub>4</sub></td><td>NaNH<sub>4</sub>HPO<sub>4</sub>.4H<sub>2</sub>O</td><td>Ti<sub>6</sub>O<sub>11</sub></td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub>.TiO<sub>2</sub></td><td>Fe<sub>2</sub>TiO<sub>4</sub></td><td>MgO</td><td>NaNO<sub>2</sub></td><td>Ti<sub>7</sub>O<sub>13</sub></td></tr><tr><td>Al<sub>2</sub>SiO<sub>5</sub></td><td>Fe<sub>2</sub>TiO<sub>5</sub></td><td>MgO<sub>2</sub></td><td>NaNO<sub>3</sub></td><td>Ti<sub>8</sub>O<sub>15</sub></td></tr><tr><td>Al<sub>4</sub>CO<sub>4</sub></td><td>Fe<sub>3</sub>(CO)<sub>12</sub></td><td>MgSiO<sub>3</sub></td><td>NaO<sub>2</sub></td><td>Ti<sub>9</sub>O<sub>17</sub></td></tr><tr><td>AlO</td><td>Fe<sub>3</sub>Al<sub>2</sub>Si<sub>3</sub>O<sub>12</sub></td><td>MgSO<sub>3</sub></td><td>NaO<sub>3</sub></td><td>TiNO</td></tr><tr><td>AlPO<sub>4</sub></td><td>Fe<sub>3</sub>O<sub>4</sub></td><td>MgSO<sub>3</sub>.3H<sub>2</sub>O</td><td>NaOH</td><td>TiO</td></tr><tr><td>Ca(PO<sub>3</sub>)<sub>2</sub></td><td>FeAl<sub>2</sub>O<sub>4</sub></td><td>MgSO<sub>3</sub>.6H<sub>2</sub>O</td><td>NaOH.H<sub>2</sub>O</td><td>TiO<sub>2</sub></td></tr><tr><td>Ca<sub>2</sub>P<sub>2</sub>O<sub>7</sub></td><td>FeCO<sub>3</sub></td><td>MgSO<sub>4</sub></td><td>NaPO<sub>3</sub></td><td>V<sub>2</sub>O<sub>3</sub></td></tr><tr><td>Ca<sub>2</sub>SiO<sub>4</sub></td><td>FeCr(VO<sub>4</sub>)<sub>2</sub></td><td>MgSO<sub>4</sub>.2H<sub>2</sub>O</td><td>NaV<sub>2</sub>O<sub>5</sub></td><td>V<sub>2</sub>O<sub>4</sub></td></tr><tr><td>Ca<sub>2</sub>V<sub>2</sub>O<sub>7</sub></td><td>FeO</td><td>MgSO<sub>4</sub>.4H<sub>2</sub>O</td><td>NaVO<sub>3</sub></td><td>V<sub>2</sub>O<sub>5</sub></td></tr><tr><td>Ca<sub>3</sub>(Al<sub>2</sub>Si<sub>2</sub>O<sub>8</sub>)<sub>3</sub>.CaCO<sub>3</sub></td><td>FeO</td><td>MgSO<sub>4</sub>.5H<sub>2</sub>O</td><td>NH<sub>2</sub>OH</td><td>V<sub>2</sub>O<sub>5</sub>.H<sub>2</sub>O</td></tr><tr><td>Ca<sub>3</sub>(PO<sub>4</sub>)<sub>2</sub></td><td>FeO.SiO<sub>2</sub></td><td>MgSO<sub>4</sub>.6H<sub>2</sub>O</td><td>NaH<sub>4</sub>.H<sub>2</sub>PO<sub>4</sub></td><td>V<sub>3</sub>O<sub>5</sub></td></tr><tr><td>Ca<sub>3</sub>Cr<sub>2</sub>(SiO<sub>4</sub>)<sub>3</sub></td><td>FeO.TiO<sub>2</sub></td><td>MgSO<sub>4</sub>.7H<sub>2</sub>O</td><td>NaH<sub>4</sub>HSO<sub>4</sub></td><td>V<sub>4</sub>O<sub>7</sub></td></tr><tr><td>Ca<sub>3</sub>Fe<sub>2</sub>Si<sub>3</sub>O<sub>12</sub></td><td>FePO<sub>4</sub></td><td>MgSO<sub>4</sub>.H<sub>2</sub>O</td><td>NaH<sub>4</sub>NO<sub>2</sub></td><td>V<sub>3</sub>O<sub>9</sub></td></tr><tr><td>Ca<sub>3</sub>Si<sub>2</sub>O<sub>7</sub></td><td>FeSiO<sub>3</sub></td><td>MgTi<sub>2</sub>O<sub>5</sub></td><td>NaH<sub>4</sub>NO<sub>3</sub></td><td>V<sub>6</sub>O<sub>11</sub></td></tr><tr><td>Ca<sub>3</sub>Si<sub>2</sub>O<sub>7</sub>.2CaCO<sub>3</sub></td><td>FeSO<sub>4</sub></td><td>MgTiO<sub>3</sub></td><td>NaH<sub>4</sub>OH</td><td>V<sub>6</sub>O<sub>13</sub></td></tr><tr><td>Ca<sub>3</sub>SiO<sub>5</sub></td><td>FeTi<sub>2</sub>O<sub>5</sub></td><td>MgV<sub>2</sub>O<sub>6</sub></td><td>NaH<sub>4</sub>VO<sub>3</sub></td><td>V<sub>7</sub>O<sub>13</sub></td></tr><tr><td>Ca<sub>3</sub>Ti<sub>2</sub>O<sub>7</sub></td><td>FeTiO<sub>3</sub></td><td>Mn(NO<sub>3</sub>)<sub>2</sub></td><td>Ni(NO<sub>3</sub>)<sub>2</sub></td><td>V<sub>8</sub>O<sub>15</sub></td></tr><tr><td>Ca<sub>3</sub>V<sub>2</sub>O<sub>8</sub></td><td>FeV<sub>2</sub>O<sub>4</sub></td><td>Mn(OH)<sub>2</sub></td><td>Ni(NO<sub>3</sub>)<sub>2</sub>.6H<sub>2</sub>O</td><td>VO</td></tr><tr><td>CaAl<sub>2</sub>SiO<sub>6</sub></td><td>FeVO<sub>4</sub></td><td>Mn(VO<sub>3</sub>)<sub>2</sub></td><td>Ni(OH)<sub>2</sub></td><td>VO<sub>2</sub></td></tr><tr><td>CaC<sub>2</sub>O<sub>4</sub></td><td>H<sub>2</sub>O</td><td>Mn<sub>2</sub>(SO<sub>4</sub>)<sub>3</sub></td><td>Ni(OH)<sub>3</sub></td><td>VOCO<sub>3</sub></td></tr><tr><td>CaCO<sub>3</sub></td><td>H<sub>2</sub>Si<sub>2</sub>O<sub>5</sub></td><td>Mn<sub>2</sub>O<sub>3</sub></td><td>Na<sub>2</sub>O<sub>3</sub>.3H<sub>2</sub>O</td><td>VOSO<sub>4</sub></td></tr><tr><td>CaCrO<sub>4</sub></td><td>H<sub>2</sub>SiO<sub>3</sub></td><td>Mn<sub>2</sub>O<sub>3</sub>.3H<sub>2</sub>O</td><td>Na<sub>2</sub>SiO<sub>4</sub></td><td>VOSO<sub>4</sub>.4H<sub>2</sub>O</td></tr><tr><td>CaFe(SiO<sub>3</sub>)<sub>2</sub></td><td>H<sub>2</sub>SO<sub>4</sub></td><td>Mn<sub>2</sub>SiO<sub>4</sub></td><td>Na<sub>2</sub>TiO<sub>4</sub></td><td>VOSO<sub>4</sub>.5H<sub>2</sub>O</td></tr><tr><td>CaFe<sub>3</sub>O<sub>5</sub></td><td>H<sub>2</sub>SO<sub>4</sub>.2H<sub>2</sub>O</td><td>Mn<sub>3</sub>O<sub>4</sub></td><td>Na<sub>3</sub>Mn<sub>2</sub>O<sub>4</sub></td><td>VOSO<sub>4</sub>.6H<sub>2</sub>O</td></tr><tr><td>CaFe<sub>5</sub>O<sub>7</sub></td><td>H<sub>2</sub>SO<sub>4</sub>.3H<sub>2</sub>O</td><td>MnO</td><td>NiO</td><td>Zn(NO<sub>3</sub>)<sub>2</sub></td></tr><tr><td>CaFeSiO<sub>3</sub></td><td>H<sub>2</sub>SO<sub>4</sub>.4H<sub>2</sub>O</td><td>MnO.OH</td><td>NiO.OH</td><td>Zn(NO<sub>3</sub>)<sub>2</sub>.2H<sub>2</sub>O</td></tr><tr><td>CaFeSiO<sub>4</sub></td><td>H<sub>2</sub>SO<sub>4</sub>.6.5H<sub>2</sub>O</td><td>MnO.TiO<sub>2</sub></td><td>NiSO<sub>4</sub></td><td>Zn(NO<sub>3</sub>)<sub>2</sub>.4H<sub>2</sub>O</td></tr><tr><td>CaO</td><td>H<sub>2</sub>SO<sub>4</sub>.H<sub>2</sub>O</td><td>MnO<sub>2</sub></td><td>NiSO<sub>4</sub>.4H<sub>2</sub>O</td><td>Zn(NO<sub>3</sub>)<sub>2</sub>.6H<sub>2</sub>O</td></tr><tr><td>CaO.2Al<sub>2</sub>O<sub>3</sub></td><td>H<sub>3</sub>PO<sub>4</sub></td><td>MnSiO<sub>3</sub></td><td>NiSO<sub>4</sub>.6H<sub>2</sub>O</td><td>Zn(OH)(NO<sub>3</sub>).H<sub>2</sub>O</td></tr><tr><td>CaO.6Al<sub>2</sub>O<sub>3</sub></td><td>H<sub>3</sub>PO<sub>4</sub>.0.5H<sub>2</sub>O</td><td>MnSO<sub>4</sub></td><td>NiSO<sub>4</sub>.7H<sub>2</sub>O</td><td>Zn(OH)<sub>2</sub></td></tr><tr><td>CaO.Al<sub>2</sub>O<sub>3</sub></td><td>HNO<sub>3</sub></td><td>MnSO<sub>4</sub>.4H<sub>2</sub>O</td><td>NiSO<sub>4</sub>.H<sub>2</sub>O</td><td>Zn<sub>2</sub>SiO<sub>4</sub></td></tr><tr><td>CaO.Al<sub>2</sub>O<sub>3</sub>.2SiO<sub>2</sub></td><td>KV<sub>2</sub>O<sub>5</sub></td><td>MnSO<sub>4</sub>.5H<sub>2</sub>O</td><td>NiTi<sub>2</sub>O<sub>5</sub></td><td>Zn<sub>2</sub>TiO<sub>4</sub></td></tr><tr><td>CaO.Al<sub>2</sub>O<sub>3</sub>.SiO<sub>2</sub></td><td>K<sub>2</sub>H<sub>2</sub>P<sub>2</sub>O<sub>7</sub></td><td>MnSO<sub>4</sub>.7H<sub>2</sub>O</td><td>NiTiO<sub>3</sub></td><td>Zn<sub>3</sub>(OH)<sub>4</sub>(NO<sub>3</sub>)<sub>2</sub></td></tr><tr><td>CaO.Cr<sub>2</sub>O<sub>3</sub></td><td>K<sub>2</sub>HPO<sub>4</sub></td><td>MnSO<sub>4</sub>.H<sub>2</sub>O</td><td>NiZnTiO<sub>4</sub></td><td>Zn<sub>3</sub>(PO<sub>4</sub>)<sub>2</sub></td></tr><tr><td>CaO.Fe<sub>2</sub>O<sub>3</sub></td><td>K<sub>2</sub>Mg(SO<sub>4</sub>)<sub>2</sub>.4H<sub>2</sub>O</td><td>MnTi<sub>2</sub>O<sub>5</sub></td><td>P<sub>2</sub>O<sub>5</sub></td><td>Zn<sub>3</sub>(VO<sub>4</sub>)<sub>2</sub></td></tr><tr><td>CaO.TiO<sub>2</sub></td><td>K<sub>2</sub>Mg(SO<sub>4</sub>)<sub>2</sub>.6H<sub>2</sub>O</td><td>MnTiO<sub>2</sub></td><td>P<sub>4</sub>O<sub>10</sub></td><td>Zn<sub>5</sub>(OH)<sub>8</sub>(NO<sub>3</sub>)<sub>2</sub>.2H<sub>2</sub>O</td></tr><tr><td>CaO.V<sub>2</sub>O<sub>5</sub></td><td>K<sub>2</sub>O</td><td>MnTiO<sub>3</sub></td><td>P<sub>4</sub>O<sub>8</sub></td><td>ZnCO<sub>3</sub></td></tr><tr><td>CaO<sub>2</sub></td><td>K<sub>2</sub>O.2SiO<sub>2</sub></td><td>N<sub>2</sub>O<sub>4</sub></td><td>PO<sub>2</sub></td><td>ZnCrO<sub>4</sub></td></tr><tr><td>CaSiO<sub>3</sub></td><td>K<sub>2</sub>O.4SiO<sub>2</sub></td><td>N<sub>2</sub>O<sub>5</sub></td><td>Sc(OH)<sub>3</sub></td><td>ZnFe<sub>2</sub>O<sub>4</sub></td></tr><tr><td>CaSO<sub>3</sub></td><td>K<sub>2</sub>O.SiO<sub>2</sub></td><td>Na<sub>2</sub>H<sub>2</sub>P<sub>2</sub>O<sub>7</sub></td><td>Sc(PO<sub>3</sub>)<sub>3</sub></td><td>ZnFeO<sub>4</sub></td></tr><tr><td>CaSO<sub>4</sub></td><td>K<sub>2</sub>O<sub>2</sub></td><td>Na<sub>2</sub>HPO<sub>4</sub></td><td>Sc<sub>2</sub>O<sub>3</sub></td><td>ZnMn<sub>2</sub>O<sub>4</sub></td></tr><tr><td>CaTiSiO<sub>5</sub></td><td>K<sub>2</sub>O<sub>3</sub></td><td>Na<sub>2</sub>HPO<sub>4</sub>.12H<sub>2</sub>O</td><td>ScPO<sub>4</sub></td><td>ZnO</td></tr><tr><td>CaV<sub>2</sub>O<sub>6</sub></td><td>K<sub>2</sub>O<sub>4</sub></td><td>Na<sub>2</sub>HPO<sub>4</sub>.2H<sub>2</sub>O</td><td>ScPO<sub>4</sub>.2H<sub>2</sub>O</td><td>ZnO.2ZnSO<sub>4</sub></td></tr><tr><td>Co<sub>2</sub>(CO)<sub>8</sub></td><td>K<sub>2</sub>S<sub>2</sub>O<sub>7</sub></td><td>Na<sub>2</sub>HPO<sub>4</sub>.7H<sub>2</sub>O</td><td>Si<sub>2</sub>N<sub>2</sub>O</td><td>ZnO.Al<sub>2</sub>O<sub>3</sub></td></tr><tr><td>Co<sub>2</sub>TiO<sub>4</sub></td><td>K<sub>2</sub>S<sub>2</sub>O<sub>8</sub></td><td>Na<sub>2</sub>O</td><td>SiO<sub>2</sub></td><td>ZnSiO<sub>3</sub></td></tr><tr><td>Co<sub>3</sub>O<sub>4</sub></td><td>K<sub>2</sub>S<sub>4</sub>O<sub>6</sub></td><td>Na<sub>2</sub>O.2SiO<sub>2</sub></td><td>SO<sub>3</sub></td><td>ZnSO<sub>4</sub></td></tr><tr><td>CoCO<sub>3</sub></td><td>K<sub>2</sub>SO<sub>3</sub>.H<sub>2</sub>O</td><td>Na<sub>2</sub>O.2TiO<sub>2</sub></td><td>Sr(NO<sub>2</sub>)<sub>2</sub></td><td>ZnSO<sub>4</sub>.2H<sub>2</sub>O</td></tr><tr><td>CoO</td><td>K<sub>2</sub>SO<sub>3</sub>.H<sub>2</sub>O</td><td>Na<sub>2</sub>O.3SiO<sub>2</sub></td><td>Sr(NO<sub>3</sub>)<sub>2</sub></td><td>ZnSO<sub>4</sub>.6H<sub>2</sub>O</td></tr><tr><td>CoO.Al<sub>2</sub>O<sub>3</sub></td><td>K<sub>2</sub>SO<sub>4</sub></td><td>Na<sub>2</sub>O.3TiO<sub>2</sub></td><td>Sr(NO<sub>3</sub>)<sub>2</sub>.4H<sub>2</sub>O</td><td>ZnSO<sub>4</sub>.7H<sub>2</sub>O</td></tr><tr><td>CoO.Cr<sub>2</sub>O<sub>3</sub></td><td>K<sub>2</sub>SO<sub>4</sub>.2MgSO<sub>4</sub></td><td>Na<sub>2</sub>O.TiO<sub>2</sub></td><td>Sr(OH)<sub>2</sub></td><td>ZnSO<sub>4</sub>.H<sub>2</sub>O</td></tr><tr><td>CoO.Fe<sub>2</sub>O<sub>3</sub></td><td>K<sub>2</sub>TiO<sub>3</sub></td><td>Na<sub>2</sub>O<sub>2</sub></td><td>Sr<sub>2</sub>CrO<sub>4</sub></td><td></td></tr><tr><td>CoO.TiO<sub>2</sub></td><td>K<sub>2</sub>Zn<sub>3</sub>(P<sub>2</sub>O<sub>7</sub>)<sub>2</sub>.3H<sub>2</sub>O</td><td>Na<sub>2</sub>P<sub>2</sub>O<sub>6</sub></td><td>Sr<sub>2</sub>Fe<sub>2</sub>O<sub>5</sub></td><td></td></tr><tr><td>CoSO<sub>4</sub></td><td>K<sub>3</sub>PO<sub>4</sub></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>5. Phosphides</td><td></td><td></td></tr><tr><td>AlP</td><td>CrP</td><td>K<sub>3</sub>P<sub>7</sub></td><td>Na<sub>3</sub>P<sub>7</sub></td><td>P<sub>4</sub>H<sub>2</sub></td></tr><tr><td>Ca<sub>3</sub>P<sub>2</sub></td><td>Cu<sub>3</sub>P</td><td>Mg<sub>3</sub>P<sub>2</sub></td><td>Ni<sub>2</sub>P</td><td>SiP</td></tr><tr><td>Co<sub>2</sub>P</td><td>CuP<sub>2</sub></td><td>Mn<sub>2</sub>P</td><td>Ni<sub>3</sub>P</td><td>Zn<sub>3</sub>P<sub>2</sub></td></tr><tr><td>CoP</td><td>Fe<sub>2</sub>P</td><td>Mn<sub>3</sub>P</td><td>Ni<sub>5</sub>P<sub>2</sub></td><td>ZnP<sub>2</sub></td></tr><tr><td>CoP<sub>3</sub></td><td>Fe<sub>3</sub>P</td><td>Mn<sub>3</sub>P<sub>2</sub></td><td>Ni<sub>6</sub>P<sub>5</sub></td><td></td></tr><tr><td>Cr<sub>12</sub>P<sub>7</sub></td><td>FeP</td><td>MnP</td><td>NiP<sub>2</sub></td><td></td></tr><tr><td>Cr<sub>3</sub>P</td><td>FeP<sub>2</sub></td><td>MnP<sub>3</sub></td><td>NiP<sub>3</sub></td><td></td></tr><tr><td></td><td></td><td>6. Sulfides</td><td></td><td></td></tr><tr><td>Al<sub>2</sub>S<sub>3</sub></td><td>Fe<sub>2</sub>S</td><td>MnS<sub>2</sub></td><td>NiS</td><td>Ti<sub>2</sub>S</td></tr><tr><td>CaS</td><td>Fe<sub>2</sub>S<sub>3</sub></td><td>Na<sub>2</sub>S</td><td>P<sub>2</sub>S<sub>3</sub></td><td>TiS</td></tr></tbody></table>
<table><caption>Table A1
(Continued)</caption><tbody><tr><td>Co<sub>3</sub>S<sub>4</sub></td><td>Fe<sub>7</sub>S<sub>8</sub></td><td>Na<sub>2</sub>S<sub>2</sub></td><td>P<sub>2</sub>S<sub>5</sub></td><td>TiS<sub>2</sub></td></tr><tr><td>Co<sub>9</sub>S<sub>8</sub></td><td>Fe<sub>9</sub>S<sub>8</sub></td><td>Na<sub>2</sub>S<sub>3</sub></td><td>P<sub>4</sub>S<sub>10</sub></td><td>TiS<sub>3</sub></td></tr><tr><td>CoS</td><td>FeS</td><td>Na<sub>2</sub>S<sub>4</sub></td><td>P<sub>4</sub>S<sub>3</sub></td><td>V<sub>2</sub>S<sub>3</sub></td></tr><tr><td>CoS<sub>2</sub></td><td>FeS<sub>2</sub></td><td>NaS</td><td>P<sub>4</sub>S<sub>5</sub></td><td>VS</td></tr><tr><td>Cr<sub>2</sub>S<sub>3</sub></td><td>K<sub>2</sub>S</td><td>Na<sub>2</sub>S<sub>2</sub></td><td>P<sub>4</sub>S<sub>6</sub></td><td>VS<sub>4</sub></td></tr><tr><td>CrAl<sub>2</sub>S<sub>4</sub></td><td>K<sub>2</sub>S<sub>2</sub></td><td>NH<sub>4</sub>HS</td><td>P<sub>4</sub>S<sub>7</sub></td><td>ZnCr<sub>2</sub>S<sub>4</sub></td></tr><tr><td>CrS</td><td>K<sub>2</sub>S<sub>3</sub></td><td>Ni<sub>3</sub>S<sub>2</sub></td><td>Sc<sub>2</sub>S<sub>3</sub></td><td>ZnS</td></tr><tr><td>Cu<sub>2</sub>S</td><td>K<sub>2</sub>S<sub>4</sub></td><td>Ni<sub>3</sub>S<sub>4</sub></td><td>ScS</td><td></td></tr><tr><td>Cu<sub>5</sub>FeS<sub>4</sub></td><td>K<sub>2</sub>S<sub>5</sub></td><td>Ni<sub>6</sub>S<sub>5</sub></td><td>SiS</td><td></td></tr><tr><td>CuFeS<sub>2</sub></td><td>MgS</td><td>Ni<sub>7</sub>S<sub>6</sub></td><td>SiS<sub>2</sub></td><td></td></tr><tr><td>CuS</td><td>MnS</td><td>Ni<sub>9</sub>S<sub>8</sub></td><td>SrS</td><td></td></tr><tr><td colspan="5">7. Nitrides</td></tr><tr><td>AlN</td><td>CuCN</td><td>KCN</td><td>NaCN</td><td>P<sub>3</sub>N<sub>5</sub></td></tr><tr><td>C<sub>2</sub>N<sub>4</sub>H<sub>4</sub></td><td>CuN<sub>3</sub></td><td>KN<sub>3</sub></td><td>NaN<sub>3</sub></td><td>ScN</td></tr><tr><td>Ca<sub>3</sub>N<sub>2</sub></td><td>Fe<sub>2</sub>N</td><td>KNH<sub>2</sub></td><td>NaNH<sub>2</sub></td><td>Si<sub>3</sub>N<sub>4</sub></td></tr><tr><td>CaCN<sub>2</sub></td><td>Fe<sub>4</sub>N</td><td>Mg<sub>3</sub>N<sub>2</sub></td><td>NaNH<sub>3</sub></td><td>Sr<sub>3</sub>N<sub>2</sub></td></tr><tr><td>Co<sub>3</sub>N</td><td>K<sub>3</sub>Co(CN)<sub>6</sub></td><td>Mn<sub>3</sub>N<sub>2</sub></td><td>NH<sub>4</sub>N<sub>3</sub></td><td>TiN</td></tr><tr><td>Cr<sub>2</sub>N</td><td>K<sub>3</sub>Fe(CN)<sub>6</sub></td><td>Mn<sub>4</sub>N</td><td>Ni(CN)<sub>2</sub></td><td>VN</td></tr><tr><td>CrN</td><td>K<sub>4</sub>Fe(CN)<sub>6</sub></td><td>Mn<sub>5</sub>N<sub>2</sub></td><td>Ni<sub>3</sub>N</td><td>Zn<sub>3</sub>N<sub>2</sub></td></tr></tbody></table>

ORCID iDs

G. S. Pekmezci https://orcid.org/0000-0001-6230-1093

References

Alessi, M., Pudritz, R. E., & Cridland, A. J. 2017, MNRAS, 464, 428

Alibert, Y., Mousis, O., & Benz, W. 2005a, ApJL, 622, L145

Alibert, Y., Mousis, O., & Benz, W. 2005b, ApJ, 622, L145

Ali-Dib, M. 2017, MNRAS, 464, 4282

Ali-Dib, M., Mousis, O., Petit, J.-M., & Lunine, J. I. 2014, ApJ, 785, 125

Artur de la Villarmois, E., Jørgensen, J. K., Kristensen, L. E., et al. 2019, A&A, 626, 19

Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481

Bézard, B., Lellouch, E., Strobel, D., Maillard, J.-P., & Drossart, P. 2002, Icar, 159, 95

Bitsch, B., & Battistini, C. 2020, A&A, 633, A10

Bond, J. C., Lauretta, D. S., & OBrien, D. P. 2010a, Icar, 205, 321

Bond, J. C., O'Brien, D. P., & Lauretta, D. S. 2010b, ApJ, 715, 1050

Booth, R. A., & Clarke, C. J. 2018, MNRAS, 473, 757

Booth, R. A., & Ilee, J. D. 2019, MNRAS, 487, 3998

Carter-Bond, J. C., O'Brien, D. P., & Raymond, S. N. 2012, ApJ, 760, 13

Cassen, P. 1996, M&PS, 31, 793

Cridland, A. J., Bosman, A. D., & van Dishoeck, E. F. 2020, A&A, 635, A68

Cridland, A. J., Eistrup, C., & van Dishoeck, E. F. 2019a, A&A, 627, A127

Cridland, A. J., Pudritz, R. E., & Alessi, M. 2016, MNRAS, 461, 3274

Cridland, A. J., Pudritz, R. E., Birnstiel, T., Cleeves, L. I., & Bergin, E. A. 2017, MNRAS, 469, 3910

Cridland, A. J., van Dishoeck, E. F., Alessi, M., & Pudritz, R. E. 2019b, A&A, 632, A63

Dodson-Robinson, S. E., Willacy, K., Bodenheimer, P., Turner, N. J., & Beichman, C. A. 2009, Icar, 200, 672

Ehrenfreund, P., & Schutte, W. A. 2000, AdSpR, 25, 2177

Eistrup, C., Walsh, C., & van Dishoeck, E. F. 2016, A&A, 595, A83

Eistrup, C., Walsh, C., & van Dishoeck, E. F. 2018, A&A, 613, A14

Elser, S., Meyer, M. R., & Moore, B. 2012, Icar, 221, 859

Fegley, B. 1988, in Workshop on the Origins of Solar Systems, ed. J. A. Nuth & P. Sylvester (Houston, TX: Lunar and Planetary Institute), 51

Fegley, B. 2000, SSRv, 92, 177

Fegley, B., & Lodders, K. 1994, Icar, 110, 117

Fletcher, L. N., Orton, G. S., Teanby, N. A., & Irwin, P. G. J. 2009, Icar, 202, 543

Fortney, J. J. 2012, ApJL, 747, L27

Fossati, L., Bagnulo, S., Elmasli, A., et al. 2010, ApJ, 720, 872

Frerkling, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 590

Gibb, E. L., Whittet, D. C. B., Schutte, W. A., et al. 2000, ApJ, 536, 347

Goettel, K. A. 1988, Mercury (A89-43751 19-91) (Tucson, AZ: Univ. of Arizona Press), 613

Guillot, T. 2008, PhyS, T130, 014023

Hebb, L., Collier-Cameron, A., Loeillet, B., et al. 2009, ApJ, 693, 1920

Helling, C., Woitke, P., Rimmer, P. B., et al. 2014, Life, 4, 142

Hersant, F., Gautier, D., & Hur, J.-M. 2001, ApJ, 554, 391

Ilee, J. D., Boley, A. C., Caselli, P., et al. 2011, MNRAS, 417, 2950

Johnson, T. V., Mousis, O., Lunine, J. I., & Madhusudhan, N. 2012, ApJ, 757, 192

Kozlovskaja, S. V., & Pechernikova, G. V. 1992, EM&P, 57, 225

Kuchner, M. J., & Seager, S. 2005, Extrasolar Carbon Planets, arXiv:astro-ph/0504214

Léger, A., Selsis, F., Sotin, C., et al. 2004, Icar, 169, 499

Lewis, J. S. 1974, Sci, 186, 440

Li, C., Ingersoll, A., Bolton, S., et al. 2020a, NatAs, 4, 609

Li, M., Huang, S., Petaev, M. I., Zhu, Z., & Steffen, J. H. 2020b, MNRAS, 495, 2543

Lide, D. R. 2002, CRC Handbook of Chemistry and Physics: A Ready-reference Book of Chemical and Physical Data (Boca Raton, FL: CRC Press)

Lin, D. N. C. 2008, SciAm, 298, 50

Lin, D. N. C., & Papaloizou, J. 1986, ApJ, 309, 846

Line, M. R., Liang, M. C., & Yung, Y. L. 2010, ApJ, 717, 496

Lodders, K. 2003, ApJ, 591, 1220

Lodders, K. 2010, in Formation and Evolution of Exoplanets, ed. R. Barnes (New York: Wiley), 157

Lodders, K., & Fegley, B. 1994, Icar, 112, 368

Lodders, K., & Fegley, B. 2002, Icar, 155, 393

Lunine, J. I., & Stevenson, D. J. 1985, ApJS, 58, 493

Madhusudhan, N. 2012, ApJ, 758, 21

Madhusudhan, N., Harrington, J., Stevenson, K. B., et al. 2011a, Natur, 469, 64

Madhusudhan, N., Mousis, O., Johnson, T. V., & Lunine, J. I. 2011b, ApJ, 743, 191

Marboeuf, U., Thiabaud, A., Alibert, Y., Cabral, N., & Benz, W. 2014, A&A, 570, A35

Miyazaki, Y., & Korenaga, J. 2020, arXiv:2004.13911

Mordasini, C., Alibert, Y., Klahr, H., & Henning, T. 2012, A&A, 547, A111

Mordasini, C., van Boekel, R., Mollière, P., Henning, T., & Benneke, B. 2016, ApJ, 832, 32

Morgan, J. W., & Anders, E. 1980, PNAS, 77, 6973

Moriarty, J., Madhusudhan, N., & Fischer, D. 2014, ApJ, 787, 10

Moses, J. I., Madhusudhan, N., Visscher, C., & Freedman, R. S. 2013, ApJ, 763, 26

Moses, J. I., Visscher, C., Fortney, J. J., et al. 2011, ApJ, 737, 37

Mousis, O., Aguichine, A., Atkinson, D. H., et al. 2020, SSRv, 216, 77

Mousis, O., & Gautier, D. 2004, P&SS, 52, 361

Mousis, O., Lunine, J. I., Madhusudhan, N., & Johnson, T. V. 2012, ApJL, 751, L7
Mousis, O., Lunine, J. I., Petit, J.-M., et al. 2011, ApJ, 727, 7

Mousis, O., Lunine, J. I., Thomas, C., et al. 2009a, ApJ, 691, 1780

Mousis, O., Lunine, J. I., Tinetti, G., et al. 2009b, A&A, 507, 1671

Mousis, O., Marboeuf, U., Lunine, J. I., et al. 2009c, ApJ, 696, 1348

Mousis, O., Ronnet, T., & Lunine, J. I. 2019, ApJ, 875, 6

Öberg, K. I., Murray-Clay, R., & Bergin, E. A. 2011, ApJL, 743, 5

Ohishi, M., Irvine, W. M., & Kaifu, N. 1992, IAUS, 150, 171

Pasek, M. A., Milsom, J. A., Ciesla, F. J., et al. 2005, Icar, 175, 1

Pekmezci, G. S., Johnson, T. V., Lunine, J. I., & Mousis, O. 2019, ApJ, 887, 10

Piso, A.-M. A., Pegues, J., & Öberg, K. I. 2016, ApJ, 833, 8

Price, E. M., Cleeves, L. I., & Öberg, K. I. 2020, ApJ, 890, 154

Prinn, R. G. 1993, in Protostars and Planets III, ed. E. H. Levy & J. I. Lunine (Tucson, AZ: Univ. of Arizona Press), 1005

Prinn, R. G., & Fegley, B. 1981, ApJ, 249, 308

Raymond, S. N., Quinn, T., & Lunine, J. I. 2005, ApJ, 632, 670

Rivoldini, A., van Hoolst, T., Verhoeven, O., Mocquet, A., & Dehant, V. 2011, Icar, 213, 451

Roine, A. 2006, HSC Chemistry 7.1 Help Menu HSC Chemistry 7.1 Commercial Package

Ronco, M. P., Thiabaud, A., Marboeuf, U., et al. 2015, BAAA, 57, 251

Santos, N. C., Adibekyan, V., Dorn, C., et al. 2017, A&A, 608, A94

Schaefer, L., Lodders, K., & Fegley, B. 2012, ApJ, 755, 41

Schneider, J. 2020, The Extrasolar Planetary Encyclopedia. Accessed on 2020 November 4 www.exoplanet.eu/catalog/

Seager, S., & Elkins-Tanton, L. 2008, ApJ, 688, 628

Stolper, E. 1980, in 11th Lunar and Planetary Science Conf. Proc. (New York: Pergamon), 235

Strom, R. G., & Sprague, A. L. 2003, Exploring Mercury: The Iron Planet (London: Springer)

Taylor, G. J. 2012, The Composition of Bulk Silicate Mars The Mantle of Mars: Insights from Theory, Geophysics, High-Pressure Studies, and Meteorites, LPI Contribution #1684, ID, 10

Taylor, S. R. 1991, Metic, 26, 267

Thiabaud, A., Marboeuf, U., Alibert, Y., Leya, I., & Mezger, K. 2015, A&A, 580, 7

Treiman, A. H. 2009, in Venus Geochemistry: Progress, Prospects, and New Missions (Houston, TX: Lunar and Planetary Institute), 47

Visscher, C., Lodders, K., & Fegley, B., Jr. 2006, ApJ, 648, 1181

Visscher, C., & Moses, J. I. 2011, ApJ, 738, 12

Visscher, C., Moses, J. I., & Saslow, S. A. 2010, Icar, 209, 602

Wahl, S. M., Hubbar, W., Militzer, B., et al. 2017, GeoRL, 44, 4649

Ward, W. R. 1997, Icar, 126, 261

Wehrstedt, M., & Gail, H.-P. 2003, A&A, 410, 917

Wong, M. H., Lunine, J. I., Atreya, S. K., et al. 2008, RvMG, 68, 219

Wong, M. H., Mahay, P. R., Atreya, S. K., Niemann, H. B., & Owen, T. C. 2004, Icar, 171, 153