Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Gradient corrections to the exchange-correlation free energy

We develop the first-order gradient correction to the exchange-correlation free energy of the homogeneous electron gas for use in finite-temperature density functional calculations. Based on this, we propose and implement a simple temperature-dependent extension for functionals beyond the local density approximation. These finite-temperature functionals show improvement over zero-temperature functionals, as compared to path-integral Monte Carlo calculations for deuterium equations of state, and perform without computational cost increase compared to zero-temperature functionals and so should be used for finite-temperature calculations. Furthermore, while the present functionals are valid at all temperatures including zero, non-negligible difference with zero-temperature functionals begins at temperatures above 10 000 K.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Ab initio calculations of polarization, piezoelectric constants, and elastic constants of InAs and InP in the wurtzite phase

We report first-principle density functional calculations of the spontaneous polarization, piezoelectric stress constants, and elastic constants for the III–V wurtzite structure semiconductors InAs and InP. Using the density functional theory implemented in the VASP code, we obtain polarization values–0.011 and–0.013 C/m{sup 2}, and piezoelectric constants e{sub 33} (e{sub 31}) equal to 0.091 (–0.026) and 0.012 (–0.081) C/m{sup 2} for structurally relaxed InP and InAs respectively. These values are consistently smaller than those of nitrides. Therefore, we predict a smaller built-in electric field in such structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Accuracy of Partial Core Corrections Using Fourier Transforms in Pseudopotential–Density Functional Theory

Partial core corrections can be important in obtaining an accurate description of nonlinear exchange-correlation functionals and improving the transferability of pseudopotentials. We show that a widely used procedure, which calculates partial core charge density, $ρ^{partial}_{core}$, in Fourier space and then converts it to real space with fast Fourier transforms, can lead to sizable numerical errors of exchange-correlation potentials in the vacuum region. Such errors occur in modeling low-dimensional materials or surfaces with supercells. The loss of accuracy originates from the slow-decaying feature of core charge density in reciprocal space. Numerical errors on the order of 1 eV in the Kohn–Sham energies of unoccupied states can occur in pseudopotential–density functional calculations. The direct calculation of the partial core charge in real space can avoid the numerical errors caused by Fourier transforms.

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

Interaction energy and isosteric heat of adsorption between hydrogen and magnesium diboride

Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).

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

First-principles study of electronic states in LiBe{sub 2}

By coupling structure prediction methods with first-principles total energy calculations we find lattice constant and bulk modulus of hcp LiBe{sub 2}. Structural and electronic properties are investigated using Linear Combination of Atomic Orbital’s method. The generalized gradient approximation based density functional calculations are attempted to look for the ground state. Thereafter ground state electronic band structure, total density of states and anisotropies in directional Compton profiles are computed. Step like behavior seen broadly in density of states (DOS) is appreciated by band structure suggesting 2D features. Anisotropic behavior of electronic states around Fermi level is also reflected in anisotropies of directional Compton profiles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Hybrid Inorganic-Organic Materials with an Optoelectronically Active Aromatic Cation: (C 7 H 7 ) 2 SnI 6 and C 7 H 7 PbI 3

Inorganic materials with organic constituents—hybrid materials—have shown incredible promise as chemically tunable functional materials with interesting optical and electronic properties. Here, the preparation and structure are reported of two hybrid materials containing the optoelectronically active tropylium ion within tin- and lead-iodide inorganic frameworks with distinct topologies. The crystal structures of tropylium tin iodide, (C 7 H 7 ) 2 SnI 6 , and tropylium lead iodide, C 7 H 7 PbI 3 , were solved using high-resolution synchrotron powder X-ray diffraction informed by X-ray pair distribution function data and high-resolution time-of-flight neutron diffraction. Tropylium tin iodide contains isolated tin(IV)-iodide octahedra and crystallizes as a deep black solid, while tropylium lead iodide presents one-dimensional chains of face-sharing lead(II)-iodide octahedra and crystallizes as a bright red-orange powder. Experimental diffuse reflectance spectra are in good agreement with density functional calculations of the electronic structure. Calculations of the band decomposed charge densities suggest that the deep black color of tropylium tin iodide is attributed to iodide ligand to tin metal charge transfer, while the bright red-orange color of tropylium lead iodide arises from charge transfer between iodine and tropylium states. Understanding the origins of the observed optoelectronic properties of these two compounds, with respect to their distinct topologies and organic–inorganic interactions, provides insight into the design of tropylium-containing compounds for potential optical and electronic applications.

Maughan, Annalise E.↗
DOE OSTI2017

Ion-gel-gating-induced oxygen vacancy formation in epitaxial L a 0.5 S r 0.5 Co O 3 – δ films from in operando x-ray and neutron scattering

