Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structure optimization”

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

Exsolution of NiCo alloys over Ruddlesden-Popper perovskite for mild electrochemical synthesis of ammonia on protonic ceramic electrochemical cells

Ammonia synthesis from renewable energies on protonic ceramic electrochemical cells (PCECs) shows great potential. The primary challenges in ammonia synthesis on PCECs include sluggish catalytic activity, competition from the hydrogen evolution reaction, and unsatisfactory durability of the cathode, which is the active site of ammonia generation. Here, in this study, we report an elaborate design of the cathode with an intended formula of Pr 4 Ni 1.79 Co 1.2R u 0.01 O 10-δ , where NiCo alloy nanoparticles are exsolved from Ruddlesden-Popper perovskite substrates after the reduction in 5 % H 2 /Ar at 400 °C for 1 h, for electrocatalysis of the nitrogen reduction reaction to ammonia. The host material Pr 4 Ni 1.8 Co 1.2 O 10-δ with embeddable layered structure and stability was deliberately chosen to fix the Ru cation and maximize the catalytic activity of NiCo. Also, density functional theory calculations suggest that Ru doping provides an optimal balance between structural stability and redox activity, facilitating the controlled exsolution of NiCo nanoparticles and enhancing catalytic performance. As a result, the composite electrode with exsolved NiCo alloy and abundant oxygen vacancies on fuel-electrode-supported PCECs achieves a superior electrochemical activity towards ammonia synthesis: a peak ammonia formation rate of 27.84 μg h −1 cm −2 and excellent Faradaic efficiencies of 62.6 % at 350 °C.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Crystal structure prediction with host-guided inpainting generation and foundation potentials

Unconditional crystal structure generation with diffusion models faces challenges in identifying symmetric crystals as the unit cell size increases. Here, we present the crystal host-guided generation (CHGGen) framework to address this challenge through conditional generation using an inpainting method, which optimizes a fraction of atomic positions within a predefined and symmetrized host structure to improve the success rate for symmetric structure generation. By integrating inpainting structure generation with a foundation potential for structure optimization, we demonstrate the method on the ZnS–P 2 S 5 and Li–Si chemical systems, where the inpainting method generates a higher fraction of symmetric structures than unconditional generation. The practical significance of CHGGen extends to enabling the structural modification of crystal structures, particularly for systems with partial occupancy or intercalation chemistry. The inpainting method also allows for seamless integration with other generative models, providing a versatile framework for accelerating materials discovery.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2025

Dark-field X-ray microscopy with structured illumination for three-dimensional imaging

Dark-field X-ray microscopy is a lens-based technique that enables real-space imaging of heterogeneous micro- and meso-scale ordered materials. However, achieving accurate three-dimensional (3D) reconstruction often requires meticulous sample alignment or rastering, requiring complex rotational setups and extended acquisition times. To address these challenges, we introduce a structured illumination technique optimized for 3D imaging of ordered materials at sub-micrometer length scales. Our approach employs a coded aperture to spatially modulate the incident X-ray beam, enabling 3D structural reconstruction from images captured at various aperture positions. Unlike current 3D imaging approaches, which often rely on rotational or rastering methods, our technique uses scanning X-ray silhouettes of the coded aperture for depth resolution along the diffraction axis. This eliminates the need for sample rotation or rastering, resulting in a highly stable and efficient imaging modality. We validated the efficacy of this approach through experimental imaging of an isolated twin domain within a bulk single crystal of an iron pnictide using a dark-field X-ray microscope. This advancement aligns with the enhanced brightness upgrades of modern synchrotron radiation facilities, unlocking new possibilities for high-resolution imaging of ordered materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Elucidating metal–organic framework structures using synchrotron serial crystallography

