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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Towards quantum turbulence in cold atomic fermionic superfluids

Fermionic superfluids provide a new realization of quantum turbulence, accessible to both experiment and theory, yet relevant to phenomena from both cold atoms to nuclear astrophysics. In particular, the strongly interacting Fermi gas realized in cold-atom experiments is closely related to dilute neutron matter in neutron star crusts. Unlike the liquid superfluids 4 He (bosons) and 3 He (fermions), where quantum turbulence has been studied in laboratory for decades, superfluid Fermi gases stand apart for a number of reasons. They admit a rather reliable theoretical description based on density functional theory called the time-dependent superfluid local density approximation that describes both static and dynamic phenomena. Cold atom experiments demonstrate exquisite control over particle number, spin polarization, density, temperature, and interaction strength. Topological defects such as domain walls and quantized vortices, which lie at the heart of quantum turbulence, can be created and manipulated with time-dependent external potentials, and agree with the time-dependent theoretical techniques. While similar experimental and theoretical control exists for weakly interacting Bose gases, the unitary Fermi gas is strongly interacting. The resulting vortex line density is extremely high, and quantum turbulence may thus be realized in small systems where classical turbulence is suppressed. Fermi gases also permit the study of exotic superfluid phenomena such as the Larkin–Ovchinnikov–Fulde–Ferrell pairing mechanism for polarized superfluids which may give rise to 3D supersolids, and a pseudo-gap at finite temperatures that might affect the regime of classical turbulence. The dynamics associated with these phenomena has only started to be explored. Finally, superfluid mixtures have recently been realized, providing experimental access to phenomena like Andreev–Bashkin entrainment predicted decades ago. Superfluid Fermi gases thus provide a rich forum for addressing phenomena related to quantum turbulence with applications ranging from terrestrial superfluidity to astrophysical dynamics in neutron stars.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Detecting deviations from Gaussianity in high-redshift CMB lensing maps

ABSTRACT While the probability density function of the cosmic microwave background (CMB) convergence field approximately follows a Gaussian distribution, primordial non-Gaussianities, and small contributions from structures at low redshifts make the overall distribution slightly non-Gaussian. Some of the late-time component can be modelled using the distribution of galaxies and subtracted off from the original CMB lensing map to produce a map of matter distribution at high redshifts. Using this high-redshift mass map, we are able to directly study the early phases of structure formation. In this work, we forecast the detectability of signatures of non-Gaussianity due to non-linear structure formation at z > 1.2. Assuming the optimal case of no systematics, we find that it is challenging to detect such signatures in current surveys, but future experiments such as the deep field of CMB-S4 will be able to make detections of ∼7σ.

Astronomy & Astrophysics↗
DOE OSTI2016

Nonlocal energy-optimized kernel: Recovering second-order exchange in the homogeneous electron gas

In order to remedy some of the shortcomings of the random phase approximation (RPA) within adiabatic connection fluctuation-dissipation (ACFD) density functional theory, we introduce a short-ranged, exchange-like kernel that is one-electron self-correlation free and exact for two-electron systems in the high-density limit. By tuning a free parameter in our model to recover an exact limit of the homogeneous electron gas correlation energy, we obtain a nonlocal, energy-optimized kernel that reduces the errors of RPA for both homogeneous and inhomogeneous solids. Using wave-vector symmetrization for the kernel, we also implement RPA renormalized perturbation theory for extended systems, and demonstrate its capability to describe the dominant correlation effects with a low-order expansion in both metallic and nonmetallic systems. Further, the comparison of ACFD structural properties with experiment is also shown to be limited by the choice of norm-conserving pseudopotential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Quantum phase transition and ferromagnetism in C o 1 + x Sn

