Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Nonstationary envelope process and first excursion probability

A definition of the envelope of nonstationary random processes is proposed. The establishment of the envelope definition makes it possible to simulate the nonstationary random envelope directly. Envelope statistics, such as the density function, joint density function, moment function, and level crossing rate, which are relevent to analyses of catastrophic failure, fatigue, and crack propagation in structures, are derived. Applications of the envelope statistics to the prediction of structural reliability under random loadings are discussed in detail.

Yang, J.↗
DOE OSTI2021

Point defect induced incommensurate dipole moments in the K Ca 2 Nb 3 O 10 Dion-Jacobson layered perovskite

Local structural distortions due to isolated atomic defects and defect complexes strongly affect the macroscopic properties of oxide ceramics. While the characterization of local defect structures is more common in simple AB O 3 perovskites, unambiguous determination of the same in layered perovskites is more difficult due to their complex crystal structures. In this study, we combined x-ray pair distribution function and density functional theory calculations to characterize the structure of cation-oxygen divacancy pairs in a Dion-Jacobson (D-J) layered perovskite. Our results indicate that local incommensurate dipole moments with polarization density in the range of ~ 0.1–17μC/cm 2 are created due to divacancy-induced structural distortions in the D-J phase K Ca 2 Nb 3 O 10 . This is comparable with defect dipole moments observed in well-known perovskite ferroelectrics. The current results imply that controlling the atomic defects can potentially lead to significant control of dielectric properties in D-J layered perovskites.

3-dimensional systems↗
DOE OSTI2021

Nuclear charge densities in spherical and deformed nuclei: Toward precise calculations of charge radii

Precise measurements of atomic transitions affected by electron-nucleus hyperfine interactions offer sensitivity to explore basic properties of the atomic nucleus and study fundamental symmetries, including the search for new physics beyond the standard model of particle physics. Specifically, such measurements, augmented by atomic and nuclear calculations, will enable extraction of the higher-order radial moments of the charge-density distribution in spherical and deformed nuclei. The new data impose higher precision requirements on a theoretical description. The nuclear charge density is composed of the proton point distribution folded with the nucleonic charge distributions. The latter induce subtle relativistic corrections due to the coupling of nucleon magnetic moments with the nuclear spin-orbit density. Additional corrections come from the effect of center-of-mass projection. We assess the precision of nuclear charge density calculations by studying the behavior of relativistic and center-of-mass motion corrections to the second and fourth charge radial moments. Special attention has been paid to the magnetic spin-orbit density associated with the local variations of the spin-orbit current. The calculations for semimagic and open-shell nuclei are performed in the framework of self-consistent mean-field theory using quantified energy density functionals and density-dependent pairing forces. We used the general expression for the spin-orbit form factor that is valid for spherical and deformed nuclei. We studied the impact of various correction terms on the charge radii, fourth radial moments, diffraction radii, and surface thickness of spherical and deformed nuclei. The spin-orbit corrections to charge radial moments and surface thickness show strong shell fluctuations which can make an appreciable effect when aiming at high-precision predictions of isotopic shifts. The inclusion of relativistic and center-of-mass corrections impacts the quality of energy density functionals optimized to charge radii data. To establish reliable constraints on the existence of new forces from isotope shift measurements, precise calculations of nuclear charge densities of deformed nuclei are needed. The proper inclusion of the spin-orbit charge density and other correction terms is essential when aiming at extraction of subtle effects which become particularly visible in isotopic trends. It is also important when developing high-quality nuclear energy density functionals optimized using heterogeneous datasets involving absolute charge radii, differential charge radii, and charge form factor properties deduced from electron-scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Computationally Guided Investigation of the Optical Spectra of Pure β-UO 3

Single-phase β-UO 3 is synthesized by flash heating UO 2 (NO 3 )·6H 2 O in air to 450 °C and annealing for 60 h under the same conditions. For the first time, we report the Raman spectra of pure β-UO 3 . To facilitate the assignment of Raman and infrared vibrational modes, we use density functional theory with density functional perturbation theory. By employing a novel analysis scheme that includes the mode frequencies as well as a quantitative analysis of the mode eigenvectors, we assign the observed spectral features to individual chemical modes. In particular, the density functional theory optimized structure, observed Raman spectrum, and eigenvector analysis suggest the presence of four crystallographically distinct uranyl ions, one more than has previously been suggested.

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

Analysis of defects in In 2 O 3 :H synthesized in presence of water vapor and hydrogen gas mixture

