# Formation of the Trappist-1 system in a dry protoplanetary disk

Antoine Schneeberger, Olivier Mousis, Magali Deleuil, Jonathan I Lunine

► To cite this version:

Antoine Schneeberger, Olivier Mousis, Magali Deleuil, Jonathan I Lunine. Formation of the Trappist-1 system in a dry protoplanetary disk. Astronomy & Astrophysics - A&A, 2024, Astronomy & Astrophysics, 682, pp.L10. 10.1051/0004-6361/202348309. hal-04452204

HAL Id: hal-04452204

https://hal.science/hal-04452204v1

Submitted on 12 Feb 2024

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
LETTER TO THE EDITOR

# Formation of the Trappist-1 system in a dry protoplanetary disk

Antoine Schneeberger<fnref n="1" />, Olivier Mousis<fnref n="1" />,<fnref n="2" />, Magali Deleuil<fnref n="1" />,<fnref n="2" />, and Jonathan I. Lunine<fnref n="3" />

Received 18 October 2023 / Accepted 18 January 2024

ABSTRACT

A key feature of the Trappist-1 system is its monotonic decrease in bulk density with growing distance from the central star, which indicates an ice mass fraction that is zero in the innermost planets, b and c, and about 10% in planets d through h. Previous studies suggest that the density gradient of this system could be due to the growth of planets from icy planetesimals that progressively lost their volatile content during their inward drift through the protoplanetary disk. Here we investigate the alternative possibility that the planets formed in a dry protoplanetary disk populated with pebbles made of phyllosilicates, a class of hydrated minerals with a water fraction possibly exceeding 10 wt%. We show that the dehydration of these minerals in the inner regions of the disk and the outward diffusion of the released vapor up to the ice-line location allow the condensation of ice onto grains. Pebbles with water mass fractions consistent with those of planets d-h would have formed at the snow-line location. In contrast, planets b and c would have been accreted from drier material in regions closer to the star than the phyllosilicate dehydration line.

**Key words.** methods: numerical – planets and satellites: composition – planets and satellites: formation – protoplanetary disks

1. Introduction

Because of its compactness and density gradient, the Trappist-1 exoplanetary system is the closest analog to the well-characterized Galilean moon system. It is also an outstanding laboratory for exploring planet formation around ultracool M dwarfs. The Trappist-1 system contains seven exoplanets in a resonance chain, three of which are within the habitable zone (Gillon et al. 2016, 2017), with semimajor axes ranging between 0.015 and 0.062 AU (Luger et al. 2017; Agol et al. 2021). This resonance chain could be a sign of a past inward drift of the whole system before orbital capture (Ormel et al. 2017; Bitsch et al. 2019; Coleman et al. 2019; Izidoro et al. 2021; Madeira et al. 2021).

The measured densities of the Trappist-1 planets are among the most accurate ones available for an exoplanetary system (Luger et al. 2017). These densities have been used by Acuña et al. (2021) to assess their interiors and compositions. The authors show evidence for an increasing water mass fraction (WMF) of the planets with increasing distance from the host star. The planets can be separated into two groups, with planets b–c containing little to zero water, and planets d–h presenting significant amounts of water (see Table 1). The significant amount of water estimated in the five outermost planets has been interpreted as a sign that they formed beyond the position of the snow line (SL), within the water-rich protoplanetary disk (PPD), prior to their inward migration (Ormel et al. 2017; Johansen et al. 2023). Such a mechanism requires a rapid growth of the planets, which could be eased by the formation of planetesimals via streaming instability triggered by a solid-to-gas ratio above unity at the SL due to ice enrichment (Ormel et al. 2017).

