Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “lead optimization”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Lewis Acid Site Assisted Bifunctional Activity of Tin Doped Gallium Oxide and Its Application in Rechargeable Zn‐Air Batteries

The enhanced safety, superior energy, and power density of rechargeable metal‐air batteries make them ideal energy storage systems for application in energy grids and electric vehicles. However, the absence of a cost‐effective and stable bifunctional catalyst that can replace expensive platinum (Pt)‐based catalyst to promote oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) at the air cathode hinders their broader adaptation. Here, in this work, it is demonstrated that Tin (Sn) doped β‐gallium oxide (β‐Ga 2 O 3 ) in the bulk form can efficiently catalyze ORR and OER and, hence, be applied as the cathode in Zn‐air batteries. The Sn‐doped β‐Ga 2 O 3 sample with 15% Sn (Sn x =0.15 ‐Ga 2 O 3 ) displayed exceptional catalytic activity for a bulk, non‐noble metal‐based catalyst. When used as a cathode, the excellent electrocatalytic bifunctional activity of Sn x =0.15 ‐Ga 2 O 3 leads to a prototype Zn‐air battery with a high‐power density of 138 mW cm −2 and improved cycling stability compared to devices with benchmark Pt‐based cathode. The combined experimental and theoretical exploration revealed that the Lewis acid sites in β‐Ga 2 O 3 aid in regulating the electron density distribution on the Sn‐doped sites, optimize the adsorption energies of reaction intermediates, and facilitate the formation of critical reaction intermediate (O*), leading to enhanced electrocatalytic activity.

Lewis acid sites↗
DOE OSTI2026

Origin of Electrochemical Activation Leading to Enhanced Cycling Stability of Li‐ and Mn‐Rich Cathodes

Electrochemical activation is a critical step for optimal functioning of Li- and Mn-rich (LMR) cathodes, yet the underlying mechanism for such activation remains elusive. Here, by using scanning/transmission electron microscopy (S/TEM) combined with the associated energy-dispersive x-ray spectroscopy (EDS) and electron energy-loss spectroscopy (EELS), we decipher the origin of the activation enhanced electrochemical properties. We reveal that activation induces the formation of a spinel-like phase within the C2/m domains of the LMR cathode, where the transition-metal ions partially occupy both the tetrahedral (8a) and octahedral (16c) sites of the $Fd\bar{3}m$ spinel lattice, distinguishing the spinel-like phase from the conventional high-voltage spinel. Systematic varying the cycling voltage reveals a critical activation voltage above which this spinel-like phase forms, while lower voltages preserve the layered bulk structure. As the spinel-like phase is a stable structure for electrochemical cycling, the present findings provide direct mechanistic insight into the voltage-dependent activation process and explain how the C2/m to spinel-like transformation upon activation contributes to the electrochemical performance of LMR cathodes, providing guidance for the rational design of Li-rich cathodes with enhanced cycling durability.

Li-rich and Mn-rich cathode↗
DOE OSTI2022

Optimal control to handle variations in moisture content and reactor in-feed rate

The variations in feedstock characteristics, such as moisture and particle size distribution, lead to an inconsistent flow of feedstock from the biomass pre-processing system to the reactor in-feed system. These inconsistencies result in low on-stream times at the reactor in-feed equipment. This research develops an optimal process control method for a biomass pre-processing system comprised of milling and densification operations to provide the consistent flow of feedstock to a reactor's throat. This method uses a mixed-integer optimization model to identify optimal bale sequencing, equipment in-feed rate, and buffer location and size in the biomass pre-processing system. This method, referred to as the hybrid process control (HPC), aims to maximize throughput over time. We compare HPC with a baseline feed forward process control. Our case study based on switchgrass finds that HPC reduces the variation of a reactor's feeding rate by up to 100% without increasing the operating cost of the biomass pre-processing system for biomass with moisture ranging from 10 to 25%. Additionally, HPC reduces the cost of processing biomass by 0.36%–2.22%, and reduces processing time by 0.35%-2.24%. Furthermore, a biorefinery can adapt HPC to achieve its design capacity.

09 BIOMASS FUELS↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2025

Consensus-Based Control and Optimization of Power System Inertia

