SEARCH · Search NASA

Results for “Non-linear Simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

842 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Aerothermodynamic CFD Analysis

This presentation provides an introduction to two Computational Fluid Dynamics (CFD) codes, LAURA and FUN3D, which have been developed and widely used at NASA Langley Research Center in Entry, Descent, and Landing applications. A workflow using the LAURA and FUN3D CFD codes to predict aerodynamic and aerothermodynamic engineering quantities is presented. The presentation will cover the basics of setting up and running simulations using both CFD codes, as well as how to post-process results obtained from each code. More advanced topics will also be presented, including a state-of-the-art uncertainty quantification approach to be included in the next LAURA release, and a new geometry-based workflow used with FUN3D focused on replacing user-defined meshing with an adaptation-based approach.

Kyle B Thompson↗

DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗

DOE OSTI2026

Used Nuclear Fuel Management Using the Next Generation System Analysis Model

The U.S. Department of Energy (DOE) is leading the National effort to manage the back end of the nuclear fuel cycle, encompassing the safe transportation, storage/staging, and/or eventual disposal of used nuclear fuel (UNF) and high-level radioactive waste. The Next Generation System Analysis Model (NGSAM) is DOE’s discrete-event, agent-based simulation tool designed to model the full life cycle of UNF from reactor discharge to final disposal. NGSAM supports the DOE Office of Spent Fuel and High-Level Waste Disposition by enabling a detailed, scenario-based analysis of logistics, infrastructure, and shipping strategies. NGSAM replaces legacy models with a modern, flexible platform built on Repast Simphony and enhanced by the Process Analysis Tool. NGSAM simulates the movement and interaction of individual fuel assemblies with system components such as canisters, casks, railcars, and facilities. The model integrates with the Java Transportation Operations Model to plan and execute transportation scenarios, supporting both constrained and unconstrained resource allocation. Key features include customizable allocation and acceptance algorithms, detailed facility-level operations, and a Quick Edit tool for rapid scenario adjustments. NGSAM supports multimodal transportation modeling (e.g. rail, road, barge) and provides comprehensive cost, schedule, and infrastructure data. NGSAM utilizes data from sources such as DOE’s STANDARDS UNF database and DOE’s Stakeholder Tool for Assessing Radioactive Transportation, while also allowing user-defined inputs for scenario customization. NGSAM enables stakeholders to evaluate complex UNF management strategies, assess system performance under varying assumptions, and inform decision making for future infrastructure investments. Its modular architecture and integration with other Integrated Waste Management System tools make it a critical asset for planning the safe and efficient disposition of the Nation’s growing UNF inventory.

Craig, Brian [Argonne National Laboratory (ANL)]↗

DOE OSTI2024

Making Plasticized Polymer Electrolytes Stable Against Sodium Metal for High‐Energy Solid‐State Sodium Batteries

Solid polymer electrolytes based on plastic crystals are promising for solid-state sodium metal (Na 0 ) batteries, yet their practicality has been hindered by the notorious Na 0 -electrolyte interface instability issue, the underlying cause of which remains poorly understood. Here, in this study, by leveraging a model plasticized polymer electrolyte based on conventional succinonitrile plastic crystals, we uncover its failure origin in Na 0 batteries is associated with the formation of a thick and non-uniform solid electrolyte interphase (SEI) and whiskery Na 0 nucleation/growth. Furthermore, we design a new additive-embedded plasticized polymer electrolyte to manipulate the Na 0 deposition and SEI formulation. For the first time, we demonstrate that introducing fluoroethylene carbonate (FEC) additive into the succinonitrile-plasticized polymer electrolyte can effectively protect Na 0 against interfacial corrosion by facilitating the growth of dome-like Na 0 with thin, amorphous, and fluorine-rich SEIs, thus enabling significantly improved performances of Na//Na symmetric cells (1,800 h at 0.5 mA cm −2 ) and Na//Na 3 V 2 (PO 4 ) 3 full cells (93.0 % capacity retention after 1,200 cycles at 1 C rate in coin cells and 93.1 % capacity retention after 250 cycles at C/3 in pouch cells at room temperature). Our work provides valuable insights into the interfacial failure of plasticized polymer electrolytes and offers a promising solution to resolving the interfacial instability issue.

