Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MO 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Basal Plane Doping to Activate Colloidal MoS 2 Nanosheets for Catalytic Hydrodeoxygenation of para- Cresol

The valorization of biomass into biofuels is a critical process for producing renewable fuels. Hydrodeoxygenation (HDO), particularly over doped molybdenum disulfide (MoS 2 ), a transition metal dichalcogenide (TMD) material, is a common representative catalytic reaction system for converting biomass-derived materials into useful hydrocarbons. However, the location and role of dopants, such as Co, in HDO is not fully understood. The effects of dopant location and oxidation state are often precluded by inhomogeneity in the ensemble properties of nanosheet size and dopant dispersion, as well as difficulty in observing the behavior of atomic site behavior directly. Using a colloidal approach to synthesize cobalt-doped MoS 2 nanosheets with controlled dopant concentration, combined with X-ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations, we determine that basal plane doped Co (25% Co:Mo mole ratio) shows peak catalytic activity in HDO of para-cresol, a model biomass-derived compound, and that basal Co sites are demonstrably more active than edge sites. By observing these doping effects in MoS 2 catalysts for HDO, we can further optimize not only the production of carbon-neutral fuels but also direct the tailoring of doped TMD catalysts toward their intended applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomic Molybdenum for Synthesis of Ammonia with 50% Faradic Efficiency

The electrochemical dinitrogen (N 2 ) reduction reaction (NRR) under ambient conditions has gained significant interest as an environmentally friendly alternative to the traditional Haber–Bosch process for the synthesis of ammonia (NH 3 ). However, up to now, most of the reported NRR electrocatalysts with satisfactory catalytic activities have been hindered by the large overpotential in N 2 activation. The preparation of highly efficient Mo-based NRR electrocatalyst in acidic electrolytes under ambient conditions is demonstrated here, consisting of stabilized single Mo atoms anchored on holey nitrogen-doped graphene synthesized through a convenient potassium-salt-assisted activation method. At -0.05 V versus a reversible hydrogen electrode (RHE), an electrode consisting of the resultant electrocatalyst immobilized on carbon fiber paper can attain an exceptional Faradaic efficiency of 50.2% and a NH3 yield rate of 3.6 µg h -1 mg cat -1 with low overpotentials. Density functional theory calculations further unveil that compared to the original graphene without holes, the edge coordinated Mo atoms and the existence of vacancies on holey graphene lower the overpotential of N 2 reduction, thereby promoting the NRR catalytic activity. Finally, this work could provide new guidelines for future designs in single-atom catalysis that would be beneficial to ambient N 2 fixation, and replacement of classical synthesis processes that are very energy-intensive.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Prediction of the structures and heats of formation of MO 2 , MO 3 , and M 2 O 5 for M = V, Nb, Ta, Pa

Structures for the mono-, di-, and tri-bridge isomers of M 2 O 5 as well as those for the MO 2 and MO 3 fragments for M = V, Nb, Ta, and Pa were optimized at the density functional theory (DFT) level. Single point CCSD(T) calculations extrapolated to the complete basis set (CBS) limit at the DFT geometries were used to predict the energetics. Here, the lowest energy dimer isomer was the di-bridge for M = V and Nb and the tri-bridge for M = Ta and Pa. The di-bridge isomers were predicted to be composed of MO 2 + and MO 3 - fragments, whereas the mono- and tri-bridge are two MO 2 + fragments linked by an O 2 - . The heats of formation of M 2 O 5 dimers, as well as MO 2 and MO 3 neutral and ionic species were predicted using the Feller–Peterson–Dixon (FPD) approach. The heats of formation of the MF 5 species were calculated to provide additional benchmarks. Dimerization energies to form the M 2 O 5 dimers are predicted to become more negative going down group 5 and range from -29 to -45 kcal mol -1 . The ionization energies (IEs) for VO 2 and TaO 2 are essentially the same at 8.75 eV whereas the IEs for NbO 2 and PaO 2 are 8.10 and 6.25 eV, respectively. The predicted adiabatic electron affinities (AEAs) range from 3.75 eV to 4.45 eV for the MO 3 species and vertical detachment energies from 4.21 to 4.59 eV for MO 3 - . The calculated M[double bond, length as m-dash]O bond dissociation energies increase from 143 kcal mol -1 for M = V to ~170 kcal mol -1 for M = Nb and Ta to ~200 kcal mol -1 for M = Pa. The M–O bond dissociation energies are all similar ranging from 97 to 107 kcal mol -1 . Natural bond analysis provided insights into the types of chemical bonds in terms of their ionic character. Pa 2 O 5 is predicted to behave like an actinyl species dominated by the interactions of approximately linear PaO 2 + groups.

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

