Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional approximation”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Quantum Monte Carlo analysis of a charge ordered insulating antiferromagnet: The Ti 4 O 7 Magneli phase

The Magneli phase Ti 4 O 7 is an important transition metal oxide with a wide range of applications because of its interplay between charge, spin, and lattice degrees of freedom. At low temperatures, it has non-trivial magnetic states very close in energy, driven by electronic exchange and correlation interactions. We have examined three low- lying states, one ferromagnetic and two antiferromagnetic, and calculated their energies as well as Ti spin moment distributions using highly accurate Quantum Monte Carlo methods. We compare our results to those obtained from density functional theory- based methods that include approximate corrections for exchange and correlation. Our results confirm the nature of the states and their ordering in energy, as compared with density-functional theory methods. However, the energy differences and spin distributions differ. Here, a detailed analysis suggests that non-local exchange-correlation functionals, in addition to other approximations such as LDA+U to account for correlations, are needed to simultaneously obtain better estimates for spin moments, distributions, energy differences and energy gaps.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Non-local thermodynamic equilibrium 1.5D modeling of red giant stars

Spectra for two-dimensional (2D) stars in the 1.5D approximation are created from synthetic spectra of one-dimensional (1D) non-local thermodynamic equilibrium (NLTE) spherical model atmospheres produced by the PHOENIX code. The 1.5D stars have the spatially averaged Rayleigh-Jeans flux of a K3-4 III star while varying the temperature difference between the two 1D component models (ΔT {sub 1.5D}) and the relative surface area covered. Synthetic observable quantities from the 1.5D stars are fitted with quantities from NLTE and local thermodynamic equilibrium (LTE) 1D models to assess the errors in inferred T {sub eff} values from assuming horizontal homogeneity and LTE. Five different quantities are fit to determine the T {sub eff} of the 1.5D stars: UBVRI photometric colors, absolute surface flux spectral energy distributions (SEDs), relative SEDs, continuum normalized spectra, and TiO band profiles. In all cases except the TiO band profiles, the inferred T {sub eff} value increases with increasing ΔT {sub 1.5D}. In all cases, the inferred T {sub eff} value from fitting 1D LTE quantities is higher than from fitting 1D NLTE quantities and is approximately constant as a function of ΔT {sub 1.5D} within each case. The difference between LTE and NLTE for the TiO bands is caused indirectly by the NLTE temperature structure of the upper atmosphere, as the bands are computed in LTE. We conclude that the difference between T {sub eff} values derived from NLTE and LTE modeling is relatively insensitive to the degree of the horizontal inhomogeneity of the star being modeled and largely depends on the observable quantity being fit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Pressure induced magnetic phase transition in RhFe{sub 3}N and IrFe{sub 3}N: An ab-initio study

The structural, electronic, elastic and magnetic properties of RhFe{sub 3}N and IrFe{sub 3}N are investigated using ab-initio calculations based on density functional theory as implemented in VASP code within the gradient generalized approximation. The non-spin polarized and spin polarized calculations are performed for these nitrides at normal and high pressures. It is found that these ternary nitrides are stable in ferromagnetic state at normal pressure. The lattice constant and bulk modulus values are calculated. The electronic structure reveals that these nitrides are metallic at normal pressure. The calculated elastic constants indicate that they are mechanically stable at ambient pressure. Ferromagnetic to nonmagnetic phase transition is observed in RhFe{sub 3}N and IrFe{sub 3}N at high pressure. Ferromagnetism is quenched in these nitrides at high pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Fractional charge and spin errors in self-consistent Green’s function theory

We examine fractional charge and spin errors in self-consistent Green’s function theory within a second-order approximation (GF2). For GF2, it is known that the summation of diagrams resulting from the self-consistent solution of the Dyson equation removes the divergences pathological to second-order Møller-Plesset (MP2) theory for strong correlations. In the language often used in density functional theory contexts, this means GF2 has a greatly reduced fractional spin error relative to MP2. The natural question then is what effect, if any, does the Dyson summation have on the fractional charge error in GF2? To this end, we generalize our previous implementation of GF2 to open-shell systems and analyze its fractional spin and charge errors. We find that like MP2, GF2 possesses only a very small fractional charge error, and consequently minimal many electron self-interaction error. This shows that GF2 improves on the critical failings of MP2, but without altering the positive features that make it desirable. Furthermore, we find that GF2 has both less fractional charge and fractional spin errors than typical hybrid density functionals as well as random phase approximation with exchange.

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