25 ENERGY STORAGE↗

NASA NTRS2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust↗

NASA NTRS2026

Performance Optimization for Lunar Extravehicular Activity Readiness (POLAR) Study: Methods Paper

To better understand which aspects of physical fitness may be most related to performance during Lunar surface operations and thereby help to inform the current NASA fitness standards, much can be learned from fields encompassing the “tactical athlete.” Other physically demanding professions such as law enforcement, military, or rescue professionals often require candidates to meet occupationally-relevant fitness standards. The determination of such standards is a multistep process, including both objective and subjective measures, to determine tasks essential to occupational performance and identify the minimal fitness profile needed to meet physical demands of the job. Notably, fitness is only one component which may contribute to the demands of astronaut selection, flight assignment, and occupational performance. Utilizing a framework to systematically determine which domains of fitness most contribute to relevant job tasks can aid in the refinement of current NASA-STD-3001 fitness standards. Therefore, NASA’s Exercise Physiology & Countermeasures Laboratory conducted the Performance Optimization for Lunar Extravehicular Activity Readiness (POLAR) study to identify and examine a comprehensive list of fitness tests (including NASA-STD-3001 assessments: 1-Repetition Maximum [1-RM] bench press and deadlift) and determine preliminary relationships between identified fitness parameters and novel Artemis-relevant tasks to help inform future investigations for the continued development of aerobic and muscular fitness standards for surface EVAs. This was accomplished through 1) a review of the literature relating fitness assessments to simulated or real EVA performance to identify fitness tests that are most correlated with simulated EVA task performance; 2) a task analysis following a modified framework for physical employment standards development to down select mission critical tasks; and 3) development and pilot testing of a novel, portable Artemis-relevant EVA task circuit to relate to a battery of fitness assessments. The current report describes the methodology used to complete the task analysis, EVA task circuit development, and the pilot study.

Nicole C Strock↗

NASA NTRS2026

Remote Sensing-Driven Hydrodynamic Modeling in Data-Scarce Regions: Integrating ICESat-2, Sentinel-2, SWOT and Re-analysis Models for Coastal Monitoring

Hydrodynamic models in coastal and estuarine systems are typically constrained by sparse bathymetry, boundary, and validation data, especially in regions where field campaigns are costly or impractical. Here we develop and test a fully satellite-driven framework for hydrodynamic modeling in South Africa’s Langebaan Lagoon without using any local in situ measurements. Bathymetry is derived by training multispectral Sentinel-2 reflectance against ICESat-2 ATL24 photon-derived depths using an XGBoost model optimized with Bayesian search. The final satellite-derived bathymetry reproduces independent ATL24 points with RMSE = 0.45 m and R 2 = 0.97. This bathymetry was used in a depth-averaged Delft3D Flexible Mesh model driven at the open boundary by TPXO tidal harmonics and by ERA5 winds. We validate modeled water surface elevation against 16 SWOT low-rate (250 m, unsmoothed) passes in 2023. SWOT–model comparisons yield an overall RMSE of 0.11 m and R 2 = 0.61, with typical point differences <0.10 m (∼7% of the 1.5 m tidal range), and showed consistent spatial gradients in water level from the offshore boundary, through Saldanha Bay, and into the lagoon. At the offshore boundary, TPXO and SWOT sea surface heights agree closely (R 2 = 0.86). A simple phase adjustment of ∼26,min between TPXO and SWOT lowers the RMSE from 0.18,m to 0.11,m, showing that phase offset accounts for some of the discrepancy, with additional errors likely linked to non-tidal signals. Our results demonstrate that combining passive optical, photon-counting LiDAR, radar interferometry, and global tidal/atmospheric models enables robust, transferrable hydrodynamic modeling in data-scarce coastal systems, offering a cost-effective pathway for monitoring.

ICESat-2↗

NASA NTRS1979

Dynamic analysis of a photovoltaic power system with battery storage capability

A photovolataic power system with a battery storage capability is analyzed. A dual battery current control concept is proposed, which enables the battery to either supply or accept power depending upon system environment and load conditions. A simulation of the power system, including the battery current control, is developed and evaluated. The evaulation demonstrate the visbility of the battery control concept of switch the battery from a charge to discharge mode and back as required by load and environmental conditions. An acceptable system operation is demonstrated over the entire insolation range. Additionally, system sensitivity, bandwidth, and damping characteristics of the battery control are shown to be acceptable for a projected hardware implementation.

