---
title: 'Hydrogen–Water Demixing: Planetary & Materials Insights'
url: https://www.emergentmind.com/topics/hydrogen-water-demixing
type: topic
---

# Hydrogen–Water Demixing: Planetary & Materials Insights

Searching arXiv for recent papers on hydrogen–water demixing and related planetary interior studies.
Hydrogen–water demixing is the phase separation of a nominally mixed $\mathrm{H_2}$–$\mathrm{H_2O}$ system into coexisting hydrogen-rich and water-rich phases under conditions where mixing is thermodynamically disfavored. In planetary science, the topic is primarily associated with the interiors of Uranus and Neptune, where the miscibility of major volatile constituents affects interior stratification, gravitational harmonics, magnetic-field interpretation, thermal evolution, and atmospheric composition constraints [2012.04166; 2410.21099; 2507.06288]. In condensed-matter and cryogenic materials research, hydrogen removal from hydrogen-filled ice provides a distinct but related manifestation of hydrogen–water separation, producing a metastable porous water framework known as ice XVII [1607.07617]. Across these settings, the central question is whether hydrogen and water remain a single phase or undergo immiscibility governed by the Gibbs free energy of mixing, chemical-potential equalities, and the geometry of binodal and spinodal boundaries [2410.21099; 2507.06288].

## 1. Thermodynamic definition and phase-separation criteria

For a binary mixture of hydrogen and water at pressure $P$, temperature $T$, and composition $x$, demixing is formulated through the Gibbs free energy of mixing. One representation writes
$$
\Delta G_{\rm mix}(x,T,P)=G_{\rm mix}-\left[x\,G_{\mathrm{H_2}}+(1-x)\,G_{\mathrm{H_2O}}\right],
$$
with $x$ taken as the $\mathrm{H_2O}$ particle fraction in one formulation [2410.21099] and as the water mole fraction $x\equiv x_{\mathrm{H_2O}}$ in another [2507.06288]. In the planetary demixing literature summarized here, phase separation occurs when the free-energy surface develops a composition range associated with an immiscibility gap, and coexistence is determined by equality of chemical potentials between the two phases [2410.21099; 2507.06288].

A detailed statement of the coexistence condition uses
$$
\mu_{\mathrm{H_2}}(A)=\mu_{\mathrm{H_2}}(B), \qquad
\mu_{\mathrm{H_2O}}(A)=\mu_{\mathrm{H_2O}}(B),
$$
where $A$ is water-poor and $B$ is water-rich [2410.21099]. The same condition is equivalently expressed through a common-tangent construction in the composition-dependent free energy, formulated as solving $\mu_1(x_1)=\mu_1(x_2)$ and $\mu_2(x_1)=\mu_2(x_2)$ [2507.06288]. The spinodal is defined by the vanishing of the curvature of the Gibbs free energy with respect to composition, $\partial^2 G/\partial x^2|_{P,T}=0$ [2507.06288]. One summary also states that the immiscibility-gap criterion is $\partial^2G_{\rm mix}/\partial x^2>0$ over some range in $x$, with the locus where $\partial^2G_{\rm mix}/\partial x^2=0$ defining the spinodal and the maximum of the coexistence curve constituting a critical point $(x_c,T_c,P_c)$ [2410.21099]. This discrepancy in sign convention reflects differing definitions and summaries of the free-energy construction rather than a settled controversy in the supplied material.

A simplified thermodynamic parameterization for warm sub-Neptunes adopts a Flory–Huggins–style expression,
$$
G_{\rm mix}(x,T,P)=RT\bigl[x\ln x+(1-x)\ln(1-x)\bigr]+x(1-x)\chi(T,P),
$$
with $x\equiv x_{\rm H_2}$ and $\chi(T,P)$ an interaction parameter [2512.01805]. In that representation, absolute instability is written as
$$
\frac{\partial^2 G_{\rm mix}}{\partial x^2}
= \frac{RT}{x(1-x)} - 2\chi(T,P) < 0,
$$
or equivalently $\chi(T,P) > RT/[2x(1-x)]$ [2512.01805]. The same source states that, in practice, it does not adopt the toy Flory–Huggins parameterization directly, but instead uses ab-initio-derived coexistence curves.

## 2. Phase diagrams, critical curves, and experimental-computational disagreements

