Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic Interface”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Unveiling interfaces between In-rich and Ga-rich GaInP vertical slabs of laterally composition modulated structures

Here, we report changes at the interface between Ga-rich/In-rich GaInP vertical slabs in laterally composition modulated (LCM) GaInP as a function of the V/III ratio. The photoluminescence exhibits satellite peaks, indicating that the parasitic potential between the GaInP vertical slabs disappears as the V/III ratio decreases. However, a high V/III ratio leads to an abrupt interface, increasing the parasitic potential because of the phosphorus-amount-dependent diffusion of group-III atoms during growth. These results suggest that the V/III ratio is an important parameter that must be wisely chosen in designing optoelectronic devices incorporating LCM structure.

36 MATERIALS SCIENCE↗
DOE OSTI2014

The importance of anharmonicity in thermal transport across solid-solid interfaces

Understanding interfacial thermal transport is of great importance for applications like energy devices and thermal management of electronics. Despite the significant efforts in the past few decades, thermal transport across solid-solid interfaces is still not fully understood and cannot be accurately predicted. Anharmonicity is often ignored in many prediction models, such as the mismatch models, the wave-packet method, and the Atomic Green's function. In this paper, we use molecular dynamics to systematically study the role of anharmonicity in thermal transport across solid-solid interfaces. The interatomic interactions are modeled using force constants up to the third order. This model allows controlling the anharmonicity independently by tuning the cubic force constants. The interfacial thermal conductance as a function of anharmonicity inside the materials and that at the interface is studied. We found that the anharmonicity inside the materials plays an important role in the interfacial thermal transport by facilitating the energy communication between different phonon modes. The anharmonicity at the interface has much less impact on the interfacial thermal transport. These results are important to the modification of traditional models to improve their prediction power.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Surface-Directed Fabrication of Integrated Membrane-Electrode Interfaces

In the past five years, we have established a research program to investigate the integration of electrode/catalyst support, gas diffusion pathways, and ionomer into a single network based on a bottom-up, surface-directed approach. Here, the porous electrode surface provides the common foundation for both catalyst deposition and ionomer growth to develop a highly controlled process. Our efforts have focused on (1) the preparation of surface-initiated ionomer that can be grown from electrodes of essentially any geometry and is chemisorbed to the active electrode surface, and (2) the integration of this ionomer with both planar and 3-D electrodes that are atomically modified by Pt catalysts. Figure 1 shows a schematic of a single pore of such an integrated cathode where a catalytic Pt atomic layer is supported along the pores of a nanoporous gold electrode. The catalyst layer is electronically accessible from the metal cathode and in contact with proton-conducting ionomer to yield extremely high utilization efficiencies of the catalyst (approaching 100%) with the ultralow loadings expected from atomic overlayers. We have introduced a new class of ionomer film that is grown from the catalyst/electrode interface to develop proton conducting pathways throughout the 3-D electrode and to the adsorbed catalyst particles. This ionomer is a comb polymer designed to contain sulfonate-lined channels for rapid proton conduction as well as hydrophobic domains (fluorocarbon or hydrocarbon) to promote O 2 solubility and transfer. This interfacial ionomer is designed to interface with a bulk membrane (such as Nafion) and an anode that could be engineered in a similar way. A key conclusion thus far is that ionomer and Pt catalyst can be integrated at a single interface (both 2-D and 3-D) to provide a high level of control over the construction of electrodes. To expand our efforts to date, we seek to utilize the materials chemistry we have developed to investigate how the depth-dependent structure of an ionomer film affects its performance in transferring protons and oxygen to accelerate the rate of the oxygen reduction reaction. By correlating structure and composition to performance and identifying barriers that limit efficiency, we seek to design ultra high performance ionomer films that are customized for the oxygen reduction reaction.

08 HYDROGEN↗
DOE OSTI2020

Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data

Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Mixed Conducting Electrodes for Better AMTEC Cells

