Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “performance bound”

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 199 records · Page 11

Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-like Calculations

Here, this Letter introduces a unified emulation framework for studying continuum physics in finite quantum systems. Using a reduced basis method, we construct powerful emulators for the inhomogeneous Schrödinger equation that operate in a combined parameter space of complex energy (𝐸) and other inputs (𝜽). Within the space, the emulators simultaneously perform analytical continuation in 𝐸—extracting continuum physics from numerically simpler bound-state-like calculations—and interpolate this entire process across 𝜽. This yields a small, non-Hermitian system whose properties (e.g., resonances and scattering observables) can be rapidly predicted for any 𝜽. Crucially, the complex-𝐸 emulation provides a pathway to compute continuum observables for complex systems where advanced bound-state methods exist but direct continuum calculations are yet to be developed, while the 𝜽 emulation enables rapid parameter-space exploration and can be adapted to accelerate other existing continuum calculations. Demonstrations with two- and three-body systems highlight the method’s effectiveness and suggest its connection to (near-)optimal rational approximation. This Letter presents the key results, with further details reserved for a companion paper.

ab initio calculations↗

Julia as a unifying end-to-end workflow language on the Frontier exascale system

We evaluate Julia as a single language and ecosystem paradigm powered by LLVM to develop workflow components for high-performance computing. We run a Gray-Scott, 2-variable diffusion-reaction application using a memory-bound, 7-point stencil kernel on Frontier, the US Department of Energy’s first exascale supercomputer. We evaluate the performance, scaling, and trade-offs of (i) the computational kernel on AMD’s MI250x GPUs, (ii) weak scaling up to 4,096 MPI processes/GPUs or 512 nodes, (iii) parallel I/O writes using the ADIOS2 library bindings, and (iv) Jupyter Notebooks for interactive analysis. Results suggest that although Julia generates a reasonable LLVM-IR, a nearly 50% performance difference exists vs. native AMD HIP stencil codes when running on the GPUs. As expected, we observed near-zero overhead when using MPI and parallel I/O bindings for system-wide installed implementations. Consequently, Julia emerges as a compelling high-performance and high-productivity workflow composition language, as measured on the fastest supercomputer in the world.

Godoy, William↗

Frameshifting Stimulatory Sequence Induces Large Structural Change of Ribosomal Proteins When Bound to E. coli Ribosomes

Biological macromolecular machines occupy a continuum of structural conformations to perform cellular tasks. Mapping this conformational space provides an insight into its functionality. While the cryo-electron microscopy resolution revolution has expanded our ability to characterize the conformational continuums, there are obstacles in structurally characterizing regions of high flexibility. These technical barriers have impeded characterization of flexible ribosomal proteins when the ribosome is interacting with mRNA stem-loop structures such as a frameshifting stimulatory sequence (FSS). Small-angle neutron/X-ray scattering and electron microscopy were used to study ribosomal samples and compared structural differences between a ribosome that is bound to an FSS stem-loop compared to a ribosome bound to linear mRNA. This comparison shows that a large protein stalk elongates by 22% when the 70S interacts with an mRNA stem-loop. Finally, our results suggest that ribosomal proteins have extensive flexibility and may influence important ribosomal mechanisms, such as those that involve FSS.

36 MATERIALS SCIENCE↗

Improving Progressive Retrieval for HPC Scientific Data using Deep Neural Network

As the disparity between compute and I/O on high-performance computing systems has continued to widen, it has become increasingly difficult to perform post-hoc data analytics on full-resolution scientific simulation data due to the high I/O cost. Error-bounded data decomposition and progressive data retrieval framework has recently been developed to address such a challenge by performing data decomposition before storage and reading only part of the decomposed data when necessary. However, the performance of the progressive retrieval framework has been suffering from the over-pessimistic error control theory, such that the achieved maximum error of recomposed data is significantly lower than the required error. Therefore, more data than required is fetched for recomposition, incurring additional I/O overhead. In order to tackle this issue, we propose a DNN-based progressive retrieval framework that can better identify the minimum amount of data to be retrieved. Our contributions are as follows: 1) We provide an in-depth investigation of the recently developed progressive retrieval framework; 2) We propose two designs of prediction models (named D-MGARD and E-MGARD) to estimate the amount of retrieved data size based on error bounds. 3) We evaluate our proposed solutions using scientific datasets generated by real-world simulations from two domains. Evaluation results demonstrate the effectiveness of our solution in accurately predicting the amount of retrieval data size, as well as the advantages of our solution over the traditional approach to reducing the I/O overhead. Based on our evaluation, our solution is shown to read significantly less data (5% - 40% with D-MGARD, 20% - 80% with E-MGARD).