The integration of distributed energy resources (DERs), such as solar, wind, and energy storage systems, into power grids through inverter-based resources reduces power system inertia, leading to faster frequency dynamics and potential grid instability. To address this challenge, this paper proposes a distributed, consensus-based approach for the real-time control and optimization of inertia sources (synchronous generators and/or DERs) during system disturbances, enhancing both system stability and economic performance. The distributed control and optimization approach assumes each inertia source exchanges information solely with its neighbor ones, making it easily scalable to large power grid networks. The impacts of the communication connectivity among the inertia sources as well as their generation capacity limits on the distributed approach are investigated. We also demonstrate the approach’s robustness in scenarios involving communication time delays and packet losses, validating its effectiveness through numerical simulations on a 4-bus test system and a two-area 8-bus test system.

Power system inertia↗
DOE OSTI2022

Physics makes the difference: Bayesian optimization and active learning via augmented Gaussian process

Abstract Both experimental and computational methods for the exploration of structure, functionality, and properties of materials often necessitate the search across broad parameter spaces to discover optimal experimental conditions and regions of interest in the image space or parameter space of computational models. The direct grid search of the parameter space tends to be extremely time-consuming, leading to the development of strategies balancing exploration of unknown parameter spaces and exploitation towards required performance metrics. However, classical Bayesian optimization (BO) strategies based on the Gaussian process (GP) do not readily allow for the incorporation of the known physical behaviors or past knowledge. Here we explore a hybrid optimization/exploration algorithm created by augmenting the standard GP with a structured probabilistic model of the expected system’s behavior. This approach balances the flexibility of the non-parametric GP approach with a rigid structure of physical knowledge encoded into the parametric model. The fully Bayesian treatment of the latter allows additional control over the optimization via the selection of priors for the model parameters. The method is demonstrated for a noisy version of a standard univariate test function used to evaluate optimization algorithms and further extended to physical lattice models. This methodology is expected to be universally suitable for injecting prior knowledge in the form of physical models and past data in the BO framework.

42 ENGINEERING↗
DOE OSTI2024

High-throughput combinatorial approach expedites the synthesis of a lead-free relaxor ferroelectric system

Developing novel lead-free ferroelectric materials is crucial for next-generation microelectronic technologies that are energy efficient and environment friendly. However, materials discovery and property optimization are typically time-consuming due to the limited throughput of traditional synthesis methods. In this work, we use a high-throughput combinatorial synthesis approach to fabricate lead-free ferroelectric superlattices and solid solutions of (Ba 0.7 Ca 0.3 )TiO 3 (BCT) and Ba(Zr 0.2 Ti 0.8 )O 3 (BZT) phases with continuous variation of composition and layer thickness. High-resolution x-ray diffraction (XRD) and analytical scanning transmission electron microscopy (STEM) demonstrate high film quality and well-controlled compositional gradients. Ferroelectric and dielectric property measurements identify the “optimal property point” achieved at the composition of 48BZT–52BCT. Displacement vector maps reveal that ferroelectric domain sizes are tunable by varying {BCT–BZT} N superlattice geometry. This high-throughput synthesis approach can be applied to many other material systems to expedite new materials discovery and properties optimization, allowing for the exploration of a large area of phase space within a single growth.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Standard Operating Procedure for Optimal Deployment of Meteorological Instrumentation Within the Solar Radiation Research Laboratory (BMS and RMS SOP)

The objective of the National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) is to collect and use high-quality solar radiation data sets for research leading to the widespread adoption of solar technologies. To appropriately populate and track the diverse array of instruments at the NREL-SRRL, NREL has established a Standard Operating Procedure (SOP) for optimal instrument deployment within the SRRL for both the Baseline Measurement System (BMS) and the Research Measurement System (RMS). Using best practices methodologies, the NREL-SRRL maintains a varied and extensive array of solar monitoring equipment to test, evaluate, and characterize the solar sensors used by federal and international agencies as well as the solar industry to determine the solar resource. The SOP provides the industry with guidance for solar resource assessment and is used for procedures in the long-term continuous monitoring of legacy instruments alongside state-of-the-art instruments. Based on the SOP, instruments are annually evaluated for continued deployment. Instruments that do not meet the SOP criteria are decommissioned, and new instruments that meet the criteria are deployed. Streamlining and optimizing the use of this facility ensures that the lab continues to be a world-leading solar calibration and measurement facility.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Optimal operation of solid-oxide electrolysis cells considering long-term chemical degradation