Electrode materials that exhibit mixed conductivity (that is, both electronic and ionic conductivity) have been investigated in a continuing effort to improve the performance of the alkali metal thermal-to-electric converter (AMTEC). These electrode materials are intended primarily for use on the cathode side of the sodium-ion-conducting solid electrolyte of a sodium-based AMTEC cell. They may also prove useful in sodium-sulfur batteries, which are under study for use in electric vehicles. An understanding of the roles played by the two types of conduction in the cathode of a sodium-based AMTEC cell is prerequisite to understanding the advantages afforded by these materials. In a sodium-based AMTEC cell, the anode face of an anode/solid-electrolyte/cathode sandwich is exposed to Na vapor at a suitable pressure. Upon making contact with the solid electrolyte on the anode side, Na atoms oxidize to form Na+ ions and electrons. Na+ ions then travel through the electrolyte to the cathode. Na+ ions leave the electrolyte at the cathode/electrolyte interface and are reduced by electrons that have been conducted through an external electrical load from the anode to the cathode. Once the Na+ ions have been reduced to Na atoms, they travel through the cathode to vaporize into a volume where the Na vapor pressure is much lower than it is on the anode side. Thus, the cathode design is subject to competing requirements to be thin enough to allow transport of sodium to the low-pressure side, yet thick enough to afford adequate electronic conductivity. The concept underlying the development of the present mixed conducting electrode materials is the following: The constraint on the thickness of the cathode can be eased by incorporating Na+ -ionconducting material to facilitate transport of sodium through the cathode in ionic form. At the same time, by virtue of the electronically conducting material mixed with the ionically conducting material, reduction of Na+ ions to Na atoms can take place throughout the thickness of the cathode. The net effect is to reduce the diffusion and flow resistance to sodium through the electrode while reducing the electronic resistance by providing shorter conduction paths for electrons. Reduced resistance to both sodium transport and electronic conductivity results in an increase in electric power output.

Ryan, Margaret↗
DOE OSTI2014

Influence of oxygen pressure and aging on LaAlO{sub 3} films grown by pulsed laser deposition on SrTiO{sub 3} substrates

The crystal structures of LaAlO{sub 3} films grown by pulsed laser deposition on SrTiO{sub 3} substrates at oxygen pressure of 10{sup −3} millibars or 10{sup −5} millibars, where kinetics of ablated species hardly depend on oxygen background pressure, are compared. Our results show that the interface between LaAlO{sub 3} and SrTiO{sub 3} is sharper when the oxygen pressure is lower. Over time, the formation of various crystalline phases is observed while the crystalline thickness of the LaAlO{sub 3} layer remains unchanged. X-ray scattering as well as atomic force microscopy measurements indicate three-dimensional growth of such phases, which appear to be fed from an amorphous capping layer present in as-grown samples.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗
DOE OSTI2019

Defects‐Induced In‐Plane Heterophase in Cobalt Oxide Nanosheets for Oxygen Evolution Reaction

Abstract Cobalt oxides as efficient oxygen evolution reaction (OER) electrocatalysts have received much attention because of their rich reserves and cheap cost. There are two common cobalt oxides, Co 3 O 4 (spinel phase, stable but poor intrinsic activity) and CoO (rocksalt phase, active but easily be oxidatized). Constructing Co 3 O 4 /CoO heterophase can inherit both characteristic features of each component and form a heterophase interface facilitating charge transfer, which is believed to be an effective strategy in designing excellent electrocatalysts. Herein, an atomic arrangement engineering strategy is applied to improve electrocatalytic activity of Co 3 O 4 for the OER. With the presence of oxygen vacancies, cobalt atoms at tetrahedral sites in Co 3 O 4 can more easily diffuse into interstitial octahedral sites to form CoO phase structure as revealed by periodic density functional theory computations. The Co 3 O 4 /CoO spinel/rocksalt heterophase can be in situ fabricated at the atomic scale in plane. The overpotential to reach 10 mA cm −2 of Co 3 O 4 /CoO is 1.532 V, which is 92 mV smaller than that of Co 3 O 4 . Theoretical calculations confirm that the excellent electrochemical activity is corresponding to a decline in average p ‐state energy of adsorbed‐O on the Co 3 O 4 /CoO heterophase interface. The reaction Gibbs energy barrier has been significantly decreased with the construction of the heterophase interface.