Goal-Oriented Optimal Approximations of Bayesian Linear Inverse Problems

We propose optimal dimensionality reduction techniques for the solution of goal-oriented linear-Gaussian inverse problems, where the quantity of interest (QoI) is a function of the inversion parameters. These approximations are suitable for large-scale applications. In particular, we study the approximation of the posterior covariance of the QoI as a low-rank negative update of its prior covariance and prove optimality of this update with respect to the natural geodesic distance on the manifold of symmetric positive definite matrices. Assuming exact knowledge of the posterior mean of the QoI, the optimality results extend to optimality in distribution with respect to the Kullback--Leibler divergence and the Hellinger distance between the associated distributions. We also propose the approximation of the posterior mean of the QoI as a low-rank linear function of the data and prove optimality of this approximation with respect to a weighted Bayes risk. Both of these optimal approximations avoid the explicit computation of the full posterior distribution of the parameters and instead focus on directions that are well informed by the data and relevant to the QoI. These directions stem from a balance among all the components of the goal-oriented inverse problem: prior information, forward model, measurement noise, and ultimate goals. We illustrate the theory using a high-dimensional inverse problem in heat transfer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Deep Learning without Global Optimization by Random Fourier Neural Networks

Here we introduce a new training algorithm for deep neural networks that utilize random complex exponential activation functions. Our approach employs a Markov chain Monte Carlo sampling procedure to iteratively train network layers, avoiding global and gradient-based optimization while maintaining error control. It consistently attains the theoretical approximation rate for residual networks with complex exponential activation functions, determined by network complexity. Additionally, it enables efficient learning of multiscale and high-frequency features, producing interpretable parameter distributions. Despite using sinusoidal basis functions, we do not observe Gibbs phenomena in approximating discontinuous target functions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2011

The Evolution of Ly-alpha Emitting Galaxies Between z = 2.1 and z = 3.l

We describe the results of a new, wide-field survey for z= 3.1 Ly-alpha emission-line galaxies (LAEs) in the Extended Chandra Deep Field South (ECDF-S). By using a nearly top-hat 5010 Angstrom filter and complementary broadband photometry from the MUSYC survey, we identify a complete sample of 141 objects with monochromatic fluxes brighter than 2.4E-17 ergs/cm^2/s and observers-frame equivalent widths greater than ~ 80 Angstroms (i.e., 20 Angstroms in the rest-frame of Ly-alpha). The bright-end of this dataset is dominated by x-ray sources and foreground objects with GALEX detections, but when these interlopers are removed, we are still left with a sample of 130 LAE candidates, 39 of which have spectroscopic confirmations. This sample overlaps the set of objects found in an earlier ECDF-S survey, but due to our filter's redder bandpass, it also includes 68 previously uncataloged sources. We confirm earlier measurements of the z=3.1 LAE emission-line luminosity function, and show that an apparent anti-correlation between equivalent width and continuum brightness is likely due to the effect of correlated errors in our heteroskedastic dataset. Finally, we compare the properties of z=3.1 LAEs to LAEs found at z=2.1. We show that in the approximately 1 Gyr after z approximately 3, the LAE luminosity function evolved significantly, with L * fading by approximately 0.4 mag, the number density of sources with L greater than 1.5E42 ergs/s declining by approximately 50%, and the equivalent width scalelength contracting from 70^{+7}_{-5} Angstroms to 50^{+9}_{-6} Angstroms. When combined with literature results, our observations demonstrate that over the redshift range z approximately 0 to z approximately 4, LAEs contain less than approximately 10% of the star-formation rate density of the universe.

Ciardullo, Robin↗
DOE OSTI2015

A comparative density functional study on electrical properties of layered penta-graphene

