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,297 records · Page 72

Accurate atomic quantum defects from particle–particle random phase approximation

The accuracy of calculations of atomic Rydberg excitations cannot be judged by the usual measures, such as mean unsigned errors of many transitions. We show how to use quantum defect (QD) theory to (a) separate errors due to approximate ionisation potentials, (b) extract smooth QDs to compare with experiment, and (c) quantify those defects with a few characteristic parameters. The particle–particle random phase approximation (pp-RPA) produces excellent Rydberg transitions that are an order of magnitude more accurate than those of time-dependent density functional theory with standard approximations. Here, we even extract reasonably accurate defects from the lithium Rydberg series, despite the reference being open-shell. Our methodology can be applied to any Rydberg series of excitations with four transitions or more to extract the underlying threshold energy and characteristic QD parameters. Our pp-RPA results set a demanding challenge for other excitation methods to match.

74 ATOMIC AND MOLECULAR PHYSICS↗

Bosonic field digitization for quantum computers

Quantum simulation of quantum field theory is a flagship application of quantum computers that promises to deliver capabilities beyond classical computing. The realization of quantum advantage will require methods that can accurately predict error scaling as a function of the resolution and parameters of the model and that can be implemented efficiently on quantum hardware. In this paper, we address the representation of lattice bosonic fields in a discretized field amplitude basis, develop methods to predict error scaling, and present efficient qubit implementation strategies. A low-energy subspace of the bosonic Hilbert space, defined by a boson occupation number cutoff, can be represented with exponentially good accuracy by a low-energy subspace of a finite-size Hilbert space. The finite representation construction and the associated errors are directly related to the accuracy of the Nyquist-Shannon sampling and the finite Fourier transforms of the boson number states in the field and the conjugate-field bases. We analyze the relation between the boson mass, the discretization parameters used for wave function sampling, and the finite representation size. Numerical simulations of small size Φ 4 problems demonstrate that the boson mass optimizing the sampling of the ground state wave function is a good approximation to the optimal boson mass yielding the minimum low-energy subspace size. However, we find that accurate sampling of general wave functions does not necessarily result in accurate representation. Finally, we develop methods for validating and adjusting the discretization parameters to achieve more accurate simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

TU-F-CAMPUS-I-05: Parameterization of the Noise Power Spectrum in X-Ray Computed Tomography

Purpose: The purpose of this work was to develop a method so that the noise power spectrum (NPS) can be approximated for arbitrary levels of mAs, from a single determination in CT. Methods: The NPS is factorized into 2 components, 1) a parameterized function representing the 1D normalized spatial frequency distribution and 2) a function to scale the magnitude of 1) for arbitrary values of mAs. The 1D NPS, normalized by image variance (NNPS), was determined for 2 FBP reconstruction kernels (smoothing and edge enhancing) for 400 mAs. The NNPS were fit to the parameterized function and a scaling function was established to approximate the image variance at arbitrary values of mAs. Using the root mean square error normalized by the maximum value (NRMSE), the NPS approximated with the factorization method was compared to the NPS determined at 5 different mAs levels. Results: The factorization resulted in a set of 7 coefficients that can be used to approximate the 1D NPS, for arbitrary levels of mAs, for the convolution kernels studied in this work. The approximated NPS (factorization) agreed well with the determined NPS for all mAs levels. The greatest NRMSE was 0.02 and was observed for the edge enhancing kernel. Conclusion: The proposed factorization method has been demonstrated as applicable for FBP reconstruction. It can be used to approximate the 1D NPS for arbitrary levels of mAs, from a single NPS determination. Furthermore, approximations of the 1D NPS can conveniently be distributed since the factorization method only used 7 coefficients in the approximation.

60 APPLIED LIFE SCIENCES↗

Exchange-correlation magnetic fields in spin-density-functional theory

In spin-density-functional theory for noncollinear magnetic materials, the Kohn-Sham system features exchange-correlation (xc) scalar potentials and magnetic fields. The significance of the xc magnetic fields is not very well explored; in particular, they can give rise to local torques on the magnetization, which are absent in standard local and semilocal approximations. We obtain exact benchmark solutions for two electrons on four-site extended Hubbard lattices over a wide range of interaction strengths, and compare exact xc potentials and magnetic fields with approximations obtained from orbital-dependent xc functionals. Here, the xc magnetic fields turn out to play an increasingly important role as systems becomes more and more correlated and the electrons begin to localize; the effects of the xc torques, however, remain relatively minor. The approximate xc functionals perform overall quite well, but tend to favor symmetry-broken solutions for strong interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Multi-Resolution Nonlinear Mapping Technique for Design and Analysis Application

