SEARCH · Search NASA

Results for “charge patterning”

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Two-frequency RF fields induced multipactor in coaxial transmission lines

This study presents a comprehensive investigation of two-surface multipactor discharge in coaxial transmission lines under two-frequency radio frequency (RF) excitation using one-dimensional Monte Carlo simulations validated against three-dimensional particle-in-cell simulations and experimental data. The results show that introducing a second carrier mode can suppress multipactor by reshaping and shrinking the susceptibility region, with the extent and location of suppression strongly dependent on the device's aspect ratio and the relative phase of the second carrier mode. Distinct suppression patterns are observed across different fd regimes, while in some cases, susceptibility expansion also occurs under two-frequency operation. A key outcome is the identification and delineation of pure and mixed multipactor modes in coaxial geometry, where analytical mode boundaries are not readily available. Unlike planar geometries, pure-mode regions in coaxial systems overlap with mixed-mode domains, complicating mode identification. Additionally, image charge forces are found to have negligible influence on susceptibility thresholds but strongly affect electron growth rates. These findings offer valuable insights into the use of waveform engineering for controlling multipactor in high-power RF systems.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Tethered from the Head and from the Tail: The Structure of Hydroxyl-Functionalized Ionic Liquids

Ionic liquids with special functionalities are synthesized with the specific purpose of creating new patterns of interaction in the condensed phase. This Letter discusses the case of alcohol-functionalized ILs, the so-called HFILs, which are part of the larger cohort of task-specific ionic liquids. We find that this small chemical modification can cause massive changes in the liquid landscape when the cationic tails are longer. For prototypical ionic liquids, larger alkyl tails act as separators of charge networks, but in the case of HFILs these become physical charge network linkers. The OH functionality adds a large repertoire of interactions and correlations that were mostly unavailable to traditional ILs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Modeling Electrodeposition in 3D Porous Architectures for Solid-State Li-Metal Batteries

Li-metal storage in three-dimensional (3D) electrodes is considered a potential dendrite-mitigation strategy. The large surface area and high porosity of these electrodes result in reduced local Li-plating current densities. The porous topology provides a scaffold for Li-deposition and stripping, maintaining both mechanical integrity and Li accessibility. The goal of this study is to understand how characteristics, such as geometry and material properties, affect the current distribution and deposition pattern. To this end, we developed a computational method to track material growth driven by electrodeposition within a complex geometry. This method ensures that the finite-element discretization remains conforming to the moving boundary while preserving an adequate mesh quality, and thus maintains solution accuracy. Using this new computational tool, we analyze the conditions under which porous anode architectures effectively expand the surface area of the charge-transfer interface, and self-regulate current density and dendrite growth.

3D electrode architectures↗

DOE OSTI2025

Revealing rotational symmetry breaking charge density wave order in the kagome superconductor (Rb, K)⁢V 3 ⁢Sb 5 by ultrafast pump-probe experiments

The recently discovered Kagome superconductor 𝐴⁢V 3 ⁢Sb 5 (where 𝐴 refers to K, Rb, Cs) has stimulated widespread research interest due to its interplay of nontrivial topology and unconventional correlated physics including charge-density waves (CDW) and superconductivity. The essential prerequisite to understanding the microscopic mechanisms of this complex electronic landscape is to unveil the configuration and symmetry of the charge-density wave order. As to now, little consensus has been made on what symmetry is broken. Herein, we clarify the microscopic structure and symmetry breaking of the CDW phase in RbV 3 ⁢Sb 5 and KV 3 ⁢Sb 5 by ultrafast time-resolved reflectivity. Our approach is based on extracting coherent phonon spectra induced by three-dimensional CDW and comparing them to calculated phonon frequencies via density-functional theory. The combination of these experimental results and calculations provides compelling evidence that the CDW structure of both compounds prevailing up to 𝑇 CDW is the 2 × 2 × 2 staggered inverse Star-of-David pattern with interlayer 𝜋 phase shift, in which the sixfold rotational symmetry is broken. Finally, these observations thus corroborate sixfold rotational symmetry breaking throughout the CDW phase of RbV 3 ⁢Sb 5 and KV 3⁢ Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Solid-State Side-Chain Functionalization of Conjugated Polymers for Expanded Chemical and Functional Versatility

