Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hybrid density functional 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Ab Initio Study of Atomic and Electronic Structure of Promising Ba4XMn3O12 (X = Nb, Ce, Pr) Oxides for Solar Thermochemical Hydrogen Production

The two-step metal oxide water-splitting cycle is one of the most viable approach for Solar Thermochemical Hydrogen (STCH) production. Challenges exist in finding suitable oxides that can satisfy thermodynamics of the STCH redox cycle under viable range of temperatures and partial pressures. Recently, BXM family members, i.e., Ba4NbMn3O12 (BNM), Ba4CeMn3O12 (BCM), Ba4PrMn3O12 (BPM) are discovered to exhibit promising STCH performance. However, the magnetic degrees of freedom of their experimental crystal structures is not characterized, and only hypothetical models exist for their electronic structure. We performed Monte-Carlo sampling for the magnetic spin arrangement of Mn atoms to provide explicit atomic structure models for their ground state 12R polytype and determined the most stable spin configuration for each member after structure relaxation using density functional theory (DFT) calculations. We also elucidate the most stable charge configuration (Mn3+/Mn4+ arrangement) for BNM among the existing hypothesis in literature. We investigate these strongly correlated oxides using different levels of theory (DFT+U, Hybrid-DFT) to develop better understanding of their electronic structure. These calculations allowed us to develop a model based on molecular orbital interactions that help explain the observed changes in their experimental x-ray spectroscopy measurements. Thereby, help link differences in electronic structure of BXM family members to their varying water-splitting behavior. These structure models will further facilitate future in-depth defect studies to model STCH redox cycle of these oxides.

atomic and magnetic structure↗
DOE OSTI2022

Machine learning for impurity charge-state transition levels in semiconductors from elemental properties using multi-fidelity datasets

Quantifying charge-state transition energy levels of impurities in semiconductors is critical to understanding and engineering their optoelectronic properties for applications ranging from solar photovoltaics to infrared lasers. While these transition levels can be measured and calculated accurately, such efforts are time-consuming and more rapid prediction methods would be beneficial. Here, we significantly reduce the time typically required to predict impurity transition levels using multi-fidelity datasets and a machine learning approach employing features based on elemental properties and impurity positions. We use transition levels obtained from low-fidelity (i.e., local-density approximation or generalized gradient approximation) density functional theory (DFT) calculations, corrected using a recently proposed modified band alignment scheme, which well-approximates transition levels from high-fidelity DFT (i.e., hybrid HSE06). Further, the model fit to the large multi-fidelity database shows improved accuracy compared to the models trained on the more limited high-fidelity values. Crucially, in our approach, when using the multi-fidelity data, high-fidelity values are not required for model training, significantly reducing the computational cost required for training the model. Our machine learning model of transition levels has a root mean squared (mean absolute) error of 0.36 (0.27) eV vs high-fidelity hybrid functional values when averaged over 14 semiconductor systems from the II–VI and III–V families. As a guide for use on other systems, we assessed the model on simulated data to show the expected accuracy level as a function of bandgap for new materials of interest. Finally, we use the model to predict a complete space of impurity charge-state transition levels in all zinc blende III–V and II–VI systems.

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

Molecular Level Understanding of Polyethylene Terephthalate (PET) Depolymerization in Base/Alcohol Hybrid Systems