Optimizing the performance of solid oxide electrolysis cells (SOECs) for long-term hydrogen (H 2 ) production at high temperatures is crucial, as prolonged operation leads to efficiency losses and shorter cell lifespans due to chemical degradation. Here, in this work, we adopt a quasi-steady state approach for dynamic optimization over extended operational periods to address the disparity in timescales between cell operation and degradation. Integrating a 2-D non-isothermal SOEC model with balance-of-plant (BOP) equipment, we explore three optimization objectives: minimizing terminal degradation, maximizing integral efficiency, and minimizing the levelized cost of H 2 (LCOH). Our dynamic optimization algorithm reduces LCOH by 9.5% and 16% compared to strategies focusing solely on terminal degradation and integral efficiency, respectively. For electricity prices of 0.03 $\$$/mWh and 0.3 $\$$ mWh optimal replacement schedules range from 5 to 2 years, depending on the operational mode. Furthermore, a flexible operational mode yields additional improvements in LCOH over traditional galvanostatic and potentiostatic modes.

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

Assessment of vertical stability for negative triangularity pilot plants

Abstract Negative triangularity (NT) tokamak configurations may be more susceptible to magneto-hydrodynamic instability, posing challenges for recent reactor designs centered around their favorable properties, such as improved confinement and operation free of edge-localized modes. In this work, we assess the vertical stability of plasmas with NT shaping and develop potential reactor solutions. When coupled with a conformal wall, NT equilibria are confirmed to be less vertically stable than equivalent positive triangularity (PT) configurations. Unlike PT, their vertical stability is degraded at higher poloidal beta. Furthermore, improvements in vertical stability at low aspect ratio do not translate to the NT geometry. NT equilibria are stabilized in PT vacuum vessels due to the increased proximity of the plasma and the wall on the outboard side, but this scenario is found to be undesirable due to reduced vertical gaps which give less spatial margin for control recovery. Instead, we demonstrate that informed positioning of passively conducting plates can lead to improved vertical stability in NT configurations on par with stability metrics expected in PT scenarios. An optimal setup for passive plates in highly elongated NT devices is presented, where plates on the outboard side of the device reduce vertical instability growth rates to 16% of their baseline value. For lower target elongations, integration of passive stabilizers with divertor concepts can lead to significant improvements in vertical stability. Plates on the inboard side of the device are also uniquely enabled in NT geometries, providing opportunity for spatial separation of vertical stability coils and passive stabilizers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Influence of counterion substitution on the properties of imidazolium-based ionic liquid clusters

Due to their unique physiochemical properties that may be tailored for specific purposes, ionic liquids (ILs) have been investigated for various applications, including chemical separations, catalysis, energy storage, and space propulsion. The different cations and anions comprising ILs may be selected to optimize a range of desired properties, such as thermal stability, ionic conductivity, and volatility, leading to the designation of certain ILs as designer “green” solvents. The effect of counterions on the properties of ILs is of both fundamental scientific interest and technological importance. Herein, we report a systematic experimental and theoretical investigation of the size, charge, stability toward dissociation, and geometric/electronic structure of 1-ethyl-3-methyl imidazolium (EMIM)-based IL clusters containing two different atomic counterions (i.e., bromide [Br − ] and iodide [I − ]). This work extends our studies of EMIM + cations with atomic chloride (Cl − ) and molecular tetrafluoroborate (BF 4 − ) anions reported previously by Baxter et al. [Chem. Mater. 34, 2612 (2022)] and Zhang et al . [J. Phys. Chem. Lett. 11, 6844 (2020)], respectively. Distributions of anionic IL clusters were generated in the gas phase using electrospray ionization and characterized by high mass resolution mass spectrometry, energy-resolved collision-induced dissociation, and negative ion photoelectron spectroscopy experiments. The experimental results reveal anion-dependent trends in the size distribution, relative abundance, ionic charge state, stability toward dissociation, and electron binding energies of the IL clusters. Complementary global optimization theory provides molecular-level insights into the bonding and electronic structure of a selected subset of clusters, including their low energy structures and electrostatic potential maps, and how these fundamental characteristics are influenced by anion substitution. Collectively, our findings demonstrate how the fundamental properties of ILs, which determine their suitability for many applications, may be tuned by substituting counterions. These observations are critical in the sub-nanometer cluster size regime where phenomena do not scale predictably to the bulk phase, and each atom counts toward determining behavior.

