Search NASA⌕ Search

SEARCH · Search NASA

Results for “integrated gradients”

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 73 records · Page 4

Differentiating densities on smooth manifolds

Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. Furthermore, we highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands.

97 MATHEMATICS AND COMPUTING↗

Mo-100 Disk Structural Integrity Tests I

The NorthStar pressed and sintered 29 mm diameter, 0.5 mm thick Mo-100 disks will be subjected to large temperature gradients when in beam during production. The beam preferentially heats the disk center, while the aggressive cooling by the helium keeps areas outside the beam spot relatively cool. The result is a temperature gradient of up to 900°C over the radial distance of 14.5 mm. To test the disk integrity with this high temperature gradient, a special holder for a single disk was made such that proper cooling can be applied as in the design target and the disk heating can be replicated with induction heating. NorthStar press and sintered disks were used for these tests.

36 MATERIALS SCIENCE↗

Effective kinetics driven by dynamic concentration gradients under coupled transport and reaction

Biogeochemical reaction kinetics are generally established from batch reactors where concentrations are uniform. In natural systems, many biogeochemical processes are characterized by spatially and temporally variable concentration gradients that often occur at scales which are not resolved by field measurements or biogeochemical and reactive transport models. Yet, it is not clear how these sub-scale chemical gradients affect reaction kinetics compared to batch kinetics. Here we investigate this question by studying the paradigmatic case of localized pulses of solute reacting with a solid or a dissolved species in excess. Additionally, we consider non-linear biogeochemical reactions, representative of mineral dissolution, adsorption and redox reactions, which we quantify using simplified power-law kinetics. The combined effect of diffusion and reaction leads to effective kinetics that differ quantitatively and qualitatively from the batch kinetics. Depending on the nonlinearity (reaction order) of the local kinetics, these effects lead to either enhancement or decrease of the overall reaction rate, and result in a rich variety of reaction dynamics. We derive analytical results for the effective kinetics, which are validated by comparison to direct numerical simulations for a broad range of Damköhler numbers and reaction order. Our findings provide new insights into the interpretation of imperfectly mixed lab experiments, the effective kinetics of field systems characterized by intermittent reactant release and the integration of sub-scale concentration gradients in reactive transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mountain biodiversity and ecosystem functions: interplay between geology and contemporary environments

Although biodiversity and ecosystem functions are strongly shaped by contemporary environments, such as climate and local biotic and abiotic attributes, relatively little is known about how they depend on long-term geological processes. Here, along a 3000-m elevational gradient with tectonic faults on the Tibetan Plateau (that is, Galongla Mountain in Medog County, China), we study the joint effects of geological and contemporary environments on biological communities, such as the diversity and community composition of plants and soil bacteria, and ecosystem functions. We find that these biological communities and ecosystem functions generally show consistent elevational breakpoints at 2000–2800 m, which coincide with Indus-Yalu suture zone fault and are similar to the elevational breakpoints of soil bacteria on another mountain range 1000 km away. Mean annual temperature, soil pH and moisture are the primary contemporary determinants of biodiversity and ecosystem functions, which support previous findings. However, compared with the models excluding geological processes, inclusion of geological effects, such as parent rock and weathering, increases 67.9 and 35.9% of the explained variations in plant and bacterial communities, respectively. Such inclusion increases 27.6% of the explained variations in ecosystem functions. The geological processes thus provide additional links to ecosystem properties, which are prominent but show divergent effects on biodiversity and ecosystem functions: parent rock and weathering exert considerable direct effects on biodiversity, whereas indirectly influence ecosystem functions via interactions with biodiversity and contemporary environments. Thus, the integration of geological processes with environmental gradients could enhance our understanding of biodiversity and, ultimately, ecosystem functioning across different climatic zones.

54 ENVIRONMENTAL SCIENCES↗

Non-classical ion-indiffusion processes in quasi-one-dimensional ionic conductor KTiOPO4

