Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electron correlation calculations”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

The bright knots at the tops of soft X-ray flare loops: Quantitative results from Yohkoh

Soft X-ray Telescope (SXT) observations from the Japanese Yohkoh spacecraft have shown that confined bright regions are common features at the tops of flare loops throughout most of the duration of the flares. In this paper we present quantitative results for these flare knots, in relation to other flare regions, for four relatively 'simple' flares. Emission measure distributions, electron temperatures, and electron densities are derived from SXT and Yohkoh Bragg Crystal Spectrometer (BCS) observations. The four flares selected are dominated by what appear to be single-loop structures, with bright knots at the loop tops. The flares are neither long-duration nor impulsive events. The spatial distributions of brightness and emission measure in the flares are found to be quite similar for all four events, even though there are significant differences in dynamical behavior between at least two of the events. Temperatures and densities calculated for these flares are consistent with previous results from many solar experiments. An investigation of intensity correlations between adjacent pixels at the tops of the loops suggests the existence of local disturbances in the magnetic loops that occur on spatial scales less than the radii of the loops.

Doschek, G. A.↗
DOE OSTI2024

Polarization-driven band topology evolution in twisted MoTe 2 and WSe 2

Motivated by recent experimental observations of opposite Chern numbers in R-type twisted MoTe 2 and WSe 2 homobilayers, we perform large-scale density-functional-theory calculations with machine learning force fields to investigate moiré band topology across a range of twist angles in both materials. We find that the Chern numbers of the moiré frontier bands change sign as a function of twist angle, and this change is driven by the competition between moiré ferroelectricity and piezoelectricity. Our large-scale calculations, enabled by machine learning methods, reveal crucial insights into interactions across different scales in twisted bilayer systems. The interplay between atomic-level relaxation effects and moiré-scale electrostatic potential variation opens new avenues for the design of intertwined topological and correlated states, including the possibility of mimicking higher Landau level physics in the absence of magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An efficient random-sampling method for calculating double occupancy of Gutzwiller wave function in single-band 1D and 2D lattices

In this paper, we report a random sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for under- standing the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbor hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbor hopping and 2D lattices. In addition, constraints on real-space con gurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Selectivity of Electrochemical Ion Insertion into Manganese Dioxide Polymorphs

The ion insertion redox chemistry of manganese dioxide has diverse applications in energy storage, catalysis, and chemical separations. Unique properties derive from the assembly of Mn–O octahedra into polymorphic structures that can host protons and nonprotonic cations in interstitial sites. Despite many reports on individual ion-polymorph couples, much less is known about the selectivity of electrochemical ion insertion in MnO 2 . In this work, we use density functional theory to holistically compare the electrochemistry of A x MnO 2 (where A = H + , Li + , Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ , Al 3+ ) in aqueous and nonaqueous electrolytes. We develop an efficient computational scheme demonstrating that Hubbard-U correction has a greater impact on calculating accurate redox energetics than choice of exchange-correlation functional. Using PBE+U, we find that for nonprotonic cations, ion selectivity depends on the oxygen coordination environments inside a polymorph. When H + is present, however, the driving force to form hydroxyl bonds is usually stronger. In aqueous electrolytes, only three ion-polymorph pairs are thermodynamically stable within water’s voltage stability window (Na + and K + in α-MnO 2 , and Li + in λ-MnO 2 ), with all other ion insertion being metastable. We find Al 3+ may insert into the δ, R, and λ polymorphs across the full 2-electron redox of MnO 2 at high voltage; however, electrolytes for multivalent ions must be designed to impede the formation of insoluble precipitates and facilitate cation desolvation. We also show that small ions coinsert with water in α-MnO 2 to achieve greater coordination by oxygen, while solvation energies and kinetic effects dictate water coinsertion in δ-MnO 2 . Finally, taken together, these findings explain reports of mixed ion insertion mechanisms in aqueous electrolytes and highlight promising design strategies for safe, high energy density electrochemical energy storage, desalination batteries, and electrocatalysts.

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