Two lines of evidence for possible hydrogen–water immiscibility in ice-giant interiors were identified in work on thermodynamically governed interior models of Uranus and Neptune [2012.04166]. The first arises from crude extrapolation of the experimental hydrogen–water critical curve to $\sim 3$ GPa using data obtained for an impure system containing silicates, as reported by Bali et al. (2013); the same source notes that Uranus and Neptune could also be “dirty” [2012.04166]. The second invokes reasoning based on the gravitational and magnetic fields [2012.04166]. That work also states that current ab initio models disagree and cites Soubiran and Militzer (2015), while remarking that hydrogen and water are difficult to model from first-principles quantum mechanics with the necessary precision [2012.04166].

A more systematic planetary treatment constructs seven $\mathrm{H_2}$–$\mathrm{H_2O}$ phase diagrams from available experimental and computational data [2410.21099]. These are the SFB-linear-3 GPa, SFB-linear-4 GPa, and SFB-linear-5 GPa extrapolations of Seward and Franck (1981) and Bali et al. (2013); the V23 flat, V23 conv–1800 K, and V23 conv–2000 K extensions of Vlasov et al. (2023); and Berg24, based on ab initio DFT-MD results from Bergermann et al. (2024) up to $\sim 12$ GPa [2410.21099]. For each critical curve $T_n(P)\equiv T_{\rm dmx}(P,x=0.5)$, the low-pressure U-shape of the $0.2$ GPa binodal observed by Seward and Franck (1981) is shifted vertically by $\Delta T=T_n(P)-T_n(0.2\ {\rm GPa})$, yielding approximate isobaric boundaries $x_A(P,T)$ and $x_B(P,T)$ [2410.21099]. Clausius–Clapeyron,
$$
\frac{dP}{dT}=\frac{\Delta S}{\Delta V},
$$
is given but not explicitly required in the interpolation procedure [2410.21099].

A later study of Uranus, Neptune, K2-18 b, and TOI-270 d uses recent ab initio calculations and an analytic fit to the demixing temperature:
$$
T_{\rm mix}(P,x)=a\cdot
\left[
\frac{1}{\pi}\cdot
\frac{0.5\,(b+cP)}
{(x-d)^2+(0.5\,b)^2}
\right]
\cdot(eP^3+fP^2+gP+h)+iP,
$$
with coefficients $a=1.2035\times10^{-4}$, $b=0.5501$, $c=1.9163\times10^{-2}$, $d=0.4498$, $e=-6.2253\times10^{-2}$, $f=74.5041$, $g=-3.1495\times10^{-4}$, $h=5.0828\times10^6$, and $i=4.0719$ [2507.06288]. The same work introduces constant temperature offsets,
$$
T'_{\rm mix}(P,x)=T_{\rm mix}(P,x)+T_{\rm offset},
$$
with explored values $T_{\rm offset}=0$, $+600$ K, and $+1100$ K, and for exoplanets up to $+500$ K, in order to bracket uncertainties in the miscibility gap [2507.06288].

For warm sub-Neptunes, a merged low-pressure and high-pressure critical-curve construction is described as
$$
T_{\rm crit}(P)\approx
\begin{cases}
T_{\rm Seward}(P),&P\lesssim 200\ {\rm bar}\\
T_{\rm Gupta}(P),&1\ {\rm kbar}\lesssim P\lesssim 1\ {\rm Mbar}\\
\text{(linear interp.)},&200\ {\rm bar}<P<1\ {\rm kbar},
\end{cases}
$$
and full $P$–$T$ coexistence curves are computed for metallicities from $1\times$ to $1000\times$ solar [2512.01805]. At $Z\sim0.3\,(100\times)$ solar, the dome of immiscibility peaks at $T\sim1600$ K around $P\sim10$ kbar, whereas at $Z\sim0.75\,(300\times)$ solar the dome peaks above $T\sim3000$ K near $P\sim30$ kbar [2512.01805].

## 3. Demixing in Uranus and Neptune

Interior models that assume discrete layers are only directly justified if the major constituents are immiscible; otherwise, diffuse interfaces may arise from accretion that centrally concentrates the least volatile and most dense constituents, with resulting compositional gradients likely inhibiting convection [2012.04166]. Within that framework, hydrogen–water immiscibility has been treated as a candidate explanation for the contrasting internal properties of Uranus and Neptune [2012.04166; 2410.21099].

