Search NASA⌕ Search

SEARCH · Search NASA

Results for “First Principles”

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 217 records · Page 12

Boiler Health Monitoring Using a Hybrid First Principles-Artificial Intelligence Model

Due to increased penetration of the intermittent renewables to the grid, pulverized coal (PC) plants are being forced to cycle their load frequently and rapidly, operate at low load condition for sustained period, and start up and shut down several hundred times in a year in the worst case. These severe operations are causing substantial damage to the boiler components compromising the reliability of PC plants. An online health monitoring tool can be instrumental in understanding the impacts of load-following and can eventually help PC plants to develop advanced process control strategies for improved flexibility without compromising safety nor reliability.

20 FOSSIL-FUELED POWER PLANTS↗

Search for active and inactive ion insertion sites in organic crystalline materials

The position of mobile active and inactive ions, specifically ion insertion sites, within organic crystals, significantly affects the properties of organic materials used for energy storage and ionic transport. Identifying the positions of these atomic (and ionic) sites in an organic crystal is challenging, especially when the element has low X-ray scattering power, such as lithium (Li) and hydrogen, which are difficult to detect by powder X-ray diffraction. First-principles calculations, exemplified by density functional theory (DFT), are very practical for confirming the relative stability of ion positions in materials. However, the lack of effective strategies to identify ion sites in these organic crystalline frameworks renders this task extremely challenging. This work presents two algorithms: the (i) efficient location of ion insertion sites from extrema in electrostatic local potential and charge density (ELIISE), and the (ii) ElectRostatic InsertioN (ERIN), which leverage charge density and electrostatic potential fields accessed from first-principles calculations, combined with the Simultaneous Ion Insertion and Evaluation (SIIE) workflow. SIIE inserts all ions simultaneously—to determine ion positions in organic crystals. We demonstrate that these methods accurately reproduce known ion positions in 16 organic materials and identify previously overlooked low-energy sites in tetralithium 2,6-naphthalenedicarboxylate (Li4NDC), an organic electrode material, highlighting the importance of inserting all ions simultaneously, as in the SIIE workflow. These algorithms are also integrated with off-the-shelf machine learning potentials, yielding promising results comparable to first-principles findings.

