Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional coupling”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probing plexciton emission from 2D materials on gold nanotrenches

Probing strongly coupled quasiparticle excitations at their intrinsic length scales offers unique insights into their properties and facilitates the design of devices with novel functionalities. In this work, we investigate the formation and emission characteristics of plexcitons, arising from the interaction between surface plasmons in narrow gold nanotrenches and excitons in monolayer WSe 2 . We study this strong plasmon–exciton coupling in both the far-field and the near-field. Specifically, we observe a Rabi splitting in the far-field reflection spectra of about 80 meV under ambient conditions, consistent with our theoretical modeling. Using a custom-designed near-field probe, we find that plexciton emission originates predominantly from the lower-frequency branch, which we can directly probe and map its local field distribution. We precisely determine the plexciton's spatial extension, similar to the trench width, with nanometric precision by collecting spectra at controlled probe locations. Our work opens exciting prospects for nanoscale mapping and engineering of plexcitons in complex nanostructures with potential applications in nanophotonic devices, optoelectronics, and quantum electrodynamics in nanoscale cavities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamical origins of weakly coupled relaxor behavior in Sn-doped (Ba,Ca)TiO 3 -BiScO 3

The peculiar characteristics of relaxors, viz., a frequency-dependent dielectric permittivity peak and good functional properties (dielectric, electromechanical, electrocaloric, etc.), are attributed to nanoscale regions with correlated dipoles, or polar nanoregions (PNRs). However, the exact nature of PNRs and their contribution to relaxor behavior remains debatable. In recent years, solid solutions of BaTiO3-BiMeO3 (where Me is a metal), have emerged as an interesting system with characteristics in between that of relaxors and dipole glasses. Here, we have examined the atomistic origins of weakly coupled relaxor behavior, specifically with regard to formation of PNRs, in Sn-doped (1-x)(Ba,Ca)TiO 3-x BiScO 3 using macroscopic polarization and neutron dynamic pair distribution function measurements. We show that the short-range atomic correlations observed within the PNRs dynamically fluctuate with frequencies of the order of THz. Furthermore the composition-dependent dielectric and polarization behaviors are critically influenced by the relative stability of the atomic correlations near ~1 THz, while the instantaneous atomic correlations are largely independent of x. The current results are discussed based on a model of intrinsic local modes distributed in a dielectrically soft matrix.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Reduction of electron-phonon coupling in warm dense iron

The electron-ion relaxation dynamics in warm dense iron were investigated by time-resolved x-ray absorption near-edge spectroscopy (XANES). A novel analysis, combining ab initio density functional theory (DFT) and two temperature model (TTM) simulations, was developed to calculate the x-ray absorption spectra as a function of delay time. Here we present experimental evidence of changes at the XAS L3 edge of iron that are consistent with the reduction of the electron-phonon coupling under warm dense matter conditions. The experimental results are in agreement with the model that takes both the electron ($T_e$) and the ion temperature ($T_i$) dependence of the thermophysical properties into consideration, while models where either constant electron-phonon coupling factor (G) or only $T_e$- dependent G are used do not agree with the observed relaxation dynamics of iron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Improving and testing machine learning methods for benchmarking soil carbon dynamics representation of land surface models

Representation of soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon climate feedbacks. The magnitude of this uncertainty can be reduced by accurate representation of environmental controllers of SOC stocks in ESMs. In this study, we used data of environmental factors, field SOC observations, ESM projections and machine learning approaches to identify dominant environmental controllers of SOC stocks and derive functional relationships between environmental factors and SOC stocks. Our derived functional relationships predicted SOC stocks with similar accuracy as the machine learning approach. We used the derived relationships to benchmark the coupled model intercomparison project phase six ESM representation of SOC stocks. We found divergent environmental control representation in ESMs in comparison to field observations. Representation of SOC in ESMs can be improved by including additional environmental factors and representing their functional relationships with SOC consistent with observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method