Adiabatic structure calculations compare planetary adiabats $T_{\rm ad}(P;Z)$ to the various demixing boundaries to infer the onset and depth of phase separation [2410.21099]. The onset occurs where the planetary adiabat intersects the phase boundary, stated as $P\approx4$–$11$ GPa and $T\approx260$–$310$ K [2410.21099]. Rain-out then proceeds until the adiabat of the depleted outer layer just grazes the binodal at a single transition pressure $P_Z$ [2410.21099]. In this formalism, $Z_{\rm atm}$ is the outer-envelope water mass fraction after demixing, $Z_{\rm deep}$ is the deep-interior water mass fraction, and $P_Z$ is the pressure of the sharp transition between the water-poor and water-rich layers [2410.21099].

For Uranus with $T_{1{\rm bar}}\approx76$ K, the inferred upper limit is $Z_{\rm atm}\le 0.21$ with $P_Z\approx4$–$11$ GPa [2410.21099]. For Neptune with $T_{1{\rm bar}}\approx72$ K, the inferred upper limit is $Z_{\rm atm}\le 0.16$ with the same $P_Z\approx4$–$11$ GPa range [2410.21099]. The sensitivity to phase-diagram choice is explicit: SFB-linear models are shallow with $P_Z\lesssim5$ GPa, V23 extensions give $P_Z\approx4$–$6$ GPa, and Berg24 yields the deepest $P_Z\approx11$ GPa and the widest tangential region, $4$–$11$ GPa [2410.21099].

An earlier thermodynamic interior study drew a different asymmetry between the two planets. It found that Neptune models with envelopes containing a substantial water mole fraction, as much as $\chi \gtrsim 0.1$ relative to hydrogen, can satisfy observations, whereas Uranus models appear to require $\chi \lesssim 0.01$, potentially suggestive of fully demixed hydrogen and water [2012.04166]. The same study argued that different hydrogen–water demixing states could account for the different heatflows of Uranus and Neptune [2012.04166]. This suggests that the sharper depletion inferred for Uranus in some models is linked not only to present-day composition but also to a specific thermodynamic and evolutionary pathway.

## 4. Interior structure, gravitational harmonics, and layered envelopes

The adiabatic-structure approach for the ice giants combines equations of state for an H/He mixture, water, and rock with hydrostatic equilibrium and mass continuity:
$$
\frac{dP}{dr}=-\rho(r)\,g(r)=-\rho(r)\frac{Gm(r)}{r^2},
\qquad
\frac{dm}{dr}=4\pi r^2\rho(r),
$$
and an adiabatic gradient $\nabla_{\rm ad}\equiv(\partial\ln T/\partial\ln P)_S$ computed from the equation of state mixture including ideal-mixing entropy [2410.21099]. The H/He equation of state is given as the SCvH-like EoS by Chabrier and Debras (2021), the water equation of state as AQUA EoS, and the rock core as the Hubbard and Marley (1989) silicate/iron mixture [2410.21099].

The resulting density structures are evaluated against the observed gravitational harmonics. The even zonal harmonics are written as
$$
J_{2n}=-\frac{1}{M\,R^{2n}}
\int \rho(r,\theta)\,r^{2n}\,P_{2n}(\cos\theta)\,d^3r.
$$
In a case with a water-only deep envelope adjusted to fit the observed $J_2$, the comparison with $J_4$ discriminates among demixing prescriptions [2410.21099]. SFB-linear-3 GPa models with $P_Z\sim3$ GPa yield $|J_4|$ larger than observed, including dynamic wind correction, and are therefore excluded [2410.21099]. SFB-linear-5 GPa and V23 conv–2000 K models with $P_Z\approx5$–$6$ GPa can just match Neptune’s $J_4$ but overpredict Uranus’s $|J_4|$ unless $P_Z\gtrsim10$ GPa [2410.21099]. Berg24 with $P_Z\approx11$ GPa matches both $J_2$ and $J_4$ within uncertainties [2410.21099].

The inferred deep compositions depend on whether rocks are permitted below a “rock-cloud” level. In a water-only deep envelope, $Z_{\rm deep}(\mathrm{H_2O})\approx0.68$–$0.87$ for Uranus and $0.73$–$0.90$ for Neptune [2410.21099]. If a rock fraction is allowed below a “rock-cloud” level at $T\approx2000$ K and $20$–$30$ GPa, with the ice-to-rock ratio fixed to $0.5\times$ solar ($I:R\approx1.35$), the enhanced central condensation lowers $|J_4|$ and improves agreement [2410.21099]. In that case, $Z_{\rm deep}(\mathrm{H_2O})\approx0.44$–$0.48$ and $Z_{\rm deep}(\mathrm{rocks})\approx0.32$–$0.36$ for Uranus, while Neptune yields $Z_{\rm deep}(\mathrm{H_2O})\approx0.45$–$0.64$ and $Z_{\rm deep}(\mathrm{rocks})\approx0.34$–$0.24$ [2410.21099].

