Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “exponential integrators”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Physics‐Informed Machine Learning for Inverse Design of Optical Metamaterials

Optical metamaterials manipulate light through various confinement and scattering processes, offering unique advantages like high performance, small form factor and easy integration with semiconductor devices. However, designing metasurfaces with suitable optical responses for complex metamaterial systems remains challenging due to the exponentially growing computation cost and the ill‐posed nature of inverse problems. To expedite the computation for the inverse design of metasurfaces, a physics‐informed deep learning (DL) framework is used. A tandem DL architecture with physics‐based learning is used to select designs that are scientifically consistent, have low error in design prediction, and accurate reconstruction of optical responses. The authors focus on the inverse design of a representative plasmonic device and consider the prediction of design for the optical response of a single wavelength incident or a spectrum of wavelength in the visible light range. The physics‐based constraint is derived from solving the electromagnetic wave equations for a simplified homogenized model. The model converges with an accuracy up to 97% for inverse design prediction with the optical response for the visible light spectrum as input, and up to 96% for optical response of single wavelength of light as input, with optical response reconstruction accuracy of 99%.

Sarkar, Sulagna↗
DOE OSTI2020

Quantum complexity of time evolution with chaotic Hamiltonians

We study the quantum complexity of time evolution in large-N chaotic systems, with the SYK model as our main example. This complexity is expected to increase linearly for exponential time prior to saturating at its maximum value, and is related to the length of minimal geodesics on the manifold of unitary operators that act on Hilbert space. Using the Euler-Arnold formalism, we demonstrate that there is always a geodesic between the identity and the time evolution operator e -iHt whose length grows linearly with time. This geodesic is minimal until there is an obstruction to its minimality, after which it can fail to be a minimum either locally or globally. We identify a criterion — the Eigenstate Complexity Hypothesis (ECH) — which bounds the overlap between offdiagonal energy eigenstate projectors and the k-local operators of the theory, and use it to argue that the linear geodesic will at least be a local minimum for exponential time. We show numerically that the large-N SYK model (which is chaotic) satisfies ECH and thus has no local obstructions to linear growth of complexity for exponential time, as expected from holographic duality. In contrast, we also study the case with N = 2 fermions (which is integrable) and find short-time linear complexity growth followed by oscillations. Our analysis relates complexity to familiar properties of physical theories like their spectra and the structure of energy eigenstates and has implications for the hypothesized computational complexity class separations PSPACE * BQP/poly and PSPACE * BQSUBEXP/subexp, and the “fast-forwarding” of quantum Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Extracting scattering amplitudes for arbitrary two-particle systems with one-particle left-hand cuts via lattice QCD

We derive a general formalism that relates the spectrum of two-particle systems in a finite volume to physical scattering amplitudes, taking into account the presence of any left-hand branch cuts due to single-particle exchanges. The method first relates the finite-volume spectrum to an infinite-volume short-range quantity, denoted ${\mathcal{M}}_0$, and then relates the latter to the physical scattering amplitudes via known integral equations. The derivation of both relations is performed using all-orders perturbation theory and is exact up to neglected exponentially suppressed volume dependence. The relations hold for arbitrary two-particle systems with any number of coupled channels, non-identical and non-degenerate particles, and any intrinsic spin.

algorithms↗
DOE OSTI2023

Challenging the Curse of Dimensionality in Multidimensional Numerical Integration by Using a Low-Rank Tensor-Train Format

Numerical integration is a basic step in the implementation of more complex numerical algorithms suitable, for example, to solve ordinary and partial differential equations. The straightforward extension of a one-dimensional integration rule to a multidimensional grid by the tensor product of the spatial directions is deemed to be practically infeasible beyond a relatively small number of dimensions, e.g., three or four. In fact, the computational burden in terms of storage and floating point operations scales exponentially with the number of dimensions. This phenomenon is known as the curse of dimensionality and motivated the development of alternative methods such as the Monte Carlo method. The tensor product approach can be very effective for high-dimensional numerical integration if we can resort to an accurate low-rank tensor-train representation of the integrand function. In this work, we discuss this approach and present numerical evidence showing that it is very competitive with the Monte Carlo method in terms of accuracy and computational costs up to several hundredths of dimensions if the integrand function is regular enough and a sufficiently accurate low-rank approximation is available.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Exponential vs Gaussian Correlation Functions in the Characterization of Block Copolymer Grain Structure by Depolarized Light Scattering