Vacancy formation energies and migration barriers in multi-principal element alloys

Multi-principal element alloys (MPEAs) continue to garner interest as structural and plasma-facing materials due to their structural (phase) stability and increased resistance to radiation damage. Despite sensitivity of mechanical behavior to irradiation and point-defect formation, there has been scant attention on understanding vacancy stability and diffusion in refractory-based MPEAs. Using density-functional theory, we examine vacancy stability and diffusion barriers in body-centered cubic (Mo 0.95 W 0.05 ) 0.85 Ta 0.10 (TiZr) 0.05 . The results in this MPEA show strong dependence on environment, originating from local lattice distortion associated with charge-transfer between neighboring atoms that vary with different chemical environments. We find a correlation between degree of lattice distortion and migration barrier: (Ti, Zr) with less distortion have lower barriers, while (Mo, W) with larger distortion have higher barriers, depending up local environments. Under irradiation, our findings suggest that (Ti, Zr) are significantly more likely to diffuse than (Mo, W) while Ta shows intermediate effect. Finally, material degradation caused by vacancy diffusion can be controlled by tuning composition of alloying elements to enhance creep strength at extreme operating temperatures and harsh conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-Performance Iridium–Molybdenum Oxide Electrocatalysts for Water Oxidation in Acid: Bayesian Optimization Discovery and Experimental Testing

Ir oxides are costly and scarce catalysts for oxygen evolution reaction (OER) in acid. There has been extensive interest in developing alternatives that are either Ir-free or require smaller amounts of Ir to drive the reactions at acceptable rates. One design strategy is to identify Ir-based mixed oxides that achieve similar performance while requiring smaller amounts of Ir. The obstacle to this strategy has been a very large phase space of the Ir-based mixed metal oxides, in terms of the metals combined with Ir and the different crystallographic structures of the mixed oxides, which prevents a thorough exploration of possible materials. In this work, we developed a workflow that uses machine-learning-aided Bayesian optimization in combination with density functional theory to make the exploration of this phase space plausible. This screening identified Mo as a promising dopant for forming acid-tolerant Ir-based oxides for the OER. We synthesized and characterized the Ir–Mo mixed oxides in the form of thin-film electrocatalysts with a known surface area. We show that these mixed oxides exhibited overpotentials ~30 mV lower than a pure Ir control while maintaining 24% lower Ir dissolution rates than the Ir control. Furthermore, these findings suggest that Mo is a promising dopant and highlight the promise of machine learning to guide the experimental exploration and optimization of catalytic materials.

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

Optimal parameterizations for observational constraints on thawing dark energy

Time-varying dark energy is often modeled in observational analyses through generic parameterizations of its equation of state w(z), which typically use two free parameters {w 0 , w a } to span a broad range of behaviors as a function of redshift. However, this broad range of behaviors can only approximately capture the dynamics of any given microphysical theory of dark energy. A complementary approach is to use targeted parameterizations designed to model specific classes of dynamical dark energy with greater precision. Focusing on the class of thawing dark energy, we quantify and compare the precision with which nineteen generic and targeted parameterizations can capture the dynamics of physically motivated thawing quintessence theories. We find that a targeted parameterization derived from a Padé expansion of w is the most reliable of these, producing accurate reconstructions of w(z), the expansion history H(z), and cosmological parameters such as H 0 and Ω m for a broad range of microphysical theories.

dark energy experiments↗
DOE OSTI2022

Xenon mobility in γ-uranium and uranium–molybdenum alloys

Diffusion in bcc uranium and U–Mo alloys is of great interest because fission gas and other fission products impact the performance of nuclear fuels. We investigate the mobility of xenon and molybdenum in bcc uranium (⁠$γ$-U) and metallic U–Mo alloys by calculating the migration energies of xenon and molybdenum for various local compositions using density functional theory. We also calculate the solute–vacancy binding energies of different solutes to vacancies in bcc uranium. We find that the solute–vacancy binding energy in bcc uranium is significantly higher than it is in other bcc metals (e.g., Fe and W). We also find that the migration energy of molybdenum is substantially higher than the migration energy of xenon, indicating that xenon is much more mobile than molybdenum in bcc uranium. The presence of molybdenum in the nearest-neighbor shell around a xenon atom typically increases the migration energy of xenon, which indicates a reduction of xenon mobility in U–Mo alloys compared to pure bcc uranium.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Manipulating the Geometric and Electronic Structures of Manganese Nitrides for Ammonia Synthesis