This report describes a nonlinear mapping technique where the unknown static or dynamic system is approximated by a sum of dimensionally increasing functions (one-dimensional curves, two-dimensional surfaces, etc.). These lower dimensional functions are synthesized from a set of multi-resolution basis functions, where the resolutions specify the level of details at which the nonlinear system is approximated. The basis functions also cause the parameter estimation step to become linear. This feature is taken advantage of to derive a systematic procedure to determine and eliminate basis functions that are less significant for the particular system under identification. The number of unknown parameters that must be estimated is thus reduced and compact models obtained. The lower dimensional functions (identified curves and surfaces) permit a kind of "visualization" into the complexity of the nonlinearity itself.

Phan, Minh Q.↗

A Multi-Resolution Nonlinear Mapping Technique for Design and Analysis Applications

This report describes a nonlinear mapping technique where the unknown static or dynamic system is approximated by a sum of dimensionally increasing functions (one-dimensional curves, two-dimensional surfaces, etc.). These lower dimensional functions are synthesized from a set of multi-resolution basis functions, where the resolutions specify the level of details at which the nonlinear system is approximated. The basis functions also cause the parameter estimation step to become linear. This feature is taken advantage of to derive a systematic procedure to determine and eliminate basis functions that are less significant for the particular system under identification. The number of unknown parameters that must be estimated is thus reduced and compact models obtained. The lower dimensional functions (identified curves and surfaces) permit a kind of "visualization" into the complexity of the nonlinearity itself.

Phan, Minh Q.↗

Kernel learning backward SDE filter for data assimilation

In this paper, we develop a kernel learning backward SDE filter method to estimate the state of a stochastic dynamical system based on its partial noisy observations. A system of forward backward stochastic differential equations is used to propagate the state of the target dynamical model, and Bayesian inference is applied to incorporate the observational information. Further, to characterize the dynamical model in the entire state space, we introduce a kernel learning method to learn a continuous global approximation for the conditional probability density function of the target state by using discrete approximated density values as training data. Numerical experiments demonstrate that the kernel learning backward SDE is highly effective.

97 MATHEMATICS AND COMPUTING↗

Time-dependent potential-functional embedding theory

We introduce a time-dependent potential-functional embedding theory (TD-PFET), in which atoms are grouped into subsystems. In TD-PFET, subsystems can be propagated by different suitable time-dependent quantum mechanical methods and their interactions can be treated in a seamless, first-principles manner. TD-PFET is formulated based on the time-dependent quantum mechanics variational principle. The action of the total quantum system is written as a functional of the time-dependent embedding potential, i.e., a potential-functional formulation. By exploiting the Runge-Gross theorem, we prove the uniqueness of the time-dependent embedding potential under the constraint that all subsystems share a common embedding potential. We derive the integral equation that such an embedding potential needs to satisfy. As proof-of-principle, we demonstrate TD-PFET for a Na{sub 4} cluster, in which each Na atom is treated as one subsystem and propagated by time-dependent Kohn-Sham density functional theory (TDDFT) using the adiabatic local density approximation (ALDA). Our results agree well with a direct TDDFT calculation on the whole Na{sub 4} cluster using ALDA. We envision that TD-PFET will ultimately be useful for studying ultrafast quantum dynamics in condensed matter, where key regions are solved by highly accurate time-dependent quantum mechanics methods, and unimportant regions are solved by faster, less accurate methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding ferromagnetism and optical absorption in 3d transition metal-doped cubic ZrO{sub 2} with the modified Becke-Johnson exchange-correlation functional