Wang, Jinzhen↗

Universal Optimization and Tighter Fidelity Bounds for Approximate Quantum Error Correction

Approximate quantum error correction (AQEC) not only dictates the performance of discrete- and continuous-variable quantum error correction codes but also serves as a unifying framework across various physical disciplines. Identifying the optimal recovery channel to maximize the entanglement fidelity via standard semidefinite programming is computationally bottlenecked by the exponentially growing number of Kraus operators with system size, rendering large-scale optimization prohibitive. While analytical near-optimal maps exist, they typically work only when the Knill-Laflamme conditions are nearly satisfied. In this Letter, we establish an efficient framework by leveraging the duality between recovery and environment decoupling. This framework yields a tighter analytical lower bound on entanglement fidelity than the conventional limit set by the transpose channel. Furthermore, by exploiting the decayed weights of noise Kraus operators, we introduce a framework based on principal component analysis to reduce the dimension. In thermal loss channels where the weights decay exponentially, our approach achieves a 33-fold computational speedup while maintaining rigorous accuracy. Our approach enables high-precision optimization for AQEC codes that were previously intractable due to the curse of dimensionality.

Wu, Jing [Fermilab] (ORCID:0000000249460732)↗

SPRUCE Climate Warming and Elevated CO2 Rapidly Alter Peatland Soil Carbon Sources and Stability: Supporting Data

This data set reports a suite of complementary biogeochemical analyses of peat samples from the SPRUCE (Spruce and Peatland Responses Under Changing Environments) experiment. Results were collected using quantitative molecular analysis of bulk soil carbon to assess the stability of soil organic carbon following whole-ecosystem warming and exposure to elevated carbon dioxide concentrations (eCO2). Targeted soil organic carbon components include solvent-extractable compounds (alkanoic acids, alkanols, alkanes, steroids, and terpenoids), ester-bound hydrolysable biopolymers (cutin and suberin markers), lignin phenols, and pyrogenic carbon. Bulk peat samples were analysed by Soxhlet extraction and solid phase separation for solvent-extractable compounds, alkaline hydrolysis to extract hydrolysable biopolymers, copper (II) oxide oxidation to extract lignin phenols and benzene polycarboxylic acids (BPCAs) as an approximation of pyrogenic carbon. Samples were analysed by gas chromatography (GC) equipped with a flame ionization detector (GC-FID) and compound identification was performed on GC coupled to mass selective detector (MS) for solvent-extractable compounds, ester-bound hydrolysable biopolymers and lignin phenols, and high-performance liquid chromatograph (HPLC) for pyrogenic carbon. Results are presented in Ofiti et al. (accepted). The experimental work was conducted on samples collected in August 2018 at the SPRUCE climate manipulation experiment in northern Minnesota, 40 km north of Grand Rapids, in the USDA Forest Service Marcell Experimental Forest (MEF). Samples were collected and later analysed in a 10 cm increments over 0 to 50 cm depth and 25 cm intervals from 50 to 75 cm. Samples were analyzed for lignin phenols over 0 to 30 cm depth. This data set contains one file in comma separate (*.csv) format. This dataset contains data used to produce: Ofiti, N.O.E., Schmidt, M.W.I., Abiven, S., Hanson, P.J., Iversen, C.M., Wilson, R.M., Kostka, J.E., Wiesenberg, G.L.B., Malhotra, A. 2023. Climate warming and elevated CO2 rapidly alter peatland soil carbon sources and stability. Nat Commun 14, 7533. https://doi.org/10.1038/s41467-023-43410-z.

SPRUCE experiment, Marcell Experimental Forest, so↗

Computing Bottleneck Structures at Scale for High-Precision Network Performance Analysis

