Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional approximations”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

The muscle spindle as a feedback element in muscle control

The muscle spindle, the feedback element in the myotatic (stretch) reflex, is a major contributor to muscular control. Therefore, an accurate description of behavior of the muscle spindle during active contraction of the muscle, as well as during passive stretch, is essential to the understanding of muscle control. Animal experiments were performed in order to obtain the data necessary to model the muscle spindle. Spectral density functions were used to identify a linear approximation of the two types of nerve endings from the spindle. A model reference adaptive control system was used on a hybrid computer to optimize the anatomically defined lumped parameter estimate of the spindle. The derived nonlinear model accurately predicts the behavior of the muscle spindle both during active discharge and during its silent period. This model is used to determine the mechanism employed to control muscle movement.

Andrews, L. T.↗
DOE OSTI2016

Effects of the preparation processes on structural, electronic, and optical properties of LaHoO{sub 3}

Highlights: • The preparation of LaHoO{sub 3} by combustion and hydrothermal techniques was performed. • The crystal structure and characterization of the synthesized oxides were carried out. • The lattice distortions have effects on Raman, luminescence and optical band gap. • Electronic properties and band structure have been carried out by DFT method. • Sample prepared by hydrothermal process gives higher emission intensities. - ABSTRACT: Lanthanum holmium oxide (LaHoO{sub 3}) was synthesized using combustion and hydrothermal processes in order to compare the effects of preparation method on structural, electronic, and optical properties. For this purpose, the as-prepared samples were characterized with infrared spectroscopy, Raman scattering, X-ray diffraction, and scanning electronic microscopy. The Rietveld refinements structure of as-prepared samples show that lattice distortions and volume expansion were obtained for LaHoO{sub 3} synthesized with hydrothermal method compared to combustion process. The luminescence properties were investigated upon 488 and 750 nm laser excitations. The emission curves show the presence of characteristic Ho{sup 3+} transitions with a higher emission intensity for a sample prepared by hydrothermal method than the one obtained by combustion process. The structures properties effects on the band structure and the density of states were realized via generalized gradient approximation functional with Hubbard correction. The lattice distortions were found to have effects on Raman, luminescence and optical band gap.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
DOE OSTI2016

DFT calculations of electronic and optical properties of SrS with LDA, GGA and mGGA functionals

The theoretical investigations of electronic and optical properties of SrS are made using the first principle DFT calculations. The calculations are performed for the local-density approximation (LDA), generalized gradient approximation (GGA) and for an alternative form of GGA i.e. metaGGA for both rock salt type (B1, Fm3m) and cesium chloride (B2, Pm3m) structures. The band structure, density of states and optical spectra are calculated under various available functional. The calculations with LDA and GGA functional underestimate the values of band gaps with all functional, however the values with mGGA show reasonably good agreement with experimental and those calculated by using other methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Cross correlations from back reaction on stochastic magnetic fields

The induction equation induces non trivial correlations between the primordial curvature mode and the magnetic mode which is a non linear effect. Assuming a stochastic, gaussian magnetic field the resulting power spectra determining the two point cross correlation functions between the primordial curvature perturbation and the magnetic energy density contrast as well as the magnetic anisotropic stress are calculated approximately. The corresponding numerical solutions are used to calculate the angular power spectra determining the temperature anisotropies and polarization of the cosmic microwave background, C{sub l}. It is found that the resulting C{sub l} are sub-leading in comparison to those generated by the compensated mode for a magnetic field which only redshifts with the expansion of the universe. The main focus are scalar modes, however, vector modes will also be briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Approximate DFT-based methods for generating diabatic states and calculating electronic couplings: models of two and more states

Four types of density functional theory (DFT)-based approaches are assessed in this work for the approximate construction of diabatic states and the evaluation of electronic couplings between these states. These approaches include the constrained DFT (CDFT) method, the constrained noninteracting electron (CNE) model to post-process Kohn–Sham operators, the approximate block-diagonalization (BD) of the Kohn–Sham operators, and the generalized Mulliken–Hush method. It is shown that the first three approaches provide a good description for long-distance intermolecular electron transfer (ET) reactions. On the other hand, inconsistent results were found when applying these approaches to intramolecular ET in strongly coupled, mixed-valence systems. Model analysis shows that this discrepancy is caused by the inappropriate use of the two-state model rather than the defects of the approaches themselves. The situation is much improved when more states are included in the model electronic Hamiltonian. The CNE and BD approaches can thus serve as efficient and robust alternatives for building ET models based on DFT calculations.

