Search NASA⌕ Search

SEARCH · Search NASA

Results for “Alloy theory”

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

Machine-learning and first-principles investigation of lightweight medium-entropy alloys for hydrogen-storage applications

The transition to a low-carbon economy demands efficient and sustainable energy-storage solutions, with hydrogen emerging as a promising clean-energy carrier and with metal hydrides recognized for their hydrogen-storage capacity. Here, we leverage machine learning (ML) to predict hydrogen-to-metal (H/M) ratios and solution energy by incorporating thermodynamic parameters and local lattice distortion (LLD) as key features. Our best-performing ML model provides improvements to H/M ratios and solution energies over a broad class of medium-entripy alloys (easily extendable to multi-principal-element alloys), such as Ti–Nb-X (X = Mo, Cr, Hf, Ta, V, Zr) and Co–Ni-X (X = Al, Mg, V). Ti–Nb–Mo alloys reveal compositional effects in H-storage behavior, in particular Ti, Nb, and V enhance H-storage capacity, while Mo reduces H/M and hydrogen weight percent by 40–50 %. We attributed results in molybdenum-rich alloys to slow hydrogen kinetics, as validated by our pressure-composition-temperature (PCT) isotherm experiments on pure Ti and Ti 5 Mo 95 alloys. Density functional theory (DFT) and molecular dynamics (MD) simulations also confirm that Ti and Nb promote H diffusion, whereas Mo hinders it, highlighting the interplay between electronic structure, lattice distortions, and hydrogen uptake. Notably, our Gradient Boosting Regression model identifies LLD as a critical factor in H/M predictions. Here, to aid material selection, we present two periodic tables illustrating elemental effects on (a) H 2 wt% and (b) solution energy, derived from ML, and provide a reference for identifying alloying elements that enhance hydrogen solubility and storage.

08 HYDROGEN↗

First-principles calculation of lattice distortion, electronic structure, and bonding properties of GeTe-based and PbSe-based high-entropy chalcogenides

The massive amount of wasted heat energy from industry has pushed the development of thermoelectric (TE) materials that directly convert heat into electricity to a new level of concern. Recently, multicomponent alloys such as GeTe-based and PbSe-based high-entropy (HE) chalcogenides have attracted a great deal of attention due to their potential application as TE materials. The nature of the interatomic bonding, lattice distortion (LD), and the electronic structure in this class of materials is not fully understood. Herein, we report a comprehensive computational investigation of nine GeTe-based HE alloys with eight metallic elements (Ag, Pb, Sb, Bi, Cu, Cd, Mn, and Sn) with large supercells of 1080 atoms each; seven PbSe-based HE solid solutions: Pb0.99−ySb0.012SnySe1−2xTexSx (x = 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, with y = 0) with supercells of 1000 atoms each; and five Pb0.99−ySb0.012SnySe1−2xTexSx (y = 0.05, 0.1, 0.15, 0.2, 0.25 with x = 0.25) solid solutions with supercells of 1000 atoms each. All these HE models are theoretically investigated for the first time. The electronic structure, interatomic bonding, charge transfer, and lattice distortion (LD) are investigated by first-principles calculations based on density functional theory. Multicomponent HE alloys can cause a significant LD, which affects their mechanical, thermal, and TE properties. The calculations for the GeTe-based HE chalcogenides showed that they are semiconductors with a narrow bandgap, except for m8, which has a semi-metallic characteristic, and this makes them good candidates for TE applications. For most of these models, the Fermi level shifts upward and locates deeply in the conduction bands, resulting in the enhancement of the electrical conductivity (σ). The bonding properties showed that most bonds in m5 are more dispersed, indicating highest LD and lower lattice thermal conductivity. For PbSe-based HE solid solutions, the LD calculations showed that the models Pb0.99Sb0.012Se0.5Te0.25S0.25 and Pb0.89Sb0.012Sn0.1Se0.5Te0.25S0.25 have the higher LD, and thus a lower lattice thermal conductivity. Such investigations are in high demand since it enables us to design new HE chalcogenides for TE applications. We use the novel concept of total bond order density as a single quantum mechanical metric to characterize the internal cohesion of these HE alloys and correlate with calculated properties, especially the mechanical properties. This work provides a solid database for HE chalcogenides and a road map for many potential applications. Moreover, the computational procedure we developed can be used to design new HE chalcogenides for specific TE applications.