The same work emphasizes that a sharp transition at $\sim4$–$11$ GPa emerges because the binodal is nearly vertical in $x$, and that gradual layering would require more complex binodal shapes at these pressures [2410.21099]. This is significant because it links a microscopic phase boundary directly to the macroscopic legitimacy of “classical few-layer models” for Uranus and Neptune [2410.21099].

## 5. Thermal evolution, rain-out energetics, and planetary consequences

Hydrogen–water demixing is not only a structural effect but also an energy source. In thermodynamically governed interior models of Uranus and Neptune, enough gravitational potential energy is available from gradual hydrogen–water demixing to supply Neptune’s present-day heatflow for roughly ten solar system lifetimes [2012.04166]. The same study states that hydrogen–water demixing could slow Neptune’s cooling rate by an order of magnitude [2012.04166]. Within the scope of the supplied material, these statements are presented as consequences of gradual phase separation and settling rather than of a transient catastrophic event.

A later evolutionary treatment couples phase separation directly into a 1D structure–energy solver, CEPAM, using
$$
\frac{\partial P}{\partial m}=-\frac{Gm}{4\pi r^4},
\qquad
\frac{\partial r}{\partial m}=\frac{1}{4\pi r^2\rho},
$$
and the internal-energy equation
$$
\frac{\partial u}{\partial t}+P\,\frac{\partial(1/\rho)}{\partial t}
=-\frac{\partial L}{\partial m}+\sum_i \mu_i\,\frac{\partial X_i}{\partial t},
$$
where the last term accounts for chemical work associated with composition changes [2507.06288]. At each timestep, layers satisfying $T(P)<T'_{\rm mix}(P,x_{\rm local})$ are flagged as unstable; the local water mass fraction is reduced to its saturation value, and the excess water is instantaneously redeposited below, maintaining a smooth and monotonic water-versus-depth profile [2507.06288]. The chemical potential difference appears as a positive source term and physically includes both latent heat release and gravitational potential energy as water sinks [2507.06288]. The additional energy raises the intrinsic luminosity and can slow contraction or even cause transient radius inflation [2507.06288].

The following summary organizes the planetary outcomes explicitly stated for the evolutionary models [2507.06288]:

| Planet | $T_{\rm offset}$ | Stated consequence |
|---|---:|---|
| Uranus | $0$ K | no demixing; $Z_{\rm atm}$ remains $0.40\rightarrow0.40$ |
| Uranus | $+600$ K | outer $3\%$ mass fully depleted; $\Delta R/R_0=+3\%$ |
| Uranus | $+1100$ K | outer $16\%$ mass fully depleted; $\Delta R/R_0=+20\%$ |
| Neptune | $0$ K | no demixing; $Z_{\rm atm}$ remains $0.50\rightarrow0.50$ |
| Neptune | $+600$ K | outer $13\%$ depleted; $\Delta R/R_0=+3\%$ |
| Neptune | $+1100$ K | outer $13\%$ fully depleted; $\Delta R/R_0=+20\%$ |

For Uranus, a $+600$ K offset gives onset at $P\approx0.24\times10^2$ kbar at $2$ Gyr and complete depletion of the outer $3\%$ mass, while a $+1100$ K offset yields onset at $P\approx0.24\times10^2$ kbar at $0.5$ Gyr and full demixing of the outer $16\%$ by mass, with a radius increase of nearly $20\%$ [2507.06288]. For Neptune, a $+600$ K offset gives onset at $P\approx1$–$3$ kbar at $1.5$ Gyr with the outer $13\%$ depleted, and a $+1100$ K offset gives onset at $P\approx1$–$20$ kbar at $0.5$ Gyr with full depletion of the outer envelope and again nearly $20\%$ radius increase [2507.06288].

A plausible implication is that thermodynamic uncertainty in the phase diagram propagates directly into uncertainty in the inferred luminosity history, contraction rate, and present-day stratification. That implication is explicitly anticipated in calls for coupled thermal-evolution models including latent heat release and gravitational energy from rain-out to refine cooling-time predictions and address Uranus’s anomalously low luminosity [2410.21099].

## 6. Materials manifestation: hydrogen removal from filled ice and ice XVII