Yang, Chou-Hsun↗
DOE OSTI2014

Derivative couplings between TDDFT excited states obtained by direct differentiation in the Tamm-Dancoff approximation

Working within the Tamm-Dancoff approximation, we calculate the derivative couplings between time-dependent density-functional theory excited states by assuming that the Kohn-Sham superposition of singly excited determinants represents a true electronic wavefunction. All Pulay terms are included in our derivative coupling expression. The reasonability of our approach can be established by noting that, for closely separated electronic states in the infinite basis limit, our final expression agrees exactly with the Chernyak-Mukamel expression (with transition densities from response theory). Finally, we also validate our approach empirically by analyzing the behavior of the derivative couplings around the T{sub 1}/T{sub 2} conical intersection of benzaldehyde.

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

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

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

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

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

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗
DOE OSTI2013

Time-dependent density functional theory quantum transport simulation in non-orthogonal basis

Basing on the earlier works on the hierarchical equations of motion for quantum transport, we present in this paper a first principles scheme for time-dependent quantum transport by combining time-dependent density functional theory (TDDFT) and Keldysh's non-equilibrium Green's function formalism. This scheme is beyond the wide band limit approximation and is directly applicable to the case of non-orthogonal basis without the need of basis transformation. The overlap between the basis in the lead and the device region is treated properly by including it in the self-energy and it can be shown that this approach is equivalent to a lead-device orthogonalization. This scheme has been implemented at both TDDFT and density functional tight-binding level. Simulation results are presented to demonstrate our method and comparison with wide band limit approximation is made. Finally, the sparsity of the matrices and computational complexity of this method are analyzed.

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

Tunable noninteracting free-energy density functionals for high-energy-density physics applications

In this work, we introduce the concept of a tunable noninteracting free-energy density functional and present two examples realized: (i) via a simple one-parameter convex combination of two existing functionals and (ii) via the construction of a generalized gradient approximation (GGA) enhancement factor that contains one free parameter and is designed to satisfy a set of incorporated constraints. Functional (i), constructed as a combination of the local Thomas–Fermi and a pseudopotential-adapted GGA for the noninteracting free-energy, has already demonstrated its practical usability for establishing the high temperature end of the equation of state of deuterium [Phys. Rev. B 104, 144104 (2021)] and CHON resin [Phys. Rev. E 106, 045207 (2022)] for inertial confinement fusion applications. Hugoniot calculations for liquid deuterium are given as another example of how the application of computationally efficient orbital-free density functional theory (OF-DFT) can be utilized with the employment of the developed functionals. Once the functionals have been tuned such that the OF-DFT Hugoniot calculation matches the Kohn–Sham solution at some low-temperature point, agreement with the reference Kohn–Sham results for the rest of the high temperature Hugoniot path is very good with relative errors for compression and pressure on the order of 2% or less.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

First-principles studies of the structural, electronic, and optical properties of a novel thorium compound Rb{sub 2}Th{sub 7}Se{sub 15}

The structural, electronic, and optical properties of a recently synthesized thorium compound Rb{sub 2}Th{sub 7}Se{sub 15} have been calculated in the density functional theory framework for the first time. The calculated direct band gap was 1.471 eV (generalized gradient approximation) and 1.171 eV (local density approximation), with both results being close to the experimental result of 1.83 eV. High covalency/iconicity of the Th–Se/Rb–Se bonds was demonstrated by calculating effective Mulliken charges of all ions. The polarized calculations of the complex dielectric function are presented; dependence of the calculated index of refraction was fitted to the Sellmeyer equation in the wavelength range from 500 to 2500 nm. - Graphical abstract: Calculated band structure of Rb{sub 2}Th{sub 7}Se{sub 15}. - Highlights: • The first theoretical analysis of the Rb{sub 2}Th{sub 7}Se{sub 15} properties is reported. • Structural, electronic and optical properties of Rb{sub 2}Th{sub 7}Se{sub 15} were calculated. • An indirect character of Rb{sub 2}Th{sub 7}Se{sub 15} band gap was confirmed. • Dependence of the refractive index on the wavelength was calculated.

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

