Search NASA⌕ Search

SEARCH · Search NASA

Results for “activity descriptors”

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 181 records · Page 10

In-situ IR Spectroscopy Study of Reactions of C3 Oxygenates on Heteroatom (Sn, Mo, and W) doped BEA Zeolites and the Effect of Co-adsorbed Water

The reactions of acetone and hydroxyacetone over heteroatom doped BEA zeolites (Sn, Mo, and W) in the presence and absence of H 2 O vapor are investigated using infrared spectroscopy. Acetone is converted to mesityl oxide over Sn-BEA exclusively. At higher temperatures, larger oxygenates such as phorones, aromatics, and coke form. The presence of co-adsorbed water in Sn-BEA suppresses tautomerization. H 2 O vapor is also beneficial for minimizing coke formation at high temperatures. Hydroxyacetone is converted into 2-hydroxypropanal over Sn-BEA, exhibiting high affinity to Sn sites up to 400 °C. Sn-BEA catalyzes conversion of hydroxyacetone into the enol in the absence of H 2 O, but exposure to H 2 O induces the formation of 2-hydroxypropanal and subsequent conversion to acrolein. Furthermore, the Lewis acid descriptors are used to rationalize the reaction pathways. For the isomerization of hydroxyacetone into 2-hydroxypropanal, the hardness of acid sites influences the reaction and correlates with the overall Lewis acidity of the catalysts, respectively. However, the size of the exchanged metal significantly affects aldol condensation, where keto and enol forms of acetone adsorb to active sites simultaneously.

36 MATERIALS SCIENCE↗

An active learning high-throughput microstructure calibration framework for solving inverse structure–process problems in materials informatics

Determining a process–structure–property relationship is the holy grail of materials science, where both computational prediction in the forward direction and materials design in the inverse direction are essential. Problems in materials design are often considered in the context of process–property linkage by bypassing the materials structure, or in the context of structure–property linkage as in microstructure-sensitive design problems. However, there is a lack of research effort in studying materials design problems in the context of process–structure linkage, which has a great implication in reverse engineering. In this paper, given a target microstructure, we propose an active learning high-throughput microstructure calibration framework to derive a set of processing parameters, which can produce an optimal microstructure that is statistically equivalent to the target microstructure. The proposed framework is formulated as a noisy multi-objective optimization problem, where each objective function measures a deterministic or statistical difference of the same microstructure descriptor between a candidate microstructure and a target microstructure. Furthermore, to significantly reduce the physical waiting wall-time, we enable the high-throughput feature of the microstructure calibration framework by adopting an asynchronously parallel Bayesian optimization to exploit high-performance computing resources. Case studies in additive manufacturing and grain growth are used to demonstrate the applicability of the proposed framework, where kinetic Monte Carlo (kMC) simulation is used as a forward predictive model, such that for a given target microstructure, the target processing parameters that produced this microstructure are successfully recovered.

36 MATERIALS SCIENCE↗

The electron-transfer intermediates of the oxygen evolution reaction (OER) as polarons by in situ spectroscopy

The conversion of diffusive forms of energy (electrical and light) into short, compact chemical bonds by catalytic reactions regularly involves moving a carrier from an environment that favors delocalization to one that favors localization. While delocalization lowers the energy of the carrier through its kinetic energy, localization creates a polarization around the carrier that traps it in a potential energy minimum. The trapped carrier and its local distortion—termed a polaron in solids—can play a role as a highly reactive intermediate within energy-storing catalytic reactions but is rarely discussed as such. Here, we present this perspective of the polaron as a catalytic intermediate through recent in situ and time-resolved spectroscopic investigations of photo-triggered electrochemical reactions at material surfaces. The focus is on hole-trapping at metal–oxygen bonds, denoted M–OH*, in the context of the oxygen evolution reaction (OER) from water. The potential energy surface for the hole-polaron defines the structural distortions from the periodic lattice and the resulting “active” site of catalysis. This perspective will highlight how current and future time-resolved, multi-modal probes can use spectroscopic signatures of M–OH* polarons to obtain kinetic and structural information on the individual reaction steps of OER. A particular motivation is to provide the background needed for eventually relating this information to relevant catalytic descriptors by free energies. Lastly, the formation of the O–O chemical bond from the consumption of M–OH*, required to release O 2 and store energy in H 2 , will be discussed as the next target for experimental investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Operando leaching of pre-incorporated Al and mechanism in transition-metal hybrids on carbon substrates for enhanced charge storage