Outside planetary interiors, hydrogen–water demixing also appears in cryogenic solid-state systems. A hydrogen-filled crystalline water compound called C$_0$ filled ice can be emptied to produce a new porous form of ice, termed ice XVII, while retaining the water-lattice framework [1607.07617]. The precursor is synthesized by exposing finely powdered $\mathrm{H_2O}$ ice to $\mathrm{H_2}$ gas at $P\approx430$ MPa and $T\approx255$ K; the sample adsorbs hydrogen above $\approx360$ MPa, and the pressure is maintained for several days to assure full conversion to the C$_0$ phase [1607.07617]. Room-pressure X-ray diffraction at $\approx100$ K is fitted by the C$_0$-II structural model in space group $P3_112$ with lattice constants $a=6.3313(2)$ Å and $c=6.1058(2)$ Å [1607.07617].

Hydrogen release is monitored by Raman spectroscopy in four spectral regions: lattice phonons, H–O–H stretch, $\mathrm{H_2}$ rotational lines, and $\mathrm{H_2}$ vibrons [1607.07617]. As temperature is slowly raised under vacuum, the intensity of the $\mathrm{H_2}$ rotational lines decreases until they vanish after $1$–$2$ h at $T\approx120$ K, with no abrupt shifts or splittings in lattice phonon or OH bands, implying no change of the water-lattice framework [1607.07617]. Hydrogen content is quantified using
$$
I_{\rm rot}=I_0/\sigma_0+I_1/\sigma_1,
\qquad
R\equiv I_{\rm rot}/I_{\rm phon}\propto X\equiv {\rm mol\ H_2/mol\ H_2O},
$$
with calibration giving $X\simeq0.233\,R$ [1607.07617]. Freshly synthesized C$_0$ samples have $R\approx0.108$–$0.112$, hence $X\approx25\%$ [1607.07617].

The equilibrium condition for guest hydrogen is stated as
$$
\mu_{\mathrm{H_2}}^{\rm filled}(T,p_{\rm eq})=\mu_{\mathrm{H_2}}^{\rm gas}(T,p_{\rm eq}),
$$
with ideal-gas chemical potential
$$
\mu_{\mathrm{H_2}}^{\rm gas}(T,p)=\mu^0_{\mathrm{H_2}}(T,p_0)+RT\ln(p/p_0).
$$
From adsorption isotherms at fixed uptake $X$, the adsorption enthalpy is obtained through
$$
\frac{\Delta h}{R}=\left(\frac{\partial \ln(p/p_0)}{\partial(1/T)}\right)_X,
$$
and experimentally $-\Delta h$ decreases from $\approx5$ kJ/mol at $X\approx10\%$ to $\approx2$ kJ/mol at $X\approx40\%$ [1607.07617]. The corresponding Gibbs free-energy change is given as $\Delta G=\Delta H-T\Delta S=RT\ln(p/p_{\rm eq})$ [1607.07617].

Neutron powder diffraction on deuterated ice XVII yields an empty-lattice structure in space group $P6_122$, with lattice constants at $T\approx25$ K of $a\approx6.15$ Å and $c\approx6.12$ Å [1607.07617]. The framework contains helical channels parallel to $c$, with free bore of $\approx5.3$ Å and channel diameter $\approx6.1$ Å [1607.07617]. The same work describes the demixing mechanism as smooth, diffusion-mediated desorption under vacuum, with guest molecules leaving the spiraling channels while the host framework persists [1607.07617]. Upon emptying, the OH-stretch mode downshifts by $\approx20\ {\rm cm^{-1}}$ and phonons upshift by $\approx7\ {\rm cm^{-1}}$, indicating loss of guest-induced strain [1607.07617].

Adsorption–desorption isotherms show strong hysteresis at $15$ K, two kinetic regimes at $40$ K, and nearly reversible behavior for $T\ge65$ K [1607.07617]. The emptied crystal can adsorb hydrogen again and release it repeatedly, and ice XVII can be refilled to $X\gg25\%$ at $T\le40$ K and $p\le0.3$ bar within minutes, with no detectable loss of crystallinity or capacity [1607.07617]. In this materials context, hydrogen–water demixing refers not to liquid immiscibility in a planetary envelope but to removal of guest $\mathrm{H_2}$ from a host water lattice while preserving a metastable porous ice framework.

## 7. Extensions to sub-Neptunes, observational implications, and open questions