We present a comparative study of the influence of the number of layers, the biaxial strain in the range of −3% to 3%, and the stacking misalignments on the electronic properties of a new 2D carbon allotrope, penta-graphene (PG), based on hybrid-functional method within the density functional theory (DFT). In comparison with local exchange-correlation approximation in the DFT, the hybrid-functional provides an accurate description on the degree of p{sub z} orbitals localization and bandgap. Importantly, the predicted bandgap of few-layer PG has a weak layer dependence. The bandgap of monolayer PG is 3.27 eV, approximately equal to those of GaN and ZnO; and the bandgap of few-layer PG decreases slowly with the number of layers (N) and converge to 2.57 eV when N ≥ 4. Our calculations using HSE06 functional on few-layer PG reveal that bandgap engineering by stacking misalignment can further tune the bandgap down to 1.37 eV. Importantly, there is no direct-to-indirect bandgap transition in PG by varying strain, layer number, and stacking misalignment. Owing to its tunable, robustly direct, and wide bandgap characteristics, few-layer PG is promising for optoelectronic and photovoltaic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Water dissociation on Ni(100) and Ni(111): Effect of surface temperature on reactivity

Water adsorption and dissociation on Ni(100) and Ni(111) surfaces are studied using density functional theory calculations. Water adsorbs on top site on both the surfaces, while H and OH adsorb on four fold hollow and three fold hollow (fcc) sites on Ni(100) and Ni(111), respectively. Transition states (TS) on both surfaces are identified using climbing image-nudged elastic band method. It is found that the barrier to dissociation on Ni(100) surface is slightly lower than that on Ni(111) surface. Dissociation on both the surfaces is exothermic, while the exothermicity on Ni(100) is large. To study the effect of lattice motion on the energy barrier, TS calculations are performed for various values of Q (lattice atom coordinate along the surface normal) and the change in the barrier height and position is determined. Calculations show that the energy barrier to reaction decreases with increasing Q and increases with decreasing Q on both the surfaces. Dissociation probability values at different surface temperatures are computed using semi-classical approximation. Results show that the influence of surface temperature on dissociation probability on the Ni(100) is significantly larger compared to that of Ni(111). Moreover, on Ni(100), a dramatic shift in energy barrier to lower incident energy values is observed with increasing surface temperature, while the shift is smaller in the case of Ni(111)

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

Trapped energetic electron driven modes with second stable regime in tokamak plasmas

Features of fishbone-like trapped energetic electrons driven modes (TEEMs), which are a particular manifestation of the fishbone gap modes, are investigated taking into account model radial profile and the pitch angle distribution of the energetic electrons (EEs). The TEEMs are found unstable only when the beta values of EEs β{sub h} (=pressure of the energetic electrons/pressure of magnetic field) are higher than certain critical values which are proportional to perturbed energy δW{sup ^}{sub c} of background plasma and much lower than that for traditional fishbone modes. In addition, TEEMs become stable again and enter a second stable regime when β{sub h} values are higher than second critical values. Furthermore, the modes can only be excited in plasmas which are stable for MHD instabilities. The real frequency and growth rate of the modes are approximately linear and parabolic functions of β{sub h}, respectively. The real frequency is very low but not zero in the vicinity of lower beta region, whereas it is comparable to the toroidal precession frequency of the EEs in higher critical beta region. The numerical results show that the second stable regime is easy to form when q = 1 flux surface is close to the magnetic axis. Suitable density gradient of EEs and magnetic shear are other two necessary conditions for excitation of the TEEMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Design of Distributed Controllers Seeking Optimal Power Flow Solutions Under Communication Constraints

This paper focuses on power distribution networks featuring inverter-interfaced distributed energy resources (DERs), and develops feedback controllers that drive the DER output powers to solutions of time-varying AC optimal power flow (OPF) problems. Control synthesis is grounded on primal-dual-type methods for regularized Lagrangian functions, as well as linear approximations of the AC power-flow equations. Convergence and OPF-solution-tracking capabilities are established while acknowledging: i) communication-packet losses, and ii) partial updates of control signals. The latter case is particularly relevant since it enables asynchronous operation of the controllers where DER setpoints are updated at a fast time scale based on local voltage measurements, and information on the network state is utilized if and when available, based on communication constraints. As an application, the paper considers distribution systems with high photovoltaic integration, and demonstrates that the proposed framework provides fast voltage-regulation capabilities, while enabling the near real-time pursuit of solutions of AC OPF problems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Intensity Distribution of the Three-Wave Diffraction from Dislocation Epitaxial Layers in the Reciprocal Space

