Search NASASearch

SEARCH · Search NASA

Results for “Oxygen from Regolith”

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.

873 records · Page 45

Comparative neutron-irradiation effects on thermal conductivity degradation and dimensional stability of TiC, TiB 2 , and ZrB 2 at 200–1000 °C

Ultra-high-temperature ceramics (UHTCs), including TiC, Ti 11 B₂, and Zr 11 B₂, show great potential for plasma-facing components due to their excellent high-temperature properties prior to irradiation. However, their response to neutron irradiation remains insufficiently understood, limiting robust assessment of their viability for fusion energy applications. Here, this study examines the thermal conductivity, dimensional stability and microstructure of TiC, TiB₂, and ZrB₂ following neutron irradiation at temperatures of 200–1000 °C and fast neutron fluences of 2.0 × 10 25 to 1.1 × 10 26 n/m 2 (E > 0.1 MeV). Lattice swelling measured by synchrotron X-ray diffraction in all three UHTCs was maximized at 200 °C and decreased with increasing irradiation temperature, with no evidence of amorphization observed at 200 °C. Above 600 °C, significant macroscopic volume swelling was observed in irradiated Ti 11 B₂ and Zr 11 B₂, but not in TiC, likely due to cavity formation in the diborides. The post-irradiation thermal conductivity, measured at the irradiation temperature, ranged from 28 to 45 W/m·K, representing a 34–45% reduction relative to the unirradiated material. Notably, neutron-irradiated UHTCs exhibit recoverable thermal conductivity at elevated temperatures, comparable to ferritic–martensitic steels and potentially superior to W when transmutation effects are considered, highlighting promise for shielding or armor plasma-facing components. At 600 °C, both thermal conductivity degradation and lattice swelling saturated at doses exceeding 2–4 dpa.

fusion materials

Validation Data for Benchmarking Wire Arc Additive Manufacturing Process Simulations

Residual stresses cause geometric distortion and affect mechanical performance of additively manufactured structures, yet they are notoriously difficult to assess and predict. Distortion (warpage) can drive parts outside dimensional tolerance limits, leading to part rejection or rework. For parts that meet tolerance, locked-in residual stress fields can affect structural integrity during operation, particularly subcritical cracking by fatigue, creep, or corrosion. This work develops benchmark data for a common additive manufacturing process (Wire Arc Additive Manufacturing) that can be applied for calibration and validation of physical process models that predict residual stress fields. The work includes design of two different samples of differing geometry, detailed manufacturing records for a set of physical samples, and an extensive set of residual stress measurement data developed using two diverse techniques (the contour method and neutron diffraction). An initial application of the work is also reported, where a modeling challenge was issued to secure residual stress model predictions from two independent laboratories that were blind to residual stress measurement data. These initial blind residual stress predictions show significant discrepancies relative to the measurement data, illustrating the potential value of the underlying validation data. An open repository for this work, including the sample designs, manufacturing process records, and the residual stress data, is also provided for future application in non-blind validation efforts.

36 MATERIALS SCIENCE

Comparing the fluorination of Pb+ and Po+: insights into p-block relativistic effects for superheavy elements

It is expected that the chemical properties of the heaviest of the superheavy elements (SHEs, = 113 – 118) do not align with what is suggested by their current positions on the periodic table. Specifically, the onset of significant relativistic effects, including increased spin–orbit splitting of the p-orbitals may lead to enhanced stability of low-oxidation states. Notably, it is predicted that flerovium ( = 114) may exhibit pseudo-noble gas behavior from a electron configuration that acts as an ‘inert pair’. Even though it is expected that the – splitting becomes pronounced toward the end of the p-block’s sixth row, there is currently limited experimental evidence to confirm its extent or impact. Here, the production of gas-phase lead ( = 82) and polonium ( = 84) fluoride cations (PbF and PoF ) were compared to elucidate differences in accessible oxidation states. The PoF and PbF species were produced and identified with the FIONA spectrometer at the Lawrence Berkeley National Laboratory 88-Inch Cylctron Facility. Polonium showed notably different fluorination chemistry than lead, producing PoF and PoF as primary products compared to lead’s PbF and PbF . The distribution of PoF products observed offer insights as to the role of spin–orbit splitting for polonium. Similar studies of superheavy elements would elucidate the accessibility of their low-oxidation states as well as to inform chemical predictions for eighth-row elements not yet discovered.

