Search NASA⌕ Search

SEARCH · Search NASA

Results for “STABILITY CALCULATIONS”

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 127 records · Page 7

Revisiting the Role of Entropy for Charge Separation in 1D Pi-Conjugated Semiconductors

Free carrier generation in organic donor/acceptor heterojunctions and redox-doped organic semiconductors is poorly understood, since assumed tight electron-hole binding conflicts with observed high free carrier yields. Cornerstone analyses that have guided the field for over 15 years predict that entropy can stabilize free charges in 2D and 3D pi-conjugated semiconductors but not in 1D systems. Here, the impact of entropy on charge generation in 1D pi-conjugated semiconductors is revisited by exploiting a greatly simplified system where enthalpy considerations alone should not allow for free charge generation. Noncontact solution-phase microwave conductivity is used to investigate the carrier density-dependent conductivity and dielectric constant in isolated chemically doped semiconducting single-walled carbon nanotubes in a low-dielectric solvent. Dopant chemical structure dramatically influences the carrier density-dependent complex conductivity, with bulky dopants facilitating carrier escape even at carrier densities below one carrier per nanotube. Three distinct numerical calculations show that entropic stabilization dramatically lowers the Gibbs energy barrier for free charge generation, explaining the high yield of free carriers, even in 1D. This renewed understanding of entropy's role in carrier generation has important implications for designing organic electronic devices-such as solar cells and thermoelectric energy harvesters-for enhanced carrier yield, conductivity, and performance.

36 MATERIALS SCIENCE↗

Computationally aided, entropy-driven synthesis of highly efficient and durable multi-elemental alloy catalysts

Multi-elemental alloy nanoparticles (MEA-NPs) hold great promise for catalyst discovery in a virtually unlimited compositional space. However, rational and controllable synthesize of these intrinsically complex structures remains a challenge. Here, we report the computationally aided, entropy-driven design and synthesis of highly efficient and durable catalyst MEA-NPs. The computational strategy includes prescreening of millions of compositions, prediction of alloy formation by density functional theory calculations, and examination of structural stability by a hybrid Monte Carlo and molecular dynamics method. Selected compositions can be efficiently and rapidly synthesized at high temperature (e.g., 1500 K, 0.5 s) with excellent thermal stability. We applied these MEA-NPs for catalytic NH 3 decomposition and observed outstanding performance due to the synergistic effect of multi-elemental mixing, their small size, and the alloy phase. We anticipate that the computationally aided rational design and rapid synthesis of MEA-NPs are broadly applicable for various catalytic reactions and will accelerate material discovery.

36 MATERIALS SCIENCE↗

Impact of Stabilizing Cations on Lithium Intercalation in Tunneled Manganese Oxide Cathodes

Stabilizing cations such as K + , Ba 2+ , and Ag + are known to provide charge neutrality and enhance structural stability in low-cost tunneled manganese dioxide (MnO 2 ) cathodes for Li ion batteries. However, a fundamental understanding of the role of these cations in the electrochemical performance of tunneled MnO 2 cathodes remains unclear, especially at low stabilizing cation concentrations. Here, we employ density functional theory (DFT + U) calculations to reveal the impact of stabilizing potassium cation (K + ) concentration on the structural stability, electronic properties, and kinetics of lithium transport in 2 x 2 tunneled manganese oxide (α-K y Mn 8 O 16 , at y = 0, 1, and 2) battery cathodes during lithium intercalation. Specifically, we provide insights into the effect of K + ions on several critical factors governing the electrochemical storage performance of tunneled MnO2 cathodes, including (a) energetically favorable Li+ host sites, (ii) Li + and electron transport capabilities, (iii) optimal intercalation pathways, crystal distortion, microstructural stability, and tunneled-to-layer phase transformation as a function of lithium content, and (iv) cell output voltage profile. Interestingly, we find that low K + concentrations (y ≤ 1) yield partially cation-deficient tunnels in the MnO 2 cathode. Such unique tunnel structures in the cathode enable (a) low kinetic barriers for Li transport, (b) excellent thermodynamic stability of the tunneled structure even at a high Li + loading (up to ~ 0.625 Li/Mn), and (c) good electronic conductivity facilitated by Jahn-Teller distortions; all of which are critical for achieving high capacity batteries with enhanced rate capability. Additionally, these results provide perspectives to design low-cost transition metal oxide cathodes for high-performance Li-ion batteries with excellent cycle life.

