Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation model”

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

Machine learning predictions of diffusion in bulk and confined ionic liquids using simple descriptors

Ionic liquids have many intriguing properties and widespread applications such as separations and energy storage. However, ionic liquids are complex fluids and predicting their behavior is difficult, particularly in confined environments. We introduce fast and computationally efficient machine learning (ML) models that can predict diffusion coefficients and ionic conductivity of bulk and nanoconfined ionic liquids over a wide temperature range (350–500 K). The ML models are trained on molecular dynamics simulation data for 29 unique ionic liquids as bulk fluids and confined in graphite slit pores. This model is based on simple physical descriptors of the cations and anions such as molecular weight and surface area. Here, we also demonstrate that accurate results can be obtained using only descriptors derived from SMILES (simplified molecular-input line-entry system) codes for the ions with minimal computational effort. This offers a fast and efficient method for estimating diffusion and conductivity of nanoconfined ionic liquids at various temperatures without the need for expensive molecular dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Deciphering the morphology of transition metal carbonate cathode precursors

The performance and life of Li-ion battery cathode materials are determined by both the composition (crystal structure and transition metal ratio) and the morphology (particle size, size distribution, and surface area). Careful control of these two aspects is the key to long lasting, high-energy batteries that can undergo fast charge. Developing such cathodes requires manipulation of the synthesis conditions, namely the coprecipitation process to develop the precursor and a calcination step to lithiate and convert it to a transition metal oxide. In this paper, we utilize a combination of controlled synthesis, microscopic and spectroscopic characterization, and multi-scale mathematical modeling to shed light on the synthesis of cathode precursors. The complex interplay between the various chemical reactions in the co-precipitation process is studied to provide experimentalists with guidance on achieving composition control during synthesis. Further, the formation of a variety of morphologies of the primary particles and the driving force for agglomeration is mathematically described, for the first time, based on an energy minimization approach. Results suggest that the presence of Ni and/or Co significantly lowers the reaction rate constant compared to Mn, resulting in agglomerated growth in the former and single crystal growth in the latter. Modeling studies are used to provide a phase map describing the synthesis conditions needed to control the secondary particle size and corresponding size distribution. Finally, this paper represents an important step in developing a computationally guided approach to the synthesis of battery cathode materials.

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

Large Ensemble Diagnostic Evaluation of Hydrologic Parameter Uncertainty in the Community Land Model Version 5 (CLM5)

Abstract Land surface models such as the Community Land Model version 5 (CLM5) seek to enhance understanding of terrestrial hydrology and aid in the evaluation of anthropogenic and climate change impacts. However, the effects of parametric uncertainty on CLM5 hydrologic predictions across regions, timescales, and flow regimes have yet to be explored in detail. The common use of the default hydrologic model parameters in CLM5 risks generating streamflow predictions that may lead to incorrect inferences for important dynamics and/or extremes. In this study, we benchmark CLM5 streamflow predictions relative to the commonly employed default hydrologic parameters for 464 headwater basins over the conterminous United States (CONUS). We evaluate baseline CLM5 default parameter performance relative to a large (1,307) Latin Hypercube Sampling‐based diagnostic comparison of streamflow prediction skill using over 20 error measures. We provide a global sensitivity analysis that clarifies the significant spatial variations in parametric controls for CLM5 streamflow predictions across regions, temporal scales, and error metrics of interest. The baseline CLM5 shows relatively moderate to poor streamflow prediction skill in several CONUS regions, especially the arid Southwest and Central U.S. Hydrologic parameter uncertainty strongly affects CLM5 streamflow predictions, but its impacts vary in complex ways across U.S. regions, timescales, and flow regimes. Overall, CLM5's surface runoff and soil water parameters have the largest effects on simulated high flows, while canopy water and evaporation parameters have the most significant effects on the water balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Paleotribological models using preserved fossil tissue properties reveal functional significance of Eurasian mammoth dental evolution

Eurasian mammoths (Mammuthus) underwent substantial modifications in molar morphology as later-diverging species evolved progressively thinner enamel and increased enamel crest complexity. These features have been hypothesized to reduce whole-tooth wear and extend dental longevity as increasingly graze-dominated diets evolved within the lineage. This hypothesis has yet to be directly tested. Here, in this study, we developed an in-silico wear model using experimentally derived wear rates from fossil and extant proboscidean dental tissues. The models revealed that shifts in tissue topology do not affect whole-tooth wear rate, as inverse trends in lamellar frequency and enamel thickness preserve a consistent surface enamel area fraction; the determining factor of wear. Rather, topological shifts produce a wear-emergent secondary occlusal surface with greater numbers of triturating crests that create a regular, low-relief, file-like shearing pavement. These changes in occlusal architecture likely directly impacted the mastication capacity of Mammuthus dentitions, facilitating their dietary expansion to incorporate fibrous, lower-nutrient graze.