Liu, Zhijuan↗
DOE OSTI2018

Diffusion-driven Fluid Dynamics in Ideal Gases and Plasmas

The classical transport theory based on Chapman-Enskog methods provides self-consistent approximations for the kinetic flux of mass, heat, and momentum in a fluid limit characterized with a small Knudsen number. The species mass fluxes relative to the center of mass, or “diffusive fluxes,” are expressed as functions of known gradient quantities with kinetic coefficients evaluated using similar analyses for mixtures of gases or plasma components. The sum over species of the diffusive mass fluxes is constrained to be zero in the Lagrange frame, and thus results in a non-zero molar flux leading to a pressure perturbation. At an interface between two species initially in pressure equilibrium, the pressure perturbation driven by the diffusive molar flux induces a center of mass velocity directed from the species of greater atomic mass towards the lighter atomic mass species. As the ratio of the species particle masses increases, this center of mass velocity carries an increasingly greater portion of the mass across the interface and for a particle mass ratio greater than about two, the center of mass velocity carries more mass than the gradient driven diffusion flux. Early time transients across an interface between two species in a 1D plasma regime and initially in equilibrium are compared using three methods; a fluid code with closure in a classical transport approximation, a particle in cell simulation, and an implicit Fokker-Planck solver for the particle distribution functions. The early time transient phenomenology is shown to be similar in each of the computational simulation methods, including a pressure perturbation associated with the stationary “induced” component of the center of mass velocity which decays to pressure equilibrium during diffusion. At early times, the diffusive process generates pressure and velocity waves which propagate outward from the interface and are required to maintain momentum conservation. The energy in the outgoing waves dissipates as heat in viscous regions, and it is hypothesized that these diffusion driven waves may sustain fluctuations in less viscid finite domains after reflections from the boundaries. Finally, these fluid dynamic phenomena are similar in gases or plasmas and occur in flow transients with a moderate Knudsen number. The analysis and simulation results show how the kinetic flux, represented in the fluid transport closure, directly modifies the mass averaged flow described with the Euler equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dynamic Behavior of Spatially Confined Sn Clusters and Its Application in Highly Efficient Sodium Storage with High Initial Coulombic Efficiency

Advanced battery electrodes require a cautious design of microscale particles with built-in nanoscale features to exploit the advantages of both micro- and nano-particles relative to their performance attributes. Herein, the dynamic behavior of nanosized Sn clusters and their host pores in carbon nanofiber) during sodiation and desodiation is revealed using a state-of-the-art 3D electron microscopic reconstruction technique. For the first time, the anomalous expansion of Sn clusters after desodiation is observed owing to the aggregation of clusters/single atoms. Pore connectivity is retained despite the anomalous expansion, suggesting inhibition of solid electrolyte interface formation in the sub-2-nm pores. Taking advantage of the built-in nanoconfinement feature, the CNF film with nanometer-sized interconnected pores hosting Sn clusters (≈2 nm) enables high utilization (95% at a high rate of 1 A g –1 ) of Sn active sites while maintaining an improved initial Coulombic efficiency of 87%. Finally, the findings provide insights into electrochemical reactions in a confined space and a guiding principle in electrode design for battery applications.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Thickness-based adaptive mesh refinement methods for multi-phase flow simulations with thin regions