Block copolymer (BCP) grain structure affects the mechanical, optical, and electrical properties of BCP materials, making the accurate characterization of this grain structure an important goal. In this study, improved BCP grain parameters were obtained by employing an exponentially decaying correlation function within the ellipsoidal grain model, instead of the Gaussian correlation function that was used in previous work. The exponential correlation function provides a better fit to the experimental depolarized light scattering data, which outweighs the disadvantage that it requires numerical integration to obtain the model scattered intensity.

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

The rate of microtubule breaking increases exponentially with curvature

Microtubules, cylindrical assemblies of tubulin proteins with a 25 nm diameter and micrometer lengths, are a central part of the cytoskeleton and also serve as building blocks for nanobiodevices. Microtubule breaking can result from the activity of severing enzymes and mechanical stress. Breaking can lead to a loss of structural integrity, or an increase in the numbers of microtubules. We observed breaking of taxol-stabilized microtubules in a gliding motility assay where microtubules are propelled by surface-adhered kinesin-1 motor proteins. We find that over 95% of all breaking events are associated with the strong bending following pinning events (where the leading tip of the microtubule becomes stuck). Furthermore, the breaking rate increased exponentially with increasing curvature. These observations are explained by a model accounting for the complex mechanochemistry of a microtubule. The presence of severing enzymes is not required to observe breaking at rates comparable to those measured previously in cells.

42 ENGINEERING↗
DOE OSTI2026

Integrating AI Data Centers with the Power Grid

The rapid expansion of artificial intelligence (AI) has triggered an unprecedented surge in electricity demand, with US data center energy use projected to double or triple 2023 levels by 2028. This exponential growth places strain on grid infrastructure, which can hinder timely construction of desired computing capacity. To bridge this supply-demand gap, utilities and AI developers are increasingly turning to demand flexibility, a strategy that incentivizes shifting or reducing power use during peak periods of grid stress. Data centers are uniquely equipped for flexible operations due to their digital workloads, built-in redundancy, and onsite energy assets. This article outlines four primary mechanisms to enable data center flexibility: computational load flexibility (shifting tasks temporally or geographically), flexible use of core facility infrastructure adjustments, energy storage utilization, and onsite electricity generation. To encourage adoption, utilities are deploying new tariff designs, including voluntary interruptible service riders, mandated flexibility requirements, and streamlined interconnection processes for flexible loads. For the highly capitalized and rapidly growing AI industry, the primary motivators for embracing these strategies are expediting facility interconnection, satisfying emerging regulatory mandates, and mitigating community resistance. While demand flexibility cannot substitute the long-term need for new bulk power generation, it serves as an essential, immediate solution for enabling near-term deployment. By transforming data centers from grid stressors into stabilizing assets, flexible operations can ensure reliable grid integration, ease market pressures, and support a resilient power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

High-Accuracy Mesh-Free Quadrature for Trimmed Parametric Surfaces and Volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theorem-based numerical quadrature scheme for integrating functions over trimmed parametric surfaces and volumes. The algorithm relies on two fundamental steps: (1) We iteratively reduce the dimensionality of integration using the generalized Stokes theorem to line integrals over trimming curves, and (2) we employ numerical antidifferentiation in the generalized Stokes theorem using high-order quadrature rules. Additionally, the scheme achieves exponential convergence up to trimming curve approximation error and has applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order curvilinear meshes, and initialization of multi-material simulations. We compare the quadrature scheme to commonly-used quadrature schemes in the literature and show that our scheme is much more efficient in terms of number of quadrature points used. We provide an open-source implementation of the scheme in MATLAB as part of QuaHOG, a software package for Quadrature of High-Order Geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spectral scheme for atomic structure calculations in density functional theory

