Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium core”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An experimental study on gas-liquid phase fluid migration in hydrate-bearing sediments during hydrate dissociation

Natural gas hydrate production tests face problems such as severe sand blockage, poor gas-liquid phase separation, and significant land subsidence. This is because of the insufficient understanding of the complex phase transition and gas-liquid multi-phase fluid migration during hydrate dissociation. In hydrate-bearing sediment systems, hydrate phase transition couples with gas-liquid fluid migration. The phase transition causes changes in pore structure, which in turn modifies porous infiltration parameters and fluid flow capacity. Meanwhile, alterations in phase interfaces affect key parameters like surface tension and wettability. Gas-liquid fluid migration influences heat and mass transfer, thus affecting phase equilibrium and dissociation rates. To bridge the gap in describing gas-liquid fluid migration during hydrate dissociation in experiments, this research innovatively integrated an unsteady-state gas displacement by water and a quantitative hydrate dissociation process, independently developed a multi-phase seepage experimental system suitable for hydrate dissociation and determined the relationship between seepage parameters and hydrate saturation under different porosity. The results are as follows: a) Core samples with higher initial porosity show a greater recovery rate of fluid flow capacity. b) The retarding effect of multi-phase fluid has a more significant impact on the migration of the wetting phase fluid (water) than that of the non-wetting phase fluid (methane). c) During hydrate dissociation, the evolution of absolute permeability shows an “S-shaped” pattern, and the evolution of relative permeability shows a “wiring-harness” pattern. In conclusion, the findings can provide a theoretical basis for preventing geological disasters and for geotechnical engineering design during hydrate production.

58 GEOSCIENCES↗
DOE OSTI2025

Vortex-Controlled Quasiparticle Multiplication and Self-Growth Dynamics in Superconducting Resonators

Even in the quantum limit, non-equilibrium quasiparticle (QP) populations induce QP poisoning that irreversibly relaxes the quantum state and significantly degrades the coherence of transmon qubits. A particularly detrimental yet previously unexplored mechanism arises from QP multiplication facilitated by vortex trapping in superconducting quantum circuits, where a high-energy QP relaxes by breaking additional Cooper pairs and amplifying the QP population due to the locally reduced excitation gap and enhanced quantum confinement within the vortex core. Here we directly resolve this elusive QP multiplication process by revealing vortex-controlled QP self-generation in a highly nonequilibrium regime preceding the phonon bottleneck of QP relaxation. At sufficiently low fluence, femtosecond-resolved magneto-reflection spectroscopy directly reveals a continuously increasing QP population that is strongly dependent on magnetic-field-tuned vortex density and absent at higher excitation fluences. Quantitative analysis of the emergent QP pre-bottleneck dynamics further reveals that, although the phonon population saturates within $\simeq$10~ps, both free and trapped QPs continue to grow in a self-sustained manner--hallmarks of the long-anticipated QP-vortex interactions in nonequilibrium superconductivity. We estimate a substantial increase of $\sim$34% in QP density at vortex densities of $\sim$ 100 magnetic flux quanta per $\mathrm{μm^{2}}$. Our findings establish a powerful spectroscopic tool for uncovering QP multiplication and reveal vortex-assisted QP relaxation as a critical materials bottleneck whose mitigation will be essential for resolving QP poisoning and enhancing coherence in superconducting qubits.

Park, Joong M. [Ames Lab]↗
NASA NTRS1990

The earth's forced nutations - Geophysical implications

Theories of nutation are examined critically to develop an extension of the theoretical description based on an earth in nonhydrostatic equilibrium. The estimation of nutation amplitude as a function of frequency is reviewed with reference to contributing variables such as the oceans, mantle anelasticity, and nonhydrostatic structure. Theoretical results by Wahr (1981) are compared to those of VLBI observations, and other VLBI data are used to discuss Wahr and Bergen's (1986) theoretical nutation admittances for an anelastic earth. The contributions associated with the tilt-over mode and with free-core nutations (FCN) are discussed, and the contribution of anelasticity is examined. The VLBI data show that nutation data are important sources of information regarding the earth's interior. Nutation amplitudes are derived with a model for a nonhydrostatically prestressed earth to determine the FCN contributions.

