Shock compression of FeOOH and implications for iron-water interactions in super-earth magma oceans

Domain: earth_science
Published: 2025-12-27
Editorial outcome: major_revision

Abstract
Iron(Fe)-water reactions in a magma ocean can influence water storage and density of planets. These reactions can form Fe-O-H phases, whose density, melting, and electronic properties at planetary interior conditions are important for informing planetary models. Here, we study natural goethite (α-FeOOH) that is shock-compressed along its principal Hugoniot. Analysis of our velocity interferometer system for any reflector (VISAR) results extends the equation of state to over 800 GPa. X-ray diffraction and VISAR reflectivity results indicate the onset of melting occurs at ~95 GPa with complete melting by 166 GPa, which may be relevant to low seismic velocity anomalies observed above the core-mantle boundary. Analysis of X-ray emission spectroscopy results up to 285 GPa shows the spin crossover of Fe, with dominantly low spin Fe above ~265 GPa in the melt, supporting formation of dense basal magma oceans in terrestrial planets. Using our measured FeOOH densities, we model planetary interiors up to 10 Earth masses. Assuming FeOOH forms via iron-water reactions, the radius decreases by up to 28%, while the density increases by up to 165% compared to the unreacted case, providing an avenue to investigate water storage and evolution in super-Earths and sub-Neptunes.

Introduction
The fate of water in super-Earths is critical to understanding their potential habitability. If present in sufficient amount, water in the interiors of super-Earths can also affect their interior structures and mass-radius relationship1. While some sub-Neptunes may accrete tens of percents of water by mass2, most super-Earths have lower water mass fractions. The majority of super-Earths with well-characterized masses and radii are inferred to have less than 3% of water by mass3, especially with Earth-like core mass fractions. Models that combined formation and evolution predict the population of super-Earths to be relatively dry4,5. However, these models considering thermodynamics suggest that water could be common (~1 wt%) in super-Earths6. For example, redox reactions could produce water endogenically from the oxidation of primordial hydrogen by reducing oxidized magma7. Hence, water is likely a common although limited component of super-Earths that formed during the lifetime of the proto-planetary disk. The majority of observed super-Earths are warm and hot worlds, allowing for the presence of magma-oceans even after Gyrs of evolution1,8,9.

The distribution of water within super-Earth interiors influences not only their observed densities but also the evolution of their thermal profile, redox state, and atmospheric composition10,11. Therefore, determining how water is distributed between different depths, and how it chemically interacts with other compounds is essential for interpreting planetary data. These include both their mass-radius data from missions like Kepler, Transiting Exoplanet Survey Satellite (TESS), characterizing ExOPlanet Satellite (CHEOPS), and PLAnetary Transits and Oscillations of stars (PLATO) with radial velocity or transit-timing-variation-follow-ups as well as spectral data from missions like James Webb Space Telescope (JWST), Extremely Large Telescope (ELT), or proposed Habitable Worlds Observatory (HWO) and Large Interferometer For Exoplanets (LIFE) mission.

Conventional super-Earth interior models assume a distinct water layer atop the rocky mantle12. When water partitions into the mantle or metallic core, it can reduce planetary radii at a given mass13. Additionally, the presence of water and the oxidation of iron are closely associated in magma oceans, with implications on the metallic core size and mantle composition14. Given the diverse fates of initial water inventories, it is crucial to investigate how iron and water together affect magma oceans under extreme pressure conditions. Previous high-pressure experiments have shown that metallic iron reacts with water to form FeO and FeHx15–17. This scenario typically occurs in water undersaturated conditions in the present-day Earth’s mantle. In water-rich environments, iron can be further oxidized and hydrogenated to form FeOOH phases18,19. The influence of these reactions on water storage and bulk densities of planets depends on the physical properties of iron hydroxide phases, which remain poorly understood under super-Earth interior conditions.

FeOOH crystallizes in an orthorhombic structure named goethite (α-FeOOH) under ambient conditions. It transforms into ε-FeOOH at ~6 GPa20, and then dehydrogenates to form the pyrite-type FeO2Hx (x=0–1) at approximately 74 GPa and 1600 K21. The melting conditions of FeOOH are key information for the planetary interior model but have only been estimated by first-principles calculations22. Additionally, the predicted equation of state (EOS) of liquid FeOOH is only available up to 140 GPa and 3000–6000 K22, which is insufficient for modeling the interior structure of super-Earths with interior pressures of 0.5–1 TPa. Reaching such extreme pressures is a challenge with current static compression techniques, such as laser-heated diamond-anvil cells; and the long duration of laser heating can locally alter the system’s stoichiometry due to fast hydrogen diffusion23. Shock wave compression presents a promising alternative, because it can not only achieve pressure-temperature conditions relevant to super-Earth interiors, but also likely preserve hydrogen in the sample due to the fast dynamic loading24. However, recent gas gun experiments on FeOOH have been limited to the solid-phase region (below 90 GPa and 2000 K)25, leaving the melting pressure and the densities of liquid FeOOH unknown.

Another important property in constraining the magma ocean behavior is the iron spin state, which has been proposed to affect the iron partitioning between crystallizing solid and residual liquid26,27. Low-spin iron preferentially partitions into the melt at high pressures, potentially leading to the melt being denser than the coexisting solid and affecting magma ocean dynamics26,27. A high to low spin transition can itself reduce the volume of the host phase, resulting in a density increase. Iron in FeOOH undergoes the spin transition at ~50 GPa at 300 K based on previous X-ray emission spectroscopy (XES) and Mössbauer spectroscopy data in diamond anvil cells28,29 and previous theoretical calculations30,31. This spin transition is accompanied by a ~10.9% increase in density29. However, the spin state of iron in liquid FeOOH under high-pressure and high-temperature magma ocean conditions is currently unknown, leaving a critical gap in our understanding of silicate/oxide melts containing iron and water, and consequently their influence on planetary differentiation.

In this study, we performed laser shock compression experiments on natural goethite starting materials. Velocity interferometer system for any reflector (VISAR) measurements at the facility at Laboratoire pour l’utilisation des lasers intenses (LULI) allow us to determine the equation of state on the Hugoniot to over 800 GPa. In-situ X-ray diffraction (XRD) and XES measurements at the Matter in Extreme Conditions (MEC) end-station of the Linear Coherent Linac Source (LCLS) at SLAC National Accelerator Laboratory revealed the melting pressure (~95 GPa) and a spin transition in iron along the Hugoniot. Density data were incorporated into an interior model where iron reacts with water to form FeOOH, providing an avenue for understanding water storage and formation histories of super-Earths that formed under water-poor versus water-rich conditions.

Results
Hugoniot equation of state
The starting sample is natural goethite with the purity of 99.2 wt% (0.8 wt% SiO2 + Al2O3 impurities) and the bulk density, ρ0 = 4.236(81) g/cm3 (see Methods for details; Text S1, Fig. S1, Table S1). The LULI experiments were conducted up to over 800 GPa, while the MEC experiments were conducted below 100 GPa. Step targets were used for both experiments, impedance match methods were used for data analysis, and Monte Carlo simulations were used to evaluate uncertainties. Details of target design, experiments, and data analysis can be found in the Methods and the supporting information (Texts S2, S3; Figs. S2–S4).

Combined LULI and MEC data show that US increases linearly as UP increases (Fig. 1; Table S2). A linear fit to all data yields US = 3.9(3) + 1.6(1) UP. Our data at UP < 3 km/s is consistent with previous gas gun data which show a plateau at US =~2 km/s25 (inset in Fig. 1c). However, we are unable to discriminate between a potential plateau and a linear trend, particularly at UP < 2 km/s; the resulting fit shows negligible variation whether the full dataset is used or restricted to UP > 2 km/s (Fig. S5). Using this relation, pressure was calculated using the equation of P = P0 + ρ0 UP US, where P0 is the ambient pressure (Table S2). We reached the maximum pressure of 819(72) GPa in our experiments, which is relevant to the interior pressure of super-Earths.

