Search NASA⌕ Search

SEARCH · Search NASA

Results for “Real Functions”

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 163 records · Page 9

Shortcomings of meta-GGA functionals when describing magnetism

Several recent studies have shown that SCAN, a functional belonging to the meta-generalized gradient approximation (MGGA) family, leads to significantly overestimated magnetic moments in itinerant ferromagnetic metals. However, this behavior is not inherent to the MGGA level of approximation since TPSS, for instance, does not lead to such severe overestimations. In order to provide a broader view of the accuracy of MGGAfunctionals for magnetism, we extend the assessment to more functionals but also to antiferromagnetic solids. The results show that to describe magnetism there is overall no real advantage in using a MGGA functional compared to GGAs. For both types of approximation, an improvement in ferromagnetic metals is necessarily accompanied by a deterioration (underestimation) in antiferromagnetic insulators, and vice versa. Furthermore, we also provide some analysis in order to understand in more detail the relation between the mathematical form of the functionals and the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Interpolation of computed gamma-ray detector response functions

Gamma-ray spectra measured by traditional detectors contain features that result from a combination of the effects of detector materials/geometry, the incident gamma-ray energy, and the angle of entry. The features, such as the full-energy photopeak, Compton continuum, annihilation peak, and escape peaks, are governed by simple relationships depending on incident energy and have been known for a long time. Monte Carlo computer simulations of gamma rays interacting with a detector will show these features, and with a resolution function applied, the results should look similar to real measurements. The traditional approach to creating a detector response function requires many separate simulations of monoenergetic gamma rays striking the detector. This paper presents a new approach to developing computed detector response functions. The new approach involves a much smaller number of monoenergetic gamma-ray simulations and uses interpolation to quickly generate the responses of gamma rays that were not simulated. During the interpolation process, the underlying physics equations are used to accurately compute the response of a given energy gamma ray from the small set of simulations. Such work enables accelerated generation of synthetic radiation detector data.

Detector response↗

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-space atomic dynamics in liquid gallium studied by inelastic neutron scattering

Gallium is a prototypical liquid metal and has gained renewed attention because of its unique properties. Characterizing and elucidating its atomic dynamics remains elusive despite numerous studies, primarily because of the challenges in quantifying atomic-scale dynamics in liquids. Recent developments in inelastic neutron scattering enable us to measure the Van Hove correlation function that describes the real-space motion of liquid atoms. Here, in this work, we use this approach to reveal the dynamics in gallium liquids and find the co-existence of two dynamical medium-range orders (MROs), which have a dynamical behavior distinct from that of the short-range order (SRO). We propose that these MROs are driven by global forces in the form of two density waves, as a direct consequence of the underlying competition between ionic core repulsion and valence electron cohesion. We suggest that the density wave approach is not only applicable to other metallic liquids exhibiting similar structural anomalies, but also offers a promising direction for elucidating the dynamics of complex liquids and glasses by linking electronic-state fluctuations to atomic dynamics.