Polyethylene terephthalate (PET) depolymerization in base/alcohol hybrid systems represents a promising low-energy approach for chemically recycling PET waste into valuable monomers. This study investigates the mechanistic pathways of PET depolymerization in NaOH/alcohol solutions, emphasizing the competing roles of hydroxide and alkoxide species. Utilizing a combination of experimental techniques, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations, we explore how factors such as base concentration, alcohol chain length, and pKa values of alcohols influence PET depolymerization efficiency and pathways. Our findings indicate that alkoxide ions (RO⁻) exhibit notably higher reactivity than hydroxide ions (HO⁻), favoring an alcoholysis pathway in the base/alcohol hybrid system. Experimental results across a series of C1 to C5 alcohols show that longer-chain alcohols, particularly 1-butanol, achieve higher PET conversion, although this does not align solely with simple nucleophilicity trends of alkoxides. While DFT calculations reveal comparable activation energies for various alkoxides in PET depolymerization, MD simulations underscore the significant role of alcohol chain length, with longer-chain alcohols forming more stable or frequent interactions with PET. Additionally, the alkoxide concentration, influenced by the alcohol’s pK a , directly impacts PET conversion. These suggest that PET depolymerization is governed by a balance between alkoxide concentration and alkoxide-PET interactions, rather than activation energies or nucleophilicity alone. From a practical perspective, incorporating long-chain alcohols as cosolvents may enhance process efficiency but increases raw material costs by approximately 30%. However, long-chain alcohols present a safer and more sustainable alternative to hazardous cosolvents such as dichloromethane. This work offers a molecular-level understanding of PET depolymerization in base/alcohol systems and provides insights into optimizing these systems for more efficient and sustainable PET recycling processes.

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

Magnetic exchange interactions in binuclear and tetranuclear iron(III) complexes described by spin‐flip DFT and Heisenberg effective Hamiltonians

Abstract Low‐energy spectra of single‐molecule magnets (SMMs) are often described by Heisenberg Hamiltonians. Within this formalism, exchange interactions between magnetic centers determine the ground‐state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants ( J ) for a set of iron (III) binuclear and tetranuclear complexes from all‐electron calculations using non‐collinear spin‐flip time‐dependent density functional theory (NC‐SF‐TDDFT). For 12 binuclear complexes with J ‐values ranging from −6 to −132 cm −1 , our benchmark calculations using the short‐range hybrid ω PBEh functional and 6‐31G(d,p) basis set agree well with the experimentally derived values (mean absolute error of 4.7 cm −1 ). For the tetranuclear SMMs, the computed J constants are within 6 cm −1 from the experimentally derived values. We explore the range of applicability of the Heisenberg model by analyzing bonding patterns in these Fe(III) complexes using natural orbitals (NO), their occupations, and the number of effectively unpaired electrons. The results illustrate the efficiency of the spin‐flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.

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

Discovering Dinuclear Dioxygen-Bridged Cobalt(III) Complexes for Selective Binding of O 2 from Air

The design and development of dioxygen activation in porous crystalline materials is a useful avenue for exploring selective adsorption of O 2 that shows significant potential to enable separations of O 2 and N 2 from air. Porous materials featuring redox-active metal centers have received attention regarding selective O 2 adsorption via chemisorption. Drawing inspiration from a dinuclear cobalt material ([(Co(III) 2 (bpbp)O 2 ) 2 bdc](PF 6 ) 4 (CSD code: GAMVIB; bpbp – = 2,6-bis(N,N-bis(2-pyridylmethyl)aminomethyl)-4-tert-butylphenolato; bdc 2– = 1,4-benzenedicarboxylato)) that displays reversible and selective O 2 adsorption, we focus on searching for potential O 2 -selective materials with dinuclear cobalt clusters that have dioxygen-bridged Co(III) complexes. We combine structure screening with a high-level hybrid periodic density functional theory (DFT) workflow to investigate O 2 and N 2 adsorption in materials from validated crystal structure databases (e.g., the CSD database). These calculations identify multiple materials that are predicted to have superior O 2 binding capability relative to GAMVIB. Grand Canonical Monte Carlo (GCMC) simulations based on DFT-developed force fields were performed for selected candidates to estimate the adsorption performance for O 2 /N 2 mixtures.

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

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗
DOE OSTI2025

Strain Effects in SrHfO 3 Films Grown by Hybrid Molecular Beam Epitaxy