The Theory of Bottleneck Structures is a recently-developed framework for studying the performance of data networks. It describes how local perturbations in one part of the network propagate and interact with others. This framework is a powerful analytical tool that allows network operators to make accurate predictions about network behavior and thereby optimize performance. Previous work implemented a software package for bottleneck structure analysis, but applied it only to toy examples. In this work, we introduce the first software package capable of scaling bottleneck structure analysis to production-size networks. Here, we benchmark our system using logs from ESnet, the Department of Energy's high-performance data network that connects research institutions in the U.S. Using the previously published tool as a baseline, we demonstrate that our system achieves vastly improved performance, constructing the bottleneck structure graphs in 0.21 s and calculating link derivatives in 0.09 s on average. We also study the asymptotic complexity of our core algorithms, demonstrating good scaling properties and strong agreement with theoretical bounds. These results indicate that our new software package can maintain its fast performance when applied to even larger networks. They also show that our software is efficient enough to analyze rapidly changing networks in real time. Overall, we demonstrate the feasibility of applying bottleneck structure analysis to solve practical problems in large, real-world data networks.

benchmark↗

Communication Lower Bounds and Optimal Algorithms for Multiple Tensor-Times-Matrix Computation

Multiple tensor-times-matrix (Multi-TTM) is a key computation in algorithms for computing and operating with the Tucker tensor decomposition, which is frequently used in multidimensional data analysis. Here, we establish communication lower bounds that determine how much data movement is required (under mild conditions) to perform the Multi-TTM computation in parallel. The crux of the proof relies on analytically solving a constrained, nonlinear optimization problem. We also present a parallel algorithm to perform this computation that organizes the processors into a logical grid with twice as many modes as the input tensor. We show that, with correct choices of grid dimensions, the communication cost of the algorithm attains the lower bounds and is therefore communication optimal. Finally, we show that our algorithm can significantly reduce communication compared to the straightforward approach of expressing the computation as a sequence of tensor-times-matrix operations when the input and output tensors vary greatly in size.

HBL-inequalities↗

Bose–Einstein condensation of a two-magnon bound state in a spin-1 triangular lattice

In ordered magnets, the elementary excitations are spin waves (magnons), which obey Bose–Einstein statistics. Similarly to Cooper pairs in superconductors, magnons can be paired into bound states under attractive interactions. The Zeeman coupling to a magnetic field is able to tune the particle density through a quantum critical point, beyond which a ‘hidden order’ is predicted to exist. Here, in this work, we report direct observation of the Bose–Einstein condensation of the two-magnon bound state in Na 2 BaNi(PO 4 ) 2 . Comprehensive thermodynamic measurements confirmed the two-dimensional Bose–Einstein condensation quantum critical point at the saturation field. Inelastic neutron scattering experiments were performed to establish the microscopic model. An exact solution revealed stable two-magnon bound states that were further confirmed by electron spin resonance and nuclear magnetic resonance experiments, demonstrating that the quantum critical point is due to the pair condensation, and the phase below the saturation field is likely the long-sought-after spin nematic phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct functionalization of C-H bonds by electrophilic anions

Bond formation between two of the most inert molecules, alkanes and [B12X12]2- (X = Cl, Br), is performed in a two step process. Fragmentation of [B12X12]2- in the gas phase generates highly reactive [B12X11]- ions which spontaneously react with alkanes. The reaction mechanism was investigated using tandem mass spectrometry and gas-phase vibrational spectroscopy combined with electronic structure calculations. This effort revealed the unprecedented and counterintuitive electrophilic substitution of a proton in an alkane by a negatively charged ion resulting in a B-C bond formation. The product is a dianionic [B12X11CnH2n+1]2- species, in which H+ is electrostatically bound to the dianion. High flux ion soft-landing was performed to co-deposit [B12X11]- and complex organic molecules (phthalates) in thin layers on surfaces. Molecular structure analysis of the product films revealed that C-H functionalization by [B12X11]- occured in the presence of other more reactive functional groups. This observation demonstrates the utility of highly reactive fragment ions for defined bond formation processes and may pave the way to the use of gas phase ion chemistry for chemical synthesis in the condensed phase.

Warneke, Jonas↗

Precision three-dimensional imaging of nuclei using recoil-free jets