Conjugated polymers that integrate diverse chemical functionalities with high electronic performance are essential for advanced organic electronic, optoelectronic, and biointerfaced technologies. Achieving such multifunctionality typically requires grafting functional units onto polymer side chains. However, conventional solution-phase approaches remain constrained by solubility limitations, side-chain-induced packing disruptions, and challenges in preserving charge transport. Here, we introduce a solid-state side-chain functionalization strategy that exploits a swellable polar side-chain architecture and azide–alkyne click chemistry to enable efficient molecular diffusion and reaction throughout predeposited polymer thin films. This approach substantially broadens the chemical compatibility of graftable units and mitigates the adverse effects on the electrical property. Moreover, the method supports spatially selective functionalization within a continuous film, enabling the patterned incorporation of chemically distinct groups to produce spatially defined optical properties. Here, this solid-state strategy thus provides a versatile platform for constructing conjugated polymers with expanded chemical versatility, preserved electronic performance, and programmable spatial functionality.

Conjugated polymers↗

DOE OSTI2024

Directing Charge Carriers and Ferroelectric Domains at Lateral Interfaces in van der Waals Heterostructures

Emergent phenomena in traditional ferroelectrics are frequently observed at heterointerfaces. Accessing such functionalities in van der Waals ferroelectrics requires the formation of layered heterostructures, either vertically stacked (similar to oxide ferroelectrics) or laterally stitched (without equivalent in 3D-crystals). Here, we investigate lateral heterostructures of the ferroelectric van der Waals semiconductors SnSe and SnS. A two-step process produces ultrathin crystals comprising an SnSe core laterally joined to an SnS edge-band, as confirmed by Raman spectroscopy, transmission electron microscopy (TEM) imaging, and electron diffraction. TEM shows a moiré pattern across the SnSe core due to coverage by an ultrathin SnS layer. The ability of the lateral interface (IF) to direct excited carriers, probed by cathodoluminescence, shows electron transfer over 560 nm diffusion length from the SnS edge-band. Large, thin flakes supporting ferroelectricity allow investigating domains and domain wall interactions in uniform crystals and lateral heterostructures. Polarized optical microscopy of sub-20 nm flakes consistently shows ⟨110⟩ oriented stripe domains with mirror-twin domain walls. Heterostructures adopt two domain configurations, with domains either constrained to the SnSe core or propagating across the entire SnSe–SnS flakes. Furthermore, the combined results demonstrate multifunctional van der Waals heterostructures with high-quality IFs presenting extraordinary opportunities for manipulating carrier flows and ferroelectric domain patterns.

2D ferroelectrics↗

DOE OSTI2026

Time-Resolved Neutron Imaging for Hydrogen Uptake in Subsurface Lithologies

Geologic hydrogen production and underground storage are increasingly important for meeting rising energy demands while providing clean-combustion advantages. However, hydrogen’s high diffusivity and propensity for leakage through porous media necessitate direct evaluation of its transport behavior in subsurface materials. Whereas X-ray microcomputed tomography (μCT) studies often employ contrast agents or surrogate gases, this study leverages neutron transmission radiography/CT to observe hydrogen migration in situ. This work represents the first demonstration of real-time neutron radiography of hydrogen migration in reservoir and caprock lithologies. Cylindrical cores of Indiana limestone, Amherst Gray sandstone, and Tumey shale were subjected to constant-pressure hydrogen charging and scanned in real time using high-resolution neutron radiography. Results indicate immediate hydrogen infiltration in sandstone and limestone, with homogeneous distribution detected throughout their pore structure. In contrast, hydrogen remained largely absent from fine-grained shale under the same pressure, except in an apparently localized fracture zone, where neutron signatures confirmed the presence of hydrogen. Subsequent neutron CT of the sandstone sample, using image subtraction against an uncharged reference, corroborated hydrogen distribution patterns. Even under lowpressure, single-phase conditions, distinct neutron imaging signatures of hydrogen were achieved. These preliminary findings underscore the potential of neutron imaging for advancing subsurface hydrogen migration research.