Hydrogen–water demixing has been extended from Solar System ice giants to sub-Neptunes. One study finds that demixing may occur in Uranus, Neptune, K2-18 b, and TOI-270 d and could lead to complete depletion of water in the outermost regions of Uranus and Neptune [2507.06288]. For K2-18 b, a temperature offset of $500$ K is required to obtain complete depletion of water in the atmosphere, and the model is proposed as an explanation for the absence of water features in its JWST spectrum [2507.06288]. For TOI-270 d, the same offset yields partial atmospheric depletion, consistent with JWST’s detection of water [2507.06288].

A later atmosphere–interior inference framework, ATHENAIA, argues for “a window for demixing” on warm metal-rich sub-Neptunes [2512.01805]. The atmosphere is modeled with SCARLET and the interior with one-dimensional structure models following Thorngren et al. (2016, 2019), linked by minimizing
$$
\delta_{\rm TPR}
=\sqrt{\bigl[T_{\rm atm}(P_{\rm link})-T_{\rm int}(P_{\rm link})\bigr]^2
+\bigl[R_{\rm atm}(P_{\rm link})-R_{\rm int}(P_{\rm link})\bigr]^2}
$$
at $P_{\rm link}=\min(P_{\rm RCB},1\,{\rm kbar})$ [2512.01805]. For solar-type irradiation levels equivalent to TOI-270 d, the region in which demixing first appears is approximately
$$
Z_{\rm env}\sim0.5\mbox{--}0.95 \quad (\text{100--700}\times\text{solar}),
\qquad
T_{\rm eq}\sim330\mbox{--}500\ {\rm K},
$$
and the window broadens for lower $T_{\rm eq}$, higher mass, or larger envelope mass fraction [2512.01805].

That framework emphasizes the role of adiabatic gradients. The dry adiabatic gradient is
$$
\nabla_{\!ad}=\left(\partial\ln T/\partial\ln P\right)_s
=\frac{\alpha_T P}{\rho c_P},
$$
and water-rich mixtures have systematically smaller $\nabla_{\rm ad}$ than pure H$_2$/He mixtures [2512.01805]. The reported values are $\nabla_{\rm ad}\sim0.28$–$0.30$ for pure H$_2$/He, $\sim0.20$–$0.25$ for $50\%$ by mass H$_2$O, and $\sim0.15$–$0.20$ for $80\%$ H$_2$O in the $1$–$10$ kbar region [2512.01805]. Because a shallower adiabat heats up more slowly with depth, a water-rich envelope can remain below $T_{\rm crit}(P)$ and therefore enter the demixing region [2512.01805].

For TOI-270 d specifically, the inferred posterior is bimodal: either a thin ($1$–$3$ wt %) solar-metallicity envelope or a thick ($30$–$60$ wt %) water-rich one [2512.01805]. The JWST water abundance, stated as $Z_{\rm atm}\sim200\pm100\times$, lies inside the demixing window $Z_{\rm env}\simeq0.5$–$0.95$ and $f_{\rm env}\sim5$–$30\%$, so the planet is argued likely to host compositional gradients, with the true bulk $Z_{\rm env}$ possibly $2$–$3\times$ higher than the photospheric $Z_{\rm atm}$ [2512.01805]. The same models place the envelope–mantle boundary at $P\gtrsim20$ kbar and $T\lesssim2600$ K for $T_{\rm int}=25$ K and $f_{\rm env}\le0.5$, which is below high-pressure silicate melting curves of $3000$–$4000$ K; on that basis, no molten magma ocean is predicted for TOI-270 d [2512.01805].

Several unresolved issues recur across the literature. Consensus remains lacking on the phase boundary itself, because extrapolated experiments, newer experiments, and ab initio calculations do not yield a unique miscibility curve [2012.04166; 2410.21099]. The predicted atmospheric water abundance limits are upper limits for water alone, while real atmospheres also contain CH$_4$, NH$_3$, and He [2410.21099]. If $\mathrm{H_2}$–$\mathrm{H_2O}$ demixing does not occur, alternative explanations for low outer-envelope water include cloud-inhibited convection and progressive planetesimal enrichment during formation [2410.21099]. Other demixing processes—H–He at Mbar pressures, H–C forming diamonds, and MgO–H$_2$O in the deep mantle—may superpose additional layering [2410.21099]. Reduced uncertainty in $J_4$, better estimates of dynamic corrections, atmospheric water-abundance measurements from an orbiter-plus-probe mission, and laboratory and ab initio studies of multicomponent phase behavior remain the stated priorities for testing the hydrogen–water demixing hypothesis [2410.21099].

Source: https://www.emergentmind.com/topics/hydrogen-water-demixing