Perovskite oxide heterostructures host a large number of interesting phenomena such as ferroelectricity, which are often driven by octahedral distortions in the crystal that may induce polarization. SrHfO 3 (SHO) is a perovskite oxide with a pseudocubic lattice parameter of 4.08 Å that previous density functional theory (DFT) calculations suggest can be stabilized in a ferroelectric P4mm phase when stabilized with sufficient compressive strain. Additionally, it is insulating and possesses a large band gap and a high dielectric constant, making it an ideal candidate for oxide electronic devices. Here, to test the viability of epitaxial strain as a driver of ferroic phase transitions, SHO films were grown by hybrid molecular beam epitaxy (hMBE) with a tetrakis(ethylmethylamino)hafnium(IV) source on GdScO 3 and TbScO 3 substrates. Strained SHO phases were characterized using X-ray diffraction, X-ray absorption spectroscopy, and scanning transmission electron microscopy to determine the space group of the strained films, with the results compared to those of DFT-optimized models of phase stability versus strain. Contrary to past reports, we find that compressively strained SrHfO 3 undergoes octahedral tilt distortions without associated ferroelectric polarization and most likely takes on the I4/mcm phase with the a 0 a 0 c – tilt pattern.

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

Enhanced Guest@MOF Interaction via Stepwise Thermal Annealing: TCNQ@Cu 3 (BTC) 2

Confinement of guest molecules in porous materials such as metal organic frameworks (MOFs) promises to deliver emergent properties separate from those of the individual components. Understanding the confinement mechanism is therefore important for the development of new synthesis routes that adjust MOF properties for specific applications. In this work, we developed a new synthetic method to confine guest molecules into MOF pores through a stepwise thermal annealing process, wherein the confinement of 7,7,8,8-tetracyanoquinodimethane (TCNQ) guest molecules into Cu 3 (BTC) 2 (BTC = benzene-1,3,5-tricarboxylic acid) MOF host is used as an example of how novel materials can be created with new physical properties. The stepwise thermal annealing process includes 1) an activation process of pristine Cu 3 (BTC) 2 MOF to maximizes the TCNQ guest loading in the MOF host by effectively removing the residual solvents and 2) post-annealing of the TCNQ infiltrated MOF to enhances the interaction of the confined guest molecules with the MOF host. Obtained experimental results based on thermogravimetric analysis, N 2 gas adsorption, electron microscopy, X-ray diffraction and infrared absorption, combined with density functional theory calculations provide evidence that the use of a stepwise thermal annealing process yields enhancements in the guest loading, packing and interaction between the TCNQ guest and the MOF host. The new hybrid TCNQ@Cu 3 (BTC) 2 system is stable and shows no significant signs of structural degradation even after submersion in water. This is due to the presence of significantly stronger interactions of TCNQ with the frame-work metal ions compared to those of the water molecules competing for the same framework binding sites. It was also found that TCNQ@Cu 3 (BTC) 2 system maintains a significant CO 2 and CH 4 adsorption potential compared to the pristine MOF. The synthetic route developed in this work yields novel guest@MOF hybrid systems that will be useful for many MOF-based applications such as gas separations and chemical sensors performed under humid conditions.

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

Accelerated coupled Monte Carlo-Thermal hydraulic calculations using a hybrid GTF-diffusion-based prediction block: first results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections' variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint Eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Accelerated Coupled Monte Carlo-Thermal Hydraulic Calculations using a Hybrid GTF-Diffusion-based Prediction Block: First Results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work [1, 2] investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections’ variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Synthetic Control of Water-Stable Hybrid Perovskitoid Semiconductors