In numerical simulations of multi-scale, multi-phase flows, grid refinement is required to resolve regions with small scales. A notable example is liquid-jet atomization and subsequent droplet dynamics. It is essential to characterize the detailed flow physics with variable length scales with high fidelity, in order to elucidate the underlying mechanisms. In this paper, two thickness-based mesh refinement schemes are developed based on distance- and topology-oriented criteria for thin regions with confining wall/plane of symmetry and in any situation, respectively. Both techniques are implemented in a general framework with a volume-of-fluid formulation and an adaptive-mesh-refinement capability. The distance-oriented technique compares against a critical value, the ratio of an interfacial cell size to the distance between the mass center of the cell and a reference plane. The topology-oriented technique is developed from digital topology theories to handle more general conditions. The requirement for interfacial mesh refinement can be detected swiftly, without the need of thickness information, equation solving, variable averaging or mesh repairing. The mesh refinement level increases smoothly on demand in thin regions. The schemes have been verified and validated against several benchmark cases to demonstrate their effectiveness and robustness. These include the dynamics of colliding droplets, droplet motions in a microchannel, and atomization of liquid impinging jets. Overall, the thickness-based refinement technique provides highly adaptive meshes for problems with thin regions in an efficient and fully automatic manner.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Emerging Electrochemical Techniques for Probing Site Behavior in Single-Atom Electrocatalysts

Single-atom catalysts (SACs) have aroused tremendous interest over the past decade, particularly in the community of energy and environment-related electrocatalysis. A rapidly growing number of recent publications have recognized it as a promising candidate with maximum atomic utilization, distinct activity, and selectivity in comparison to bulk catalysts and nanocatalysts. However, the complexity of localized coordination environments and the dispersion of isolated sites lead to significant difficulties when it comes to gaining insight into the intrinsic behavior of electrocatalytic reactions. Furthermore, the low metal loadings of most SACs make conventional ensemble measurements less likely to be accurate on the subnanoscale. Thus, it remains challenging to probe the activity and properties of individual atomic sites by available commercial instruments and analytical methods. In spite of this, continuing efforts have lately focused on the development of advanced measurement methodologies, which are very useful to the fundamental understanding of SACs. There have recently been a number of in situ/operando techniques applied to SACs, such as electron microscopy, spectroscopy, and other analysis methods, which support relevant functions to identify the active sites and reaction intermediates and to investigate the dynamic behavior of localized structures of the catalytic sites. This Account aims to present recent electrochemical probing techniques which can be used to identify single-atomic catalytic sites within solid supports. First, we describe the basic principles of molecular probe methods for the study and analysis of electrocatalytic site behavior. In particular, the in situ probing technique enabled by surface interrogation scanning electrochemical microscopy (SISECM) can measure the active site density and kinetic rate with high resolution. An alternative electrochemical probing technique is further demonstrated on the basis of single-entity electrochemistry, which allows the unique electrochemical imaging of the size and catalytic rate of single atoms, molecules, and clusters. Next, the merits and limitations of different electrochemical techniques are then discussed, along with perspectives for future prospects. Apart from this, we further showcase the powerful capability of emerging electrochemical probing techniques for determining significant effects and properties of SACs for various electrocatalytic reactions, including oxygen reduction and evolution, hydrogen evolution, and nitrate reduction. Overall, electrochemical techniques with atomic resolution have greatly increased opportunities for observing, measuring, and understanding the surface and interface chemistry during energy conversion. In the future, it is anticipated that the development of electrochemical probing techniques will be advanced with innovative perspectives on the behavior and features of SACs. We hope that this Account can contribute in several ways to promoting the fundamental knowledge and technical progress of emerging electrochemical measurements for studying SACs.

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

Feld-induced modulation of two-dimensional electron gas at LaAlO 3 /SrTiO 3 interface by polar distortion of LaAlO 3