Gopidi, Harshan Reddy [Argonne National Laboratory↗

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cost analysis of alternative large-scale high-temperature solid oxide electrolysis hydrogen production facilities

We extend our past cost analysis of gigawatt-scale solid oxide electrolysis (SOE) facilities that produce high purity hydrogen gas from water by estimating construction and operating costs for three new alternative design cases: (1) offsite feed steam generation; (2) near-atmospheric pressure (NAP) stack; and (3) onsite electric boiler feed steam generation. Pressure effects on hydrogen electrode-(cathode-)supported SOE cell (SOEC) stack performance are estimated for the same assumed cell and stack construction and used to determine facility-wide stack capital costs for achieving a fixed H2 production at different pressures. Modular balance of plant (BOP) process equipment capital costs are estimated for each new alternative design case using our past equipment sizing, design, and cost data and scaling relationships. Furthermore, we update BOP equipment sizing and design for the NAP case using Aspen®. Vendor quotes for electric boilers are used to estimate costs for the electric boiler design case. Factory and onsite assembly and installation costs for SOEC stacks and BOP equipment are calculated using our past simplified first-principles approach. First-of-a-kind (FOAK) and N th -of-a-kind (NOAK) production maturity cost estimates are included for all results. The case with NAP stacks offers the lowest facility total capital cost (TCC, ~23% lower than base) while use of small electric boilers requires the highest TCC (~3% higher than base). H 2 production prices decrease from the base of ~$\$2.17$ /kgH 2 to ~$\$1.92$/kgH 2 for 1 GW e DC SIP facilities utilizing NAP stacks supplied by offsites steam situated in large modules and blocks for $\$0.030$/kWh e and $0.009/kWh t prices for electricity and thermal energy, respectively. We report all costs in 2021 US dollars.

Balance of plant (BOP) process equipment↗

Near-room-temperature ferromagnetic ordering in the pressure-induced collapsed-tetragonal phase in SrC⁢o 2 ⁢P 2

We present high-pressure electrical transport, magnetization, and single-crystal x-ray diffraction data on SrC⁢o 2 ⁢P 2 single crystals, revealing a transition from the uncollapsed tetragonal (ucT) to the collapsed tetragonal (cT) phase >∼10 GPa. This structural transition is accompanied by ferromagnetic (FM) ordering with a Curie temperature up to 260 K. First-principles-based density functional theory calculations show that there is a first-order transition between ucT and cT phases, with an onset ∼10 GPa as well as the appearance of FM ordering in the cT phase. Above ∼30 GPa, the experimental signatures of the magnetic ordering vanish in a first-order-like manner, consistent with the theoretical calculation results, indicating that SrC⁢o 2 ⁢P 2 is another example of the avoidance of quantum criticality in FM intermetallic compounds with almost room-temperature critical temperature. SrC⁢o 2 ⁢P 2 provides clear evidence that the structural, electronic, and magnetic properties associated with the cT transition are strongly entangled and are not only qualitatively captured by our first-principles-based calculations but are quantitatively reproduced as well.

density functional theory↗

Is Ba 3 In 2 O 6 a high- T c superconductor?

It has been suggested that Ba 3 In 2 O 6 might be a high-Tc superconductor. Experimental investigation of the properties of Ba 3 In 2 O 6 was long inhibited by its instability in air. Recently epitaxial Ba 3 In 2 O 6 with a protective capping layer was demonstrated, which finally allows its electronic characterization. The optical bandgap of Ba3In2O6 is determined to be 2.99 eV in-the (001) plane and 2.83 eV along the c-axis direction by spectroscopic ellipsometry. First-principles calculations were carried out, yielding a result in good agreement with the experimental value. Various dopants were explored to induce (super-)conductivity in this otherwise insulating material. Neither A- nor B-site doping proved successful. The underlying reason is predominately the formation of oxygen interstitials as revealed by scanning transmission electron microscopy and first-principles calculations. Additional efforts to induce superconductivity were investigated, including surface alkali doping, optical pumping, and hydrogen reduction. To probe liquid-ion gating, Ba 3 In 2 O 6 was successfully grown epitaxially on an epitaxial SrRuO 3 bottom electrode. So far none of these efforts induced superconductivity in Ba 3 In 2 O 6 , leaving the answer to the initial question of whether Ba 3 In 2 O 6 is a high-T c superconductor to be 'no' thus far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of interparticle fields and radiation reaction on beam dynamics

The dynamics of relativistic particles in an intense electromagnetic field can be described by the Landau-Lifshitz (LL) equation, where the radiation reaction (RR) is accounted for via a self-force, and interparticle fields are often neglected as an approximation. However, the inclusion of interparticle fields is necessary to ensure energy-momentum conservation, particularly during coherent emission. Here we present (i) an analytical proof showing that the energy-momentum conservation law of the Hamilton-Rohrlich-Dirac action, which is divergence free and describes a generic system of interacting charges, respects causality and provides physically sensible results; (ii) a simple generalization of the LL equation for many particles evaluated as a function of the total field, i.e., the sum of the external and interparticle fields. By performing first-principles numerical simulations of a neutral, relativistic bunch of electrons and positrons (e − /e + ) colliding with a laser pulse, this theory is shown to satisfy energy-momentum conservation when interparticle fields and RR are simultaneously taken into account; and (iii) the combined effect of interparticle fields and RR primarily affects the tail of the particle energy distribution. Additionally, our first-principles simulations show that the effect of interparticle fields on beam energy loss becomes smaller when most of the radiated energy is incoherent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fermiology and transport properties of the candidate topological crystalline insulator SrAg 4 ⁢Sb 2

Compared to time-reversal symmetry-protected ℤ 2 topological insulators and Dirac/Weyl semimetals, there are significantly fewer candidates for topological crystalline insulators. SrAg 4 ⁢Sb 2 is predicted to exhibit topological crystalline insulator behavior when considering spin-orbit coupling. In this study, we systematically investigate single crystals of SrAg 4 ⁢Sb 2 using electrical transport and magnetic torque measurements, along with first-principles calculations. Our transport data reveals its compensated semimetal nature with a magnetoresistance up to around 700% at 2 K and 9 T. Analysis of de Haas–van Alphen oscillations uncovers a Fermi surface consisting of three distinct Fermi pockets with light effective masses. Comparison between the three-dimensional fermiology obtained from our oscillation data and the first-principles calculations demonstrates excellent agreement. This confirms the accuracy of the calculations, which indicate a band inversion centered at the 𝑇 point and identify the existence of nontrivial tube and needle hole Fermi pockets at Γ, alongside one trivial diamond electron pocket at the 𝐹 point in the Brillouin zone. Furthermore, symmetry and topology analysis results in two potential sets of topological invariants, suggesting the emergence of two-dimensional gapless Dirac surface states either on the 𝑎⁢𝑏 planes or on both the 𝑎⁢𝑏 planes and mirror planes, protected by crystal symmetries. Therefore, SrAg 4 ⁢Sb 2 emerges as a promising candidate topological crystalline insulator.

36 MATERIALS SCIENCE↗

Chromium segregation-induced oxide evolution in Ni-10Cr alloys during high-temperature oxidation

The oxidation behavior of a Ni-10(wt%)Cr alloy under high-temperature O 2 conditions is investigated using transmission electron microscopy and first-principles calculations. Results reveal that chromium segregation plays a central role in driving the evolution of complex oxide phase structures during oxidation. At low Cr concentrations, Cr preferentially segregates to NiO grain boundaries or internal pores, substituting for Ni atoms and forming Ni(Cr)O solid solutions. As Cr content increases, enhanced diffusion promotes Cr penetration into the NiO lattice, leading to the formation of multiphase oxide structures. First-principles modeling corroborates these findings: at low Cr concentrations, Cr atoms favor surface and grain-boundary segregation, while higher concentrations lead to Cr aggregation within the NiO bulk. Furthermore, the integrated experimental-theoretical approach provides atomistic insights into Cr-mediated mass transport mechanisms during alloy oxidation and offers valuable guidance for controlling oxide growth kinetics and phase stability in Ni-Cr alloys, with implications for improving oxidation resistance in high-temperature structural applications.

36 MATERIALS SCIENCE↗

Ferroelectric Semimetals with α-Bi/SnSe van der Waals Heterostructures and Their Topological Currents

We show that proximity effects can be utilized to engineer van der Waals heterostructures (vdWHs) displaying semimetallic spin-ferroelectricity locking, where ferroelectricity and semimetallic spin states are confined to different layers, but are correlated by means of proximity effects. Our findings are supported by first principles calculations involving α-Bi/SnSe bilayers. We show that such systems support ferroelectrically switchable nonlinear anomalous Hall effect originating from large Berry curvature dipoles as well as direct and inverse spin Hall effects with giant bulk spin-charge interconversion efficiencies. In conclusion, the giant efficiencies are consequences of the proximity-induced semimetallic nature of low energy electron states, which are shown to behave as two-dimensional pseudo-Weyl fermions by means of symmetry analysis and first principles calculations as well as direct angle-resolved photoemission spectroscopy measurements.

2-dimensional systems↗

Computationally guided experimental validation of divacancy defect formation in 4H-SiC

Recent research into solid-state qubits for quantum information science has focused on optically addressable spin defects such as the negatively charged nitrogen-vacancy center in diamond and the neutrally charged divacancy (VV) in 4H-SiC as scalable quantum sensors and networking qubits. Within this context, direct investigations of the structural origin and defect formation dynamics of a sub-set of the VV center in 4H-SiC remain lacking. Here, we take a systematic experimental approach guided by predictions from first-principles simulations to gain a thorough mechanistic understanding of the VV defect formation and control in 4H-SiC. We study the effect of annealing time and temperature on VV formation in high-purity semi-insulating 4H-SiC samples following electron irradiation. Three different temperatures (1123, 1273, and 1473 K) and annealing duration (from 0.5 to 72 h) are chosen to explore VV formation in different regions. We find that samples annealed at 1273 K give the highest VV-related photoluminescence (PL) intensities, in agreement with the prediction from first-principles calculations. Furthermore, the logarithmic dependence of VV-related PL intensities on the annealing duration at 1273 K indicates that 1273 K provides sufficient thermal energy for silicon vacancy migration but not for VV migration. Together, these results suggest that efficient VV formation occurs above the V Si migration temperature and below the VV migration threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗

Direct and indirect spin current generation and spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers

Spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers are studied using a combination of symmetry analysis, circuit theory, semiclassical simulations, and first-principles calculations using the nonequilibrium Green's function method with supercell disorder averaging. Here, we focus on unconventional processes involving the interplay between the two ferromagnetic layers, which are classified into direct and indirect mechanisms. The direct mechanism involves spin current generation by one ferromagnetic layer and its subsequent absorption by the other. In the indirect mechanism, the in-plane spin-polarized current from one ferromagnetic layer “leaks” into the other layer, where it is converted into an out-of-plane spin current and reabsorbed by the original layer. The direct mechanism results in a predominantly dampinglike torque, which damps the magnetization towards a certain direction 𝐬 𝑑 . The indirect mechanism results in a predominantly fieldlike torque with respect to a generally different direction 𝐬 𝑓 . Similarly to the current-in-plane giant magnetoresistance, the indirect mechanism is only active if the thickness of the nonmagnetic spacer is smaller than or comparable to the mean free path. Numerical calculations for a semiclassical model based on the Boltzmann equation confirm the presence of both direct and indirect mechanisms of spin current generation. First-principles calculations reveal sizable unconventional spin-orbit torques in Co/Cu/Co, Py/Cu/Py, and Co/Pt/Co trilayers and provide strong evidence of indirect spin current generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cation Vacancies Enable Anion Redox in Li Cathodes

Conventional Li-ion battery intercalation cathodes leverage charge compensation that is formally associated with redox on the transition metal. Employing the anions in the charge compensation mechanism, so-called anion redox, can yield higher capacities beyond the traditional limitations of intercalation chemistry. Here, we aim to understand the structural considerations that enable anion oxidation and focus on processes that result in structural changes, such as the formation of persulfide bonds. Using a Li-rich metal sulfide as a model system, we present both first-principles simulations and experimental data that show that cation vacancies are required for anion oxidation. First-principles simulations show that the oxidation of sulfide to persulfide only occurs when a neighboring vacancy is present. To experimentally probe the role of vacancies in anion redox processes, we introduce vacancies into the Li 2 TiS 3 phase while maintaining a high valency of Ti. When the cation sublattice is fully occupied and no vacancies can be formed through transition metal oxidation, the material is electrochemically inert. Upon introduction of vacancies, the material can support high degrees of anion redox, even in the absence of transition metal oxidation. The model system offers fundamental insights to deepen our understanding of structure–property relationships that govern reversible anion redox in sulfides and demonstrates that cation vacancies are required for anion oxidation, in which persulfides are formed.

25 ENERGY STORAGE↗

Density functional theory-based surrogate kinetic models for heterogeneous reactions of hydrocarbon intermediates on silicon carbide

The increasing demand for high-performance materials in advanced technologies highlights the importance of achieving a fundamental understanding and potential control of silicon carbide (SiC) deposition processes. However, existing models often lack sufficient theoretical detail, relying heavily on empirical data and offering limited predictive capability. In particular, the complex surface chemistry governing SiC growth remains poorly understood. This study addresses these challenges by employing density functional theory (DFT) to investigate key heterogeneous reactions involving hydrocarbon intermediates on SiC surfaces, including dehydrogenation, hydrogenation, and carbon deposition. Transition state searches were conducted to identify reaction pathways and energy barriers. While first-principles calculations offer high accuracy, they are computationally intensive. To extend the utility of these first-principles results, vibrational analyses were performed using phonon-based statistical thermochemistry to compute temperature-dependent reaction rates which were used to develop Arrhenius-type surrogate kinetic models. Furthermore, the resulting framework provides a more rigorous, physically grounded basis for integrating atomistic insights into continuum-scale modeling, ultimately enabling improved prediction and optimization of SiC film growth in high-performance material systems.

Density Functional Theory↗

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗

TiV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys titanium-vanadium (Ti-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ti and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