08 HYDROGEN↗

DOE OSTI2026

Towards dislocation-driven quantum interconnects

A central problem in the deployment of quantum technologies is the realization of robust architectures for quantum interconnects. We propose to engineer interconnects in semiconductors and insulators by patterning spin qubits at dislocations, thus forming quasi one-dimensional lines of entangled point defects. To gain insight into the feasibility and control of dislocation-driven interconnects, we investigate the optical cycle and coherence properties of nitrogen-vacancy (NV) centers in diamond, in proximity of dislocations, using a combination of advanced first-principles calculations. We show that one can engineer spin defects with properties similar to those of their bulk counterparts, including charge stability and a favorable optical cycle, and that NV centers close to dislocations have much improved coherence properties. Finally, we predict optically detected magnetic resonance spectra that may facilitate the experimental identification of specific defect configurations. Our results provide a theoretical foundation for the engineering of one-dimensional arrays of spin defects in the solid state.

Materials science↗

DOE OSTI2024

Proton rings from late-forming ballistic sheath fields

Many laser-driven ion experiments have seen ring-like patterns in the proton angular distribution across a wide array of laser and target parameters. These rings can impede measurement due to the small acceptance angle of detectors and often inhibit potential applications. A myriad of explanations for their formation have been proposed, yet most studies attribute them to some aspect of the laser–plasma interaction. Using 3D particle-in-cell simulations, we show that late-forming strong radial electric fields can arise due to charge separation while the beam is in flight, long after the laser–plasma interaction. These fields can accelerate ions to significant divergences (≈ 10°) as they propagate away from the target. We compare our results to a recent experiment where a high intensity, short pulse laser (I 0 ≈ 10 21 W/cm 2 , τ ≈ 30 fs) was incident upon thin (≈ 1 μm) liquid crystal targets. Our simulations capture all the main features of the experimental results—namely, robust ring formation and larger rings for higher energy protons. In addition, we show that rings do not form for sufficiently short preplasma scale lengths. Finally, we develop a phenomenological model to describe the spatiotemporal structure of the radial electric field and use this to explain the proton rings' energy and preplasma dependencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Compositionally Complex Spinel Oxides as Conversion Anodes for Lithium-Ion Batteries

Four different compositionally complex multicomponent M 3 O 4 spinels containing 5–8 distinct metals were prepared by a rapid combustion synthesis method or solvothermal synthesis. High resolution synchrotron X-ray diffraction patterns show that the materials consist primarily of spinel phases with small amounts of rock salt impurities, and, in several samples, a minor amount of contracted spinel phase. Materials were investigated as conversion anodes in lithium half-cells and delivered significantly higher capacities than two-component MgFe 2 O 4 made by combustion synthesis. X-ray absorption near-edge structure (XANES) was used to estimate the oxidation states of the metals in the pristine, lithiated (discharged) and delithiated (charged) materials to better understand the redox processes in half cells that led to the improvement. Co, Ni, and Zn are reduced to low oxidation states during lithiation (cell discharge) but are only partially oxidized. The presence of a conductive metallic network that forms after lithiation is thought to account for the improved electrochemical characteristics. Interestingly, in most of the samples, iron is not fully reduced during initial lithiation unlike what happens with a set of related high entropy spinel ferrites studied previously. Finally, the improved electrochemical properties of these materials illustrates both the advantages of complexity and the difficulties in predicting their behavior.