The onset of ferromagnetism in cobalt-tin alloys is investigated experimentally and theoretically. The Co 1+x Sn alloys were prepared by rapid quenching from the melt and form a modified hexagonal NiAs-type crystal structure for 0.45 ≤ x ≤ 1. The magnetic behavior is described analytically and by density-functional theory using supercells and the coherent-potential approximation. The excess of Co concentration x, which enters the interstitial 2d sites in the hypothetical NiAs-ordered parent alloy CoSn, yields a Griffiths-like phase and, above a quantum critical point (x c ≈ 0.65), a quantum phase transition to ferromagnetic order. Quantum critical exponents are determined on the paramagnetic and ferromagnetic sides of the transition and related to the nature of the magnetism in itinerant systems with different types of chemical disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Evaluating transition metal oxides within DFT-SCAN and SCAN + U frameworks for solar thermochemical applications

Using the strongly constrained and appropriately normed (SCAN) and SCAN+U approximations for describing electron exchange-correlation (XC) within density functional theory, we investigate the oxidation energetics, lattice constants, and electronic structure of binary Ce-, Mn-, and Fe-oxides, which are crucial ingredients for generating renewable fuels using two-step, oxide-based, solar thermochemical reactors. Unlike other common XC functionals, we find that SCAN does not over-bind the O 2 molecule, based on direct calculations of its bond energy and robust agreement between calculated formation enthalpies of main group oxides versus experiments. However, in the case of transition-metal oxides (TMOs), SCAN systematically overestimates (i.e., yields too negative) oxidation enthalpies due to remaining self-interaction errors in the description of their ground-state electronic structure. Adding a Hubbard U term to the transition-metal centers, where the magnitude of U is determined from experimental oxidation enthalpies, significantly improves the qualitative agreement and marginally improves the quantitative agreement of SCAN+U-calculated electronic structure and lattice parameters, respectively, with experiments. Importantly, SCAN predicts the wrong ground-state structure for a few oxides, namely, Ce 2 O 3 , Mn 2 O 3 , and Fe 3 O 4 , while SCAN+U predicts the right polymorph for all systems considered in this work. Furthermore, the SCAN+U framework, with an appropriately determined U, will be required to accurately describe ground-state properties and yield qualitatively consistent electronic properties for most transition-metal and rare earth oxides.

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

Negative thermal expansion in the Ruddlesden-Popper calcium titanates

Materials exhibiting negative thermal expansion (NTE) are important for the fabrication and operation of microelectronic devices and optical systems. In this paper, combining density-functional-theory calculations and the self-consistent quasiharmonic approximation method, we investigate the variation of the NTE mechanisms active in layered calcium titanate oxides. We find a NTE trend with layer number for the orthorhombic materials comprising the RP series, but the monoclinic polymorph is an outlier.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Surface fitting with biharmonic and harmonic models

The physical and geometric interpretation of a surface fitting technique, multiquadric equations (MQ), is discussed. The relation of this method and its reciprocal form to biharmonic and harmonic techniques respectively, is pointed out. In both forms the solutions may be viewed as being numerical approximations of an integral equation in which an unknown density function is the physical source for disturbing potential (reciprocal-MQ) or elastic displacement (MQ).

Hardy, R. L.↗
NASA NTRS1982

A numerical study of the effects of ambipolar diffusion on the collapse of magnetic gas clouds

The gravitational collapse of isothermal, nonrotating magnetic gas clouds have been calculated numerically, including the effects of ambipolar diffusion. The fractional ionization in the clouds is approximated by a power-law function of the gas density, f = K/n to the q-power, where K and q are adjustable parameters. Eleven numerical experiments were run, and the results indicate that the asymptotic character of collapse is determined mainly by the value of q and is largely independent of the other parameters characterizing a cloud (e.g., K, cloud mass). In particular, there is nearly a one-to-one correspondence between q and the slope, x, of the central magnetic field strength-gas density relationship. If q is no more than 0.8, a cloud collapses asymptotically, as though the magnetic field were 'frozen' to the neutral matter. The magnetic field strength at the center of a collapsing cloud is strongly amplified during collapse even for values of q of about 1, despite extremely low values of fractional ionization. A discussion of the theoretical basis for this unexpected behavior is given. Possible implications of our results for the problems of magnetic braking of rotating protostars and star formation in general are also presented.

Black, D. C.↗
NASA NTRS1990