Merrill, W. C.↗

DOE OSTI2026

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality↗

DOE OSTI2024

Thermally unstable roosts influence winter torpor patterns in a threatened bat species

Abstract Many hibernating bats in thermally stable, subterranean roosts have experienced precipitous declines from white-nose syndrome (WNS). However, some WNS-affected species also use thermally unstable roosts during winter that may impact their torpor patterns and WNS susceptibility. From November to March 2017–19, we used temperature-sensitive transmitters to document winter torpor patterns of tricolored bats (Perimyotis subflavus) using thermally unstable roosts in the upper Coastal Plain of South Carolina. Daily mean roost temperature was 12.9 ± 4.9°C SD in bridges and 11.0 ± 4.6°C in accessible cavities with daily fluctuations of 4.8 ± 2°C in bridges and 4.0 ± 1.9°C in accessible cavities and maximum fluctuations of 13.8 and 10.5°C, respectively. Mean torpor bout duration was 2.7 ± 2.8 days and was negatively related to ambient temperature and positively related to precipitation. Bats maintained non-random arousal patterns focused near dusk and were active on 33.6% of tracked days. Fifty-one percent of arousals contained passive rewarming. Normothermic bout duration, general activity and activity away from the roost were positively related to ambient temperature, and activity away from the roost was negatively related to barometric pressure. Our results suggest ambient weather conditions influence winter torpor patterns of tricolored bats using thermally unstable roosts. Short torpor bout durations and potential nighttime foraging during winter by tricolored bats in thermally unstable roosts contrasts with behaviors of tricolored bats in thermally stable roosts. Therefore, tricolored bat using thermally unstable roosts may be less susceptible to WNS. More broadly, these results highlight the importance of understanding the effect of roost thermal stability on winter torpor patterns and the physiological flexibility of broadly distributed hibernating species.

Biodiversity & Conservation↗

NASA NTRS2006

Battery Control Boards for Li-Ion Batteries on Mars Exploration Rovers

Rechargeable Lithium-ion batteries have been operating successfully on both Spirit and Opportunity rovers for the last two years, which includes six months of Assembly Launch and Test Operations (ATLO), seven months of cruise and about eleven months of surface operations. The Battery Control Boards designed and fabricated in-house would protect cells against overcharge and over-discharge and provide cell balance. Their performance has thus far been quite satisfactory. The ground data o the mission simulation battery project little capacity loss of less than 3% during cruise and 180 sols. Batteries are poised to extend the mission beyond six months, if not a couple of years.

Ewell, R.↗

DOE OSTI2026

Rinse-Free, Sodium-Efficient Synthesis of O3-Type Layered Oxide Materials Enabled by Acetate Precursors

Sodium-ion batteries (SIBs) are a sustainable alternative to lithium-ion systems for global electrification, with O3-type layered oxide cathodes offering high specific capacity and feasibility of scalable synthesis. Industrial co-precipitation synthesis of these cathodes typically uses transition metal sulfates, requiring extensive water rinsing to remove Na 2 SO 4 impurities, a process that consumes significant water and risks residual inactive phases if incomplete. Here, this work introduces a rinse-free, resource-efficient approach using metal acetate precursors. Residual sodium acetate in non-rinsed precursors decomposes during sintering to generate Na 2 CO 3 in situ, partially substituting an external sodium resource (e.g., NaOH and Na 2 CO 3 ) and reducing its consumption by ∼18–20%. Phase-pure O3-Na 1.0 Ni 1/3 Fe 1/3 Mn 1/3 O 2 (NFM111) cathodes synthesized via this method exhibit microstructure and electrochemical performance comparable to rinsed sulfate-derived counterparts, with initial capacities of 141 mAh g −1 at C/20 (7.5 mA g −1 ). By eliminating rinsing and minimizing sodium reagent use, this acetate-based route enhances sustainability and scalability of layered oxide production for SIBs.

acetate vs. sulfate↗

DOE OSTI2025

Adsorptive behavior of poly (vinylidene fluoride) membranes for the recovery of lignin-derived hydrophobic deep eutectic solvents