Abstract An interface-modified reproducing kernel particle method (IM-RKPM) is introduced in this work to allow for a direct model construction from image pixels of heterogeneous polycrystalline Li-ion battery microstructures. The interface-modified reproducing kernel (IM-RK) approximation is constructed through scaling of a kernel function by a regularized distance function in conjunction with strategic placement of interface node locations. This leads to RK shape functions with either weak or strong discontinuities across material interfaces, suitable for modeling various interface mechanics. With the placement of a triple junction node and distance-based scaling of kernel functions, the resulting IM-RK shape function also possesses proper discontinuities at the triple junctions. This IM-RK approximation effectively remedies the well-known Gibb’s oscillation in the smooth approximation of discontinuities. Different from the conventional meshfree approaches for interface discontinuities, this IM-RK approach is done without additional degrees of freedom associated with the enrichment functions, and it is formulated with the standard procedures in the RK shape function construction. This work focuses on identifying the accuracy and convergence properties of IM-RKPM for modeling the coupled electro-chemo-mechanical system. A linear patch test is formulated and numerically tested for the electro-chemo-mechanical coupled problem with a Butler–Volmer boundary condition representing the physical conditions in Li-ion battery microstructures. This is followed by verification of the optimal rates of convergence of IM-RKPM for solving the coupled problem with higher order solutions. The image-based modeling of Li-ion battery microstructures in the numerical examples demonstrates the applicability of the proposed method to realistic Li-ion battery materials modeling.

25 ENERGY STORAGE↗
DOE OSTI2026

Influence of cellular redox reactions on the structure and function of light harvesting and photosystems

Photosynthesis enables the conversion of one of the most abundant and free forms of energy, sunlight, into chemical bonds through the utilization of highly tailored protein complexes. These enzymes work in unison to absorb, convert, and transform light into high-energy electrons which are used for various functions important to metabolism and cellular protection. Over the last ∼50 years, photosynthetic organisms, such as cyanobacteria, have been adapted and engineered to produce valuable compounds like hydrogen and ethylene, among others. Often this is performed by removing native and/or adding in exogenous energy utilization pathways so that light energy is re-directed towards the synthesis of desired compounds. However, the interplay between primary light capture, conversion reactions, and the downstream electron utilization sinks is not fully understood. Further complicating these strategies are the plethora of compensatory mechanisms that facilitate steady electron flow and the maintenance of photosynthesis under dynamic conditions. This manifests as structural and functional plasticity of the photosynthetic machinery, often seen in modulations of oligomeric compositions or changes in protein-protein interactions and coupling with redox enzymes. Understanding these mechanisms is crucial to biotechnology applications because re-engineering electron utilization sinks has profoundly different effects on the light capture and conversion reactions of photosynthesis. Optimization requires a molecular-level understanding of the functional interrelationships between electron sinks and photosynthetic components that influence photosynthetic efficiencies to realize potential improvements in product yields. Here, we aim to highlight how perturbation of reductive reactions is revealing the functional plasticity in key components of the photosynthetic energy transduction pathway.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Revisiting Understanding of Electrochemical CO 2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory

Copper (Cu) electrodes, as the most efficacious of CO 2 reduction reaction (CO 2 RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO 2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO 2 RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO 2 RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO 2 RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at –0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Furthermore, our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

A differentiable, physics-informed ecosystem modeling and learning framework for large-scale inverse problems: demonstration with photosynthesis simulations

Photosynthesis plays an important role in carbon, nitrogen, and water cycles. Ecosystem models for photosynthesis are characterized by many parameters that are obtained from limited in situ measurements and applied to the same plant types. Previous site-by-site calibration approaches could not leverage big data and faced issues like overfitting or parameter non-uniqueness. Here we developed an end-to-end programmatically differentiable (meaning gradients of outputs to variables used in the model can be obtained efficiently and accurately) version of the photosynthesis process representation within the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) model. As a genre of physics-informed machine learning (ML), differentiable models couple physics-based formulations to neural networks (NNs) that learn parameterizations (and potentially processes) from observations, here photosynthesis rates. We first demonstrated that the framework was able to correctly recover multiple assumed parameter values concurrently using synthetic training data. Then, using a real-world dataset consisting of many different plant functional types (PFTs), we learned parameters that performed substantially better and greatly reduced biases compared to literature values. Further, the framework allowed us to gain insights at a large scale. Our results showed that the carboxylation rate at 25 °C (V c,max25 ) was more impactful than a factor representing water limitation, although tuning both was helpful in addressing biases with the default values. This framework could potentially enable substantial improvement in our capability to learn parameters and reduce biases for ecosystem modeling at large scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

