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

An atomistic theory of nucleation: Self-organization via non-equilibrium work and fluctuations

For this work, insights from non-equilibrium statistical mechanics, highlighting the role of work and fluctuations at the microscale, are applied toward the development of a fundamental, rigorous and purely atomistic theory of nucleation. Nanoscale fluctuations in order, density and heat influence the local nucleation rate by orders of magnitude, necessitating their inclusion through a modern approach. Coarse-graining over the underlying Hamiltonian dynamics allows derivation of a microscale expression for the nucleation rate in terms of a classical path integral over far from equilibrium trajectories and their associated work. Second law violating states at the microscale, as found from the dynamics of small critical nucleation clusters, contribute exponentially to the observable macroscale nucleation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Soft factorisation and exponentiation from Schwinger-space geometry

Infrared divergences in Quantum Field Theory govern the low-energy dynamics of many physical theories, and their understanding is a crucial ingredient in predicting the outcomes of collider experiments. We present a novel approach to deriving the structure of these divergences by employing the Schwinger parametrization of Feynman integrals. After using tropical geometry to identify divergent limits, we study the all-orders asymptotic properties of Feynman diagrams via matrix manipulations of graph Laplacians, which allows us to analyse their IR behaviour systematically. We explicitly demonstrate the soft-hard factorization of the integrand for a broad class of diagrams, and reveal that when written in terms of worldline distances, topologically distinct diagrams asymptote to the same integrand at leading order in the soft limit. In particular, for the case of Quantum Electrodynamics (with massive fermions), we use this fact to show how ladder-type diagrams combine in Schwinger-parameter space to yield the correct exponentiated soft anomalous dimension. This framework provides a foundation for extending these methods to more complex theories like Quantum Chromodynamics and offers a pathway towards a systematic understanding of infrared divergences in perturbative amplitudes.

Factorization↗
DOE OSTI2020

Sudden-Discharge Cycling Characteristics and Millisecond Dynamic Behaviors of a HTS Stainless-Steel Insulated Double-Pancake Coil With Thin Copper Plates

We have tested an intermediate-sized HTS stainless-steel (co-wound) insulated double pancake coil (~132 turns per coil) and discovered dynamic effects during sudden discharge on both millisecond and second scales. For soft break, the power supply is turned off suddenly but the shunt resistor (164 milliohm) is still connected; for hard break, a high voltage contactor is used to suddenly open the circuit so that the leads are completely open. Hard break is dangerous to traditionally insulated LTS magnets, but the SS-insulated HTS DPC retained integrity over ¡80 charge-discharge cycles. Voltage decay curves for soft and hard breaks are studied on short and long timescales. On long timescales, we observe the expected exponential decay and measure the time constant to calculate contact resistance. On short timescales, we observe highly dynamic effects: a hyper-exponential decay suggesting that the inductance or contact resistance is suddenly changing with time, likely an indication of the RL circuit forming. Moreover, peak voltage immediately (ms) after the sudden discharge were collected and compared. It is found that the peak voltage in soft-discharge is proportional to the applied current, while the those in hard-break has a unique pattern and tends to stay flat even though the applied current increases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Light cone dynamics in excitonic states of two-component Bose and Fermi gases

Herein we consider the non-equilibrium dynamics of two-component one dimensional quantum gases in the limit of extreme population imbalance where the minority species has but a single particle. We consider the situation where the gas is prepared in a state with a single spatially localized exciton: the single particle of the minority species is spatially localized while the density of the majority species in the vicinity of the minority particle sees a depression. Remarkably, we are able to consider cases where the gas contains on the order of $N$ = 100 particles, comparable to that studied in experiments on cold atomic gases. We are able to do so by exploiting the integrability of the gas together with the observation that the excitonic state can be constructed through a simple superposition of exact eigenstates of the gas. The number of states in this superposition, rather than being exponentially large in the number of particles, scales linearly with $N$. We study the evolution of such spatially localized states in both strongly interacting Bose and Fermi gases. The behavior of the light cones when the interaction strength and density of the gas is varied can be understood from exact results for the spin excitation spectrum in these systems. We argue that the light cone in both cases exhibits scaling collapse. However unique to the Bose gas, we show that the presence of gapped finite-momentum roton-like excitations provide the Bose gas dynamics with secondary light cones.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Effective description of sub-maximal chaos: stringy effects for SYK scrambling