Dental wear↗
DOE OSTI2025

MS25: Materials Science-Focused Benchmark Data Set for Machine Learning Interatomic Potentials

Here, we present MS25, a benchmark data set for evaluating machine learning interatomic potentials (MLIPs) across diverse materials-relevant systems including MgO surfaces, liquid water, zeolites, a catalytic Pt surface reaction, high-entropy alloys (HEAs), and disordered Zr-oxides. Five MLIP architectures (MACE, NequIP, Allegro, MTP, and Torch-ANI) are trained and tested, focusing not only on traditional metrics (energies, forces, and stresses) but also explicitly validating derived physical observables such as lattice constants, volumes, and reaction barriers. We find that most models reach comparable accuracy on standard error metrics across the simple systems, although equivariant MLIPs offer 1.5–2× improvements over nonequivariant MLIPs in energy and force error for structurally complex or compositionally disordered environments such as HEAs and Zr–O systems. Our analysis highlights that low errors in energy and force predictions do not guarantee reliable observables, emphasizing the necessity of explicit validation. We demonstrate limitations in cross-framework transferability, as models trained on one zeolite framework (CHA) fail to reliably generalize to predictions of structurally distinct frameworks (e.g., MFI). Size-extensive tests show some dependence on system size for MgO, resulting from forced periodicity. The HEA and Zr–O data sets are identified as challenging tests for future benchmarks and MLIP model architecture developments as they show significant differentiation in error between MLIP architectures and are still relatively difficult at 1000 training images. Moving forward, we recommend that benchmarking efforts shift their focus from marginal accuracy improvements in energy and force errors toward identifying and understanding model failure modes, rigorously assessing transferability, and evaluating how their errors affect observable predictions. For researchers looking to choose an MLIP architecture, we suggest selecting equivariant MLIP architectures if the complexity of the system is a challenge. For simple materials problems, auxiliary features such as integration with molecular dynamics engines, trade-offs between computational data set generation cost vs MLIP inference speed, and framework integration may play a more important decision factor than small differences in error metrics that are unlikely to matter for production-level research.

chemical structure↗
DOE OSTI2024

Machine learning potential assisted exploration of complex defect potential energy surfaces

Abstract Atomic-scale defects generated in materials under both equilibrium and irradiation conditions can significantly impact their physical and mechanical properties. Unraveling the energetically most favorable ground-state configurations of these defects is an important step towards the fundamental understanding of their influence on the performance of materials ranging from photovoltaics to advanced nuclear fuels. Here, using fluorite-structured thorium dioxide (ThO 2 ) as an exemplar, we demonstrate how density functional theory and machine learning interatomic potential can be synergistically combined into a powerful tool that enables exhaustive exploration of the large configuration spaces of small point defect clusters. Our study leads to several unexpected discoveries, including defect polymorphism and ground-state structures that defy our physical intuitions. Possible physical origins of these unexpected findings are elucidated using a local cluster expansion model developed in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tumor targeted alpha particle therapy with an actinium-225 labelled antibody for carbonic anhydrase IX

Selective antibody targeted delivery of α particle emitting actinium-225 to tumors has significant therapeutic potential. This work highlights the design and synthesis of a new bifunctional macrocyclic diazacrown ether chelator, H 2 MacropaSqOEt, that can be conjugated to antibodies and forms stable complexes with actinium-225. The macrocyclic diazacrown ether chelator incorporates a linker comprised of a short polyethylene glycol fragment and a squaramide ester that allows selective reaction with lysine residues on antibodies to form stable vinylogous amide linkages. This new H 2 MacropaSqOEt chelator was used to modify a monoclonal antibody, girentuximab (hG250), that binds to carbonic anhydrase IX, an enzyme that is overexpressed on the surface of cancers such as clear cell renal cell carcinoma. This new antibody conjugate (H 2 MacropaSq-hG250) had an average chelator to antibody ratio of 4 : 1 and retained high affinity for carbonic anhydrase IX. H 2 MacropaSq-hG250 was radiolabeled quantitatively with [ 225 Ac]Ac III within one minute at room temperature with micromolar concentrations of antibody and the radioactive complex is stable in human serum for >7 days. Evaluation of [ 225 Ac]Ac(MacropaSq-hG250) in a mouse xenograft model, that overexpresses carbonic anhydrase IX, demonstrated a highly significant therapeutic response. It is likely that H 2 MacropaSqOEt could be used to modify other antibodies providing a readily adaptable platform for other actinium-225 based therapeutics.

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