25 ENERGY STORAGE↗

The solubility of La hydroxide and stability of La 3+ and La hydroxyl complexes at acidic to mildly acidic pH from 25 to 250 °C

The mobility of rare earth elements (REE) in natural hydrothermal systems can be assessed using geochemical modeling, which requires reliable thermodynamic data of relevant aqueous species. In this study, we evaluate the controls of pH and temperature on La speciation and the role of hydroxyl complexes in REE transport at hydrothermal conditions. Batch-type hydrothermal solubility experiments were conducted using synthetic La hydroxide powders equilibrated in perchloric acid-based aqueous solutions at temperatures between 150 and 250 °C and starting pH of 2 to 5. The La hydroxide solubility is retrograde with temperature and displays a strong pH dependence with a decrease in La concentrations from acidic to mildly acidic pH spanning between 3 and 5 orders of magnitude (e.g. log La molality of –2.5 to –7.2 at 250 °C). Thermodynamic optimizations using GEMSFITS allow to retrieve the standard partial molal Gibbs energies for the La 3+ aqua ion and the formation constants for the La hydroxyl species (i.e., LaOH 2+ , La(OH) 2 + , La(OH) 3 0 ) between 25 and 250 °C. A comparison between the experimentally derived thermodynamic properties with the calculated values from the Helgeson-Kirkham-Flowers equation of state parameters indicates an increased divergence with temperature. Discrepancies in standard partial molal Gibbs energies range between ~ 1 – 12 kJ/mol and result in a predicted La hydroxide solubility differing by up to 3 orders of magnitude at 250 °C. Speciation calculations indicate a higher stability of La 3+ and LaOH 2+ over the other La hydroxyl species in the studied pH range of 3.4 to 6. Here, the optimized thermodynamic properties for La aqueous species have important implications for modeling the solubility of REE minerals such as monazite and the mobility of REE in hydrothermal systems.

58 GEOSCIENCES↗

Visualizing Oxidation Mechanisms in Few-Layered Black Phosphorus via In Situ Transmission Electron Microscopy

Layered two-dimensional (2D) black phosphorus (BP) exhibits novel semiconducting properties including a tunable bandgap and high electron mobility. However, the poor stability of BP in ambient environment severely limits potential for application in future electronic and optoelectronic devices. While passivation or encapsulation of BP using inert materials/polymers has emerged as a plausible solution, a detailed fundamental understanding of BP's reaction with oxygen is imperative to rationally advance its use in applications. Here, we use in situ environmental transmission electron microscopy to elucidate atomistic structural changes in mechanically exfoliated few-layered BP during exposure to varying partial pressures of oxygen. An amorphous oxide layer is seen on the actively etching BP edges, and the thickness of this layer increases with increasing oxygen partial pressure, indicating that oxidation proceeds via initial formation of amorphous P x O y species which sublime to result in the etching of the BP crystal. We observe that while few-layered BP is stable under the 80 kV electron beam (e-beam) in vacuum, the lattice oxidizes and degrades at room temperature in the presence of oxygen only in the region under the e-beam. The oxidative etch rate also increases with increasing e-beam dosage, suggesting the presence of an energy barrier for the oxidation reaction. Preferential oxidative etching along the [0 0 1] and [0 0 1] crystallographic directions is observed, in good agreement with density functional theory calculations showing favorable thermodynamic stability of the oxidized BP (0 0 $\bar{1}$) planes compared to the (1 0 0) planes. Finally, we expect the atomistic insights and fundamental understanding obtained here to aid in the development of novel approaches to integrate BP in future applications.

2D materials↗

Di- and Tetrameric Molybdenum Sulfide Clusters Activate and Stabilize Dihydrogen as Hydrides