Abstract Recently, membrane technology has gained significant traction as an energy-efficient alternative to traditional thermal processes for solvent recovery. Deep eutectic solvents (DESs) have emerged as sustainable alternatives to conventional organic solvents, yet a systematic methodology for selecting compatible membrane materials for their recovery remains underdeveloped. This study established a predictive framework for membrane material selection in hydrophobic DES applications using Hansen Solubility Parameters (HSP) with inverted criteria targeting materials with relative energy difference (RED) values greater than 1.0. Flat sheet membranes were fabricated via the non-solvent induced phase separation (NIPS) technique. Four NIPS fabricated polymer membranes were evaluated: polysulfone, cellulose acetate, polyvinylidene fluoride (PVDF) fabricated with polyethylene glycol (PEG) as a pore-forming agent, and polybenzimidazole (PBI). The HSP approach successfully predicted membrane-solvent compatibility, with polysulfone (RED = 0.6) and cellulose acetate (RED = 0.9) dissolving completely within 24 h, while PVDF (RED = 1.9) and PBI (RED = 1.1) maintained structural integrity throughout a 7-day exposure period. Furthermore, PVDF demonstrated superior performance with minimal weight gain (3.0%), hydrophobic surface characteristics (122° water contact angle), and enhanced mechanical properties following DES exposure. Comprehensive chemical and morphological characterization confirmed PVDF’s chemical stability and revealed a surface-selective interaction mechanism involving simultaneous PEG (pore-forming agent) extraction and DES component adsorption. Adsorption kinetics followed pseudo-first-order behavior with reversible characteristics, best described by the Temkin isotherm model (R² = 0.9987). PVDF membranes-maintained separation functionality with average lignin rejection (75.2 ± 7.69%) and demonstrated filtration permeability of 2.0 ± 0.34 LMH/bar. This methodology provides a rational approach for membrane selection in emerging solvent systems, contributing to the advancement of sustainable separation technologies for DES-based biomass processing applications.

Science & Technology - Other Topics↗

NASA NTRS1992

Flight Testing

It is considered by the Flight Mechanics Panel of AGARD as vitally important that the NATO flight test community meet regularly so that new techniques for flight test, instrumentation and data analysis and lessons learned from past and on-going programs be disseminated to ensure that safe efficient cost-effective and timely testing is accomplished. There are many new systems being tested or planned for testing in the near future. These include programmable signal processor radars, integrated flight, fire and propulsion control systems, thrust vectoring, low observable technologies, multifunction pilot displays and multisensor integration. Acquisition and processing of large quantities of avionics multiplex data are challenges that must be met. There is a need for greater use of simulators and other hardware-in-the-loop ground test facilities.

Source record↗

NASA NTRS2026

Highly Porous Polyimide Gel for Use as Battery Separator with Room Temperature Ionic Liquid Electrolytes

Advanced aerospace vehicular concepts require advances in many existing technologies, including space power and energy storage systems. Batteries represent one of the major areas in need of improvement, both in terms of energy density and safety, with growing concerns over the fire safety of commercial lithium-ion batteries. This has prompted efforts to develop nonflammable battery components, namely the electrolyte and separator. Existing commercial lithium-ion batteries utilize polyolefin microporous membranes as separators with an electrolyte consisting of a lithium salt dissolved in a mixture of cyclic carbonate solvents. This separator/electrolyte combination has ionic conductivities in the range of 10 −2 to 10 −3 S/cm. However, the cyclic carbonate solvents are inherently flammable. Room-temperature ionic liquids (RTILs) appear to be a safer alternative. They offer good ionic conductivities and are inherently nonvolatile and nonflammable, giving them a safety advantage. However, many promising RTILs for battery electrolytes are not compatible with commercial polyolefin separator materials. Alternative separator materials, such as polyimides, are non-flammable and are capable of accepting RTILs into their structure. Polyimide gels, with a composition of 4,4′-oxydianiline, 3,3′,4,4′-tetracarboxylic dianhydride and cross-linked with Desmodur N3300A, possess an open-porous, fibrillar network architecture which offers a high degree of porosity (typically greater than 85% porosity) for lithium-ion transport and conduction, as well as good mechanical properties. Furthermore, these polyimide gels are compatible with selected imidazolium-based RTILs. Nonflammable separator/electrolyte systems with room-temperature conductivities in the range of 10 −3 S/cm have been evaluated. It has been demonstrated that 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide was the most promising among six RTILs screened, in terms of both ionic conductivity and constant current cycling.

