Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster ions”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗
DOE OSTI2023

Tandem Mass-Selective Cryogenic Digital Ion Traps for Enhanced Cluster Formation

We present the implementation of tandem mass-selective cryogenic ion traps, designed to enhance the range of ion processing capabilities that can be performed prior to spectroscopic interrogation. We show that both the formation of ion clusters and mass filtering steps can be combined in a single cryogenic linear quadrupole ion trap driven by RF square waves. Mass filtering and mass isolation can be achieved by manipulation of the RF frequency and duty cycle. Very importantly, this scheme circumvents the need for high-amplitude RF voltages that can be incompatible with typical cryogenic ion processing conditions. In addition, proper adjustment of the stability boundaries during the clustering process allows for the preferential formation of a specific cluster size rather than a broad distribution of sizes. Lastly, we show that a specific cluster size can be formed, mass-selected, and then transferred to another ion trap for a second completely separate ion processing step. Furthermore, the instrumentation and modular design developed here expand the scope of ionic species and clusters that can be accessed by processing electrosprayed ions.

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

Microsolvation in V + (H 2 O) n Clusters Studied with Selected-Ion Infrared Spectroscopy

Gas phase ion-molecule clusters of the form V + (H 2 O) n (n = 1-30) are produced by laser vaporization in a supersonic expansion. These ions are analyzed and mass-selected with a time-of-flight mass spectrometer and investigated with infrared laser photodissociation spectroscopy. The small clusters (n ≤ 7) are studied with argon tagging, while the larger clusters are studied via the elimination of water molecules. The vibrational spectra for the small clusters include only free O-H stretching vibrations, while larger clusters exhibit red-shifted hydrogen bonding vibrations. The spectral patterns reveal that the coordination around V + ions is completed with four water molecules. A symmetric square-planar structure forms for the n = 4 ion, and this becomes the core ion in larger structures. Clusters up to n = 8 have mostly two-dimensional structures, but hydrogen bonding networks evolve to three-dimensional structures in larger clusters. The free O-H vibration of AAD (acceptor-acceptor-donor) coordinated surface molecules converges to a frequency near that of bulk water by the cluster size of n = 30. However, the splitting of this vibration for AAD- versus AD-coordinated molecules is still different compared to other singly-charged or doubly-charged cation-water clusters. Finally, this indicates that cation identity and charge-site location in the cluster can produce discernible spectral differences for clusters in this size range.

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

Anion-dependent phase behavior of methylimidazolium-based ionic liquids mixed with water: Correlation between local molecular structure and mesoscale behaviors

The mesoscopic phase behavior of decylmethylimidazolium (C 10 mim) ionic liquids (ILs) bearing three monovalent anions—thiocyanate (SCN − ), nitrate (NO 3 − ), and chloride (Cl − )—mixed with water at relatively high IL contents (50–95 wt%) was investigated. Small-angle and wide-angle X-ray scattering (SAXS/WAXS) were employed to follow the evolution of both local and mesoscale structures across this composition range. In the absence of water, C 10 mimSCN and C 10 mimCl behaved as disordered liquids, whereas C 10 mimNO 3 spontaneously formed a hexagonally ordered cylindrical mesophase and displayed a sticky-solid macroscopic appearance. Upon addition of water, C 10 mimSCN remained a viscous liquid and only weakly ordered lamellar domains were observed. This limited ordering is attributed to the lack of hydrogen bonding and weak interaction energy between SCN − anions and C 10 mim + cations. In contrast, the trigonal-planar NO 3 − and point-like Cl − anions promoted the formation of well-defined hexagonal mesophases up to 35–45 wt% water. The formation of hydrogen bonding of the two anions with imidazolium ring protons likely enabled the creation of compact ion clusters that effectively exclude water molecules from the immediate vicinity of the IL aggregates. These findings demonstrate that the interaction energy between ion pairs dictate IL-water interactions and therefore control the transition from disordered liquids to ordered mesophases in IL/water mixtures with high IL contents. In conclusion, the combined SAXS/WAXS analysis reveals a correlation between local intermolecular structure and the emergence of mesoscopic order, providing a systematic framework for tailoring mesoscale structures in alkylimidazolium-based IL/water systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Elucidating the Molecular Origins of the Transference Number in Battery Electrolytes Using Computer Simulations

