Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffuse fraction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Soot formation and precursor chemistry in Counterflow flames of aviation fuel surrogates

To meet market demands, the aviation sector is interested in utilizing drop-in Synthetic Aviation Turbine Fuels (SATF), either as neat fuels or in blends with conventional Jet A. SATF currently approved in standard specifications may have lower aromatic content with significant fractions of normal, branched, and cyclo-alkanes. Fundamental studies on soot formation from aviation fuels (Jet A, SATF) and their surrogate components are essential to understand how fuel composition influences soot and soot precursor formation. Here, this study reports new measurements of polycyclic aromatic hydrocarbons (PAH) and soot in counterflow diffusion flames (CDFs) of aviation fuel surrogates. Both intrusive and non-intrusive diagnostics are employed to determine the profiles of temperature, gas phase species, PAHs (up to C16), and soot volume fraction (SVF) in CDFs of iso-octane and surrogate mixtures. These measurements shed light on the transition of soot precursors to primary soot particles. In addition to serving as a common surrogate component in Jet A surrogate mixtures, iso-octane is a template species for larger, less volatile branched alkanes found in SATF mixtures. The newly developed Lawrence Livermore National Laboratory (LLNL) PAH and soot model successfully captures temperature, precursor species, and SVF profiles for the mixtures and conditions discussed in this work. Finally, a high-fidelity surrogate for Jet A is proposed that matches targeted physical and chemical properties well, while leveraging the wide range of candidate fuel molecules available in the LLNL detailed chemical model. The proposed surrogate formulation is validated against newly acquired measurements of the surrogate and literature measurements of Jet A. These new experiments and simulations provide critical insights into the PAH and soot formation from aviation fuels. Reaction pathways which require further investigation are highlighted, such that future work may bridge the remaining quantitative gaps in predicting soot formation from aviation fuel surrogates and surrogate components.

Aviation fuels↗
DOE OSTI2025

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗
DOE OSTI2025

A Route to Design Novel Functional Peptides by Applying a Denoising Diffusional Model to mRNA Display Libraries

In vitro directed evolution techniques, such as mRNA display, enable peptide ligand discovery and optimization. However, physical libraries that rely on a genetic code can only search a small fraction of sequence space due to inherent biases in the genetic code and experimental limitations. To address this challenge, denoising diffusion implicit models (DDIMs) are applied to generate novel peptide ligands against B‐cell lymphoma extra‐large (Bcl‐x L ), a key cancer target. Starting with high‐throughput sequencing data from previous selections, a DDIM is trained to produce novel sequences with high affinity binding. Experimental validation confirms that most generated sequences are functionally equivalent to the original library members for Bcl‐x L binding and demonstrated comparable binding kinetics and affinity relative to the wildtype and nearest original neighbors. Importantly, this approach generated rare sequences not easily accessible via mutation and directed evolution. These results indicate that DDIMs can complement and expand directed evolution data, efficiently exploring underrepresented regions of sequence space. This approach provides a broadly applicable framework for accelerating ligand discovery and optimizing molecular properties across diverse targets.

Qi, Pearl [Mork Family Department of Chemical Engi↗
DOE OSTI2020

On the solid-state-bonding mechanism in friction stir welding

This study is a critical assessment of various solid-state-bonding mechanisms is established for friction stir welding (FSW) processes of engineering alloys. The commonly assumed sintering-like diffusional-bonding hypothesis is criticized in this work as not the dominant mechanism. For the wide spectrum of material constitutive laws and FSW processing conditions examined and employed in realistic applications, the thermomechanical history on the workpiece–workpiece interface traverses in the creep-dominated regime for the growth/shrinkage of interfacial cavities. The evolution of the bonding fraction relies mainly on the creep strain rate in the adjourning workpieces, weakly on stress triaxiality, and negligibly on interfacial diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Integrated simulation of U-10Mo monolithic fuel swelling behavior

Here, a separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to take into account various processes occurring in RTR fuels during irradiation. In order to accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions vs. fission density) predicted using the phase-field method. The remainder of fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed in order to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Measurement of Three Transport Coefficients and the Thermodynamic Factor in Block Copolymer Electrolytes with Different Morphologies

