Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT calculation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nuclear–electronic orbital approach to quantization of protons in periodic electronic structure calculations

The nuclear–electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born–Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems, such as two-dimensional materials and liquid–semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide–water interface. The zero-point energy effects of the protons as well as electron–proton correlation are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

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

S K-edge XAS of Cu Ⅱ , Cu Ⅰ , and Zn Ⅱ oxidized Dithiolene complexes: Covalent contributions to structure and the Jahn-Teller effect

We report reduced dithiolene ligands are bound to high valent Mo centers in the active site of the oxotransferase family of enzymes. Related model complexes have been studied with great insight by Prof. Holm and his colleagues. This study focuses on the other limit of dithiolene chemistry: an investigation of the 2-electron oxidized dithiolene bound to low-valent late transition metal (TM) ions (Zn II , Cu I , and Cu II ). The bonding descriptions of the oxidized dithiolene [N,N-dimethyl piperazine 2,3-dithione (Me 2 Dt 0 )] complexes are probed using S K-edge X-ray absorption spectroscopy (XAS) and the results are correlated to density functional theory (DFT) calculations. These experimentally supported calculations are then extended to explain the different geometric structures of the three complexes. The Zn II (Me 2 Dt 0 ) 2 complex has only ligand-ligand repulsion so it is stabilized at the D 2d symmetry limit. The Cu I (Me 2 Dt 0 ) 2 complex has additional weak backbonding thus distorts somewhat from D 2d toward D 2h symmetry. The Cu II (Me 2 Dt 0 ) 2 complex has a strong σ donor bond that leads to both a large Jahn-Teller stabilization to D 2h and an additional covalent contribution to the geometry. The combined strong stabilization results in the square planar, D 2h structure. This study quantifies the competition between the ligand-ligand repulsion and the change in electronic structures in determining the final geometric structures of the oxidized dithiolene complexes, and provides quantitative insights into the Jahn-Teller stabilization energy and its origin.

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

Novel ceramic capacitors with ultrahigh energy density and efficiency (Final Technical Report)

Antiferroelectric ceramics are a special class of material that have shown great potential as the dielectric in electrical capacitors due to their high energy- and power-density. During each charge-discharge cycle, the ceramic undergoes transformation to a ferroelectric phase and resumes its antiferroelectric phase. The hysteresis associated with the transitions leads to a mediocre energy efficiency and service lifetime of antiferroelectric capacitors and, hence, their almost absence in commercial products. Under the support of this research project, we first formulated a universal lattice-compatibility theory that included electrostatic polarization energy along with elastic energy and thermal energy to understand the origin of the hysteresis in antiferroelectric oxides. Guided by this compatibility theory, we conducted high-throughput density functional theory (DFT) calculations to assess chemical modifiers and their effect on crystal structures of 400+ PbZrO 3 -based compositions. Down-selected compositions were experimentally validated for their suppressed hysteresis and higher energy efficiency. The verified low-hysteresis compositions were then expanded to an antiferroelectric ceramic library with nearly 500 new compositions (more than 1,500 samples) using high-throughput experiments involving ceramic synthesis and property screening. The large quantity of data generated (theory and experimental) in these tasks were processed by machine-learning techniques and identified trends were fed to the next iteration. In the end, we successfully discovered four compositions with near-zero hysteresis, yielding a world-record energy efficiency of 98.2% at an energy density of 3.0 J/cm 3 . Furthermore, our antiferroelectric ceramic capacitor reaches 79.5 million charge-discharge cycles lifetime, a factor of 80 enhancement over previous antiferroelectric ceramics with large hysteresis. These research accomplishments have not only met the milestones set in the SOPO, but also led to two patent filings, three journal publications (one of them was in Advanced Materials, impact factor 29.4), and nine oral presentations at various venues. Through the course of the project, three postdocs, four Ph.D. students, and one M.S. student were trained. In short, our project established a new methodology in searching next-generation functional ceramics on the fundamental side and discovered several high-efficiency antiferroelectric compositions for capacitors on the applied side. Once fabricated into the multilayer form for commercial applications, these ceramic capacitors can potentially enable the high temperature high power density DC-link capacitors that are critical for the next generation inverters in electric vehicles. The project also significantly contributed to the nation’s workforce development in the STEM fields.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Discrete Fourier Transform Analysis in a Complex Vector Space