42 ENGINEERING↗
DOE OSTI2025

Coupling Interfacial Redox‐Reactions with In Situ Proton Generation for the Photoelectrochemical Separation of Rare‐Earth Elements

Abstract A dual‐function photoelectrochemical (PEC) separation system is demonstrated for rare‐earth element (REE) recovery. The sustainable release of the captured REEs is promoted through the synergistic integration of a redox‐reaction for electrostatic repulsion, and in situ proton generation for ion‐exchange, all driven by photoelectrochemistry. The platform consists of a redox‐copolymer, poly(ferrocenylpropyl methacrylamide‐ co ‐methacrylic acid) (P(FPMAm‐ co ‐MAA)) (PFM), conjugated with carbon nanotubes (CNTs) and coated onto titanium dioxide nanorods (TNRs). The (PFM‐CNT)/TNR spontaneously adsorbs up to 214.2 mg of Yttrium/g PFM by ion‐exchange, and demonstrates broad applicability for other REEs. The adsorbed REEs are released through the PEC oxidation of ferrocene (Fc) to ferrocenium (Fc + ), and the simultaneous PEC water splitting reactions at the TNRs that protonate the carboxylate binding groups. This dual photoelectrochemically‐driven mechanism for REE release is investigated by in situ pH measurements, as well as vibrational and X‐ray photoelectron spectroscopy. Through PEC approaches, a 68.8% reduction in energy consumption during REE recovery has been achieved compared to purely electrochemical systems, with a regeneration efficiency close to 100%. For NdFeB magnets from waste hard disk drives, Nd and Dy recovery efficiencies of 59.2 and 61.1% are achieved. The dual‐functionality of these copolymer PEC systems offers a sustainable platform for modulating critical element recovery.