Hua, Chengyun [Oak Ridge National Laboratory (ORNL↗

NeuralCubes: Deep Representations for Visual Data Exploration

Visual exploration of large multi-dimensional datasets has seen tremendous progress in recent years, allowing users to express rich data queries that produce informative visual summaries, all in real time. Techniques based on data cubes are some of the most promising approaches. However, these techniques usually require a large memory footprint for large datasets. To tackle this problem, we present NeuralCubes: neural networks that predict results for aggregate queries, similar to data cubes. NeuralCubes learns a function that takes as input a given query, for instance, a geographic region and temporal interval, and outputs the result of the query. The learned function serves as a real-time, low-memory approximator for aggregation queries. Our models are small enough to be sent to the client side (e.g. the web browser for a web-based application) for evaluation, enabling data exploration of large datasets without database/network connection. Here, we demonstrate the effectiveness of NeuralCubes through extensive experiments on a variety of datasets and discuss how NeuralCubes opens up opportunities for new types of visualization and interaction.

97 MATHEMATICS AND COMPUTING↗

Predicting fatigue from heart rate signatures using functional logistic regression

Physical fatigue can have adverse effects on humans in extreme environments. Therefore, being able to predict fatigue using easy to measure metrics such as heart rate (HR) signatures has potential to have an impact in real-life scenarios. We apply a functional logistic regression model that uses HR signatures to predict physical fatigue, where physical fatigue is defined in a data-driven manner. Data were collected using commercially available wearable devices on 47 participants hiking the 20.7-mile Grand Canyon rim-to-rim trail in a single day. Fitted model provides good predictions and interpretable parameters for real-life application.

60 APPLIED LIFE SCIENCES↗

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low-rank approximation of the dielectric matrix

We present an efficient preconditioning technique for accelerating the fixed-point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low-rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed-point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous self-consistent field iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from ~100 to 1100 atoms (~500–20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized d orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20–30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed-point iteration in some cases. Lastly, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4-fold reduction in computational cost for the total ground-state calculation compared to other preconditioners.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum effects on the dynamics and properties of soft materials

The quantum effects of nuclear and electronic motion play an important role in the structure, dynamics, and function of soft materials, yet they are difficult to capture with conventional classical simulations or static electronic–structure methods. In this work several complementary approaches for treating quantum effects in polymeric and soft–matter systems are demonstrated, with a focus being on the hydrogen-bonded networks, ion and charge transport, and photoactive chromophores. The proton transfer, tunneling, and isotope effects are captured within the reduced-dimensionality models by implementing grid-based nuclear quantum dynamics in terms of the discrete variable and Fourier bases. The nuclear quantum dynamics is extended to larger systems by employing the quantum trajectories and quantum–thermal bath schemes combined with on-the-fly electronic structure, enabling the description of high-dimensional polymeric environments at feasible cost. The dynamics in the electronic degrees of freedom, simulating the optical response in large chromophores such as chlorophylls, is performed using the real-time time-dependent density functional theory implemented in the real-space multigrid (RMG) code. These approaches are demonstrated on case studies of the proton and hydroxide transport in hydrated polymer membranes, charge transfer in conjugated polymers, and the optical spectra of chlorophyll chromophores relevant to polymerized chlorophyll materials and chlorophyll–polymer hybrids. The reviewed methods and applications highlight practical routes of including quantum effects in simulations of soft functional materials.

Garashchuk, Sophya [Univ. of South Carolina, Colum↗

Oxidative Deboronation of Boronic Acids by Hydrogen Peroxide in Planta Generates Borate for Cross-Linking of Rhamnogalacturonan II

Vascular plants require boron to cross-link the rhamnogalacturonan-II (RG-II) domain of pectin to form functional cell walls. Boronic acids, which form reversible esters with cis-diols like borate, have been proposed to influence RG-II cross-linking, though the mechanism remains unclear. We used suspension-cultured rose cells adapted to grow without boron to investigate the effect of boronic acids on RG-II dimerization. When grown with phenylboronic acid (PBA) as the sole boron source, nearly all RG-II was crosslinked, whereas methylboronic acid (MBA) only partially restored cross-linking. In contrast, in vitro assays showed that homogeneous RG-II monomers did not dimerize with alkyl or aryl boronic acids unless supplemented with hydrogen peroxide (H2O2), which oxidatively converts boronic acids to boric acid. Real-time NMR spectroscopy and density functional theory calculations provided insight into the reaction mechanism and energetics of oxidation respectively. Together, our data show that exogenous boronic acids are a source of boric acid for plants, and that the deboronation reaction generates aryl or alkyl alcohol byproducts that can undergo further chemical modification in planta. The fate and potential roles of these byproducts in planta remain to be determined.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Soft pions and transport near the chiral critical point

During the expansion of a heavy ion collision, the system passes close to the $\textit{O}$(4) critical point of QCD, and thus the fluctuations of the order parameter $(σ, \overrightarrow{π}$) are expected to be enhanced. In this work, our goal is to compute how these enhanced fluctuations modify the transport coefficients of QCD near the pseudocritical point. We also make a phenomenological estimate for how chiral fluctuations could effect the momentum spectrum of soft pions. We first formulate the appropriate stochastic hydrodynamic equations close to the $\textit{O}$(4) critical point. Then, working in mean field, we determine the correlation functions of the stress tensor and the currents which result from this stochastic real-time theory, and use these correlation functions to determine the scaling behavior of the transport coefficients. The hydrodynamic theory also describes the propagation of pion waves, fixing the scaling behavior of the dispersion curve of soft pions. We present scaling functions for the shear viscosity and the charge conductivities near the pseudocritical point and estimate the absolute magnitude of the critical fluctuations to these parameters and the bulk viscosity. Using the calculated pion dispersion curve, we estimate the expected critical enhancement of soft pion yields, and this estimate provides a plausible explanation for the excess seen in experiment relative to ordinary hydrodynamic computations. Our results motivate further phenomenological and numerical work on the implications of chiral symmetry on real-time properties of thermal QCD near the pseudocritical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

ForSE: A GAN-based Algorithm for Extending CMB Foreground Models to Subdegree Angular Scales

We present ForSE (Foreground Scale Extender), a novel Python package that aims to overcome the current limitations in the simulation of diffuse Galactic radiation, in the context of cosmic microwave background (CMB) experiments. ForSE exploits the ability of generative adversarial neural networks (GANs) to learn and reproduce complex features present in a set of images, with the goal of simulating realistic and non-Gaussian foreground radiation at subdegree angular scales. This is of great importance in order to estimate the foreground contamination to lensing reconstruction, delensing, and primordial B-modes for future CMB experiments. We applied this algorithm to Galactic thermal dust emission in both total intensity and polarization. Our results show how ForSE is able to generate small-scale features (at 12') having as input the large-scale ones (80'). The injected structures have statistical properties, evaluated by means of the Minkowski functionals, in good agreement with those of the real sky and which show the correct amplitude scaling as a function of the angular dimension. Furthermore, the obtained thermal dust Stokes Q and U full-sky maps as well as the ForSE package are publicly available for download.

79 ASTRONOMY AND ASTROPHYSICS↗

Origin of charge density wave in topological semimetals SrAl 4 and EuAl 4

Topological semimetals in BaAl 4 -type structure show many interesting behaviors, such as charge density wave (CDW) in SrAl 4 and EuAl 4 , but not the isostructural and isovalent BaAl 4 , SrGa 4 , and BaGa 4 . Here using Wannier functions based on density functional theory, we calculate the susceptibility functions with millions of k-points to reach the small q-vector and study the origin and driving force behind the CDW. Our comparative study reveals that the origin of the CDW in SrAl 4 and EuAl 4 is the strong electron-phonon coupling interaction for the transverse acoustic mode at small q-vector along the Γ-Z direction besides the maximum of the real part of the susceptibility function from the nested Fermi surfaces of the Dirac-like bands, which explains well the absence of CDW in the other closely related compounds in a good agreement with experiment. We also connect the different CDW behaviors in the Al compounds to the macroscopic elastic properties.

36 MATERIALS SCIENCE↗

Probing the Chemical Kinetics of Minimalist Functional Group Gasoline Surrogates

We report surrogate mixtures are routinely used for understanding gasoline fuel combustion in engine simulations. The general trend in surrogate formulation has been to increase the number of fuel components in a mixture to better emulate real fuel properties. Recently, a new surrogate design strategy based on functional group analysis of real gasolines was proposed using a minimal number of species [minimalist functional group (MFG)-approach]. MFG surrogates (having just one or two components) could experimentally capture the ignition delay time (IDT), threshold sooting index, and smoke point of different gasoline fuels with hundreds of components. However, other combustion characteristics were not explored, and kinetic modeling of MFG surrogates was not reported. These aspects are addressed in this paper, where the combustion behavior of MFG surrogates for various gasolines was assessed by simulating IDT, jet-stirred reactor oxidation, and premixed laminar flame speeds using chemical kinetic modeling. MFG simulations were compared with experimental data of the real gasolines as well as with the more complex multicomponent (five to nine species) surrogates. This study reveals that binary MFG surrogate mixtures are capable of accurately simulating the combustion behavior of more complex gasoline fuels with hundreds of components.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MPC-based double-layer real-time conditional cSelf-restoration for interconnected microgrids

In this report, we propose a novel model predictive control(MPC)-based real-time conditional self-restoration energy management system (CSR-EMS) for interconnected microgrids (IMGs) integrated with renewable energy sources (RESs) and energy storage systems (ESSs). Superior to the existing IMG self-restoration methods, the “conditionality” of the proposed CSR-EMS can economically realize self-restoration and grid-assisted restoration during energy deficiency or faults, in both islanded and grid-connected modes. Cost minimization is implemented as the objective function to judge in real-time which restoration mode is economically preferred. The proposed CSR-EMS comprises two layers–the lower layer operates locally to eliminate electricity fluctuations created by RESs and ensure economic effectiveness within an MG, whereas the upper layer oversees the real-time operational status of the IMG system and determines power exchange among microgrids (MGs) during abnormalities. In detail, when a microgrid inside the IMG system experiences an energy deficiency, the CSR-EMS, on an MPC basis, intelligently optimizes power production from each dispatchable distributed generator (DG), ESS, power imported from the main grid, and power exchange among the IMGs to maintain the demand–supply balance, while considering system recovery cost, state of charge (SoC) of ESSs and operation modes of the IMGs (i.e., grid-connected or islanded mode). Simulation results and comparisons with existing IMG self-healing EMSs demonstrate the economic efficacy of the proposed CSR-EMS strategy during normal and abnormal operations, which can be used as an energy control framework for modern power systems with multiple interconnected microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Scalable implementation of polynomial filtering for density functional theory calculation in PARSEC

In this work, we present an efficient implementation of polynomial filtering methods in PARSEC, a real-space pseudopotential based Kohn–Sham density functional theory solver. The implementation described here improves upon a Chebyshev-filtered subspace iteration algorithm used in the previous version of PARSEC. We present a hybrid polynomial filtering scheme that combines Chebyshev-filtered subspace iteration and a spectrum slicing method that partitions the spectrum into several spectral slices and uses bandpass-filtered subspace iteration to compute approximate eigenpairs within each interior slice simultaneously. We describe a procedure to partition a spectrum and construct polynomial filters. We also discuss a number of practical issues such as the use of appropriate data layouts for carrying out the computation on a two-dimensional process grid and how to achieve good load balance by allocating an appropriate number of process groups to each spectral slice. Numerical examples are presented to demonstrate the effectiveness of the hybrid polynomial filtering method as well as the superior parallel scalability of spectrum slicing in comparison to that of Chebyshev-filtered subspace iteration.

97 MATHEMATICS AND COMPUTING↗

DERMS-RT (Distributed Energy Resource Management Solution using Real-Time Optimization) [SWR-20-46]

DERMS-RT provides a distributed energy resource management solution using real-time optimization to control different types of distributed energy resources (DERs) — such as rooftop PVs, battery energy storage systems, HVAC loads, electric water heater loads, and EV charging loads, and use these DERs to provide distribution voltage regulation and peak demand management services. The control approach implemented by DERM-RT is developed based on the real-time optimal power flow and distributed control previously developed by NREL researchers. The software codes are designed to provide modular DER management solution that can be easily applied to control heterogeneous DERs in distribution grids in coordination with centralized utility management systems. Also, DERMS-RT provides application programming interface (API) to extract useful grid model information from an open-source power system simulation software, use the information to solve the optimization problem, and apply the optimal set-points to the controlled DERs. In this way, DERMS-RT provides a plug-and-play function for implementing the real-time DER control on different utility system models, and it can be used as a simulation testbed to demonstrate the impact of real-time DER optimization on improving grid operations.

Ding, Fei↗

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