Using hybrid functional-based density functional theory calculations, we analyze the structure and kinetics of defects formed in two competing synthesis routes to prepare hydrogen-doped In 2 O 3 films, using a hydrogen and oxygen gas mixture vs synthesis in the presence of water vapor. For both of these synthesis routes, we find that H + is the dominant defect species: when the Fermi level is close to the conduction band, H + has a lower formation energy than other intrinsic or extrinsic defects. Our results also suggest that water molecules spontaneously split into H + (which occupies octahedral voids) and OH – interstitials (which occupies vacant oxygen lattice sites or oxygen vacancies). From the analysis of the binding energies between these different defects, we conclude that these defects do not cluster and are most likely to stay spatially distributed throughout the films. In addition, the sum of formation energies of an oxygen (i.e., O$^{2–}_i$) and a H + interstitial is close to the formation energy of a OH – interstitial, meaning that water molecules are completely split into 2H + and O$^{2–}_i$ at the synthesis conditions. Further, in the presence of H 2 + O 2 gas mixture, oxygen interstitials occupy oxygen vacancies while hydrogen interstitials occupy vacant oxygen lattice sites and form bonds with lattice oxygens. Our analysis of the defect equilibria suggests that the hydrogen content in films synthesized in the presence of water vapor is higher than films synthesized in the presence of a hydrogen gas mixture. At high dopant concentrations, a hydrogen bond network is formed in the system and this leads to large distortions in the lattice.

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

The role of water vapor during the synthesis of hydrogen doped In 2 O 3

Hydrogen-doped In 2 O 3 synthesized using magnetron sputtering or atomic layer deposition in the presence of water vapor has high transparency in the near infra-red region and mobility values that are almost three to four times of those in commercially used In 2 O 3 :Sn at a much lower carrier concentration. However, simple questions like what happens when water molecules enter In 2 O 3 or why the dominant charge carrier changes from an oxygen vacancy in In 2 O 3 :Sn to H + in In 2 O 3 are not clear. Using hybrid functional based density functional theory calculations, in this study we show that water molecules spontaneously split into H + and OH – . The H + interstitial bonds with an oxygen in In 2 O 3 , and the OH – interstitial occupies a vacant lattice site or is trapped by an oxygen vacancy. Thus, in agreement with experimental findings, our results suggest that the concentration of oxygen vacancies decreases as the partial pressure of water increases. The charge transition level, ϵ ( – / + ) , of a hydrogen interstitial and the Fermi level at which the formation energies of OH – and H + are equal lie above the conduction band minimum, suggesting that the H + interstitial has a lower formation energy. In addition, water molecules do not exhibit any propensity to form clusters, meaning that H + ions can diffuse and remain spatially separated from each other.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculating and Characterizing the Charge Distributions in Solids

Accurate estimation of the partial atomic charges on metal centers is useful for understanding electronic and catalytic properties of materials. However, different methods of calculating these charges may give quite different results; this issue has been more widely studied for molecules than for solids. Here we study the charges on the metal centers of a test set of 18 solids containing transition metals by using density functional theory with several density functionals (PBE, PBE+U, TPSS, revTPSS, HLE17, revM06-L, B3LYP, B3LYP*, and other exchange-modified B3LYP functionals) and four charge models (Bader, Hirshfeld, CM5, and DDEC6). The test set contains 12 systems with nonmagnetic metal centers (eight metal oxides (MO 2 ), two metal sulfides (MS 2 ), and two metal selenides (MSe 2 )) and six ferromagnetic transition metal complexes. Our study shows that, among the four types of charges, Bader charges are the highest and Hirshfeld charges are the lowest for all the systems, regardless of the functional being used. The CM5 charges are bigger than DDEC6 charges for MX 2 with M = Ti or Mo and X = S or Se, but for the other 14 cases they are lower. We found that the most of the systems are sensitive to the Hubbard U parameters in PBE+U and to the percentage X of Hartree–Fock exchange in exchange-modified B3LYP; as we increase U or X, the charges on the metal atoms in MX 2 increase steadily. Testing different density functionals shows charges calculated with higher Hubbard U parameters in PBE+U are comparable to B3LYP (with 20% Hartree–Fock exchange). Among four meta-GGA functionals studied, the charges with HLE17 have the closest agreement with B3LYP. In conclusion, the variation of charges with choice of charge model is greater than the variation with choice of density functional.

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

Effects of baryonic feedback on the cosmic web

Here, we study the effect of baryons on the cosmic web—halos, filaments, walls, and voids. To do so, we apply a modified version of NEXUS, a cosmic web morphological analysis algorithm, to the IllustrisTNG simulations. We find that halos lose more than 10% of their mass due to baryons, mostly to filaments and a small portion to walls and voids. However, the mass transfer does not significantly shift the boundaries of structures, leaving the volume fractions of the cosmic structures largely unaffected. We quantify the effects of baryonic feedback on the power spectrum and the probability density function of the density field for individual cosmic structures. For the power spectrum, most suppression due to feedback can be accounted for by including M ≥ 10 12 M ⊙ /h halos, without considering other cosmic structures. However, when examining the probability density function of the density field, we find nearly 100% suppression of the emptiest regions and 10%-level effects (boost or suppression) in the remaining regions of filaments, walls, and voids. Our results indicate the importance of modeling the effects of baryons in the whole cosmic web, not just halos, for cosmological analysis beyond two-point statistics or field-based inferences.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Improved first-principles equation-of-state table of deuterium for high-energy-density applications