To evaluate potential hydrogen preservation during shock compression, we consider the rise time of the drive laser power, t = 10-10 s, in our laser shock experiments. The hydrogen diffusion coefficient, DH, is unknown in the liquid FeOOH under our experimental conditions. Instead, we used DH of 10−3 cm2/s in superionic liquid water at 400–600 GPa and ~5000 K as an upper bound approximation, where hydrogen is unbounded to any other atoms and moves freely across the entire sample23. The diffusion length is thus at most L = √DH t = 10-6.5 cm. Our sample is 10−1 to 100 cm in length and width and 4 to 5 × 10-3 cm in thickness. Therefore, even if superionic behavior is taken into account32, hydrogen is likely to stay in the atomic structure of the sample during the laser shock. Furthermore, the full width at half maximum of the X-ray beam is 3 × 10−3 cm in diameter. Thus, hydrogen is likely to stay within the XRD probing area, although it is difficult to completely rule out the possibility of hydrogen loss because DH in our sample at the relevant pressure-temperature conditions is not constrained.

Melting pressure and temperature
In-situ XRD patterns and VISAR reflectivity measurements were used to evaluate the melting pressure of FeOOH along the Hugoniot. Two criteria were applied to melting. The first is the emergence of diffuse scattering, which may coexist with diffraction peaks, and is indicative of the onset of melting33. The second criterion involves a change in reflectivity at the rear surface of the window, commonly associated with the solid-to-liquid transition of the ejected material34.

To focus on detection of diffuse scattering, which is much weaker than crystalline diffraction peaks, we masked diffraction peaks from solid phases (Fig. 2a and b). At 1 atm, masked diffraction patterns show relatively flat intensity above 15°, compatible with background (Fig. 2b). The signal below 15° results from kapton35. At 95 GPa, the masked pattern starts to show two broad scattering features at 16° and 22° 2θ angles, indicating that a small amount of melt is beginning to form while the shocked FeOOH is still predominantly solid (Fig. 2b). The intensity of the diffuse scattering increases as pressure increases to 138 GPa where the diffraction spots are still observed (Fig. 2b). This observation indicates coexistence of the solid and melt at this pressure. Dehydration has been reported in previous diamond-anvil cell experiments21,36,37, but not observed here. Laser shock compression has a compression rate in the order of 1012 GPa/s. Such ultrafast compression may lead to solid-melt coexistence and suppress dehydration due to kinetic limitations. At 166 GPa, crystalline diffraction peaks were no longer observed, and only diffuse scattering is present (Fig. 2b), indicating the melting completes by 166 GPa (Fig. 2c).

The loss of fringe motion at the sample-vacuum interface serves as additional evidence of melting, likely resulting from the loss of strength upon melting of FeOOH, as suggested by previous studies on shock-compressed SiO234. The VISAR images at 95 GPa (Fig. S6) show a gradual disappearance of fringe motion from the edge to the center. Since the drive laser is focused at the image center, the laser power decreases from the center to the edge. Hence, this change in fringe motion is indicative of onset melting at the center which is consistent with our XRD data (Fig. 2a and b). This onset pressure is lower than that predicted from previous first-principles calculations22 (Fig. 2c). At 138 GPa, the fringe motion vanishes in the whole VISAR images (Fig. S6), indicating that the sample is within the melting region.

The temperature along the Hugoniot can be estimated using the equation of dT = (T/CV) γ dV + (V0 - V) dP + (P - P0) dV, where γ is the Grüneisen parameter and CV is the heat capacity at a constant volume38. γ is assumed to have a form of γ = γ0 (ρ/ρ0)q, where the subscripted ‘0’ denotes the ambient conditions and q is a volume-independent constant. CV can be calculated using the Debye model38 with Debye temperature at ambient conditions, θ0. The γ0, q, and θ0 of goethite were taken from Gleason et al.20 as 0.91, 1, and 740 K, respectively. As temperature approaches to the melting point, the electronic and anharmonic effect on γ and CV may be significant, but experimental constraints are not available. For this reason, in the pressure range of 70–95 GPa we used parameters from theoretical calculations22, specifically CV of 130 J/mol/K and γ of 1.84 for pyrite-type FeO2H. Our calculated temperature for solid FeOOH increases from 300 K to approximately 1855 K at 95 GPa (Fig. 2c), consistent with a previous gas gun study25. Furthermore, previous theoretical calculations suggest that CV remains nearly constant at 139.6 J/mol/K while γ increases with increasing pressure (γ0 = 0.876 and q = −0.914) in liquid FeOOH. Our calculated temperature for liquid FeOOH increases from approximately 3620 K at 166 GPa to approximately 7543 K at 300 GPa (Fig. 2c).

Previous static compression experiments in diamond-anvil cells have suggested that the pyrite-type FeO2Hx phase could account for the anomalously low seismic velocities observed in ultralow-velocity zones above the core-mantle boundary18. Our study reveals that the melt and solid phases may coexist at the core-mantle boundary. The presence of the melt would further enhance the velocity decrease so less FeOOH phase is needed to explain ultralow-velocity zones.

Electronic spin transition
Electronic properties along the Hugoniot were investigated by in-situ XES measurements using an XFEL. Our emission spectra at ambient conditions consist of the main peak Kβ1,3 and the satellite peak Kβ’ (Fig. 3a), consistent with previous FeOOH data in the high spin state29. At high pressures along the Hugoniot, the main peak shifts to a lower energy and the intensity of satellite peak reduces, indicating that a fraction of the iron in the shock-compressed FeOOH transforms to a low-spin state (Figs. 3a and S7). We calculated the integrated absolute difference (IAD) between reference and shot spectra39. As the low spin spectrum remains unknown, we used the high spin emission spectrum at ambient conditions as a reference. Our emission measurements have high accuracy due to experimental geometry, X-ray energy, and probing time (see methods). The precision was evaluated similar to the method in Pardo et al.40 The noise oscillation of our spectra is approximately 3% relative to the height of the main peak. Using a Monte Carlo method, we generated 1000 emission spectra with 3% noise and calculated the IAD for each spectrum. The obtained 1000 IAD values display a normal distribution from which the mean value and the standard deviation were obtained (Fig. S8; Table S3). The obtained IAD values and their uncertainties are plotted in Fig. 3b. The applied 3% variation in the emission spectra results in 3–7% uncertainties in the IAD values, indicating robustness of using IAD values to indicate spin state.

The IAD increases from zero at ambient conditions to ~0.3 at 81 GPa. The emission spectrum at 81 GPa shows the weak satellite peak Kβ’ (Fig. 3), indicating that the fully low spin state has not been reached, i.e., mixed spin state. This observation is different from a previous study in diamond-anvil cells which shows the low spin Fe in goethite above 45 GPa at ambient temperature29. The difference indicates that the high temperature in our shock experiments broadens the spin transition region, as observed in other materials such as magnesiowüstite41.

The IAD values and the main peak position do not change much between 81 to 165 GPa. At this pressure region, melting was identified as discussed above (Fig. 2). As such, the measured emission spectra in this pressure range contain information from both solid and liquid phases. By extending the IAD trend of the solid phase (below 81 GPa) to higher pressures and the liquid phase (above 165 GPa) to lower pressures, we find that the plateau region could be caused by a mixture of signals from the low-spin-dominated solid and the high-spin-dominated liquid phases. That is, the melting may lead to the spin “recovery” from low to high. The loss of structural order upon melting weakens Fe-O interactions which reduces covalency, leading to a diminished crystal field splitting as a consequent “recovery” of the high-spin state.