Wahr, J. M.↗
DOE OSTI2025

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Heavy dark matter in rapidly evolving massive stars

We study the impact of heavy dark matter (DM) captured in massive stars via scattering(s) with the star constituents. We focus on the first stars and use stellar evolution simulations to track down how DM capture evolves over time from the zero-age main sequence to the late metal-rich stages of stellar evolution. During the early hydrogen-helium-dominated phase, the capture process is well described by scattering with two targets. As a star evolves, metal production leads to the formation of a dense core surrounded by a lighter envelope. The core significantly enhances the capture of ultra-heavy DM; in this case, three distinct nuclear species are required to accurately describe multiple-scattering capture. We use the Eddington inversion method to obtain a realistic DM velocity distribution, better suited when the star is near the center of a halo, than the widely used Maxwell-Boltzmann distribution. We find that heavy DM would be able to thermalize and achieve capture-annihilation equilibrium within a massive star's lifetime for regions of the parameter space not excluded by direct detection. For non-annihilating DM, because of the high amount of targets available for capture and despite massive stars being short-lived, it would even be possible for DM to achieve self-gravitation and collapse to a black hole, which eventually could swallow the star from within before the expected end of the star's life, for non-excluded regions of the parameter space. Our results highlight the dependence of DM capture on the stellar evolutionary stage, composition, and halo location, demonstrating that accurate modeling of massive stars is essential for constraining heavy DM with primordial stellar populations.

dark matter theory↗
NASA NTRS1978

Microcrater and solar-flare track maturation of the lunar regolith

The microcrater, pancake and cosmic ray track record in the 100-200 microns feldspar fraction of a variety of surface soils and Apollo and Luna 16-24 cores have been studied. In Luna 24, a study of the glassy agglutinates plus breccias in 9 levels revealed this core was deposited as discrete layers, most of which are immature. Luna 16 and 20 cores were found respectively mature and immature. Glassy agglutinates and impact microcraters, absent to rare in immature soils such as 12033 and 14141, are present in similar abundances on crystals from submature to mature soils. In submature and mature soils, pancakes are present on 60 to 100 percent of crystals and craters on 50 to 70 percent. Microcrater, track and solar wind gas data suggest feldspar crystals in the size range considered reach an equilibrium maturity grade faster than bulk soils. This is interpreted as due to the limited lifetime of these crystals in the lunar regolith.

Poupeau, G.↗
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗
NASA NTRS2020

Mantle-melt Partitioning of the Highly Siderophile Elements: New Results and Application to Mars

Trace elements and extant and extinct isotopic attributes in martian meteorites have been used to argue that Mars accreted quickly, differentiated into core and mantle, and established several mantle reservoirs, possibly within 10 Ma of T0. The partitioning of trace elements in the deep mantle has been relatively unstudied, despite the need for such knowledge in understanding magma ocean crystallization and the origin of depleted and enriched mantle reservoirs. The siderophile element composition of the martian mantle, and lithophile isotopic systems such as Sr, Hf, and Nd, are thought to record evidence for early metal-silicate equilibrium and deep magma ocean at an intermediate depth and pressure of 800 km or 14 GPa. We have carried out experiments across this pressure range to better understand the mineral/melt partitioning of a wide range of elements. These new data are used to evaluate differentiation models for Mars and to help interpret the available isotopic data. The relatively incompatible nature of Re compared to mildly compatible Os means that the crystallization of a deep magma ocean will lead to residual liquids with super chondritic Re/Os, and solids with sub-chondritic Re/Os. Such material available in the mantle could be the source of enriched isotopic reservoir that produced shergottites with + Os values. On the other hand, slightly sub-chondritic Re/Os ratios in the crystallizing solids would provide a reservoir that could produce - Os values. Melting of mixtures of these two enriched and depleted end members could explain the Nd-Os isotopic correlations and systematics of shergottites.

K Righter↗
NASA NTRS1975

Binding of an adatom to a simple metal surface

