Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Configuration interaction”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An over-approximation of entropy for elemental multiconfigurational ground state electronic structures

The ability to predict the chemical and physical properties of a material is directly related to the structure and interactions of its electrons. For materials comprised of f-block elements (the lanthanides and actinides found in the last two rows of the periodic table), the complexity of electronic structure has presented great difficulty in understanding, modeling, and predicting material properties. The complexity of multiconfigurational ground state electronic structures is illustrated herein by the combinatorics of electron permutations within individual and cumulative occupancy configurations. A non-integer orbital occupancy representation of multiconfigurational ground states is described for superposition mixing between multiple near-energy degenerate occupancy configurations and generalized in such a way that established ground states are returned by approximation for elements with less-complex electronic structures. By considering the occupancy configurations as statistical mechanics macrostates, and the permutations of electrons as statistical mechanics microstates within those macrostates, an over-approximation of entropy for multiconfigurational elemental ground state electronic structures has been calculated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

The Energy Exascale Earth System Model Version 3: 2. Overview of the Coupled System

The Energy Exascale Earth System Model version 3 (E3SMv3) represents the latest advancement in Earth system modeling developed by the U.S. Department of Energy (DOE). Building upon previous versions, E3SMv3 introduces significant updates across its coupled components to enhance capability and improve fidelity. The atmosphere component incorporates advancements in chemistry, aerosol-cloud interactions, convection, and microphysics. The ocean features a new time-stepping scheme and a higher-resolution unstructured mesh with sub-ice-shelf cavities, while the sea ice model integrates advanced snow and ice physics for more realistic cryospheric simulations. The land model introduces prognostic vegetation dynamics and a new sub-grid topographic treatment of solar radiation. A new tri-grid configuration harmonizes the horizontal grids of the land and river components for improved process coupling. It is enabled by a new non-linear remapping between the atmosphere and land. E3SMv3 underwent extensive testing through a comprehensive simulation campaign, including pre-industrial control, idealized CO 2 experiments, and historical simulations spanning 1850–2024. The model demonstrates significant improvements in simulating the evolution of the historical surface temperature, particularly addressing the “pothole cooling” bias in earlier versions. Reduced aerosol-related forcing contributes to more realistic radiative forcing and better alignment with the observational record. Ocean heat content (OHC) and sea ice trends are also improved as a result.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Collective neural network behavior in a dynamically driven disordered system of superconducting loops

Collective properties of complex systems composed of many interacting components such as neurons in our brain can be modeled by artificial networks based on disordered systems. We show that a disordered neural network of superconducting loops with Josephson junctions can exhibit computational properties like categorization and associative memory in the time evolution of its state in response to information from external excitations. Superconducting loops can trap multiples of fluxons in many discrete memory configurations defined by the local free energy minima in the configuration space of all possible states. A memory state can be updated by exciting the Josephson junctions to fire or allow the movement of fluxons through the network as the current through them surpasses their critical current thresholds. Simulations performed with a lumped element circuit model of a 4-loop network show that information written through excitations is translated into stable states of trapped flux and their time evolution. Experimental implementation on a high-Tc superconductor YBCO-based 4-loop network shows dynamically stable flux flow in each pathway characterized by the correlations between junction firing statistics. Neural network behavior is observed as energy barriers separating state categories in simulations in response to multiple excitations, and experimentally as junction responses characterizing different flux flow patterns in the network. The state categories that produce these patterns have different temporal stabilities relative to each other and the excitations. This provides strong evidence for time-dependent (short-to-long-term) memories, that are dependent on the geometrical and junction parameters of the loops, as described with a network model.

Josephson junctions↗
DOE OSTI2021

The DAQ system of the 12,000 channel CMS high granularity calorimeter prototype