At pressures above 165 GPa, the IAD increases and the main peak position shifts to lower energies, indicating that the low-spin fraction increases with increasing pressures in the melt (Fig. 3a). Evidence that the iron in the melt may reach the fully low spin state above 265 GPa include (1) the emission spectra showing no observable satellite peak Kβ’, (2) the IAD reaching to 0.5–0.6 that has been linked to low-spin iron in iron compounds42, and (3) the IAD and main peak positions remaining unchanged within experimental uncertainties for the highest pressures (Fig. 3b). We also simulated Fe3+ Kβ XES spectra using the crystal field multiplet theory43,44 in the CTM4XAS program45. The spin transition from the high-spin to low-spin state occurs at a crystal field splitting energy of about 3.0 eV. We used 10Dq = 1.8 eV for the HS state and 10Dq = 3.6 eV for the LS state. The theoretical Fe3+ XES show well agreement with the experimental observations. Upon the high-spin to low-spin transition, the main Kβ1,3 line exhibits a negative energy shift, while the Kβ’ feature shifts to higher energy (Fig. 3c). The reduced Kβ’-Kβ1,3 splitting reflects an increase of covalency, associated with the increase of crystal field. In addition, the ratio of Kβ’/Kβ1,3 reflects S/(S+1). All these changes originate from the strong 3p-3d exchange interaction, which plays a dominant role in the Kβ emission channel.

These results indicate that iron- and hydrogen-rich melts, which are likely to form during the late stages of fractional crystallization, can exhibit varying proportions of low-spin iron depending on depth and the size of the super-Earth46. In larger super-Earths with deeper magma oceans, a greater proportion of low-spin iron may be stabilized. If low-spin iron preferentially partitions into the melt26,27, it could become denser than the coexisting solid, potentially leading to the formation of a basal magma ocean in a super-Earth interior.

Density of FeOOH along the Hugoniot
Densities along the Hugoniot can be calculated using the equation of ρ = ρ0 / (1 – UP / US) (Fig. S9; Table S2)38. This equation is derived from the mass conservation which is universal regardless of phase transitions. Thus, we used ρ0 of our starting sample for all different phases/structures along the Hugoniot. The density increases significantly from 4.236 g/cm3 at ambient conditions to approximately 6.8 g/cm3 at 196 GPa, and then slowly increases to approximately 8.3 g/cm3 at 819 GPa (Fig. S9). The different slope is due to the pressure effect on compressibility and/or different solid versus liquid compressibility. Using the EOSfit7 software47, we fit the P-ρ data of FeOOH melt data at pressure above 166 GPa to the second-order Birch-Murnaghan EOS. We used a reference pressure of 196 GPa where the density is 6.8(4) g/cm3 from our LULI data (Table S2). The bulk modulus at this reference pressure was obtained to be 1407(423) GPa. The density information is further used in the planetary interior model in the next section.

Planetary interior model
We investigate the implications of the formation of FeOOH on interior models of planets. Specifically, we are interested in the effect on planet bulk density, which we will quantify in the following, and on the potential gas speciation from outgassing, which we address speculatively in the discussion.

For investigating the effect on planet bulk density, we employ the interior model from Dorn et al.48 with recent updates in Luo et al.13 (refer to Methods for details). For the choice of the interior model, we are guided by the chemical reactions based on previous mineral physics experiments15,16,18,19:

3Fe + H2O → FeO + 2FeH (P < 86 GPa) (1)

4Fe + 2H2O → FeOOH + 3FeH (P > 86 GPa) (2)

Our reference case (case A) is according to commonly used interior models in exoplanetary science that assume pure iron cores and iron-free silicate mantles with water only added to the surface (case A in Fig. 4, and green curves in Fig. 5a). The largest effect of reactions (1) and (2) can be studied focusing on a hypothetical end-member case (case C in Fig. 4), where all the iron in the mantle takes part in reaction (1) and (2) depending on the pressure range. If all the iron is reacting with water assuming an Earth-like interior (i.e., a core mass fraction fcore of 0.325), the amount of water is constrained to values below 5 wt%, depending on the planet mass. The upper limit of 1.3 wt% can be understood from considering only reaction (2) and is derived from the molar ratio of reactants and the molar mass ratios of water (MH2O) and iron (MFe). Specifically, since only one out of four Fe atoms is oxidized to FeOOH in reaction (2), we assume that 25% of the total Fe participates in the reaction. Accordingly, the water mass fraction fH2O equals to 0.25 * fcore * (MH2O / MFe) * (2/4), where MH2O and MFe are equal to 18.015 and 55.845 g/mol, respectively. When we allow for the formation of FeO and FeOOH (case C), we keep the same bulk composition as in case A, but all water is now chemically bound to FeO and FeOOH in the mantle15–17 and FeH in the core. At higher planet masses, the pressure gradient in the interior is steeper and hence the mass fraction of the planet experiencing pressures above 86 GPa, where FeOOH is stable, becomes larger. The formation of 1 mole FeOOH requires 1 mole more water more than the formation of 1 mole FeO. In consequence, the amount of water required to fully react iron is 1.3 wt% for a 1 M⊕ (Earth mass) planet. Here, we depict the molten silicate mantle with a 2-layer structure, according to the stability regions of FeO and FeOOH. We simply represent the mantle to be a mixture of silicates with FeO and FeOOH, depending on pressure.

The main effect of accounting for reactions (1) and (2) in case C compared to commonly used interior models (case A) is a reduction by ~6-28% in radii which implies a density increase of ~22–165% (Fig. 5a). Compared to case A, the core in case C shrinks by about 1/3–1/4 and the mantle expands as O and H are added to the mantle (Fig. 4). This combined effect causes the radius at the top of the mantle in case C to increase by 1.5–3%. Thus, the total radii reduction in case C is dominated by the lack of any steam water layer at the surface as water stores in the interior of the planets following reactions (1) and (2).

Also, our hypothetical end-member case C yields smaller radii compared to case B when water is allowed to be dissolved as O and H in the (core and mantle) melts according to solubilities and partitioning coefficients as done in Luo et al.13 (Fig. 5). This is expected as case B includes a steam water layer by definition, which significantly increases planet radii compared to bare rocky worlds.

Our mass-radius curves span a region from super-Earths to sub-Neptunes. Here, we now focus on the end-member scenario where we allow for the maximum amount of water that is needed to react all iron according to reactions (1) and (2), which varies from 0.9 to 1.3% across the planet mass range of interest. Interestingly, this value is lower than the 3% of water that has been inferred to be an upper bound on the possible amount of water for the majority of super-Earths3 given observational data. In addition, most observed super-Earths are hot worlds that likely underwent significant atmospheric loss due to their proximity to the host star49. This implies that the amount of water available for chemical reactions early in the history of the planet may not be present at the time of observation. With case C we highlight a case when no water is present at the surface; however, not all available water may have reacted with iron early in a planet’s history and some water can be dissolved in silicates or remain present at the surface. Our hypothetical end-member scenario highlights that the treatment of water in interior models has profound implications for the interpretation of observational results. Treating interactions between water and iron (case C, red) leads to radii that are matched in the super-Earth population, however, the identical bulk composition with a commonly layered model (case A, green) leads to a mass-radius curve that lies in a transitional regime between super-Earths and sub-Neptunes where only few planets are observed. However, we note that for realistic super-Earths, not all iron may have reacted with water as the reactants may be partially separated in growing super-Earths interiors and likely not all iron participates in reactions (1) and (2) or the amount of water is simply limited well below 3%. Furthermore, we do not consider any silicates to be part of the reaction system as the experiments are restricted to iron and water only. Therefore, case C is a hypothetical end-member case.