In this study, we present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number Ζ = 1 to 83.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion sign problem is crucial, as the signal-to-noise ratio decreases exponentially with imaginary time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parametrized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. Here, we only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross sections of larger nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Particle-number distribution in large fluctuations at the tip of branching random walks

Here, we investigate properties of the particle distribution near the tip of one-dimensional branching random walks at large times t , focusing on unusual realizations in which the rightmost lead particle is very far ahead of its expected position, but still within a distance smaller than the diffusion radius ~$\sqrt{t}$. Our approach consists in a study of the generating function $G_{Δx}(λ) = Σ_n$ ${λ^n}p_n(Δx)$ for the probabilities $p_n(Δx)$ of observing $\textit{n}$ particles in an interval of given size $Δ\textit{x}$ from the lead particle to its left, fixing the position of the latter. This generating function can be expressed with the help of functions solving the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation with suitable initial conditions. In the infinite-time and large-$Δ\textit{x}$ limits, we find that the mean number of particles in the interval grows exponentially with $Δ\textit{x}$, and that the generating function obeys a nontrivial scaling law, depending on $Δ\textit{x}$ and λ through the combined variable $[Δx — f(λ)]^3 / Δx^2$, where $\textit{f}$(λ) ≡ – ln(1 – λ) – ln [– ln(1 – λ)]. From this property, one may conjecture that the growth of the typical particle number with the size of the interval is slower than exponential, but, surprisingly enough, only by a subleading factor at large Δ$\textit{x}$. The scaling we argue is consistent with results from a numerical integration of the FKPP equation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion-sign problem is crucial, as the signal-to-noise ratio decreases exponentially with Euclidean time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parameterized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. We only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross-sections of larger nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Million-scale data integrated deep neural network for phonon properties of heuslers spanning the periodic table

Existing machine learning potentials for predicting phonon properties of crystals are typically limited on a material-to-material basis, primarily due to the exponential scaling of model complexity with the number of atomic species. We address this bottleneck with the developed Elemental Spatial Density Neural Network Force Field, namely Elemental-SDNNFF. The effectiveness and precision of our Elemental-SDNNFF approach are demonstrated on 11,866 full, half, and quaternary Heusler structures spanning 55 elements in the periodic table by prediction of complete phonon properties. Self-improvement schemes including active learning and data augmentation techniques provide an abundant 9.4 million atomic data for training. Deep insight into predicted ultralow lattice thermal conductivity (<1 Wm –1 K –1 ) of 774 Heusler structures is gained by p–d orbital hybridization analysis. Additionally, a class of two-band charge-2 Weyl points, referred to as “double Weyl points”, are found in 68% and 87% of 1662 half and 1550 quaternary Heuslers, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Laser-driven ion and electron acceleration from near-critical density gas targets: Towards high-repetition rate operation in the 1 PW, sub-100 fs laser interaction regime

Ion acceleration from gaseous targets driven by relativistic-intensity lasers was demonstrated as early as the late 1990s, yet most of the experiments conducted to date have involved picosecond-duration, Nd:glass lasers operating at low repetition rate. Here, we present measurements on the interaction of ultraintense ( ≈ 10 20 W cm − 2 , 1 PW ) , ultrashort ( ≈ 70 fs ) Ti:Sa laser pulses with near-critical ( ≈ 10 20 cm − 3 ) helium gas jets, a debris-free targetry with the potential for future compatibility with high ( ≈ 1 Hz ) repetition rate operation. We provide evidence of α particles being forward accelerated up to ≈ 2.7 − MeV energy with a total flux of ≈ 10 11 sr − 1 as integrated over > 0.1 − MeV energies and detected within a 0.5 − mrad solid angle. We also report on on-axis emission of relativistic electrons with an exponentially decaying spectrum characterized by a ≈ 10 − MeV slope, i.e., five times larger than the standard ponderomotive scaling. The total charge of these electrons with energy above 2 MeV is estimated to be of ≈ 1 nC , corresponding to ≈ 0.1 % of the laser drive energy. In addition, we observe the formation of a plasma channel, extending longitudinally across the gas density maximum and expanding radially with time. These results are well captured by large-scale particle-in-cell simulations, which reveal that the detected fast ions most likely originate from reflection off the rapidly expanding channel walls. The latter process is predicted to yield ion energies in the MeV range, which compare well with the measurements. Finally, direct laser acceleration is shown to be the dominant mechanism behind the observed electron energization. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Chaotic roots of the modular multiplication dynamical system in Shor's algorithm