The CMS experiment at the CERN LHC will be upgraded to accommodate the 5-fold increase in the instantaneous luminosity expected at the High-Luminosity LHC (HL-LHC). Concomitant with this increase will be an increase in the number of interactions in each bunch crossing and a significant increase in the total ionising dose and fluence. One part of this upgrade is the replacement of the current endcap calorimeters with a high granularity sampling calorimeter equipped with silicon sensors, designed to manage the high collision rates. As part of the development of this calorimeter, a series of beam tests have been conducted with different sampling configurations using prototype segmented silicon detectors. In the most recent of these tests, conducted in late 2018 at the CERN SPS, the performance of a prototype calorimeter equipped with ≈12,000 channels of silicon sensors was studied with beams of high-energy electrons, pions and muons. This paper describes the custom-built scalable data acquisition system that was built with readily available FPGA mezzanines and low-cost Raspberry Pi computers.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Probing proton cross-shell excitations through the two-neutron removal from 38 Ca

Bound states of the neutron-deficient, near-dripline nucleus 36 Ca were populated in two-neutron removal from the ground state of 38 Ca, a direct reaction sensitive to the single-particle configurations and couplings of the removed neutrons in the projectile wave function. Final-state exclusive cross sections for the formation of 36 Ca and the corresponding longitudinal momentum distributions, both determined through the combination of particle and γ-ray spectroscopy, are compared to predictions combining eikonal reaction theory and shell-model two-nucleon amplitudes from the USDB, USDC, and ZBM2 effective interactions. The final-state cross-section ratio σ($2$$^{+}_{1}$) / σ (0 + ) shows particular sensitivity and is approximately reproduced only with the two-nucleon amplitudes from the ZBM2 effective interaction that includes proton cross-shell excitations into the pf shell. Characterizing the proton pf-shell occupancy locally and schematically, an increase of the sd – pf shell gap by 250 keV yields an improved description of this cross-section ratio and simultaneously enables a reproduction of the B(E2;$0$$^{+}_{1}$ → $2$$^{+}_{1}$) excitation strength of 36 Ca. Furthermore, this highlights an important aspect if a new shell-model effective interaction for the region was to be developed on the quest to model the neutron-deficient Ca isotopes and surrounding nuclei whose structure is impacted by proton cross-shell excitations.

20 ≤ A ≤ 38↗
DOE OSTI2025

A Versatile Simulated Data Transport Layer for in Situ Workflows Performance Evaluation

In situ processing does not only allow scientific applications to face the explosion in data volume and velocity but also to address the time constraints of many simulation-analysis workflows by providing scientists with early insights about their applications at runtime. Multiple frameworks implement the concept of a data transport layer (DTL) to enable such in situ workflows. These tools are very versatile, directly or indirectly access the data generated on the same node, another node of the same compute cluster, or a completely distinct node, and allow data publishers and subscribers to run on the same computing resources or not. This versatility puts on researchers the onus of taking key decisions related to resource allocation and how to transport data to ensure the most efficient execution of their in situ workflows. However, domain scientists and workflow practitioners lack the appropriate tools to assess the respective performance of particular design and deployment options. In this paper we introduce a versatile simulated DTL designed to provide researchers with insights on the respective performance of different execution scenarios of in situ workflows. This open-source, standalone library builds on the SimGrid toolkit and can be linked to any SimGrid-based simulator. It facilitates the evaluation of the performance behavior, at scale, of different data transport configurations and the study of the effects of resource allocation strategies. We demonstrate the scalability, versatility, and accuracy of this simulated DTL by reproducing the execution of two synthetic benchmarks and of a real-world in situ workflow composed of an MPI application and a parallel data analysis. Results of simulations run on a single core show that the proposed library can simulate the interactions of tens of thousands of simulated processes deployed on two interconnected commodity clusters in a few seconds, and the execution by a thousand simulated processes of an in situ workflow in less than three minutes.

Suter, Fred [ORNL] (ORCID:0000000319021955)↗
DOE OSTI2023

Uniform Distribution Technique for Neutrino Beam Scan Simulation

In Fermilab's neutrino facilities such as the Neutrinos at the Main Injector (NuMI) and the upcoming Long Baseline Neutrino Facility (LBNF), a proton beam strikes high-power target, producing positively and negatively charged pions and kaons. There is a need for detailed simulations in order to capture all particle interactions and beam propagation from protons on target to short-lived mesons decaying into muons and neutrinos. The generation of individual beam simulations is a resource-intensive and time-consuming process. In this paper, we describe a method through which many simulation samples with high statistics can be generated to study the effects of beam scan across a target for given beam configurations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