Alternative computational strategies for the Discrete Fourier Transform (DFT) have been developed using analysis of geometric manifolds. This approach provides a general framework for performing DFT calculations, and suggests a more efficient implementation of the DFT for applications using iterative transform methods, particularly phase retrieval. The DFT can thus be implemented using fewer operations when compared to the usual DFT counterpart. The software decreases the run time of the DFT in certain applications such as phase retrieval that iteratively call the DFT function. The algorithm exploits a special computational approach based on analysis of the DFT as a transformation in a complex vector space. As such, this approach has the potential to realize a DFT computation that approaches N operations versus Nlog(N) operations for the equivalent Fast Fourier Transform (FFT) calculation.

Dean, Bruce H.↗
DOE OSTI2023

Validation of the Cossee–Arlman mechanism for propylene oligomerization on Ni/UiO-66

Steady state rate expressions can be derived to distinguish the Cossee–Arlman and metallacycle mechanisms postulated for propylene oligomerization on nickel-based catalysts based on product selectivities, where product selectivities for the former are a function of olefin pressure because sequential coordination and insertion steps lead to independent mechanistic pathways for different hexene isomers. In contrast, the metallacycle mechanism presents pressure-independent product selectivities due to successive coordination prior to the kinetically relevant steps in each mechanism. In this work, steady state propylene oligomerization rates and selectivities were measured in the absence of an activator on nickel functionalized UiO-66 metal organic framework (MOF), Ni/UiO-66, to validate the Cossee–Arlman mechanism for light olefin oligomerization. In situ NO titrations reveal that ~5% of nickel sites were active during the reaction, and thus, not all nickel sites are relevant for catalysis. Propylene dimerization was first order in propylene pressure from 5 to 500 kPa with an apparent activation energy of ~20 kJ mol -1 from 453 to 493 K. Calculated apparent activation energies with density functional theory (DFT) calculations on cluster models of Ni/UiO-66 are in agreement with experiment to corroborate the Cossee–Arlman mechanism. Selectivities of hexene products and the ratio of hexene product selectivities on Ni/UiO-66 are in accordance with selectivity expressions derived from the Cossee–Arlman mechanism. In conclusion, analysis of product selectivities can be used more extensively to demarcate the Cossee–Arlman and metallacycle mechanisms for olefin oligomerization on metal-based catalysts.

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

Scale-invariant machine-learning model accelerates the discovery of quaternary chalcogenides with ultralow lattice thermal conductivity

We design an advanced machine-learning (ML) model based on crystal graph convolutional neural network that is insensitive to volumes (i.e., scale) of the input crystal structures to discover novel quaternary chalcogenides, AMM'Q 3 (A/M/M' = alkali, alkaline earth, post-transition metals, lanthanides, and Q = chalcogens). These compounds are shown to possess ultralow lattice thermal conductivity (κ l ), a desired requirement for thermal-barrier coatings and thermoelectrics. Upon screening the thermodynamic stability of ~1 million compounds using the ML model iteratively and performing density-functional theory (DFT) calculations for a small fraction of compounds, we discover 99 compounds that are validated to be stable in DFT. Taking several DFT-stable compounds, we calculate their κ l using Peierls–Boltzmann transport equation, which reveals ultralow κ l (<2 Wm -1 K -1 at room temperature) due to their soft elasticity and strong phonon anharmonicity. Our work demonstrates the high efficiency of scale-invariant ML model in predicting novel compounds and presents experimental-research opportunities with these new compounds.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Association Mechanisms of Unsaturated C2 Hydrocarbons with Their Cations: Acetylene and Ethylene

The ion-molecule association mechanism of acetylene and ethylene with their cations is investigated by ab initio quantum chemical methods to understand the structures, association energies, and the vibrational and electronic spectra of the products. Stable puckered cyclic isomers are found as the result of first forming less stable linear and bridge isomers. The puckered cyclic complexes are calculated to be strongly bound, by 87, 35 and 56 kcal/mol for acetylene-acetylene cation, ethylene-ethylene cation and acetylene-ethylene cation, respectively. These stable complexes may be intermediates that participate in further association reactions. There are no association barriers, and no significant inter-conversion barriers, so the initial linear and bridge encounter complexes are unlikely to be observable. However, the energy gap between the bridged and cyclic puckered isomers greatly differs from complex to complex: it is 44 kcal/mol in C4H4 +, but only 6 kcal/mol in C4H8 +. The accurate CCSD(T) calculations summarized above are also compared against less computationally expensive MP2 and density functional theory (DFT) calculations for structures, relative energies, and vibrational spectra. Calculated vibrational spectra are compared against available experiments for cyclobutadiene cation. Electronic spectra are also calculated using time-dependent DFT.

