Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oxygen vacancy”

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 289 records · Page 16

Effects of local compositional and structural disorder on vacancy formation in entropy-stabilized oxides from first-principles

Entropic stabilization has evolved into a strategy to create new oxide materials and realize novel functional properties engineered through the alloy composition. Achieving an atomistic understanding of these properties to enable their design, however, has been challenging due to the local compositional and structural disorder that underlies their fundamental structure-property relationships. Here, we combine high-throughput atomistic calculations and linear regression algorithms to investigate the role of local configurational and structural disorder on the thermodynamics of vacancy formation in (MgCoNiCuZn)O-based entropy-stabilized oxides (ESOs) and their influence on the electrical properties. We find that the cation-vacancy formation energies decrease with increasing local tensile strain caused by the deviation of the bond lengths in ESOs from the equilibrium bond length in the binary oxides. The oxygen-vacancy formation strongly depends on structural distortions associated with the local configuration of chemical species. Vacancies in ESOs exhibit deep thermodynamic transition levels that inhibit electrical conduction. By applying the charge-neutrality condition, we determine that the equilibrium concentrations of both oxygen and cation vacancies increase with increasing Cu mole fraction. Our results demonstrate that tuning the local chemistry and associated structural distortions by varying alloy composition acts an engineering principle that enables controlled defect formation in multi-component alloys.

36 MATERIALS SCIENCE↗

Investigation and Suppression of Oxygen Release by LiNi 0.8 Co 0.1 Mn 0.1 O 2 Cathode under Overcharge Conditions

The safety issue of lithium-ion batteries is a crucial factor limiting their large-scale application. Therefore, it is of practical significance to evaluate the impact of their overcharge behavior because of the severe levels of oxygen release of cathode materials during this process. Here, by combining a variety of in situ techniques of spectroscopy and electron microscopy, this work studies the structural degradation of LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM811) accompanying the oxygen release in the overcharge process. It is observed that a small amount of O 2 evolves from the initial surface at ≈4.7V. When charging to a higher voltage (≈5.5V), a large amount of O 2 evolves on the newly formed surface due to the occurrence of microcracks. Based on experimental results and theoretical calculations, it is determined that the oxygen release mainly occurs in the near-surface regions, where the remaining oxygen vacancies accumulate to create voids. To suppress the oxygen release, single-crystalline NCM811 with integrated structure is introduced and serves as a cathode, which can effectively inhibit morphology destruction and reduce the activation of lattice oxygen in the surface region. These findings provide a theoretical basis and effective strategy for improving the safety performance of Ni-rich cathode materials in practical applications.

25 ENERGY STORAGE↗

Defect Thermodynamics and Transport Properties of Proton Conducting Perovskite Electrode and Electrolyte Materials Evaluated Based on Density Functional Theory Modeling

Both electron-rich and electron-poor perovskite oxides have been used in solid oxide cell applications as electrode and electrolyte materials. The rich oxygen defect chemistry and its coupling to temperature, hydrogen-steam or oxygen-steam gas pressure, or to the applied potentials creates enormous complexities for modeling performance and degradation of the materials. Herein, density functional theory-based thermodynamic modeling was carried out to describe the defect chemistry and transport properties of the proton-conducting electrolyte BaZr1-xYxO3-δ (x≤0.1) and of the triple-conducting perovskite (La,Ba)(Fe,M)O3-δ (M=Y and Zr). The defect thermodynamics of intrinsic point defects and the hydrogen-related defect reactions were solved in integrated defect models and further used to predict the Brouwer diagram and the transport properties of the functional perovskites. For the electron-poor electrolytes BaZr0.9Y0.1O3-δ, the developed model has been used to describe the experimental transport properties in the SOC operating conditions. Specifically, the roles played by the acceptor-bound holes and the intrinsic and hydrogen point defects upon the conductivities of holes, protons, and oxygen vacancies under the hydrogen-rich and oxygen-rich conditions at various humidity levels were demonstrated. A defect modeling tool was also developed for the triple-conducting perovskite (La,Ba)(Fe,M)O3-δ (M=Y and Zr) to examine magnetic effects and hydride defects in defect equilibria.

defect thermodynamics↗

Redox-Based Chemical Looping Large-Scale Air Separation Unit Designs Using Perovskite Material