Where does the water come from? Although water can be added primordially when building blocks from outside the water snow line are accreted50, water can also be formed endogenously by redox reactions between FeO and any primordial hydrogen envelope, which has been estimated to lead to water amounts in the range of 1%, so roughly similar to our end-member case C. If only 1 wt% of water is available, almost all of the available iron in Earth-like interiors can react to form FeOOH, ignoring any interaction of silicates. In principle, O and H can also partition to the metal phase during core formation, which we have not taken into account here. Chemical thermodynamic approaches are needed to take this consistently into account.

The exoplanet community is eager to learn which super-Earths evolved from the sub-Neptune population via atmospheric loss and which super-Earths were born rocky51. Differences in bulk composition and their effect on bulk densities are significant (Fig. 5). The chemical imprint of hydrogen envelopes or the lack thereof on the chemistry of super-Earth atmospheres is another promising avenue for distinguishing the planetary evolution histories6 thanks to missions like JWST, ELT, HWO, and LIFE.

The imprint of FeOOH formation on the gas speciation of any outgassing secondary atmosphere is therefore of broad interest. FeOOH-rich mantles likely outgas more oxidized species compared to mantles containing FeO. Clearly, the higher oxidation power of FeOOH compared to FeO provides the possibility to distinguish dry super-Earths from water-rich worlds.

Methods
Sample characterization and target design
Our starting samples are natural goethite from the Restormel Royal Iron Mine, Cornwall U.K. They exhibit an acicular or prismatic crystal habit and a size of a few millimeters in width and height and approximately a centimeter in length (Fig. S1). Single-crystal XRD analysis52 reveals that the sample has orthorhombic symmetry (space group Pnma) with lattice parameters of a = 9.9506(6) Å, b = 3.0190(2) Å, c = 4.6001(3) Å, V = 138.191(15) Å3, and ρ0 = 4.2707(5) g/cm3, consistent with the goethite structural parameters at ambient conditions (Text S1; Table S1)53. Density measurements using the Archimedes method give a density of 4.236(81) g/cm3, based on the mean and standard deviation calculated from seven analyses. The consistent results for densities derived from XRD and the Archimedes method indicate negligible porosity within measurement uncertainties. Energy-dispersive X-ray spectroscopy (EDS) shows a homogeneous FeOOH composition with approximately 0.8 wt.% SiO2 + Al2O3 impurities (Text S1; Fig. S1).

The samples were double-sided polished to a thickness of 35–45 µm using 3M diamond films, and the thickness was measured using a micrometer. To enhance the reflectivity in the VISAR measurements, a 300 nm thick Al coating was deposited on one side of each sample. Step target designs were used in both LULI and MEC experiments. The LULI step target consists of a 10 or 50 µm thick kapton ablator, a 40 µm thick Al foil serving as a pusher to create a uniform shock wave, a 44–200 µm thick SiO2 or LiF window material attached to one half of the Al surface, and the sample on the other half. A small gap of a few µm between the sample and the window allows for determining the breakout time of the kapton (Fig. S2a). In the MEC experiments, the Al-coated sample platelets were glued to one half of the Al-coated kapton surface, with the kapton having a ~80 µm thickness (Fig. S2b).

Laser-driven shock experiments
The experiments were performed at the LULI 2000 high-energy laser facility, utilizing two laser beams at a wavelength of 527 nm with either 2 ns or 5 ns square pulses. Each laser pulse has a single-pulse energy of up to 500 J. Phase plate smoothing was employed to produce an 800 µm diameter focal spot, achieving a maximum intensity of 1014 W/cm2 on the target when both laser beams were combined. Diagnostics such as VISAR were used to measure the reflectivity and shock velocity along the Hugoniot at wavelengths of 1064 nm (ω) and 532 nm (2ω).

Laser-driven shock compression experiments have been performed using the frequency-doubled Nd-glass laser beams available at the MEC end-station of LCLS at SLAC National Accelerator Laboratory. We produced Hugoniot states inside the sample by using flat-top temporal pulses of 8 ns with focal spots of 300 µm diameter on target using continuous phase plates available at MEC. The hydrodynamics was monitored using two velocity interferometry systems for any reflector (VISAR) operating at 532 nm.

VISAR data analysis
VISAR images in the LULI experiments show fringe motion in the LiF or SiO2 window and the transit time of the shock wave inside the sample. Analysis of the fringe motion using the Neutrino software, together with the Hugoniot of LiF and SiO254,55, reveals UP in LiF or UP in SiO2. Note that the shock front in SiO2 becomes reflective at pressures above approximately 100 GPa along the Hugoniot56. The transit time and sample thickness were used to calculate US in the sample. US and P of the sample can be further constrained using an impedance matching technique from the release isentrope and the reshock Hugoniot of Al38 (Text S2). The mean and standard deviation of US and P were obtained by running 300 Monte Carlo simulations with consideration of uncertainties of Hugoniot parameters of the window and Al, the fringe motion analysis, the sample initial density, and the transit time in the sample (Text S2).

VISAR images from the MEC experiments provide the transit time of the shock wave inside the kapton. Using a previous Hugoniot of kapton and the Hugoniot of goethite FeOOH from our LULI experiments, the sample’s US and P were constrained by employing an impedance matching method38 (Text S3). Also, 1000 Monte Carlo simulations were performed to derive the frequency distribution of US and P in sample with uncertainty sources of the transit time in kapton and Hugoniot parameters of both kapton and goethite FeOOH. The mean and standard deviation of US and P were finally calculated from the frequency distribution.

XES and XRD
In-situ Kβ XES spectra or XRD patterns were collected during laser shock compression at LCLS, MEC, following the previous experimental setups33,46. X-ray pulses are approximately 50 fs wide and are quasi-monochromatic with only 0.2–0.5% energy variation. X-ray beam focuses down to 30 µm in diameter in the center of the uniformly shocked area on the sample. X-ray probing time is 0.5–1 ns prior to the sample breakout so that both emission and diffraction probe uniformly shocked samples.

X-ray energy was set to 8 keV for XES measurements. We used a multi-crystal energy dispersive spectrometer where a crystal analyzer and a position-sensitive detector (PSD) follow the von Hamos geometry57. The analyzer was placed at the 90° angle horizontally to the incident X-ray to avoid elastic scattering. The analyzer consists of arrays of bent crystals. The curvature of these crystals is used to focus X-rays emitted from the sample onto the PSD. The direction perpendicular to this curvature is responsible for dispersing X-ray energies across the PSD. By integrating the PSD signal along the focusing at different energies, an emission spectrum is obtained. This geometry allows us to collect most of the emission signal from the first 10 µm of compressed sample46.

X-ray energy was set to 17 keV for XRD measurements. The diffraction data were collected using four ePix10k detectors which were placed in the standard detector configuration at MEC, LULI33. We used LaB6 to calibrate the tilting and the distance to the sample of four detectors in the Dioptas software58. Further data analysis was performed in the same software.