Physics↗

Effect of substitutional doping and disorder on the phase stability, magnetism, and half-metallicity of Heusler alloys

Spintronics is the fast growing field that will play a key role in optimizing power consumption, memory, and processing capabilities of nanoelectronic devices. Heusler alloys are potential candidates for application in spintronics due to their room temperature (RT) half-metallicity, high Curie temperature, low lattice mismatch with most substrates, and strong control on electronic density of states at Fermi level. In this work, we investigate the effect of substitutional doping and disorder on the half-metallicity, phase stability, and magnetism of Heusler alloys using density functional theory methods. Our study shows that electronic and magnetic properties of half/full-Heusler alloys can be tuned by changing electron-count through controlled variation of chemical compositions of alloying elements. We provide a detailed discussion on the effect of substitutional doping and disorder on the tunability of half-metallic nature of Co 2 MnX and NiMnX based Heusler alloys, where X represents group 13–16 and period 3–6 elements of the periodic table. Based on the idea of electron count and disorder, we predicted a possible existence of thermodynamically stable half-metallic multicomponent bismuthides, for example, (CuNi 3 )Mn 4 Bi 4 and (ZnNi 7 )Mn 8 Bi 8 , through substitution doping at Ni site by specific Cu and Zn composition in half-Heusler NiMnBi. We believe that the design guide based on electron-counts presented for half-metals will play a key role in electronic-structure engineering of novel Heusler alloys for spintronic application, which will accelerate the development and synthesis of novel materials.

36 MATERIALS SCIENCE↗

Ab Initio Simulation of Structure and Properties in Ni-Based Superalloys: Haynes282 and Inconel740

The electronic structure, interatomic bonding, and mechanical properties of two supercell models of Ni-based superalloys are calculated using ab initio density functional theory methods. The alloys, Haynes282 and Inconel740, are face-centered cubic lattices with 864 atoms and eleven elements. These multi-component alloys have very complex electronic structure, bonding and partial-charge distributions depending on the composition and strength of the local bonding environment. We employ the novel concept of total bond order density (TBOD) and its partial components (PBOD) to ascertain the internal cohesion that controls the intricate balance between the propensity of metallic bonding between Ni, Cr and Co, and the strong bonds with C and Al. We find Inconel740 has slightly stronger mechanical properties than Haynes282. Both Inconel740 and Haynes282 show ductile natures based on Poisson’s ratio. Poisson’s ratio shows marginal correlation with the TBOD. Comparison with more conventional high entropy alloys with equal components are discussed.

36 MATERIALS SCIENCE↗

Pseudoelastic deformation in Mo-based refractory multi-principal element alloys

Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal-element alloys. We present phase and property analysis of quinary (MoW) x Zr y (TaTi) 1-x-y refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular dynamics simulations. Both CALPHAD and density-functional theory analysis of phase stability indicates a Mo-W-rich region of this quinary has a stable single-phase body-centered-cubic structure. We report first quinary composition from Mo-W-Ta-Ti-Zr family of alloy with pseudo-elastic behavior, i.e., hysteresis in stress-strain. Our analysis shows that only Mo-W-rich compositions of Mo-W-Ta-Ti-Zr, i.e., Mo+W, show reproducible hysteresis in stress-strain responsible for pseudo-elastic behavior. The (MoW) 85 Zr 7.5 (TaTi) 7.5 was down-selected based on temperature-dependent phase diagram analysis and molecular dynamics simulations predicted elastic behavior that reveals twinning-assisted pseudoelastic behavior. While mostly unexplored in body-centered-cubic crystals, twinning is a fundamental deformation mechanism that competes against dislocation slip in crystalline solids. Furthermore, this alloy shows identical cyclic deformation characteristics during uniaxial < 100 > loading, i.e., the pseudoelasticity is isotropic in loading direction. Additionally, a temperature increase from 77 to 1,500 K enhances the elastic strain recovery in load-unload cycles, offering possibly control to tune the pseudoelastic behavior.