It has been proposed that the exponential decay and subsequent power law saturation of out-of-time-order correlation functions can be universally described by collective ‘scramblon’ modes. We develop this idea from a path integral perspective in several examples, thereby establishing a general formalism. After reformulating previous work on the Schwarzian theory and identity conformal blocks in two-dimensional CFTs relevant for systems in the infinite coupling limit with maximal quantum Lyapunov exponent, we focus on theories with sub-maximal chaos: we study the large-q limit of the SYK quantum dot and chain, both of which are amenable to analytical treatment at finite coupling. In both cases we identify the relevant scramblon modes, derive their effective action, and find bilocal vertex functions, thus constructing an effective description of chaos. The final results can be matched in detail to stringy corrections to the gravitational eikonal S-matrix in holographic CFTs, including a stringy Regge trajectory, bulk to boundary propagators, and multi-string effects that are unexplored holographically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A Quantum-Classical Collaborative Training Architecture Based on Quantum State Fidelity

Recent advancements have highlighted the limitations of current quantum systems, particularly the restricted number of qubits available on near-term quantum devices. This constraint greatly inhibits the range of applications that can leverage quantum computers. Moreover, as the available qubits increase, the computational complexity grows exponentially, posing additional challenges. Consequently, there is an urgent need to use qubits efficiently and mitigate both present limitations and future complexities. To address this, existing quantum applications attempt to integrate classical and quantum systems in a hybrid framework. In this study, we concentrate on quantum deep learning and introduce a collaborative classical-quantum architecture called co-TenQu. The classical component employs a tensor network for compression and feature extraction, enabling higher-dimensional data to be encoded onto logical quantum circuits with limited qubits. On the quantum side, we propose a quantum-state-fidelity-based evaluation function to iteratively train the network through a feedback loop between the two sides. co-TenQu has been implemented and evaluated with both simulators and the IBM-Q platform. Compared to state-of-the-art approaches, co-TenQu enhances a classical deep neural network by up to 41.72% in a fair setting. Additionally, it outperforms other quantum-based methods by up to 1.9 times and achieves similar accuracy while utilizing 70.59% fewer qubits.

42 ENGINEERING↗
DOE OSTI2021

Replication Data for: The effect of frustrated rotations on the pre-exponential factor for unimolecular reactions on surfaces: a case study of alkoxy dehydrogenation

The data underlying this published work have been made publicly available in this repository as part of the IMASC Data Management Plan. This work was supported as part of the Integrated Mesoscale Architectures for Sustainable Catalysis (IMASC), an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences under Award # DE-SC0012573.

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

Lyapunov growth in quantum spin chains

The Ising spin chain with longitudinal and transverse magnetic fields is often used in studies of quantum chaos, displaying both chaotic and integrable regions in its parameter space. However, even at a strongly chaotic point this model does not exhibit Lyapunov growth of the commutator squared of spin operators, as this observable saturates before exponential growth can manifest itself (even in situations where a spatial suppression factor makes the initial commutator small). In this work, we extend this model from the spin 1/2 Ising model to higher spins, demonstrate numerically that a window of exponential growth opens up for sufficiently large spin, and extract a quantity which corresponds to a notion of a Lyapunov exponent. In the classical infinite-spin limit, we identify and compute the appropriate classical analog of the commutator squared, and show that the corresponding exponent agrees with the infinite-spin limit extracted from the quantum spin chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quark and gluon two-loop beam functions for leading-jet p T and slicing at NNLO

We compute the complete set of two-loop beam functions for the transverse momentum distribution of the leading jet produced in association with an arbitrary colour-singlet system. Our results constitute the last missing ingredient for the calculation of the jet-vetoed cross section at small veto scales at the next-to-next-to-leading order, as well as an important ingredient for its resummation to next-to-next-to-next-to-leading logarithmic order. Our calculation is performed in the soft-collinear effective theory framework with a suitable regularisation of the rapidity divergences occurring in the phase-space integrals. We discuss the occurrence of soft-collinear mixing terms that might violate the factorisation theorem, and demonstrate that they are naturally absorbed into the beam functions at two loops in the exponential rapidity regularisation scheme when performing a multipole expansion of the measurement function. As in our recent computation of the two-loop soft function, we present the results as a Laurent expansion in the jet radius R. We provide analytic expressions for all flavour channels in x space with the exception of a set of R-independent non-logarithmic terms that are given as numerical grids. We also perform a fully numerical calculation with exact R dependence, and find that it agrees with our analytic expansion at the permyriad level or better. Our calculation allows us to define a next-to-next-to-leading order slicing method using the leading-jet p T as a slicing variable. As a check of our results, we carry out a calculation of the Higgs and Z boson total production cross sections at the next-to-next-to-leading order in QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Anomalous losses of energetic particles in the presence of an oscillating radial electric field in fusion plasmas