Developments in coercive field engineering and waveguide writing via ion-exchange in the ferroelectric ionic-conductor KTiOPO4 (KTP) have led to breakthroughs within the field of nonlinear optics and quantum photonics. However, the microscopic diffusion dynamics of ion-exchange in KTP are poorly understood and show signs of non-classical diffusion dynamics. The precise control of ion-exchange gradients is critical for developing KTP into a mature, integrated nonlinear optical material but requires a level of control that, in turn, demands a deeper understanding of the ion-exchange process. We employ direct ToF-SIMS mapping to untangle the complex diffusion dynamics of the commonly used ion-exchange dopant ions: Rb, Ba, and K. We map periodic diffusion gradients of two exchanges, one dominated by Rb-ions and the other by Ba-ions. For the first time, we show the interdependent nature of the Ba- and Rb ion in the ion-exchange diffusion gradient. The low incorporation of the Ba ion is mapped and quantified, and dopant–crystal interactions are identified. This work highlights the complex diffusion dynamics of dopants in the KTP structure and lays the foundation for deterministic control of diffusion gradients for the next generation of integrated nonlinear and quantum photonic devices.

Barrett, Laura [ORNL]↗

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗

High-Nickel Cathodes with Mechanical and Interfacial Robustness via Tailored Concentration Gradients for Stable Li-Ion Batteries

Here, we have developed a versatile mathematical framework integrated with an automated reactor system to design and reify highly customizable full concentration gradient (FCG) in high-nickel cathodes for advanced Li-ion batteries. This method provides precise and independent control of the average composition, slope, and curvature of FCGs, enabling the optimization of structural and mechanical properties of the cathode materials. We have showcased this method with Ni 0.8 Co 0.1 Mn 0.1 (OH) 2 precursors of controlled FCGs, which unlocked an optimized cathode with excellent cycling stability without crack formation after repeated cycles. This work opens up new possibilities for the design and manufacturing of advanced cathode materials, enabling safer, high-performance batteries.

25 ENERGY STORAGE↗

Detection of Diversion in a Realistic Heat Pipe Microreactor Using Supervised Machine Learning

Microreactors (MRs) pose new challenges for international safeguards. Here, their small size and mass reproducibility make them ideal for deployment in greater numbers and in remote locations, making the job of safeguards inspectors more challenging. Machine learning (ML) is currently being applied to many fields to augment human performance and increase automation; in particular, ML could be used to provide insight for international inspectors to help detect the diversion of nuclear fuel from MR cores. Four ML model types (k-nearest neighbors, decision tree, random forest, and histogram-based gradient boosted ensemble) were trained on integrated flux and critical control drum angle data generated with Serpent 2 for a realistic heat pipe MR design, achieving nearly 100% binary classification accuracy of nominal and diversion core configurations by the end of 1 full power year for three of the four model types. Regression model variants were also trained, using the same input data, for predicting the number of fuel pins diverted. Root-mean-square errors below 5% of the total number of fuel pins were achieved by the 1 full power year mark for all models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Quandary

Quandary numerically simulates and optimizes the time-evolution of open quantum systems. The underlying dynamics are modelled by Lindblad's master equation, a linear ordinary differential equation (ODE) describing quantum systems interacting with the environment. Quandary solves this ODE numerically by applying a time-stepping integration scheme, and utilizes a gradient-based optimization approach to determine optimal control pulses that drive the quantum system to a desired target state. Two optimization objectives are considered: (a) Unitary gate optimization that finds controls to realize a unitary gate transformation, and (b) optimal reset that aims to drive the quantum system to the ground states. Gradient-based optimization schemes utilizing Petsc's Tao optimization package are applied to generate control pulses that minimize the respective measure. To evaluate the gradient of the objective function, the discrete adjoint method is used while leveraging techniques from Algorithmic Differentiation to produce exact and consistent gradients. To mitigate excessive execution run times, the software can be build together with the XBraid software library which provides a parallelization strategy to distribute the time-evolution of the underlying dynamics onto multiple processor.

Petersson, NilsA.↗

TEAM project review, year 1

This note summarizes our research activities within the TEAM project between October 2019 and November 2020, funded by the ASCR Advanced Research in Quantum Computing program. Three research efforts have begun during the review period. The first effort, resulting in the JuQBox software, implemented a novel algorithm for optimal control in closed quantum systems based on symplectic time integration and exact calculations of the gradient of the discretized objective function. In our second effort, resulting in a code called Quandary, we consider the optimal control problem for open quantum systems. Here the evolution of the quantum system is modeled by the Lindblad master equation. In this code, we are considering both the realization of quantum gates as well as the optimal reset problem in which a general mixed quantum state is driven towards its ground state. In our third research effort we have begun the development of computational software for what we call Lindblad learning, i.e. the process of characterizing the equations governing our quantum computing hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessment of Electrode Contamination Mitigation at 0.5 MA Scale