Equilibrium models of self-gravitating inviscid disks resulting from the collapse of rotating clouds

A method is presented for the construction of velocity and surface density profiles of infinitely thin, self-gravitating disks formed from the collapse of spherical clouds, under the condition that the angular momentum of each material parcel is conserved. Results are presented for initially uniformly rotating clouds with density distributions varying as R exp -n and n = 0, 1, 2. The resulting disk surface density distributions are well represented by power laws in the cylindrical radius over most of their extent. The form of the zeroth-order approximation provides a useful representation of the final surface density function near the center for all three cases, and for n = 2 this result holds throughout the disk. This result is not generally true for the velocity functions, however.

Stemwedel, Sally W.↗
NASA NTRS1991

Center for Modeling of Turbulence and Transition (CMOTT). Research briefs: 1990

Brief progress reports of the Center for Modeling of Turbulence and Transition (CMOTT) research staff from May 1990 to May 1991 are given. The objectives of the CMOTT are to develop, validate, and implement the models for turbulence and boundary layer transition in the practical engineering flows. The flows of interest are three dimensional, incompressible, and compressible flows with chemistry. The schemes being studied include the two-equation and algebraic Reynolds stress models, the full Reynolds stress (or second moment closure) models, the probability density function models, the Renormalization Group Theory (RNG) and Interaction Approximation (DIA), the Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS).

Povinelli, Louis A.↗
NASA NTRS1994

Center for modeling of turbulence and transition: Research briefs, 1993

This research brief contains the progress reports of the research staff of the Center for Modeling of Turbulence and Transition (CMOTT) from June 1992 to July 1993. It is also an annual report to the Institute for Computational Mechanics in Propulsion located at Ohio Aerospace Institute and NASA Lewis Research Center. The main objectives of the research activities at CMOTT are to develop, validate, and implement turbulence and transition models for flows of interest in propulsion systems. Currently, our research covers eddy viscosity one- and two-equation models, Reynolds-stress algebraic equation models, Reynolds-stress transport equation models, nonequilibrium multiple-scale models, bypass transition models, joint scalar probability density function models, and Renormalization Group Theory and Direct Interaction Approximation methods. Some numerical simulations (LES and DNS) have also been carried out to support the development of turbulence modeling. Last year was CMOTT's third year in operation. During this period, in addition to the above mentioned research, CMOTT has also hosted the following programs: an eighteen-hour short course on 'Turbulence--Fundamentals and Computational Modeling (Part I)' given by CMOTT at the NASA Lewis Research Center; a productive summer visitor research program that has generated many encouraging results; collaborative programs with industry customers to help improve their turbulent flow calculations for propulsion system designs; a biweekly CMOTT seminar series with speakers from within and without the NASA Lewis Research Center including foreign speakers. In addition, CMOTT members have been actively involved in the national and international turbulence research activities. The current CMOTT roster and organization are listed in Appendix A. Listed in Appendix B are the abstracts of the biweekly CMOTT seminar. Appendix C lists the papers contributed by CMOTT members.

Liou, William W.↗
DOE OSTI2013

Charge asymmetry in rovibrationally excited HD{sup +} determined using explicitly correlated all-particle Gaussian functions

Very accurate non-Born-Oppenheimer quantum-mechanical calculations are performed to determine the average values of the interparticle distances and the proton-deuteron density function for the rovibrationally excited HD{sup +} ion. The states corresponding to excitations to all bound vibrational states (v = 0, …, 22) and simultaneously excited to the first excited rotational state (N = 1) are considered. To describe each state up to 8000 explicitly correlated all-particle Gaussian functions are used. The nonlinear parameters of the Gaussians are variationally optimized using a procedure that employs the analytical energy gradient determined with respect to these parameters. The results show an increasing asymmetry in the electron distribution with the vibrational excitation as the electron density shifts towards deuteron and away from the proton.

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

Affordable and accurate large-scale hybrid-functional calculations on GPU-accelerated supercomputers