36 MATERIALS SCIENCE↗

Density Functional Theory Investigation of Oxidation Intermediates on Gold and Gold–Silver Surfaces

Gold and gold-silver alloys can be active and selective oxidation catalysts. Previous work has suggested that O 2 dissociation occurs at bimetallic step sites on gold-silver alloys, but the site responsible for the rest of the reaction steps has not been studied. As a first step in gaining insight into this issue, in this work we investigated the adsorption of oxygen and other oxidation intermediates on the (111) and (211) facets of gold-silver alloys using density functional theory. Oxygen and silver coverage effects were analyzed, and different model structures were compared. We also examined the energy barriers for the diffusion of atomic oxygen to gain insight into O migration and spillover. On (111) surfaces, O adsorption is much stronger at low O coverage (less than 0.22 ML), while on (211) surfaces O is strongly bound at both high and low O coverage. O diffusion across the step is faster than diffusion along the step. Ag stabilizes O, both when directly bound to it and when in an adjacent site. Ag also reduces repulsive O-O interactions at low O coverage. Our calculated reaction barriers for O-assisted CH 3 O dehydrogenation suggest that reaction is faster on steps than on terraces. Overall, our findings suggest that spillover of O from Ag-rich steps to Au-rich terraces does not occur, and that oxidation reactions on gold-silver alloys occurs on step sites. More specifically, oxidation likely occurs either on Ag-rich step sites or on Au-rich step sites that are adjacent to Ag-rich step sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elastic properties of the W 0.75 R e 0.25 alloy at high pressure up to 183 GPa

The high pressure equation of state for the W 0.75 R e 0.25 alloy is experimentally determined up to 183 GPa with synchrotron angle-dispersive powder x-ray diffraction in the diamond-anvil cell and to ∼925 GPa with density-functional theory. W-Re alloys are used in many industrial high-temperature applications and as a confining gasket material in high-pressure diamond-anvil cell research. The inclusion of 25 wt. % Re achieves the highest performance in terms of strength and ductility while also maintaining the body-centered-cubic (bcc) crystal structure, yet to date there has been no investigation into its elastic behavior at high pressure. We present the experimentally and theoretically determined volumetric and elastic pressure response and systematically compare these results to other W-Re alloys, finding that the bulk modulus of W-Re alloys varies nonlinearly with Re content and W 0.75 Re 0.25 becomes less incompressible than W at 85 GPa. Published by the American Physical Society 2025

Alloys↗

Machine learning prediction of the mechanical properties of refractory multicomponent alloys based on a dataset of phase and first principles simulation

In this work, a dataset including structural and mechanical properties of refractory multicomponent alloys was developed by fusing computations of phase diagram (CALPHAD) and density functional theory (DFT). The refractory multicomponent alloys, also named refractory complex concentrated alloys (CCAs) which contain 2–5 types of refractory elements were constructed based on Special Quasi-random Structure (SQS). The phase of alloys was predicted using CALPHAD and the mechanical property of alloys with stable and single body-centered cubic (BCC) at high temperature (over 1,500°C) was investigated using DFT-based simulation. As a result, a dataset with 393 refractory alloys and 12 features, including volume, melting temperature, density, energy, elastic constants, mechanical moduli, and hardness, were produced. To test the capability of the dataset on supporting machine learning (ML) study to investigate the property of CCAs, CALPHAD, and DFT calculations were compared with principal components analysis (PCA) technique and rule of mixture (ROM), respectively. It is demonstrated that the CALPHAD and DFT results are more in line with experimental observations for the alloy phase, structural and mechanical properties. Furthermore, the data were utilized to train a verity of ML models to predict the performance of certain CCAs with advanced mechanical properties, highlighting the usefulness of the dataset for ML technique on CCA property prediction.

