Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reproducible science”

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 991 records · Page 55

Simulating Atmospheric Processes in Earth System Models and Quantifying Uncertainties With Deep Learning Multi‐Member and Stochastic Parameterizations

Abstract Deep learning is a powerful tool to represent subgrid processes in climate models, but many application cases have so far used idealized settings and deterministic approaches. Here, we develop stochastic parameterizations with calibrated uncertainty quantification to learn subgrid convective and turbulent processes and surface radiative fluxes of a superparameterization embedded in an Earth System Model (ESM). We explore three methods to construct stochastic parameterizations: (a) a single Deep Neural Network (DNN) with Monte Carlo Dropout; (b) a multi‐member parameterization; and (c) a Variational Encoder Decoder with latent space perturbation. We show that the multi‐member parameterization improves the representation of convective processes, especially in the planetary boundary layer, compared to individual DNNs. The respective uncertainty quantification illustrates that methods (b) and (c) are advantageous compared to a dropout‐based DNN parameterization regarding the spread of convective processes. Hybrid simulations with our best‐performing multi‐member parameterizations remained challenging and crash within the first days. Therefore, we develop a pragmatic partial coupling strategy relying on the superparameterization for condensate emulation. Partial coupling reduces the computational efficiency of hybrid Earth‐like simulations but enables model stability over 5 months with our multi‐member parameterizations. However, our hybrid simulations exhibit biases in thermodynamic fields and differences in precipitation patterns. Despite this, the multi‐member parameterizations enable improvements in reproducing tropical extreme precipitation compared to a traditional convection parameterization. Despite these challenges, our results indicate the potential of a new generation of multi‐member machine learning parameterizations leveraging uncertainty quantification to improve the representation of stochasticity of subgrid effects.