Abstract Manganese (Mn) nitrides are important nitrogen (N) carriers for small‐scale intermittent ammonia (NH 3 ) synthesis. However, only 3∼8 % of lattice N are converted into NH 3 . In this study, the geometric and electronic structures of well‐defined Mn 4 N and Mn 2 N lattices were altered using transition metal heteroatoms (Cr, Fe, Co, Ni, Mo) to understand the driving force behind lattice N diffusion and extraction. Density Functional Theory (DFT) revealed that the binding of early hydrogenation product (NH) follows a linear relationship with lattice N over a wide range of close‐packed surfaces. But, the binding of NH 2 and NH 3 are more sensitive to the geometric and electronic structures. Further, the chemical bonding can be quantitatively characterized with the covalency derived from Crystal Orbital Hamiltonian Population (COHP). In the Eley Rideal‐Mars van Krevelan pathway, the overall NH 3 formation free energy ( ) and lattice N diffusion barrier ( ) are the respective thermodynamic and kinetic determining factors. Aided by a rate‐determining step (RDS) model, Mn 4 N modified with Fe, Co, Ni single‐atom dopants all show enhanced rates for NH 3 formation.

Shan, Nannan↗
DOE OSTI2024

An active, stable cubic molybdenum carbide catalyst for the high-temperature reverse water-gas shift reaction

Although technologically promising, the reduction of carbon dioxide (CO 2 ) to produce carbon monoxide (CO) remains economically challenging owing to the lack of an inexpensive, active, highly selective, and stable catalyst. We show that nanocrystalline cubic molybdenum carbide (α-Mo 2 C), prepared through a facile and scalable route, offers 100% selectivity for CO 2 reduction to CO while maintaining its initial equilibrium conversion at high space velocity after more than 500 hours of exposure to harsh reaction conditions at 600°C. The combination of operando and postreaction characterization of the catalyst revealed that its high activity, selectivity, and stability are attributable to crystallographic phase purity, weak CO-Mo 2 C interactions, and interstitial oxygen atoms, respectively. Mechanistic studies and density functional theory (DFT) calculations provided evidence that the reaction proceeds through an H 2 -aided redox mechanism.

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

The striking influence of oxophilicity differences in heterometallic Mo–Mn oxide cluster reactions with water

Mixed-metal oxides have proven to be effective catalysts for the hydrogen evolution reaction, often outperforming either of the binary metal oxides. The reactivity of Mn x MoO y - (x = 1, 2; y = 3, 4) clusters toward H 2 O was investigated via time-of-flight mass spectrometry with clear evidence of cluster oxidation and corresponding H 2 production, specifically for Mn x MoO 3 - (x = 1, 2) clusters. Unlike previously studied MoxOy- clusters, which assumed a broad distribution of stoichiometries (typically x ≤ y ≤ 3x), both MnMoOy- and Mn 2 MoO y - preferentially formed y = 3 and 4 compositions in significant quantities under our source conditions. The electronic and molecular structures of the Mn x MoO y (x = 1, 2; y = 3, 4) anion and neutral clusters were probed with anion photoelectron spectroscopy and analyzed with supporting density functional theory calculations. Our studies suggest that both metal centers are involved in initial cluster–water complex formation, while Mo is the center that undergoes oxidation; hence, reactivity terminates when Mo is saturated in its highest oxidation state of +6. Across these four clusters, Mn remains relatively reduced and is stable in a high-spin electronic configuration. The preferential reactivity of water molecules toward the Mo center rather than Mn is rationalized by the much lower relative oxophilicity of Mn.

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

Activation of Low-Valent, Multiply M–M Bonded Group VI Dimers toward Catalytic Olefin Metathesis via Surface Organometallic Chemistry

