Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function 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 775 records · Page 43

Solution of the Problem of Free Vibrations of a Nonthin Orthotropic Shallow Shell of Variable Thickness in the Refined Statement

We consider the problem of investigation of the spectrum of natural vibrations of a nonthin orthotropic shallow shell variable in two coordinate directions of thickness in the nonclassical statement. The approach to the solution of the obtained two-dimensional boundary-value problem is based on its reduction (by the method of spline-approximation of the unknown functions along one coordinate direction) to the one-dimensional problem with its subsequent solution. We study different cases of boundary conditions imposed on the contours of the shell. We also perform the comparison and analysis of the natural frequencies and modes of vibrations of orthotropic shells of constant and variable thickness.

42 ENGINEERING↗

Approaching periodic systems in ensemble density functional theory via finite one-dimensional models

Ensemble density functional theory (EDFT) is a generalization of ground-state DFT, which is based on an exact formal theory of finite collections of a system's ground and excited states. EDFT in various forms has been shown to improve the accuracy of calculated energy level differences in isolated model systems, atoms, and molecules, but it is not yet clear how EDFT could be used to calculate band gaps for periodic systems. We extend the application of EDFT toward periodic systems by estimating the thermodynamic limit with increasingly large finite one-dimensional 'particle in a box' systems, which approach the uniform electron gas (UEG). Using ensemble-generalized Hartree and local spin density approximation exchange-correlation functionals, we find that corrections go to zero in the infinite limit, as expected for a metallic system. However, there is a correction to the effective mass, with results comparable to other calculations on 1D, 2D, and 3D UEGs, which indicates promise for non-trivial results from EDFT on periodic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Lineshape theory of pigment-protein complexes: How the finite relaxation time of nuclei influences the exciton relaxation-induced lifetime broadening

In pigment-protein complexes, often the excited states are partially delocalized and the exciton-vibrational coupling in the basis of delocalized states contains large diagonal and small off-diagonal elements. This inequality may be used to introduce potential energy surfaces (PESs) of exciton states and to treat the inter-PES coupling in Markov and secular approximations. The resulting lineshape function consists of a Lorentzian peak that is broadened by the finite lifetime of the exciton states caused by the inter-PES coupling and a vibrational sideband that results from the mutual displacement of the excitonic PESs with respect to that of the ground state. So far analytical expressions have been derived that relate the exciton relaxation-induced lifetime broadening to the Redfield [T. Renger and R. A. Marcus, J. Chem. Phys. 116, 9997 (2002)] or modified Redfield [M. Schröder, U. Kleinekathöfer, and M. Schreiber, J. Chem. Phys. 124, 084903 (2006)] rate constants of exciton relaxation, assuming that intra-PES nuclear relaxation is fast compared to inter-PES transfer. Here, we go beyond this approximation and provide an analytical expression, termed Non-equilibrium Modified Redfield (NeMoR) theory, for the lifetime broadening that takes into account the finite nuclear relaxation time. In an application of the theory to molecular dimers, we find that, for a widely used experimental spectral density of the exciton-vibrational coupling of pigment-protein complexes, the NeMoR spectrum at low-temperatures (T < 150 K) is better approximated by Redfield than by modified Redfield theory. At room temperature, the lifetime broadening obtained with Redfield theory underestimates the NeMoR broadening, whereas modified Redfield theory overestimates it by a similar amount. A fortuitous error compensation in Redfield theory is found to explain the good performance of this theory at low temperatures. Since steady state spectra of PPCs are often measured at low temperatures, Redfield theory still provides a numerically efficient alternative to NeMoR theory. At higher temperatures, we suggest to use NeMoR theory, because it has the same numerical costs as modified Redfield theory, but is more accurate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine-learning Kohn–Sham potential from dynamics in time-dependent Kohn–Sham systems

Abstract The construction of a better exchange-correlation potential in time-dependent density functional theory (TDDFT) can improve the accuracy of TDDFT calculations and provide more accurate predictions of the properties of many-electron systems. Here, we propose a machine learning method to develop the energy functional and the Kohn–Sham potential of a time-dependent Kohn–Sham (TDKS) system is proposed. The method is based on the dynamics of the Kohn–Sham system and does not require any data on the exact Kohn–Sham potential for training the model. We demonstrate the results of our method with a 1D harmonic oscillator example and a 1D two-electron example. We show that the machine-learned Kohn–Sham potential matches the exact Kohn–Sham potential in the absence of memory effect. Our method can still capture the dynamics of the Kohn–Sham system in the presence of memory effects. The machine learning method developed in this article provides insight into making better approximations of the energy functional and the Kohn–Sham potential in the TDKS system.