Correlating Thermodynamic and Kinetic Hydricities of Rhenium Hydrides

The kinetics of hydride transfer from Re( R bpy)(CO) 3 H (bpy = 4,4'-R-2,2'-bipyridine; R = OMe, t Bu, Me, H, Br, COOMe, CF 3 ) to CO 2 and seven different cationic N-heterocycles were determined. Additionally, the thermodynamic hydricities of complexes of the type Re( R bpy)(CO) 3 H were established primarily using computational methods. Linear free-energy relationships (LFERs) derived by correlating thermodynamic and kinetic hydricities indicate that, in general, the rate of hydride transfer increases as the thermodynamic driving force for the reaction increases. Kinetic isotope effects range from inverse for hydride transfer reactions with a small driving force to normal for reactions with a large driving force. Hammett analysis indicates that hydride transfer reactions with greater thermodynamic driving force are less sensitive to changes in the electronic properties of the metal hydride, presumably because there is less buildup of charge in the increasingly early transition state. Bronsted α values were obtained for a range of hydride transfer reactions and along with DFT calculations suggest the reactions are concerted, which enables the use of Marcus theory to analyze hydride transfer reactions involving transition metal hydrides. It is notable, however, that even slight perturbations in the steric properties of the Re hydride or the hydride acceptor result in large deviations in the predicted rate of hydride transfer based on thermodynamic driving forces. This indicates that thermodynamic considerations alone cannot be used to predict the rate of hydride transfer, which has implications for catalyst design.

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

Drastic enhancement of magnetic critical temperature and amorphization in topological magnet EuSn 2 P 2 under pressure

High pressure is an effective tool to induce exotic quantum phenomena in magnetic topological insulators by controlling the interplay of magnetic order and topological state. This work presents a comprehensive high-pressure study of the crystal structure and magnetic ground state up to 62 GPa in an intrinsic topological magnet EuSn 2 P 2 . With a combination of high resolution X-ray diffraction, 151 Eu synchrotron Mössbauer spectroscopy, X-ray absorption spectroscopy, molecular orbital calculations, and electronic band structure calculations, it has been revealed that pressure drives EuSn 2 P 2 from a rhombohedral crystal to an amorphous phase at 36 GPa accompanied by a fourfold enhancement of magnetic ordering temperature. In the pressure-induced amorphous phase, Eu ions take an intermediate valence state. The drastic enhancement of magnetic ordering temperature from 30 K at ambient pressure to 130 K at 41.2 GPa resulting from Ruderman–Kittel–Kasuya–Yosida (RKKY) interactions likely attributes to the stronger Eu–Sn interaction at high pressure. These rich results demonstrate that EuSn 2 P 2 is an ideal platform to study the correlation of the enhanced RKKY interactions, disordered lattice, intermediate valence, and topological state.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Measurement of atom resolvability in cryo-EM maps with Q-scores

Cryogenic electron microscopy (cryo-EM) maps are now at the point where resolvability of individual atoms can be achieved. However, resolvability is not necessarily uniform throughout the map. We introduce a quantitative parameter to characterize the resolvability of individual atoms in cryo-EM maps, the map Q-score. Q-scores can be calculated for atoms in proteins, nucleic acids, water, ligands and other solvent atoms, using models fitted to or derived from cryo-EM maps. Q-scores can also be averaged to represent larger features such as entire residues and nucleotides. Averaged over entire models, Q-scores correlate very well with the estimated resolution of cryo-EM maps for both protein and RNA. Here, assuming the models they are calculated from are well fitted to the map, Q-scores can be used as a measure of resolvability in cryo-EM maps at various scales, from entire macromolecules down to individual atoms. Q-score analysis of multiple cryo-EM maps of the same proteins derived from different laboratories confirms the reproducibility of structural features from side chains down to water and ion atoms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Electrochemical control over stoichiometry via cation intercalation into Chevrel-phase sulphides (Cu x Mo 6 S 8 , x = 1-3)