However, red dwarfs produce far-ultraviolet (FUV) and extreme-ultraviolet (XUV) irradiation during their early evolution (Shkolnik & Barman 2014). This can cause the photolytic destruction of the water present in the PPD or in the material originating from the interstellar medium and accreted by the PPD (Shkolnik & Barman 2014; Loyd et al. 2018; Youngblood et al. 2021). Furthermore, repeated episodes of accretional outbursts by Trappist-1 could have produced rapid temperature rises in the PPD that would dehydrate pebbles (Houge & Krijt 2023). Finally, sequestration of oxygen in rock-and metal-forming species must be considered when the carbon-to-oxygen (C/O) ratio exceeds the solar value (0.55; Lodders 2021) in the PPD. When the C/O ratio reaches or exceeds 0.8, the depletion of oxygen is severe; it is bound up in rock, metal oxides, CO, and CO₂, hindering the production of water (Pekmezci et al. 2019).

In systems where one or both of the water loss processes occur, an interesting alternative is to consider the case where water would be embedded in phyllosilicate particles. Phyllosilicates form a large group of hydrous minerals, which are structured in layers bounded by cation-anion couples that contain OH groups in their crystalline structure. They can contain H₂O molecules between layers, and dehydration experiments show that serpentinite minerals (one of the most common forms) can hold up to ~13 wt% of water (Bezacier et al. 2010; Chauviré et al. 2021; Pettke & Bretscher 2022).

Detection of phyllosilicate grains was claimed in the interstellar medium and PPDs (Zaikowski et al. 1975; Whittet et al. 1997; Reach et al. 2009). However, this evidence has since been reinterpreted as a mixture of amorphous silicates and water ice (Gibb et al. 2004; Potapov et al. 2021). On the other hand,

<footnote n="1">¹ Aix-Marseille Université, CNRS, CNES, Institut Origines, LAM, Marseille, France
e-mail: antoine.schneeberger@lam.fr</footnote>

<footnote n="2">² Institut Universitaire de France (IUF), Paris, France</footnote>

<footnote n="3">³ Department of Astronomy, Cornell University, Ithaca, NY, USA</footnote>
phyllosilicates are found in CI and CM chondrites, which have WMFs as high as 10% (Beck et al. 2010, 2014; Alexander 2019). The current consensus is that phyllosilicates found in CM chondrites are formed by aqueous alteration of silicates in water-rich planetesimals (Bischoff 1998; Suttle et al. 2021). Such ice-rich solids can form in the vicinity of the water SL and become precursors of phyllosilicate-rich pebbles (Bischoff 1998) via collisional evolution.