97 MATHEMATICS AND COMPUTING↗

Galaxy and Mass Assembly (GAMA): Exploring the WISE Web in G12

We present an analysis of the mid-infrared Wide-field Infrared Survey Explorer (WISE) sources seen within the equatorial GAMA G12 field, located in the North Galactic Cap. Our motivation is to study and characterize the behavior of WISE source populations in anticipation of the deep multiwavelength surveys that will define the next decade, with the principal science goal of mapping the 3D large-scale structures and determining the global physical attributes of the host galaxies. In combination with cosmological redshifts, we identify galaxies from their WISE W1 (3.4 μm) resolved emission, and we also perform a star-galaxy separation using apparent magnitude, colors, and statistical modeling of star counts. The resulting galaxy catalog has ≃590,000 sources in 60 deg{sup 2}, reaching a W1 5σ depth of 31 μJy. At the faint end, where redshifts are not available, we employ a luminosity function analysis to show that approximately 27% of all WISE extragalactic sources to a limit of 17.5 mag (31 μJy) are at high redshift, z>1. The spatial distribution is investigated using two-point correlation functions and a 3D source density characterization at 5 Mpc and 20 Mpc scales. For angular distributions, we find that brighter and more massive sources are strongly clustered relative to fainter sources with lower mass; likewise, based on WISE colors, spheroidal galaxies have the strongest clustering, while late-type disk galaxies have the lowest clustering amplitudes. In three dimensions, we find a number of distinct groupings, often bridged by filaments and superstructures. Using special visualization tools, we map these structures, exploring how clustering may play a role with stellar mass and galaxy type.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Robust Cross-correlation-based Measurement of Clump Sizes in Galaxies

Stars form in molecular complexes that are visible as giant clouds (∼10{sup 5−6} M{sub ⊙}) in nearby galaxies and as giant clumps (∼10{sup 8−9} M{sub ⊙}) in galaxies at redshifts z≈1–3. Theoretical inferences on the origin and evolution of these complexes often require robust measurements of their characteristic size, which is hard to measure at limited resolution and often ill-defined due to overlap and quasi-fractal substructure. We show that maximum and luminosity-weighted sizes of clumps seen in star formation maps (e.g., Hα) can be recovered statistically using the two-point correlation function (2PCF) if an approximate stellar surface density map is taken as the normalizing random field. After clarifying the link between Gaussian clumps and the 2PCF analytically, we design a method for measuring the diameters of Gaussian clumps with realistic quasi-fractal substructure. This method is tested using mock images of clumpy disk galaxies at different spatial resolutions and perturbed by Gaussian white noise. We find that the 2PCF can recover the input clump scale at ∼20% accuracy, as long as this scale is larger than the spatial resolution. We apply this method to the local spiral galaxy NGC 5194, as well as to three clumpy turbulent galaxies from the DYNAMO-HST sample. In both cases, our statistical measurements of Hα clump size agree with previous measurements and with the estimated Jeans lengths. However, the new measurements are free from subjective choices when fitting individual clumps.

79 ASTRONOMY AND ASTROPHYSICS↗

A Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT): Part I. The clear-sky model

The solar energy industry often uses individual steps to empirically compute plane-of-array (POA) irradiance from horizontal irradiance and decompose it to narrow-wavelength bands. Conventional radiative transfer models designed for meteorological applications requires significant computing efforts in practice; however, they provide a physics-based solution of radiance and therefore are capable of computing spectral POA irradiances in a single step. In this study, we integrate the advantages of the current models and develop an innovative radiative transfer model, the Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT), to efficiently compute irradiances on inclined photovoltaics (PV) panels for 2002 narrow-wavelength bands from 0.28 to 4.0 um. This study is reported in two parts. Part I presents the methodology and performance evaluation of the new model under clear-sky conditions. The Simple Model of the Atmospheric Radiative Transfer of Sunshine (SMARTS), which was designed to compute clear-sky irradiances, is employed to rapidly provide the optical properties of a given clear-sky atmosphere. The clear-sky radiances in the narrow-wavelength bands are computed by considering three paths of photon transmission and solving the radiative transfer equation with the single-scattering approximation. The Bi-directional Transmittance Distribution Function (BTDF) of aerosols is given by their single-scattering phase function with a correction using a two-stream approximation. The validation analysis confirms that FARMS-NIT has improved accuracy compared to TMYSPEC as evaluated by both surface observations and a state-of-the-art radiative transfer model. This model substantially improves computational efficiency compared to other radiative transfer models though it uses slightly more computing time than TMYSPEC. Part II of this study addresses the model in cloud-sky conditions and will be published as a companion paper.

14 SOLAR ENERGY↗