The confinement of energetic particles in nuclear fusion devices is studied in the presence of an oscillating radial electric field and an axisymmetric magnetic equilibrium. Here, it is shown that, despite the poloidal and toroidal symmetries, initially integrable orbits turn into chaotic regions that can potentially intercept the wall of the tokamak, leading to particle losses. It is observed that the losses exhibit algebraic time decay different from the expected exponential decay characteristic of radial diffusive transport. A dynamical explanation of this behaviour is presented, within the continuous time random walk theory. The central point of the analysis is based on the fact that, contrary to the radial displacement, the poloidal angle is not bounded and a proper statistical analysis can therefore be made, showing for the first time that energetic particle transport can be super-diffusive in the poloidal direction and characterised by asymmetric poloidal displacement. The connection between poloidal and radial positions ensured by the conservation of the toroidal canonical momentum, implies that energetic particles spend statistically more time in the inner region of the tokamak than in the outer one, which explains the observed algebraic decay. This indicates that energetic particles might be efficiently slowed down by the thermal population before leaving the system. Also, the asymmetric transport reveals a new possible mechanism of self-generation of momentum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Gauss-Radau-Laguerre Discrete Variable Representation for Use in Continuum Electron Dynamics

In this work, we detail an implementation, suitable for calculations on highly correlated ionizing systems, of a modified finite element discrete variable representation (FE-DVR) appended with a Gauss-Radau-Laguerre element. The appended element includes exterior complex scaling (ECS) to impose outgoing wave boundary conditions on treatments of processes involving continuum electrons. In this “infinite range” ECS (irECS), the complications that introduce reflections from the end of the grid when the last ECS finite element has finite range are avoided by the use of the Laguerre-weighted exponentially decaying tails, while outgoing wave boundary conditions are still imposed via the ECS transformation. For highly correlated systems in the absence of strong external fields we find that accurate two-electron integrals are essential in this modified FE-DVR. To accurately compute the two-electron integrals over the entire ECS contour, we present a detailed examination of the implications from the boundary terms that arise in a solution of Poisson’s equation with the Radau-Laguerre basis. A boundary term correction is necessary, and when included, the Radau-Laguerre DVR can accurately describe highly correlated states such as the doubly excited states of helium over the entire ECS contour.

elements↗
DOE OSTI2021

Prediction of Transition-State Scaling Relationships and Universal Transition-State Vibrational and Entropic Correlations for Dehydrogenations

Linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) or transition-state scaling (TSS) relations aid with the prediction of electronic energies. However, temperature effects and pre-exponentials are often taken as constants across metal surfaces or a homologous series. Vibrational scaling relationships (VSRs) offer a way to determine such parameters. Transition-state VSRs (TSVSRs) between local minima and transition states of AH X (A = C, N, O) surface diffusions correlate with BEP relations and broaden to thermochemical property scaling. Using density functional theory, we extend TSVSRs to AHX dehydrogenation reactions on transition-metal surfaces, relating vibrational modes of local minima to transition states. We first predict the slopes of the TSS relations by incorporating bond angles using the Slater–Koster structure factors and hybridization through crystal orbital overlap population analysis and energy overlap integrals between adsorbates and metal surfaces. Additionally, we uncover universal thermochemical property scaling, enabling the estimation of entropies and temperature corrections to enthalpies across a homologous series. Here, we demonstrate both significant vibrational corrections in reactions with low intrinsic electronic barriers and considerable variation in the pre-exponential of a simple dehydrogenation reaction across metals and AH X adsorbates.

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

Optimizing Transportation Networks for E-Waste Reverse Logistics: A Multi-Modal Cost Allocation and Pricing Strategy

The exponential growth of electronic waste (e-waste) poses critical challenges for sustainable reverse logistics and transportation network optimization. This study develops a dual-channel transportation framework for e-waste logistics that integrates dynamic freight pricing, cost allocation mechanisms, and game-theoretic coordination. The model captures interactions between centralized hubs and distributed processing networks, accounting for freight rate elasticity, volume allocation, and capacity constraints. Using Stackelberg game theory and cost-sharing strategies, the framework optimizes transportation efficiency and profit distribution across logistics channels. Numerical simulations show that the dual-channel structure increases centralized hub profit by 226.8% compared to baseline single-channel operations, while boosting total transported volume by 1.2% and nearly doubling freight collector profit under cost-sharing. Scenario analyses across regional infrastructures reveal that network density, policy incentives, and logistics costs shape routing efficiency and profit allocation. These findings suggest that coordinated strategies combining dynamic pricing, targeted infrastructure investment, and strategic cost allocation are needed to design efficient, resilient, and regionally adaptable e-waste transportation systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