Oxygen production by means of air separation through redox cycle reactions of Sr 1–x Ca x FeO 3-δ perovskite particles is examined numerically in the context of modular large-scale units. Fixed-bed, multicomponent reactor designs allowing for the recovery of the energy released by the exothermic oxidation (or adsorption) to enhance the endothermic reduction (or desorption) of the perovskite material are considered and compared against a baseline cylindrical packed bed. Numerically, the gas–solid reacting system is approximated as a single-phase reacting gas flow through a porous medium of constant porosity. The redox kinetics account for the oxygen nonstoichiometry of the perovskite, which exhibits a dependence on the temperature and the oxygen partial pressure, in order to describe some mechanisms of the oxygen vacancies. Results show that the spontaneous release of the oxygen from the lattice arrangements, through the thermal management of the heat of the oxidation reaction from the geometrical design of the reactors, can be enhanced up to 4% with respect to the baseline oxygen production of 0.201 g O 2 ·h –1 ·g perovskite –1 with the cylindrical packed bed, in the condition of operations in this work, at 500 °C. Sensitivity to the operation temperature showed that the oxygen production is greatly enhanced by approximately 30% with the increase of the temperature from 500 to 550 °C. We report this result is consistent with experimental observations of faster release of the oxygen from the lattice arrangements of this perovskite at 550 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The role of defect charge, crystal chemistry, and crystal structure on positron lifetimes of vacancies in oxides

Density functional theory based positron lifetime (PL) calculations for cation and oxygen monovacancies in a range of oxides—hematite, magnetite, hercynite, and alumina—have been conducted to compare the impact of defect chemistry and crystal structure on the predicted lifetimes. The role of defect charge state has also been examined. A comparison across the same type of crystalline structure but different composition shows that oxygen vacancies only induce a slight increase in the positron-electron overlap and thus barely modify the PL as compared to the bulk. A much more substantial increase of PL is observed for cation monovacancies, regardless of crystal structure or the elemental nature of the vacancy, which we ascribe to an enhanced localization of charge density around the vacant site. The structural and compositional richness of the oxide leads to longer defect PLs, with defected hercynite exhibiting the longest PLs. The charge state of cation monovacancies modifies only by a small percentage the positron localization, relegating to secondary importance the metal defect's oxidation state in modifying the lifetime of positrons within vacancy traps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Tailored (La 0.2 Pr 0.2 Nd 0.2 Tb 0.2 Dy 0.2 ) 2 Ce 2 O 7 as a Highly Active and Stable Nanocatalyst for the Oxygen Evolution Reaction

Abstract Designing highly active and robust catalysts for the oxygen evolution reaction is key to improving the overall efficiency of the water splitting reaction. It has been previously demonstrated that evaporation induced self‐assembly (EISA) can be used to synthesize highly porous and high surface area cerate‐based fluorite nanocatalysts, and that substitution of Ce with 50% rare earth (RE) cations significantly improves electrocatalyst activity. Herein, the defect structure of the best performing nanocatalyst in the series are further explored, Nd 2 Ce 2 O 7 , with a combination of neutron diffraction and neutron pair distribution function analysis. It is found that Nd 3 + cation substitution for Ce in the CeO 2 fluorite lattice introduces higher levels of oxygen Frenkel defects and induces a partially reduced RE 1.5 Ce 1.5 O 5 +x phase with oxygen vacancy ordering. Significantly, it is demonstrated that the concentration of oxygen Frenkel defects and improved electrocatalytic activity can be further enhanced by increasing the compositional complexity (number of RE cations involved) in the substitution. The resulting novel compositionally‐complex fluorite– (La 0.2 Pr 0.2 Nd 0.2 Tb 0.2 Dy 0.2 ) 2 Ce 2 O 7 is shown to display a low OER overpotential of 210 mV at a current density of 10 mAcm −2 in 1M KOH, and excellent cycling stability. It is suggested that increasing the compositional complexity of fluorite nanocatalysts expands the ability to tailor catalyst design.

Chemistry↗

Hydrogen defects in LaBi 2 O 4 X (X = Cl, Br, and I) Sillén oxyhalide phases and their impacts on ionic transport