Pore, Jennifer L [Nuclear Science Division, Lawren

Dimensional Reduction Guides Electronic Structure Evolution in the A n Cu 4–n SnS 4 Semiconductor Series

The search for new functional materials with tunable properties remains a central challenge in chemistry, particularly for applications in energy and electronics. In this work, we present a framework for predictive crystal design in alkali metal chalcogenides that enables controlled dimensional reduction of a parent covalent motif, yielding a broad range of electronic structures, which systematically evolve from one parent to the other. We present 11 new members of the A n Cu 4–n SnS 4 family (A = alkali metal; n = 0–4), which reduce the three-dimensional (3D) covalent network of Cu 4 SnS 4 into various 3D, 2D, 1D, and 0D [Cu 4–n SnS 4 ] n− motifs through the substitution of Cu with alkali metals of various radii. The end members of the family set the range in achievable band gaps at 0.99 eV for fully covalent Cu 4 SnS 4 (n = 0) and 3.38 eV for K 4 SnS 4 (n = 4) with 0D [SnS 4 ] n− tetrahedra. As the dimensionality of [Cu 4–n SnS 4 ] n− systematically reduces within A n Cu 4–n SnS 4 (n = 1–3), a stepwise increase in band gap energy occurs through a gradual decrease in the energy of the valence band maximum and an increase in the conduction band minimum, with an increase in the effective masses of charge carriers. Furthermore, irrespective of the alkali metal, the thermal stability decreases with decreasing [Cu 4–n SnS 4 ] n− dimensionality within the quaternary members. Most importantly, we demonstrate that predictable crystal structure and property evolution for a given composition space is possible by deriving a general formula based on substituting the covalent metals of a parent structure with alkali metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Reduced instability growth and improved radiation trapping with optimized shock timing in double-shell inertial confinement fusion capsules

The double shell is a volume-burn inertial confinement fusion concept consisting of two concentric shells: a low-Z outer shell that collides with and transfers momentum to a high-Z inner shell which compresses and heats the thermonuclear fuel. The increased number of capsule interfaces and severe hydrodynamic instability of the high-density pusher during its acceleration phase provide challenges to the success of the double shell. Two-dimensional radiation-hydrodynamics simulations predict the hydrodynamic instability growth on the outer surface of the pusher can be greatly reduced through appropriate timing of two shocks that cross this interface. One of these shocks, unique to multi-shell designs, arises from radiation-driven ablation of the inner shell ahead of the main shock, the second shock of concern. The shock timing is optimized by increasing the thickness of a low-Z tamper layer exterior to the pusher, resulting in only minimal changes to the implosion timing. Reducing the instability growth on the outer surface of the high-Z pusher can dramatically decrease the modulations that feedthrough to the pusher inner surface, improving the efficacy of radiation trapping in the thermonuclear fuel and increasing the predicted thermonuclear yield by ≳20×.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Disorder-induced spin-cluster magnetism in a doped kagome spin liquid candidate

The search for new quantum spin liquid materials relies on systems with strong frustration such as spins on an ideal kagome lattice. However, lattice imperfections can have substantial effects which are as yet not well understood. In recent work, the two-dimensional kagome system YCu 3 ⁢(OH) 6 ⁢[(Cl 𝑥 ⁢Br (1−𝑥) ) 3−𝑦 ⁢(OH) 𝑦 ] has emerged as a leading candidate hosting a Dirac spin liquid which appears to survive at least for 𝑥 < 0.4, associated with alternating-bond-hexagon (ABH) disorder. Here in magnetic samples with 𝑥 = 0.58, 𝑦 = 0.1 we report unusual in-plane ferromagnetic canting (FM) of the in-plane antiferromagnet (AFM), with an unusually wide regime of short-ranged order, and propose theoretical models to explain this behavior. First, we show that Kitaev-type exchanges naturally arise on the kagome lattice to second order in the known Dzyaloshinskii-Moriya exchanges, and that these interactions can produce the unusual in-plane FM canting from antichiral AFM. Second, we propose a phenomenological model of weakly FM-canted spin clusters to describe the short-ranged regime and analyze quantum fluctuations in an ABH toy model to show how ABH disorder can stabilize this regime. Here, the combination of experimental observation and theory suggests that kagome-Kitaev interactions and ABH disorder are necessary for describing the magnetic fluctuations in this family of materials, with potential implications for the proposed proximate spin liquid phase.

Seth, Arnab [Georgia Institute of Technology, Atla

Robustness of pairwise kinematic Sunyaev–Zel’dovich effect to optical-cluster-selection bias

The pairwise kinematic Sunyaev–Zel’dovich (kSZ) effect measures both the pairwise motion between galaxy groups and clusters and the amount of gas within them, providing a tracer for cosmic growth. To interpret the cosmological information in the kSZ measurements, it is crucial to understand the optical-cluster-selection bias on the kSZ observables. Line-of-sight structures that contribute to both the optical observable (e.g., richness) and the cosmological signal can induce a correlation between these two quantities at a fixed cluster mass. The selection bias arising from this correlation is a key systematic effect for cosmological analyses. For cosmological observables such as cluster abundance and weak lensing, controlling this selection bias may help explain the tension between the DES-Y1 results and the Planck constraints. In order to test for a kSZ effect equivalent of such a bias, we adopted an alternative mock richness based on galaxy counts within cylindrical volumes along the line of sight. We applied the cylindrical count method to hydrodynamical simulations across a wide range of galaxy-selection criteria, assigning richness consistent with DES-Y1 to the mock clusters. When comparing optically selected clusters to mass-selected halos, we find no significant bias on pairwise kSZ signals, pairwise velocities, or optical depth within our uncertainty limits of approximately 16, 10, and 8%, respectively.

79 ASTRONOMY AND ASTROPHYSICS

Distributed Mafic Rock Resources for Carbon Mineralization in Arizona

Ex-situ carbon mineralization is a process by which CO2 is reacted with alkaline silicate minerals and rocks to produce stable carbonate materials, which can be used for other industrial processes. Arizona, U.S.A., hosts abundant surficial mafic rocks in three young volcanic fields, Geronimo-San Bernardino, San Francisco, and Springerville, and other distributed locations throughout the state. We created a Mafic Rock Resource Inventory (MRRI) that categorizes geochemical, physical, and textural characteristics of a diverse suite of surficial mafic rock samples and provide a benchmark reaction dataset parameterizing the temperature, pressure, and pH conditions best suited ex-situ mineralization in different rock types. MRRI data is publicly available online via a map-viewer. We establish two reaction condition sets, varied in temperature, pressure, and pH, where crystal-rich and glassy rocks reach maximum reaction extent and different carbonate phases are formed. Systematic ex-situ mineralization experiments on 21 diverse rock types show trends in geochemical, mineralogical, and reactivity behavior and establish maximum effective capture capacity. From this, scoria cones in three Arizona volcanic fields have a ~62 Gt effective CO₂ storage capacity with one of the fields having a ~42 Gt storage capacity in lava flows. Reactivity results have applications to alkaline mafic rock resources exposed globally, including producing additional effective storage capacity estimates and scaled commercialization of mafic rock ex-situ mineralization, should reaction extents be improved through advances in mineralization techniques. MRRI data were used to create a Direct Air Capture to Mineralization (DACM) systems model, technoeconomic analysis, and life-cycle assessment. These documents are presented as three appendices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY

SUDA: A SUrface Dust Analyser for Compositional Mapping of the Galilean Moon Europa

Abstract The Surface Dust Analyser (SUDA) is a mass spectrometer onboard the Europa Clipper mission for investigating the surface composition of the Galilean moon Europa. Atmosphereless planetary moons such as the Galilean satellites are wrapped into a ballistic dust exosphere populated by tiny samples from the moon’s surface produced by impacts of fast micrometeoroids. SUDA will measure the composition of such surface ejecta during close flybys of Europa to obtain key chemical signatures for revealing the satellite’s composition such as organic molecules and salts, history, and geological evolution. Because of their ballistic orbits, detected ejecta can be traced back to the surface with a spatial resolution roughly equal to the instantaneous altitude of the spacecraft. SUDA is a Time-Of-Flight (TOF), reflectron-type impact mass spectrometer, optimized for a high mass resolution which only weakly depends on the impact location. The instrument will measure the mass, speed, charge, elemental, molecular, and isotopic composition of impacting grains. The instrument’s small size of$268 ~\mathrm {mm} \times 250 ~\mathrm {mm} \times 171$ 268 mm × 250 mm × 171 $~\mathrm {mm}$ mm , radiation-hard design, and rather large sensitive area of 220 cm 2 matches well the challenging demands of the Clipper mission.

Astronomy & Astrophysics

Future Projection of Tropical Upper-Tropospheric Troughs and Implications for Tropical Cyclone Activity

Summertime tropical upper-tropospheric troughs (TUTTs) provide a unified framework to better understand how extratropical and tropical forcings jointly modulate basin-scale tropical cyclone (TC) activity. In this study, we examine future changes in TUTTs and their implications for TC activity. Multimodel ensemble-mean projections from 45 Coupled Model Intercomparison Project phase 6 (CMIP6) models suggest a contraction of the Pacific TUTT and an expansion of the Atlantic TUTT as the climate warms. Consistently, future changes in environment-based TC indices indicate that the large-scale conditions will become more favorable for TC genesis and intensification over the central North Pacific but less favorable over the tropical North Atlantic and Gulf of Mexico. Utilizing a TC-permitting large-ensemble dataset [i.e., the Database for Policy Decision-Making for Future Climate Change (d4PDF)] that adequately captures the observed interannual TUTT–TC relationships, we further confirm the impacts of projected TUTT changes on the TC activity in a warmer climate. In contrast, the TUTT–TC relationship is poorly represented in most CMIP6 High-Resolution Model Intercomparison Project (HighResMIP) models; such deficiencies call for caution when assessing future TC risk based on explicitly tracked TCs in these models. Additionally, CMIP6 projections show large intermodel spread in TUTT changes, implying uncertainty in projected TC activity, especially over the central-to-eastern Pacific and the North Atlantic. This intermodel spread is associated with interhemispheric sea surface temperature warming asymmetry, which leads to a meridional shift of the intertropical convergence zone (ITCZ) and the simultaneous weakening or strengthening of TUTTs in the North Pacific and North Atlantic. The potential contributions of anthropogenic aerosol forcing and oceanic circulation to this interhemispheric warming asymmetry are briefly examined.

Climate Change

Probing the atomic dynamics of ultrafast melting with femtosecond electron diffraction

Melting is an every-day phase transition that is determined by thermodynamic parameters like temperature and pressure. In contrast, ultra-fast melting is governed by the microscopic response to a rapid energy input and, thus, can reveal the strength and dynamics of atomic bonds as well as the energy flow rate to the lattice. Accurately describing these processes remains challenging and requires detailed insights into transient states encountered. Here, we present data from femtosecond electron diffraction measurements that capture the structural evolution of copper during the ultrafast solid-to-liquid phase transformations. At absorbed energy densities 2-4 times the melting threshold, melting begins at the surface slightly below the nominal melting point followed by rapid homogeneous melting throughout the volume. Molecular dynamics simulations reproduce these observations and reveal a weak electron-lattice energy transfer rate for the given experimental conditions. Both simulations and experiments show no indications of rapid lattice collapse when its temperature surpasses proposed limits of superheating, providing evidence that the inherent dynamics limits the speed of disordering in ultrafast melting of metals.

FOS: Physical sciences

Paleotribological models using preserved fossil tissue properties reveal functional significance of Eurasian mammoth dental evolution

Eurasian mammoths (Mammuthus) underwent substantial modifications in molar morphology as later-diverging species evolved progressively thinner enamel and increased enamel crest complexity. These features have been hypothesized to reduce whole-tooth wear and extend dental longevity as increasingly graze-dominated diets evolved within the lineage. This hypothesis has yet to be directly tested. Here, in this study, we developed an in-silico wear model using experimentally derived wear rates from fossil and extant proboscidean dental tissues. The models revealed that shifts in tissue topology do not affect whole-tooth wear rate, as inverse trends in lamellar frequency and enamel thickness preserve a consistent surface enamel area fraction; the determining factor of wear. Rather, topological shifts produce a wear-emergent secondary occlusal surface with greater numbers of triturating crests that create a regular, low-relief, file-like shearing pavement. These changes in occlusal architecture likely directly impacted the mastication capacity of Mammuthus dentitions, facilitating their dietary expansion to incorporate fibrous, lower-nutrient graze.

Dental wear

Discovery of magnetic-field-tunable density modulations and spin tilting in a layered altermagnet

Altermagnets recently emerged as a new class of magnetic materials, arising from specific spin crystal symmetries. They exhibit a spin-polarized electronic band structure similar to ferromagnets, yet possess zero net magnetization, promising exotic properties. Here we study a layered triangular lattice altermagnet, cobalt-intercalated NbSe 2 using scanning tunneling microscopy and spectroscopy (STM/S). Spectroscopic-imaging STM and spin-polarized STM reveals emergent 2 a 0 tri-directional charge and spin density modulations on the selenium surface. Density functional theory simulations suggest these modulations reflect the underlying cobalt superstructure. We discover that an out-of-plane magnetic field tunes the modulation amplitudes and the electronic density-of-states in a manner dependent on the field direction and strength. This behavior is attributed to the field-induced tilting of cobalt spins, which can have profound implications on the electronic properties of the altermagnet. Our results provide atomic-scale insights to uncover a magnetic-field tunable altermagnetic band structure, highlight the importance of understanding spin canting in altermagnets.

condensed-matter physics

Time-Resolved Neutron Imaging for Hydrogen Uptake in Subsurface Lithologies

Geologic hydrogen production and underground storage are increasingly important for meeting rising energy demands while providing clean-combustion advantages. However, hydrogen’s high diffusivity and propensity for leakage through porous media necessitate direct evaluation of its transport behavior in subsurface materials. Whereas X-ray microcomputed tomography (μCT) studies often employ contrast agents or surrogate gases, this study leverages neutron transmission radiography/CT to observe hydrogen migration in situ. This work represents the first demonstration of real-time neutron radiography of hydrogen migration in reservoir and caprock lithologies. Cylindrical cores of Indiana limestone, Amherst Gray sandstone, and Tumey shale were subjected to constant-pressure hydrogen charging and scanned in real time using high-resolution neutron radiography. Results indicate immediate hydrogen infiltration in sandstone and limestone, with homogeneous distribution detected throughout their pore structure. In contrast, hydrogen remained largely absent from fine-grained shale under the same pressure, except in an apparently localized fracture zone, where neutron signatures confirmed the presence of hydrogen. Subsequent neutron CT of the sandstone sample, using image subtraction against an uncharged reference, corroborated hydrogen distribution patterns. Even under lowpressure, single-phase conditions, distinct neutron imaging signatures of hydrogen were achieved. These preliminary findings underscore the potential of neutron imaging for advancing subsurface hydrogen migration research.

08 HYDROGEN

Kondo effect in ferromagnetic quantum critical CeRh 6 ⁢Ge 4

The mechanism of a pressure-induced quantum critical point in the heavy fermion ferromagnet CeRh 6 ⁢Ge 4 has attracted interest, as ferromagnetic quantum criticality in a clean itinerant Ce compound is typically avoided. The localized versus itinerant character of the 4⁢𝑓 electrons is a key aspect for understanding this behavior. We investigated the electronic structure of the 4⁢𝑓 shell in CeRh 6 ⁢Ge 4 using core-level photoelectron and x-ray absorption spectroscopy, demonstrating the hybridization of Ce 4⁢𝑓 with the conduction electrons. Linearly polarized x-ray absorption reveals a temperature-dependent linear dichroism consistent with the crystal-electric-field sequence as inferred from the static susceptibility. This dichroism cannot be described by an ionic full-multiplet model alone, but is reproduced by including the Kondo effect within a single-impurity Anderson model in the noncrossing approximation. The Kondo effect mixes higher-lying crystal-field states into a resulting multiorbital ground state with 4⁢𝑓 occupancy, 𝑛 𝑓 ∼ 0.9. Deviations at low temperatures between the measured linear dichroism and calculated dichroism suggest an orbital-dependent Kondo effect. A scenario in which there is a multiorbital ground state and orbital-dependent Kondo hybridization should be a starting point for a model of pressure-induced criticality in CeRh 6 ⁢Ge 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Large moiré superstructure of stacked incommensurate charge density waves

This article reports how two different charge density waves with slightly different wavevectors can exist in the same crystal and create moiré superstructure. Advances in heterostructure fabrication have opened new frontiers in moiré physics. Here we extend moiré engineering from artificially assembled thin flakes with mismatched lattice parameters to materials that host incommensurate orders, presenting a long-period moiré superlattice in a layered charge-density-wave compound, EuTe 4 . Using high-momentum-resolution X-ray diffraction, we found two coexisting incommensurate charge density waves with slightly mismatched in-plane wavevectors. The interaction between these two charge density waves leads to joint commensuration with the lattice and a moiré superstructure with a period of ~13.6 nm, offering key insights into the unique properties of EuTe4, such as the temperature-invariant incommensurate wavevectors and unconventional in-gap states. Owing to interlayer phase shifts, the moiré superstructure exhibits a clear thermal hysteresis, accounting for the large hysteresis in electrical resistivity and numerous metastable states. Our findings open new directions for moiré engineering based on incommensurate lattices and highlight the important role of interlayer ordering in stacked structures.

36 MATERIALS SCIENCE

Mapping the Microhardness of Al Matrix in U-7Mo/Al Dispersion Fuels at Medium and High Burn-up by in situ Nanoindentation

The mechanical degradation of the Al matrix in U-7Mo dispersion fuels was detected after irradiation in previous studies, but more details of the degradation mechanism are not clear yet to clarify the root cause. In this study, microhardness mapping in fresh and irradiated U-7Mo dispersion fuels coated with ZrN, with burn-up of 3.35 × 10 21 and 6.28 × 10 21 f/cm 3 , are measured to reveal the hardening of the Al matrix as a function of distance to fuel particles. Micro cubes were pulled by plasma focused ion beam scanning electron microscopy and tested by in situ nanoindenter with Berkovitch tip. Size effects are measured in medium and high burn-up samples. The higher burn-up specimen exhibits a more pronounced size effect when indent load is lower than 20 mN. Size effect correction models are also calculated by plots fitting, which predicts the true microhardness without size effect. Results shows that the hardness of the Al matrix near the ZrN coating has the highest value, then decreases moving away from the coating and become stable when the distance reaches ∼10 µm. The hardness starts increasing again getting closer to the next fuel particle, which makes the hardness distribution between two fuel particles to be “U” shape.

36 - MATERIALS SCIENCE