Herein this work presents a systematic investigation of the electrochemical intercalation of aqueous copper cations into the Chevrel phase (CP) Mo 6 S 8 and its effect on the host's electronic and structural characteristics as a function of stoichiometry. The electrochemical potentials at which copper intercalates into the CP Mo 6 S 8 host in 1 M CuSO 4 were identified for Cu x Mo 6 S 8 (x = 1–3) as 0.23 V, 0.11 V, and 0.055 V (vs. Cu/Cu 2+ ) for Cu 1 Mo 6 S 8 , Cu 2 Mo 6 S 8 , and Cu 3 Mo 6 S 8 , respectively, and correlated via ex situ PXRD to stoichiometric content. Potentials were also identified that yield non-stoichiometric phases with traceable Cu content-dependent structural changes. Ab initio DFT calculations were performed to probe thermodynamic stability of Cu x Mo 6 S 8 CP compositions and the metal-to-semiconductor transitions as a function of copper content. Lastly, we demonstrate that the SK-pre-edge intensity is a successful probe for intercalant content: as copper equivalents increase, so does charge transfer to the partially occupied S 3p orbitals, resulting in a decrease in the pre-edge intensity corresponding to the excitation of a S 1s-to-3p/Mo 4d hybridized orbital.

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

Scalar Breit interaction for molecular calculations

Variational treatment of the Dirac–Coulomb–Gaunt or Dirac–Coulomb–Breit two-electron interaction at the Dirac–Hartree–Fock level is the starting point of high-accuracy four-component calculations of atomic and molecular systems. In this work, we introduce, for the first time, the scalar Hamiltonians derived from the Dirac–Coulomb–Gaunt and Dirac–Coulomb–Breit operators based on spin separation in the Pauli quaternion basis. While the widely used spin-free Dirac–Coulomb Hamiltonian includes only the direct Coulomb and exchange terms that resemble nonrelativistic two-electron interactions, the scalar Gaunt operator adds a scalar spin–spin term. The spin separation of the gauge operator gives rise to an additional scalar orbit-orbit interaction in the scalar Breit Hamiltonian. Benchmark calculations of Au n (n = 2–8) show that the scalar Dirac–Coulomb–Breit Hamiltonian can capture 99.99% of the total energy with only 10% of the computational cost when real-valued arithmetic is used, compared to the full Dirac–Coulomb–Breit Hamiltonian. In conclusion, the scalar relativistic formulation developed in this work lays the theoretical foundation for the development of high-accuracy, low-cost correlated variational relativistic many-body theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Connected triple excitations in coupled-cluster calculations of hyperpolarizabilities: Neon

We have calculated the second hyperpolarizability gamma of neon using the CCSD(T) method. The accuracy of the CCSD(T) approach has been established by explicit comparison with the single, double and triple excitation coupled-cluster (CCSDT) method using extended basis sets that are known to be adequate for the description of gamma. Our best estimate for gamma(sub 0) of 110 +/- 3 a.u. is in good agreement with other recent theoretical values and with Shelton's recent experimental estimate of 108 +/- 2 a.u. Comparison of the MP2 and CCSD(T) hyperpolarizability values indicates that MP2 gives a very good description of the electron correlation contribution to gamma(sub 0). We have combined MP2 frequency-dependent corrections with the CCSD(T) gamma(sub 0) to yield values of gamma(-2 omega;omega,omega,0) and gamma(exp K)(-omega;omega,0,0).

Rice, Julia E.↗
NASA NTRS2008

Symmetric Phase-Only Filtering in Particle-Image Velocimetry