Sillén oxyhalides have recently emerged as promising materials for both photocatalytic and ionic transport applications, yet the role of likely-ubiquitous hydrogen-related defects in these layered compounds remains largely unexplored. Here, we employ first-principles defect calculations to investigate incorporation energetics for hydrogen- and oxygen-related defects, as well as their migration barriers in LaBi 2 O 4 X (X = Cl, Br, I) phases. We find that hydrogen interstitials, particularly protonic species (H i + ), are readily accommodated within the open Bi–O layers. Protons compete with oxygen vacancy donors (V O 2+ ) and charge-compensate with oxygen interstitial acceptors (O i 2− ). By linking hydrogen defect formation to water- and oxygen-related redox equilibria, we reveal that V O 2+ facilitates H i + incorporation, while O i 2− promotes interstitial hydroxide formation, establishing a direct connection between proton and oxide-ion transport. Calculated migration barriers indicate that ionic diffusion is confined to Bi–O layers with low barriers of 0.20–0.25 eV for H i + and 0.14–0.25 eV for V O 2+ , suggesting that the materials contain intrinsic pathways for mixed ionic conduction. These results provide a microscopic picture of hydrogen behavior in Sillén oxyhalides and point to design strategies for integrating protonic and oxide-ion transport in layered oxyhalide electrolytes. Band-edge alignment analysis shows that LaBi 2 O 4 I provides the optimal combination of hydrogen solubility, oxygen defect stability, and mixed ionic conductivity, highlighting its potential for low-temperature electrochemical and energy-conversion applications. Overall, this work establishes the defect-driven origin of hydrogen transport in Sillén oxyhalides and expands their applicability beyond photocatalysis to mixed ionic conduction and hydrogen electrochemistry.

Energy - Conversion↗

Identifying native point defect configurations in α-alumina

Intimately intertwined atomic and electronic structures of point defects govern diffusion-limited corrosion and underpin the operation of optoelectronic devices. For some materials, complex energy landscapes containing metastable defect configurations challenge first-principles modeling efforts. Here, we thoroughly reevaluate native point defect geometries for the illustrative case of α-Al 2 O 3 by comparing three methods for sampling candidate geometries in density functional theory calculations: displacing atoms near a naively placed defect, initializing interstitials at high-symmetry points of a Voronoi decomposition, and Bayesian optimization. We find symmetry-breaking distortions for oxygen vacancies in some charge states, and we identify several distinct oxygen split-interstitial geometries that help explain literature discrepancies involving this defect. We also report a surprising and, to our knowledge, previously unknown trigonal geometry favored by aluminum interstitials in some charge states. Importantly, these new configurations may have transformative impacts on our understanding of defect migration pathways in aluminum-oxide scales protecting metal alloys from corrosion. Overall, the Voronoi scheme appears most effective for sampling candidate interstitial sites because it always succeeded in finding the lowest-energy geometry identified in this study, although no approach found every metastable configuration. Finally, we show that the position of defect levels within the band gap can depend strongly on the defect geometry, underscoring the need to conduct careful searches for ground-state geometries in defect calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Atomically dispersed Ag on δ-MnO 2 via cation vacancy trapping for toluene catalytic oxidation

The regulation of metal–support interaction (MSI) and the improvement of noble metal utilization are the frontiers of volatile organic compound (VOC) catalytic oxidation. Here, we report a H 2 O 2 -assisted synthesis process for trapping Ag + onto the engineered Mn vacancies of δ-MnO 2 for the catalytic oxidation of toluene. Due to the strong MSI between Ag and δ-MnO 2 , the monoatomic Ag + anchored on cation vacancies forms longer and weakened Ag–O bonds compared with Ag 2 O, which can extra promote the formation of oxygen vacancies and the activation of adjacent lattice oxygen in Ag–O–Mn units compared with regular nanoparticle counterparts. It is beneficial to the steps of methyl dehydrogenation, demethylation, and benzene ring cleavage. The catalyst yields excellent toluene oxidation activity, enhanced turnover frequency in terms of Ag, and good stability. Finally, we propose in this study a facile design strategy for fabricating stably atomically dispersed metal sites, as well as clarify the structure–activity relationships during toluene catalytic oxidation, especially the important effect of single-atom Ag on the adjacent local structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Surface Hydroxyls and Lattice Oxygen in Governing Selectivity and Restructuring During Acetic Acid Conversion on Fe 3 O 4 (001)