25 ENERGY STORAGE↗

DOE OSTI2026

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Evolution and Degradation Patterns of Electrochemical Cells Based on the Analysis of Interfacial Phenomena at Li Metal Anode/Electrolyte Interfaces

In this work, we report the results of a theoretical–computational analysis of the solid electrolyte interphase (SEI) growth and degradation dynamics occurring in lithium metal batteries during cycling. We use ab initio-kinetic Monte Carlo simulations to generate a synthetic data set, which is analyzed by machine learning methods. We aim to determine: (i) how modifications in interfacial interaction energies between solid electrolyte interphase (SEI) blocks and between Li ions and SEI facets impact the Coulombic efficiency (CE) of the battery and (ii) what factors, including reactions, microscopic transport, and other interfacial events, may lead to cell performance “failure” during prolonged charge and discharge cycles, signaled as a sharp decay in the CE over cycling. The demonstration of our approach is done on a cell including a Li metal surface interfacing with a previously introduced state-of-the-art electrolyte, and the idea can be applied to any electrochemical system. Outcomes include the identification of the leading chemical, physical, and structural variables causing cell failure and relating them to the electrolyte formulation, thus paving the way to future more refined analysis and electrolyte design.

batteries↗

DOE OSTI2024

The power reliability event simulator tool (PRESTO): A novel approach to distribution system reliability analysis and applications

The growing interest in onsite solar photovoltaic and energy storage systems is partially motivated by customer concerns regarding grid reliability. However, accurately assessing the effectiveness of PVESS in mitigating these interruptions requires a comprehensive understanding of location-specific outage patterns and the ability to simulate realistic scenarios. To address the gap, we introduce the Power Reliability Event Simulation TOol (PRESTO), the first publicly available tool that simulates location-specific power interruptions at the county level. PRESTO allows for a more realistic assessment of system reliability by considering the unpredictability and location-specific patterns of power interruptions. We applied PRESTO in a case study of a single-family home across three U.S. counties, examining the performance of a solar photovoltaic system with 10kWh of battery storage during short-duration power interruptions. Our findings show that this system reliably met 93% of energy demand for essential non-heating and cooling loads, fully serving these loads in 84% of events, despite the constraints of daily time-of-use bill management which limits the battery's state-of-charge reserve. However, when heating and cooling loads were included, system performance decreased significantly, with only 70% of demand met and full service in 43% of events. These results highlight the challenges of using solar photovoltaic and energy storage systems for short-duration outages, emphasizing the need to consider factors like battery size and grid charging strategies to improve reliability. Our study demonstrates the practical applications of PRESTO, providing valuable insights into potential mitigation strategies including grid charging and optimizing battery size.

14 SOLAR ENERGY↗

DOE OSTI2025

Exact results in softly broken supersymmetric chiral gauge theories with flavor

We present exact results in softly broken supersymmetric SU⁡(𝑁 𝐶 ) chiral gauge theories with charged fermions in one antisymmetric, 𝑁 𝐹 fundamental, and 𝑁 𝐶 + 𝑁 𝐹 − 4 antifundamental representations. We achieve this by considering the supersymmetric version of these theories and utilizing anomaly mediated supersymmetry breaking at a scale 𝑚 ≪ Λ to generate a vacuum. The connection to nonsupersymmetric theories is then conjectured in the limit 𝑚 → ∞. For odd 𝑁 𝐶 , we determine the massless fermions and unbroken global symmetries in the infrared. For even 𝑁 𝐶 , we find global symmetries are nonanomalous and no massless fermions. In all cases, the symmetry breaking patterns differ from what the tumbling hypothesis would suggest.

Anomalies↗

DOE OSTI2025

Glancing Angle Deposition in Gas Sensing: Bridging Morphological Innovations and Sensor Performances