ethylene↗
DOE OSTI2022

Origin of ultralow phonon transport and strong anharmonicity in lead-free halide perovskites

All-inorganic lead-free halide double perovskites offer a promising avenue toward non-toxic, stable optoelectronic materials, properties that are missing in their prominent lead-containing counterparts. Their large thermopowers and high carrier mobilities also make them promising for thermoelectric applications. Here, in this work, we present a first-principles study of the lattice vibrations and thermal transport behaviors of Cs 2 SnI 6 and γ-CsSnI 3 , two prototypical compounds in this materials class. We show that conventional static zero temperature density functional theory (DFT) calculations severely underestimate the lattice thermal conductivities (κ l ) of these compounds, indicating the importance of dynamical effects. By calculating anharmonic renormalized phonon dispersions, we show that some optic phonons significantly harden with increasing temperature (T), which reduces the scattering of heat carrying phonons and enhances calculated κ l values when compared with standard zero temperature DFT. Furthermore, we demonstrate that coherence contributions to κ l , arising from wave like phonon tunneling, are important in both compounds. Overall, calculated κ l with temperature-dependent interatomic force constants, built from particle and coherence contributions, are in good agreement with available measured data, for both magnitude and temperature dependence. Large anharmonicity combined with low phonon group velocities yield ultralow values, with room temperature values of 0.26 W/m-K and 0.72 W/m-K predicted for Cs 2 SnI 6 and γ-CsSnI 3 , respectively. We further show that the lattice dynamics of these compounds are highly anharmonic, largely mediated by rotation of the SnI6 octahedra and localized modes originating from Cs rattling motion. These thermal characteristics combined with their previously computed excellent electronic properties make these perovskites promising candidates for optoelectronic and room temperature thermoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine-Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

To cope with legacy greenhouse gas emissions and to achieve net-zero emissions by 2050, the U.S. Department of Energy (DOE) is funding efforts to develop direct air capture (DAC), a method for removing CO2 directly from air. Metal organic frameworks (MOFs) have been well studied as DAC sorbent materials due to their tunable structural and compositional properties. Thermodynamic simulations using force fields are often used to provide predictions of a material’s performance in many separations. However, these force fields often make assumptions about bonds and the physics of the adsorption process. A new class of force fields called machine-learned force fields (MLFFs) use machine learning to form quantitative relationships between a material’s chemical structure and the forces and energies predicted by more accurate quantum mechanical calculations, such as dispersion-corrected density functional theory (DFT). In this work, models were developed to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption using MLFFs. These methods were parametrized based on thousands of DFT calculations of CO2 in flexible MOFs and used to predict MOF structural properties as well as CO2 adsorption in several MOFs.

Findley, John↗
DOE OSTI2023

Interpretable, Extensible Linear Regression Model for Electron Density Prediction

This code allows us to generate a model to predict ground state electron densities directly from atomic structure thereby bypassing the need to perform direct DFT calculations. It is based on many body correlation descriptors defined in https://arxiv.org/abs/2004.14442 (IM release numbers: LLNL-JRNL- 808791 and LLNL-JRNL-840698). The code reads in electron densities, in CHGCAR format used by the DFT code VASP, calculates many-body correlation descriptors based on atoms surrounding a grid point and then uses linear regression to obtain the fitting parameters.

Kumar, Shashikant↗
DOE OSTI2026

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2021

Elimination of the linearization error in APW/LAPW basis set: Dirac-Kohn-Sham equations