Material Extrusion Printing of Poly(Acrylonitrile‐Styrene‐Acrylate) and Poly(Acrylonitrile‐Butadiene‐Styrene) Structures Reinforced with Poly(Phenylene Oxide) Additives for Improved Thermomechanical Properties and their Surface Analysis by Time‐of‐Flight Secondary Ion Mass Spectrometry

Fused filament fabrication offers the ability to 3D print complex geometries made from plastic filament materials; however, these parts are mechanically outperformed by parts created by traditional fabrication methods. To overcome this challenge, a high-performance polymer poly(2,6-dimethyl-1,4-phenylene oxide) (PPO) is incorporated as an additive into two common engineering thermoplastics, poly(acrylonitrile-styrene-acrylate) (ASA) and poly(acrylonitrile-butadiene-styrene) (ABS). Structures printed from these polymer blends are more mechanically robust compared to those prepared from the parent polymers, with low loading levels (1–5 wt%) of PPO improving the elastic strength by up to ≈30% relative to the parent terpolymers. Even at higher loading levels (10 and 20 wt% PPO), there is no evidence of additive aggregation in the model thin films, which is supported by compositional analysis of the copolymers and chemical analysis via time-of-flight secondary ion mass spectrometry. The enhancements in mechanical properties of ASA and ABS blends appear to be a consequence of homogeneous incorporation of the PPO additive. In conclusion, this work explores expanding materials-property space using miscible blends of engineering thermoplastics to improve mechanical performance as a general approach to overcoming challenges with parts created by melt-based material extrusion printing.

additive manufacturing↗
DOE OSTI2022

Ligand Coverage and Exciton Delocalization Control Chiral Imprinting in Perovskite Nanoplatelets

Chiral perovskites have generated significant recent interest, yet little is known about how their chiro-optical properties arise. In this study, chiral methylammonium lead halide perovskite nanoplatelets (NPLs) with varied halide and ligand compositions are prepared using direct synthetic methods. Circular dichroism (CD) and 1 H NMR studies find a non-linear relationship between the chiroptical properties and the ratio of chiral phenylethylammonium (PEA) to achiral octylamine (OA) ligands on the NPL surface. We use density functional theory (DFT) computations and a chiral imprinted particle-in-a-box model to rationalize the experimentally observed CD spectra and find that the saturation point of the induced chirality depends on the size of the perovskite exciton relative to the size of the ligands. Temperature-dependent CD and 1 H NMR studies, combined with DFT analysis, show that both the CD intensity and sign depend strongly on the structure and orientation of the ligands. Furthermore, this work reveals the complex nature of chiral imprinting in perovskite nanostructures and establishes a simple physical model for ligand-induced chiral imprinting to guide the further development of chiral materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling Spatial Asymmetries in Teleconnected Extreme Temperatures

Abstract Combining strengths from deep learning and extreme value theory can help describe complex relationships between variables where extreme events have significant impacts (e.g., environmental or financial applications). Neural networks learn complicated nonlinear relationships from large datasets under limited parametric assumptions. By definition, the number of occurrences of extreme events is small, which limits the ability of the data-hungry, nonparametric neural network to describe rare events. Inspired by recent extreme cold winter weather events in North America caused by atmospheric blocking, we examine several probabilistic generative models for the entire multivariate probability distribution of daily boreal winter surface air temperature. We propose metrics to measure spatial asymmetries, such as long-range anticorrelated patterns that commonly appear in temperature fields during blocking events. Compared to vine copulas, the statistical standard for multivariate copula modeling, deep learning methods show improved ability to reproduce complicated asymmetries in the spatial distribution of ERA5 temperature reanalysis, including the spatial extent of in-sample extreme events.

Krock, Mitchell L.↗
DOE OSTI2020

Forced Patterns of Sea Level Rise in the Community Earth System Model Large Ensemble From 1920 to 2100