Symmetrical phase-only filtering (SPOF) can be exploited to obtain substantial improvements in the results of data processing in particle-image velocimetry (PIV). In comparison with traditional PIV data processing, SPOF PIV data processing yields narrower and larger amplitude correlation peaks, thereby providing more-accurate velocity estimates. The higher signal-to-noise ratios associated with the higher amplitude correlation peaks afford greater robustness and reliability of processing. SPOF also affords superior performance in the presence of surface flare light and/or background light. SPOF algorithms can readily be incorporated into pre-existing algorithms used to process digitized image data in PIV, without significantly increasing processing times. A summary of PIV and traditional PIV data processing is prerequisite to a meaningful description of SPOF PIV processing. In PIV, a pulsed laser is used to illuminate a substantially planar region of a flowing fluid in which particles are entrained. An electronic camera records digital images of the particles at two instants of time. The components of velocity of the fluid in the illuminated plane can be obtained by determining the displacements of particles between the two illumination pulses. The objective in PIV data processing is to compute the particle displacements from the digital image data. In traditional PIV data processing, to which the present innovation applies, the two images are divided into a grid of subregions and the displacements determined from cross-correlations between the corresponding sub-regions in the first and second images. The cross-correlation process begins with the calculation of the Fourier transforms (or fast Fourier transforms) of the subregion portions of the images. The Fourier transforms from the corresponding subregions are multiplied, and this product is inverse Fourier transformed, yielding the cross-correlation intensity distribution. The average displacement of the particles across a subregion results in a displacement of the correlation peak from the center of the correlation plane. The velocity is then computed from the displacement of the correlation peak and the time between the recording of the two images. The process as described thus far is performed for all the subregions. The resulting set of velocities in grid cells amounts to a velocity vector map of the flow field recorded on the image plane. In traditional PIV processing, surface flare light and bright background light give rise to a large, broad correlation peak, at the center of the correlation plane, that can overwhelm the true particle- displacement correlation peak. This has made it necessary to resort to tedious image-masking and background-subtraction procedures to recover the relatively small amplitude particle-displacement correlation peak. SPOF is a variant of phase-only filtering (POF), which, in turn, is a variant of matched spatial filtering (MSF). In MSF, one projects a first image (denoted the input image) onto a second image (denoted the filter) as part of a computation to determine how much and what part of the filter is present in the input image. MSF is equivalent to cross-correlation. In POF, the frequency-domain content of the MSF filter is modified to produce a unitamplitude (phase-only) object. POF is implemented by normalizing the Fourier transform of the filter by its magnitude. The advantage of POFs is that they yield correlation peaks that are sharper and have higher signal-to-noise ratios than those obtained through traditional MSF. In the SPOF, these benefits of POF can be extended to PIV data processing. The SPOF yields even better performance than the POF approach, which is uniquely applicable to PIV type image data. In SPOF as now applied to PIV data processing, a subregion of the first image is treated as the input image and the corresponding subregion of the second image is treated as the filter. The Fourier transforms from both the firs and second- image subregions are normalized by the square roots of their respective magnitudes. This scheme yields optimal performance because the amounts of normalization applied to the spatial-frequency contents of the input and filter scenes are just enough to enhance their high-spatial-frequency contents while reducing their spurious low-spatial-frequency content. As a result, in SPOF PIV processing, particle-displacement correlation peaks can readily be detected above spurious background peaks, without need for masking or background subtraction.

Wemet, Mark P.↗
DOE OSTI2021

Electronic structure and topology across T c in the magnetic Weyl semimetal Co 3 Sn 2 S 2

Co 3 Sn 2 S 2 is a magnetic Weyl semimetal, in which ferromagnetic ordering at 177 K is predicted to stabilize Weyl points. We perform temperature and spatial dependent angle--resolved photoemission spectroscopy measurements through the Curie temperature (T c ), which show large band shifts and renormalization concomitant with the onset of magnetism. We argue that Co 3 Sn 2 S 2 evolves from a Mott ferromagnet below T c to a correlated metallic state above T c . To understand the magnetism, we derive a tight-binding model of Co-3d x 2 -y 2 orbitals on the kagome lattice. At the filling obtained by first-principles calculations, this model reproduces the ferromagnetic ground state, and results in the reduction of Coulomb interactions due to cluster effects. Using a disordered local moment simulation, we show how this reduced Hubbard U leads to a collapse of the bands across the magnetic transition, resulting in a correlated state which carries associated characteristic photoemission signatures that are distinct from those of a simple lifting of exchange splitting. Finally, the behavior of topology across T c is discussed in the context of this description of the magnetism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The Energy-Energy Correlation in the back-to-back limit at N 3 LO and N 3 LL'