In this study, we explore the azimuthal angle decorrelation of lepton-jet pairs in e-p and e-A collisions as a means for precision measurements of the three-dimensional structure of bound and free nucleons. Utilizing soft-collinear effective theory, we perform the first-ever resummation of this process in e-p collisions at NNLL accuracy using a recoil-free jet axis. Our results are validated against Pythia simulations. In e-A collisions, we address the complex interplay between three characteristic length scales: the medium length L, the mean free path of the energetic parton in the medium λ, and the hadronization length L h . We demonstrate that in the thin-dilute limit, where L $\ll$ L h and L ~ λ, this process can serve as a robust probe of the three-dimensional structure for bound nucleons. We conclude by offering predictions for future experiments at the Electron-Ion Collider within this limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Groundwater and Intruder Radionuclide Screening

Solid Waste Management (SWM) operates the E-Area Low-Level Waste Facility (ELLWF) where low level on-site and off-site solid waste streams are buried. The facility has been in operation since late 1994 and is currently projected to remain in operation until 2065. This facility can accommodate a broad range of waste forms resulting from the six different types of disposal unit options (i.e., varying degrees of engineered barriers → trenches to concrete vaults). This facility is currently operating under a Performance Assessment (PA) issued back in 2008, along with several subsequent supporting Special Analyses (SA). The Savannah River National Laboratory (SRNL) developed the prior PAs and SAs and has been tasked to update the facility’s PA. For operating the E-Area facility, a Waste Inventory Tracking System (e.g., WITS) is actively employed by waste generators where every radionuclide entering the facility (to be buried in one of its many disposal units) must be either directly or indirectly accounted for. Since there is a large number of radionuclides in existence (>3,000), the International Commission on Radiological Protection (ICRP-107) has provided guidance on the subset of radionuclides requiring further assessment in landfills such as the E-Area LLWF. The ICRP-107 publication provides critical radiological information on 1,252 radionuclides of 97 elements. This database is the critical starting point for developing a consistent inventory limit system. The explicit measurement and tracking of all radionuclides are not necessary when process knowledge, burial history, and radiological aspects are factored into conservative groundwater and intruder screening processes. Across the DOE complex these screening processes have been historically performed using the methodology suggested by the National Council on Radiation Protection and Measurements (NCRP) as presented in their original report (NCRP 1984) and then refined in a later report (NCRP 1996). In the recommended screening models employed within this report the traditional NCRP models are updated to better handle progeny and better reflect the known characteristics of E-Area. This improved screening process is referred to as the “NCRP-like” method. Various upgrades to the traditional NCRP methodology have been used by others (e.g., NRC funded effort by Kennedy and Strenge 1992). The more detailed models by Kennedy and Strenge (1992) were considered in this effort and are compared to the results from the more traditional NCRP-like models. The groundwater and inadvertent intruder screening analyses presented in this report start with this 1,252 radionuclide list and reduces it down to more manageable lists that are applicable to the various disposal unit types contained within E-Area. In order to reduce this starting list, some level of exposure risk must be considered acceptable. Historically, a dose (or concentration level) has been compared with a screening criterion set to 1% of a performance measure (e.g., a beta-gamma dose not to exceed the 4 mrem/yr beta-gamma performance measure x 0.01 = 0.04 mrem/yr). Thus, if a radionuclide produced a bounding or screening-level dose (or concentrations) less than the screening criterion, it could be safely removed from further consideration. Several tiers of screening and bounding level analyses have been considered in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Reduced order models for Lagrangian hydrodynamics

It is reported, as a mathematical model of high-speed flow and shock wave propagation in a complex multimaterial setting, Lagrangian hydrodynamics is characterized by moving meshes, advection-dominated solutions, and moving shock fronts with sharp gradients. These challenges hinder the existing projection-based model reduction schemes from being practical. We develop several variations of projection-based reduced order model techniques for Lagrangian hydrodynamics by introducing three different reduced bases for position, velocity, and energy fields. A time-windowing approach is also developed to address the challenge imposed by the advection-dominated solutions. Lagrangian hydrodynamics is formulated as a nonlinear problem, which requires a proper hyper-reduction technique. Therefore, we apply the over-sampling DEIM and SNS approaches to reduce the complexity due to the nonlinear terms. Finally, we also present both a posteriori and a priori error bounds associated with our reduced order model. We compare the performance of the spatial and time-windowing reduced order modeling approaches in terms of accuracy and speed-up with respect to the corresponding full order model for several numerical examples, namely Sedov blast, Gresho vortices, Taylor-Green vortices, and triple-point problems.