Insufficient exposure and utilization of active sites often induces an inferior reactivity for transition-metal-based two-dimensional (2D) materials. In response, we for the first time propose a universal "nano-tailoring" strategy to incorporate abundant defects and active sites into low-crystallinity nanosheets by electrochemically leaching of Al species. With MnAl layered double hydroxides (LDHs) as a representative example, potassium-birnessite MnO 2 (AK-MnO 2 ) with oxygen vacancies and abundant edge sites is successfully produced. The oxygen vacancies are shown to help optimize the electron-transfer and ion-adsorption capability. These integrated advantages endow the AK-MnO 2 with a high capacitance value of 239 F g –1 at 100 A g –1 . By further combining with soft X-ray absorption spectroscopy techniques, we unravel that the reducibility of M 2+ in M 2+ Al-LDH serves as the key descriptor for the reconstruction rate. Here, this "nano-tailoring" strategy can provide some important implications and clues to manipulating 2D materials for efficient energy storage and conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deep learning of interface structures from simulated 4D STEM data: cation intermixing vs. roughening ∗

Abstract Interface structures in complex oxides remain an active area of condensed matter physics research, largely enabled by recent advances in scanning transmission electron microscopy (STEM). Yet the nature of the STEM contrast in which the structure is projected along the given direction precludes separation of possible structural models. Here, we utilize deep convolutional neural networks (DCNN) trained on simulated 4D STEM datasets to predict structural descriptors of interfaces. We focus on the widely studied interface between LaAlO 3 and SrTiO 3 , using dynamical diffraction theory and leveraging high performance computing to simulate thousands of possible 4D STEM datasets to train the DCNN to learn properties of the underlying structures on which the simulations are based. We test the DCNN on simulated data and show that it is possible (with >95% accuracy) to identify a physically rough from a chemically diffuse interface and create a DCNN regression model to predict step positions. We quantify the applicability of the model to different thicknesses and the transferability of the approach. The method shown here is general and can be applied for any inverse imaging problem where forward models are present.

42 ENGINEERING↗

Catalytic Hydrogen Production from Methane: A Review on Recent Progress and Prospect

Natural gas (Methane) is currently the primary source of catalytic hydrogen production, accounting for three quarters of the annual global dedicated hydrogen production (about 70 M tons). Steam–methane reforming (SMR) is the currently used industrial process for hydrogen production. However, the SMR process suffers with insufficient catalytic activity, low long-term stability, and excessive energy input, mostly due to the handling of large amount of CO 2 coproduced. With the demand for anticipated hydrogen production to reach 122.5 M tons in 2024, novel and upgraded catalytic processes are desired for more effective utilization of precious natural resources. In this review, we summarized the major descriptors of catalyst and reaction engineering of the SMR process and compared the SMR process with its derivative technologies, such as dry reforming with CO 2 (DRM), partial oxidation with O 2 , autothermal reforming with H 2 O and O 2 . Finally, we discussed the new progresses of methane conversion: direct decomposition to hydrogen and solid carbon and selective oxidation in mild conditions to hydrogen containing liquid organics (i.e., methanol, formic acid, and acetic acid), which serve as alternative hydrogen carriers. We hope this review will help to achieve a whole picture of catalytic hydrogen production from methane.

08 HYDROGEN↗

Covalent Functionalization of Nickel Phosphide Nanocrystals with Aryl-Diazonium Salts