The three-wave X-ray diffraction in strongly disordered epitaxial layers of GaN and ZnO is experimentally investigated. The charts of the intensity distribution in the reciprocal space are plotted in coordinates q{sub θ} and q{sub ϕ} for the most intensive three-wave combination (1010)/(1011) by means of subsequent θ- and ϕ-scanning. A nontrivial shape of the θ-sections of these contours at a distance from the ϕ center of reflection is revealed; it is different for different samples. For the θ-curves at the center of reflection, we observed a common peak that may be approximated by the Voigt function with a power-low decrease in the intensity at the wings; the decrease law (from–4.5 to–5.0) is found to be considerably greater than that for the similar curves of two-wave diffraction and not depending on the dislocation density and distribution in layers. In some films we observed a coarse-block structure; in addition, it follows from the distribution in the reciprocal space that these blocks are turned with respect to each other around a normal to the surface, which allows us to suggest the existence of low-angle boundaries between them, consisting exclusively of edge dislocations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Quinone 1 e – and 2 e – /2 H + Reduction Potentials: Identification and Analysis of Deviations from Systematic Scaling Relationships

Quinones participate in diverse electron transfer and proton-coupled electron transfer processes in chemistry and biology. To understand the relationship between these redox processes, an experimental study was carried out to probe the 1 e – and 2 e – /2 H + reduction potentials of a number of common quinones. The results reveal a non-linear correlation between the 1 e – and 2 e – /2 H + reduction potentials. This unexpected observation prompted a computational study of 134 different quinones, probing their 1 e – reduction potentials, p K a values, and 2 e – /2 H + reduction potentials. We report the density functional theory calculations reveal an approximately linear correlation between these three properties and an effective Hammett constant associated with the quinone substituent(s). However, deviations from this linear scaling relationship are evident for quinones that feature intramolecular hydrogen bonding, halogen substituents, charged substituents, and/or sterically bulky substituents. These results, particularly the different substituent effects on the 1 e – versus 2 e – /2 H + reduction potentials, have important implications for designing quinones with tailored redox properties.

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

Exploring the Random Phase Approximately for materials chemistry and physics

This proposal focuses on improved accuracy for the delicate energy differences of interest in materials chemistry with the fully nonlocal random phase approximation (RPA) in a density functional context. Could RPA or RPA-like approaches become standard methods of first-principles electronic-structure calculation for atoms, molecules, solids, surfaces, and nano-structures? Direct RPA includes the full exact exchange energy and a nonlocal correlation energy from the occupied and unoccupied Kohn-Sham orbitals and orbital energies, with an approximate but universal description of long-range van der Waals attraction. RPA also improves upon simple pair-wise interaction potentials or vdW density functional theory. This improvement is essential to capture accurate energy differences in metals and different phases of semiconductors. The applications in this proposal are challenges for the simpler approximations of Kohn-Sham density functional theory, which are part of the current “standard model” for quantum chemistry and condensed matter physics. Within this project we already applied RPA on different structural phase transitions on semiconductors, metals and molecules. Although RPA predicts accurate structural parameters, RPA has proven not equally accurate in all kinds of structural phase transitions. Therefore a correction to RPA can be necessary in many cases. We are currently implementing and testing a nonempirical, spatially nonlocal, frequency-dependent model for the exchange-correlation kernel in the adiabatic-connection fluctuation-dissipation context. This kernel predicts a nearly-exact correlation energy for the electron gas of uniform density. If RPA or RPA-like approaches prove to be reliably accurate, then expected increases in computer power may make them standard in the electronic-structure calculations of the future.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Theoretical study of the relativistic molecular rotational g-tensor

An original formulation of the relativistic molecular rotational g-tensor valid for heavy atom containing compounds is presented. In such formulation, the relevant terms of a molecular Hamiltonian for non-relativistic nuclei and relativistic electrons in the laboratory system are considered. Terms linear and bilinear in the nuclear rotation angular momentum and an external uniform magnetic field are considered within first and second order (relativistic) perturbation theory to obtain the rotational g-tensor. Relativistic effects are further analyzed by carrying out the linear response within the elimination of the small component expansion. Quantitative results for model systems HX (X=F, Cl, Br, I), XF (X=Cl, Br, I), and YH{sup +} (Y=Ne, Ar, Kr, Xe, Rn) are obtained both at the RPA and density functional theory levels of approximation. Relativistic effects are shown to be small for this molecular property. The relation between the rotational g-tensor and susceptibility tensor which is valid in the non-relativistic theory does not hold within the relativistic framework, and differences between both molecular parameters are analyzed for the model systems under study. It is found that the non-relativistic relation remains valid within 2% even for the heavy HI, IF, and XeH{sup +} systems. Only for the sixth-row Rn atom a significant deviation of this relation is found.

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