The rate at which rechargeable batteries can be charged and discharged is governed by the selective transport of the working ions through the electrolyte. Conductivity, the parameter commonly used to characterize ion transport in electrolytes, reflects the mobility of both cations and anions. The transference number, a parameter introduced over a century ago, sheds light on the relative rates of cation and anion transport. This parameter is, not surprisingly, affected by cation–cation, anion–anion, and cation–anion correlations. In addition, it is affected by correlations between the ions and neutral solvent molecules. Computer simulations have the potential to provide insights into the nature of these correlations. We review the dominant theoretical approaches used to predict the transference number from simulations by using a model univalent lithium electrolyte. In electrolytes of low concentration, one can obtain a quantitative model by assuming that the solution is made up of discrete ion-containing clusters–neutral ion pairs, negatively and positively charged triplets, neutral quadruplets, and so on. These clusters can be identified in simulations using simple algorithms, provided their lifetimes are sufficiently long. In concentrated electrolytes, more clusters are short-lived and more rigorous approaches that account for all correlations are necessary to quantify transference. Elucidating the molecular origin of the transference number in this limit remains an unmet challenge.

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

Observation of N-rich solid-electrolyte interphase by ToF-SIMS.

Formation of a stable solid electrolyte interphase (SEI) between lithium electrodes and electrolyte upon multiple charge/discharge cycles is crucial to a long-term lithium-ion battery performance. Addition of LiNO3 to lithium bis (fluorosulfonyl) imide/poly(ethylene oxide) (LiFSI/PEO) electrolyte leads to a durable SEI that is electrically insulating yet highly conductive to Li ions, chemically and electrochemically stable, physically uniform, and mechanically robust. ToF-SIMS was used here in combination with sputtering by a gaseous cluster ion beam (GCIB) to examine how the addition of a small proportion of LiNO3 to the LiFSI/PEO electrolyte affects the SEI composition. Negative ion ToF-SIMS spectra of the cycled samples display an intense m/z 26 peak associated with the SEI. Exact mass assignments and isotopic ratios indicate that this peak should be assigned as (CN-)-C-12, with little to no negative secondary ion signal arising from (LiF-)-Li-7. This CN- signal appears to arise from an N-rich portion of the SEI adjacent to the Li electrode that is depleted in LiF relative to the bulk electrolyte. The dearth of LiF- (and LiF+ from the positive ion spectra) is unexpected because LiF has been identified in the SEI in similar samples. Finally, GCIB sputtering indicates that the SEI adheres more strongly to the Li electrode than to the LiFSI/PEO electrolyte.

Shavandi, Seyedeh Reyhaneh↗
DOE OSTI2021

Gas-phase stability of large lanthanide:diglycolamide clusters evaluated using collision-induced dissociation

Efficient utilization of long-term deep geologic storage repositories is critical for the large-scale deployment of nuclear energy. As the minor actinides (americium and curium) are the largest contributor to the decay heat of used nuclear fuel after a few hundred years, their removal from used fuel prior to disposal can significantly increase the amount of waste that can be stored in a given volume. The development of ligands that can effectively separate the minor actinides from other components of used nuclear fuel are crucial to efficient utilization of geological storage repositories. N,N,N’,N’-tetraoctyl diglycolamide (TODGA) is a promising extractant for the separation of lanthanides and the minor actinides from other components of used nuclear fuel. While the use of TODGA in f-element separations has been investigated in process-based formulations, gas-phase metal ion cluster experiments enable the study of covalent interactions in reprocessing systems absent from solvent effects. Exploring fundamental differences in lanthanide-ligand covalent interactions can impact the development and implementation of actinide-lanthanide separation systems in nuclear fuels reprocessing. In this study, lanthanide-TODGA clusters were synthesized in the gas-phase and identified using mass spectrometry fragmentation experiments. Large europium and samarium clusters were identified that contained up to 8 and 10 bound TODGA ligands, respectively; this was surprising due to the size and multidentate binding that is normal for TODGA. In addition, while smaller clusters showed evidence of sequential ligand fragmentation, larger clusters displayed the loss of neutral TODGA with applied collision voltage. Interestingly, the voltage required for this removal decreased as more TODGA ligands were bound to the metal, suggesting that the metal coordination sphere was becoming more saturated and TODGA ligands were more weakly bound as the clusters got larger.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Unraveling disadhesion mechanism of epoxy/CSH interface under aggressive conditions

Almost all FRP-reinforced concrete exposed to the marine environment are susceptible to durability deficiencies due to failure of the FRP/epoxy/concrete interface, especially the vulnerable epoxy/concrete interface, which is frequently attacked by water and ions. Herein, molecular dynamics simulations of epoxy/CSH interface under different environmental exposures are utilized to specifically elucidate the degradation principles of water and ions on the bonding properties and mechanical responses of the epoxy/concrete interface. The environmental sequences associated with the interfacial bonding degradation are arranged as Na{sub 2}SO{sub 4} + NaCl>NaCl>Na{sub 2}SO{sub 4} > Water>Dry. Remarkably, water molecules weaken the interaction between the epoxy and CSH by robbing interfacial H- and Ca-O bonds, thereby reducing the energy barrier to interfacial peel and shear failure. The weakening effects are exacerbated by the presence of aggressive ions, as Na{sup +} ions attracted to the CSH surface accumulate more water molecules in the interfacial region by forming hydrated ions clusters, which ultimately accelerates the epoxy/CSH interfacial degradation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy

The ease with which the pH is routinely determined for aqueous solutions masks the fact that the cationic product of Arrhenius acid dissolution, the hydrated proton, or H + (aq), is a remarkably complex species. Here, we review how results obtained over the past 30 years in the study of H + ·(H 2 O) n cluster ions isolated in the gas phase shed light on the chemical nature of H + (aq). This effort has also revealed molecular-level aspects of the Grotthuss relay mechanism for positive-charge translocation in water. Recently developed methods involving cryogenic cooling in radiofrequency ion traps and the application of two-color, infrared–infrared (IR–IR) double-resonance spectroscopy have established a clear picture of how local hydrogen-bond topology drives the diverse spectral signatures of the excess proton. This information now enables a new generation of cluster studies designed to unravel the microscopic mechanics underlying the ultrafast relaxation dynamics displayed by H + (aq).

Chemistry↗
DOE OSTI2023

Correlations between Molecular Structure, Solvation Topology, and Transport Properties of Aqueous Organic Flow Battery Electrolyte Solutions

Aqueous organic redox flow batteries (AORFBs) are considered promising technologies for storing energy generated from renewable resources. However, designing organic electrolyte molecules is limited by gaps between fundamental understanding of coupling between solvation structure and dynamics, and macroscopic transport properties like viscosity. Herein, we used molecular dynamics simulations to understand correlations between ionic molecular structures, ion clustering, and transport properties in 2,3-dihydrophenazine (2,3-DHP), a promising AORFB anolyte. We show that experimentally measured viscosity can be reproduced from simulations at relevant concentrations and that the asymmetric structure of 2,3-DHP leads to unique inhomogeneity in the solvation topology. However, order parameters and metrics need to be developed for better correlations over spatiotemporal scales, with careful consideration of the inhomogeneity of organic anolyte molecules. In conclusion, we show that the increased size and asymmetry of the anolyte leads to breakdown of assumptions within methods for determining ion transport mechanisms previously developed for Li-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Modifying Li + and Anion Diffusivities in Polyacetal Electrolytes: A Pulsed-Field-Gradient NMR Study of Ion Self-Diffusion

Polyacetal electrolytes have been demonstrated as promising alternatives to liquid electrolytes and poly(ethylene oxide) (PEO) for rechargeable lithium-ion batteries; however, the relationship between polymer structure and ion motion is difficult to characterize. Here, we study structure-property trends in ion diffusion with respect to polymer composition for a systematic series of five polyacetals with varying ratios of ethylene oxide (EO) to methylene oxide (MO) units, denoted as P(xEO-yMO), and PEO. We first use 7 Li and 19 F pulsed-field-gradient NMR spectroscopy to measure cation and anion self-diffusion, respectively, in polymer/lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt mixtures. At 90 °C, we observe modest changes in Li + diffusivity across all polymer compositions, while anion (TFSI - ) self-diffusion coefficients decrease significantly with increasing MO content. At a given reduced temperature (T - T g ), all polyacetal electrolytes exhibit faster Li + self-diffusion than PEO. Intriguingly, P(EO-MO) and P(EO-2MO) also show slower TFSI - anion self-diffusion than PEO at a given reduced temperature. Molecular dynamics simulations reveal that shorter distances between acetal oxygen atoms (O-CH 2 -O) compared to ether oxygens (O-CH 2 -CH 2 -O) promote more diverse, often asymmetric, Li + coordination environments. Finally, Raman spectra reveal that anion-rich ion clusters in P(EO-MO) and P(EO-2MO) lead to decreased anion diffusivity, which along with increased cation diffusivity, support the viability of polyacetals as high-performance polymer electrolytes.

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

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

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

From classical to quantum dynamics of atomic and ionic species interacting with graphene and its analogue

Graphene and its analogues offer a broad range of application opportunities for (opto)-electronics, sensing, catalysis, phase separation, energy conversion and storage, etc. Engineering graphene properties often relies on its controllable functionalization, defect formation and patterning, and reactive gas etching. In this chapter, we survey the dynamics of graphene using classical and quantum-classical dynamics methods. We discuss the reactivity, scattering, and transmission of atomic and ionic species including Ar cluster ion, H/D, and H+/D+ on graphene flakes of various sizes, focusing on the atomic-scale motion and energy dissipation pathways involved in forming and breaking covalent bonding. Discussions on the nuclear quantum effects of light species, the effects of isotopic substitution, and the methodologies for such modeling are also included.