To provide context for observed sea level rise, the forced response (FR) in dynamic sea level (DSL) during the 20th and 21st centuries is examined in the Community Earth System Model Large Ensemble (LE). After accounting for simulation drift, which in the LE is sizable, the DSL FR is found to be complex, both in space and time. Its evolving character is suggested to arise from both the diversity and transient evolution of climate forcing agents and the slow adjustment timescales of the intermediate and deep oceans. Nonetheless, various intervals of spatially coherent change simulated for the recent past and near future are identified, and their characteristics and associated driving mechanisms are identified and discussed. The mid- to late-20th century DSL FR is characterized by an hemispherically asymmetric pattern of change, with depressed rates of rise in the northern oceans. There is also a dipole of change in the Southern Ocean caused by changes in near-surface zonal winds. Through the late 20th and early 21st centuries, a different pattern of rise emerges with elevated rates in the tropics and depressed rates at high latitudes. Zonal and interbasin variations characterize both intervals and involve the pattern and depth of anomalous ocean heat content storage and spatial contrasts in the expansion coefficient tied mainly to base state temperature, with greater rates of rise in warm regions per unit warming. Finally, the relative roles of surface flux and ocean convergence anomalies are examined.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Pyroelectric Heat Detection for Calibrated Measurement of Atomic Layer Deposition Reaction Heat

Atomic layer deposition (ALD) precursors and processes have been identified that enable the deposition of materials that span the periodic table. However, the detailed mechanisms, kinetics, thermodynamics, and complexities of even the most well-studied ALD half-reactions remain insufficiently resolved experimentally to test present computational hypotheses. In this work, a pyroelectric calorimeter is designed, implemented, calibrated, and electrically modeled to enable quantitative ALD reaction heat detection. The operation and response of three ALD calorimeter designs were evaluated for use as absolute ALD half-reaction heat detectors. The most effective design was calibrated in situ before being utilized to measure the heat generated by the reaction of trimethylaluminum with a hydroxylated alumina surface to produce 96 μJ/cm 2 over a reaction time of ~30 ms. With thermal and temporal resolutions down to 0.1 μJ/cm 2 and 50 ns, respectively, this sensitive and ultrafast probe is poised to improve our fundamental understanding of numerous ALD reaction mechanisms through quantitative comparison with first-principles computation predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Correlated Response of Protein Side-Chain Fluctuations and Conformational Entropy to Ligand Binding

The heterogeneous fast side-chain dynamics of proteins plays crucial roles in molecular recognition and binding. Site-specific NMR experiments quantify these motions by measuring the model-free order parameter (O axis 2 ) on a scale of 0 (most flexible) to 1 (least flexible) for each methyl-containing residue of proteins. In this work, we have examined ligand-induced variations in the fast side-chain dynamics and conformational entropy of calmodulin (CaM) using five different CaM–peptide complexes. O axis 2 of CaM in the ligand-free (O axis,U 2 ) and ligand-bound (O axis,B 2 ) states are calculated from molecular dynamics trajectories and conformational energy surfaces obtained using the adaptive biasing force (ABF) method. ΔO axis 2 = O axis,B 2 – O axis,U 2 follows a Gaussian-like unimodal distribution whose second moment is a potential indicator of the binding affinity of these complexes. The probability for the binding-induced O axis,U 2 → O axis,B 2 transition decreases with increasing magnitude of ΔO axis 2 , indicating that large flexibility changes are improbable for side chains of CaM after ligand binding. A linear correlation established between ΔOaxis2 and the conformational entropy change of the protein makes possible the determination of the conformational entropy of binding of protein–ligand complexes. The results not only underscore the functional importance of fast side-chain fluctuations but also highlight key motional and thermodynamic correlates of protein–ligand binding.

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

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗
DOE OSTI2020

Graphene-extracted membrane lipids facilitate the activation of integrin α v β 8

Despite the remarkable electrochemical properties of graphene, strong van der Waals attraction between graphene and biomolecules often causes cytotoxicity, which hinders its applications in the biomedical field. Unfortunately, surface passivation of graphene might stimulate undesired immune response as the nanomaterial triggers cytokine production through membrane receptor activation. In this work, we use all-atom Molecular Dynamics (MD) simulations to unravel the underlying mechanism of graphene-induced inside-out activation of integrin αvβ8, a prominent membrane receptor expressed in immune cells. We model the transmembrane (TM) domains of integrin αvβ8 in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and observe the structural changes in the integrin–membrane complex when interacting with a graphene nanosheet across the membrane. We find that the β8 TM domain interacts with the graphene nanosheet directly or indirectly through extracted lipids, facilitating the pulling of a β8 subunit away from an αv subunit and thus leading to the disruption of the TM domain association by breaking the hydrophobic cluster in the cytoplasmic domains of the αv and β8 subunits. Alanine substitution of two conserved phenylalanine residues on the αv subunit at this hydrophobic cluster further reveals the importance of a stable T-shaped structure in retaining integrin in its inactive state. Our results agree with previous studies on the interactions between other integrin subtypes and their endogenous activators, suggesting an intriguing role that the graphene nanosheet may play in the integrin-related signal transduction during its interaction with the membrane.

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