The Z Machine at Sandia National Laboratories uses current pulses with peaks up to 27 MA to drive target implosions and generate high energy density conditions of interest for stockpile stewardship programs pertinent to the NNSA program portfolio . Physical processes in the region near the Z Machine target create electrode plasmas which seed parasitic current loss that reduce the performance and output of a Z experiment. Electrode surface contaminants (hydrogen, water, hydrocarbons) are thought to be the primary constituent of electrode plasmas which contribute to loss mechanisms. The Sandia team explore d in situ heating and plasma discharge techniques by integrating requisite infrastructure into Sandia's Mykonos LTD accelerator, addressing potential impacts to accelerator operation, and reporting on the impact of these techniques on electrode plasma formation and shot performance. The in situ discharge cleaning utilizes the electrodes of the accelerator to excite an argon-oxygen plasma to sputter and chemically react contaminants from electrode surfaces. Insulating breaks are required to isolate the plasma in electrode regions where loss processes are most likely to occur. The shots on Mykonos validate that these breaks do not perturb experiment performance, reducing the uncertainty on the largest unknown about the in situ cleaning system. Preliminary observations with electrical and optical diagnostics suggest that electrode plasma formation is delayed, and overall inventory has been substantively reduced. In situ heating embeds cartridge heaters into accelerator electrodes and employs a thermal bakeout to rapidly desorb contaminants from electrode surfaces. For the first time, additively manufactured (AM) electrode assemblies were used on a low impedance accelerator to integrate cooling channels and manage thermal gradients. Challenges with poor supplier fabrication to specifications, load alignment, thermal expansion and hardware movement and warpage appears to have introduced large variability in observed loss, though, preventing strong assertions of loss reduction via in situ heating. At this time, an in situ discharge cleaning process offers the lowest risk path to reduce electrode contaminant inventories on Z, though we recommend continuing to develop both approaches. Additional engineering and testing are required to improve the implementation of both systems. .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization

Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite-difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 12.78 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.

17 WIND ENERGY↗

Time-varying gradient metasurface with applications in all-optical beam steering

Integrating the large, subpicosecond nonlinear optical response of epsilon-near-zero (ENZ) materials with the broad design freedoms of plasmonic metasurfaces shows potential for creating rapidly modulated optical devices with possible applications in telecommunications, sensing, and reactive beam steering. In this work, we experimentally investigate a metasurface consisting of a plasmonic gradient array on a thin layer of indium tin oxide (ITO), characterize how incident probe pulses diffract from a system as it is being dynamically modulated by a pump pulse at wavelengths near the ENZ region. Angular shifts in the diffraction orders are observed and can be principally attributed to the adiabatic wavelength conversion of the probe as it witnesses the temporal change of index induced by the pump. Of note, the asymmetric gradient metasurface, considered to be a blazed diffraction grating, shows significantly different dynamic responses for different diffraction orders. The free-space wavelength shift to +1 and -1 diffraction orders is 6 and 12 nm, resulting in steering angle changes of 0.65 and 1.5°, respectively.

42 ENGINEERING↗

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗

ZEUS: An Efficient GPU Optimization Method Integrating PSO, BFGS, and Automatic Differentiation

We introduce a novel, efficient computational method, ZEUS, for numerical optimization, and provide an open-source implementation. It has four key ingredients: (1) particle swarm optimization (PSO), (2) the use of the Broyden-Fletcher-Goldfarb-Shanno (BFGS) method, (3) automatic differentiation (AD), and (4) GPUs. Our approach addresses the computational challenges inherent in high-dimensional, non-convex optimization problems. In the first phase of the algorithm, we get a potentially good set of starting points using PSO. Thereafter, we run BFGS independently in parallel from these starting points. BFGS is one of the best-performing algorithms for numerical optimization. However, it requires the gradient of the function being optimized. ZEUS integrates automatic differentiation into BFGS thus avoiding the need for the user to calculate derivatives explicitly. The use of GPUs allows ZEUS to speed up the calculations substantially. We carry out systematic studies to explore the trade-offs between the number of PSO iterations taken, starting points, and BFGS iteration depth. We show that a handful of iterations of PSO can improve global convergence when combined with BFGS. We also present performance studies using common test functions. The source code can be found at https://github.com/fnal-numerics/global-optimizer-gpu.