Garashchuk, Sofiya V.↗
DOE OSTI2020

Consequences of product inhibition in the quantification of kinetic parameters

While the potential for product inhibition in catalytic reactions is well known, the impact of neglected inhibition on measured kinetic parameters is often overlooked. The presence of product inhibition, most often caused by the competitive adsorption of products with reactants on catalytic active sites, is difficult to determine a priori for an arbitrary catalytic system. The significance of product inhibition relies on the concentration of the products, their adsorption thermodynamics on catalytically relevant sites, and process parameters such as temperature and pressure. When inhibition is significant, however, apparent activation energies and reaction orders vary from the differential-reactor apparent activation energy by a factor of 1/(1 - δ), where δ is the total inhibition order (e.g., the factor (1 - δ) = 1.6 for a system with product inhibition of -0.6 order). This is illustrated here with the kinetics of NO oxidation over Cu ion clusters (Cu x O y ) in Cu-SSZ-13, for which the product NO 2 inhibits the forward reaction. Furthermore, in the presence of inhibition, when only reactants are fed to a flow reactor or placed in a batch reactor, there is often no practical conversion that is low enough to guarantee differential behavior. Inclusion of products in the feed solves this problem, allowing accurate determination of kinetic parameters such as apparent reaction orders and activation energies. We also demonstrate that evaluation of the necessity of co-feeding products to assure measurement of differential-reactor data in a given catalytic system is straightforward from a plot of the log of the rate (or conversion) versus the log of the space time in a flow reactor or elapsed time in a batch reactor. We encourage inclusion of this test in all kinetic analyses that are reasonably approximated by power law rate expressions.

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

Correlating Solvation Free Energy to Electrolyte Properties for Lithium Metal Batteries

The electrolyte plays a critical role in lithium metal batteries. In particular, ion solvation profoundly impacts key electrolyte properties and battery performance. Here, in this study, we systematically investigate solvation-property relationships in a series of electrolytes with different solvent-diluent ratios. Through potentiometric techniques that measure the relative solvation free energies of electrolytes, we find that weaker solvation correlates with larger ion clusters, lower ionic conductivity and diffusion coefficient, and superior electrochemical stability. Weaker solvation leads to the formation of a small number of Li clusters with large hydrodynamic radii, which lowers the Li + diffusivity and ionic conductivity of the electrolyte. Concurrently, weaker solvation leads to improved electrochemical stability at both the cathode and anode interfaces. Understanding these solvation-property relationships and trade-offs is important to designing electrolytes for optimized lithium metal battery performance.

25 ENERGY STORAGE↗
DOE OSTI2023

High-entropy electrolytes for practical lithium metal batteries

Electrolyte engineering is crucial for improving battery performance, particularly for lithium metal batteries. Recent advances in electrolytes have greatly improved cyclability by enhancing electrochemical stability at the electrode interfaces, but concurrently achieving high ionic conductivity has remained challenging. Here we report an electrolyte design strategy for enhanced lithium metal batteries by increasing the molecular diversity in electrolytes, which essentially leads to high entropy electrolytes (HEEs). We find that in weakly solvating electrolytes, the entropy effect reduces ion clustering while preserving the characteristic anion-rich solvation structures, which is characterized by synchrotron-based X-ray scattering and molecular dynamics simulations. Electrolytes with smaller- sized clusters exhibit a 2-fold improvement in ionic conductivity compared to conventional weakly- solvating electrolytes, enabling stable cycling at high current densities up to 2C (6.2 mA cm -2 ) in anode- free LiNi 0.6 Mn 0.2 Co 0.2 (NMC622)||Cu pouch cells. In conclusion, the efficacy of the design strategy is verified by performance improvements in three disparate weakly solvating electrolyte systems.

25 ENERGY STORAGE↗
DOE OSTI2021

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

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

Formation and optical properties of IR photoluminescence centres in bismuth-containing aluminophosphate glass

We have prepared aluminophosphate glass samples differing in bismuth oxide content and demonstrating broadband near-IR photoluminescence. Analysis of their photoluminescence spectra leads us to conclude that they contain two main types of emission centres, one of which seems to be a bismuth monocation and the other is a cluster ion. In addition to these luminescence centres, the aluminophosphate glasses contain a nonluminescent bismuth-containing centre responsible for the broad, strong optical absorption band peaking at 450 nm. (active media)

36 MATERIALS SCIENCE↗