Since the discovery of two-dimensional electron gas at the LaAlO 3 /SrTiO 3 interface, its intriguing physical properties have garnered significant interests for device applications. Yet, understanding its response to electrical stimuli remains incomplete. Our in-situ transmission electron microscopy analysis of a LaAlO 3 /SrTiO 3 two-dimensional electron gas device under electrical bias reveals key insights. Inline electron holography visualized the field-induced modulation of two-dimensional electron gas at the interface, while electron energy loss spectroscopy showed negligible electromigration of oxygen vacancies. Instead, atom-resolved imaging indicated that electric fields trigger polar distortion in the LaAlO 3 layer, affecting two-dimensional electron gas modulation. This study refutes the previously hypothesized role of oxygen vacancies, underscoring the lattice flexibility of LaAlO 3 and its varied polar distortions under electric fields as central to two-dimensional electron gas dynamics. These findings open pathways for advanced oxide nanoelectronics, exploiting the interplay of polar and nonpolar distortions in LaAlO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2018

Spatially controlled doping of two-dimensional SnS 2 through intercalation for electronics

Doped semiconductors are the most important building elements for modern electronic devices. In silicon-based integrated circuits, facile and controllable fabrication and integration of these materials can be realized without introducing a high-resistance interface. Besides, the emergence of two-dimensional (2D) materials enables the realization of atomically thin integrated circuits. However, the 2D nature of these materials precludes the use of traditional ion implantation techniques for carrier doping and further hinders device development10. Here, we demonstrate a solvent-based intercalation method to achieve p-type, n-type and degenerately doped semiconductors in the same parent material at the atomically thin limit. In contrast to naturally grown n-type S-vacancy SnS 2 , Cu intercalated bilayer SnS 2 obtained by this technique displays a hole field-effect mobility of ~40 cm 2 V -1 s -1 , and the obtained Co-SnS 2 exhibits a metal-like behaviour with sheet resistance comparable to that of few-layer graphene. Combining this intercalation technique with lithography, an atomically seamless p–n–metal junction could be further realized with precise size and spatial control, which makes in-plane heterostructures practically applicable for integrated devices and other 2D materials. Therefore, the presented intercalation method can open a new avenue connecting the previously disparate worlds of integrated circuits and atomically thin materials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Clustering on Magnesium Surfaces – Formation and Diffusion Energies

The formation and diffusion energies of atomic clusters on Mg surfaces determine the surface roughness and formation of faulted structure, which in turn affect the mechanical deformation of Mg. This paper reports first principles density function theory (DFT) based quantum mechanics calculation results of atomic clustering on the low energy surfaces {0001} and {$\bar{1}$011} . In parallel, molecular statics calculations serve to test the validity of two interatomic potentials and to extend the scope of the DFT studies. On a {0001} surface, a compact cluster consisting of few than three atoms energetically prefers a face-centered-cubic stacking, to serve as a nucleus of stacking fault. On a {$\bar{1}$011} , clusters of any size always prefer hexagonal-close-packed stacking. Adatom diffusion on surface {$\bar{1}$011} is high anisotropic while isotropic on surface (0001). Three-dimensional Ehrlich–Schwoebel barriers converge as the step height is three atomic layers or thicker. FInally, adatom diffusion along steps is via hopping mechanism, and that down steps is via exchange mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Facet-selective morphology-controlled remote epitaxy of ZnO microcrystals via wet chemical synthesis

We report on morphology-controlled remote epitaxy via hydrothermal growth of ZnO micro- and nanostructure crystals on graphene-coated GaN substrate. The morphology control is achieved to grow diverse morphologies of ZnO from nanowire to microdisk by changing additives of wet chemical solution at a fixed nutrient concentration. Although the growth of ZnO is carried out on poly-domain graphene-coated GaN substrate, the direction of hexagonal sidewall facet of ZnO is homogeneous over the whole ZnO-grown area on graphene/GaN because of strong remote epitaxial relation between ZnO and GaN across graphene. Atomic-resolution transmission electron microscopy corroborates the remote epitaxial relation. The non-covalent interface is applied to mechanically lift off the overlayer of ZnO crystals via a thermal release tape. The mechanism of facet-selective morphology control of ZnO is discussed in terms of electrostatic interaction between nutrient solution and facet surface passivated with functional groups derived from the chemical additives.

36 MATERIALS SCIENCE↗