Ionic-liquid/gel-based transistors have emerged as a potentially ideal means to accumulate high charge-carrier densities at the surfaces of materials such as oxides, enabling control over electronic phase transitions. Substantial gaps remain in the understanding of gating mechanisms, however, particularly with respect to charge carrier vs oxygen defect creation, one contributing factor being the dearth of experimental probes beyond electronic transport. Here we demonstrate the use of synchrotron hard x-ray diffraction and polarized neutron reflectometry as in operando probes of ion-gel transistors based on ferromagnetic La 0.5 Sr 0.5 CoO 3-δ . An asymmetric gate-bias response is confirmed to derive from electrostatic hole accumulation at negative gate bias vs oxygen vacancy formation at positive bias. The latter is detected via a large gate-induced lattice expansion (up to 1%), complementary bulk measurements and density functional calculations enabling quantification of the bias-dependent oxygen vacancy density. Remarkably, the gate-induced oxygen vacancies proliferate through the entire thickness of 30-40-unit-cell-thick films, quantitatively accounting for changes in the magnetization depth profile. In conclusion, these results directly elucidate the issue of electrostatic vs redox-based response in electrolyte-gated oxides, also demonstrating powerful approaches to their in operando investigation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GDB-9-Ex_TD-DFT-PBE0: Dataset containing Time Dependent Density Functional Theory (TDDFT) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 96,766 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Time Dependent Density Functional Theory (TDDFT) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The time-dependent density-functional theory (TDDFT) approach with the PBE0 exchange-correlation functional and ORCAs default integration grid was employed. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2021

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
DOE OSTI2016

Efficient Use of an Adapting Database of Ab Initio Calculations To Generate Accurate Newtonian Dynamics

We develop and demonstrate a method to efficiently use density functional calculations to drive classical dynamics of complex atomic and molecular systems. The method has the potential to scale to systems and time scales unreachable with current ab initio molecular dynamics schemes. It relies on an adapting dataset of independently computed Hellmann–Feynman forces for atomic configurations endowed with a distance metric. The metric on configurations enables fast database lookup and robust interpolation of the stored forces. Here, we discuss mechanisms for the database to adapt to the needs of the evolving dynamics, while maintaining accuracy, and other extensions of the basic algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography

Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

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

Electronic and thermoelectric properties of InN studied using ab initio density functional theory and Boltzmann transport calculations

The thermoelectric properties of indium nitride in the most stable wurtzite phase (w-InN) as a function of electron and hole concentrations and temperature were studied by solving the semiclassical Boltzmann transport equations in conjunction with ab initio electronic structure calculations, within Density Functional Theory. Based on maximally localized Wannier function basis set and the ab initio band energies, results for the Seebeck coefficient are presented and compared with available experimental data for n-type as well as p-type systems. Also, theoretical results for electric conductivity and power factor are presented. Most cases showed good agreement between the calculated properties and experimental data for w-InN unintentionally and p-type doped with magnesium. Our predictions for temperature and concentration dependences of electrical conductivity and power factor revealed a promising use of InN for intermediate and high temperature thermoelectric applications. The rigid band approach and constant scattering time approximation were utilized in the calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Atomic Structure and Electronic Properties of Anionic Germanium–Zirconium Clusters

This paper presents spatial structure optimization results and calculated electronic spectra for ZrGe{sub n}{sup −} (n = 8–20) anion clusters. Comparison of density-functional calculation results and available experimental data allows us to identify real spatial structures of the clusters. The formation of stable endohedral ZrGe{sub n}{sup −} clusters is possible for n ≥ 12. The clusters with a smaller number of germanium atoms predominantly have exohedral structures.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Understanding the Photocatalytic Properties of Pt/CeO x /TiO 2 : Structural Effects on Electronic and Optical Properties

Ceria-titania interfaces play a crucial role in different chemical processes but are especially promising for the photocatalytic splitting of water using light in the visible wavelength region when Pt is added to the system. Yet, the complexity of this hierarchical structure hampers the study of the origin of its outstanding properties. In this article, the structural, electronic and optoelectronic properties of CeO 2 /TiO 2 systems containing 1D, 2D, and 3D particles of ceria are analyzed by means of density functional calculations. Adsorption sites and vacancy effects have been studied to model Pt adsorption. Density of states calculations and absorption spectra simulations explain the behavior of these systems. Conclusively, these models are used for the screening of other metals that can be combined with this heterostructure to possibly find more efficient water splitting photocatalysts.

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

Impact of excess iron on the calculated electronic and magnetic properties of Co{sub 3–x}Fe{sub x}Si Heusler-compound

Ab initio density functional calculations are performed on the Co{sub 3–x}Fe{sub x}Si alloys with variable iron composition (1≤x≤3). The evolution of the structural, electronic, and magnetic properties is investigated within different levels of approximations. These alloys crystallize in cubic Heusler structures, which evolve from the regular L2₁ structure for Co₂FeSi to the inverse X structure for x≥2. Using on-site Coulomb interactions of U{sub eff}(Co)=3.07 eV and U{sub eff}(Fe)=3.4 eV is found to describe consistently the experimental properties for x≤2. A good agreement between calculated and experimental magnetic moments is found for the cubic inverse Heusler phases without the addition of Hubbard-model. The spin polarization is found to range from 1 for Fe concentrations of x≤2, indicating the half-metallic character, to –0.29 at x=3.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Reproducibility of calculations on Li species with correlation-consistent basis sets

Here we have found that two different basis set series, both labeled cc-pVxZ, are widely available and widely used for calculation on Li compounds such that calculations that should be the same are different, causing a reproducibility problem. Here we study the magnitude of the disagreement of these nominally-the-same basis sets for wave function and density functional calculations of Li clusters. We find deviations among the basis sets of as large as 22 millihartrees. In most cases one obtains lower energies with the noncanonical choices. Therefore, it is important to specify which version of the basis set is used for calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