Behrens, Gunnar [Deutsches Zentrum für Luft‐ und R↗

Neuromorphic intermediate representation: A unified instruction set for interoperable brain-inspired computing

Abstract Spiking neural networks and neuromorphic hardware platforms that simulate neuronal dynamics are getting wide attention and are being applied to many relevant problems using Machine Learning. Despite a well-established mathematical foundation for neural dynamics, there exists numerous software and hardware solutions and stacks whose variability makes it difficult to reproduce findings. Here, we establish a common reference frame for computations in digital neuromorphic systems, titled Neuromorphic Intermediate Representation (NIR). NIR defines a set of computational and composable model primitives as hybrid systems combining continuous-time dynamics and discrete events. By abstracting away assumptions around discretization and hardware constraints, NIR faithfully captures the computational model, while bridging differences between the evaluated implementation and the underlying mathematical formalism. NIR supports an unprecedented number of neuromorphic systems, which we demonstrate by reproducing three spiking neural network models of different complexity across 7 neuromorphic simulators and 4 digital hardware platforms. NIR decouples the development of neuromorphic hardware and software, enabling interoperability between platforms and improving accessibility to multiple neuromorphic technologies. We believe that NIR is a key next step in brain-inspired hardware-software co-evolution, enabling research towards the implementation of energy efficient computational principles of nervous systems. NIR is available atneuroir.org

Science & Technology - Other Topics↗

Probing the local electronic structure of isovalent Bi atoms in InP

Cross-sectional scanning tunneling microscopy (X-STM) is used to experimentally study the influence of isovalent Bi atoms on the electronic structure of InP. In this work, we map the spatial pattern of the Bi impurity state, which originates from Bi atoms down to the sixth layer below the surface, in topographic, filled-state X-STM images on the natural {110} cleavage planes. The Bi impurity state has a highly anisotropic bowtielike structure and extends over several lattice sites. These Bi-induced charge redistributions extend along the $\langle$110$\rangle$ directions, which define the bowtielike structures we observe. Local tight-binding calculations reproduce the experimentally observed spatial structure of the Bi impurity state. In addition, the influence of the Bi atoms on the electronic structure is investigated in scanning tunneling spectroscopy measurements. These measurements show that Bi induces a resonant state in the valence band, which shifts the band edge toward higher energies. Furthermore, we show that the energetic position of the Bi-induced resonance and its influence on the onset of the valence band edge depend crucially on the position of the Bi atoms relative to the cleavage plane.

74 ATOMIC AND MOLECULAR PHYSICS↗

Origins of electronic bands in the antiferromagnetic topological insulator MnBi 2 Te 4

Despite the rapid progress in understanding the first intrinsic magnetic topological insulator MnBi 2 Te 4 , its electronic structure remains a topic under debates. In this work, we perform a thorough spectroscopic investigation into the electronic structure of MnBi 2 Te 4 via laser-based angle-resolved photoemission spectroscopy. Through quantitative analysis, we estimate an upper bound of 3 meV for the gap size of the topological surface state. Furthermore, our circular dichroism measurements reveal band chiralities for both the topological surface state and quasi-2D bands, which can be well reproduced in a band hybridization model. A numerical simulation of energy-momentum dispersions based on a four-band model with an additional step potential near the surface provides a promising explanation for the origin of the quasi-2D bands. Our study represents a solid step forward in reconciling the existing controversies in the electronic structure of MnBi 2 Te 4 , and provides an important framework to understand the electronic structures of other relevant topological materials MnBi 2n Te 3n+1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Synthetic Control of Exciton Dynamics in Bioinspired Cofacial Porphyrin Dimers

Understanding how the complex interplay among excitonic interactions, vibronic couplings, and reorganization energy determines coherence-enabled transport mechanisms is a grand challenge with both foundational implications and potential payoffs for energy science. We use a combined experimental and theoretical approach to show how a modest change in structure may be used to modify the exciton delocalization, tune electronic and vibrational coherences, and alter the mechanism of exciton transfer in covalently linked cofacial Zn-porphyrin dimers (meso-beta linked AB m-β and meso-meso linked AA m-m ). While both AB m-β and AA m-m feature zinc porphyrins linked by a 1,2-phenylene bridge, differences in the interporphyrin connectivity set the lateral shift between macrocycles, reducing electronic coupling in AB m-β and resulting in a localized exciton. Pump-probe experiments show that the exciton dynamics is faster by almost an order of magnitude in the strongly coupled AA m-m dimer, and two-dimensional electronic spectroscopy (2DES) identifies a vibronic coherence that is absent in AB m-β . Theoretical studies indicate how the interchromophore interactions in these structures, and their system-bath couplings, influence excitonic delocalization and vibronic coherence-enabled rapid exciton transport dynamics. Real-time path integral calculations reproduce the exciton transfer kinetics observed experimentally and find that the linking-modulated exciton delocalization strongly enhances the contribution of vibronic coherences to the exciton transfer mechanism, and that this coherence accelerates the exciton transfer dynamics. These benchmark molecular design, 2DES, and theoretical studies provide a foundation for directed explorations of nonclassical effects on exciton dynamics in multiporphyrin assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Recent Developments and Applications of Quantum Well Infrared Photodetector Focal Plane Arrays

There are many applications that require long wavelength, large, uniform, reproducible, low cost, stable, and radiation-hard infrared (IR) focal plane arrays (FPAs). For example, the absorption lines of many gas molecules, such as ozone, water, carbon monoxide, carbon dioxide, and nitrous oxide occur in the wavelength region from 3 to 18 micron. Thus, IR imaging systems that operate in the long wavelength IR (LWIR) region (6 - 18 micron) are required in many space borne applications such as monitoring the global atmospheric temperature profiles, relative humidity profiles, cloud characteristics, and the distribution of minor constituents in the atmosphere which are being planned for future NASA Earth and planetary remote sensing systems. Due to higher radiation hardness, lower 1/f noise, and larger array size the GaAs based Quantum Well Infrared Photodetector (QWIP) FPAs are very attractive for such space borne applications compared to intrinsic narrow band gap detector arrays. In this presentation we will discuss the optimization of the detector design, material growth and processing that has culminated in realization of large format long-wavelength QWIP FPAs, portable and miniature LWIR cameras, holding forth great promise for myriad applications in 6-18 micron wavelength range in science, medicine, defense and industry. In addition, we will present some system demonstrations using broadband, two-color, and high quantum efficiency long-wavelength QWIP FPAs.

Gunapala, S. D.↗

Acquisition of a Biomedical Database of Acute Responses to Space Flight during Commercial Personal Suborbital Flights

There is currently too little reproducible data for a scientifically valid understanding of the initial responses of a diverse human population to weightlessness and other space flight factors. Astronauts on orbital space flights to date have been extremely healthy and fit, unlike the general human population. Data collection opportunities during the earliest phases of space flights to date, when the most dynamic responses may occur in response to abrupt transitions in acceleration loads, have been limited by operational restrictions on our ability to encumber the astronauts with even minimal monitoring instrumentation. The era of commercial personal suborbital space flights promises the availability of a large (perhaps hundreds per year), diverse population of potential participants with a vested interest in their own responses to space flight factors, and a number of flight providers interested in documenting and demonstrating the attractiveness and safety of the experience they are offering. Voluntary participation by even a fraction of the flying population in a uniform set of unobtrusive biomedical data collections would provide a database enabling statistical analyses of a variety of acute responses to a standardized space flight environment. This will benefit both the space life sciences discipline and the general state of human knowledge.

Charles, John B.↗

Whole-organism 3D quantitative characterization of zebrafish melanin by silver deposition micro-CT

We previously described X-ray histotomography, a high-resolution, non-destructive form of X-ray microtomography (micro-CT) imaging customized for three-dimensional (3D), digital histology, allowing quantitative, volumetric tissue and organismal phenotyping (Ding et al., 2019). Here, we have combined micro-CT with a novel application of ionic silver staining to characterize melanin distribution in whole zebrafish larvae. The resulting images enabled whole-body, computational analyses of regional melanin content and morphology. Normalized micro-CT reconstructions of silver-stained fish consistently reproduced pigment patterns seen by light microscopy, and further allowed direct quantitative comparisons of melanin content across wild-type and mutant samples, including subtle phenotypes not previously noticed. Silver staining of melanin for micro-CT provides proof-of-principle for whole-body, 3D computational phenomic analysis of a specific cell type at cellular resolution, with potential applications in other model organisms and melanocytic neoplasms. Advances such as this in whole-organism, high-resolution phenotyping provide superior context for studying the phenotypic effects of genetic, disease, and environmental variables.

59 BASIC BIOLOGICAL SCIENCES↗

Keck Infrared Transient Survey. I. Survey Description and Data Release 1

Abstract We present the Keck Infrared Transient Survey, a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roman SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. Here, we describe the first release of data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients <17 mag in a red optical band (usually ZTF r or ATLAS o bands); a volume-limited sample including all transients within redshift z < 0.01 ( D ≈ 50 Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete to z = 0.005. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events, luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction using PypeIt , which requires minimal human interaction to ensure reproducibility.

Tinyanont, S. (ORCID:0000000214814676)↗

Geometrical optics without singularities: using the ray time as the coordinate space

Geometrical optics (GO) is widely used for reduced modelling of waves in plasmas, but it fails near reflection points, where it predicts a spurious singularity of the wave amplitude. We show how to avoid this singularity by adopting a different representation of the wave equation. Instead of the physical coordinate 𝑥 and the wavevector 𝑘, we use the ray time 𝜏 as the new canonical coordinate and the ray energy ℎ as the associated canonical momentum. To derive the envelope equation in the 𝜏-representation, we construct the Weyl symbol calculus on the (𝜏,ℎ) space and show that the corresponding Weyl symbols are related to their (𝑥,𝑘) counterparts by the Airy transform. This allows us to express the coefficients in the envelope equation through the known properties of the original dispersion operator. When necessary, solutions of this equation can be mapped to the 𝑥-space using a generalised metaplectic transform. However, the field per se might not even be needed in practice. Instead, knowing the corresponding Wigner function usually suffices for linear and quasilinear calculations. As a Weyl symbol itself, the Wigner function can be mapped analytically, using the aforementioned Airy transform. We show that the standard Airy patterns that form in regions where conventional GO fails are successfully reproduced within metaplectic GO (MGO) simply by remapping the field from the 𝜏-space to the 𝑥-space. An extension to mode-converting waves is also presented. This formulation, which we call generalised MGO, can be particularly useful, for example, for reduced modelling of the O–X conversion in inhomogeneous plasma near the critical density, an effect that is important for fusion applications and also occurs in the ionosphere. Overall, MGO can replace GO for any practical purposes, because it better handles cutoffs and is similar otherwise.

plasma waves↗

Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures

Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.

Environmental Sciences & Ecology↗

Helicon edge heating in Proto-MPEX: Azimuthally resolved RF modeling

The Material Plasma Exposure eXperiment (MPEX) is a linear device under construction at ORNL. MPEX employs a variety of heating systems, including a helicon antenna, which is the subject of this work. Lacking a resonant absorption mechanism, as tokamaks have, modeling the RF power absorption in such antennas requires detailed consideration of cold collisional physics. Previous modeling of the heating in this antenna was done exclusively in 2D, assuming axisymmetry and a single azimuthal mode number m = 1⁠. In this study, we introduce a fully 3D model. By including the full 3D physics, with its self-consistent geometry-imposed azimuthal spectrum, we are able to reproduce the experimentally observed transition from edge heating to core heating.

Tierens, Wouter [Oak Ridge National Laboratory (OR↗

Correlating transcription and protein expression profiles of immune biomarkers following lipopolysaccharide exposure in lung epithelial cells

Universal and early recognition of pathogens occurs through recognition of evolutionarily conserved pathogen associated molecular patterns (PAMPs) by innate immune receptors and the consequent secretion of cytokines and chemokines. The intrinsic complexity of innate immune signaling and associated signal transduction challenges our ability to obtain physiologically relevant, reproducible and accurate data from experimental systems. One of the reasons for the discrepancy in observed data is the choice of measurement strategy. Immune signaling is regulated by the interplay between pathogen-derived molecules with host cells resulting in cellular expression changes. However, these cellular processes are often studied by the independent assessment of either the transcriptome or the proteome. Correlation between transcription and protein analysis is lacking in a variety of studies. In order to methodically evaluate the correlation between transcription and protein expression profiles associated with innate immune signaling, we measured cytokine and chemokine levels following exposure of human cells to the PAMP lipopolysaccharide (LPS) from the Gram-negative pathogen Pseudomonas aeruginosa . Expression of 84 messenger RNA (mRNA) transcripts and 69 proteins, including 35 overlapping targets, were measured in human lung epithelial cells. We evaluated 50 biological replicates to determine reproducibility of outcomes. Following pairwise normalization, 16 mRNA transcripts and 6 proteins were significantly upregulated following LPS exposure, while only five (CCL2, CSF3, CXCL5, CXCL8/IL8, and IL6) were upregulated in both transcriptomic and proteomic analysis. This lack of correlation between transcription and protein expression data may contribute to the discrepancy in the immune profiles reported in various studies. The use of multiomic assessments to achieve a systems-level understanding of immune signaling processes can result in the identification of host biomarker profiles for a variety of infectious diseases and facilitate countermeasure design and development.

59 BASIC BIOLOGICAL SCIENCES↗

Dataset for 'Ombadi et al. (2023). A warming-induced reduction in snow fraction amplifies rainfall extremes, Nature'

This package contains the main codes, sample input data and main result files to reproduce the analysis and results presented in the article: “Ombadi et al. (2023), A warming-induced reduction in snow fraction amplifies rainfall extremes, Nature”. The folder consists of the following: (1) “Raw data”: a folder that contains sample input data which is used in some of the codes for demonstration purposes. It also contains data that was not pre-processed such as Elevation data; (2) “Results”: this folder contains files of the main results presented in the paper including: “Annual-Max-Series”, “Change-rainfall-extremes”, “Change-snow-fraction”, “Warming levels_By scenario_model_year” and “Masks”. Description of these folders is detailed in the "Readme.rtf" file; (3) Python jupyter notebooks (Extract_Annual Max Series (AMS).ipynb, Elevation-dependent amplification of rainfall extremes.ipynb, Sensitivity_to_global_warming.ipynb) demonstrate the main steps of analysis. Further description of those notebooks is provided in the "Readme.rtf" file; (4) R code for extreme value analysis (Extreme_Value_Analysis.R). The sample and pre-processed dataset in "Raw data" is obtained from publicly available repositories of CMIP6 and ERA5 datasets; see Methods for more detail. This research was supported by Office of Science, Office of Biological and Environmental Research of the US Department of Energy under contract no. DE-AC02-05CH11231 for the CASCADE Scientific Focus (funded by the Regional and Global Model Analysis Program area within the Earth and Environmental Systems Modeling Program) and the iNAIADS Early Career Research Project (funded by the Environmental Systems Science program).

54 ENVIRONMENTAL SCIENCES↗

Multigenerational Independent Colony for Extraterrestrial Habitation, Autonomy, and Behavior Health (MICEHAB): An Investigation of a Long Duration, Partial Gravity, Autonomous Rodent Colony

The path from Earth to Mars requires exploration missions to be increasingly Earth-independent as the foundation is laid for a sustained human presence in the following decades. NASA pioneering of Mars will expand the boundaries of human exploration, as a sustainable presence on the surface requires humans to successfully reproduce in a partial gravity environment independent from Earth intervention. Before significant investment is made in capabilities leading to such pioneering efforts, the challenges of multigenerational mammalian reproduction in a partial gravity environment need be investigated. The Multi-generational Independent Colony for Extraterrestrial Habitation, Autonomy, and Behavior health is designed to study these challenges. The proposed concept is a conceptual, long duration, autonomous habitat designed to house rodents in a partial gravity environment with the goal of understanding the effects of partial gravity on mammalian reproduction over multiple generations and how to effectively design such a facility to operate autonomously while keeping the rodents healthy in order to achieve multiple generations. All systems are designed to feed forward directly to full-scale human missions to Mars. This paper presents the baseline design concept formulated after considering challenges in the mission and vehicle architectures such as: vehicle automation, automated crew health management/medical care, unique automated waste disposal and hygiene, handling of deceased crew members, reliable long-duration crew support systems, and radiation protection. This concept was selected from an architectural trade space considering the balance between mission science return and robotic and autonomy capabilities. The baseline design is described in detail including: transportation and facility operation constraints, artificial gravity system design, habitat design, and a full-scale mock-up demonstration of autonomous rodent care facilities. The proposed concept has the potential to integrate into existing mission architectures in order to achieve exploration objectives, and to demonstrate and mature common capabilities that enable a range of destinations and missions.

Rodgers, Erica M.↗

Many-Body Benchmark of Electronic Charge and Spin Densities for Li 1–x NiO 2

Accurate benchmarks are particularly important for highly correlated oxides as mean-field approximations often fail to describe the subtle balance of charge transfer and magnetism in these materials with an accuracy comparable to experimental needs. Here we present accurate diffusion Monte Carlo (DMC) results of the electronic charge and spin densities for the tunable highly correlated oxide Li 1–x NiO 2 for x = 0, 1/2, and 1. To enable quantitative comparisons, we introduce a robust density-partitioning scheme, extending Voronoi analysis to assign atomic charges from spatially noisy DMC densities. We then benchmark common approximations used in density functional theory (DFT). Comparison against DMC shows that r 2 SCAN delivers the most balanced performance across charge, spin, and radial density descriptors, nearly reproducing DMC results for LiNiO 2 and apical Ni sites in Li 0.5 NiO 2 . Hybrid functionals (PBE0, SCAN0) perform unexpectedly poorly, and PBE + U + V yields inconsistent trends between charge and spin densities. Therefore, the r 2 SCAN functional minimizes errors relative to DMC while capturing the variable valence of the Ni ion and also retaining the computational efficiency of DFT for large-scale simulations of the tunable structural and electronic phases of Li1−xNiO2. Our study highlights the importance of accurate benchmarking of the fundamental quantities involved in DFT to select appropriate DFT approximations in order to advance the predictive modeling of charge-transfer-driven phenomena in correlated electron systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-dimensional overdamped fluctuations of the soft perovskite lattice in CsPbBr 3

Lead halide perovskites exhibit structural instabilities and large atomic fluctuations thought to impact their optical and thermal properties, yet detailed structural and temporal correlations of their atomic motions remain poorly understood. Here, these correlations are resolved in CsPbBr 3 crystals using momentum-resolved neutron and X-ray scattering measurements as a function of temperature, complemented with first-principles simulations. We uncover a striking network of diffuse scattering rods, arising from the liquid-like damping of low-energy Br-dominated phonons, reproduced in our simulations of the anharmonic phonon self-energy. Here, these overdamped modes cover a continuum of wave vectors along the edges of the cubic Brillouin zone, corresponding to two-dimensional sheets of correlated rotations in real space, and could represent precursors to proposed two-dimensional polarons. Further, these motions directly impact the electronic gap edge states, linking soft anharmonic lattice dynamics and optoelectronic properties. These results provide insights into the highly unusual atomic dynamics of halide perovskites, relevant to further optimization of their optical and thermal properties.

36 MATERIALS SCIENCE↗

Bottom-up fabrication of 2D Rydberg exciton arrays in cuprous oxide

Solid-state platforms provide exceptional opportunities for advancing on-chip quantum technologies by enhancing interaction strengths through coupling, scalability, and robustness. Cuprous oxide (Cu 2 O) has recently emerged as a promising medium for scalable quantum technology due to its high-lying Rydberg excitonic states, akin to those in hydrogen atoms. To harness these nonlinearities for quantum applications, the confinement dimensions must match the Rydberg blockade size, which can reach several microns in Cu 2 O. Using a CMOS-compatible growth technique, this study demonstrates the bottom-up fabrication of site-selective arrays of Cu 2 O microparticles. We observed Rydberg excitons up to the principal quantum number n = 5 within these Cu 2 O arrays on a quartz substrate and analyzed the spatial variation of their spectrum across the array, showing robustness and reproducibility on a large chip. These results lay the groundwork for the deterministic growth of Cu 2 O around photonic structures, enabling substantial light-matter interaction on integrated photonic platforms and paving the way for scalable, on-chip quantum devices.

36 MATERIALS SCIENCE↗