The density functional formalism of Hohenberg and Kohn is used to investigate the energies, charge densities and forces which hold an adatom on the surface of a simple metal. The valence wavefunction of the adatom is fitted to the Herman-Skillman solutions at large distance and is simplified somewhat in the core region. The field of the ion is represented by the Ashcroft pseudopotential. For the metal the jellium model is used. Detailed calculations are carried out for a sodium adatom on a sodium surface. Simply juxtaposing adatom and surface gives a binding energy of about 1/3 eV. This value is approximately twice the surface energy per atom in the close-packed plane. Charge redistributions as determined variationally increase the binding energy by about 10%. The equilibrium distance for the adatom turns out to be 1.66 A from the surface, as compared with 1.52 A, the observed value for one-half the distance between the close-packed planes.

Huntington, H. B.↗
NASA NTRS2018

High-resolution SOFIA/EXES Spectroscopy of SO2 Gas in the Massive Young Stellar Object MonR2 IRS3: Implications for the Sulfur Budget

Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the ν3 rovibrational band of SO2 obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO2 gas in the massive YSO MonR2 IRS3. This yields a SO2/H abundance lower limit of 5.6 ± 0.5 × 10(exp -7), or >4% of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of b < 3.20 km s(exp -1). The small line widths and high temperature (Tex = 234 ± 15 K) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov–Gerasimenko), which is followed by SO2 formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO2 from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.

Infrared ISM↗
NASA NTRS1987

Investigation of turbulent flow in highly curved ducts with application to turbomachinery components

Numerical predictions were performed using a semi-elliptic calculation procedure for the case of turbulent flow in passage through a 90 deg bend of square cross section. Two versions of the isotropic turbulent viscosity two equation k-epsilon model were used. The first, the wall function model (WFM), employs the logarithmic law-of-the-wall relation and the notion of equilibrium flow to set all the necessary boundary conditions at the first grid point adjacent to a solid wall. The second model, the Van Driest model (VDM), employs Prandtl's original mixing length formulation, in conjunction with Van Driest's semi-empirical relation for the mixing length, to calculate the turbulent viscosity in the near wall regions of the flow. In this case, boundary conditions for k and epsilon, required to calculate these quantities in the core of the flow, are obtained by matching the mixing length and Reynolds number model formulations in an overlapping region of the flow near the walls. In both cases the results obtained show an improvement over earlier calculations using an elliptic numerical procedure. This is attributed to the finer grids possible in the present work. Of the two models, the VDM formulation shows better overall conformity with the mean flow measurements. Neither model reproduces well the details of the stress distribution as a result of the implied isotropic turbulent viscosity.

Humphrey, J. A. C.↗
NASA NTRS2006

Internal Structure of Rhea

We model the interior of Rhea on the basis of observational constraints and the results from geodynamical models available in the literature. Ten main types of models are defined, depending on the presence or absence of a high-pressure ice layer (ice II), and the extent of separation of the rock component from the volatiles. The degree-two gravity coefficients are computed for each of these models in order to assess which properties of the interior are likely to be inferred from Cassini radio science measurements scheduled on 26 November 2005. C22 greater than 2.5 x 10(exp -4) indicates that the satellite is undifferentiated, except for a slight increase in density with depth resulting from material self-compression. C22 between 1.67 x 10(exp -4)(lower bound) and 1.90 x 10(exp -4) indicates the presence of a rocky core, whose radius can be determined from the satellite's mass and ices densities, for a given temperature profile. For other values, most of the ten models cannot be distinguished from each other. However, assumptions on the density of the rock phase, presence or absence of ice II, and the degree of differentiation could allow a unique model to be determined in many cases. While the calculation presented in this work assumes that Rhea is in hydrostatic equilibrium, it is likely that Rhea' gravity field is partly affected by nonhydrostatic anomalies.

gravity science↗
NASA NTRS1974

Violation of the Vogt-Russell theorem for homogeneous nondegenerate stars