From oversimplified to overlooked: The case for exploring rich dark sectors

The Standard Model (SM) of particle physics provides a very successful description of fundamental particles and their interactions but it is incomplete, as neutrino masses, dark matter and the baryon asymmetry of the Universe indicate. In addition, the origin of masses and of the approximate fundamental symmetries call out for deeper explanations. The quest for a New SM Theory, that extends the SM to a more general theory, is ongoing. For decades the main focus has been on the TeV scale, but despite an impressive theoretical and experimental effort, no hints of new physics at such scale has been found in experiments. Dark sectors provide an interesting alternative to TeV scale extensions of the SM to explain the open questions in particle and astroparticle physics. Going beyond minimal models, rich dark sectors extend the SM to a complex theory with multiple particles and interactions, in analogy to the SM itself. They have a wealth of theoretical and astrophysical/cosmological consequences and can lead to phenomenological signatures that can be markedly different to that of minimal ones. These include short-lived particles and semi-visible decay signatures, as opposed to minimal models where new states are typically long-lived and purely visible or invisible resonances. Given the experimental configurations and analysis strategies, current dark sector searches might miss such signatures. We advocate a dedicated programme of searches for rich dark sectors that overcomes the assumptions on minimality and on the long lifetime of particles and encompasses a broader range of possibilities. Here, we discuss a prototype model that includes a complex structure akin to the SM: multiple generations of fermions charged under a new spontaneously-broken gauge symmetry.

Dark matter↗
DOE OSTI2023

Role of Anions in Stabilizing the [Zn–Al] Layered Double Hydroxides: A Thermodynamic Study

Here, room-temperature acid solution calorimetry, high-temperature oxide melt solution calorimetry, and low-temperature heat capacity measurements were employed to calculate the thermodynamic stabilities of the [Zn–Al–X] layered double hydroxides (LDH) containing different anions (X = Cl – , CO 3 2– , and SO 4 2– ). Cryogenic heat capacity measurements demonstrated a Schottky-type anomaly in the heat capacity of all three LDHs below 11 K. This anomaly is attributed to the tunneling of protons between adjacent oxygen atoms in the LDH interlayer as this creates an energy system similar to a two-level system modeled with a Schottky term. These heat capacity measurements were also used to determine vibrational entropies which, when combined with configurational entropies, provide standard entropies of these LDHs. Enthalpies of formation of LDHs from binary components were determined and combined with the entropies of formation to calculate Gibbs free energies. Based on these values, the order of stability is [Zn–Al–SO 4 ] > [Zn–Al–CO 3 ] > [Zn–Al–Cl]. This trend results from a combination of the interlayer spacing, amount of water in the interlayer, interactions among the interlayer species, and interactions between the metal hydroxide layer and the interlayer.

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

Heterostructured Lepidocrocite Titanate-Carbon Nanosheets for Electrochemical Applications

Lepidocrocite-type titanates that reversibly intercalate sodium ions at low potentials (~0.6 V vs Na/Na + ) are promising anode candidates for sodium-ion batteries. However, large amounts of carbon additives are often used to improve their electrical conductivity and overcome poor cycling performance in the electrode composites. To ameliorate electronic transport issues of lepidocrocite titanate (K 0.8 Ti 1.73 Li 0.27 O 4 , KTL) in sodium-ion batteries, we have designed and synthesized heterostructures of exfoliated lepidocrocite-type titanium oxide (LTO) nanosheets with alternating carbon layers via a solution-based self-assembly approach. Positively charged dopamine (Dopa) was used as the carbon precursor and intercalated between negatively charged exfoliated titania nanosheets through electrostatic interaction. Dopa-intercalated LTO was then annealed under argon to form conductive carbon layers between titania sheets. The carbon content in the heterostructures was controlled by modifying the self-assembly conditions (i.e., pH, stirring duration, and Dopa-to-LTO ratio). Electrodes were prepared using carbonized heterostructures (LTO-C) without adding more carbon to the composites and tested in sodium half-cell configurations. Further, higher capacities and improved capacity retention over 250 cycles and lower impedance were observed, as the carbon content of LTO-C heterostructures was increased from 0% (LTO nanosheets) to 30%. These results indicate that the self-assembly approach for 2D heterostructured electrode materials is a promising strategy to overcome electronic transport limitations of layered transition-metal oxides and improve their electrochemical performance for next-generation energy storage applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Relaxation of a single defect spin by the low-frequency gyrotropic mode of a magnetic vortex

