Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical mechanics models”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling noise in global Mølmer-Sørensen interactions applied to quantum approximate optimization

Many-qubit Mølmer-Sørensen (MS) interactions applied to trapped ions offer unique capabilities for quantum information processing, with applications including quantum simulation and the quantum approximate optimization algorithm (QAOA). Here, we develop a physical model to describe many-qubit MS interactions under four sources of experimental noise: vibrational mode frequency fluctuations, laser power fluctuations, thermal initial vibrational states, and state preparation and measurement errors. The model parametrizes these errors from simple experimental measurements, without free parameters. Herein, we validate the model in comparison with experiments that implement sequences of MS interactions on two 171 Yb + ions. The model shows reasonable agreement after several MS interactions as quantified by the reduced chi-squared statistic $χ^{2}_{red}$≈2. As an application we examine MaxCut QAOA experiments on three and six ions. The experimental performance is quantified by approximation ratios that are 91% and 83% of the optimal theoretical values. Our model predicts $0.93^{+0.03}_{-0.02}$ and $0.95^{+0.04}_{-0.03}$, respectively, with disagreement in the latter value attributable to secondary noise sources beyond those considered in our analysis. With realistic experimental improvements to reduce measurement error and radial trap frequency variations, the model achieves approximation ratios that are 99% of the optimal. Incorporating these improvements into future experiments is expected to reveal new aspects of noise for future modeling and experimental improvements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Nonequilibrium universality of the nonreciprocally coupled 𝑂⁡(𝑛 1 ) × 𝑂⁡(𝑛 2 ) model

Nonequilibrium dynamics play an important role in all contexts of physics, both classical and quantum as well as living and nonliving, so it is crucial to develop a foundational understanding of nonequilibrium phase transitions. In this work we investigate an important class of nonequilibrium dynamics in the form of nonreciprocal interactions. In particular we study how nonreciprocal coupling between two 𝑂⁡(𝑛𝑖) order parameters (with 𝑖 = 1,2) affects the universality at a multicritical point, extending the analysis of J. T. Young et al. [Phys. Rev. X 10, 011039 (2020)], which considered the case 𝑛 1 = 𝑛 2 = 1, i.e., a ℤ 2 × ℤ 2 model. We show that nonequilibrium fixed points (NEFPs) emerge for a broad range of 𝑛 1 ,𝑛 2 and exhibit intrinsically nonequilibrium critical phenomena, namely a violation of fluctuation-dissipation relations at all scales and underdamped oscillations near criticality in contrast to the overdamped relaxational dynamics of the corresponding equilibrium models. Furthermore, the NEFPs exhibit an emergent discrete scale invariance in certain physically relevant regimes of 𝑛 1 ,𝑛 2 , but not others, depending on whether the critical exponent 𝜈 is real or complex. The boundary between these two regions is described by an exceptional point in the renormalization group (RG) flow, leading to distinctive features in correlation functions and the phase diagram. Another contrast with the previous work is the number and stability of the NEFPs as well as the underlying topology of the RG flow. Lastly, we investigate an extreme form of nonreciprocity where one order parameter is independent of the other order parameter but not vice versa. Unlike the ℤ 2 × ℤ 2 model, which becomes nonperturbative in this case, we identify a distinct nonequilibrium universality class whose dependent field similarly violates fluctuation-dissipation relations but does not exhibit discrete scale invariance or underdamped oscillations near criticality.

Critical phenomena↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2020

Comparative differential cuproproteomes of Rhodobacter capsulatus reveal novel copper homeostasis related proteins

Abstract Copper (Cu) is an essential, but toxic, micronutrient for living organisms and cells have developed sophisticated response mechanisms towards both the lack and the excess of Cu in their environments. In this study, we achieved a global view of Cu-responsive changes in the prokaryotic model organism Rhodobacter capsulatus using label-free quantitative differential proteomics. Semi-aerobically grown cells under heterotrophic conditions in minimal medium (∼0.3 μM Cu) were compared with cells supplemented with either 5 μM Cu or with 5 mM of the Cu-chelator bathocuproine sulfonate. Mass spectrometry based bottom-up proteomics of unfractionated cell lysates identified 2430 of the 3632 putative proteins encoded by the genome, producing a robust proteome dataset for R. capsulatus. Use of biological and technical replicates for each growth condition yielded high reproducibility and reliable quantification for 1926 of the identified proteins. Comparison of cells grown under Cu-excess or Cu-depleted conditions to those grown under minimal Cu-sufficient conditions revealed that 75 proteins exhibited statistically significant (p < 0.05) abundance changes, ranging from 2- to 300-fold. A subset of the highly Cu-responsive proteins was orthogonally probed using molecular genetics, validating that several of them were indeed involved in cellular Cu homeostasis.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1990