NaY zeolite-encapsulated dimeric (Mo 2 S 4 ) and tetrameric (Mo 4 S 4 ) molybdenum sulfide clusters stabilize hydrogen as hydride binding to Mo atoms. Density functional theory (DFT) calculations and adsorption measurements suggest that stabilization of hydrogen as sulfhydryl (SH) groups, as typical for layered MoS 2 , is thermodynamically disfavored. Competitive adsorption of H 2 and ethene on Mo was probed by quantifying adsorbed CO on partly hydrogen and/or ethene covered samples with IR spectroscopy. During hydrogenation, experiment and theory suggest that Mo is covered predominately with ethene and sparsely with hydride. DFT calculations further predict that under reaction conditions, each Mo x S y cluster can activate only one H 2 , suggesting that the entire cluster (irrespective of its nuclearity) acts as one active site for hydrogenation. The nearly identical turnover frequencies (24.7 ± 3.3 mol ethane ·h -1 ·mol cluster -1 ), apparent activation energies (31-32 kJ·mol -1 ), and reaction orders (~0.5 in ethene and ~1.0 in H 2 ) show that the active sites in both clusters are catalytically indistinguishable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stability of Pt Skin Intermetallic Core Catalysts and Adsorption Properties for the Oxygen Reduction Reaction

Density functional theory calculations were used to determine the stability of metal slabs consisting of a Pt surface monolayer and intermetallic supporting layers made from combinations of six transition metal elements (Pt, Fe, Co, Ni, Cu, and Ag), as a model system for Pt skin intermetallic core nanoparticle catalysts. The stability of the slabs is largely determined by strain at the interface of the Pt skin and the subsurface intermetallic, which was described by a lattice matching parameter (r). The surface charge on the Pt skin was found to be correlated with the average electronegativity (EN) of the intermetallic core, so this average EN was used as a descriptor for how the electronic coupling (or ligand effect) affects adsorption energies. A total of 46 slabs were investigated in terms of their stability, from which 10 stable slabs were selected for further studies of adsorbate binding (OOH*, O*, and OH*) that are intermediates in the oxygen reduction reaction (ORR). The correlation between all three adsorption energies and descriptors r and EN was found. Using a linear fit between our descriptors and the calculated adsorption energies, the overpotential for the ORR was obtained as a function of r and EN, from which a volcano plot was produced. The volcano peak was found at r = 0.96 or at EN = 2.025. Interestingly, neither r nor EN was a sufficient single reactivity descriptor as the data points were well off the general trend in both linear fits; this implies that both the strain effect and the ligand effect influence the adsorption energies, although they are partly correlated. The (r, EN) target peak parameters were used to screen over 241 intermetallic combinations of transition metal elements as active ORR activity. In conclusion, this analysis identified 11 intermetallic compounds which can support a Pt skin to have a high predicted ORR activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synergistic Co-Ir/Ru Composite Electrocatalysts Impart Efficient and Durable Oxygen Evolution Catalysis in Acid

Exploring highly active and robust catalysts, which have low precious metal content, to boost the kinetically sluggish oxygen evolution reaction (OER) is a key concern for hydrogen production via proton exchange membrane water electrolysis (PEMWE). Here, in this work, rational engineering of the morphology and the local geometric ligand environment of Ir and Ru catalysts are presented by using defect-rich, lanthanum- and lithium-doped Co 3 O 4 nanofiber (LLCF) as substrate that promotes the electrocatalytic OER. Two catalysts, IrCoOx@LLCF and RuCoOx@LLCF, achieve mass activities of 1013.5 A g Ir –1 and 1911.4 A g Ru –1 in 0.1 M HClO 4 at 300 mV overpotential, respectively, which are 26 and 50 times higher than that of commercial IrO 2 and RuO 2 . Operando X-ray absorption spectroscopy unveils the reversible structure of IrCoOx during the OER and the suppression of over-oxidation of Co and Ir, giving rise to high stability. Density functional theory calculations reveal that the local geometric ligand engineering optimizes the binding of oxygenated species to the active sites, resulting in strongly enhanced catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis, Electronic Structure, and Reactivity of a Planar Four‐Coordinate, Cobalt–Imido Complex