36 MATERIALS SCIENCE↗

Mapping Thermoelectric Transport in a Multicomponent Alloy Space

Abstract Interest in high entropy alloy thermoelectric materials is predicated on achieving ultralow lattice thermal conductivity κ L through large compositional disorder. However, here it is shown that for a given mechanism, such as mass contrast phonon scattering, κ L will be minimized along the binary alloy with highest mass contrast, such that adding an intermediate mass atom to increase atomic disorder can increase thermal conductivity. Only when each component adds an independent scattering mechanism (such as adding strain fluctuation to an existing mass fluctuation) is there a benefit. In addition, both charge carriers and heat‐carrying phonons are known to experience scattering due to alloying effects, leading to a trade‐off in thermoelectric performance. Analytic transport models are applied, based on perturbation and effective medium theories, to predict how alloy scattering will affect the thermal and electronic transport across the full compositional range of several pseudo‐ternary and pseudo‐quaternary alloy systems. To do so, a multicomponent extension is demonstrated to both thermal and electronic binary alloy scattering models based on the virtual crystal approximation. Finally, it is shown that common functional forms used in computational thermodynamics can be applied to this problem to further generalize the scattering behavior that is modeled.

36 MATERIALS SCIENCE↗

Helium interaction with solutes and impurities in neutron-irradiated nanostructured ferritic alloys: A first principles study

Density functional theory calculations are performed to explore the binding between He and alloying solutes, impurities, and transmutation products expected in neutron irradiated nanostructured ferritic alloys (NFAs), here 14YWT is taken as an example. Elements that exhibit significant binding (attraction) with an interstitial He are Y (binding energy = 0.46 eV), Mg (0.32), O (0.33), Ti (0.16), and C (0.15). Those that provide significant binding to a substitutional He are O (1.44), Y (1.24), N (0.73), H (0.56), Mg (0.52), Ti (0.34), Si (0.34), C (0.33), Al (0.32), Ni (0.26), Ta (0.23), and Mn (0.16). The presence of these elements in Fe matrix could reduce the transport of He towards oxide particles, dislocations, and internal boundaries, and could promote He bubble nucleation in the matrix. For convenience, we compile existing binding energy data of He with He n and He n V (He-vacancy) clusters. Dissociation pathway analysis reveals that, in general, the most likely dissociation of a He n V cluster is by a sequential emission of individual He atoms. Furthermore, larger bubbles are more prone to dissociation than smaller ones. In addition, we estimate the binding energy (segregation energy) of He in bulk Y 2 Ti 2 O 7 (YTO) single crystal, YTO/Fe interface, and YTO particle embedded in Fe, with respect to interstitial He in Fe, from existing formation energies of He in these structures. We also compile available data of He binding with Fe self-interstitial atom (SIA), SIA clusters, and edge and screw dislocations. Note that given the absence of DFT data, the binding with SIA clusters and dislocations are gathered from simulations with empirical potentials. Finally, the data presented in this paper is important to inform multiscale simulations of He bubble accumulation.

36 MATERIALS SCIENCE↗

Theory Guided Fine‐Tune of Strain Effects in Pt Ternary Alloy via Rare Earth Templating: Achieving High Performance PEMFCs Catalysts