The correction of infrasound signals for upper atmospheric winds

Infrasound waves propagate in the atmosphere by a well known mechanism produced by refraction of the waves, return to earth, and reflection at the surface into the atmosphere for subsequent bounces. A figure illustrates this phenomenon with results from a ray trace model. In this instance three rays are returned to earth from a region centered at about 50 kilometers in altitude and two from a region near 110 kilometers in altitude. The control of the wave refraction is largely dominated by the temperature-height profile and inversions; however, a major influence is also produced by the atmospheric wind profile. Another figure illustrates the considerable ray differences for rays moving in the wind direction (to the right) and in the counter direction (to the left). It obviously can be expected that infrasonic signal amplitudes will be greatly influenced by the winds in the atmosphere. The seasonal variation of the high altitude atmospheric winds is well documented. A third figure illustrates this with average statistics on the observed zonal wind in the region of 50 plus or minus 5 kilometers in altitude. The results are based upon a survey by Webb; Webb terms this parameterization the Stratospheric Circulation Index (SCI). The very strong seasonal variation has the ability to exert a major seasonal influence on infrasonic signals. The purpose here is to obtain a method for the correction of this effect.

Mutschlecner, J. Paul↗
DOE OSTI2025

Widening of Wind Stress Anomalies Amplifies ENSO in a Warming Climate

Abstract Climate change simulations generally indicate the strengthening of El Niño–Southern Oscillation (ENSO) sea surface temperature (SST) variability through the twenty-first century, yet a robust physical mechanism explaining this change across different models is still lacking, and the projections for ENSO amplitude exhibit a large spread. Most commonly, changes in the background state of the tropical Pacific are invoked to explain these changes of ENSO. Here, we show that changes in the structure of wind stress anomalies associated with ENSO are potentially as important as these background state changes. Specifically, changes in the magnitude, meridional width, and zonal structure of wind stress anomalies can explain approximately 53% of the intermodel variance in the projected change of ENSO magnitude through the twenty-first century as well as 43% in ENSO periodicity changes. Among these changes in the wind structure, the meridional widening of wind anomalies plays the most important role. To demonstrate that these changes are indeed critical, we develop a hybrid model of ENSO based on the Community Earth System Model, version 2, which incorporates a dynamical ocean coupled to a simplified statistical atmosphere within the tropical Pacific. In the absence of external forcing and corresponding mean-state changes, the imposed changes in wind stress anomalies in this hybrid model result in an increase of ENSO amplitudes of nearly 10% along the equator. Our results are also theoretically supported by a recharge-oscillator model that incorporates the meridional wind structure. Thus, changes in the structure of wind stress anomalies, together with changes in the mean state, likely play a critical role in the projected strengthening of ENSO.