Polyimide↗

NASA NTRSDate not supplied

Evaluation of Low-Energy Hydrogen Separation Method Using Metal-Organic Frameworks (MOFs) for Closed-Loop ECLSS Air Revitalization (CLEAR)

The State-of-the-Art (SOA) air revitalization architecture onboard the International Space Station (ISS) recovered approximately 50% of the oxygen (O 2 ) from metabolic carbon dioxide (CO 2 ) via the Sabatier process from 2011 to 2017. O 2 recovery is currently constrained by the limited availability of reactant hydrogen (H 2 ) preventing complete conversion of CO 2 to H 2 O. Increasing O 2 recovery within Closed-Loop ECLSS is essential to reducing resupply mass for long-duration manned missions; specifically focusing on water (H 2 O) which supplies H 2 for Sabatier via water electrolysis. Past ground test endeavors at Marshall Space Flight Center (MSFC) have attempted to recover H 2 from Sabatier-produced CH 4 using technologies such as carbon vapor deposition (CVD) and plasma pyrolysis. The byproducts of these technologies can act as a catalyst poison or reactor deadload to the Sabatier reactor. Hydrogen separation techniques must be utilized to maintain the Sabatier catalyst during gas recycling and must be scalable, non-energy intensive, and safe to operate in a habitation setting. Research indicated that metal-organic frameworks (MOFs) could meet these criteria and were tested for their capability to capture the various carbon-based gaseous products of CVD and plasma pyrolysis such as acetylene (C 2 H 2 ), ethylene (C 2 H 4 ), ethane (C 2 H 6 ), and carbon monoxide (CO) which would purify the hydrogen gas stream passing through the MOF. A sub-scale adsorption column was developed by Marshall Space Flight Center to test three MOF candidates against a synthetic gas mixture comprised of process-relevant carbonous gases and hydrogen to evaluate the separation capability of the MOFs. The results of the hydrogen separation capability, isothermal desorption capability, and demonstrated cyclic reuse of the MOF are presented in this paper.

Kagen Crawford↗

DOE OSTI2023

Temperature Dependence of Low‐Frequency Noise Characteristics of NiO x /β‐Ga 2 O 3 p–n Heterojunction Diodes

Temperature dependence of the low-frequency electronic noise in NiO x /β-Ga 2 O 3 p–n heterojunction diodes is reported. The noise spectral density is of the 1/f-type near room temperature but shows signatures of Lorentzian components at elevated temperatures and at higher current levels (f is the frequency). It is observed that there is an intriguing non-monotonic dependence of the noise on temperature near T = 380 K. The Raman spectroscopy of the device structure suggests material changes, which results in reduced noise above this temperature. The normalized noise spectral density in such diodes is determined to be on the order of 10 −14 cm 2 Hz −1 (f = 10 Hz) at 0.1 A cm −2 current density. In terms of the noise level, NiO x /β-Ga 2 O 3 p–n diodes perform excellently for new technology and occupy an intermediate position among devices of various designs implemented with different ultra-wide-bandgap semiconductors. The obtained results are important for understanding the electronic properties of NiO x /β-Ga 2 O 3 heterojunctions and contribute to the development of noise spectroscopy as the quality assessment tool for new electronic materials and device technologies.

1/f noise↗

NASA NTRS2004

Energetics of Single Substitutional Impurities in NiTi

Shape-memory alloys are of considerable current interest, with applications ranging from stents to Mars rover components. In this work, we present results on the energetics of single substitutional impurities in B2 NiTi. Specifically, energies of Pd, Pt, Zr and Hf impurities at both Ni and Ti sites are computed. All energies are computed using the CASTEP ab initio code, and, for comparison, using the quantum approximate energy method of Bozzolo, Ferrante and Smith. Atomistic relaxation in the vicinity of the impurities is investigated via quantum approximate Monte Carlo simulation, and in cases where the relaxation is found to be important, the resulting relaxations are applied to the ab initio calculations. We compare our results with available experimental work.

Binary Alloys↗