We excite the gyrotropic mode of a magnetic vortex and observe the resulting effect on the spin state of a nearby nitrogen-vacancy (NV) defect in diamond. Thin permalloy disks fabricated on a diamond sample are magnetized in a vortex state in which the magnetization curls around a central core. The magnetization dynamics of this configuration are described by a discrete spectrum of confined magnon modes as well as a low-frequency gyrotropic mode in which the vortex core precesses about its equilibrium position. Despite the spin transition frequencies being far-detuned from the modes of the ferromagnet, we observe enhanced relaxation of the NV spin when driving the gyrotropic mode. Moreover, we map the spatial dependence of the interaction between the vortex and the spin by translating the vortex core within the disk with an applied magnetic field, resulting in steplike motion as the vortex is pinned and de-pinned. Strong spin relaxation is observed when the vortex core is within approximately 250 nm of the NV center defect. We attribute this effect to the higher frequencies in the spectrum of the magnetic fringe field arising from the soliton-like nature of the gyrotropic mode when driven with sufficiently large amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High-throughput native mass spectrometry as experimental validation for in silico drug design

In this project, we developed automated workflows for both experimental validation and computational prediction of protein-ligand interactions. The ultimate goal is to establish an integrated pipeline for high throughput design of inhibitors to enzymes relevant to all areas of biological research. Our experimental approach is based on native mass spectrometry (native MS), which measures accurate masses and quantify the relative abundance of protein-ligand complexes to define binding affinity. We set up an in-house built autosampler with highly flexible configurations to minimize the manual steps for high throughput native MS. In parallel, we also performed manual native MS to characterize the binding of substrates and inhibitors of SARS-Cov-2 nonstructural protein nsp10/16 in order to optimize the experimental parameters for future automation. On the computational side, we streamlined the pipeline to achieve minimal manual intervention for predicting enzyme inhibitors via simulation, using the same nsp10/16 system as an example. Using the native MS method we examined 8 top-ranked designed compounds, 2 of which showed weak binding of ~50 µM. The information from native MS experiment provided critical insights and the foundation for a fully integrated workflow for enzyme inhibitor design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

GLEAM: Galaxy Line Emission & Absorption Modeling