Stuivenvolt-Allen, Jacob [Yale University, New Hav↗
NASA NTRS2017

Network Analysis of Rodent Transcriptomes in Spaceflight

Network analysis methods leverage prior knowledge of cellular systems and the statistical and conceptual relationships between analyte measurements to determine gene connectivity. Correlation and conditional metrics are used to infer a network topology and provide a systems-level context for cellular responses. Integration across multiple experimental conditions and omics domains can reveal the regulatory mechanisms that underlie gene expression. GeneLab has assembled rich multi-omic (transcriptomics, proteomics, epigenomics, and epitranscriptomics) datasets for multiple murine tissues from the Rodent Research 1 (RR-1) experiment. RR-1 assesses the impact of 37 days of spaceflight on gene expression across a variety of tissue types, such as adrenal glands, quadriceps, gastrocnemius, tibalius anterior, extensor digitorum longus, soleus, eye, and kidney. Network analysis is particularly useful for RR-1 -omics datasets because it reinforces subtle relationships that may be overlooked in isolated analyses and subdues confounding factors. Our objective is to use network analysis to determine potential target nodes for therapeutic intervention and identify similarities with existing disease models. Multiple network algorithms are used for a higher confidence consensus.

genomics↗
DOE OSTI2024

Surface Pressure Fluctuations Induced by a Hypersonic Turbulent Boundary Layer on a Sharp Cone at Angle of Attack

High-fidelity simulations are performed to characterize the turbulence-induced wall pressure fluctuations on a sharp cone at a 5.5° angle-of-attack in a Mach 8 flow. Wall-resolved large-eddy simulation (LES) and wall-modeled large-eddy simulation (WMLES) results are compared to measurements at several locations on the cone body, where comparisons with high-frequency PCB sensors are good, while comparison with low-frequency kulite sensors varies depending on the location. For a given streamwise location, simulation results show significant azimuthal variation in the wall pressure fluctuation statistics. Comparisons between LES and WMLES results indicate that WMLES can accurately reproduce the low-frequency component of the autospectra quite well, while slight sensitivity is apparent in the coherence. A grid sensitivity study for WMLES further demonstrated sensitivity in the wall pressure fluctuation coherence, but not in the autospectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Topology and many-body localization

We discuss the problem of localization in two dimensional electron systems in the quantum Hall (single Landau level) regime. After briefly summarizing the well-studied problem of Anderson localization in the non-interacting case, we concentrate on the problem of disorder induced many-body localization (MBL) in the presence of electron–electron interactions using numerical exact diagonalization and eigenvalue spacing statistics as a function of system size. We provide evidence showing that MBL is not attainable in a single Landau level with short range (white noise) disorder in the thermodynamic limit. We then study the interplay of topology and localization, by contrasting the behavior of topological and nontopological subbands arising from a single Landau level in two models — (i) a pair of extremely flat Hofstadter bands with an optimally chosen periodic potential, and (ii) a Landau level with a split-off nontopological impurity band. Both models provide convincing evidence for the strong effect of topology on the feasibility of many-body localization as well as slow dynamics starting from a nonequilibrium state with charge imbalance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Electrical Responses of Modified Mineral Surfaces as Observed With Spectral Induced Polarization and Atomic Force Microscopy

Abstract Atomic force microscopy (AFM) and spectral induced polarization (SIP) are widely used to investigate the electrical properties of mineral surfaces at vastly different scales of measurement. We compare AFM and SIP measurements made on two different materials (glass beads and silica gel) subjected to etching, deposition of iron oxide particles, and inclusion of calcite grains. We found that the treatments produced qualitatively consistent behaviors in the AFM and SIP data. Direct AFM measurements of surface charge density for silica and calcite surfaces were quantitatively compared to values estimated from the SIP results using a grain polarization model. No statistically significant difference (at a 95% confidence level) was found between the surface charge density of silica estimated by AFM (2.3 ± 6.6 mC/m 2 for glass beads and 1.6 ± 0.1 mC/m 2 for silica gel) versus SIP (5.4 ± 4.4 mC/m 2 for glass beads and 1.6 ± 0.5 mC/m 2 for silica gel). The surface charge density for calcite determined by AFM (43.5 ± 12.9 mC/m 2 ) was approximately 19 times higher than that found for silica. While the charge density of calcite surfaces determined by SIP was also generally higher than that found for silica, different treatments produced significantly different values between 4.7 and 258 mC/m 2 (with a maximum 95% CI of ±8.7 mC/m 2 ). Several possible explanations exist for the range of the observed SIP measurements, including aging of the calcite surfaces. Overall, this study suggests the potential for the complementary use of AFM and SIP measurements to constrain future investigations of polarization mechanisms in porous media.

Geochemistry & Geophysics↗
DOE OSTI2021

Searching for fat tails in CRISPR-Cas systems: Data analysis and mathematical modeling

Understanding CRISPR-Cas systems—the adaptive defence mechanism that about half of bacterial species and most of archaea use to neutralise viral attacks—is important for explaining the biodiversity observed in the microbial world as well as for editing animal and plant genomes effectively. The CRISPR-Cas system learns from previous viral infections and integrates small pieces from phage genomes called spacers into the microbial genome. The resulting library of spacers collected in CRISPR arrays is then compared with the DNA of potential invaders. One of the most intriguing and least well understood questions about CRISPR-Cas systems is the distribution of spacers across the microbial population. Here, using empirical data, we show that the global distribution of spacer numbers in CRISPR arrays across multiple biomes worldwide typically exhibits scale-invariant power law behaviour, and the standard deviation is greater than the sample mean. We develop a mathematical model of spacer loss and acquisition dynamics which fits observed data from almost four thousand metagenomes well. In analogy to the classical ‘rich-get-richer’ mechanism of power law emergence, the rate of spacer acquisition is proportional to the CRISPR array size, which allows a small proportion of CRISPRs within the population to possess a significant number of spacers. Our study provides an alternative explanation for the rarity of all-resistant super microbes in nature and why proliferation of phages can be highly successful despite the effectiveness of CRISPR-Cas systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Uncertainties in the production of iron-group nuclides in core-collapse supernovae from Monte Carlo variations of reaction rates

Core-collapse supernovae, occurring at the end of massive star evolution, produce heavy elements, including those in the iron peak. Although the explosion mechanism is not yet fully understood, theoretical models can reproduce optical observations and observed elemental abundances. However, many nuclear reaction rates involved in explosive nucleosynthesis have large uncertainties, impacting the reliability of abundance predictions. To address this, we have previously developed a Monte Carlo-based nucleosynthesis code that accounts for reaction rate uncertainties and has been applied to nucleosynthesis processes beyond iron. Our framework is also well suited for studying explosive nucleosynthesis in supernovae. In this paper, we investigate 1D explosion models using the ‘PUSH method’ , focusing on progenitors with varying metallicities and initial masses around $M_{\rm ZAMS} = 16\, {\rm M}_{\odot }$. Detailed post-process nucleosynthesis calculations and Monte Carlo analyses are used to explore the effects of reaction rate uncertainties and to identify key reaction rates in explosive nucleosynthesis. We find that many reactions have little impact on the production of iron-group nuclei, as these elements are primarily synthesized in the nuclear statistical equilibrium. However, we identify a few ‘key reactions’ that significantly influence the production of radioactive nuclei, which may affect astrophysical observables. In particular, for the production of ${}^{44}{\rm Ti}$, we confirm that several traditionally studied nuclear reactions have a strong impact. However, determining a single reaction rate is insufficient to draw a definitive conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Dynamic microbial assembly processes correspond to soil fertility in sustainable paddy agroecosystems

1. Microorganisms underpin the sustainability of agroecosystems, but we have a poor understanding of how the assembly processes that govern microbial community structure interact through time to maintain soil fertility. 2. Ancient rice paddies have a history of cultivation spanning hundreds or thousands of years and provide an ideal opportunity to explore relationships between community assembly processes and long-term agroecosystem sustainability. We used statistical null models to evaluate temporal patterns of bacterial community assembly processes across rice cultivation years and their associations with increases in soil fertility in paddy agroecosystems. 3. We found that stochasticity was a persistent force structuring paddy soil microbial communities, putatively due to frequent flooding that facilities dispersal. However, determinism gradually increased over centuries of agriculture and corresponded to increases in soil fertility. Accordingly, multifunctionality, calculated using enzymes catalyzing key soil functions in carbon, nitrogen, phosphorous and sulfur cycling significantly increased along rice cultivation years. We therefore propose that balanced assembly processes preserve a diverse array of traits that can dominate under various circumstances (stochasticity) while allowing for traits relevant to soil fertility to be selected for by soil fertility (determinism), thus supporting multifunctionality that underlies sustainable agroecosystems. 4. In total, our results highlight that knowledge on temporal microbial assembly processes, community composition, and ecological function in context of soil fertility is key to understanding ecological mechanisms maintaining agroecosystems sustainability.

Liu, Wenjing↗
DOE OSTI2021

Molecular dynamics simulation of metallic Al–Ce liquids using a neural network machine learning interatomic potential

Al-rich Al-Ce alloys have the possibility of replacing heavier steel and cast-irons for use in high-temperature applications. Knowledge about the structures and properties of Al-Ce alloys at liquid state is vital for optimizing the manufacture process to produce desired allows. However, reliable molecular dynamics simulation of Al-Ce alloy systems remains a great challenge due to the lack of accurate Al-Ce interatomic potential. In this work, an artificial neural network (ANN) deep machine learning (ML) method is used to develop a reliable interatomic potential for Al-Ce alloy. Ab initio molecular dynamics (AIMD) simulation data on Al-Ce liquid with small unit cell (~200 atoms) and on the known Al-Ce crystalline compounds are collected to train the interatomic potential using ANN-ML. The obtained ANN-ML model reproduces well the energies, forces, and atomic structure of Al 90 Ce 10 liquid and crystalline phases of Al-Ce compounds in comparison with ab initio results. The developed ANN-ML potential is applied in molecular dynamics simulations to study the structures and properties of metallic Al 90 Ce 10 liquid, which would provide useful insight for guiding experimental process to produce desired Al-Ce allows.

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

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Specific features of fluorescence transfer in multiply scattering randomly inhomogeneous layers under intense laser pumping

Based on the analysis of experimental data on the effect of the pulsed laser pump intensity on the spectral properties and the size of the fluorescent response zone in randomly inhomogeneous fluorescent layers, we found that the amplification of spontaneous and stimulated emission significantly affects the statistical properties of the propagation lengths of the fluorescent field partial components in the layers. The experiments are performed with layers of SiO{sub 2} and TiO{sub 2} nanoparticles saturated with rhodamine 6G, pumped by 532-nm laser radiation in the intensity range corresponding to the transient regime from excitation of spontaneous fluorescence to random lasing in the layer. The experimental data are compared with the results of statistical modelling of fluorescence transfer. It is shown that, even at a pump intensity below the random lasing threshold, the spontaneous emission amplification in a layer leads to a significant increase in the contributions to the fluorescence response from partial components with propagation lengths much larger than the layer thickness. This can be interpreted as a manifestation of the quasi-waveguide effect, in which the probability of propagation of diffuse fluorescence components along the layer over distances many times greater than its thickness and the size of the pumped region increases significantly with a decrease in the characteristic radiation amplification length in the layer. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Bayesian-based response expansion and uncertainty quantification using sparse measurement sets

Systems subjected to dynamic loads often require monitoring of their vibrational response, but limitations on the total number and placement of the measurement sensors can hinder the data-collection process. Here, we present an indirect approach to estimate a system’s full-field dynamic response, including all uninstrumented locations, using response measurements from sensors sparsely located on the system. This approach relies on Bayesian inference that utilizes a system model to estimate the full-field response and quantify the uncertainty in these estimates. By casting the estimation problem in the frequency domain, this approach utilizes the modal frequency response functions as a natural, frequency-dependent weighting scheme for the system mode shapes to perform the expansion. This frequency-dependent weighting scheme enables an accurate expansion, even with highly correlated mode shapes that may arise from spatial aliasing due to the limited number of sensors, provided these correlated modes do not have natural frequencies that are closely spaced. Furthermore, the inherent regularization mechanism that arises in this Bayesian-based procedure enables the utilization of the full set of system mode shapes for the expansion, rather than any reduced subset. This approach can produce estimates when considering a single realization of the measured responses, and with some modification, it can also produce estimates for power spectral density matrices measured from many realizations of the responses from statistically stationary random processes. A simply supported beam provides an initial numerical validation, and a cylindrical test article excited by acoustic loads in a reverberation chamber provides experimental validation.

42 ENGINEERING↗
DOE OSTI2021

Learning-based approach to plasticity in athermal sheared amorphous packings: Improving softness

The plasticity of amorphous solids undergoing shear is characterized by quasi-localized rearrangements of particles. While many models of plasticity exist, the precise relationship between the plastic dynamics and the structure of a particle’s local environment remains an open question. Previously, machine learning was used to identify a structural predictor of rearrangements called “softness.” Although softness has been shown to predict which particles will rearrange with high accuracy, the method can be difficult to implement in experiments where data are limited and the combinations of descriptors it identifies are often difficult to interpret physically. Here, we address both of these weaknesses, presenting two major improvements to the standard softness method. First, we present a natural representation of each particle’s observed mobility, allowing for the use of statistical models that are both simpler and provide greater accuracy in limited datasets. Second, we employ persistent homology as a systematic means of identifying simple, topologically informed, structural quantities that are easy to interpret and measure experimentally. We test our methods on two-dimensional athermal packings of soft spheres under quasi-static shear. We find that the same structural information that predicts small variations in the response is also predictive of where plastic events will localize. We also find that an excellent accuracy is achieved in athermal sheared packings using simply a particle’s species and the number of nearest neighbor contacts.

36 MATERIALS SCIENCE↗