Here, a detailed account of the implementation of equations of the relativistic density functional theory (RDFT) using basis sets of APW/LAPW type with flexible extensions provided by local orbitals is given. Earlier discoveries of the importance of the high derivative local orbital (HDLO) extension of the APW/LAPW basis set for enhancing the accuracy of DFT calculations are confirmed using a fully relativistic approach and α – U as an example. High energy local orbitals (HELO's), however indispensable for GW calculations, are considerably less efficient in enhancing the accuracy of DFT applications. It is shown that a simplified approach to the relativistic effects, namely considering them only inside the muffin-tin (MT) spheres, produces basically identical results (as compared to fully relativistic approach) for the electronic free energy of the five materials considered in this work. By comparing the effect of the simplified approach on the electronic free energy with its effect on the electronic kinetic energy we conclude that the insensitivity of the free energy to the way we describe the relativistic effects in the interstitial region is related to the variational property of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Surfactant inhibition mechanisms of carbonate mineral dissolution in shale

Surfactants are common additives to hydraulic fracturing and enhanced oil recovery (EOR) fluids, and are under consideration for amendment to supercritical carbon dioxide for geological carbon sequestration (GCS). The effect of a common anionic surfactant, internal olefin sulfonate (IOS), on mineral dissolution from shale into brine was evaluated. When added to brine at concentrations exceeding the critical micelle concentration (94 mg/L), IOS inhibited carbonate mineral dissolution in an Eagle Ford shale, as well as dissolution of optical quality calcite (the dominant carbonate in the shale). Laser profilometry images provide spatial resolution across > 3 orders of magnitude, and indicate that IOS addition to brine both enhances the formation of new etch pits in calcite, and impedes their further growth. Time-of-flight secondary ion mass spectrometry surface profiles show for the first time that IOS preferentially adsorbs at calcite pit edges versus flat calcite surfaces (i.e., terraces). Surface pressure calculations, sulfur K-edge near edge X-ray absorption fine structure (NEXAFS) spectroscopy results, and density functional theory (DFT) calculations support this observation; the DFT results indicate that the sulfonate head group of the IOS molecule binds strongly to the calcite step site as compared to the terrace site. The S K-edge NEXAFS results indicate that IOS adsorbed more to etched calcite surfaces compared to smooth calcite surfaces. Overall, the results indicate that weak adsorption on flat calcite surfaces (i.e., terraces) disrupts water structure and enhances mass transfer of dissolution, while strong adsorption on calcite pit edges displaces adsorbed water and inhibits further etch pit growth. This work provides the first direct evidence of preferential adsorption of IOS to etched calcite surfaces and links it to macroscopic dissolution kinetics. Finally, this work has implications for surfactant-containing fluids used in hydraulic fracturing, EOR and potentially GCS for subsurface injection into carbonate rich reservoirs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unexpected features in the optical vibrational spectra of δ-UO 3

Uranium trioxide displays a complex chemical phase space, with at least six structurally distinct polymorphs accessible via different synthetic routes. Remarkably, despite its technological importance, full structural and electronic characterization of these polymorphs remains an open area of study. δ-UO 3 in particular has attracted significant theoretical attention due to its high point group and space group symmetries, having U (VI) in octahedral coordination with polyhedra interconnected through corner-sharing to build a 3-D cubic lattice with space group symmetry Pm-3m and Z = 1. Critical experimental information, such as its optical vibrational spectra, are not known. Here, we study the Raman and infrared (IR) spectra of δ-UO 3 together with the support of density functional theory (DFT) calculations for spectral interpretation. A symmetry analysis of the DFT-predicted phonon eigenmodes indicates that δ-UO 3 should have two IR active modes and no Raman active modes. Experimental results, however, indicate significant Raman scattering from δ-UO 3 . We therefore propose four potential explanations for this apparent contradiction: a possible tetragonal distortion to the cubic cell, the existence of a surface impurity layer, vacancy scattering, and structural activation of Raman signal. We use powder X-ray diffraction and confocal Raman spectroscopy with depth profiling to investigate these possibilities and suggest future experiments to explore this phenomenon in more detail. Understanding the lattice dynamics of δ-UO 3 is important for identification of technogenic U phases via Raman and infrared spectroscopy and our results indicate that the simple understanding of δ-UO 3 as a high-symmetry cubic structure should be reconsidered.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Single crystal structure, vibrational spectroscopy, gas sorption and antimicrobial properties of a new inorganic acidic diphosphates material (NH 4 ) 2 Mg(H 2 P 2 O 7 ) 2 • 2 H 2 O

We report on the successful synthesis of diammonium magnesium dihydrogendiphosphate (V) dihydrate compound (NH 4 ) 2 Mg(H 2 P 2 O 7 ) 2 •2H 2 O using a wet chemical route. Single crystal X-ray diffraction analysis and micro Raman spectroscopy are employed to characterize the compound. We demonstrate, using a multidisciplinary approach, that this compound is ideal for carbon dioxide (CO 2 ) capture in addition to other anthropogenic gasses. We show here -from both an experimental as well as from a density functional theory (DFT) calculations routes- the potential for adopting this compound into domestic air-conditioning units (ACUs). From these experiments, the resistance to bacterial growth is also investigated, which is critical for the adoption of this compound in ACUs. Our compound exhibits a higher methane (CH 4 ) sorptivity as compared to CO 2 at 25 °C and 45 °C under pressures up to 50 bars. Furthermore, DFT electronic structure calculations are used to compute the main structural and electronic properties of the compound, taking into consideration the characteristics of the identified pores as a function of the progressive CO 2 vs. CH 4 loadings. Finally, the antibacterial assay reveals a strong antibacterial activity against the tested Gram-positive and Gram-negative bacteria, with a large zone of inhibition against the tested E. Coli, S. Aureus and K. Pneumonia

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

Rhenium-Selenido Corroles: Reflections on 5d Metalloporphyrins and Metallocorroles as Triplet Emitters and Photosensitizers

Following the successful synthesis of rhenium-oxido, rhenium-imido, and rhenium-sulfido corroles, a series of para-X-substituted rhenium-selenido triarylcorroles (X = OCH 3 , CH 3 , H, F, CF 3 ) have been prepared in 70–76% yields. A one-pot, two-step procedure was used, analogous to that used for ReS corroles, involving rhenium insertion via high-temperature reaction of Re 2 (CO) 10 and free-base corroles in refluxing decalin, followed by exposure to PCl 3 (which is thought to deoxygenate ReO corroles formed under the reaction) and elemental selenium. An analogous procedure, however, failed to yield rhenium-tellurido corroles, presumably reflecting, according to DFT calculations, the limited stability of these species. Unlike ReO corroles, ReS and ReSe corroles were found not to exhibit phosphorescence in the NIR-I regime (600–1000 nm); nor did they sensitize singlet oxygen formation. In the hope of obtaining a broader perspective of this negative result, we also examined ReN porphyrins and found them not to phosphoresce or to sensitize singlet oxygen formation. These results were explained by DFT and TDDFT calculations in terms of low-energy triplet states, which are not energetic enough to excite molecular oxygen to its lowest singlet state. Whether some of the new Re corroles exhibit phosphorescence in the NIR-II regime remains an interesting question for future investigation.

Aromatic compounds↗
DOE OSTI2023

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simple self-consistent method for excited states in density functional theory to characterize defect-derived behavior in wide-band-gap-based microelectronic materials

This final report summarizes the results of the Laboratory Direct Research and Development (LDRD) Project Number 229740. Wide band gap semiconductors such as gallium nitride (GaN) have features highly desirable for multiple mission electronic applications. Realization of their potential requires atomic-scale understanding of electronic behavior. The principal experimental tools for electronically probing defects in GaN are chemically undifferentiating and lack a practical theoretical counterpart needed to identify and characterize specific defects. This project investigated whether a simple idea for modeling defect excited states and their associated photoluminescence (PL) energies is viable, as a path to accelerate the understanding of defect behavior and gain valuable insights into engineering new electronic materials and devices. The research implemented a non-self-consistent total-energy evaluation of a Koopmans-type estimation of an excited electronic state energy in density functional theory (DFT) calculations, and proceeded to design, implement, and assess a self-consistent method for computing excited states based upon an OCcupation-Constrained-DFT (occ-DFT). The occ-DFT was verified in test calculations of defect excited states and validated against well-characterized PL data for 3d transition metal defects in GaN. The method proved stable and robust in computing excited states and gave accurate predictions compared to experimental PL data. The combined ground state/excited-state capability proved capable of chemically differentiating defect species in GaN. In application to 3d dopants in GaN, we reinterpreted extensive experimental literature, proposed new defects as prospective candidates for use in quantum information applications, and outlined design strategies to create and exploit these potentially useful functional defects in GaN.

36 MATERIALS SCIENCE↗