A systematic study is made of the number and types of solutions of the equilibrium equations of stellar structure, in the case of homogeneous stars of Population I over the mass range 2-1000 solar masses, with four different opacity representations. A variant of the usual fitting method permits the simultaneous investigation of convergence and tendency toward multiplicity of the solutions. Quadratic interpolation and extrapolation of Carson's new opacity tables produces a very large opacity at low temperatures that greatly affects the loose outer layers of massive stars, while leaving the cores practically unaffected. As a result, over a small mass range, well above 100 solar masses, triple solutions exist, always near an effective temperature of log Te = 4.73. A simple classification of the known exceptions to the Vogt-Russell theorem on the uniqueness of stellar structure is given.

Stothers, R.↗
NASA NTRS1992

Evidence for a high temperature differentiation in a molten earth: A preliminary appraisal

If the earth were molten during its later stages of accretion as indicated by the present understanding of planetary accretion process, the differentiation that led to the formation of the core and mantle must have occurred at high temperatures in the range of 3000-5000 K because of the effect of pressure on the temperature of melting in the interior of the earth. This calls into question the use of low-temperature laboratory measurements of partition coefficients of trace elements to make inferences about earth accretion and differentiation. The low temperature partition coefficients cannot be directly applied to high temperature fractionations because partition coefficients refer to an equilibrium specific to a temperature for a given reaction, and must change in some proportion to exp 1/RT. There are no laboratory data on partition coefficients at the high temperatures relevant to differentiation in the interior of the earth, and an attempt to estimate high temperature distribution coefficients of siderophile elements was made by considering the chemical potential of a given element at equilibrium and how this potential changes with temperature, under some specific assumptions.

Murthy, V. Rama↗
NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 13

Contents include the following: A Fast, Non-Destructive Method for Classifying Ordinary Chondrite Falls Using Density and Magnetic Susceptibility. An Update on Results from the Magnetic Properties Experiments on the Mars Exploration Rovers, Spirit and Opportunity. Measurement Protocols for In Situ Analysis of Organic Compounds at Mars and Comets. Piping Structures on Earth and Possibly Mars: Astrobiological Implications. Uranium and Lead in the Early Planetary Core Formation: New Insights Given by High Pressure and Temperature Experiments. The Mast Cameras and Mars Descent Imager (MARDI) for the 2009 Mars Science Laboratory. MGS MOC: First Views of Mars at Sub-Meter Resolution from Orbit. Analysis of Candor Chasma Interior Layered Deposits from OMEGA/MEX Spectra. Analysis of Valley Networks on Valles Marineris Plateau Using HRSC/MEX Data. Solar Abundance of Elements from Neutron-Capture Cross Sections. Preliminary Evaluation of the Secondary Ion/Accelerator Mass Spectrometer, MegaSIMS. Equilibrium Landforms in the Dry Valleys of Antarctica: Implications for Landscape Evolution and Climate Change on Mars. Continued Study of Ba Isotopic Compositions of Presolar Silicon Carbide Grains from Supernovae. Paleoenviromental Evolution of the Holden-Uzboi Area. Stability of Magnesium Sulfate Minerals in Martian Environments. Tungsten Isotopic Constraints on the Formation and Evolution of Iron Meteorite Parent Bodies. Migration of Dust Particles and Volatiles Delivery to the Inner Planets. On the Sitting of Trapped Noble Gases in Insoluble Organic Matter of Primitive Meteorites. Trapping of Xenon Upon Evaporation-Condensation of Organic Matter Under UV Irradiation: Isotopic Fractionation and Electron Paramagnetic Resonance Analysis. Stability of Water on Mars. A Didactic Activity. Analysis of Coronae in the Parga Chasma Region, Venus. Photometric and Compositional Surface Properties of the Gusev Crater Region, Mars, as Derived from Multi-Angle, Multi-Spectral Investigation of Mars Express HRSC Data. Mapping Compositional Diversity on Mars: Spatial Distribution and Geological Implications. A New Simulation Chamber for Studying Planetary Environments. Folded Structure in Terra Sirenum. Mars. Nitrogen-Noble Gas Static Mass Spectrometry of Genesis Collector Materials. Neon Isotope Heterogeneity in the Terrestrial Mantle: Implication for the Acquisition of Volatile Elements in Terrestrial Planets. The Cosmic Clock, the Cycle of Terrestrial Mass Extinctions.