We present the analytic formula for the Energy-Energy Correlation (EEC) in electron-positron annihilation computed in perturbative QCD to next-to-next-to-next-to-leading order (N 3 LO) in the back-to-back limit. In particular, we consider the EEC arising from the annihilation of an electron-positron pair into a virtual photon as well as a Higgs boson and their subsequent inclusive decay into hadrons. Our computation is based on a factorization theorem of the EEC formulated within Soft-Collinear Effective Theory (SCET) for the back-to-back limit. We obtain the last missing ingredient for our computation — the jet function — from a recent calculation of the transverse-momentum dependent fragmentation function (TMDFF) at N 3 LO. We combine the newly obtained N 3 LO jet function with the well known hard and soft function to predict the EEC in the back-to-back limit. The leading transcendental contribution of our analytic formula agrees with previously obtained results in N = 4 supersymmetric Yang-Mills theory. We obtain the N = 2 Mellin moment of the bulk region of the EEC using momentum sum rules. Finally, we obtain the first resummation of the EEC in the back-to-back limit at N 3 LL' accuracy, resulting in a factor of ~ 4 reduction of uncertainties in the peak region compared to N 3 LL predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

An energetic link between order and strength in metals: A nanocrystalline strength limit in high-entropy alloys and intermetallic compounds

The metallurgy and materials communities have long understood and exploited fundamental links between chemical and structural ordering in metallic solids to tailor their mechanical properties. We extend these ideas to include prediction of the nanocrystalline strength limit in high-entropy alloys and intermetallic compounds, where a breakdown occurs in the classical Hall-Petch strengthening behavior. The highest reported strength achievable through alloying has rapidly climbed and given rise to new classifications of materials with extraordinary properties, with a notable case being nanocrystalline metals. High-entropy alloys (chemically disordered, concentrated solid solutions) and intermetallic compounds are two boundary cases of how tailored order can be used to manipulate mechanical behavior. Here, we show that the complex electronic-structure mechanisms governing the peak strength of alloys and pure metals can be reduced to a few physically meaningful parameters based on their atomic arrangements and used – with no fitting parameters – to predict the maximum strength of these materials. This includes a generalized energy-based accounting for the degree of structural and chemical ordering that allows for rapid and reasonably accurate prediction of peak strength (validated in the nanocrystalline limit) as a function of temperature. Predictions of maximum strength based on the activation energy (with all materials properties derived from DFT calculations or experiments) for a stress-driven transition to an amorphous state is shown to accurately describe the breakdown in Hall-Petch behavior at the smallest crystallite sizes for pure metals, intermetallic compounds, high-entropy alloys, and metallic glasses. Further, this activation energy is also shown to be directly proportional to interstitial electronic charge density, which is a good predictor of ductility, stiffness (moduli), and phase stability in high-entropy alloys and solid metals generally. The proposed framework suggests the possibility of coupling ordering and intrinsic strength to mechanisms like dislocation nucleation, hydrogen embrittlement, and transport properties, such as through correlations between the activation energies for amorphization with stacking-fault and grain boundary energies. It additionally opens the prospect for greatly accelerated structural materials design and development to address materials challenges limiting more sustainable and efficient use of energy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantifying the effect of late bombardment on terrestrial zircons