Covalent functionalization of Ni 2 P nanocrystals was demonstrated using aryl-diazonium salts. Spontaneous adsorption of aryl functional groups was observed, with surface coverages ranging from 20 to 96% depending on the native reactivity of the salt as determined by the aryl substitution pattern. Increased coverage was possible for low reactivity species using a sacrificial reductant. Functionalization was confirmed using thermogravimetric analysis, Fourier transform infrared spectroscopy, and X-ray photoelectron spectroscopy. The structure and energetics of this nanocrystal electrocatalyst system, as a function of ligand coverage, were explored with density functional theory calculations. The Hammett parameter of the surface functional group was found to linearly correlate with the change in Ni and P core-electron binding energies and the nanocrystal's experimentally and computationally determined work function. The electrocatalytic activity and stability of the functionalized nanocrystals for hydrogen evolution were also improved when compared to the unfunctionalized material, but a simple trend based on electrostatics was not evident. Density functional theory was used to understand this discrepancy, revealing that H adsorption energies on the covalently functionalized Ni 2 P also do not follow the electrostatic trend and are predictive descriptors of the experimental results.

14 SOLAR ENERGY↗

Cluster-Graph Fingerprinting: A Framework for Quantitative Analysis of Machine-Learned Interatomic Model Training and Simulation Data

Machine-learned interatomic models represent a significant advancement in simulation methods, extending the predictive ability of first-principles methods to previously inaccessible length and time scales. However, the data-driven nature of these models can lead to difficult-to-detect errors that can compromise prediction accuracy. To address this challenge, we introduce a novel fingerprinting approach based on the Chebyshev Interaction Model for Efficient Simulation (ChIMES) ML-IAM graph-based descriptor. Our strategy enables efficient and statistically rigorous analysis of system configurations used in ML-IAM training and those generated by their application, e.g., in molecular dynamics simulations. We demonstrate that these fingerprints can effectively assess novelty of a configuration relative to an existing data set and determine dissimilarity among individual configurations, which are two key tasks in workflows for active learning-based ML-IAM training, data set curation, and on-the-fly uncertainty quantification.

36 MATERIALS SCIENCE↗

Automated AI-driven Molecular Design for Therapeutic Discovery

In recent years, artificial intelligence and machine learning (AI/ML) approaches have revolutionized the process of designing new therapeutics, enabling scientists to rapidly respond to emerging threats from various pathogens. A prime example is the SARS-CoV-2 main protease, a key target for the development of antiviral inhibitors. In this study, we employed a novel, integrated approach that combines AI-driven iterative design of inhibitor candidates, screening based on physio-chemical properties and toxicity, physics-based computational modeling of protein-inhibitor interactions, and AI-assisted analysis of Native MS biophysical assay and characterization of designed candidates. Our deep learning 3D-scaffold model, which uses an input scaffold as a starting point, generated tens of thousands of compounds while preserving the key scaffold. To optimize these candidates, we calculated a comprehensive set of 136 descriptors, including both 2D and 3D molecular features, for compounds targeting the SARS-CoV-2 Main protease (Mpro) and a neurodegenerative disease-associated protein, cyclophilin (Cyp). The generated compounds were initially filtered based on their properties and then ranked according to their predicted binding affinity using our automated modeling and ML methods. Experimental validation of the Mpro candidates showing inhibitory activity demonstrates that our workflow can expedite the therapeutic discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-destructive electrochemical diagnosis of failure mechanisms in aqueous zinc batteries

The early detection of secondary reactions that affect the life and performance of zinc manganese oxide batteries requires a shift from conventional time-consuming and often destructive procedures to rapid lifetime-predictive techniques. In this work, an electrochemical approach is employed to elucidate independent signatures for four common types of failure mechanisms in zinc manganese dioxide (Zn||MnO2) batteries—namely, the loss of zinc inventory, the loss of active material at the cathode, electrolyte depletion, and increased cell impedance. Our findings, specific to coin cell configurations, reveal that each induced failure mechanism can be distinctively modeled and identified based on responses from the rest voltage and columbic-efficiency data for prompt detection. For instance, electrolyte depletion response manifests a distinctive abrupt (>80 %) decrease in columbic efficiency (CE) and charge-rest voltage (Vc) while the discharge-rest voltage remained constant at ~1.3 V. Furthermore, electrolyte rejuvenation of the cell increased the CE to >95 % and restored Vc from ~0.3 to >1.7 V. Recovery experiments and reference performance tests demonstrated consistency between electrochemical descriptors and their associated failure mechanisms. Further, the outcomes of this work provide valuable insights and data models for some of the dominant failure mechanisms present in zinc manganese battery chemistries, which are beneficial to accelerated early-lifetime diagnosis and advancement of Zn batteries development.