Cho, Ki‐Hyun [Department of Chemical and Biomolecu↗
DOE OSTI2021

Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities

Here, we present a framework for calibration of parameters in elastoplastic constitutive models that is based on the use of automatic differentiation (AD). The model calibration problem is posed as a partial differential equation-constrained optimization problem where a finite element (FE) model of the coupled equilibrium equation and constitutive model evolution equations serves as the constraint. The objective function quantifies the mismatch between the displacement predicted by the FE model and full-field digital image correlation data, and the optimization problem is solved using gradient-based optimization algorithms. Forward and adjoint sensitivities are used to compute the gradient at considerably less cost than its calculation from finite difference approximations. Through the use of AD, we need only to write the constraints in terms of AD objects, where all of the derivatives required for the forward and inverse problems are obtained by appropriately seeding and evaluating these quantities. We present three numerical examples that verify the correctness of the gradient, demonstrate the AD approach's parallel computation capabilities via application to a large-scale FE model, and highlight the formulation's ease of extensibility to other classes of constitutive models.

42 ENGINEERING↗
DOE OSTI2023

Search for light long-lived neutral particles that decay to collimated pairs of leptons or light hadrons in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for light long-lived neutral particles with masses in the O(MeV–GeV) range is presented. The analysis targets the production of long-lived dark photons in the decay of a Higgs boson produced via gluon–gluon fusion or in association with a W boson. Events that contain displaced collimated Standard Model fermions reconstructed in the calorimeter or muon spectrometer are selected in 139 fb -1 of $\sqrt{s}$ = 13 TeV $pp$ collision data collected by the ATLAS detector at the LHC. Background estimates for contributions from Standard Model processes and instrumental effects are extracted from data. The observed event yields are consistent with the expected background. Exclusion limits are reported on the production cross-section times branching fraction as a function of the mean proper decay length cτ of the dark photon, or as a function of the dark-photon mass and kinetic mixing parameter that quantifies the coupling between the Standard Model and potential hidden (dark) sectors. A Higgs boson branching fraction above 1% is excluded at 95% CL for a Higgs boson decaying into two dark photons for dark-photon mean proper decay lengths between 10 mm and 250 mm and dark photons with masses between 0.4 GeV and 2 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Search for an axion-like particle with forward proton scattering in association with photon pairs at ATLAS

A search for forward proton scattering in association with light-by-light scattering mediated by an axion-like particle is presented, using the ATLAS Forward Proton spectrometer to detect scattered protons and the central ATLAS detector to detect pairs of outgoing photons. Proton-proton collision data recorded in 2017 at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV were analysed, corresponding to an integrated luminosity of 14.6 fb –1 . A total of 441 candidate events were selected. A search was made for a narrow resonance in the diphoton mass distribution, corresponding to an axion-like particle (ALP) with mass in the range 150–1600 GeV. No excess is observed above a smooth background. Upper limits on the production cross section of a narrow resonance are set as a function of the mass, and are interpreted as upper limits on the ALP production coupling constant, assuming 100% decay branching ratio into a photon pair. The inferred upper limit on the coupling constant is in the range 0.04–0.09 TeV –1 at 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Effective entropy of quantum fields coupled with gravity

Entanglement entropy, or von Neumann entropy, quantifies the amount of uncertainty of a quantum state. For quantum fields in curved space, entanglement entropy of the quantum field theory degrees of freedom is well-defined for a fixed background geometry. In this paper, we propose a generalization of the quantum field theory entanglement entropy by including dynamical gravity. The generalized quantity named effective entropy, and its Renyi entropy generalizations, are defined by analytic continuation of a replica calculation. The replicated theory is defined as a gravitational path integral with multiple copies of the original boundary conditions, with a co-dimension-2 brane at the boundary of region we are studying. We discuss different approaches to define the region in a gauge invariant way, and show that the effective entropy satisfies the quantum extremal surface formula. When the quantum fields carry a significant amount of entanglement, the quantum extremal surface can have a topology transition, after which an entanglement island region appears. Our result generalizes the Hubeny-Rangamani-Takayanagi formula of holographic entropy (with quantum corrections) to general geometries without asymptotic AdS boundary, and provides a more solid framework for addressing problems such as the Page curve of evaporating black holes in asymptotic flat spacetime. We apply the formula to two example systems, a closed two-dimensional universe and a four-dimensional maximally extended Schwarzchild black hole. We discuss the analog of the effective entropy in random tensor network models, which provides more concrete understanding of quantum information properties in general dynamical geometries. We show that, in absence of a large boundary like in AdS space case, it is essential to introduce ancilla that couples to the original system, in order for correctly characterizing quantum states and correlation functions in the random tensor network. Using the superdensity operator formalism, we study the system with ancilla and show how quantum information in the entanglement island can be reconstructed in a state-dependent and observer-dependent map. We study the closed universe (without spatial boundary) case and discuss how it is related to open universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Quantifying microbial control of soil organic matter dynamics at macrosystem scales

Soil organic matter (SOM) stocks, decomposition and persistence are largely the product of controls that act locally. Yet the controls are shaped and interact at multiple spatiotemporal scales, from which macrosystem patterns in SOM emerge. Theory on SOM turnover recognizes the resulting spatial and temporal conditionality in the effect sizes of controls that play out across macrosystems, and couples them through evolutionary and community assembly processes. For example, climate history shapes plant functional traits, which in turn interact with contemporary climate to influence SOM dynamics. Selection and assembly also shape the functional traits of soil decomposer communities, but it is less clear how in turn these traits influence temporal macrosystem patterns in SOM turnover. Here, we review evidence that establishes the expectation that selection and assembly should generate decomposer communities across macrosystems that have distinct functional effects on SOM dynamics. Representation of this knowledge in soil biogeochemical models affects the magnitude and direction of projected SOM responses under global change. Yet there is high uncertainty and low confidence in these projections. To address these issues, we make the case that a coordinated set of empirical practices are required which necessitate (1) greater use of statistical approaches in biogeochemistry that are suited to causative inference; (2) long-term, macrosystem-scale, observational and experimental networks to reveal conditionality in effect sizes, and embedded correlation, in controls on SOM turnover; and (3) use of multiple measurement grains to capture local- and macroscale variation in controls and outcomes, to avoid obscuring causative understanding through data aggregation. Here, when employed together, along with process-based models to synthesize knowledge and guide further empirical work, we believe these practices will rapidly advance understanding of microbial controls on SOM and improve carbon cycle projections that guide policies on climate adaptation and mitigation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Quasiparticle description of angle-resolved photoemission spectroscopy for SrCuO 2

SrCuO 2 has long been considered a near-archetypal realization of a quasi-one-dimensional (1D) system of interacting electrons with short-range interactions. Within this framework, experimental observations - interpreted through the lens of the 1D Hubbard model - suggest that electron and hole excitations decay into two types of (unphysical) collective bosonic modes: "spinons", which carry the spin degree of freedom, and "holons", which carry the charge degree of freedom. This model, known as spin-charge separation, is most directly evidenced by angle-resolved photoemission spectroscopy (ARPES), where a photo-induced hole decays into a continuum of these excitations. Here we present an alternative perspective grounded in first-principles, self-consistent, and parameter-free many-body perturbation theory. In this revised quasiparticle framework, ARPES can be understood as a one-body effect arising from mild disorder in a long-range antiferromagnetic ground state. the emergence of the so-called spinon branch arises naturally from spin disorder, the anomalous linewidths are accurately captured, and we provide a compelling explanation for the spectral weight observed at the non-magnetic zone boundary. This reinterpretation provides a unified explanation for key experimental signatures previously attributed to spin-charge separation, including features observed in optical conductivity. Additionally, we show that SrCuO 2 exhibits a nontrivial interchain coupling that significantly influences both its one-particle and two-particle spectral functions. By comparing the spectral features of SrCuO 2 with those of La 2 CuO 4 , we argue that SrCuO 2 shares notable similarities with the two-dimensional cuprates-both being rooted in a common CuO 4 plaquette-based molecular orbital framework.

1D cuprates↗
DOE OSTI2022

Magnetic phase diagram of A 2 [FeCl 5 (H 2 O)] ( A = K, Rb, NH 4 )

Here, erythrosiderites with the formula A 2 FeX 5 ∙H 2 O, where A = Rb, K, and (NH 4 ) and X = Cl and Br are intriguing systems that possess various magnetic and electric phases, as well as multiferroic phases in which magnetism and ferroelectricity are coupled. In this report, we study the magnetic phase diagram of erythrosiderites as a function of superexchange interactions and magnetic anisotropies. To this end, we perform classical Monte Carlo simulations on magnetic Hamiltonians that contain five different superexchange interactions with single-ion anisotropies. Our phase diagram contains all magnetic ground states that have been experimentally observed in these materials. We argue that the ground states can be explained by varying the ratio of $\frac{J_4}{J_2}$. For $\frac{J_4}{J_2}$ > 0.95 a cycloidal spins structure is stabilized as observed in (NH 4 ) 2 FeCl 5 ∙H 2 O and otherwise a collinear spin structure is stabilized as observed in (K,Rb) 2 FeCl 5 ∙H 2 O. We also show that the difference in the single-ion anisotropy along a- and c- axes is essential to stabilize the intermediate state observed in (NH) 2 FeCl 5 ∙H 2 O.

36 MATERIALS SCIENCE↗