Abstract A four‐coordinate cobalt–imido complex, ( t Bu mPNP)Co=NMes ( t Bu mPNP=modified PNP pincer ligand) has been synthesized from addition of 2,4,6‐trimethylphenylazide (Mes–N 3 ) to the corresponding dinitrogen complex. The solid‐state structure determined by X‐ray diffraction established a rare, idealized planar geometry with a Co=N bond distance of 1.716(2) Å. Magnetic measurements revealed an S =1 ground state with CAS‐SCF calculations supporting radical character on the imide nitrogen. Thermolysis of the cobalt–imido compound induced selective insertion of the imido group into a Co−P bond and yielded a three‐coordinate cobalt complex with a distorted T‐shaped geometry. Transition state analysis conducted with DFT calculations established the thermodynamic stability of the P–N coupled product and provided insight into the exclusive selectivity.

Park, Yoonsu↗

Transition Metal Nitrides as Promising Catalyst Supports for Tuning CO/H 2 Syngas Production from Electrochemical CO 2 Reduction

The electrochemical carbon dioxide reduction reaction (CO 2 RR) to produce synthesis gas (syngas) with tunable CO/H 2 ratios has been studied by supporting Pd catalysts on transition metal nitride (TMN) substrates. Combining experimental measurements and density functional theory (DFT) calculations, Pd-modified niobium nitride (Pd/NbN) is found to generate much higher CO and H 2 partial current densities and greater CO Faradaic efficiency than Pd-modified vanadium nitride (Pd/VN) and commercial Pd/C catalysts. In-situ X-ray diffraction identifies the formation of PdH in Pd/NbN and Pd/C under CO 2 RR conditions, whereas the Pd in Pd/VN is not fully transformed into the active PdH phase. DFT calculations show that the stabilized *HOCO and weakened *CO intermediates on PdH/NbN are critical to achieving higher CO 2 RR activity. This work suggests that NbN is a promising substrate to modify Pd, resulting in an enhanced electrochemical conversion of CO 2 to syngas with a potential reduction in precious metal loading.

25 ENERGY STORAGE↗

Synthesis, Electronic Structure, and Reactivity of a Planar Four–Coordinate, Cobalt–Imido Complex

A four-coordinate cobalt imido complex, ( tBu mPNP)Co=NMes ( tBu mPNP = modified PNP pincer ligand) has been synthesized from addition of 2,4,6-trimethylphenylazide (Mes–N 3 ) to the corresponding dinitrogen complex. The solid-state structure determined by X-ray diffraction established a rare, idealized planar geometry with a Co=N bond distance of 1.716(2) Å. Magnetic measurements revealed an S = 1 ground state with CAS-SCF calculations supporting radical character on the imide nitrogen. Here, thermolysis of the cobalt-imido compound induced selective insertion of the imido group into a Co–P bond and yielded a three-coordinate cobalt complex with a distorted T-shaped geometry. Transition state analysis conducted with DFT calculations established the thermodynamic stability of the P–N coupled product and provided insight into the exclusive selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tailoring the d-band center on Ru 1 Cu single-atom alloy nanotubes for boosting electrochemical non-enzymatic glucose sensing

The development of cost-effective and highly efficient electrocatalysts is critical to help electrochemical non-enzymatic sensors achieve high performance. Here, in this work, a new class of catalyst, Ru single atoms confined on Cu nanotubes as a single-atom alloy (Ru 1 Cu NTs), with a unique electronic structure and property, was developed to construct a novel electrochemical non-enzymatic glucose sensor for the first time. The Ru 1 Cu NTs with a diameter of about 24.0 nm showed a much lower oxidation potential (0.38 V) and 9.0-fold higher response (66.5 μA) current than Cu nanowires (Cu NWs, oxidation potential 0.47 V and current 7.4 μA) for glucose electrocatalysis. Moreover, as an electrochemical non-enzymatic glucose sensor, Ru1Cu NTs not only exhibited twofold higher sensitivity (54.9 μA mM −1 cm −2 ) and wider linear range (0.5–8 mM) than Cu NWs, but also showed a low detection limit (5.0 μM), excellent selectivity, and great stability. According to theoretical calculation results, the outstanding catalytic and sensing performance of Ru 1 Cu NTs could be ascribed to the upshift of the d-band center that helped promote glucose adsorption. This work presents a new avenue for developing highly active catalysts for electrochemical non-enzymatic sensors.