25 ENERGY STORAGE↗

Understanding and design of interstitial oxygen conductors

Highly efficient oxygen-active materials that react with, absorb, and transport oxygen is essential for fuel cells, electrolyzers and related applications. While vacancy-mediated oxygen-ion conductors have long been the focus of research, they are limited by high migration barriers at intermediate temperatures (400–600 °C), which hinder their practical applications. In contrast, interstitial oxygen conductors exhibit significantly lower migration barriers enabling higher ionic conductivity at lower temperatures. This review systematically examines both well-established and recently identified families of interstitial oxygen-ion conductors, focusing on how their unique structural motifs such as corner-sharing polyhedral frameworks, isolated polyhedral, and cage-like architectures, facilitate low migration barriers through interstitial and/or interstitialcy diffusion mechanisms. A central discussion of this review focuses on the evolution of design strategies, from targeted donor doping, element screening, to physical-intuition descriptor material screening and machine learning approach, which leverage computational tools to explore vast chemical spaces in search for new interstitial conductors. The success of these strategies demonstrates that a significant, largely unexplored space remains for discovering high-performing interstitial oxygen conductors. Crucial features enabling high-performance interstitial oxygen diffusion include the availability of electrons for oxygen reduction and sufficient structural flexibility with accessible volume for interstitial accommodation and migration. This review concludes with a forward-looking perspective, proposing a knowledge-driven methodology that integrates current understanding with data-centric approaches to identify promising interstitial oxygen conductors outside traditional search paradigms. These approaches are expected to significantly accelerate the development of high-performance interstitial oxygen conductors for a variety of oxygen-active applications, ultimately paving the way for more efficient and sustainable energy technologies.

Interstitial oxygen conductors↗

Decoupling plasma, catalyst, and gaseous mechanisms for non-oxidative methane conversion

Direct non-oxidative methane (CH 4 ) conversion to value-added hydrogen (H 2 ) and C 2 products remains hindered by fundamental catalytic scaling constraints and rapid surface deactivation at elevated temperatures. Plasma-enabled catalysis offers a promising route to overcome the thermodynamic and kinetic barriers of direct non-oxidative methane upgrading at mild conditions, yet control over C–C product selectivity and catalyst stability remains elusive. Here, we establish a unified mechanistic framework including Langmuir–Hinshelwood (L–H) and Langmuir–Rideal (L–R) mechanisms that disentangles the roles of plasma excitation (including vibrationally activated species and radicals), surface temperature (T sur ), and catalyst binding energy in steering CH 4 conversion to H 2 and C 2 hydrocarbons. Through a combination of density functional theory (DFT) informed microkinetic modeling, in situ and ex situ surface characterization, and product quantification under dielectric barrier discharge conditions, we show that vibrationally excited CH 4 lowers activation barriers selectively for dissociative chemisorption, enabling surface activation across a wide range of transition metal catalysts at low thermal energy input. We find that once CH 4 is dissociatively chemisorbed, the branching between C 2 H 2 , C 2 H 4 , and C 2 H 6 is governed by surface properties (carbon binding energy, T sur , etc), regardless of plasma excitation. The DFT informed microkinetic model decouples the effects of molecular activation from surface properties and indentifies operating windows that maximize target yields while suppressing carbon accumulation and subsequent catalytic inactivation. Experiments on polycrystalline Cu/Al 2 O 3 , Ni/Al 2 O 3 , and Pt/Al 2 O 3 validate these predictions, revealing catalyst-dependent branching toward ethane or ethylene and distinct deactivation profiles. We unify these trends into a generalized three-dimensional plasma-thermal-catalytic design space, from which reduced descriptors such as T vib /T sur in the limit of vibrationally excited L–H pathways emerge as predictive metrics. These results enable rational tuning of methane conversion pathways and unlock selective C 2 formation using earth-abundant metals under mild plasma conditions.

catalyst inactivation↗

Non-bonding interaction of dual atom catalysts for enhanced oxygen reduction reaction