XES spectra simulation
We simulate the Fe3+ Kβ XES spectra using the crystal field multiplet theory43,44 implemented in the CTM4XAS program45. Fluorescence emission was considered part of the x-ray inelastic scattering and can be described by the Kramers-Heisenberg equation59. Kβ XES multiplets are fully considered, including 3d-3d, 1s-3d, 3p-3d interactions, 3d spin-orbit coupling (SOC), and crystal field splitting. The Slater integrals Fdd and Fpd were scaled to 80% of their Hartree-Fock values, while Gpd was scaled to 56%. Temperature effects use the Boltzmann weighting factor. The core-hole lifetime broadening and instrumental resolution were modeled using a Lorentzian function with a half-width at half-maximum (HWHM) of 1.5 eV and a Gaussian function with a HWHM of 1.2 eV, respectively. Due to the symmetry distortion in FeO(OH), a distorted D4h symmetry was considered in the XES calculations, where Ds and Dt considered as 0.1 eV and 0 eV, respectively. These specific numerical values are not critical for this study. The simulated high-spin spectrum was energy-shifted to align with the experimental spectrum at 1 atm. The Fe3+ ion (3d5) has a high-spin ground state corresponding to the 6A1 term. In the low-spin configuration, the ground state becomes the 2T2g term. A crystal field splitting energy (10 Dq) of approximately 3.0 eV is required to induce the spin transition from the high-spin to the low-spin state.

Interior model
The interior model from Dorn et al.48 with recent updates in Dorn & Lichtenberg1 and Luo et al.13 solves hydrostatic equilibrium, mass conservation, and equations of state. Here, we consider an iron dominated core, a silicate dominated mantle and a water (steam) layer. For the core, we consider Fe, H, and O, if allowed. For solid Fe, we use the equations of state for hexagonal close packed (hcp) iron from literature60,61. For liquid iron with, and without H and O, we use Luo et al.13. The core thermal profile is assumed to be adiabatic throughout the core. At the core-mantle boundary (CMB), there is a temperature jump as the core can be hotter than the mantle due to the residual heat released during core formation. We follow Stixrude et al.62 and add a temperature jump at the CMB such that the temperature at the top of the core is at least as high as the melting temperature of the silicates.

The mantle is assumed to be made up of three major constituents, i.e., MgO, SiO2, and FeO/FeOOH. For the solid mantle, we use the thermodynamical model Perple_X63, to compute stable mineralogy and density for a given composition, pressure, and temperature, employing the database of Stixrude et al.64. For pressures higher than 125 GPa, we define stable minerals a priori and use their respective equation of states from various sources65–69. For the liquid mantle, which is the focus of this work, we calculate its density assuming an ideal mixture of main components (Mg2SiO4, SiO2, and FeO)66,70–72 and add them using the additive volume law. Note that we use Mg2SiO4 instead of MgO since the data for forsterite is available in a high-pressure temperature regime, which is not available for MgO to our knowledge. The temperature gradient is assumed to be adiabatic, and for multi-component mixtures, we use the harmonic mean of the single-component adiabatic gradients.

Water can be added to the mantle melts, while the solid mantle is assumed to be dry. The addition of water reduces the density, for which we follow Bajgain and coauthors73 and decrease the melt density per wt% water by 0.036 g/cm3. For small water mass fractions, this reduction is nearly independent of pressure and temperature.

The melting curve of mantle material is calculated for dry and pure MgSiO3 to which the addition of water74 and iron75 can lower the melting temperatures. The presence of H and O within the core will also lower its melting temperature, for which we follow Luo et al.13. The partitioning between mantle melts and the water layer is determined by a modified Henry’s law, for which we use the fitted solubility function of Dorn & Lichtenberg1. For the partitioning of water between iron and silicates, we follow Luo et al.13. For the equilibration pressure of water to partition between iron and silicates, we use half of the core-mantle boundary pressure, which is roughly within typically discussed values for Earth (0.3–0.6). Varying this pressure would introduce overall small changes in the distribution of water.

Above the mantle, a water layer can be present. For most scenarios that we investigate here, water at the surface is in steam phase. We employ the EOS compilation of different water phases from Haldemann et al.76. We further assume an adiabatic temperature profile in the steam phase except for pressures below 0.1 bar, where we assume an isothermal profile following the equilibrium temperature. The radius of the planet is defined at 1 mbar. For our interior model, we use EOS for Fe13,60,61 in the core instead of FeH. The effect of H on the density of FeH is minor because hydrogen atoms are much smaller than iron atoms, and thus the volume change from Fe to FeH is minimal. The effect of using Fe versus FeH in the core on the radius is less than ~0.5%, which we quantify by density differences between the employed Fe EOS and the FeH EOS77.