Here, we present a first-principles equation-of-state (EOS) table of deuterium aimed at improving the previously established first-principles equation-of-state table (FPEOS). The EOS table presented here, referred to as iFPEOS, introduces (1) a universal density functional theory (DFT) treatment of all density and temperature conditions, (2) a fully consistent treatment of exchange-correlation (XC) thermal effects across the entire range of temperatures covered, and (3) quantum treatment of ions. Based on ab initio molecular dynamics driven by thermal density functional theory, iFPEOS includes density points in the range 1×10 –3 $\leqslant$ ρ $\leqslant$ 1.6×10 3 g/cm 3 and temperature points in the range 800 K $\leqslant$ T $\leqslant$ 256 MK, thus covering the challenging warm dense matter (WDM) regime. For an improved description of the electronic structure, iFPEOS employs an advanced free-energy XC density functional with explicit temperature dependence, which is at the metageneralized gradient approximation level of DFT. We use the latest orbital-free free-energy density functional for the high-temperature regime where it shows excellent agreement with standard Mermin-Kohn-Sham DFT. For quantum treatment of ions we use path-integral molecular dynamics in order to take into account nuclear quantum effects. Results are compared to other EOS models and most recent experimental measurements of deuterium properties such as the molecular-toatomic fluid transition, the principal and reshock Hugoniot, and sound speed. We find that iFPEOS provides an improved agreement with experimental data compared to other first-principles EOS models in the WDM regime for pressures up to 200 GPa and temperatures up to 60 000 K. For higher pressures and temperatures, however, iFPEOS is in agreement with other models in predicting lower compressibility and higher sound speed along the Hugoniot, compared to experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗
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↗
DOE OSTI2021

Audacity of huge: overcoming challenges of data scarcity and data quality for machine learning in computational materials discovery

Machine learning (ML)-accelerated discovery requires large amounts of high-fidelity data to reveal predictive structure–property relationships. For many properties of interest in materials discovery, the challenging nature and high cost of data generation has resulted in a data landscape that is both scarcely populated and of dubious quality. Data-driven techniques starting to overcome these limitations include the use of consensus across functionals in density functional theory, the development of new functionals or accelerated electronic structure theories, and the detection of where computationally demanding methods are most necessary. When properties cannot be reliably simulated, large experimental data sets can be used to train ML models. In the absence of manual curation, increasingly sophisticated natural language processing and automated image analysis are making it possible to learn structure–property relationships from the literature. Finally, models trained on these data sets will improve as they incorporate community feedback.

36 MATERIALS SCIENCE↗
DOE OSTI2020

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Gaussian windows: A tool for exploring multivariate data

Presented here is a method for interactively exploring a large set of quantitative multivariate data, in order to estimate the shape of the underlying density function. It is assumed that the density function is more or less smooth, but no other specific assumptions are made concerning its structure. The local structure of the data in a given region may be examined by viewing the data through a Gaussian window, whose location and shape are chosen by the user. A Gaussian window is defined by giving each data point a weight based on a multivariate Gaussian function. The weighted sample mean and sample covariance matrix are then computed, using the weights attached to the data points. These quantities are used to compute an estimate of the shape of the density function in the window region. The local structure of the data is described by a method similar to the method of principal components. By taking many such local views of the data, we can form an idea of the structure of the data set. The method is applicable in any number of dimensions. The method can be used to find and describe simple structural features such as peaks, valleys, and saddle points in the density function, and also extended structures in higher dimensions. With some practice, we can apply our geometrical intuition to these structural features in any number of dimensions, so that we can think about and describe the structure of the data. Since the computations involved are relatively simple, the method can easily be implemented on a small computer.

Jaeckel, Louis A.↗
DOE OSTI2022

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

Density functional theory calculations play a central role in understanding chemical and solid-state systems. Progress depends on density functionals that accurately reproduce both energies, for thermochemistry, and properly describe ground states and other properties that are of interest. The Cr dimer, benzene and graphene are especially important benchmark systems for quantum chemistry and condensed matter physics. The Strongly Constrained and Appropriately Normed (SCAN) functional, which is an advanced meta-generalized gradient approximation functional that significantly improves molecular energies is shown to perform poorly for the Cr dimer. This is connected with its poor performance for itinerant solid-state magnets and is a consequence of over localization of electrons, thus illustrating an analogy between the Cr dimer and itinerant magnets. The Cr dimer is a notoriously difficult system for density functionals. However, we additionally find that SCAN predicts an incorrect symmetry broken ground state for 2D graphene and for the benzene molecule, which is surprising considering that ground states of these are known to be well described even by the simplest local density approximation. We show that SCAN overly favors localized spin polarized states, which is a serious deficiency of this approach. Thus, the challenge of finding density functionals that accurately treat both localized and delocalized electronic systems remains.

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

ABINIT: Overview and focus on selected capabilities

ABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperature-dependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. ABINIT has strong links with many other software projects that are briefly mentioned.

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

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, which is a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Further, our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