Metal organic frameworks (MOFs) are porous crystalline materials that display a wide variety of physical and chemical properties. Their single crystal structure determination is often challenging because in most cases micro- or nano-sized crystals spontaneously form upon MOF synthesis, which cannot be recrystallized. The production of larger single crystals for structure determination involves optimizing, and thus modifying, the conditions of synthesis, in which success cannot be guaranteed. Failure to produce crystals suitable for single-crystal X-ray diffraction leaves the 3D structure of the MOF compound unknown, and scientists must resort to more challenging structure solution methods based on X-ray powder or electron diffraction data. These laborious tasks can be avoided by using serial crystallography techniques which merge data collected on many micro-crystals. Here, we report the application of three synchrotron serial crystallography methods. We call these “mesh”, “grid” and “mesh&collect” scans. “Still” images (no rotation) are collected in the mesh scan approach, whereas small rotational wedges are collected in the grid scan method. The third protocol, mesh&collect, combines the acquisition of still images and rotational wedges. Using these means, we determine the ab initio structure of benchmark MOFs, MIL-100(Fe) and ZIF-8, that differ largely in unit cell size. These methods are expected to be widely applicable and facilitate structure determination of many MOF microcrystalline systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Heterovalent Substitution of K 2 SrP 2 O 7 :Cr 3+ to Achieve Anti-Thermal-Quenching Broadband Near-Infrared Luminescence

Broadband near-infrared (NIR) light sources based on phosphor-converted light-emitting diodes are highly desirable for biochemical analysis and medical diagnosis applications. However, thermal quenching remains a demanding challenge for developing efficient NIR phosphors. Herein, we report the enhancement of both quantum efficiency and thermal stability in Cr 3+ -activated K 2 SrP 2 O 7 phosphors through a heterovalent substitution strategy by replacing Sr 2+ with Al 3+ in K 2 Sr 1–x Al x P 2 O 7 (0.05 ≤ x ≤ 0.2) to obtain optimized broadband NIR emission. Structural modulation via Al 3+ substitution leads to the optimized composition, K 2 Sr 0.88 Al 0.1 P 2 O 7 :0.02Cr 3+ , which emits across a broad NIR range of 650–1100 nm peaking at 807 nm with a full width at half-maximum of ∼130 nm under 448 nm excitation. Remarkably, its emission intensity at 150 °C remains 120% of the initial value at room temperature, demonstrating a rare antithermal-quenching behavior. Temperature-dependent XRD studies further reveal that Al 3+ substitution effectively suppresses lattice expansion at elevated temperatures, indicating enhanced lattice stability under thermal excitation. Detailed structural and spectral analyses show that the substitution enhances local site symmetry, reduces electron–phonon coupling, increases thermally induced absorption probability, and fortifies energetic barriers against nonradiative transitions. These synergistic effects collectively endow this NIR phosphor with a superior thermal stability. Furthermore, NIR light-emitting diodes fabricated with this phosphor exhibit strong potential for applications in information identification, nondestructive detection, and night vision technologies. This study demonstrates a local structure engineering strategy for designing thermally robust Cr 3+ -activated NIR phosphors, offering valuable insights into material discovery and NIR spectroscopy device development.

Cr3+↗
DOE OSTI2026

From Powder to Power: Tailored Pre-Milling Strategy that Optimizes Microstructure for Efficient Hydrogen Evolution

The hydrogen evolution reaction (HER) plays a critical role in enabling large-scale electrolytic hydrogen production and advancing future technologies and fuel production. Among non-precious metals, NiMo-based catalysts are particularly attractive due to their capability to promote both water dissociation and hydrogen adsorption in alkaline media. However, conventional NiMo catalysts often suffer from incomplete alloying, particle aggregation, weak metal–support interactions, and surface oxidation, which significantly limit active site utilization, electronic conductivity, and long-term stability. Herein, we report a scalable, solid-state, two-step ball milling strategy for constructing an efficient NiMo/C catalyst for the HER. Pre-milling the metal precursors prior to carbon incorporation promotes efficient solid-state activation, leading to the formation of uniformly alloyed, ultrafine, and defect-rich Ni–Mo nanoparticles with robust metal–carbon interfacial anchoring. Benefiting from this integrated structural design, the optimized NiMo/C catalyst exhibits low overpotentials and Tafel slopes, reduced charge-transfer resistance, and excellent durability under alkaline conditions. Comprehensive structural and surface analysis reveal that the enhanced HER performance can be attributed to the synergistic interplay of homogeneous Ni–Mo alloying, uniform nanoparticle dispersion on a defect-rich carbon scaffold, enhanced accessibility of catalytically active sites, and optimized electronic coupling that stabilizes the catalyst surface. This work highlights the two-step solid-state ball milling strategy as a simple, robust, and scalable route for the preparation of effective and durable non-precious metal HER electrocatalysts, offering practical insights toward large-scale and inexpensive hydrogen production.