Hybrid metal-halide perovskites and their derived materials have emerged as the next-generation semiconductors with a wide range of applications, including photovoltaics, light-emitting devices, and other optoelectronics. Over the past decade, numerous single-crystalline perovskite derivatives have been synthesized and developed. However, the synthetic methods for these derivatives mainly rely on acidic crystallization conditions. This approach leads to crystals comprising metal halide building blocks, which show problematic stability when directly exposed to water. In this study, a methodology is developed for synthesizing hybrid metal-halide compounds using lead iodide and the zwitterionic bifunctional molecule cysteamine (CYS), to form various perovskitoid structures under a broad pH range. Interestingly, the different pH conditions alter the coordination environment of lead halides, leading to lead-sulfide and lead-nitride covalent bond formation. This modification significantly enhances their stability when in direct contact with water, lasting for months. Photoluminescence measurements and first principal density functional theory (DFT) calculations reveal that the perovskitoids synthesized under basic and acidic pH conditions exhibit a direct bandgap nature, while those synthesized under neutral conditions display an indirect bandgap. This approach opens new avenues for manipulating synthetic methods to develop water-stable hybrid semiconductors suitable for a wide range of applications, such as solid-state light emitters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spin-polarized imaging of strongly interacting fermions in the ferrimagnetic state of the Weyl candidate CeBi

CeBi has an intricate magnetic phase diagram whose fully polarized state has recently been suggested as a Weyl semimetal, though the role of f states in promoting strong interactions has remained elusive. Here we focus on the less-studied but also time-reversal symmetry-breaking ferrimagnetic phase of CeBi, where our density functional theory (DFT) calculations predict additional Weyl nodes near the Fermi level E F . We use spin-polarized scanning tunneling microscopy and spectroscopy to image the surface ferrimagnetic order on the itinerant Bi p states, indicating their orbital hybridization with localized Ce f states. We observe suppression of this spin-polarized signature at E F , coincident with a Fano line shape in the conductance spectra, suggesting the Bi p states partially Kondo screen the f magnetic moments, and this p – f hybridization causes strong Fermi-level band renormalization. Furthermore, the p-band flattening is supported by our quasiparticle interference measurements, which also show band splitting in agreement with DFT, painting a consistent picture of a strongly interacting magnetic Weyl semimetal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

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

Version 2.0.0 - M-SPARC: Matlab-Simulation Package for Ab-initio Real-space Calculations

M-SPARC is a Matlab code for performing ab initio Kohn–Sham Density Functional Theory simulations. Version 2.0.0 of the software further extends its capability to include relativistic effects, dispersion interactions, and advanced semilocal/nonlocal exchange–correlation functionals. These features significantly increase the fidelity of first principles calculations that can be performed using M-SPARC.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A numerical Poisson solver with improved radial solutions for a self-consistent locally scaled self-interaction correction method

Abstract The universal applicability of density functional approximations is limited by self-interaction error made by these functionals. Recently, a novel one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities was proposed. The locally scaled SIC (LSIC) method is exact for the one-electron densities, and unlike the well-known Perdew–Zunger SIC (PZSIC) method recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to PZSIC method as a special case when isoorbital indicator is set to the unity. Here, we present a numerical scheme that we have adopted to evaluate the Coulomb potential of the electron density scaled by the iso-orbital indicator required for the self-consistent LSIC calculations. After analyzing the behavior of the finite difference method (FDM) and the green function solution to the radial part of the Poisson equation, we adopt a hybrid approach that uses the FDM for the Coulomb potential due to the monopole and the GF for all higher-order terms. The performance of the resultant hybrid method is assessed using a variety of systems. The results show improved accuracy than earlier numerical schemes. We also find that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Dopant levels in large nanocrystals using stochastic optimally tuned range-separated hybrid density functional theory

In this work, we apply a stochastic version of an optimally tuned range-separated hybrid functional to provide insight on the electronic properties of P- and B- doped Si nanocrystals of experimentally relevant sizes. We show that we can use the range-separation parameter for undoped systems to calculate accurate results for dopant activation energies. We apply this strategy for tuning functionals to study doped nanocrystals up to 2.5 nm in diameter at the hybrid functional level. In this confinement regime, the P- and B- dopants have large activation energies and have strongly localized states that lie deep within the energy gaps. Structural relaxation plays a greater role for B-substituted dopants and contributes to the increase in activation energy when the B dopant is near the nanocrystal surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