Understanding the reactivity of carboxylic acids on metal oxide surfaces is critical for elucidating ketonization mechanisms relevant to biomass upgrading. Here, we investigate the adsorption and thermal decomposition of acetic acid (CH 3 COOH) on Fe 3 O 4 (001) using scanning tunneling microscopy (STM), temperature-programmed reaction spectroscopy (TPRS), and X-ray photoelectron spectroscopy (XPS). At room temperature, acetic acid adsorbs dissociatively to form ordered bidentate acetate (CH 3 COO) overlayer that lifts the (√2 × √2)R 45° surface reconstruction. TPRS reveals ketene (CH 2 CO) as the dominant product, along with CO, CO 2 , and H 2 O, the latter evolving via a Mars–van Krevelen (MvK) mechanism. Isotopic labeling shows preferential CO 2 formation from the carboxyl carbon and a more balanced CO/CO 2 ratio from the methyl carbon, suggesting distinct oxidation pathways. STM imaging reveals embedded acetate intermediates filling surface oxygen vacancies created in MvK steps. Upon product formation completion (~700 K), extensive surface etching is observed, with pits elongated along the octahedral Fe rows. Approximately 20% of the surface oxygen is removed, consistent with vacancy formation stoichiometry inferred from product distributions. These findings demonstrate that carboxylate-induced restructuring of Fe 3 O 4 (001) involves both surface healing and reduction processes, offering mechanistic insights relevant to ketonization and broader carboxylic acid chemistry on metal oxides.

acetic acid↗

Machine learning informed rational design of high entropy double perovskite oxide universal air/steam electrodes for solid oxide electrochemical cells

Due to their high efficiency and versatility, solid oxide electrochemical cells (SOCs) are poised to play a significant role in future energy conversion and storage applications. In recent years, SOCs have bifurcated into two distinct categories: traditional oxygen-ion conducting SOCs that typically operate from ∼650—850 °C and the more recent proton-conducting ceramic (PCC) SOCs that typically operate from ∼400—650 °C. Current performance and lifetime of both oxygen-ion conducting SOCs and PCCs is primarily limited by the air/steam electrode, which facilitates the oxygen reduction reaction (ORR) during fuel cell operation and must also facilitate the oxygen evolution reaction (OER) during electrolysis operation. Here, we present a newly designed high-entropy double perovskite oxide suitable as a universal ORR/OER electrode for both oxygen-ion conducting SOCs and PCCs. Machine learning methods are applied to identify chemical descriptors for highly catalytic high-entropy double perovskite oxides (AA’B 2 O 6 ) across a large compositional space. Based on the machine-learning guidance, we ultimately converge on Ba 0.9 Cs 0.1 (Ca 0.2 Gd 0.2 La 0.2 Pr 0.2 Sr 0.2 )Co 1.5 Fe 0.5 O 6 (CsBaHEO) as a universal air/steam electrode. Structure stabilization is accomplished by an equimolar five-cation high-entropy composition on the A’-site, while cesium substitution on the A-site enhances the electrical conductivity and leads to a higher oxygen vacancy concentration. This material exhibits versatility and high performance in reversible oxygen-ion SOCs, reversible PCCs, and also large-scale tubular PCCs. For example, the CsBaHEO-based PCC reaches 1018 mW∙cm −2 at 600°C, while a large-scale tubular PCC using CsBaHEO for electrolysis achieves a hydrogen production rate of 21.314 ML∙min −1 at 600 °C.

Cell↗

Machine-Learning-Driven Discovery of Water Splitting BaFe 2 O 4 and Human-in-the-Loop Improvement via Al-Substitution for Increased Thermal Stability