Yang, Xiaoxuan [Oak Ridge National Laboratory (ORN↗
DOE OSTI2026

Developing cathode infiltration processes for all-solid-state bilayer LLZO cells

To realize the bilayer architecture of lithium lanthanum zirconate (LLZO) for application in solid-state batteries (SSBs), the scaffold structure must be optimized, and effective cathode infiltration strategies must be established. In this study, we fabricate a modified bilayer LLZO using a sacrificial layer to enhance surface porosity, and systematically investigate various cathode infiltration techniques to fill the scaffold with oxide cathode active materials (CAM). Structural characterizations showed that the sacrificial layer significantly increased open surface porosity, enabling the surface of the scaffold to be filled with CAM. To further increase infiltration depth, applying vacuum or vibration was compared, with the full-depth infiltration achieved using a sonicator-based vibration. Full cells prepared using the modified bilayer LLZO and vibration-assisted technique demonstrated successful operation. This work demonstrates a practical and scalable approach for engineering bilayer LLZO structures and integrating oxide cathodes into porous scaffolds, offering a promising pathway toward high-performance solid-state batteries.

Bilayer↗
DOE OSTI2024

Active Learning for Metamaterial Optimization on HPC and QC Integrated Systems

Active learning algorithms, integrating machine learning, quantum computing and optics simulation in an iterative loop, offer a promising approach to optimizing metamaterials. However, these algorithms can face difficulties in optimizing highly complex structures due to computational limitations. High-performance computing (HPC) and quantum computing (QC) integrated systems can address these issues by enabling parallel computing. In this study, we develop an active learning algorithm working on HPC-QC integrated systems. We evaluate the performance of optimization processes within active learning (i.e., training a machine learning model, problem-solving with quantum computing, and evaluating optical properties through wave-optics simulation) for highly complex metamaterial cases. Our results showcase that utilizing multiple cores on the integrated system can significantly reduce computational time, thereby enhancing the efficiency of optimization processes. Therefore, we expect that leveraging HPC-QC integrated systems helps effectively tackle large-scale optimization challenges in general.

Kim, Seongmin↗
DOE OSTI2025

Micropolar deep material network

This study extends the Deep Material Network (DMN), a physics-informed machine learning framework, to predict the homogenized mechanical response of composite materials with micropolar (Cosserat-type) constitutive behavior. This extension incorporates microstructure-dependent size effects, enabling accurate, efficient, and size-aware predictions for composites with complex internal architectures. While traditional, direct numerical simulation micropolar models effectively capture size effects by introducing extra local degrees of freedom, they bring significant computational challenges, particularly for multiscale analyses relevant to engineering applications. The micropolar DMN developed in this paper achieves high accuracy while significantly reducing computation time compared to micropolar direct numerical simulations. This advancement enables multiscale analyses and parameter studies that were previously impractical, such as high-cycle fatigue simulations and comprehensive investigations of internal length scale effects notably in size-dependent plastic response and the optimization of lattice structures. By uniting microstructure-sensitive modeling, physics-driven learning, and scalable surrogate modeling, the micropolar DMN paves the way for accelerated material design, large-scale parametric studies, and the reliable incorporation of size-dependent effects across a wide range of engineering applications, including optimization and next-generation composite design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling-Based Design and Optimization of a Gradient Composite Transition Joint

An innovative additively manufactured gradient composite transition joint (AM-GCTJ) has been designed to join dissimilar metals, to address the pressing issue of premature failure observed in conventional dissimilar metal welds (DMWs) when subjected to increased cyclic operating conditions of fossil fuel power plants. The transition design, guided by computational modeling, developed a gradient composite material distribution, facilitating a smooth transition in material volume fraction and physical properties between different alloys. This innovative design seeks to alleviate structural challenges arising from distinct material properties, including high thermal stress and potential cracking issues resulting from the thermal expansion mismatch typically observed in conventional DMWs. In this study, we investigated the creep properties of transition joints comprising Grade 91 steel and 304 stainless steel through a combination of simulations and creep testing experiments. The implementation of a gradient composite design in the plate transition joint resulted in a significant enhancement of creep resistance when compared to the baseline conventional DMW. For instance, the creep rupture life of the transition joint was improved by > 400% in a wide range of temperature and stress testing conditions. Meanwhile, the failure location shifted to the base material of Grade 91 steel. Such enhancement can be primarily attributed to the strong mechanical constraint facilitated by the gradient composite design, which effectively reduced the stresses on the less creep-resistant alloy in the transition zone. Beyond examining plate joints, it is crucial to assess the deformation response of tubular transition joints under pressure loading and transient temperature conditions to substantiate and demonstrate the effectiveness of the design. The simulation results affirm that the tubular transition joint demonstrates superior resistance compared to its counterpart DMW when subjected to multiaxial stresses in tubular structures. In addition, optimization of the transition joint’s geometry dimensions has been conducted to diminish the accumulated deformation and enhance the service life. Lastly, the scalability and potential of the innovative transition joints for large-diameter pipe applications are addressed.

Zhang, Wei↗
DOE OSTI2025

High‐Performance Al–S Batteries by Spin Polarization Modulation via Catalytic Ni‐MoS 2 Nanosheets

Aluminum–sulfur (Al–S) batteries catalysts with adsorption and catalytic capabilities can effectively improve the slow redox kinetics, but the current research often ignores the effect of optimizing the electronic structure of the catalyst on improving charge transfer and adsorption. Here, in this study, Ni-doped monolayer MoS 2 nanosheets are synthesized and used as a catalytic additive for the sulfur cathode. The addition of Ni promotes spin splitting of 4d orbital of Mo, thereby affecting polarization degree of the basal plane sulfur and making it change from a low spin state to a high spin one. This high spin configuration raises the electron energy level and provides an active electron state to react with aluminum polysulfides (AlPSs), which optimizes the adsorption energy. At the same time, it accelerates electron transfer and lowers the energy barrier for the overall conversion of the polysulfides. Benefiting from these features, Al–S batteries based on rationally designed S@Ni-MoS 2 /C cathodes exhibit a high initial capacity (1603.0 mAh g −1 at 0.5 A g −1 ) and extraordinary cycling stability (0.035% capacity decay rate during 2000 cycles). This study showcases a spin-polarized electronic structure control strategy to enhance catalytic activity, providing a viable approach for developing efficient catalysts for practical Al–S batteries.

Aluminum-sulfur batteries↗
DOE OSTI2024

TuFF internal WRAP for Rapid Pipeline Repair (TuFF iWRAP)

The goal of “TuFF internal WRAP for Rapid Pipeline Repair” (TuFF iWRAP) program was to develop a novel material system and placement process to fabricate structural pipe within the existing deteriorated pipelines without disruption of gas delivery. The team (University of Delaware – Center for Composite Materials (UD-CCM) and Plitzie Inc.) addressed this challenge by developing a new material feedstock and pipe in pipe (PIP) repair strategy. This allows the potential for significant cost reduction and has minimum operational impact on gas customers. The new robotic based placement design allows discontinuous placement of pipe sections creating a stand-alone structural liner within the legacy pipeline without the need for pipe shutdown. Here, the material is supplied using a tethered material feeding system and is placed and UV cured with the internal Wound Rapid Automated Placement (iWRAP) system. This provides maximum placement efficiency capable of traversing 90 angle bends in 12-inch pipe and overall design customization to meet pipe repair requirements (e.g., variable wall thickness, bridging gaps, etc.). UV-curable fiber reinforced composite material has been optimized to meet structural performance and placement/cure times. Superior strength and fatigue life has been demonstrated by improving fiber-matrix adhesion using new fiber sizing for UV resins. Rapid cure approaches using new liner and resins have been evaluated with industry. The appropriate design of the section joints has been developed and tested. We estimate coating time to be ~100 hours per mile enabling typical pipe repair within 1 week.

03 NATURAL GAS↗
DOE OSTI2025

Orthorhombic Structures as Inorganic Halide Perovskite Models for High-Throughput DFT Investigations

Although high-symmetry Pm-3m space group cubic models are computationally efficient for high-throughput density functional theory (DFT) calculations of inorganic ternary (ABX 3 ) halide perovskites (HPs), they frequently predict band gaps (E g ) that disagree with experiment. Conversely, while low-symmetry cubic polymorphous networks (PN) comprised of 160 to 320 atoms incur significantly greater computational cost, they predict E g ’s that are more closely aligned with experiment. Here, in this study, we compare the DFT total energies and E g ’s predicted by four high-symmetry structure models (Pnma orthorhombic, R3m trigonal, P4/mbm tetragonal, and Pm-3m cubic) to cubic PNs for 5 experimentally characterized ternary HPs and find that the orthorhombic model computes Eg’s with the smallest MAD of 0.23 eV relative to the PNs. Pair distribution functions and DFT-computed total energies show that octahedral tilting, which is present in the 20-atom orthorhombic and 160-atom cubic PN models but not in the 5-atom cubic models, stabilizes all 5 compositions in our benchmarking set. We also find that imposing PN constraints when generating and optimizing these orthorhombic structures by fixing the unit cell lattice vectors and displacing the atoms prior to ionic relaxation with DFT, which we call the orthorhombic surrogate model (OSM), lowers the MAD of E g predictions to 0.09 eV. Our OSM predicts the PN band gaps of an additional 95 theoretical inorganic ternary HPs with MAD of 0.08 eV, supporting its usage in high-throughput DFT investigations to closely estimate PN band gaps with much less computational expense.

14 SOLAR ENERGY↗
DOE OSTI2025

Quasiaperiodic grain boundary phases of Σ⁢5 tilt grain boundaries in refractory metals

We report ground-state structures and phase transitions in Σ⁢5⁢[001] tilt grain boundaries (GBs) in body-centered-cubic (bcc) refractory metals Nb, Ta, Mo, and W. Σ⁢5 tilt GBs have been extensively investigated over the past several decades, with their ground-state structure—composed of kite-shaped structural units—previously thought to be well understood. By performing a rigorous GB structure search that optimizes the number of atoms in the boundary core, we predict different quasiaperiodic “split kite” phases analogous to those previously found in GBs in face-centered-cubic metals. Furthermore, our results suggest that complex aperiodic phases of GBs appear to be a general phenomenon, as validated through density functional theory calculations. Moreover, the atoms in the split kite phase demonstrate distinct collective diffusion dynamics. Phase-contrast image simulations of split kites show better agreement with experimental observations, offering an alternative explanation for previous microscopy results and motivating future atomically resolved imaging of the GB structure.

Body-centered cubic↗
DOE OSTI2024

Decomposition and Algorithmic Approaches for Solving Large-Scale Process Family Design Problems

Our most recent work expands the water desalination case study from 76 variants to 10,897 variants using the equation-oriented model built in Pyomo as part of the PARETO project. Using the discretization formulation presented in Stinchfield (2024a), rather than solving for all 10,897 variants simultaneously, we decompose the formulation into subproblems containing subsets of variants from the process family. We solve the overall problem with Progressive Hedging (PH) deployed in parallel on a distributed HPC cluster using the open-source Python package mpi-sppy (Knueven et al., 2023). This approach allowed us to solve this process family design problem to ~1.5% relative optimality gap in about 5 hours; in comparison, Gurobi reached ~50% relative optimality gap in about 6 hours (Stinchfield et al., 2024b). However, this approach still requires discretization of the common unit module design ranges; additionally, PH acts as a heuristic for MILP’s with gap-closing capabilities. Ideally, we would not have to use ML surrogates or discretization to solve this problem, instead solving the process family design problem with the equation-oriented model directly to achieve the most accurate results. However, recall that we did not consider solving the MINLP directly due to complexity and size. In this work, we aim to decompose and solve this large-scale MINLP using a Structured Nonlinear Global Optimization algorithm presented by Cao and Zavala (2019).

Stinchfield, Georgia↗
DOE OSTI2024

Durable and High-Performance SOECs Based on Proton Conductors for Hydrogen Production

Proton-conducting solid oxide electrolysis cells (P-SOECs) are a promising technology for cost-effective and efficient production of green hydrogen. Breakthroughs in materials development, optimization of cell structure, and achievement of high performance and durability are essential to significantly increase the commercial competitiveness of these technologies. The main objective of this project is to gain scientific knowledge for the rational design, fabrication, and demonstration of a robust, highly efficient, and low-cost SOEC technology based on a proton-conducting electrolyte membrane for hydrogen production. We focused on better understanding the degradation mechanisms of proton-conducting electrolytes, air electrodes, and catalyst materials under electrolysis mode to develop an effective strategy for rationalizing new materials that are vital for enhancing cell performance and durability. The scope includes enhancing the performance and durability of the electrolyte and electrode materials under realistic operating conditions, developing highly active and robust catalysts to minimize electrode losses while improving tolerance to contaminant poisoning, revealing the mechanism of enhanced activity and stability of the catalyst, and understanding the underlying degradation mechanisms. In addition, various characterization techniques were employed to gain a fundamental understanding of the materials’ behavior and their impact on cell performance, providing vital information to guide materials discovery and cell design. After defect chemistry engineering, the optimized donor and acceptor co-doped electrolytes BaMo/W 0.03 Ce 0.71 Yb 0.26 O 3-δ (BM/W03) showed substantially improved chemical stability against high concentrations of CO 2 and H 2 O compared to the state-of-the-art electrolyte (BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ , BZCYYb1711) while maintaining comparable ionic conductivity and ionic transference number. To bypass the inherent trade-off between conductivity and chemical stability, we fabricated a bi-layer electrolyte composed of BZCYYb1711 coated with a highly-stable thin layer of BaHf 0.83 Yb 0.17 O 3-δ (BHYb). This bi-layer electrolyte displayed excellent chemical stability against high concentration CO 2 ; there was no detectable formation of BaCO 3 after exposure to 97% CO 2 (with 3% H 2 O) at 500 °C for 1000 hours and the rate of degradation in resistance was about 0.4% per 1,000 hours (kh). In contrast, the same BZCYYb1711 electrolyte without a BHYb coating degraded significantly under the same testing conditions; the degradation rate was increased to 5.1%/kh. In addition, a triple conducting air electrode Ba 0.9 Pr 0.1 Hf 0.1 Y0.1Co 0.8 O 3-δ (BPHYC) was developed by heavily doping transition metal ions into a proton-conducting material. This air electrode material, composed of 3 distinct phases, exhibits superior electrocatalytic activity due to the synergistic effect from the three component phases. Moreover, an active and durable catalyst, La 2 Ni 0.5 Fe 0.5 O 4+δ (LNF), was developed, showing excellent catalytic activity and contaminant tolerance, with a degradation rate of only 0.49%/kh when exposed to high concentrations of steam and Cr. Finally, single cells were constructed from the best electrolytes, electrodes, and catalyst coatings developed in this project. These cells demonstrated superior high current density at a given cell voltage, high roundtrip efficiency, and remarkable durability (up to 1000 hours of operation).

08 HYDROGEN↗
DOE OSTI2025

Bulk‐Boundary Correspondence of Semimetal Ru 3 Sn 7 and Topological Surface States on Chemically Realistic Terminations

Ru 3 Sn 7 is experimentally demonstrated as an efficient catalyst, with potential utilization of topological surface states for hydrogen evolution reaction. Despite its promising catalytic performance, the topological nature of Ru 3 Sn 7 remains uncertain. Particularly, the bulk-boundary correspondence has not yet been established, hence hindering a rigorous justification of its topologically-protected surface states. In this work, the bulk topology of Ru 3 Sn 7 is detailed using first-principles calculations and the topological quantum chemistry formalism. Ru3Sn7 turns out to be an enforced semimetal possessing symmetry-protected crossings within a set of bands near the Fermi level, which are enforced and prescribed by the violations of symmetry-prescribed compatibility relations. Moreover, the surface states and the associated origin from the same set of entangled bands are identified, thereby establishing the bulk-boundary correspondence. To evaluate the effects of chemical modifications, the response of topological surface states to various surface terminations, stoichiometry, and oxidation is examined. The surface structures are globally optimized, and the phase diagrams for various experimental conditions are built. It is shown that, due to changes in the local chemical environment, the original surface states are significantly altered. Modified surface bands can be observed near the Fermi level on surface terminations that preserve the C 4v symmetry.

36 MATERIALS SCIENCE↗