The design and engineering of composite materials is one strategy to satisfy the materials needs of systems with multiple orthogonal property requirements. In the case of rechargeable batteries with lithium metal anodes, the system requires a separator with fast lithium ion transport and good mechanical strength. In this work, we focus on the system polystyrene-block-poly(ethylene oxide) (SEO) with bis(trifluoromethane)sulfonimide lithium salt (LiTFSI). Ion transport occurs in the salt-containing poly(ethylene oxide)-rich domains. Mechanical rigidity arises due to the glassy nature of polystyrene (PS). If we assume that the salt does not interact with the PS-rich domains, we can describe ion transport in the electrolyte by three transport parameters (ionic conductivity, κ, salt diffusion coefficient, D, and cation transference number, t + 0 ) and a thermodynamic factor, T f . By systematically varying the volume fraction of the conducting phase, Φ c between 0.29 and 1.0, and chain length, N between 80 and 8000, we elucidate the role of morphology on ion transport. We find that κ is the strongest function of morphology, varying by three full orders of magnitude, while D is a weaker function of morphology. To calculate t + 0 and T f , we measure the current fraction, ρ + , and the open circuit potential, U, of concentration cells. We find that ρ + and U follow universal trends as a function of salt concentration, regardless of chain length, morphology, or Φ c , allowing us to calculate t + 0 for any SEO/LiTFSI or PEO/LiTFSI mixture when κ and D are known. Finally, the framework developed in this paper enables predicting the performance of any block copolymer electrolyte in a rechargeable battery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Tailoring Charge Transport and Magnetism in Complex Half-Heusler/Full-Heusler Nanocomposites (Final Report)

This project exploits (1) the facile co-crystallization of the half-Heusler (HH) and full-Heusler (FH) structures, (2) the high coherency and stability of the HH/FH interfaces, and (3) the solid-state inter-conversion between HH and FH structures through local atomic diffusion, to investigate, for selected compositions of the HH and FH phases, the effect of size, dimensionality, dispersion and volume fraction of FH nanostructures on the electronic, phonon and magnetic behavior of the resulting HH/FH nanocomposites. Our goal is to identify and control key material and growth parameters governing phase formation, microstructural evolution, and to understand the mechanism by which these changes of the internal structure modify the material’s performance. To elucidate the mechanism by which the FH nanostructures regulate electronic charge transport and induce ferromagnetism within non-magnetic and magnetic semiconducting HH matrices, suitable series of HH/0D-FH and HH/2D-FH nanocomposites will be fabricated using synthesis conditions suggested by theoretical calculations of the formation energy and stability of various FH phases (varying M) in HH matrices with known chemical composition. The structure of the HH/FH interfaces (size, chemical composition, volume fraction and dispersion) will be investigated along with the evaluation of the electronic (carrier density, mobility, thermopower, electrical conductivity, effective mass), thermal and magnetic (susceptibility, magnetization etc.) properties of the resulting HH/FH nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A time-dependent chloride diffusion model for predicting initial corrosion time of reinforced concrete with slag addition

The effect of granulated blast furnace slag (GBFS) addition on the threshold chloride concentration (TCC) for rebar corrosion was investigated. A modified diffusion model, coupled with a time-dependent effective diffusion coefficient and surface chloride concentration, was proposed to predict the chloride profile. The corrosion initiation time was estimated based on the model predictions and the measured TCC. The results indicate that adding GBFS decreases the TCC by lowering the pH value of the pore solution. The evolution of corrosion potential and current density is found to obey a 3-parameter Weibull distribution. MIP tests show that adding GBFS contributes to refinement of pore structure by decreasing the fraction of large capillary pores. The time-dependent model exhibits good predictive strength and helps understand how GBFS addition delays the corrosion initiation by retarding the chloride diffusion, though a lower TCC is obtained.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Enhancement of Ni–(Y 2 O 3 ) 0.08 (ZrO 2 ) 0.92 fuel electrode performance by infiltration of Ce 0.8 Gd 0.2 O 2–δ nanoparticles

This paper addresses the use of Ce 0.8 Gd 0.2 O 2–δ (GDC) infiltration into the Ni–(Y 2 O 3 ) 0.08 (ZrO 2 ) 0.92 (YSZ) fuel electrode of solid oxide cells (SOCs) for improving their electrochemical performance in fuel cell and electrolysis operation. Although doped ceria infiltration into Ni–YSZ has recently been shown to improve the electrode performance and stability, the mechanisms defining how GDC impacts electrochemical characteristics are not fully delineated. Furthermore, the electrochemical characteristics have not yet been determined over the full range of conditions normally encountered in fuel cell and electrolysis operation. Here we present a study of both symmetric and full cells aimed at understanding the electrochemical mechanisms of GDC-modified Ni–YSZ over a wide range of fuel compositions and temperatures. Single-step GDC infiltration at an appropriate loading substantially reduced the polarization resistance of Ni–YSZ electrodes in electrolyte-supported cells, as measured using electrochemical impedance spectroscopy (EIS) at various temperatures (600–800 °C) in a range of H 2 O–H 2 mixtures (3–90 vol% H 2 O). Fuel-electrode-supported cells had significant concentration polarization due to the thick Ni–YSZ supports. Here, a distribution of relaxation times approach is used to develop a physically-based electrochemical model; the results show that GDC reduces the reaction resistance associated with three-phase boundaries, but also appears to improve oxygen transport in the electrode. Increasing the H 2 O fraction in the H 2 –H 2 O fuel mixture reduced both the three-phase boundary resistance and the gas diffusion resistance for Ni–YSZ; with GDC infiltration, the electrode resistance showed less variation with fuel composition. GDC infiltration improved the performance of fuel-electrode-supported full cells, which yielded a maximum power density of 2.28 W cm –2 in fuel cell mode and an electrolysis current density at 1.3 V of 2.22 A cm –2 , both at 800°C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Simulation of the SPARC plasma boundary with the UEDGE code