electrochemistry↗

Critical mineral substitutions in IN617: A combined computational and experimental approach to performance evaluation and feasibility

Addressing the escalating demand for critical minerals (CMs) driven by global climate change initiatives, this study explores compositional modifications to Inconel 617 (IN617) by substituting cobalt (Co) with manganese (Mn) across various atomic percentages. Here, we conducted a computational feasibility study employing Molecular Dynamics (MD) simulations to provide strategic guidance for experimental validation. The simulations analyzed tensile strength and corrosion resistance for five modified compositions (M1 to M5) to identify optimal properties. Tensile tests on cubic simulation cells were performed to generate stress-strain curves, revealing the impact of Co replacement with Mn on tensile strength—a metric correlated with hardness. Oxygen penetration simulations were conducted to evaluate corrosion resistance, indicating that reduced oxygen penetration depth corresponds to enhanced resistance. Promising compositions underwent phase diagram calculations for assessing phase stability. The optimal composition (M1), characterized by high tensile strength and minimal oxygen penetration, was chosen for experimental validation using induction melting and friction stir consolidation techniques. The materials were further characterized using SEM-EDS, XRD, and Vickers hardness testing. Our findings suggest that Mn substitution in IN617 can yield mechanical performance at par with high-Co alloys in energy-critical applications.

36 MATERIALS SCIENCE↗

A Six-Membered Concerted Mechanism for CO 2 Capture by Amines Studied under Charged Microdroplet Reaction Conditions

Carbon dioxide (CO 2 ) capture and storage represents an important technological challenge. A mechanistic understanding of interactions involved in the capture process is necessary not only for technological development but also for efficient conversion of captured CO 2 into value-added materials. Herein, we present a novel contained secondary electrospray ionization platform for studying the interactions of gaseous amines and CO 2 gas under microdroplet reaction conditions, which enables mass spectrometry (MS) characterization of CO 2 capture products and intermediates in real time. We detected [2 M + CO 2 + H] + species, which corresponds to a six-membered intermediate. DFT calculations confirmed the high stability of the protonated six-membered ring intermediate. This finding provides a plausible concerted mechanism in the microdroplet environment that excludes the involvement of thermodynamically disfavored ionic species. The carbamic acid counterpart of the final product/salt was readily characterized by tandem MS. The carbamic acid/amine salt was also isolated and characterized by Fourier transform infrared spectroscopy. By virtue of the fact that headspace vapors of amines are sampled, we were able to establish a high-throughput platform that enabled the CO 2 capture capacity of five different amines to be studied in under 2 min. The same device also enabled the absolute quantification of capture capacity.

Amines↗

Exploring Ca–Ce–M–O (M = 3d Transition Metal) Oxide Perovskites for Solar Thermochemical Applications

Solar thermochemical (STC) processes hold promise as efficient ways to generate renewable fuels, fuel precursors, or chemical feedstocks using concentrated sunlight. Specifically, one actively researched approach is the two-step STC cycle, which uses a redox-active, off-stoichiometric, transition-metal oxide material to split water and/or CO 2 , generating H 2 and/or CO, respectively, or syngas (a combination of H 2 and CO). Identifying novel metal oxides that yield larger reduction extents (practically achievable off-stoichiometries) than the state-of-the-art CeO 2 is critical. Here, we explore the chemical space of Ca–Ce–M–O (M = 3d transition metal, except Cu and Zn) metal oxide perovskites, with Ca and/or Ce occupying the A site and M occupying the B site within an ABO 3 framework, as potential STC candidates. We use density functional theory (DFT)-based calculations and systematically evaluate the oxygen vacancy (VaO) formation energy (≈ enthalpy of reduction in an STC cycle), electronic properties, thermodynamic stability of CaMO 3 , CeMO 3 , and Ca 0.5 Ce 0.5 MO 3 perovskites, and the VaO formation energy within Ca 0.5 Ce 0.5 Ti 0.5 Mg 0.5 O 3 perovskite. We consider only Ca and/or Ce on the A site because of their similar size and the potential redox activity of Ce 4+ . If both Ce and M exhibit simultaneous reduction with Va O formation, the resulting perovskite could exhibit a larger entropy of reduction than a single cation reduction. The increased entropy produces increased reduction for fixed temperature, partial pressure of oxygen, and reduction enthalpy, and therefore increased STC efficiency. Importantly, we identify Ca 0.5 Ce 0.5 MnO 3 , Ca 0.5 Ce 0.5 FeO 3 , and Ca 0.5 Ce 0.5 VO 3 to be promising candidates based on their Va O formation energy and thermodynamic (meta)stability. Moreover, based on our calculated on-site magnetic moments, electron density of states, and electron density differences between pristine and defective structures, we find Ca 0.5 Ce 0.5 MnO 3 to exhibit simultaneous reduction of both Ce 4+ (A-site) and Mn 3+ (B-site), highlighting a particularly promising candidate for STC applications with a predicted higher entropy of reduction than CeO 2 . Lastly, we extract metrics that govern the trends in Va O formation energies, such as standard reduction potentials, and provide pointers for further experimental and theoretical studies, which will enable the design of improved materials for the STC cycle.