The electronic structure, magnetic, and optical properties in cubic crystalline phase of Zr{sub 1−x}TM{sub x}O{sub 2} (TM = V, Mn, Fe, and Co) at x = 6.25% are studied using density functional theory with the Generalized Gradient Approximation and the modified Becke-Johnson of the exchange-correlation energy and potential. In our calculations, the zirconia is a p-type semiconductor and has a large band gap. We evaluated the possibility of long-range magnetic order for transition metal ions substituting Zr. Our results show that ferromagnetism is the ground state in V, Mn, and Fe-doped ZrO{sub 2} and have a high value of energy in Mn-doped ZrO{sub 2}. However, in Co-doped ZrO{sub 2}, antiferromagnetic ordering is more stable than the ferromagnetic one. The exchange interaction mechanism has been discussed to explain the responsible of this stability. Moreover, it has been found that the V, Mn, and Fe transition metals provide half-metallic properties considered to be the leading cause, responsible for ferromagnetism. Furthermore, the optical absorption spectra in the TM -doped cubic ZrO{sub 2} are investigated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Should muffin tin radius vary in different structures of a material?: A case study

Quantum mechanical calculations based on density functional theory and a generalized gradient approximation (GGA) have been used to study the structural properties of YbN. Its predicted unit cell lattice parameter in NaCl (B1) structure is 4.7810Å and in CsCl (B2) structure it is 2.8685Å. In the determination of lattice parameter the muffin tin radius (R{sub MT}) of constituent atoms play important role. In both the structures the muffin tin radius for Yb and N converges to 2.3 and 1.4 a.u., respectively.

36 MATERIALS SCIENCE↗

Magnetism in undoped ZnS studied from density functional theory

The magnetic property induced by the native defects in ZnS bulk, thin film, and quantum dots are investigated comprehensively based on density functional theory within the generalized gradient approximation + Hubbard U (GGA + U) approach. We find the origin of magnetism is closely related to the introduction of hole into ZnS systems. The relative localization of S-3p orbitals is another key to resulting in unpaired p-electron, due to Hund's rule. For almost all the ZnS systems under study, the magnetic moment arises from the S-dangling bonds generated by Zn vacancies. The charge-neutral Zn vacancy, Zn vacancy in 1− charge sate, and S vacancy in the 1+ charge sate produce a local magnetic moment of 2.0, 1.0, and 1.0 μ{sub B}, respectively. The Zn vacancy in the neutral and 1− charge sates are the important cause for the ferromagnetism in ZnS bulk, with a Curie temperature (T{sub C}) above room temperature. For ZnS thin film with clean (111) surfaces, the spins on each surface are ferromagnetically coupled but antiferromagnetically coupled between two surfaces, which is attributable to the internal electric field between the two polar (111) surfaces of the thin film. Only surface Zn vacancies can yield local magnetic moment for ZnS thin film and quantum dot, which is ascribed to the surface effect. Interactions between magnetic moments on S-3p states induced by hole-doping are responsible for the ferromagnetism observed experimentally in various ZnS samples.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Explicit approximations to estimate the perturbative diffusivity in the presence of convectivity and damping. III. Cylindrical approximations for heat waves traveling inwards

In this paper, a number of new explicit approximations are introduced to estimate the perturbative diffusivity (χ), convectivity (V), and damping (τ) in cylindrical geometry. For this purpose, the harmonic components of heat waves induced by localized deposition of modulated power are used. The approximations are based on the heat equation in cylindrical geometry using the symmetry (Neumann) boundary condition at the plasma center. This means that the approximations derived here should be used only to estimate transport coefficients between the plasma center and the off-axis perturbative source. If the effect of cylindrical geometry is small, it is also possible to use semi-infinite domain approximations presented in Part I and Part II of this series. A number of new approximations are derived in this part, Part III, based upon continued fractions of the modified Bessel function of the first kind and the confluent hypergeometric function of the first kind. These approximations together with the approximations based on semi-infinite domains are compared for heat waves traveling towards the center. The relative error for the different derived approximations is presented for different values of the frequency, transport coefficients, and dimensionless radius. Moreover, it is shown how combinations of different explicit formulas can be used to estimate the transport coefficients over a large parameter range for cases without convection and damping, cases with damping only, and cases with convection and damping. The relative error between the approximation and its underlying model is below 2% for the case, where only diffusivity and damping are considered. If also convectivity is considered, the diffusivity can be estimated well in a large region, but there is also a large region in which no suitable approximation is found. This paper is the third part (Part III) of a series of three papers. In Part I, the semi-infinite slab approximations have been treated. In Part II, cylindrical approximations are treated for heat waves traveling towards the plasma edge assuming a semi-infinite domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding density functional theory (DFT) and completing it in practice