Peer review rounds
[{"round": 1, "reviewer_comments": {"reviewer_1": "The authors conducted a detailed investigation of goethite (α-FeOOH) under extreme pressure–temperature (P–T) conditions relevant to the interiors of super-Earths, utilizing state-of-the-art laser shock compression techniques. This comprehensive study explored key physical properties of FeOOH—including density, melting behavior, and iron spin crossover—across an unprecedented pressure range reaching up to 800 GPa. The results offer novel insights into iron–water interactions in terrestrial planetary environments and carry important implications for magma ocean dynamics and the chemical evolution of exoplanetary interiors.\nWhile the work represents a notable technical advancement in dynamic high-pressure experimentation and provides valuable new constraints for planetary modeling, several critical issues must be addressed before the manuscript can be considered for publication.\nComments and Suggestions:\n1. The Gan et al (PRB 2023)’s study delivered very detailed shock compression data up to ~90 GPa, including both experimental and theoretical results. In comparison, the current work significantly extends the pressure range to 800 GPa and ventures into the melted regime of FeOOH, which is a substantial advancement. However, the dataset above ~100 GPa appears considerably scattered. While this may reflect the inherent experimental difficulty of collecting data under such extreme P–T conditions, the conclusions drawn from these data are not sufficiently substantiated.\nThe melting pressure determinations also present inconsistencies: XRD and VISAR suggest differing melt onset pressures, and the interpretation of XES-IAD results remains speculative. Moreover, the reported observation of a low-spin FeOOH melt phase above ~260 GPa is a significant result for understanding Fe-rich melts in planetary interiors. However, the density data in Figure S8 show minimal change across the spin transition, which deviates from the typical expectation of a pronounced density increase during spin crossover. In this context, density functional theory (DFT) simulations could offer valuable support to validate these experimental observations.\n- In Figure 2C, it is essential to clearly distinguish between experimental and theoretical data. For instance, Deng’s dataset originates from DFT predictions and should be explicitly labeled as such.\n- Spin crossover in solid FeO₂H is computationally predicted to occur around 60 GPa, while their results indicate a significantly delayed transition in the melt. Authors speculated that it is due to a potential breakdown in crystal field splitting within the disordered liquid structure. While plausible, this interpretation would benefit from stronger theoretical backing—either through new simulations or by referencing analogous studies in the literature.\n2. The authors observe a broad pressure range for solid-melt coexistence along the Hugoniot. Could the authors elaborate on the factors contributing to this extended coexistence region? Specifically, given that hydrous minerals typically undergo dehydration at elevated temperatures, is there evidence that dehydration reactions (e.g., FeOOH → Fe₂O₃ + H₂O) may be occurring in these experiments? Could the observed solid-melt coexistence be partially attributed to such dehydration processes, where the \"melt\" phase might actually represent a mixture of FeOOH/Fe2O3 and liberated water?\n3. The Hugoniot path appears to intersect potential dehydration boundaries for FeOOH (CaCl2 or pyrite structures). it would be valuable to include experimental or calculated dehydration curves for FeOOH at relevant pressures. In addition, the temperature values along the axis are currently difficult to distinguish. Please consider Using a linear scale if appropriate for the temperature range.", "reviewer_2": "This is a strong, interdisciplinary paper which was a pleasure to review. The authors present novel high-pressure data on a mineral (FeOOH) which was previously unconsidered in planet interior models, but which was shown here to significantly impact the fate of water in hot planets. The fate of water in the exoplanet population is of significant interest to the exoplanet community and the resulting decrease in planetary radii from a steam atmosphere but decrease in density in the rock portion of the planet would change the interpretation of many planets. It presents clear gaps in our knowledge of planet interiors for future interdisciplinary work. While I cannot comment in depth on the high-pressure experiments which the work hinges on, I have a few minor concerns that should be clarified about the interior model.\n1.) Density of the material in the interior model is treated in several different ways between ideal mixing and new measurements. I would be interested in what the average density in the mantle and core is in the three cases for a selected planet. It would be most helpful for comparisons if the authors include a plot of density with radius/mass.\n2.) Line 513 does some heavy lifting on how the temperature gradient is chosen. How is the adiabat chosen for multi-component mixtures?\n3.) Using Fe in the core rather than FeH is explained in line 533 as having a small effect. What EOS is being used for FeH here? Can the authors further clarify why including such a light element in the core would have small effects?\n4.) In Line 498, the authors refer to using liquid Iron from Luo et al., but I don’t believe that work presented new equations of state for liquids. The authors probably mean they use the same prescription as in Luo, but it reads as if new measurements were used for the EoS from Luo.\nMinor/Stylistic comments:\n5.) The introduction presents some exoplanet demographics as decided facts versus debated astrophysics. The second sentence has grammatical issues. The 3% water in Line 48 should be noted to be for planets with Earth-like core mass fractions. Lastly, Paragraph 1 and Paragraph 3 feel contradictory at points in that P1 discusses how planets are relatively dry while P3 discusses diverse water inventories.\n6.) The text in Figure 5 is very small.\n7.) In Fig. 5 caption Line 328 there’s a quick smile-inducing typo."}, "author_response": "Point-by-point response to the reviewers’ comments:\nReviewers’ comments are presented in black bold font, while our point-by-point responses are\nin blue. Revisions to the manuscript are highlighted in yellow.\nReviewer #1 (Remarks to the Author):\nThe authors conducted a detailed investigation of goethite (α-FeOOH) under extreme\npressure–temperature (P–T) conditions relevant to the interiors of super-Earths, utilizing\nstate-of-the-art laser shock compression techniques. This comprehensive study explored\nkey physical properties of FeOOH—including density, melting behavior, and iron spin\ncrossover—across an unprecedented pressure range reaching up to 800 GPa. The results\noffer novel insights into iron–water interactions in terrestrial planetary environments and\ncarry important implications for magma ocean dynamics and the chemical evolution of\nexoplanetary interiors.\nWhile the work represents a notable technical advancement in dynamic high-pressure\nexperimentation and provides valuable new constraints for planetary modeling, several\ncritical issues must be addressed before the manuscript can be considered for publication.\n[Author] We thank the reviewer for recognizing the novelty and values of our work, and for\nconstructive comments that improve our manuscript. We have carefully considered and\naddressed these comments below.\nComments and Suggestions:\n1. The Gan et al (PRB 2023)’s study delivered very detailed shock compression data up to\n~90 GPa, including both experimental and theoretical results. In comparison, the current\nwork significantly extends the pressure range to 800 GPa and ventures into the melted\nregime of FeOOH, which is a substantial advancement. However, the dataset above ~100\nGPa appears considerably scattered. While this may reflect the inherent experimental\ndifficulty of collecting data under such extreme P–T conditions, the conclusions drawn\nfrom these data are not sufficiently substantiated.\n[Author] We appreciate the reviewer’s acknowledgment of the technical challenges in extending\nmeasurements to such extreme pressure-temperature conditions. We used different experimental\ntechniques with Gan et al.1 We used laser-driven shock compression which has a micron-second\ntimescale and an extended pressure range compared to the gas gun used in Gan et al.1 We have\nrevised the manuscript to make our conclusions more substantiated according to the specific\ncomments below.\ndiffering melt onset pressures, and the interpretation of XES-IAD results remains\nspeculative. Moreover, the reported observation of a low-spin FeOOH melt phase above\n~260 GPa is a significant result for understanding Fe-rich melts in planetary interiors.\n[Author] To better discuss melting processes under laser-driven shock compression, we have\nadded raw VISAR 1 and 2 images at 95 and 138 GPa in Figure S6. The drive laser irradiates the\ncenter of the image, where the laser power is the highest and gradually decreases towards the\nedge. At 95 GPa, the VISAR images show fringe motion on the left side that progressively\ndisappears towards the center. This indicates the onset of melting at the center, consistent with\nthe XRD measurements, which probe the same region. At 138 GPa, no fringe motion is observed\nfrom the center to the left side, suggesting that the sample is within the melting regime. We have\nmodified statements in the main text and Figure S6.\nLines 175-183: The loss of fringe motion at the sample-vacuum interface serves as additional\nevidence of melting, likely resulting from the loss of strength upon melting of FeOOH, as\nsuggested by previous studies on shock-compressed SiO 2. The VISAR images at 95 GPa (Fig.\nlaser is focused at the image center, the laser power decreases from the center to the edge. Hence,\nthis change in fringe motion is indicative of onset melting at the center which is consistent with\nour XRD data (Figs. 2a and b). This onset pressure is slightly lower than that predicted from\nprevious first-principles calculations3,4 (Fig. 2c). At 138 GPa, the fringe motion vanishes in the\nwhole VISAR images (Fig. S6), indicating that the sample is within the melting region.\nFigure S6:\nof X-ray beam. In (a), the fringe motion after the breakout gradually vanishes from the left to the\ncenter, indicating the center starts to melt at 95 GPa. In (b), the fringe motion vanishes in the\nentire images, indicating the whole sample is in the molten state.\nThe comment on XES-IAD results is addressed in our response to a more specific comment\nbelow.\nHowever, the density data in Figure S8 show minimal change across the spin transition,\nwhich deviates from the typical expectation of a pronounced density increase during spin\ncrossover. In this context, density functional theory (DFT) simulations could offer valuable\nsupport to validate these experimental observations.\n[Author] Our IAD results reveal a broad spin transition from approximately 50 to 265 GPa,\nlikely due to the temperature effect during the laser-driven shock compression. High temperature\ncan smooth out the expected density discontinuity (Trautner et al.5). Moreover, pressure and\ntemperature increase simultaneously along the Hugoniot, which is different from diamond-anvil\ncell experiments where the density jump is observed at a constant temperature. Additionally, the\nmelting could counterbalance the density increase upon spin crossover. Therefore, we do not\nexpect to see a pronounced density increase during spin crossover. This is also the case in\nprevious Hugoniot data of FeO and Fe O (Fig. S9). Conducting DFT calculations would indeed\n2 3\nbe valuable. However, accurately accounting for spin transitions introduces considerable\n3\ncomplexity and challenges. We also believe that such calculations would not substantially\nimprove the present study. Therefore, DFT calculations along the Hugoniot are beyond the scope\nof this work and could be pursued in future studies.\n- In Figure 2C, it is essential to clearly distinguish between experimental and theoretical\ndata. For instance, Deng’s dataset originates from DFT predictions and should be explicitly\nlabeled as such.\n[Author] We have labeled “DFT” for Deng’s dataset and “Exp” for experimental work in the\nfigure legend and noted them in the figure caption.\nFigure 2. … Exp, experiment; DFT, density functional theory; α, goethite; ε, CaCl -type\n- Spin crossover in solid FeO₂H is computationally predicted to occur around 60 GPa, while\ntheir results indicate a significantly delayed transition in the melt. Authors speculated that\nit is due to a potential breakdown in crystal field splitting within the disordered liquid\nstructure. While plausible, this interpretation would benefit from stronger theoretical\nbacking—either through new simulations or by referencing analogous studies in the\nliterature.\n[Author] We simulated Fe3+ Kβ X-ray emission spectra using the crystal field multiplet theory\nimplemented in the CTM4XAS program6–10. The spin transition from the high-spin (HS) to low-\nspin (LS) state occurs at a crystal field splitting energy of 10 Dq = 3.0 eV. In our calculations, we\nused 10 Dq = 1.8 eV for the HS state and 10 Dq = 3.6 eV for the LS state. The theoretical Fe³⁺\nthe main Kβ₁,₃ line exhibits a negative energy shift, while the Kβ′ feature shifts to higher energy.\nThe reduced Kβ′–Kβ₁,₃ splitting reflects an increase covalency, associated with the enhanced\ncrystal field strength. We have added panel (c) in Figure 3 and some discussion in the main text.\nFigure 3. ... (c) Simulated Fe3+ Kβ XES spectra with high-spin (HS, 10 Dq = 1.8 eV) and low-\nspin (LS, 10 Dq = 3.6 eV) states. The vertical dashed line denotes the main peak position of the\nHS spectrum.\nDiscussion at lines 257-274: That is, the melting may lead to the spin “recovery” from low to\nhigh. The loss of structural order upon melting weakens Fe-O interactions which reduces\ncovalency, leading to a diminished crystal field splitting as a consequent “recovery” of the high-\nspin state.\nAt pressures .... We also simulated Fe3+ Kβ XES spectra using the crystal field multiplet\ntheory7,8 in the CTM4XAS program9. The spin transition from the high-spin to low-spin state\noccurs at a crystal field splitting energy of about 3.0 eV. We used 10 Dq = 1.8 eV for the HS state\nand 10 Dq = 3.6 eV for the LS state. The theoretical Fe3+ XES show well agreement with the\nexperimental observations. Upon the high-spin to low-spin transition, the main Kβ line\n1,3\nexhibits a negative energy shift, while the Kβ’ feature shifts to higher energy (Fig. 3c). The\nreduced Kβ’-Kβ splitting reflects an increase of covalency, associated with the increase of\n1,3\ncrystal field. In addition, the ratio of Kβ’/Kβ reflects S/(S + 1). All these changes originate\n1,3\nfrom the strong 3p-3d exchange interaction, which plays a dominant role in the Kβ emission\nchannel.\nMethodology at lines 480-496:\nXES spectra simulation\nthe CTM4XAS program9. Fluorescence emission was considered part of the x-ray inelastic\nscattering and can be described by the Kramers-Heisenberg equation10. Kβ XES mutiplets are\nfully considered, including 3d-3d, 1s-3d, 3p-3d interactions, 3d spin-orbit coupling (SOC), and\ncrystal field splitting. The Slater integrals Fdd and Fpd were scaled to 80% of their Hartree-Fock\nvalues, while Gpd was scaled to 56%. Temperature effects use the Boltzmann weighting factor.\nThe core-hole lifetime broadening and instrumental resolution were modeled using a Lorentzian\nfunction with a half-width at half-maximum (HWHM) of 1.5 eV and a Gaussian function with a\nHWHM of 1.2 eV, respectively. Due to the symmetry distortion in FeO(OH), a distorted D\n4h\nsymmetry was considered in the XES calculations, where Ds and Dt considered as 0.1 eV and 0\neV, respectively. These specific numerical values are not critical for this study. The simulated\nhigh-spin spectrum was energy-shifted to align with the experimental spectrum at 1 atm. The\nFe3+ ion (3d5) has a high-spin ground state corresponding to the 6A term. In the low-spin\n1\nconfiguration, the ground state becomes the 2T term. A crystal field splitting energy (10 Dq) of\n2g\napproximately 3.0 eV is required to induce the spin transition from the high-spin to the low-spin\nstate.\n2. The authors observe a broad pressure range for solid-melt coexistence along the\nHugoniot. Could the authors elaborate on the factors contributing to this extended\ncoexistence region? Specifically, given that hydrous minerals typically undergo dehydration\nat elevated temperatures, is there evidence that dehydration reactions (e.g., FeOOH →\nFe₂O₃ + H₂O) may be occurring in these experiments? Could the observed solid-melt\ncoexistence be partially attributed to such dehydration processes, where the \"melt\" phase\nmight actually represent a mixture of FeOOH/Fe2O3 and liberated water?\n[Author] We did not observe diffraction peaks of Fe O or other dehydrated phases (such as\n2 3\npyrite phase) from our XRD patterns. We observed hot and compressed goethite coexisting with\nmelt, and thus, the melt is likely composed of FeOOH (Fig. 2). This may reflect kinetics of\ndehydration and melting. Laser-driven shock compression has a compression rate of more than\ncompression may cause the solid-melt coexistence and prevent dehydration from occurring. We\nhave added a few sentences to elaborate this point at lines 170-173: Dehydration has been\nreported in previous diamond-anvil cell experiments11–13, but not observed here. Laser shock\ncompression has a compression rate in the order of 1012 GPa. Such ultrafast compression may\nlead to solid-melt coexistence and suppress dehydration due to kinetic limitations.\n3. The Hugoniot path appears to intersect potential dehydration boundaries for FeOOH\n(CaCl2 or pyrite structures). it would be valuable to include experimental or calculated\n6\ndehydration curves for FeOOH at relevant pressures. In addition, the temperature values\nalong the axis are currently difficult to distinguish. Please consider Using a linear scale if\nappropriate for the temperature range.\n[Author] We have added the CaCl -type FeOOH (ε-FeOOH) dehydration curve from Gan et al.14\n2\nand Koemets et al.15, and the pyrite phase dehydration curve from Nishi et al.12 Now we use a\nlinear scale for the temperature axis. Please refer to our response to a previous comment for the\nupdated figure.\nReviewer #2 (Remarks to the Author):\nThis is a strong, interdisciplinary paper which was a pleasure to review. The authors\npresent novel high-pressure data on a mineral (FeOOH) which was previously\nunconsidered in planet interior models, but which was shown here to significantly impact\nthe fate of water in hot planets. The fate of water in the exoplanet population is of\nsignificant interest to the exoplanet community and the resulting decrease in planetary\nradii from a steam atmosphere but decrease in density in the rock portion of the planet\nwould change the interpretation of many planets. It presents clear gaps in our knowledge of\nplanet interiors for future interdisciplinary work. While I cannot comment in depth on the\nhigh-pressure experiments which the work hinges on, I have a few minor concerns that\nshould be clarified about the interior model.\n[Author] We thank the reviewer for recognizing the significance and novelty of our work. We\nhave addressed those minor concerns below.\n1.) Density of the material in the interior model is treated in several different ways between\nideal mixing and new measurements. I would be interested in what the average density in\nthe mantle and core is in the three cases for a selected planet. It would be most helpful for\ncomparisons if the authors include a plot of density with radius/mass.\n[Author] We have added Fig. S10 to show a plot of density with radius for three cases.\n2.) Line 513 does some heavy lifting on how the temperature gradient is chosen. How is the\nadiabat chosen for multi-component mixtures?\n[Author] The temperature gradient is assumed to be adiabatic. For multi-component mixtures we\nuse the harmonic mean of the single component adiabatic gradients. We have clarified this point\nat lines 517-519: The temperature gradient is assumed to be adiabatic, and for multi-component\nmixtures, we use the harmonic mean of the single-component adiabatic gradients.\n3.) Using Fe in the core rather than FeH is explained in line 533 as having a small effect.\nWhat EOS is being used for FeH here? Can the authors further clarify why including such\na light element in the core would have small effects?\n[Author] The effect of H on the density of FeH is minor, because H is very small compared to Fe\nand hence the volume change from Fe to FeH is minor. Comparision between EOS of FeH\n(Tagawa et al. 2022)16 and Fe (Hakim et al. 2018, Miozzi et al. 2020, Luo et al. 2024)17–19 shows\nabout a 10% density difference at maximum, corresponding to ~0.5% change in planetary radius.\nWe therefore treat FeH as Fe in our calculations. Changes have been made at lines 538-542: The\neffect of H on the density of FeH is minor because hydrogen atoms are much smaller than iron\natoms, and thus the volume change from Fe to FeH is minimal. The effect of using Fe versus\nFeH in the core on the radius is less than ~0.5%, which we quantify by density differences\nbetween the employed Fe EOS and the FeH EOS16.\n8\n4.) In Line 498, the authors refer to using liquid Iron from Luo et al., but I don’t believe\nthat work presented new equations of state for liquids. The authors probably mean they\nuse the same prescription as in Luo, but it reads as if new measurements were used for the\nEoS from Luo.\n[Author] The Luo et al.19 paper does in fact include new equations of state for both dry liquid Fe\nand hydrous liquid Fe.\nMinor/Stylistic comments:\n5.) The introduction presents some exoplanet demographics as decided facts versus debated\nastrophysics. The second sentence has grammatical issues.\n[Author] We have changed it at lines 43-43: If present in sufficient amounts, …\nThe 3% water in Line 48 should be noted to be for planets with Earth-like core mass\nfractions.\n[Author] We have revised the sentence at line 48: …, especially with Earth-like core mass\nfractions.\nLastly, Paragraph 1 and Paragraph 3 feel contradictory at points in that P1 discusses how\nplanets are relatively dry while P3 discusses diverse water inventories.\n[Author] We have added some changes in the P3 paragraph that takes care of it. Lines 68-70:\nAdditionally, the presence of water and the oxidation of iron are closely associated in magma\noceans, with implications on the metallic core size and mantle composition. Given the diverse\nfates of initial water inventories, …\n6.) The text in Figure 5 is very small.\n[Author] We have made the text larger in Figure 5.\n9\n7.) In Fig. 5 caption Line 328 there’s a quick smile-inducing typo.\n[Author] It should be “Mass-radius”. Modification made.\nReferences:\nproperty of iron oxyhydroxide across the spin crossover under extreme conditions,” Phys. Rev. B\n107(6), 64106 (2023).\nG.W. Collins, J.H. Eggert, and R. Jeanloz, “Shock compression of stishovite and melting of\n10\nsilica at planetary interior conditions,” Science (80-. ). 347(6220), 418–420 (2015).\nand Melt System at High Pressures: Implications for Ultralow-Velocity Zones,” J. Geophys. Res.\nSolid Earth 124(5), 4566–4575 (2019).\nEarth’s deep interior,” (2018).\nL. Cobden, and H. Marquardt, “Iron Content ‐ Dependence of Ferropericlase Elastic Properties\nAcross the Spin Crossover From Novel,” Geophys. Reseach Lett. 51, e2024GL111276 (2024).\nX-ray Emission from a Systematic Investigation of Iron Compounds,” Inorg. Chem. 59(17),\n12518–12535 (2020).\nMetal Oxides, Universiteit Utrecht, 2025.\nanalysis of transition metal L edges,” Micron 41(7), 687–694 (2010).\nmetal ions: A comparison via Kα and Kβ detection,” Phys. Rev. B - Condens. Matter Mater.\nPhys. 56(8), 4553–4564 (1997).\nunder deep lower-mantle conditions and Earth’s oxygen-hydrogen cycles,” Nature 534(7606),\n241–244 (2016).\nof FeOOH,” Nature 547(7662), 205–208 (2017).\nChariton, G. Aprilis, I. Koemets, K. Glazyrin, H.P. Liermann, M. Hanfland, E. Ohtani, N.\nDubrovinskaia, C. McCammon, and L. Dubrovinsky, “Chemical Stability of FeOOH at High\nPressure and Temperature, and Oxygen Recycling in Early Earth History**,” Eur. J. Inorg.\nChem. 2021(30), 3048–3053 (2021).\nDeoxygenation and Dehydration of Ferric Oxyhydroxide in Earth’s Subducting Slabs,” Geophys.\nRes. Lett. 48(17), 1–9 (2021).\nChariton, G. Aprilis, I. Koemets, K. Glazyrin, H.P. Liermann, M. Hanfland, E. Ohtani, N.\nDubrovinskaia, C. McCammon, and L. Dubrovinsky, “Chemical Stability of FeOOH at High\nPressure and Temperature, and Oxygen Recycling in Early Earth History,” Eur. J. Inorg. Chem.\n2021(30), 3048–3053 (2021).\n11\nImplications for Hydrogen in Earth’s Inner Core,” Geophys. Res. Lett. 49(5), 1–9 (2022).\nNeumann, “A new ab initio equation of state of hcp-Fe and its implication on the interior\nstructure and mass-radius relations of rocky super-Earths,” Icarus 313, 61–78 (2018).\nthermal equation of state of iron,” Minerals 10(2), (2020).\nand sub-Neptunes,” Nat. Astron., (2024).\n12"}, {"round": 2, "reviewer_comments": {"reviewer_1": "Authors fully addressed my comments and now it is clear enough to publish it with this revised version.", "reviewer_2": "I recommend for this interdisciplinary work on planet interiors for publication. This work continues to evolve our understanding of water in the interior of planets and the complex chemical reactions that do have sizeable impacts on our inferences from observations and the habitability of planets. The authors have addressed my concerns well, and I appreciate their detailed comments and consideration of my review. I do not require seeing any further changes, but I do have one comment on the new Figure S10.\nFigure S10 has one or two density inversions. These are very small, but likely unphysical. I believe they are caused by the inclusion of Fe in the melt while the solid region of the mantle remains less dense than the densest melt. I don't find this second order effect to change the results, but any inversions should be mentioned in at least the caption."}, "author_response": "Point-by-point response to the reviewers’ comments:\nReviewers’ comments are presented in black bold font, while our point-by-point responses are\nin blue. Revisions to the manuscript are highlighted in yellow.\nREVIEWERS' COMMENTS\nReviewer #1 (Remarks to the Author):\nAuthors fully addressed my comments and now it is clear enough to publish it with this\nrevised version.\n[Author] Thank you for your recommendation of publication.\nReviewer #2 (Remarks to the Author):\nI recommend for this interdisciplinary work on planet interiors for publication. This work\ncontinues to evolve our understanding of water in the interior of planets and the complex\nchemical reactions that do have sizeable impacts on our inferences from observations and\nthe habitability of planets. The authors have addressed my concerns well, and I appreciate\ntheir detailed comments and consideration of my review. I do not require seeing any\nfurther changes, but I do have one comment on the new Figure S10.\n[Author] Thank you for your recommendation of publication.\nFigure S10 has one or two density inversions. These are very small, but likely unphysical. I\nbelieve they are caused by the inclusion of Fe in the melt while the solid region of the\nmantle remains less dense than the densest melt. I don't find this second order effect to\nchange the results, but any inversions should be mentioned in at least the caption.\n[Author] We have mentioned the inversions in the caption of Figure S10.\nAt lines 231-233 in SI:\nFigure S10. Density profiles of three cases in our planetary interior model. The very small\ndensity inversion is unphysical and is caused by our assumption that we switch from FeO to the\nless dense FeOOH phase for pressures above 86 GPa."}]

Decision letter