Thermochemical hydrogen (TCH) production offers a promising method for converting thermal energy into hydrogen fuel through heat-driven redox cycles of metal oxides. Here, in this work a defect graph neural network (dGNN) was used to predict oxygen vacancy formation energies ΔH V O combined with Materials Project predictions of oxygen chemical potential stability to screen candidate oxides via high-throughput database analysis. BaFe 2 O 4 was identified as a promising material for experimental validation based on its predicted ΔH V O , oxygen chemical potential stability range, and potential for tunable substitutions to improve thermal properties. Experimental validation using thermogravimetric analysis (TGA), stagnation flow reactor (SFR), X-ray diffraction (XRD), and electron microscopy confirmed positive water-splitting behavior but also revealed limitations in thermal stability under aggressive reduction conditions. To address this, a human-in-the-loop modification strategy was employed introducing Al substitution in BaFe 2–x Al x O 4 ; this modification improves thermal stability, alters the crystal structure and enhances overall performance. These results demonstrate a combined computational and experimental workflow in which machine learning accelerates identification of promising candidates, while targeted experimental design enables optimization of functional performance. This approach advances the development of robust, cost-effective TCH materials and highlights the importance of integrating data-driven discovery with human-guided materials design in paving the way for scalable hydrogen production technologies.

organic↗

Thermal stability and coalescence dynamics of exsolved metal nanoparticles at charged perovskite surfaces

Exsolution reactions enable the synthesis of oxide-supported metal nanoparticles, which are desirable as catalysts in green energy conversion technologies. It is crucial to precisely tailor the nanoparticle characteristics to optimize the catalysts’ functionality, and to maintain the catalytic performance under operation conditions. We use chemical (co)-doping to modify the defect chemistry of exsolution-active perovskite oxides and examine its influence on the mass transfer kinetics of Ni dopants towards the oxide surface and on the subsequent coalescence behavior of the exsolved nanoparticles during a continuous thermal reduction treatment. Nanoparticles that exsolve at the surface of the acceptor-type fast-oxygen-ion-conductor SrTi 0.95 Ni 0.05 O 3–δ (STNi) show a high surface mobility leading to a very low thermal stability compared to nanoparticles that exsolve at the surface of donor-type SrTi 0.9 Nb 0.05 Ni 0.05 O 3–δ (STNNi). Our analysis indicates that the low thermal stability of exsolved nanoparticles at the acceptor-doped perovskite surface is linked to a high oxygen vacancy concentration at the nanoparticle-oxide interface. For catalysts that require fast oxygen exchange kinetics, exsolution synthesis routes in dry hydrogen conditions may hence lead to accelerated degradation, while humid reaction conditions may mitigate this failure mechanism.

25 ENERGY STORAGE↗

Engineering Thermally Resilient and Kinetically Active Reversible Protonic Ceramic Cells via Interfacial Design

Achieving concurrent fast electrode kinetics and long-term thermo-mechanical durability remains a critical challenge for reversible protonic ceramic electrochemical cells (R-PCECs). Herein, we report an interfacial engineering strategy that integrates a perovs.kite-type PrBaRu0.1Co1.9O5+δ (PBRC) nanoparticle layer onto a PrBa0.5Sr0.5Co1.5Fe0.5O5+δ (PBSCF) substrate (PBRC-PBSCF), together with a modified pellet-assisted sintering approach to fabricate dense BaZr0.4Ce0.4Y0.1Yb0.1O3-δ (BZCYYb4411) electrolytes. The in situ reconstructed heterointerface enhances oxygen reduction/evolution reaction (ORR/OER) kinetics, promotes H2O adsorption/dissociation, and improves steam tolerance, as verified by electrochemical measurements and interfacial microstructural analyses. Density functional theory reveals that Ru-induced electronic modulation at the PBRC-PBSCF interface lowers the energy of oxygen vacancy formation and optimizes the position of the O 2p band center, thereby accelerating oxygen redox kinetics and stabilizing the interface. The resulting R-PCECs deliver an excellent peak power density of 1.112 W cm−2 and an electrolysis current density of −1.257 A cm−2 at 1.3 V in 3% H2O wet air at 600°C, with a reasonable faradaic efficiency. Furthermore, the cells demonstrate excellent stability, sustaining 100 h of thermal cycling (400–600°C, 200°C h−1) in both fuel cell and electrolysis modes, with 600 h of stability in electrolysis mode (600°C, −0.5 to −2 A cm−2).

30 DIRECT ENERGY CONVERSION↗

Synergistic catalysis by Mn promoted ceria for molecular oxygen assisted epoxidation