Here, we demonstrate the design of graphene-supported dual atom catalysts (DACs) for the four-electron oxygen reduction reaction (ORR), by utilizing the non-bonding interaction of counterpart metals (M) that synergistically tune the electronic properties and catalytic activity of the Fe active site in FeMN 6 -DAC and FeMN 8 -DAC systems, where M stands for Fe, Co, Ni, Cu, and Zn. More specifically, for Fe-M distances below 15Å, the non-bonding interaction is significant, making the system act as the DAC. We predicted that FeNiN 6 -DAC and FeNiN 8 -DAC exhibit a low ORR overpotential (η ORR ) of 0.28 V and 0.47 V, respectively, which are at the summits of volcano plots. This low η ORR originates from the high Bader charge transfer coupled with high spin density at the Fe site in both the FeNiN 6 -DAC and FeNiN 8 -DAC systems, which weakens the adsorption of OH* intermediate while enhancing its desorption to H 2 O. Guided by these density functional theory (DFT) computational results, we synthesized FeCoN 8 -DAC and FeNiN 8 -DAC along with N-doped graphene and confirmed their structures with scanning transmission electron microscopy (STEM), X-ray photoelectron spectroscopy (XPS), X-ray absorption near-edge structure (XANES), extended X-ray absorption fine structure (EXAFS), and electron spin resonance (ESR). We verify experimentally the catalytic activities and find that FeNiN 8 -DAC has the low experimental overpotential of 0.39 V with a Tafel slope of 47 mVdec -1 . Based on these results, we propose a DFT-guided strategy to tune the charge transfer and spin population of the active site toward designing DACs for electrochemical ORR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Remote sensing: Earth's surface and atmosphere; Proceedings of Workshop X and the Topical Meeting of the 26th COSPAR Plenary Meeting, Toulouse, France, June 30-July 11, 1986

The present conference on space-based remote sensing of the earth's surface and atmosphere addresses the two broad issues of remote sensing activities of interest to developing countries and the results obtained to date by the International Satellite Cloud Climatology Project, the Earth Radiation Budget Experiment, and the Stratospheric Aerosol and Gas Experiment (SAGE). Attention is given to the remote sensing of environmental factors affecting health, applications of satellite microwave radiometry, earth science missions for the NASA Space Station, and digitally produced Landsat map images. Also discussed are time-accumulated visible and IR histograms used as cloud cover descriptors, the estimation of the radiation budget's sensitivity to cloud variations, monitoring global surface temperature variations using cloud data sets, and an analysis of preliminary SAGE II data on ozone and NO2.

Carter, W. D.↗

Predictions of Leukemia Risks to Astronauts from Solar Particle Events