Shor's factoring algorithm, believed to provide an exponential speedup over classical computation, relies on finding the period of an exactly periodic quantum modular multiplication operator. This exact periodicity is the hallmark of an integrable system, which is paradoxical from the viewpoint of quantum chaos, given that the classical limit of the modular multiplication operator is a highly chaotic system that occupies the “maximally random” Bernoulli level of the classical ergodic hierarchy. In this work, we approach this apparent paradox from a quantum dynamical systems viewpoint, and consider whether signatures of ergodicity and chaos may indeed be encoded in such an “integrable” quantization of a chaotic system. We show that Shor's modular multiplication operator, in specific cases, can be written as a superposition of quantized 𝐴-baker's maps exhibiting more typical signatures of quantum chaos and ergodicity. This work suggests that the integrability of Shor's modular multiplication operator may stem from the interference of other “chaotic” quantizations of the same family of maps, and paves the way for deeper studies on the interplay of integrability, ergodicity, and chaos in and via quantum algorithms.

Dynamical systems↗
DOE OSTI2023

HQ-Sim: High-performance State Vector Simulation of Quantum Circuits on Heterogeneous HPC Systems

Quantum circuit simulations are applied in more and more circum- stances as the quantum computing community becomes broader. It helps researchers to evaluate quantum algorithms and relieve the burden of limited quantum computing resources. However, most of the state-of-the-art quantum simulators utilize either CPU or GPU to store and calculate the state vector, which results in resource starvation. Moreover, the maximum number of qubits supported by the simulator is bounded by the memory, since the memory utilization increases exponentially with the number of qubits. In this study, we leverage Heterogeneous computing to utilize both CPU and GPU to store and update state vectors. We also integrate lossy data compression to reduce memory requirements. Specifically, we develop a heterogeneous framework that has a dynamic scheduler to fully utilize the computing resources. We apply lossy compression to chunked state vector to make the maximum number of qubits higher than the regular simulators, the compression also benefits the data movement between CPU and GPU.

Zhang, Boyuan↗
DOE OSTI2023

Integrating fine root diameter and watershed mapping to characterize rhizosphere hydrology

Root morphology and soil hydraulic characteristics were integrated using watershed distance mapping to show water distribution and uptake across the plant-soil interface. Poplar (Populus deltoides, P. trichocarpa), maize (Zea mays), juniper (Juniperus virginiana), grape (Vitis rotundifolia) and maple (Acer saccharum) seedlings were grown in sand, after which root diameter and soil water dynamics were assessed via sequential neutron radiography. Three local soil regions (root-soil interface or edge, rhizosphere, bulk soil) were classified based on both radial distance from the root surface and diameter of the nearest root, from which changes in water content and distribution were characterized using digital image processing. Water content dynamics across the rhizosphere showed two different species-independent processes: a consistently elevated water content at the root-soil edge interface which increased with root diameter, and hysteresis as the rhizosphere transitioned to bulk soil (~0.5 cm from the root), independent of root diameter. Water uptake per unit root surface area declined exponentially with root diameter, independent of species. Results highlight the species-independent hydrologic characteristics of the rhizosphere and the potential for evaluating them in a local spatially connected soil context. Avenues for improved integration of soil and root characteristics are discussed.

59 BASIC BIOLOGICAL SCIENCES↗