Performing high accuracy hybrid functional calculations for condensed matter systems containing a large number of atoms is at present computationally very demanding or even out of reach if high quality basis sets are used. We present a highly optimized multiple graphics processing unit implementation of the exact exchange operator which allows one to perform fast hybrid functional density-functional theory (DFT) calculations with systematic basis sets without additional approximations for up to a thousand atoms. With this method hybrid DFT calculations of high quality become accessible on state-of-the-art supercomputers within a time-to-solution that is of the same order of magnitude as traditional semilocal-GGA functionals. The method is implemented in a portable open-source library.

density functional theory↗
NASA NTRS2004

Electron Density and Two-Channel Neutron Emission Measurements in Steady-State Spherical Inertial-Electrostatically Confined Plasmas, with Review of the 1-D Kinetic Model

Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the 1-D Poisson equation for EC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with non-fusing hydrogen-1. The deuterium rates are consistent with predictions from the model.

Dobson, Chris C.↗
DOE OSTI2022

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations

The thermal neutron scattering law (TSL) for crystalline beryllium hydride (BeH 2 ) is developed from first-principles ab initio lattice dynamics calculations and the impact of neutron thermalization in this material on critical mass is estimated. BeH 2 has a body-centered orthorhombic crystal structure with 12 molecules per unit cell and a theoretical density of 0.755 g/cm 3 . The vibrational (phonon) densities of states for H and Be bound in BeH 2 are determined using VASP density functional theory and PHONON lattice dynamics calculations. The TSLs for H bound in BeH 2 , H(BeH 2 ), and Be bound in BeH 2 , Be(BeH 2 ), are then evaluated in the incoherent approximation from the calculated H and Be partial phonon density of states using FLASSH. Finally, critical mass as a function of 235 U loading density for bare and reflected BeH 2 moderated spheres is predicted from MC21 Monte Carlo neutron transport calculations using ENDF/B-VIII.0 cross sections and the H(BeH 2 ) and Be(BeH 2 ) TSL evaluations. Comparisons are made to water (H 2 O), polyethylene (CH 2 ), and beryllium oxide (BeO) as moderators. These critical mass predictions are a refinement upon the prior work by Rao and Srinivasan that neglected thermal neutron scattering effects. The minimum critical mass of a BeH 2 moderated assembly is estimated to be 0.207 kg 235 U for a 0.20 m thick BeO reflected sphere and 0.178 kg 235 U for a 0.40 m thick BeO reflected sphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Progress towards understanding ultranonlocality through the wave-vector and frequency dependence of approximate exchange-correlation kernels

In the framework of time-dependent density functional theory (TDDFT), the exact exchange-correlation (xc) kernel f xc (n, q, ω) determines the ground-state energy, excited-state energies, lifetimes, and the time-dependent linear density response of any many-electron system. The recently developed MCP07 xc kernel f xc (n, q, ω) of A. Ruzsinszky et al. [Phys. Rev. B 101, 245135 (2020)] yields excellent uniform electron gas (UEG) ground-state energies and plausible plasmon lifetimes. As MCP07 is constructed to describe f xc of the UEG, it cannot capture optical properties of real materials. To verify this claim, we follow Nazarov et al. [Phys. Rev. Lett. 102, 113001 (2009)] to construct the long-range, dynamic xc kernel, lim q→0 f xc (n, q, ω) = -α(ω)e 2 /q 2 , of a weakly inhomogeneous electron gas, using MCP07 and other common xc kernels. The strong wavevector and frequency dependence of the “ultranonlocality” coefficient α(ω) is demonstrated for a variety of simple metals and semiconductors. We examine how imposing exact constraints on an approximate kernel shapes α(ω). Comparisons to kernels derived from correlated-wavefunction calculations are drawn.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Electron Density and Two-Channel Neutron Emission Measurements in Steady-State Spherical Inertial-Electrostatically Confined Plasmas, with Review of the One-Dimensional Kinetic Model

Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7 x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the one-dimensional Poisson equation for IEC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with nonfusing hydrogen- 1. The deuterium rates are consistent with predictions from the model.

Dobson, Chris C.↗