Mesoporous manganese doped cerium oxide catalysts were prepared by an inverse micelle method with different manganese loadings. These materials exhibited superior activity for epoxidation of alkenes with molecular oxygen. Highest active 10MnCe material showed an activity of 80 % conversion and >99 % selectivity for cyclooctene epoxidation. The fluorite crystal structure which facilitates the high oxygen mobility, low temperature reducibility, homogeneous distribution of manganese on ceria, nanoparticle nature, higher number of oxygen vacancies, low valent states of manganese and cerium, and the ability to generate reactive oxygen species of the manganese doped cerium oxide could be correlated to the high catalytic activity. The reaction followed first order kinetics with a rate constant of 0.31 h –1 at 100 °C, and the activation energy turned out to be 10.5 kJ/mol. In conclusion, the mechanistic study revealed that the reaction proceeds via a unique mechanism which involves in situ generated superoxide and singlet oxygen species..

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pretreatment Effects on the Surface Chemistry of Small Oxygenates on Molybdenum Trioxide

Understanding surface reactions of biomass-derived oxygenates on metal oxides is important for designing catalysts for valorization of biomass. As such, This work elucidated the effect of different pretreatments on molybdenum trioxide (MoO 3 ) to understand how surface reactivity is controlled by the surface oxidation state. The catalyst was pretreated in oxidative, inert, and reducing environments. The inert and reducing pretreatments created oxygen vacancies on the catalyst surface that acted as active sites for the adsorption of oxygenated molecules, with the reducing pretreatment yielding a higher density of these active sites. Exposing the catalyst to an alcoholic solvent such as methanol also led to a partial reduction similar to the inert pretreatment. After pretreatment, the catalyst was exposed to ethanol, acetaldehyde, and crotonaldehyde with subsequent characterization by diffuse reflectance infrared spectroscopy (DRIFTS), temperature-programmed desorption (TPD), X-ray absorption near edge spectroscopy (XANES), and X-ray photoelectron spectroscopy (XPS). Density functional theory (DFT) was also used to determine adsorption configurations and energies of ethanol, acetaldehyde, and crotonaldehyde. Reduced surfaces were shown to have a stronger affinity for carbonyls, leading to a higher activity for the aldol condensation of acetaldehyde and ethanol to C 4 molecules. Catalysts pretreated in an oxidative environment were completely inactive toward chemisorption and reaction of acetaldehyde.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-principles study of defects and doping limits in CaO

Calcium oxide (CaO) is a promising host for quantum defects because of its ultrawide bandgap and potential for long spin coherence times. Using hybrid functional calculations, we investigate the intrinsic point defects and how they limit Fermi-level positions and doping in CaO. We find calcium and oxygen vacancies to be the most common intrinsic defects, acting as compensating acceptors and donors, respectively. Oxygen interstitials are also prevailing under O-rich conditions and act as compensating donors. Due to compensation by these defects, O-poor conditions are required to dope CaO n-type, while O-rich conditions are required for p-type doping. We find that, at room temperature, intrinsic CaO can only achieve Fermi-level positions between 1.76 eV above the valence-band maximum (VBM) and 1.73 eV below the conduction-band minimum (CBM). If suitable shallow dopants are found, the allowed range of Fermi levels would increase to between VBM + 0.53 eV and CBM − 0.27 eV and is set by the compensating intrinsic defects. Additionally, we study hydrogen impurities, and show that hydrogen will not only limit p-type doping but can also act as shallow donor when substituting oxygen (HO defects).

Physics↗

Excited State Properties of Point Defects in Semiconductors and Insulators Investigated with Time-Dependent Density Functional Theory

Here, we present a formulation of spin-conserving and spin-flip hybrid time-dependent density functional theory (TDDFT), including the calculation of analytical forces, which allows for efficient calculations of excited state properties of solid-state systems with hundreds to thousands of atoms. We discuss an implementation on both GPU- and CPU-based architectures along with several acceleration techniques. We then apply our formulation to the study of several point defects in semiconductors and insulators, specifically the negatively charged nitrogen-vacancy and neutral silicon-vacancy centers in diamond, the neutral divacancy center in 4H silicon carbide, and the neutral oxygen-vacancy center in magnesium oxide. Our results highlight the importance of taking into account structural relaxations in excited states in order to interpret and predict optical absorption and emission mechanisms in spin defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