We present Galaxy Line Emission & Absorption Modeling (gleam), a Python tool for fitting Gaussian models to emission and absorption lines in large samples of 1D extragalactic spectra. gleam is tailored to work well in batch mode without much human interaction. With gleam, users can uniformly process a variety of spectra, including galaxies and active galactic nuclei, in a wide range of instrument setups and signal-to-noise regimes. gleam also takes advantage of multiprocessing capabilities to process spectra in parallel. With the goal of enabling reproducible workflows for its users, gleam employs a small number of input files, including a central, user-friendly configuration in which fitting constraints can be defined for groups of spectra and overrides can be specified for edge cases. For each spectrum, gleam produces a table containing measurements and error bars for the detected spectral lines and continuum and upper limits for nondetections. For visual inspection and publishing, gleam can also produce plots of the data with fitted lines overlaid. In the present paper, we describe gleam’s main features, the necessary inputs, expected outputs, and some example applications, including thorough tests on a large sample of optical/infrared multi-object spectroscopic observations and integral field spectroscopic data. gleam is developed as an open-source project hosted at https://github.com/multiwavelength/gleam and welcomes community contributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The NuMI Flux Prediction at ICARUS

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long-baseline neutrino oscillation experiment seeking to probe fundamental symmetries within the structure of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix, and perform precision measurements its parameters including the neutrino mass ordering via the sign of $\Delta m^2_{31}$, and the charge-parity violating phase, $\delta_{CP}$. To make these measurements with high precision, DUNE will require external $\nu$-Ar scattering cross section data as a crucial input to the oscillation fit. Imaging Cosmic And Rare Underground Signals (ICARUS) is a 476 t liquid argon neutrino detector located at Fermi National Accelerator Laboratory (FNAL) where it is serving as the far detector for the Short-Baseline Neutrino (SBN) program along the Booster Neutrino Beam (BNB) axis. ICARUS additionally lies 795 m downstream and 100.1 mrad off-axis of the Neutrinos at the Main Injector (NuMI) neutrino beam. From this position, ICARUS is exposed to a large flux of NuMI (anti-)electron and (anti-)muon neutrinos, and poses a unique opportunity to provide high-statistics measurements of quasi-elastic and single pion-production cross sections for four neutrino flavors ($\nu_{\mu}$, $\nu_{e}$, $\bar{\nu}_{\mu}$, $\bar{\nu}_{e}$). This dissertation is centered around accurately characterizing the models and estimating their precision for use in making these measurements. This includes identifying major sources of uncertainty in the models such that they can be properly propagated to the cross section measurements. Specifically, this work focused on the model of the NuMI beamline and its impact on the neutrino fluxes, but also delved into the detector response model and its impact on reconstructed observables in the detector. Significant efforts were made to improve the characterization to enhance precision and thus reduce the level of propagated uncertainty. In particular, the NuMI flux was determined to be composed of 57% $\nu_{\mu}$, 38% $\bar{\nu}_{\mu}$, 3% $\nu_{e}$, and 2% $\bar{\nu}_{e}$ while the horns are operating in the positive-particle focusing configuration. The total uncertainty on the $\nu_{\mu} + \bar{\nu}_{\mu}$ ($\nu_{e} + \bar{\nu}_{e}$) flux while operating in the forward horn operating mode was determined to be 10.84% (9.04%). Compared to the on-axis flux, mesons that eventually decay to neutrinos more frequently reinteract within the NuMI structure, resulting in elevated uncertainty as these processes are not well-constrained by existing hadron interaction cross section measurements. Covariance matrices were calculated to propagate the flux uncertainty characterization to NuMI analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Electronic structure of mononuclear Cu-based molecule from density-functional theory with self-interaction correction

In this paper, we investigate the electronic structure of a planar mononuclear Cu-based molecule [Cu(C 6 H 4 S 2 ) 2 ] z in two oxidation states ( z = –2, –1) using density-functional theory (DFT) with Fermi–Löwdin orbital (FLO) self-interaction correction (SIC). The dianionic Cu-based molecule was proposed to be a promising qubit candidate. Self-interaction error within approximate DFT functionals renders severe delocalization of electron and spin densities arising from 3 d orbitals. The FLO-SIC method relies on optimization of Fermi–Löwdin orbital descriptors (FODs) with which localized occupied orbitals are constructed to create SIC potentials. Starting with many initial sets of FODs, we employ a frozen-density loop algorithm within the FLO-SIC method to study the Cu-based molecule. We find that the electronic structure of the molecule remains unchanged despite somewhat different final FOD configurations. In the dianionic state (spin S = 1/2), FLO-SIC spin density originates from the Cu d and S p orbitals with an approximate ratio of 2:1, in quantitative agreement with multireference calculations, while in the case of SIC-free DFT, the orbital ratio is reversed. Overall, FLO-SIC lowers the energies of the occupied orbitals and, in particular, the 3 d orbitals unhybridized with the ligands significantly, which substantially increases the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) compared to SIC-free DFT results. The FLO-SIC HOMO–LUMO gap of the dianionic state is larger than that of the monoanionic state, which is consistent with experiment. Our results suggest a positive outlook of the FLO-SIC method in the description of magnetic exchange coupling within 3 d -element-based systems.

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

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