97 MATHEMATICS AND COMPUTING↗

Leveraging interpolation models and error bounds for verifiable scientific machine learning

Effective verification and validation techniques for modern scientific machine learning workflows are challenging to devise. Statistical methods are abundant and easily deployed, but often rely on speculative assumptions about the data and methods involved. Error bounds for classical interpolation techniques can provide mathematically rigorous estimates of accuracy, but often are difficult or impractical to determine computationally. Here, in this work, we present a best-of-both-worlds approach to verifiable scientific machine learning by demonstrating that (1) multiple standard interpolation techniques have informative error bounds that can be computed or estimated efficiently; (2) comparative performance among distinct interpolants can aid in validation goals; (3) deploying interpolation methods on latent spaces generated by deep learning techniques enables some interpretability for black-box models. We present a detailed case study of our approach for predicting lift-drag ratios from airfoil images. Code developed for this work is available in a public Github repository.

97 MATHEMATICS AND COMPUTING↗

Surprising olefin/paraffin separation performance recovery of highly aged carbon molecular sieve hollow fiber membranes by a super-hyperaging treatment

Carbon molecular sieve (CMS) hollow fibers derived from 6FDA 1 :BPDA 1 /DAM 2 , a well-studied polyimide precursor with attractive propylene/propane separation performance, were purposely treated to simulate long time aging under storage conditions. A simple, efficient post-aging thermal-treatment regeneration termed “super-hyperaging” is shown here to recover C 3 H 6 permeance with satisfactory C 3 H 6 /C 3 H 8 selectivity. A super-hyperaging temperature around 280 °C in air is shown to provide an optimum tradeoff in permeance and selectivity for highly aged samples. The regenerated CMS membrane displays performance well beyond the polymer permeance-selectivity tradeoff upper bound, and a hypothetical mechanism for the super-hyperaging recovery of aged CMS fibers is discussed in this work.

Recovery↗

Quantifying the generalization error in deep learning in terms of data distribution and neural network smoothness

We report the accuracy of deep learning, i.e., deep neural networks, can be characterized by dividing the total error into three main types: approximation error, optimization error, and generalization error. Whereas there are some satisfactory answers to the problems of approximation and optimization, much less is known about the theory of generalization. Most existing theoretical works for generalization fail to explain the performance of neural networks in practice. To derive a meaningful bound, we study the generalization error of neural networks for classification problems in terms of data distribution and neural network smoothness. We introduce the cover complexity (CC) to measure the difficulty of learning a data set and the inverse of the modulus of continuity to quantify neural network smoothness. A quantitative bound for expected accuracy/error is derived by considering both the CC and neural network smoothness. Although most of the analysis is general and not specific to neural networks, we validate our theoretical assumptions and results numerically for neural networks by several data sets of images. The numerical results confirm that the expected error of trained networks scaled with the square root of the number of classes has a linear relationship with respect to the CC. We also observe a clear consistency between test loss and neural network smoothness during the training process. In addition, we demonstrate empirically that the neural network smoothness decreases when the network size increases whereas the smoothness is insensitive to training dataset size.

97 MATHEMATICS AND COMPUTING↗

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗

Primary observables for gauge boson collider signals

In this paper, we determine a basis for the on-shell four-point amplitudes V V V V for Standard Model gauge bosons V = W ± , Z , γ , g . Following previous work, this completes the analysis of three- and four-point amplitudes for the Standard Model and could be used for model-independent searches for beyond the Standard Model physics at colliders. Our results include a Lagrangian parametrization for the “primary” amplitudes, where including additional derivatives leads to the Mandelstam “descendant” amplitudes, and upper bounds on the coupling constants from imposing unitarity. We also perform an estimate for the sensitivity for new Z decays at the High-Luminosity LHC, finding that Z → γ ℓ ¯ ℓ could be searched for, but that other decay modes, like Z → ( γ γ γ , γ g g ) , are too small to be discovered after imposing unitarity constraints. Published by the American Physical Society 2024

Astronomy & Astrophysics↗