However, hydrous mineral inclusions within chondritic meteorites exhibit a distinct D/H ratio and hydration heterogeneity, favoring a primitive origin (Metzler et al. 1992; Ciesla et al. 2003; Piani et al. 2021). This supports models that show that micrometric phyllosilicate grains can result from the alteration of silicate grains by water vapor over a time period of ~100 kyr in the PPD (D'Angelo et al. 2019; Thi et al. 2020). This is also consistent with the idea that phyllosilicate nanocrystals could form in the wake of shock waves crossing the PPD (Ciesla et al. 2003). These shock waves would produce a temperature rise of up to 2000 K (Ciesla et al. 2003; Burkhardt et al. 2019), which would evaporate water ice and melt silicates. As the molten silicate recrystallized, it would react with the surrounding water vapor to form phyllosilicate nanocrystals within a few days of the shock wave. Both processes require the presence of water vapor, which implies that they occurred early in the history of the PPD.

The possible dearth of water as ice in red dwarf PPDs suggests an alternative scenario where (i) the Trappist system formed in a dry PPD devoid of volatiles, and (ii) water was only delivered to the inner disk in the form of phyllosilicate minerals. Figure 1 qualitatively describes this scenario in which three distinct regions coexist in the disk: an inner region occupied by dehydrated phyllosilicates that is located inside the phyllosilicate dehydration line (PDL), a second region containing hydrated phyllosilicates and bracketed by the PDL and the SL, and an outer region populated by phyllosilicates coated with ice and extending from the SL. In this scenario, phyllosilicates dehydrate in the inner part of the disk due to increasing temperature and release water vapor at the PDL, which viscously diffuses inward and outward, eventually crossing the SL, where it condenses onto grains. This formation mechanism has been successfully applied to the Galilean moons to account for their observed density gradient (Mousis et al. 2023).

We have investigated the evolution of the WMF of phyllosilicate-rich pebbles during their transport and growth in a dry PPD. We used the one-dimensional volatile evolution model in a PPD from Aguichine et al. (2020, 2022) and Schneeberger et al. (2023), to which we added the dehydration of phyllosilicates at the PDL location. We tracked the WMF evolution of pebbles at the SL location, assuming a starting WMF in the phyllosilicate grains in the 1–10% range. We then studied the impact of the proposed scenario on the formation history of the Trappist-1 system.

## 2. Disk model

Our PPD model is based on the PPD model presented in Aguichine et al. (2020, 2022), Mousis et al. (2020), and Schneeberger et al. (2023). To size the model to the Trappist-1 system, the central star mass was scaled down to 0.0802 $M_\odot$ (Gillon et al. 2016). We also opted to linearly scale the disk mass and accretion rate from $0.1 M_\odot$ and $10^{-7.6} M_\odot \text{yr}^{-1}$ for the PPD (Hartmann et al. 1998) to $0.1 M_\star$ and $10^{-7.6} M_\star \text{yr}^{-1}$ for the Trappist-1 system, with $M_\star$ being the mass of Trappist-1.

<table><caption>**Table 1.** Estimated masses and WMFs for Trappist-1 planets (Acuña et al. 2021).</caption><thead><tr><th>Planet</th><th>Mass ($M_⊕$)</th><th>WMF</th></tr></thead><tbody><tr><td>b</td><td>1.375 ± 0.041</td><td>(3.1<sup>+0.5</sup><sub>-3.1</sub>) × 10<sup>-5</sup></td></tr><tr><td>c</td><td>1.300 ± 0.036</td><td>(0.0<sup>+4.4</sup><sub>-0.0</sub>) × 10<sup>-6</sup></td></tr><tr><td>d</td><td>0.388 ± 0.007</td><td>0.084 ± 0.071</td></tr><tr><td>e</td><td>0.699 ± 0.013</td><td>0.094 ± 0.067</td></tr><tr><td>f</td><td>1.043 ± 0.019</td><td>0.105 ± 0.073</td></tr><tr><td>g</td><td>1.327 ± 0.024</td><td>0.119 ± 0.080</td></tr><tr><td>h</td><td>0.327 ± 0.012</td><td>0.081<sup>+0.059</sup><sub>-0.081</sub></td></tr></tbody></table>

Our PPD model follows the differential equation (Lynden-Bell & Pringle 1974)

$$ \frac{\partial \Sigma_g}{\partial t} = \frac{3}{r} \frac{\partial}{\partial r} \left[ r^{1/2} \frac{\partial}{\partial r} \left( r^{1/2} \Sigma_g v \right) \right], \quad (1) $$

where $\Sigma_g$ is the gas surface density and $v$ is the gas viscosity. Our model assumes invariance in the orbital direction and hydrostatic equilibrium in the azimuthal direction. The viscosity is given by (Shakura & Sunyaev 1973)

$$ v = \alpha \frac{C_s^2}{\Omega_k}, \quad (2) $$

where $\alpha$ is the viscosity parameter set to the canonical value of $10^{-3}$ (Shakura & Sunyaev 1973; Lynden-Bell & Pringle 1974), $C_s$ the sound of speed, and $\Omega_k$ the Keplerian pulsation.

The initial disk profile was computed from the self-similar solution derived by Lynden-Bell & Pringle (1974):

$$ \Sigma_{g,0} = \frac{\dot{M}_{\text{acc},0}}{3\pi v} \exp\left[-\left(\frac{r}{r_c}\right)^{0.5}\right], \quad (3) $$

with $\dot{M}_{\text{acc},0}$ the initial accretion rate onto the PPD and $r_c$ the centrifugal radius of the disk. The centrifugal radius was self-consistently computed at its initialization as described in Aguichine et al. (2020).

Figure 2 shows the surface density, midplane pressure, and temperature profiles derived from our Trappist-1 PPD model and from the PPD model used by Schneeberger et al. (2023) after 0.01, 0.1, and 1 Myr of disk evolution. The comparison shows that the Trappist-1 PPD is cooler and less dense by a factor of ~5–10 in the 0.5–10 AU region of the disk. The midplane pressure of the Trappist-1 PPD is also ~30–70 times lower than in the PPD model of Schneeberger et al. (2023).

The evolution of the different water phases follows the advection-diffusion equation (Birnstiel et al. 2012; Desch et al. 2017):

$$ \frac{\partial \Sigma_i}{\partial t} + \frac{1}{r} \frac{\partial}{\partial r} \left[ r \left( \Sigma_i v_i - D_i \Sigma_g \frac{\partial}{\partial r} \left( \frac{\Sigma_i}{\Sigma_g} \right) \right) \right] - \dot{Q}_i = 0, \quad (4) $$

where $\Sigma_i$ is the water surface density of phase i, either in vapor or ice form, $v_i$ and $D_i$ are the drift velocity and diffusivity, respectively, and $Q_i$ is the source or sink term of phase i, which takes the condensation and vaporization of water into account. The evolution of grains is governed by the two-population algorithm, which provides their mean drift speed and diffusivity (Birnstiel et al. 2012).
<figure><img src="image_2.png" /><figcaption>Fig. 1. Structure of Trappist-1's PPD. Phyllosilicate grains (purple dots) drift inward in the disk. When these grains cross the PDL (orange bar), they dehydrate and release water as vapor. Only silicate grains remain inside the PDL (pink dots). The vapor viscously diffuses inward and outward and eventually crosses the SL (blue bar), where it condenses and coats phyllosilicate grains with ice (blue shell).</figcaption></figure>

The water source and sink terms are governed by its equilib- rium pressure (Wagner et al. 2011):

$$ P_{\text{eq}} = P_{\text{tp}} \exp\left(\frac{T_{\text{tp}}}{T}\right) \left[ -21.2144006\left(\frac{T_{\text{tp}}}{T}\right)^{1/300} \right. \\ \left. \qquad + 27.3203819\left(\frac{T_{\text{tp}}}{T}\right)^{2.10666667} \right. \\ \left. \qquad -6.10598130\left(\frac{T_{\text{tp}}}{T}\right)^{1.70333333} \right], \quad (5) $$

where $P_{\text{tp}}$ and $T_{\text{tp}}$ are the water triple point pressure (611.65 Pa) and temperature (273.16 K), respectively. The equilibrium pres- sure is compared with the partial pressure of water:

$$ P_v = \frac{\Sigma_v RT}{\sqrt{2\pi\mu_{H_2O}}H}, \quad (6) $$

where $\Sigma_v$ is the water vapor surface density, $R$ is the perfect gas constant, $T$ is the disk temperature, $\mu_{H_2O}$ is the water molar mass, and $H$ is the disk scale height. If the vapor partial pressure is higher than the equilibrium one, then it condenses onto the peb- bles and forms an icy shell. Condensation then results in a source term for water ice of

$$ \dot{Q}_{\text{ice}} = \min\left(\left(P_v - P_{\text{eq}}\right) \frac{2H\mu_{H_2O}}{RT\Delta t}; \frac{\Sigma_{v,i}}{\Delta t}\right). \quad (7) $$

Conversely, if the vapor partial pressure is lower than the equilibrium pressure, ice vaporizes, resulting in a sink term for water ice of

$$ \dot{Q}_{\text{ice}} = -\min\left(\sqrt{\frac{8\pi\mu_i}{RT}} \frac{3}{\pi\bar{a}\bar{\rho}} P_{\text{eq}}\Sigma_{\text{ice}}; \frac{\Sigma_{\text{ice}}}{\Delta t}\right), \quad (8) $$

where $\bar{a}$ and $\bar{\rho}$ are the grains' mean size and mean density as given by the two-population algorithm, respectively.

We added a module that takes the dehydration of phyllosili- cates into account. We considered that phyllosilicates dehydrate when they reach a temperature of 500 K in the disk (Mousis et al. 2023). Phyllosilicate-rich particles release water vapor at the PDL location in the PPD, which corresponds to the following source term for water vapor:

$$ \dot{Q} = \text{WMF}_{\text{phyllo}} \frac{\Sigma_p}{dt}, \quad (9) $$

where $\text{WMF}_{\text{phyllo}}$ is the WMF in the phyllosilicates, which is var- ied in the 1–10% range, and $\Sigma_p$ is the surface density of phyllosil- icates. We assumed that the phyllosilicate dehydration timescale is much shorter than the drift timescale (Ciesla & Lauretta 2005; Mousis et al. 2023). All the water stored in the phyllosilicate dust or pebbles is released when they cross the PDL location. The PPD model was initialized with phyllosilicate grains of 0.1 µm in size. The initial phyllosilicate abundance, $x_p^0 = \Sigma_p^0/\Sigma_g^0$, was set to $3.212 \times 10^{-3}$ (Lodders et al. 2009).

### 3. Results

Figure 3 displays water ice and vapor enrichment profiles in the PPD, with respect to the initial water abundance, after 0.01, 0.1, and 1 Myr of disk evolution. The initial water abundance is defined as

$$ x_{H_2O}^0 = \text{WMF}_{\text{phyllo}} x_p^0. \quad (10) $$

Since the disk only contains phyllosilicate at $t = 0$, the starting water abundance only depends on the water contained in phyl- losilicate particles. $\text{WMF}_{\text{phyllo}}$ of 1%, 5%, and 10% were used in the simulations. The distance between the SL and the central star increases from 0.67 to 0.85 AU, 0.53 to 0.66 AU, and 0.35 to 0.43 AU when the $\text{WMF}_{\text{phyllo}}$ value is ranged between 1% and 10% at 0.01, 0.1, and 1 Myr of PPD evolution, respectively.

The increase in this distance range is the direct consequence of the decrease in the amount of water vapor released from the
<figure><img src="image_3.png" /><figcaption>Fig. 2. From top to bottom: Temperature, pressure, and surface density profiles of the Trappist-1 PPD (solid lines) superimposed onto those of the PPD taken from Schneeberger et al. (2023; dashed lines) after 0.01, 0.1, and 1 Myr of evolution.</figcaption></figure>

pebbles into the PPD. Vapor accumulates at the PDL because its diffusion rate in the PPD is much lower than its production rate. A fraction of this vapor diffuses outward and crosses the SL.

Figure 3 shows that the water ice enrichment reaches a maximum at ~0.4, 1.3, and 2 times the initial water abundance before 0.01 Myr of PPD evolution, assuming a $WMF_{phyllo}$ of 1, 5, and 10%, respectively. Between 0.01 and 0.1 Myr of PPD evolution, the water enrichment decreases by a factor of ~20–30 down to 0.02, 0.04, and 0.07 times the initial water abundance. After 0.1 Myr, the ice enrichment remains almost constant up to 1 Myr, which corresponds to the end of our simulations.

Figure 4 represents the time evolution of the pebbles’ WMF at the SL location when $WMF_{phyllo}$ varies between 1 and 10%. The dark orange area corresponds to the range of nominal WMFs (8.1–11.9%) computed by Acuña et al. (2021) for Trappist-1 d, e, f, g, and h. The light orange area brackets the $1\sigma$ error bars (0–19.9%) associated with those measurements. The WMF of pebbles takes both the fraction of water in ice form and water

trapped in phyllosilicates into account:

$$WMF = \frac{WMF_{phyllo}\Sigma_p + \Sigma_{ice}}{\Sigma_p + \Sigma_{ice}}, \quad (11)$$

with $\Sigma_{ice}$ the surface density of ice in the PPD. Figure 4 shows that even in the absence of ice in the PPD at initialization, the pebbles’ WMF is 2–3 times larger than $WMF_{phyllo}$ after only 0.01 Myr of evolution. The pebbles’ WMF ranges between 2 and 28% when $WMF_{phyllo}$ is varied between 1 and 10%. As the PPD depletes, the amount of water vapor released from the phyllosilicates decreases, reducing the amount of ice condensed at the SL. Consequently, the pebbles’ WMF converges toward $WMF_{phyllo}$ after 0.07 Myr of PPD evolution. This implies that only cases with $WMF_{phyllo}$ in the 8–10% range match the WMF of planets d–h after 0.07 Myr of PPD evolution. Interestingly, cases with pebbles that assume a $WMF_{phyllo}$ in the 4–7% range also match the WMF of planets d–h before 0.07 Myr of PPD evolution. However, after this epoch, water ice is too depleted in the PPD.

## 4. Discussion

The scenario we have explored for the Trappist-1 system divides the PPD into three distinct regions: (i) an inner region interior to the PDL in which solids consist of pebbles made of dehydrated silicates, (ii) an intermediary region located between the PDL and the SL populated with pebbles assembled from hydrated phyllosilicates, and (iii) a region located beyond the SL populated with pebbles consisting of hydrated phyllosilicate pebbles coated with water ice.

In this scenario, the building blocks of planets b and c would have formed in regions interior to the PDL, while those of planets d–h would have formed at the SL. Planets b and c would have grown close to the star, possibly at their present-day location. On the other hand, planets d–h would have grown farther from the star than their current location (Ormel et al. 2017; Coleman et al. 2019; Huang & Ormel 2022; Ogihara et al. 2022). They could have migrated inward to their present-day location after formation. The inward migration of all planets could have led to the observed mean motion resonance chain formed by the seven planets (Ormel et al. 2017; Bitsch et al. 2019; Coleman et al. 2019; Izidoro et al. 2021).

The dehydration temperature of phyllosilicates at pressure conditions relevant to those of PPDs is not well established. In our study we chose a dehydration temperature of 500 K to explore its effect on the pebbles’ WMF (Mousis et al. 2023). Varying the dehydration temperature between 400 K and 600 K has a significant impact on the timescales and extents of ice enrichment in the PPD. A higher dehydration temperature results in a PDL located closer to the star than in cases with lower temperatures. Therefore, colder dehydration temperatures reduce the distance between the PDL and SL, enhance the amount of vapor diffusing throughout the SL, and lead to higher abundances of solid water. A dehydration temperature of 400 K induces an ice enrichment up to ~20 times the initial water abundance after 0.01 Myr of PPD evolution, assuming $WMF_{phyllo}$ = 10% in phyllosilicates. In contrast, only an ice enrichment of ~0.5 is generated in the PPD, assuming a dehydration temperature of 600 K. These results show that our simulations are highly sensitive to the dehydration temperature. This variability impacts the evolution of the pebbles’ WMF at the SL. Consequently, a lower dehydration temperature enables the formation of pebbles with WMFs consistent with those of planets d–h, assuming they
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** From left to right: Enrichment profiles of water vapor (orange lines) and ice (blue lines) normalized with the initial abundance of water in phyllosilicates at *t* = 0.01, 0.1, and 1 Myr. The initial WMFs in phyllosilicates here are 1%, 5%, and 10%. The PDL and SL are represented as brown and blue bars, respectively.</figcaption></figure>

start with WMF<sub>phyllo</sub> lower than the nominal value (8%) after 1 Myr of PPD evolution. Conversely, a higher dehydration temperature requires higher WMF<sub>phyllo</sub> values to reproduce the WMFs of planets d-h after 1 Myr of PPD evolution.

To compare the sensitivity of our results against the PPD
parameters, we reduced the PPD’s initial mass from 10⁻¹ to
10⁻² M<sub>★</sub>. A less massive PPD has a lower density, resulting in
reduced temperature and diffusion factor profiles. As a conse-
quence, despite a reduced distance between the SL and the PDL
of ~0.3 AU, there is a decrease in the amount of vapor cross-
ing the SL and condensing into ice. The water enrichment at the
SL is reduced by a factor of ~200 after 0.1 Myr of PPD evo-
lution compared with the nominal enrichment, assuming a 10%
WMF<sub>phyllo</sub>. After 1 Myr of PPD evolution, the disk dissipates
almost completely and very few phyllosilicate grains remain,
resulting in an even lower water vapor production. This indi-
cates that the PPD has a somewhat significant mass (allowing it
to generate pebbles) that is consistent with the WMFs of planets
d-h.

The proposed scenario relies on water loss by photolysis
driven by FUV/XUV radiation from Trappist-1, whose inten-
sity and time evolution remain unknown. If the FUV/XUV radi-
ation levels were very low, it is likely that primordial water
ice would remain abundant in the PPD. On the other hand,
an extreme FUV/XUV radiation episode of sufficient duration
would destroy the vapor released by phyllosilicate grains and
pebbles, thereby preventing any ice formation at the SL. Con-
sequently, a dry disk model represents an alternative scenario to
the formation of the Trappist planets from volatile-rich material
(Ormel et al. 2017; Bitsch et al. 2019). An observational test to
investigate a dry disk scenario would be the discovery of a planet
more distant than Trappist-1 h but with a density close to those
of Trappist-1 b-c. This would imply that this planet formed from

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Time evolution of the pebbles’ WMFs at the SL location. The phyllosilicate WMF is varied between 1 and 10%. The range of the WMF central values of the different Trappist-1 planets is bracketed by the dark orange box (Acuña et al. 2021), and the 1σ errors are shown in light orange.</figcaption></figure>

phyllosilicate-rich building blocks beyond the SL. In contrast,
the same planet formed in a wet disk would have a density sim-
ilar to those of planets d-h. Both scenarios could have been at
play during the Trappist system formation. In this case, the inner planets would have formed from pebbles made of dehydrated phyllosilicates. The outer planets would have accreted from ice-rich pebbles sourced from (i) water vapor condensation after outward diffusion from the inner regions and (ii) pristine ice delivered from the outer regions.

**Acknowledgements.** The project leading to this publication has received funding from the Excellence Initiative of Aix-Marseille Université-A*Midex, a French "Investissements d'Avenir program" AMX-21-IET-018. This research holds as part of the project FACOM (ANR-22-CE49-0005-01_ACT) and has benefited from a funding provided by l'Agence Nationale de la Recherche (ANR) under the Generic Call for Proposals 2022. We acknowledge Maxime Pineau for helpful discussions about phyllosilicate properties and experiments.

## References

Acuña, L., Deleuil, M., Mousis, O., et al. 2021, A&A, 647, A53

Agol, E., Dorn, C., Grimm, S. L., et al. 2021, Planet. Sci. J, 2, 1

Aguichine, A., Mousis, O., Devouard, B., et al. 2020, ApJ, 901, 97

Aguichine, A., Mousis, O., & Lunine, J. I. 2022, Planet. Sci. J, 3, 141

Alexander, C. M. O. 2019, Geochim. Cosmochim. Acta, 254, 246

Beck, P., Quirico, E., Montes-Hernandez, G., et al. 2010, Geochim. Cosmochim. Acta, 74, 4881

Beck, P., Garenne, A., Quirico, E., et al. 2014, Icarus, 229, 263

Bezacier, L., Reynard, B., Bass, J. D., et al. 2010, Earth Planet Sci. Lett., 289, 198

Birstiel, T., Klahr, H., & Ercolano, B. 2012, A&A, 539, A148

Bischoff, A. 1998, Meteor. Planet. Sci., 33, 1113

Bitsch, B., Raymond, S. N., & Izidoro, A. 2019, A&A, 624, A109

Burkhardt, A. M., Shingledecker, C. N., Le Gal, R., et al. 2019, ApJ, 881, 32

Chauviré, B., Pineau, M., Quirico, E., et al. 2021, Earth Planet. Sci. Lett., 576, 117239

Ciesla, F., & Lauretta, D. 2005, Earth Planet. Sci. Lett., 231, 1

Ciesla, F. J., Lauretta, D. S., Cohen, B. A., et al. 2003, Science, 299, 549

Coleman, G. A. L., Leleu, A., Alibert, Y., et al. 2019, A&A, 631, A7

D'Angelo, M., Cazaux, S., Kamp, I., et al. 2019, A&A, 622, A208

Desch, S. J., Estrada, P. R., Kalyaan, A., et al. 2017, ApJ, 840, 86

Gibb, E. L., Whittet, D. C. B., Boogert, A. C. A., et al. 2004, ApJS, 151, 35

Gillon, M., Jehin, E., Lederer, S. M., et al. 2016, Nature, 533, 221

Gillon, M., Triaud, A. H. M. J., Demory, B.-O., et al. 2017, Nature, 542, 456

Hartmann, L., Calvet, N., Gullbring, E., et al. 1998, ApJ, 495, 385

Houge, A., & Krijt, S. 2023, MNRAS, 521, 5826

Huang, S., & Ormel, C. W. 2022, MNRAS, 511, 3814

Izidoro, A., Bitsch, B., Raymond, S. N., et al. 2021, A&A, 650, A152

Johansen, A., Ronnet, T., Schiller, M., et al. 2023, A&A, 671, A74

Lodders, K. 2021, Space Sci. Rev., 217, 44

Lodders, K., Palme, H., & Gail, H.-P. 2009, Landolt-Börnstein, 4B, 712

Loyd, R. O. P., Shkolnik, E. L., Schneider, A. C., et al. 2018, ApJ, 867, 70

Luger, R., Sestovic, M., Kruse, E., et al. 2017, Nat. Astron., 1, 0129

Lynden-Bell, D., & Pringle, J. E. 1974, MNRAS, 168, 603

Madeira, G., Izidoro, A., & Giuliatti Winter, S. M. 2021, MNRAS, 504, 1854

Metzler, K., Bischoff, A., & Stoeffler, D. 1992, Geochim. Cosmochim. Acta, 56, 2873

Mousis, O., Aguichine, A., Helled, R., et al. 2020, Philos. Trans. R. Soc. London Ser. A, 378, 20200107

Mousis, O., Schneeberger, A., Lunine, J. I., et al. 2023, ApJ, 944, L37

Ogihara, M., Kokubo, E., Nakano, R., et al. 2022, A&A, 658, A184

Ormel, C. W., Liu, B., & Schoonenberg, D. 2017, A&A, 604, A1

Pettke, T., & Bretscher, A. 2022, Earth Sci. Rev., 225, 103896

Piani, L., Marrocchi, Y., Vacher, L. G., et al. 2021, Earth Planet. Sci. Lett., 567, 117008

Pekmezci, G. S., Johnson, T. V., Lunine, J. I., et al. 2019, ApJ, 887, 3

Potapov, A., Bouwman, J., Jäger, C., et al. 2021, Nat. Astron., 5, 78

Reach, W. T., Faied, D., Rho, J., et al. 2009, ApJ, 690, 683

Schneeberger, A., Mousis, O., Aguichine, A., et al. 2023, A&A, 670, A28

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337

Shkolnik, E. L., & Barman, T. S. 2014, AJ, 148, 64

Suttle, M. D., King, A. J., Schofield, P. F., et al. 2021, Geochim. Cosmochim. Acta, 299, 219

Thi, W. F., Hocuk, S., Kamp, I., et al. 2020, A&A, 635, A16

Wagner, W., Riethmann, T., Feistel, R., et al. 2011, J. Phys. Chem. Ref. Data, 40, 043103

Whittet, D. C. B., Boogert, A. C. A., Gerakines, P. A., et al. 1997, ApJ, 490, 729

Youngblood, A., Pineda, J. S., & France, K. 2021, ApJ, 911, 112

Zaikowski, A., Knacke, R. F., & Porco, C. C. 1975, Ap&SS, 35, 97