Olefin metathesis is a broadly employed reaction with applications that range from fine chemicals to materials and petrochemicals. The design and investigation of olefin metathesis catalysts have been ongoing for over half a century, with advancements made in terms of activity, stability, and selectivity. Immobilization of organometallic complexes onto solid supports such as silica or alumina is a promising strategy for catalyst heterogenization, often resulting in increased activity and stability. Consequently, a broad range of early transition metal catalysts bearing alkyl, oxide/alkoxide, and amide ligands have been grafted onto silica and their reactivities investigated. In this paper, we report a series of silica-supported tungsten and molybdenum dimers (X 3 M≡MX 3, where M = W and Mo; X = neopentyl, tert-butoxide, and dimethyl amide) and their reactivities toward catalytic olefin metathesis. Dynamic nuclear polarization (DNP)-enhanced solid-state nuclear magnetic resonance (SSNMR), diffuse reflectance infrared Fourier transform (DRIFT), UV resonance Raman, and X-ray absorption (XAS) spectroscopies suggest that upon heterogenization the dimers bind to the surface in a monopodal fashion, with the M≡M triple bond remaining intact. These structural assignments were further corroborated by density functional theory (DFT) calculations. While the homogeneous dimer counterparts are inert, the supported low-valent alkyl W and Mo dimers become active for the disproportionative self-metathesis of propylene to ethylene and butenes and 4-nonene to 4-octene and 5-decene under mild conditions. The lack of activity observed for the free and supported tert-butoxide and dimethyl amide dimers likely suggests that the neopentyl groups are necessary for the formation of a putative alkylidene active species. The difference in reactivity between the free and supported dimers could be explained either by the lowering of the activation barrier of the complex through the electronic effects of the surface or by site isolation of catalytically relevant reactive intermediates.

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

Multi-objective materials bayesian optimization with active learning of design constraints: Design of ductile refractory multi-principal-element alloys

We report Bayesian Optimization (BO) has emerged as a powerful framework to efficiently explore and exploit materials design spaces. To date, most BO approaches to materials design have focused on the materials discovery problem as if it were a single expensive-to-query ‘black box’ in which the target is to optimize a single objective (i.e., material property or performance metric). Also, such approaches tend to be constraint agnostic. Here, we present a novel multi-information BO framework capable of actively learning materials design as a multiple objectives and constraints problem. We demonstrate this framework by optimally exploring a Refractory Multi-Principal-Element Alloy (MPEA) space, here specifically, the system Mo-Nb-Ti-V-W. The MPEAs are explored to optimize two density-functional theory (DFT) derived ductility indicators (Pugh’s Ratio and Cauchy pressure) while learning design constraints relevant to the manufacturing of high-temperature gas-turbine components. Alloys in the BO Pareto-front are analyzed using DFT to gain an insight into fundamental atomic and electronic underpinning for their superior performance, as evaluated within this framework.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unraveling Electronic Trends in O* and OH* Surface Adsorption in the MO 2 Transition-Metal Oxide Series

Understanding the bond strength of O* and OH* intermediates to metal-oxide surfaces is key to predicting the catalytic activity in oxygen-based electrochemistry. Here, we uncover highly non-linear trends in O* and OH* adsorption energies across the 3d, 4d, and 5d series of MO 2 transition-metal (TM) oxide surfaces computed within Hubbard- U corrected density functional theory (DFT + U ). Investigating the electronic structure with crystal orbital Hamiltonian populations (COHP) of the relevant metal–oxygen bonds reveals that the spin-dependent coupling strength between metal-d and oxygen-2p atomic orbitals together with the extent of filling of bonding and anti-bonding orbitals are the primary contributors to the adsorption energy. Importantly, we show that the integrated COHP obtained purely from bulk calculations is a highly accurate descriptor for surface adsorption energetics that captures trends across the group 5–12 TM oxide series within 0.19–0.36 eV. Our results suggest a pathway to prediction of adsorption energies for an arbitrary metal–ligand catalyst system.

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

First-Principles Modeling of the Repassivation of Corrosion Resistant Alloys: Part I. O and Cl Adsorption Energy

To understand the mechanism of repassivation, a key process that controls the performance of corrosion resistant alloys (CRAs), it is useful to obtain some fundamental parameters that describe the atomistic processes. Oxygen and chlorine adsorption energies were systematically calculated with density functional theory for Ni-22Cr alloys with seven different solute elements: Co, Cu, Fe, Mn, Mo, Ru and W. All elements other than Cu enhance both O and Cl adsorption to alloy surfaces. Cr has the strongest effect among these elements. Surface configurations that give stronger O adsorption also tend to give stronger Cl adsorption. While this makes it difficult to differentiate the effect of these solute elements on the oxidation/repassivation tendency of alloys, these fundamental properties provide the foundation for higher scale models that could be used to study the oxidation process during corrosion and potentially inform the design of CRAs. Bader charge analysis reveals a strong relationship between the adatom adsorption energy and the Bader charges of both O and Cl adatom. In conclusion, the work function of alloy surfaces is increased by O while decreased by Cl, suggesting that Cl adsorption can have an intrinsic negative impact on corrosion resistance by decreasing the barrier to electron transfer.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Alloying Effects on the Transport Properties of Refractory High-entropy Alloys