IER 501: Pulsed Neutron Die-Away Experiments at LLNL [Slides]

Pulsed neutron die away experiments occur as follows: (1) a pulse of neutrons is injected into a target material via a D-T neutron generator, (2) neutrons thermalize, (3) neutrons spatially equilibrate, and (4) exponential decay is measured in fundamental mode.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Validation of the two-region Rossi-alpha model for reflected assemblies

Nuclear criticality safety, nonproliferation and safeguards, emergency response, and stockpile stewardship utilize estimates of the k eff multiplication factor. The value of k eff cannot be directly measured, but it can be inferred from the prompt neutron period. One modality of measuring the prompt neutron period is the Rossi-alpha method, which is predicated on fitting a histogram of detection time differences due to the nonrandom temporal distribution of same-fission-chain neutrons. Recent works have developed the motivation and theory to expand traditional one-region point kinetic models resulting in one-exponential histogram fits to two-region models that result in two-exponential fits. This paper validates the new two-region model using organic scintillator measurements of copper-reflected weapons-grade plutonium (0:83 ≤ k eff ≤ 0:94) and high-density- polyethylene-reflected highly enriched uranium (0:73 ≤ k eff ≤ 0:95). Furthermore, the results show that more thermal systems have shorter prompt neutron periods in the core region due to increased induced-fission probabilities for moderated neutrons. Here, a new parameter introduced by the two- region model is also shown to be correlated to the amount of reflection, and may be used to infer assembly properties such as type and amount of reflector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A Low-Rank QTT-based Finite Element Method for Elasticity Problems

We present an efficient and robust numerical algorithm for solving the linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. By integrating Z-ordering and subdomain concatenation, our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This efficiency is maintained while still guaranteeing exponential convergence with respect to the number of degrees of freedom. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented. This includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method’s performance in the numerical approximation of three representative test cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Hidden quantum criticality and entanglement in quench dynamics

Entanglement exhibits universal behavior near the ground-state critical point where correlations are long ranged and the thermodynamic entropy is vanishing. On the other hand, a quantum quench imparts extensive energy and results in a build up of entropy, hence no critical behavior is expected at long times. In this work, we present a new paradigm in the quench dynamics of integrable spin chains which exhibit a ground-state order-disorder phase transition at a critical line. Specifically, we consider a quench along the critical line which displays a volume-law behavior of the entropy and exponentially decaying correlations; however, we show that quantum criticality is hidden in higher-order correlations and becomes manifest via measures such as the mutual information and logarithmic negativity. Furthermore, we showcase the scale invariance of the Rényi mutual information between disjoint regions as further evidence for genuine critical behavior. We attribute the emergent quantum criticality to the soft mode not getting excited in spite of the quench. Moreover, the results presented here are universal to models whose low-energy or long-wavelength dynamics are well described by a free-fermionic field theory. Our results are amenable to an experimental realization on different quantum simulator platforms, particularly the Rydberg simulators. Published by the American Physical Society 2024

Paul, Sanku↗
DOE OSTI2025

Open string amplitudes: singularities, asymptotics and new representations

Open string amplitudes at tree level have been studied for over fifty years. However, there is no known analytic form for general n-point amplitudes, and their conventional representation in terms of worldsheet integrals does not make many of their most basic physical properties manifest. Recently, a formulation of these amplitudes exposing the underlying “binary geometry” via the use of “u” variables has given us many insights into their basic features. In this paper, we initiate a systematic exploration of fundamental aspects of open string amplitudes from this new point of view. We begin by finding explicit expressions for the factorization of amplitudes at general massive levels, which are seen to be determined by products of lower-point massless amplitudes with shifted kinematics. We then study the asymptotic behavior when subsets of kinematic variables become large, delineating regimes with exponential (generalized hard scattering) and power-law (generalized Regge) behavior. We also give precise expressions for the asymptotics, which reveal another example of the recently observed property of factorization away from poles. We derive new recursion relations for the amplitude, which when repeatedly applied reduce to infinite series representations with a wider domain of convergence than the usual integral representations. For the five-point case, we present a new closed-form expression for the amplitude that for the first time gives its analytic continuation to all of kinematic space. We also discuss novel relations between amplitudes at different kinematic points following from the recently observed “split” factorizations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