The first 500 million years of Earth history is thought to be a period of intense planetary bombardment, but the timing and flux of this meteorite bombardment is poorly understood. In particular, on the basis of an inferred lunar impact history, some workers have hypothesized a ~3.9 Ga terminal cataclysm (TC) in which there was marked increase in the impact flux affecting the Moon, the Earth and possibly other terrestrial planets. Minerals that survived this enigmatic period offer a way to test early planetary bombardment models as they may contain telltale micro- to nanoscale shock features. Here, we present results from a numerical modeling calculation that assesses the probability that a zircon residing in the crust would escape shock melting or shock deformation during a TC bombardment event. Even with conservative pressure estimates for zircon shock deformation and intermediate bombardment intensities, we find that only ~6% of ≥4.0 Ga crust would be expected to survive a 3.9 Ga cataclysm without experiencing either complete melting or zircon shock metamorphism. In this study, we couple this modeling with a search for shock effects in the oldest zircons from the Acasta Gneiss Complex, which would have been present in the Earth's crust during a putative 3.9 Ga TC. Spatially correlated electron and NanoSIMS ion microscopy of 4.02 Ga igneous zircons from Acasta reveals no evidence of ancient shock. These data, together with similar results from other Hadean zircon suites, confirm that a post-Hadean TC is unlikely to have occurred. We suggest that the dearth of pre-3.9 Ga terrestrial crust and zircons is instead best explained by endogenic processes related to the mechanisms of early crust formation. Our modeling allows us to evaluate bombardment scenarios from the terrestrial zircon record by applying probabilistic interpretations to zircon shock deformation data. This approach will be valuable for other planetary bodies, allowing broader conclusions to be drawn from geographically limited datasets.

58 GEOSCIENCES↗
DOE OSTI2021

Direct observation of kink evolution due to Hund’s coupling on approach to metal-insulator transition in NiS 2- x Se x

Understanding characteristic energy scales is a fundamentally important issue in the study of strongly correlated systems. In multiband systems, an energy scale is affected not only by the effective Coulomb interaction but also by the Hund’s coupling. Direct observation of such energy scale has been elusive so far in spite of extensive studies. Here, we report the observation of a kink structure in the low energy dispersion of NiS 2- x Se x and its characteristic evolution with x , by using angle resolved photoemission spectroscopy. Dynamical mean field theory calculation combined with density functional theory confirms that this kink originates from Hund’s coupling. We find that the abrupt deviation from the Fermi liquid behavior in the electron self-energy results in the kink feature at low energy scale and that the kink is directly related to the coherence-incoherence crossover temperature scale. Our results mark the direct observation of the evolution of the characteristic temperature scale via kink features in the spectral function, which is the hallmark of Hund’s physics in the multiorbital system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exchange-correlation functional challenges in modeling quaternary chalcogenides

The development of next-generation quaternary chalcogenides, such as Cu2ZnSnS4 (CZTS) and Cu2ZnGeS4 (CZGS), for solar energy and thermoelectric applications hinges upon both careful experimentation and accurate quantum mechanical modeling. To address the latter, many have turned to density functional theory (DFT), which offers several choices for the approximate treatment of electron exchange and correlation (XC). Popular XC functionals include the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) and the recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. Extensions of DFT functionals, such as adding a Hubbard U correction and introducing a fraction of the Fock exchange (hybrid functionals), have been used widely to model systems containing 3d metal ions. However, no studies yet have compared comprehensively PBE(+U) and SCAN(+U) in the quality of their predictions of the bulk and defect thermodynamics of quaternary chalcogenides, which play a critical role in device fabrication and performance. Hence, here we calculate the (i) 0 K formation energies of bulk Ge compounds and (ii) neutral defect formation energies including charge-balanced (e.g., CuZn + ZnCu) and charge-imbalanced (e.g., CuSn) combinations of antisites and vacancies in CZTS and CZGS using the PBE, PBE +U, SCAN, SCAN +U, and the hybrid Heyd-Scuseria-Ernzerhof XC frameworks. We find that the formation energies of charge-imbalanced defects are more sensitive to the choice of the XC functional than those of charge-balanced defects, which can be explained by the differences in the extent of penalization of defect-generated delocalized electrons/holes by PBE, PBE +U, SCAN, and SCAN +U. Additionally, our results show that SCAN systematically underbinds Ge-containing compounds, thus highlighting the need for even further improvement of XC functionals. Based on our findings, we recommend the use of SCAN for modeling quaternary chalcogenides because its errors are systematic, and it has the firmest theoretical underpinning. Our work provides guidance for future modeling of quaternary chalcogenides.

14 SOLAR ENERGY↗