cluster↗
DOE OSTI2023

Ferroelectric Domain Wall Engineering Enables Thermal Modulation in PMN–PT Single Crystals

Acting like thermal resistances, ferroelectric domain walls can be manipulated to realize dynamic modulation of thermal conductivity ( k ), which is essential for developing novel phononic circuits. Despite the interest, little attention has been paid to achieving room-temperature thermal modulation in bulk materials due to challenges in obtaining a high thermal conductivity switching ratio ( k high / k low ), particularly in commercially viable materials. Here, room-temperature thermal modulation in 2.5 mm-thick Pb(Mg 1/3 Nb 2/3 )O 3 – x PbTiO 3 (PMN– x PT) single crystals is demonstrated. With the use of advanced poling conditions, assisted by the systematic study on composition and orientation dependence of PMN– x PT, a range of thermal conductivity switching ratios with a maximum of ≈1.27 is observed. Further, simultaneous measurements of piezoelectric coefficient (d 33 ) to characterize the poling state, domain wall density using polarized light microscopy (PLM), and birefringence change using quantitative PLM reveal that compared to the unpoled state, the domain wall density at intermediate poling states (0< d 33 33,max ) is lower due to the enlargement in domain size. At optimized poling conditions (d 33,max ), the domain sizes show increased inhomogeneity that leads to enhancement in the domain wall density. This work highlights the potential of commercially available PMN– x PT single crystals among other relaxor-ferroelectrics for achieving temperature control in solid-state devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Influence of microstructure and temperature on impact toughness of H13 steel produced by binder jet additive manufacturing

Binder Jet Additive Manufacturing (BJAM) is a promising manufacturing pathway to produce H13 steel dies and tooling with complex geometries for applications in high pressure aluminum die casting, hot stamping, and injection molding. While fully dense H13 coupons produced using BJAM have been subjected to detailed microstructure characterization, properties which are critical and relevant to the aforementioned applications, such as impact toughness, have not been reported. Here, this work evaluated the influence of microstructural characteristics and test temperature on the impact toughness of H13 produced by BJAM. Coupons were produced from three different powder size distributions (PSDs), with nominal powder size ranges of −22 μm, 10-32 μm, and 15-53 μm. Coupons were printed, sintered, hot isostatically pressed (HIPed), and heat treated by quenching and tempering. After HIPing, the measured porosities of the different PSDs were all less than 0.015 vol%. The −22 μm PSD BJAM material exhibited the best impact toughness of all three PSDs across the entire test temperature range from 25 to 400 °C, and also exhibited an ambient temperature impact toughness of 10.9 J at a hardness of 46 HRC, equivalent to minimum threshold requirements for premium grade H13. However, all three BJAM PSDs exhibited moderately lower impact toughness than premium grade wrought H13 from RT to 400 °C. We attribute this in part due to the very large prior-austenite grain (PAG) sizes near 400 μm and segregation from former liquid channels leading to greater amounts of carbide precipitation on PAG boundaries. Technical pathways for optimizing processing and chemistry of BJAM H13 to further improve impact toughness are discussed.

Binder jet additive manufacturing↗
DOE OSTI2021

Modeling and experimental validation of a W f /W-fabrication by chemical vapor deposition and infiltration

Tungsten (W) has a unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall material in future fusion devices. However, its intrinsic brittleness and its further embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) with extrinsic toughening is being developed. A possible synthesis route is chemical vapor deposition (CVD) using heated W fabrics as substrate. The challenge is that the growing CVD-W can isolate domains from precursor access leading to strength-reducing pores. To deepen the process understanding and to optimize the CVD parameters, models were developed with COMSOL Multiphysics and validated experimentally. W deposition rate equations as function of the temperature and the partial pressures of the precursors H 2 and WF 6 were experimentally validated in previous work. In the present article, these equations are applied to obtain partial pressures within the CVD reactor. The results are taken as input for transient simulations in the microscale, in which W coatings, growing onto multiple adjacent W fibers, were simulated via mesh deformation and remeshing. The surface-to-surface contact of the W coatings and the corresponding potential pore formation were simulated by implementing sophisticated deposition rate stop conditions. Within the measuring uncertainties of ≃ ±1%, the models are validated successfully by experimental comparison regarding the deposition rate, pore structure, and relative densities ranging from 0.6 to 0.9.

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