Here in this work, the UEDGE edge transport code is used to examine conditions in the SPARC divertor and edge plasma for various levels of carbon impurity and power from the core (P SOL ). A double-null magnetic configuration is simulated assuming up-down symmetry in geometry and physics. The anomalous heat and particle transport coefficients are tuned to match empirical predictions for SPARC's midplane density profiles, target plate heat flux profiles, and inner/outer divertor power sharing. Convective transport is included on the low-field side, while on the high-field side the transport is modeled as purely diffusive. Hydrogen neutrals are modeled as a fluid with inertial effects, and a carbon impurity is included using the fixed-fraction model. We find that detachment induced by impurity seeding could significantly reduce the heat flux to the divertor surfaces in the SPARC tokamak. At P SOL = 28 MW (the value predicted for SPARC's full-power H-mode scenario) cases with both divertor legs detached were obtained with a carbon impurity fraction between 0.3%–1.4%, far below Z eff limits for SPARC. When the plasma in the outer leg is detached, the peak heat flux density perpendicular to the target plate is below 1 MW m -2 , electron and ion temperatures are less than 1.5 eV, and momentum detachment is observed. However, the detachment state is found to be sensitive to the side-wall boundary conditions, the level of neutral pumping, and the target plate tilt. Finally, a broadly similar SOLPS simulation of SPARC is used to assess the appropriateness of the simpler impurity and neutral models used in UEDGE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The effects of turbulence on the flame structure and NO formation of ammonia turbulent premixed combustion at various equivalence ratios

Ammonia is carbon-free and is regarded as a potential fuel to address global warming issues. Here, in this work, three-dimensional direct numerical simulations (DNS) of ammonia/air turbulent premixed flames were performed to explore the influence of turbulence and equivalence ratio on the flame structure and NO formation characteristics. Two equivalence ratios were considered, i.e. Φ=0.9 and Φ=1.1. The general flame structures were presented and species distributions were examined. The NO mass fraction was found to be the highest in the product for the lean case and in the reaction zone for the rich case. The conditional mean values of species mass fractions and reaction rates were compared with those of the unstrained and strained laminar premixed flames to explore how well the laminar flame structures can approach those of the turbulent flame. The budget analysis of the species transport equations showed that turbulent diffusion plays an important role in species transport. The turbulent diffusivity D T was estimated using the gradient hypothesis based on the DNS data. Various laminar flame simulations with different diffusivities were carried out. It was shown that the conditional means of the DNS agree well with those of the laminar flames including D T in the transport property calculation. The global and local NO formation characteristics were investigated. It was found that the mean NO production conditioned on the progress variable is lower compared with the corresponding laminar flame in the rich case. However, the relative contributions from various NO pathways are rarely affected by turbulence. The NO mass fraction is higher in negative curvature regions compared with positive curvature regions of the flame surface for the rich case, which is due to the preferential diffusion of H 2 and other radicals and the enhanced NO pathways in negatively curved regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Pectin-associated mannans and xylans play distinct roles in cell-cell adhesion in pine and poplar wood

Wood particles of pine required longer treatments of acidic chlorite and dilute alkali to achieve similar proportions of cell separation as of poplar. Chemical extraction and enzymatic digestion studies indicated that rhamnogalacturonan-I (RG-I), mannan and xylan contribute to cell-cell adhesion in pine wood. Mannan associated with lignin in an acidic chlorite fraction and with xylan in a dilute alkali fraction. RG-I and xylan were extracted in both chlorite and dilute alkali in poplar. Scanning electron microscopy and a probe diffusion study using Fluorescence Recovery After Photobleaching showed that lignin, RG-I and mannan initiated cell separation at cell corners and dilute alkali extraction of xylan was required for complete cell separation in pine. Cell-cell separation increased enzymatic saccharification within 24 h by 2-fold in poplar and 3.4-fold in pine. Finally, our results inform a strategy to reduce particle size in woody biomass and provide insights into different molecular bases of cell-cell adhesion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Emergent Properties from Three-Dimensional Assemblies of (Nano)particles in Confined Spaces