Leukemias consisting of acute and chronic myeloid leukemia and acute lymphatic lymphomas represent the earliest cancers that appear after radiation exposure, have a high lethality fraction, and make up a significant fraction of the overall fatal cancer risk from radiation for adults. Several considerations impact the recommendation of a preferred model for the estimation of leukemia risks from solar particle events (SPE's): The BEIR VII report recommends several changes to the method of calculation of leukemia risk compared to the methods recommended by the NCRP Report No. 132 including the preference of a mixture model with additive and multiplicative components in BEIR VII compared to the additive transfer model recommended by NCRP Report No. 132. Proton fluences and doses vary considerably across marrow regions because of the characteristic spectra of primary solar protons making the use of an average dose suspect. Previous estimates of bone marrow doses from SPE's have used an average body-shielding distribution for marrow based on the computerized anatomical man model (CAM). We have developed an 82-point body-shielding distribution that faithfully reproduces the mean and variance of SPE doses in the active marrow regions (head and neck, chest, abdomen, pelvis and thighs) allowing for more accurate estimation of linear- and quadratic-dose components of the marrow response. SPE's have differential dose-rates and a pseudo-quadratic dose response term is possible in the peak-flux period of an event. Also, the mechanistic basis for leukemia risk continues to improve allowing for improved strategies in choosing dose-rate modulation factors and radiation quality descriptors. We make comparisons of the various choices of the components in leukemia risk estimates in formulating our preferred model. A major finding is that leukemia could be the dominant risk to astronauts for a major solar particle event.

Cucinotta, F. A.↗

Data Science-Driven Analysis of Substrate-Permissive Diketopiperazine Reverse Prenyltransferase NotF: Applications in Protein Engineering and Cascade Biocatalytic Synthesis of (-)-Eurotiumin A

Prenyltransfer is an early-stage carbon-hydrogen bond (C-H) functionalization prevalent in the biosynthesis of a diverse array of biologically active bacterial, fungal, plant, and metazoan diketopiperazine (DKP) alkaloids. Toward the development of a unified strategy for biocatalytic construction of prenylated DKP indole alkaloids, we sought to identify and characterize a substrate-permissive C2 reverse prenyltransferase (PT). As the first tailoring event within the biosynthesis of cytotoxic notoamide metabolites, PT NotF catalyzes C2 reverse prenyltransfer of brevianamide F. Solving a crystal structure of NotF (in complex with native substrate and prenyl donor mimic dimethylallyl S-thiolodiphosphate (DMSPP)) revealed a large, solvent-exposed active site, intimating NotF may possess a significantly broad substrate scope. To assess the substrate selectivity of NotF, we synthesized a panel of 30 sterically and electronically differentiated tryptophanyl DKPs, the majority of which were selectively prenylated by NotF in synthetically useful conversions (2 to > 99%). Quantitative representation of this substrate library and development of a descriptive statistical model provided insight into the molecular origins of NotF's substrate promiscuity. This approach enabled the identification of key substrate descriptors (electrophilicity, size, and flexibility) that govern the rate of NotF-catalyzed prenyltransfer, and the development of an "induced fit docking (IFD)-guided" engineering strategy for improved turnover of our largest substrates. We further demonstrated the utility of NotF in tandem with oxidative cyclization using flavin monooxygenase, BvnB. This one-pot, in vitro biocatalytic cascade enabled the first chemoenzymatic synthesis of the marine fungal natural product, (-)-eurotiumin A, in three steps and 60% overall yield.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Yield Room-Temperature Aryl 211 At Radiolabeling via Redox-Activated Column-Shipped 211 At

Efficient late-stage formation of aryl– 211 At bonds is relevant to expanding the scope of targeted α-therapy radiopharmaceuticals. In this work, a catalyst-enabled astatodeboronation approach is evaluated using wet chemically regenerated 211 At recovered from shipped 3-octanone-impregnated columns. A series of aryl boronic esters bearing electron-donating and electron-withdrawing substituents were radiolabeled under mild conditions, with reaction parameters examined through time-dependent studies. Radiochemical conversions ranging from 68–97% were observed across the substrate set. Density functional theory calculations were employed as a qualitative framework to contextualize relative B–C and At–C bonding descriptors. Overall, this study demonstrates that wet chemically regenerated 211 At can be applied to the late-stage labeling of aryl boronic esters under mild conditions and provides guidance for future exploration of 211 At radiolabeling strategies.

211At-aryl labeling↗

A Computational Framework to Accelerate the Discovery of Perovskites for Solar Thermochemical Hydrogen Production: Identification of Gd Perovskite Oxide Redox Mediators

A high-throughput computational framework to identify novel multinary perovskite redox mediators is presented, and this framework is applied to discover the Gd-containing perovskite oxide compositions Gd 2 BB'O 6 , GdA'B 2 O 6 , and GdA'BB'O 6 that split water. The computational scheme uses a sequence of empirical approaches to evaluate the stabilities, electronic properties, and oxygen vacancy thermodynamics of these materials, including contributions to the enthalpies and entropies of reduction, ΔH TR and ΔS TR . This scheme uses the machine-learned descriptor τ to identify compositions that are likely stable as perovskites, the bond valence method to estimate the magnitude and phase of BO 6 octahedral tilting and provide accurate initial estimates of perovskite geometries, and density functional theory including magnetic- and defect-sampling to predict STCH-relevant properties. Eighty-three promising STCH candidate perovskite oxides down-selected from 4392 Gd-containing compositions are reported, three of which are referred to experimental collaborators for characterization and exhibit STCH activity. Our results demonstrate that the high-throughput computational scheme described herein—which is used to evaluate Gd-containing compositions but can be applied to any multinary perovskite oxide compositional space(s) of interest—accelerates the discovery of novel STCH active redox mediators with reasonable computational expense.

36 MATERIALS SCIENCE↗