Source record↗
NASA NTRS1983

The dynamics of rich clusters of galaxies. II - The Perseus cluster

The dynamics of the Perseus cluster are analyzed using self-consistent equilibrium analytical models. Using existing data in the literature plus new radial velocities reported here, composite surface density and velocity dispersion profiles are derived. These profiles have been compared with dynamical models described by Kent and Gunn (1982). The best fit suggest the presence of a significant degree of anisotropy in the velocity distribution: galaxy orbits are constrained to pass within a radius of seven cluster radii, or 1.3 deg of the cluster center. For Hubble constant = 50, a core radius of 340 kpc 11 arcmin and a mass to visual light ratio M/L = 300 are found. Using these results, X-ray observations of the intracluster medium in Perseus are reanalyzed. A previously noted discrepancy between the observed temperature of the hot gas and the cluster velocity dispersion is reduced but not eliminated. A cooling accretion flow previously deduced to exist in this cluster is shown to extend to only about one-third of the cluster radius.

Kent, S. M.↗
NASA NTRS2012

Effects of Real-Time NASA Vegetation Data on Model Forecasts of Severe Weather

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed a Greenness Vegetation Fraction (GVF) dataset, which is updated daily using swaths of Normalized Difference Vegetation Index data from the Moderate Resolution Imaging Spectroradiometer (MODIS) data aboard the NASA-EOS Aqua and Terra satellites. NASA SPoRT started generating daily real-time GVF composites at 1-km resolution over the Continental United States beginning 1 June 2010. A companion poster presentation (Bell et al.) primarily focuses on impact results in an offline configuration of the Noah land surface model (LSM) for the 2010 warm season, comparing the SPoRT/MODIS GVF dataset to the current operational monthly climatology GVF available within the National Centers for Environmental Prediction (NCEP) and Weather Research and Forecasting (WRF) models. This paper/presentation primarily focuses on individual case studies of severe weather events to determine the impacts and possible improvements by using the real-time, high-resolution SPoRT-MODIS GVFs in place of the coarser-resolution NCEP climatological GVFs in model simulations. The NASA-Unified WRF (NU-WRF) modeling system is employed to conduct the sensitivity simulations of individual events. The NU-WRF is an integrated modeling system based on the Advanced Research WRF dynamical core that is designed to represents aerosol, cloud, precipitation, and land processes at satellite-resolved scales in a coupled simulation environment. For this experiment, the coupling between the NASA Land Information System (LIS) and the WRF model is utilized to measure the impacts of the daily SPoRT/MODIS versus the monthly NCEP climatology GVFs. First, a spin-up run of the LIS is integrated for two years using the Noah LSM to ensure that the land surface fields reach an equilibrium state on the 4-km grid mesh used. Next, the spin-up LIS is run in two separate modes beginning on 1 June 2010, one continuing with the climatology GVFs while the other uses the daily SPoRT/MODIS GVFs. Finally, snapshots of the LIS land surface fields are used to initialize two different simulations of the NU-WRF, one running with climatology LIS and GVFs, and the other running with experimental LIS and NASA/SPoRT GVFs. In this paper/presentation, case study results will be highlighted in regions with significant differences in GVF between the NCEP climatology and SPoRT product during severe weather episodes.

Case, Jonathan L.↗
NASA NTRS1986

A possible deuterium anomaly: Implications of the CH3D/CH4 mixing ratios in the atmospheres of Jupiter, Saturn, and Uranus

Observations of CH3D in the atmospheres of the outer planets provide a test of the theory of deuterium fractionation equilibrium in the formation and evolution of these planets. Recent measurements of the CH3D/CH4 mixing ratios made for Saturn and Uranus are presented and intercompared with current values of Jupiter, illustrating large differences between the planets. Their implied D/H ratios are compared to D/H ratios derived from measurements of HD/H2; and, in the cases of Jupiter and Saturn, they may be incompatible. Implications of these comparisons are discussed in terms of the deuterium fractionation chemistry and possible enrichments of deuterium in the core ices of the planets.

Lutz, Barry L.↗