The assembly of (nano)particles into compact hierarchical structures yields emergent properties not found in the individual constituents. The formation of these structures relies on a profound knowledge of the nanoscale interactions between (nano)particles, which are often designed by researchers aided by computational studies. These interactions have an effect when the (nano)particles are brought into close proximity, yet relying only on diffusion to reach these closer distances may be inefficient. Recently, physical confinement has emerged as an efficient methodology to increase the volume fraction of (nano)particles, rapidly accelerating the time scale of assembly. Specifically, the high surface area of droplets of one immiscible fluid into another facilitates the controlled removal of the dispersed phase, resulting in spherical, often ordered, (nano)particle assemblies. In this review, we discuss the design strategies, computational approaches, and assembly methods for (nano)particles in confined spaces and the emergent properties therein, such as trigger-directed assembly, lasing behavior, and structural photonic color. Finally, we provide a brief outlook on the current challenges, both experimental and computational, and farther afield application possibilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Martensitic transformation in superlattices of two non-transforming metals

A number of studies have shown that the coherent integration of a second non-transforming phase in martensitic materials results in unprecedented thermo-mechanical properties. We demonstrate the possibility of martensitic transformation by epitaxially combining two non-transforming materials, Mg (hexagonal closed packed) and MgLi (body centered cubic). Density functional theory calculations predict that MgLi/Mg superlattices undergo a transformation crystallographically equivalent to MgSc alloys and can be expected to result in finite temperature shape memory behavior. We predict the zero-temperature energy of the martensite phase to be 9.3 meV/atom lower than austenite for a 50–50 at. % MgLi/Mg nanolaminate, which is expected to result in a martensitic transformation temperature of approximately 150 K. This energy difference can be tuned between 4.9 and 18.1 meV/atom by varying the fraction of the pure Mg phase in the superlattice, potentially enabling the transformation from cryogenic to room and higher temperatures. While Li diffusion might limit the applications of the systems proposed, this work opens the door to a new approach to create martensitic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Tracer particle transport dynamics in the diffusive sandpile cellular automaton

The confinement properties of the diffusive running sandpile are characterized by tracking the motion of a population of marked grains of sand. It is found that, as the relative strength of the avalanching to the diffusive transport channel is varied, a point is reached at which the particle global confinement time and the probability density functions of the jump-sizes and waiting-times of the tracked grains experience a sudden change, thus revealing a dynamical transition, that is consistent with previous studies (Newman DE et al., Phys Rev Lett 2002;88(20):204304). Across this transition, the sandpile moves from a regime characterized by self-similarity and memory, where avalanches of all possible sizes dominate transport across the system, to another regime where transport is taken over by near system-size, quasi-periodic avalanches. Values for the fractional transport exponents that quantify effective transport across the sandpile prior to the transition are also obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Detailed mixing measurements from single-hole converging ECN injectors using Rayleigh scattering

Time resolved liquid and vapor fields of dodecane and oxymethylene ethers are measured from Spray A-3 and Spray D using high speed Rayleigh scattering and diffuse back illumination at the Engine Combustion Network (ECN) Spray A condition of 900 K and 22.8 kg/m 3 . Global quantities including mixture fraction, vapor and liquid penetration, as well as spreading angle are measured. The mixture fraction fields and vapor penetration profiles are well predicted by the 1-D Musculus-Kattke model. The mixture fraction field and vapor penetration from Spray A-3 are similar to those measured from Spray A in previous works. Spray D exhibits higher mixture fraction fields and vapor penetration due to its larger nozzle diameter. The quasi-steady mixture fraction fields from these injectors scale well when distance from the injector is normalized by the nozzle diameter. The turbulent dissipation structures were also analyzed based on the orientation, thickness, and magnitude of the mixing layers. Furthermore, the orientation and thickness are similar to other measurements in atmospheric gas jets, while the magnitude is lower. The thickness and magnitude are subject to uncertainties due to limitations in the imaging resolution of the system but still provide an order of magnitude as a useful reference for comparison against computational fluid dynamic simulations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Phase behaviours of superionic water at planetary conditions

Most water in the Universe may be superionic, and its thermodynamic and transport properties are crucial for planetary science but difficult to probe experimentally or theoretically. We use machine learning and free-energy methods to overcome the limitations of quantum mechanical simulations and characterize hydrogen diffusion, superionic transitions and phase behaviours of water at extreme conditions. Furthermore, we predict that close-packed superionic phases, which have a fraction of mixed stacking for finite systems, are stable over a wide temperature and pressure range, whereas a body-centred cubic superionic phase is only thermodynamically stable in a small window but is kinetically favoured. Our phase boundaries, which are consistent with existing—albeit scarce—experimental observations, help resolve the fractions of insulating ice, different superionic phases and liquid water inside ice giants.

36 MATERIALS SCIENCE↗