14 SOLAR ENERGY↗

Butyltin Keggin Ion with a Rare Four-Coordinate Ca Center

Alkyltin clusters are exploited in nanolithography for the fabrication of microelectronics. The alkyltin Keggin family is unique among Keggin clusters across the periodic table; its members appear to favor the lower-symmetry β and γ isomers rather than the highly symmetrical α and ε isomers. Therefore, the alkyltin Keggin family may provide important fundamental information about the formation and isomerization of Keggin clusters. In this study, we have synthesized and structurally characterized a new butyltin Keggin cluster with a tetrahedral Ca 2+ center, fully formulated [(BuSn) 12 (CaO 4 )(OCH 3 ) 12 (O) 4 (OH) 8 ] 2+ (β-CaSn 12 ). The synthesis is a simple one-step process. Extensive solution characterization including electrospray ionization mass spectrometry, small-angle X-ray scattering, and multinuclear ( 1 H, 13 C, and 119 Sn) nuclear magnetic resonance shows β-CaSn 12 is essentially phase-pure and stable. This differs from the previously reported Na-centered analogues that always form a mixture of β and γ isomers, with facile interconversion. Therefore, this study has clarified prior confusion over complex spectroscopic and crystallographic characterization of the Na-centered analogues. Density functional theory calculations showed the following stability order: γ-CaSn 12 < γ-NaSn 12 < β-CaSn 12 < β-NaSn 12 . The β analogue is always more stable than the γ analogue, consistent with experiment. Notable outcomes of this study include a rare tetrahedral Ca coordination, a Na-free alkyltin cluster (important for microelectronics manufacturing), and a better understanding of Keggin families built of different metal cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Constitutive Model of Radiation Aging Effects in Filled Silicone Elastomers under Strain

Filled silicone elastomers, an essential component in many technological applications, are often subjected to controlled or unintended radiation for a variety of reasons. Radiation exposure can lead to permanent mechanical and structural changes in the material, which is manifested as altered mechanical response, and in some cases, a permanent set. For unfilled elastomers, network theories developed and refined over decades can explain these effects in terms of chain-scission and cross-link formation and a hypothesis involving independent networks formed at different strain levels of the material. In this work, we expose a filled silicone rubber to gamma radiation while being under finite elongational strain and show that the observed mechanical and structural changes can be quantitatively modeled within the same theoretical framework developed for unfilled elastomers as long as nuances associated with the Mullins effect are accounted for in a consistent manner. In this work, we employ Ogden’s incompressible hyperelastic model within the framework of Tobolsky’s two-network scheme to describe the observed permanent set and mechanical modulus changes as a function of radiation dosage. In the process, we conclude that gamma radiation induces both direct cross-linking at chain crossings (H-links) and main-chain-scission followed by cross-linking (Y-links). We provide an estimate of the ratio of chain-scission to cross-linking rates, which is in reasonable agreement with previous experimental estimate from Charlesby–Pinner analysis. We use density functional theory (DFT)-based quantum mechanical calculations to explore the stability of –Si and –SiO radicals that form upon a radiation-induced chain-scission event, which sheds light on the relative rates of Y-linking and H-linking processes.

36 MATERIALS SCIENCE↗