Additive Manufacturing (AM) has opened new frontiers for the design of refractory high-entropy alloys (HEAs) for high-temperature applications. The thermal conductivity of the AM feedstock is among the most important thermo-physical properties that control the melting and solidification process. Despite its significance, there remains a notable gap in both computational and experimental research concerning the thermal conductivity of HEAs. Here, we use density functional theory (DFT) to systematically investigate the alloying effects on the transport properties of Ti-Cr-Mo-W-V-Nb-Ta RHEAs, including electrical and thermal conductivities and Seebeck coefficient. The relaxation time of charge carriers is a key underlying parameter determining thermal conductivity that is exceedingly challenging to predict from first principles alone, and we thus follow the approach by Mukherjee, Satsangi, and Singh [Chem Mater 32, 6507 (2022)] to optimize the relaxation time for RHEAs. Here we validated thermal conductivity predictions on elemental solids, binary and ternary alloys, and RHEAs and compared them against thermodynamic (CALPHAD) predictions and our experiments with good correlations. To understand observed trends in thermal conductivity, we assessed the phase stability, electronic structure, phonon, and intrinsic- and tensile strength of down-selected RHEAs. Our electronic structure and phonon results connect well with the observed compositional trends for thermal transport in RHEAs. Our DFT assessment and CALPHAD predictions provide a unique design guide for RHEAs with tailored thermal conductivity, a critical consideration for AM and thermal-management applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bioinspired Activation of N 2 on Asymmetrical Coordinated Fe Grafted 1T MoS 2 at Room Temperature

Inspired by the nitrogen fixation process on MoFe nitrogenase, asymmetrical coordinated Fe grafted onto 1T MoS 2 were successfully synthesized. The unique electron -rich structure with asymmetrical coordination made the 1T Fe 0.1 Mo 0.9 S 2 layered material actively react with water and dinitrogen at room temperature and atmosphere pressure. Subsequently, ammonia can be produced with a yield of 800 pmol (NH 4 + ) g -1 (12.5% yield in mole). The activation, fixation and reduction of dinitrogen were confirmed by isotopically labeled experiments. The location and the specific coordination environment of grafted Fe in Fe-Mo-S were further determined by X-ray absorption spectroscopy analysis. Our work demonstrates that the nitrogen fixation and reduction for ammonia at room temperature without any chemical and electrochemical assistance is distinctly different from traditional bionic -inspired nitrogen fixation process. The mechanism of the activation and reduction of N 2 was further investigated by density functional theory calculation and Raman spectra. Compared with 1T MoS 2 , the enriched electron nature and asymmetrical coordination of Fe in Fe-Mo-S materials play a critical role in the bioinspired activation of N 2 at ambient condition.

1T MoS2↗
DOE OSTI2021

Not all platinum surfaces are the same: Effect of the support on fundamental properties of platinum adlayer and its implications for the activity toward hydrogen evolution reaction

The adsorption of atomic hydrogen on a platinum monolayer supported on orthorhombic Mo 2 C(100) surface has been investigated, considering different hydrogen surface coverages. Calculations have been performed using density functional theory with the Perdew-Burke-Ernzerhof exchange correlation functional and a D3 Van der Waals corrections. The theoretical insight has been gained into atomic hydrogen interaction with Pt monolayer, supported on both molybdenum and well-studied tungsten carbide, and considering hydrogen surface coverage. Fundamental properties of Pt adlayer depend on the support, affecting hydrogen evolution activity of the resulting systems. At low hydrogen coverage all systems, with the exception of Pt supported on the molybdenum-terminated Mo 2 C, adsorb H comparably to a pristine Pt(111) surface and their high activity for the hydrogen evolution reaction is predicted. Additionally, at higher coverages supported Pt monolayers interact with atomic hydrogen unlike the Pt(111), suggesting that the activity of the supported and unsupported platinum toward hydrogen evolution reaction have different origins. Furthermore, the position of the supported platinum monolayers on the volcano curve is a function of the surface coverage, more so than for extended metal surfaces. Therefore, hydrogen surface coverage is a key variable to understand the catalytic potential, approaching towards an improved model for screening of electrocatalytic systems.

36 MATERIALS SCIENCE↗