Detection of multi-reference character imbalances enables a transfer learning approach for virtual high throughput screening with coupled cluster accuracy at DFT cost

Appropriately identifying and treating molecules and materials with significant multi-reference (MR) character is crucial for achieving high data fidelity in virtual high-throughput screening (VHTS). Despite development of numerous MR diagnostics, the extent to which a single value of such a diagnostic indicates the MR effect on a chemical property prediction is not well established. We evaluate MR diagnostics for over 10 000 transition-metal complexes (TMCs) and compare to those for organic molecules. We observe that only some MR diagnostics are transferable from one chemical space to another. By studying the influence of MR character on chemical properties (i.e., MR effect) that involve multiple potential energy surfaces (i.e., adiabatic spin splitting, ΔE H–L , and ionization potential, IP), we show that differences in MR character are more important than the cumulative degree of MR character in predicting the magnitude of an MR effect. Motivated by this observation, we build transfer learning models to predict CCSD(T)-level adiabatic ΔE H–L and IP from lower levels of theory. By combining these models with uncertainty quantification and multi-level modeling, we introduce a multi-pronged strategy that accelerates data acquisition by at least a factor of three while achieving coupled cluster accuracy (i.e., to within 1 kcal mol –1 MAE) for robust VHTS.

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

Impacts of phase state and water content on secondary organic aerosol formation and partitioning (Final Technical Report)

Secondary organic aerosols (SOA) account for a major fraction of particulate matter in the atmosphere, affecting climate, air quality, and public health. SOA formation and evolution are highly complex processes involving both chemical reactions and transport of molecules in air, at particle surfaces, and within particles. The state of the matter, or “phase state,” comprising SOA can vary from a liquid, over an amorphous semi-solid, to a glassy solid, depending on chemical composition, water content, relative humidity (RH), and temperature. The occurrence of glassy and amorphous semi-solid states can pose limitations on the rate of transport of molecules, affecting gas-particle interactions and challenging the treatment of SOA in atmospheric models. The objective of this project was to improve fundamental understanding of the interplay of the phase state of particles and water content on the evolution of SOA formation. We have developed a method to estimate the glass transition temperature (T g ) of organic compounds to predict the phase state and viscosity of SOA. This method has been applied to high resolution mass spectrometry data of various types of SOA, achieving a good agreement with viscosity measurements. In addition, the viscosity estimation method was implemented into a regional model CMAQ to simulate diel and seasonal variations of SOA phase state over the U.S. Applying kinetic multilayer modeling, we analyzed laboratory experiments and ARM field measurements (e.g., HI-SCALE) on kinetic limitations of amine uptake. We have also conducted numerical simulations to estimate equilibration timescales of SOA partitioning. Overall, we have evaluated the impacts of phase state and water content on SOA lifecycle that should contribute to reducing the uncertainty of SOA representation in regional climate and air quality predictions. The project has yielded 17 publications, in which DE-SC0018349 is explicitly acknowledged. Please see the list of publications as below, followed by summary of research activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Zentropy Theory for Transformative Functionalities of Magnetic and Superconducting Materials

The proposed research developed the zentropy theory through applications to complex magnetic materials and superconductors under the hypothesis that the emergent properties of complex magnetic materials and superconductors can be predicted by statistical mechanics of ergodic microstates with their partition functions computed from DFT-predicted free energies. The key objective is to develop approaches to systematically determine the types and number of microstates and the supercell size in DFT-based calculations through convergency of macroscopic functionalities, with the incorporation of our mixed-space approach accounting for the interactions between periodic supercells. In addition to use scientific intuitions to guide the design of important microstates, the key innovation of the proposed research is to integrate the domain knowledge and the material-property-descriptor database (MPDD) with 4 million microstates, which is supported by our deep neural network machine learning models (SIPFENN: structure-informed prediction of formation energy using neural networks) and integrated with our high throughput DFT Tool Kit (DFTTK). For complex magnetic materials, one of the objectives is to develop approaches to calculate short-range ordering from the statistical distribution of each microstate. For superconductors, the divergency of quasiparticle effective mass at a quantum critical point will be investigated, and the superconducting and non-superconducting microstates will be delineated through analysis of electronic band structure, density of states, charge density, and Fermi surface.

36 MATERIALS SCIENCE↗