We review some salient points in the derivation of density functional theory (DFT) and of the local density approximation (LDA) of it. We then articulate an understanding of DFT and LDA that seems to be ignored in the literature. We note the well-established failures of many DFT and LDA calculations to reproduce the measured energy gaps of finite systems and band gaps of semiconductors and insulators. We then illustrate significant differences between the results from self consistent calculations using single trial basis sets and those from computations following the Bagayoko, Zhao, and Williams (BZW) method, as enhanced by Ekuma and Franklin (BZW-EF). Unlike the former, the latter calculations verifiably attain the absolute minima of the occupied energies, as required by DFT. These minima are one of the reasons for the agreement between their results and corresponding, experimental ones for the band gap and a host of other properties. Further, we note predictions of DFT BZW-EF calculations that have been confirmed by experiment. Our subsequent description of the BZW-EF method ends with the application of the Rayleigh theorem in the selection, among the several calculations the method requires, of the one whose results have a full, physics content ascribed to DFT. This application of the Rayleigh theorem adds to or completes DFT, in practice, to preserve the physical content of unoccupied, low energy levels. Discussions, including implications of the method, and a short conclusion follow the description of the method. The successive augmentation of the basis set in the BZW-EF method, needed for the application of the Rayleigh theorem, is also necessary in the search for the absolute minima of the occupied energies, in practice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Interdependence of spin structure, anion height and electronic structure of BaFe{sub 2}As{sub 2}

Superconducting as well as other electronic properties of Fe-based superconductors are quite sensitive to the structural parameters specially, on anion height which is intimately related to z{sub As}, the fractional z co-ordinate of As atom. Due to presence of strong magnetic fluctuation in these Fe-based superconductors, optimized structural parameters (lattice parameters a, b, c) including z{sub As} using density functional theory (DFT) under generalized gradient approximation (GGA) does not match experimental values accurately. In this work, we show that the optimized value of z{sub As} is strongly influenced by the spin structures in the orthorhombic phase of BaFe{sub 2}As{sub 2} system. We take all possible spin structures for the orthorhombic BaFe{sub 2}As{sub 2} system and then optimize z{sub As}. Using these optimized structures we calculate electronic structures like density of states, band structures etc., for each spin configurations. From these studies we show that the electronic structure, orbital order which is responsible for structural as well as related to nematic transition, are significantly influenced by the spin structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exact solutions to the fermion propagator Schwinger-Dyson equation in Minkowski space with on-shell renormalization for quenched QED

With the introduction of a spectral representation, the Schwinger-Dyson equation (SDE) for the fermion propagator is formulated in Minkowski space in QED. After imposing the on-shell renormalization conditions, analytic solutions for the fermion propagator spectral functions are obtained in four dimensions with a renormalizable version of the Gauge Technique anzatz for the fermion-photon vertex in the quenched approximation in the Landau gauge. Despite the limitations of this model, having an explicit solution provides a guiding example of the fermion propagator with the correct analytic structure. The Padé approximation for the spectral functions is also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MAFT: The Multicomputer Architecture for Fault-Tolerance

Multicomputer Architecture for Fault-Tolerance (MAFT) is a loosely coupled multiprocessor system designed to achieve an unreliability of less than 10(exp -10)/hr in flight-critical real time applications. The MAFT design objectives and architecture are presented. The fault-tolerance implementation of major functions in MAFT is also presented, including communication; task scheduling; reconfiguration; clock synchronization; and data handling and voting. The need for Byzantine agreement or approximate agreement in various functions is discussed. Different methods were selected to achieve agreement in various subsystems. These methods are illustrated by a more detailed description of the task scheduling and error handling subsystems.

Kieckhafer, Roger M.↗

Time-independent hybrid enrichment for finite element solution of transient conduction–radiation in diffusive grey media

We investigate the effectiveness of the partition-of-unity finite element method for transient conduction–radiation problems in diffusive grey media. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary diffusion approximation to the radiation in grey media. The coupled equations are integrated in time using a semi-implicit method in the finite element framework. We show that for the considered problems, a combination of hyperbolic and exponential enrichment functions based on an approximation of the boundary layer leads to improved accuracy compared to the conventional finite element method. It is illustrated that this approach can be more efficient than using h adaptivity to increase the accuracy of the finite element method near the boundary walls. The performance of the proposed partition-of-unity method is analyzed on several test examples for transient conduction–radiation problems in two space dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