The sluggish kinetics and insufficient durability of platinum-based catalysts remain crucial barriers limiting proton-exchange-membrane fuel cells (PEMFCs) deployment. Here, we report a theory-guided synthesis combined with rare-earth templating to realize a previously inaccessible Pt 5 Co-like phase with tailored atomic-scale strain. Guided by density functional theory (DFT) calculations, we identified that a Pt 5 Co-like sublayer can induce a unique mild compressive strain (−1.24%) to the Pt(111) shell and an optimal *OH binding energy shift (ΔE ≈ 0.11 eV). This shift positions the alloy catalyst near the apex of the oxygen reduction reaction activity volcano. This prediction guided the synthesis of ternary alloy Pt 5 (Ce)Co@Pt multilayer nanoparticles, featuring a Ce-stabilized core, a Pt 5 Co-like sublayer, and a Pt-rich shell. This catalyst demonstrates both exceptionally high activity and durability, achieving a mass activity of 2.6 A∙mg Pt −1 in rotating disk electrode testing. In fuel cell membrane electrode assembly tests, Pt 5 (Ce)Co@Pt achieves a current density of 1.9 A∙cm −2 at 0.7 V under heavy-duty vehicle conditions. Remarkably, it maintains 1.2 A∙cm −2 after 1 80 000 AST cycles, doubling the U.S. DOE 2025 target. This work demonstrates a rational design strategy that DFT-guided strain engineering integrates with rare-earth templating to advance Pt-based catalysts for fuel cell applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Understanding Twinning and Deformation in High Entropy Alloys

A combination of high strength and high ductility has been observed in multi-principal element alloys due to twin formation attributed to low stacking fault energy (SFE). In the pursuit of low SFE alloys, a key bottleneck is the lack of understanding of the composition–SFE cor- relations that would guide tailoring SFE via alloy composition. Using density functional theory (DFT), we show that dopant radius, which have been postulated as a key descriptor for SFE in dilute alloys, does not fully explain SFE trends across different host metals. Instead, charge density is a much more central descriptor. It allows us to (1) explain contrasting SFE trends in Ni and Cu host metals due to various dopants in dilute concentrations, (2) explain the large SFE variations observed in the literature even within a given alloy composition due to the nearest neighbor environments in “model” concentrated alloys, and (3) develop a machine learning model that can be used to predict SFEs in multi-elemental alloys. This model opens a possibility to use charge density as a descriptor for predicting SFE in alloys. Furthermore, a descriptor-less machine learning (ML) model based only on charge density images extracted from density functional theory (DFT) is developed to predict stacking fault energies (SFE) in concentrated alloys. The model is based on convolutional neural networks (CNNs) as one of the promising ML techniques for dealing with complex images and data. Identification of correct descriptors is a key bottleneck to develop ML models for predicting materials properties. Often, in most ML models, textbook physical descriptors such as atomic radius, valence charge and electronegativity are used as descriptors which have limitations because these properties change in concentrated alloys when multiple elements are mixed to form a solid solution. We illustrate that, within the scope of DFT, the search for descriptors can be circumvented by electronic charge density, which is the backbone of the Kohn-Sham DFT and describes the system completely. The performance of our model is demonstrated by predicting SFE of concentrated alloys with an RMSE and R2 of 6.18 mJ/m2 and 0.87, respectively, validating the accuracy of the proposed approach.

36 MATERIALS SCIENCE↗

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Density Functional Theory Analysis of Candidate Tritium Permeation Barriers for Vanadium Alloys as a Fusion Structural Material

Vanadium alloys are being considered as a candidate structural material for the breeder blanket component of future fusion reactors due to their good thermal conductivity, low activation, and chemical compatibility with FLiBe. However, vanadium has a high hydrogen solubility, resulting in hydrogen diffusion and trapping in the vanadium structural alloy and lower than expected tritium breeding ratios. To mitigate the permeation of hydrogen into the structural alloy, a tritium permeation barrier will be coated onto the vanadium alloy. Multiple materials have been proposed as candidate tritium permeation barriers, but the extent to which each of these materials impacts hydrogen diffusion and trapping is still unknown. In this project, we perform atomic-scale calculations to narrow down the list of potential tritium permeation barrier candidates and to understand the mechanisms through which hydrogen diffusion is affected in these materials.