Conversion efficiency improvement of ELO GaAs solar cell, deposited on water soluble sacrificial buffer

Here, we demonstrate the improvement of power conversion efficiency of an epitaxially-lifted-off single junction GaAs solar cell deposited on a water-soluble sacrificial buffer architecture by introducing an additional germanium (Ge) interlayer between GaAs and the fluoride buffer. The epitaxial lift-off (ELO) technique has been extensively used to separate III-V device layers from their single crystal GaAs substrates. However, conventional ELO requires the use of concentrated hydrofluoric acid (HF) for extended times to etch out the sacrificial layer, subsequently degrading the surface roughness of the parent wafer. As a result, the wafer has to undergo expensive and intensive chemical mechanical polishing (CMP) processes, costing about 25% of a pristine 6-inch GaAs substrate. In our previous work, we demonstrated a method to eliminating the need for CMP post-processing by using water-assisted ELO (H 2 O-ELO). A water-soluble, 3-layer buffer architecture was developed using alkaline earth compounds. However, devices suffered low performance due to a high defect density in the GaAs active layers. A Ge interlayer was introduced to provide a more favorable surface energy for GaAs growth, which leads to an improvement of GaAs crystal quality. The Ge deposition conditions were optimized to achieve a high-quality Ge layer on triple-layer fluoride buffer. Single junction GaAs devices fabricated on Ge/(Ca,Sr)F 2 /BaF 2 /(Ca,Sr)F 2 showed improvement of solar cell performance parameters featuring increases in V oc by 23.7%, and fill factor (F.F.) by 4.9%, which results in an overall improvement of power conversion efficiency from 10.3% to 12.69%.

14 SOLAR ENERGY↗
DOE OSTI2024

Predicting the X-ray Absorption Spectrum of Ozone with Single Configuration State Functions

X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn–Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids can be approximated with individual CSFs, despite the presence of multireference character in the ground state. Here, we also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O 3 . This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.

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

Arylammonium-Assisted Reduction of the Open-Circuit Voltage Deficit in Wide-Bandgap Perovskite Solar Cells: The Role of Suppressed Ion Migration

Surface treatment using large alkyl/aryl ammonium cations has demonstrated reduced open-circuit voltage (V OC ) deficits in perovskite solar cells (PSCs), but the origin of the improvements has been vaguely attributed to defect passivation. Here, we combine the microscopic probing of the local electrical properties, thermal admittance spectroscopic analysis, and first principles calculations to elucidate the critical role of arylammonium interface layers on suppressing the ion migration in wide-bandgap (WBG) PSCs. Our results reveal that arylammonium surface treatment using phenethylammonium iodide increases the activation energy barrier for ion migration on the surface, which suppresses the accumulation of charge defects at surface and grain boundaries (GBs), leading to reduced dark saturation current density in WBG PSCs. Furthermore, with device optimization, our champion 1.73-eV PSC delivers a power conversion efficiency of 19.07% with a V OC of 1.25 V, achieving a V OC -deficit of 0.48 V.

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

Hierarchal structures tuned electrocaloric and electromechanical performance in PVDF-based tetrapolymers

Ferroelectrics with multifunctionalities are gaining increased interest in self-actuated electrocaloric effect (ECE) refrigerators. However, achieving high ECE and electromechanical (EM) coupling concomitantly for maximum heat transfer remains challenging. Here we present the structure-property relationship for poly(vinylidene fluoride-co-trifluoroethylene-co-chlorofluoroethylene-co-double bond), P(VDF-TrFE-CFE-DB), tetrapolymer, which exhibited a high ECE entropy change of 66.5 J Kg −1 K −1 and EM strain of −6.1%. We show that thermal treatment can be a key factor influencing multifunctional properties. High-temperature annealing incorporates DB and CFE units into crystalline grains to form extended-chain crystals, enabling CFE units to induce relaxor behavior and DB units to induce large structural changes at low electric fields. This synergy leads to an enhancement in both ECE and EM performances. Furthermore, at an optimized temperature of 50 °C, the annealed films exhibit giant cross-energy coupling, achieving ECE and EM performances of 100.8 J Kg −1 K −1 and −7.6%. This study provides insights into developing new ferroelectric polymers with electroactive multifunctionalities.

36 MATERIALS SCIENCE↗