Glancing Angle Deposition (GLAD) has emerged as a versatile and powerful nanofabrication technique for developing next-generation gas sensors by enabling precise control over nanostructure geometry, porosity, and material composition. Through dynamic substrate tilting and rotation, GLAD facilitates the fabrication of highly porous, anisotropic nanostructures, such as aligned, tilted, zigzag, helical, and multilayered nanorods, with tunable surface area and diffusion pathways optimized for gas detection. This review provides a comprehensive synthesis of recent advances in GLAD-based gas sensor design, focusing on how structural engineering and material integration converge to enhance sensor performance. Key materials strategies include the construction of heterojunctions and core–shell architectures, controlled doping, and nanoparticle decoration using noble metals or metal oxides to amplify charge transfer, catalytic activity, and redox responsiveness. GLAD-fabricated nanostructures have been effectively deployed across multiple gas sensing modalities, including resistive, capacitive, piezoelectric, and optical platforms, where their high aspect ratios, tailored porosity, and defect-rich surfaces facilitate enhanced gas adsorption kinetics and efficient signal transduction. These devices exhibit high sensitivity and selectivity toward a range of analytes, including NO2, CO, H2S, and volatile organic compounds (VOCs), with detection limits often reaching the parts-per-billion level. Emerging innovations, such as photo-assisted sensing and integration with artificial intelligence for data analysis and pattern recognition, further extend the capabilities of GLAD-based systems for multifunctional, real-time, and adaptive sensing. Finally, current challenges and future research directions are discussed, emphasizing the promise of GLAD as a scalable platform for next-generation gas sensing technologies.

Chemistry↗

DOE OSTI2025

Thermodynamic Theory of Proximity Ferroelectricity

Proximity ferroelectricity has recently been reported as a new design paradigm for inducing ferroelectricity, where a nonferroelectric polar material becomes a ferroelectric one by interfacing with a thin ferroelectric layer. Strongly polar materials, such as AlN and ZnO, which were previously unswitchable with an external field below their dielectric breakdown fields, can now be switched with practical coercive fields when they are in intimate proximity to a switchable ferroelectric. Here, we develop a general Landau-Ginzburg theory of proximity ferroelectricity in multilayers of nonferroelectrics and ferroelectrics to analyze their switchability and coercive fields. The theory predicts regimes of both “proximity switching,” where the multilayers collectively switch, and “proximity suppression,” where they collectively do not switch. The mechanism of the proximity ferroelectricity is an internal electric field determined by the polarization of the layers and their relative thickness in a self-consistent manner that renormalizes the double-well ferroelectric potential to lower the steepness of the switching barrier. Further reduction in the coercive field emerges from charged defects in the bulk that act as nucleation centers. The application of the theory to proximity ferroelectricity in Al x−1⁢ Sc x ⁢N/AlN and Zn 1−x ⁢Mg x ⁢O/ZnO bilayers is demonstrated. The theory further predicts that dielectric-ferroelectric and paraelectric-ferroelectric multilayers can potentially lead to induced ferroelectricity in the dielectric or paraelectric layers, resulting in the entire stack being switched, an exciting avenue for new discoveries. This thawing of “frozen ferroelectrics,” paraelectrics, and potentially dielectrics with high dielectric constants promises a large class of new ferroelectrics with exciting prospects for previously unrealizable domain-patterned optoelectronic and memory technologies.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Improving ICARUS Track Reconstruction Algorithms