36 MATERIALS SCIENCE↗

Discovery of single-atom alloy catalysts for CO 2 -to-methanol reaction by density functional theory calculations

The transformations of CO 2 molecules into valuable products are of increasing interest due to the negative impact of anthropogenic CO 2 emissions on global warming. The CO 2 -to-methanol hydrogenation is an economically profitable reaction of carbon fixation, but it still steps away from widespread industrialization because of the lack of efficient and selective catalysts. Recently, single-atom alloy (SAA) catalysts have been developed to work remarkably in CO 2 hydrogenation reactions. Doping isolated single atoms into metallic catalyst can dramatically alter the catalytic performance of the host. Here, we have performed a screening discovery on Ru and 6 RuX (X = Fe, Co, Ni, Cu, Ir and Pt) SAAs using density functional theory (DFT) computations. We considered 13 possible elementary reactions in 4 possible reaction pathways on Ru and all RuX surfaces. In the computed mechanisms, we found that the formation of *H 2 COOH and *HCOO intermediates plays a critical role in determining catalysts’ activities. Doping Co and Pt isolated single atoms into Ru surface can thermodynamically and kinetically facilitate these intermediates formation processes, eventually promoting the production of methanol. The combination of weak binding and enhanced charge redistribution on RuCo and RuPt surfaces give them improved catalytic activities over pure Ru. This work will ultimately facilitate the discovery and development of SAAs for CO 2 to methanol, serving as guidance to experiments and theoreticians alike.

25 ENERGY STORAGE↗

Coupling computational thermodynamics with density-function-theory based calculations to design L12 precipitates in Fe Ni based alloys

Achieving a high-volume fraction of thermodynamically stable L12-type precipitates that are resistant to coarsening is of great importance for the development of low-cost FeNi based austenitic steels. With the aid of computational thermodynamics, this work designed two model alloys: Fe-37.4Ni-6.1Al-2.9Ti (FNAT) and Fe-45.2Ni-5.9Al-8.5Si (FNAS). Both alloys were designed to contain a similar amount of L1 2 precipitate in Fe-Ni matrix without forming other precipitates. Density-Function-Theory (DFT) calculation was coupled with computational thermodynamics to predict the critical radius at which the precipitates change shape from spherical to cuboidal. The calculation results suggest that critical radius for the FNAT alloy is about two orders of magnitude larger than that for the FeNiAlSi alloy. Phase stability and morphology of the L1 2 precipitates in these two alloys were experimentally investigated through X-ray diffraction, atom probe tomography, and scanning and transmission electron microscopy. The L1 2 precipitates in the FeNiAlSi system were found to be cuboidal and rod shaped, with much larger size than the spherical ones in the FeNiAlTi system, agreeing with the calculation results. This work suggested that coupling computational thermodynamics with DFT calculations can be reliably used to design L1 2 precipitates in FeNi based alloys.

36 MATERIALS SCIENCE↗

Heterostructural Alloy Phase Diagram for (Cd 1-x Zn x ) 3 As 2

Alloying the topological semimetal Cd 3 As 2 with Zn 3 As 2 provides a potential route for controlling the electronic properties. We predict the alloy phase diagram from first-principles calculations, considering that both end members have a crystal structure derived from the antifluorite lattice, but with different arrangements of the unoccupied cation sites. To overcome the limitations of the regular solution approximation and to include short-range order effects, we perform Monte Carlo simulations, parameterize the temperature dependence of the mixing enthalpy ΔH m , and perform thermodynamic integration of the free energy. The resulting phase diagram exhibits features that are unique to heterostructural alloy systems and provides computational predictions of solubility limits and composition ranges that are stable against spinodal decomposition.

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↗