Soos, Dominik [Old Dominion U.]↗

Enabling Long-range Exploration in Minimization of Multimodal Functions

We consider the problem of minimizing multi-modal loss functions with a large number of local optima. Since the local gradient points to the direction of the steepest slope in an infinitesimal neighborhood, an optimizer guided by the local gradient is often trapped in a local minimum. To address this issue, we develop a novel nonlocal gradient to skip small local minima by capturing major structures of the loss’s landscape in black-box optimization. The nonlocal gradient is defined by a directional Gaussian smoothing (DGS) approach. The key idea of DGS is to conducts 1D long-range exploration with a large smoothing radius along d orthogonal directions in Rd, each of which defines a nonlocal directional derivative as a 1D integral. Such long-range exploration enables the nonlocal gradient to skip small local minima. The d directional derivatives are then assembled to form the nonlocal gradient. We use the Gauss-Hermite quadrature rule to approximate the d 1D integrals to obtain an accurate estimator. The superior performance of our method is demonstrated in three sets of examples, including benchmark functions for global optimization, and two real-world scientific problems.

Zhang, Jiaxin↗

Evaluating the Ratio of Electron and Hole Mobilities from a Single Bulk Sample Using Photo-Seebeck Effect

When a semiconductor is under photoexcitation, the voltage response to a temperature gradient is the photo-Seebeck effect. In this work, we study this effect, focusing on the contribution from transport of photo-excited carriers. We demonstrate that by combining photo-Seebeck with photoconductivity measurements, one can determine the ratio between electron and hole mobilities, and hence both of them when one is known. This is found for the case of defect-free samples, where no detail on the absorbance, carrier lifetime or recombination is necessary. Our method reported here does not require chemical doping, which could introduce defects and is often not feasible. It applies to both thin film and bulk samples. Experiment wise, photo-Seebeck effect is relatively easy to implement, or added to existing systems. In a broader context, for semiconductors with significant influence from defects, our result suggests that the photo-Seebeck behavior can still be understood. In this case another photo-transport property is necessary, in order to identify the mobilities of carriers and information regarding the defects. This framework integrates the information from photoexcitation and thermal gradients to provide a general method to determine fundamental electronic properties of materials.

36 MATERIALS SCIENCE↗

Colour gradients of low-redshift galaxies in the DESI Legacy Imaging Survey

ABSTRACT Radial colour gradients within galaxies arise from gradients of stellar age, metallicity, and dust reddening. Large samples of colour gradients from wide-area imaging surveys can complement smaller integral-field spectroscopy data sets and can be used to constrain galaxy formation models. Here, we measure colour gradients for low-redshift galaxies (z < 0.1) using photometry from the DESI Legacy Imaging Survey DR9. Our sample comprises ∼93 000 galaxies with spectroscopic redshifts and ∼574 000 galaxies with photometric redshifts. We focus on gradients across a radial range 0.5Reff to Reff, which corresponds to the inner disc of typical late-type systems at low redshift. This region has been the focus of previous statistical studies of colour gradients and has recently been explored by spectroscopic surveys such as MaNGA. We find that the colour gradients of most galaxies in our sample are negative (redder towards the centre), consistent with the literature. We investigate empirical relationships between colour gradient, average g − r and r − z colour, Mr, M⋆, and sSFR. Trends of gradient strength with Mr (M⋆) show an inflection around Mr ∼ −21 ($\log _{10} \, M_\star /\mathrm{M_\odot }\sim 10.5$). Below this mass, colour gradients become steeper with increasing M⋆, whereas colour gradients in more massive galaxies become shallower. We find that positive gradients (bluer stars at smaller radii) are typical for galaxies of $M_{\star }\sim 10^{8}\, \mathrm{M_\odot }$. We compare our results to age and metallicity gradients in two data sets derived from fits of different stellar population libraries to MaNGA spectra, but find no clear consensus explanation for the trends we observe. Both MaNGA data sets seem to imply a significant contribution from dust reddening, in particular, to explain the flatness of colour gradients along the red sequence.

Astronomy & Astrophysics↗