Generation, dynamics, and correlations of the fission fragments' angular momenta

The generation of angular momentum in fissioning nuclei is not well understood. The predictions of different models disagree, particularly concerning the correlation between the fragments' angular momenta. In this article, a time-dependent collective Hamiltonian model is proposed to treat the generation of the angular momentum in the fission fragments due to the quantum uncertainty principle as well as the dynamics of the collective wave function during and after scission. The model is constructed in the framework of the frozen Hartree-Fock approximation using a Skyrme energy functional to extract deformations of the fission fragments as well as the interactions in a derived collective Hamiltonian. The fission reactions studied are $^{240}$Pu $\rightarrow$ $^{132}$Sn+$^{108}$Ru and $^{240}$Pu $\rightarrow$ $^{144}$Ba+$^{96}$Sr. The model can account for a large part of the angular momentum found in experimental data. The orientation of the angular momentum of each fragment is found to be mainly in the plane perpendicular to the fission axis, in agreement with the experiment. The magnitudes of the angular momenta in the two fragments are nearly uncorrelated, in agreement with the recent experimental data of Wilson et al., Nature (London) 590, 566 (2021). Some of the conclusions of the traditional collective vibration model are supported by the present model but some are not. Surprisingly, it is found that the angular momenta of the fragments are slightly correlated positively as in a wriggling mode. As a result, it is also found that the presence of an octupole deformation in a fragment can significantly increase the generated angular momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Inference in infinite-dimensional inverse problems - Discretization and duality

Many techniques for solving inverse problems involve approximating the unknown model, a function, by a finite-dimensional 'discretization' or parametric representation. The uncertainty in the computed solution is sometimes taken to be the uncertainty within the parametrization; this can result in unwarranted confidence. The theory of conjugate duality can overcome the limitations of discretization within the 'strict bounds' formalism, a technique for constructing confidence intervals for functionals of the unknown model incorporating certain types of prior information. The usual computational approach to strict bounds approximates the 'primal' problem in a way that the resulting confidence intervals are at most long enough to have the nominal coverage probability. There is another approach based on 'dual' optimization problems that gives confidence intervals with at least the nominal coverage probability. The pair of intervals derived by the two approaches bracket a correct confidence interval. The theory is illustrated with gravimetric, seismic, geomagnetic, and helioseismic problems and a numerical example in seismology.

Stark, Philip B.↗
DOE OSTI2014

Results for aliovalent doping of CeBr{sub 3} with Ca{sup 2+}

Despite the outstanding scintillation performance characteristics of cerium tribromide (CeBr{sub 3}) and cerium-activated lanthanum tribromide, their commercial availability and application are limited due to the difficulties of growing large, crack-free single crystals from these fragile materials. This investigation employed aliovalent doping to increase crystal strength while maintaining the optical properties of the crystal. One divalent dopant (Ca{sup 2+}) was used as a dopant to strengthen CeBr{sub 3} without negatively impacting scintillation performance. Ingots containing nominal concentrations of 1.9% of the Ca{sup 2+} dopant were grown, i.e., 1.9% of the CeBr{sub 3} molecules were replaced by CaBr{sub 2} molecules, to match our target replacement of 1 out of 54 cerium atoms be replaced by a calcium atom. Precisely the mixture was composed of 2.26 g of CaBr{sub 2} added to 222.14 g of CeBr{sub 3}. Preliminary scintillation measurements are presented for this aliovalently doped scintillator. Ca{sup 2+}-doped CeBr{sub 3} exhibited little or no change in the peak fluorescence emission for 371 nm optical excitation for CeBr{sub 3}. The structural, electronic, and optical properties of CeBr{sub 3} crystals were studied using the density functional theory within the generalized gradient approximation. Calculated lattice parameters are in agreement with the experimental data. The energy band structures and density of states were obtained. The optical properties of CeBr{sub 3}, including the dielectric function, were calculated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