Cosmological parameters and evolution of the galaxy luminosity function

The relationship between the observed distribution of discrete sources of a flux limited sample, the luminosity function of these sources, and the cosmological model is discussed. It is stressed that some assumptions about the form and evolution of the luminosity function must be made in order to determine the cosmological parameters from the observed distribution of sources. Presented is a method to test the validity of these assumptions using the observations. It is shown how, using higher moments of the observed distribution, one can determine, independently of the cosmological model, all parameters of the luminosity function except those describing evolution of the density and the luminosity of the luminosity function. These methods are applied to the sample of approximately 1000 galaxies recently used by Loh and Spillar to determine a value of the cosmological density parameter Omega approx = 1. It is shown that the assumptions made by Loh and Spillar about the luminosity function are inconsistent with the data, and that a self-consistent treatment of the data indicates a lower value of Omega approx = 0.2 and a flatter luminosity function. It should be noted, however, that incompleteness in the sample could cause a flattening of the luminosity function and lower the calculated value of Omega and that uncertainty in the values of these parameters due to random fluctuations is large.

Caditz, David↗
NASA NTRS1989

Cosmological parameters and evolution of the galaxy luminosity function

The relationship between the observed distribution of discrete sources of a flux limited sample, the luminosity function of these sources, and the cosmological model is discussed. It is stressed that some assumptions about the form and evolution of the luminosity function must be made in order to determine the cosmological parameters from the observed distribution of sources. Presented is a method to test the validity of these assumptions using the observations. It is shown how, using higher moments of the observed distribution, one can determine, independently of the cosmological model, all parameters of the luminosity function except those describing evolution of the density and the luminosity of the luminosity function. These methods are applied to the sample of approximately 1000 galaxies recently used by Loh and Spillar to determine a value of the cosmological density parameter Omega approx = 1. It is shown that the assumptions made by Loh and Spillar about the luminosity function are inconsistent with the data, and that a self-consistent treatment of the data indicates a lower value of Omega approx = 0.2 and a flatter luminosity function. It should be noted, however, that incompleteness in the sample could cause a flattening of the luminosity function and lower the calculated value of Omega and that uncertainty in the values of these parameters due to random fluctuations is large.

Caditz, David↗
DOE OSTI2014

Exchange integrals in magnetoelectric hexagonal ferrite (SrCo{sub 2}Ti{sub 2}Fe{sub 8}O{sub 19}): A density functional study

The exchange integrals in magnetoelectric hexagonal ferrite SrCo{sub 2}Ti{sub 2}Fe{sub 8}O{sub 19} have been calculated by using density functional theory. To get 10 inter-sublattice and 3 intra-sublattice exchange integrals, the electronic structures and total energies of 20 spin arrangements have been calculated with General Gradient Approximation (GGA) + U method. The dependence of exchange integrals on U has been studied. The comparison between the exchange integrals in SrFe{sub 12}O{sub 19} and those in SrCo{sub 2}Ti{sub 2}Fe{sub 8}O{sub 19} shows that substitution of Co and Ti decreases the most interactions involving the 12 k sites. The investigation based on our exchange integrals indicates that magnetic interaction between R and S blocks reduces significantly in SrCo{sub 2}Ti{sub 2}Fe{sub 8}O{sub 19}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

First principle investigation of structural and electronic properties of bulk ZnSe

Electronic and structural properties of ZnSe are investigated using plane-wave self-consistent field method within the framework of density functional theory. The pseudopotential method within the local density approximation is used for the exchange-correlation potential. The equilibrium lattice parameter, static bulk modulus and its pressure derivative are calculated. The electronic band structure, partial density of states and density of states are also obtained. The results are compared with available theoretical calculations and experimental results.

36 MATERIALS SCIENCE↗