The ICARUS experiment is part of the Short-Baseline Neutrino (SBN) program at Fermilab. The main goal of the experiment is to investigate the possibility of sterile neutrinos in the O(1 eV) mass region and provide clarification of the anomaly detected from the Liquid Scintillator Neutrino Detector (LSND) and MiniBooNE experiments. The ICARUS-T600 detector is a Liquid Argon Time Projection Chamber (LAr-TPC), that can provide excellent 3D imaging and calorimetric reconstruction of any ionizing particles. This detection technique allows a detailed study of neutrino interactions, spanning a wide energy spectrum (from a few keV to several hundreds of GeV). The detector consists of two identical adjacent modules, filled with a total of 760 tons of ultra-pure liquid argon. Each module houses two LAr-TPCs separated by a common cathode with a maximum drift distance of 1.5 m, equivalent to about 1 ms drift time for the nominal $500$ V/m electric drift field. The anode is made of three parallel wire planes positioned 3 mm apart, where the stainless-steel wires are oriented on each plane at a different angle with respect to the horizontal direction ($+60^\degree$,$-60^\degree$,$0^\degree$). The first two planes (Induction 1 and Induction 2) provide a non-destructive charge measurement, whereas the ionization charge is fully collected by the last collection plane. In total, 53248 wires with a 3 mm pitch and length up to 9 m are installed in the detector. In the first stage of the reconstruction, segments of waveforms corresponding to physical signals (hits) are searched for in the deconvolved wire waveform with a threshold-based hit-finding algorithm. Each hit is then fitted with a Gaussian, whose area is proportional to the number of drift electrons generating the signal. In the second stage of the reconstruction, hits are passed as input to Pandora, a framework software composed of different pattern recognition algorithms, that performs a 3D reconstruction of the full image recorded in the collected event, including the identification of interaction vertices and tracks and showers inside the TPC. These are organized into a hierarchical structure (called slice) of particles generated starting from a primary interaction vertex. In some cases, related to the inefficiencies in the hit detection or excessive deflection of the particle trajectory, Pandora breaks the particle's track into two or more smaller pieces and considers each piece as an independent track. We studied this phenomenon focusing on primary muons from ν_μ CC interactions contained in a single module with a track at least 20 cm long, to exclude delta rays. The study determined that about $7-8\%$ of the muon tracks are broken. Approximately $80\%$ of the times, Pandora assigns all segments of the track to the same slice (intra-slice track split), while in the remaining $20\%$ of the cases, one of the segments is associated with another slice (extra-slice track split). To mitigate this phenomenon, we designed an algorithm that detects and stitches the tracks broken by Pandora for the intra-slice split. In Monte Carlo simulations, the algorithm showed an efficiency exceeding $80\%$ and a purity exceeding $93\%$.

Ricci, Alessandro Maria [Pisa U.; INFN, Pisa] (ORC↗

DOE OSTI2026

Mechanism-Informed Breakdown: Understanding Degradation by Controlling Voltage-Hold Patterns in Proton Exchange Membrane Water Electrolyzers

Low catalyst loadings pose challenges to performance stability in proton exchange membrane (PEM) water electrolysis over extended operation. To study the impact of degradation mechanisms and voltage loss rates, different stress tests are applied to membrane electrode assemblies. Potential cycling conditions were observed to induce higher degrees of iridium (Ir) oxide crystallization, ionomer degradation, and catalyst layer (CL) thinning, which likely contributed to higher kinetic loss rates. On the other hand, while Ir migrating into the PEM (Ir band) generally impairs performance, the interconnected and more uniform Ir band formed under a constant 2 V hold may allow for Ir at the catalyst/membrane interface to remain electronically connected and kinetically accessible, as well as indicate greater Ir site access during the applied stressor. The 2 V hold also demonstrates improved kinetic durability through a lower Tafel slope, faster polarization kinetics, and reduced charge transfer resistance. In contrast, potential cycling caused the migration of disconnected Ir agglomerates into the membrane bulk and created a steady increase in charge transfer resistance, a more dramatic decrease in capacitance (46.7% loss), and significant damage to the surrounding ionomer, indicating a decline in both the quality and quantity of active sites in the anode CL. This work underscores the distinct degradation pathways associated with load holds versus cycling, highlighting the role of catalyst-ionomer interactions in kinetic performance and long-term stability. These insights can inform operational strategies for PEM electrolyzers powered by intermittent energy sources, aiming to minimize efficiency losses over extended operation.

36 MATERIALS SCIENCE↗