Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theoretical and computational techniques”

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

Data-driven discovery of dynamics from time-resolved coherent scattering

Coherent X-ray scattering (CXS) techniques are capable of interrogating dynamics of nano- to mesoscale materials systems at time scales spanning several orders of magnitude. However, obtaining accurate theoretical descriptions of complex dynamics is often limited by one or more factors—the ability to visualize dynamics in real space, computational cost of high-fidelity simulations, and effectiveness of approximate or phenomenological models. In this work, we develop a data-driven framework to uncover mechanistic models of dynamics directly from time-resolved CXS measurements without solving the phase reconstruction problem for the entire time series of diffraction patterns. Our approach uses neural differential equations to parameterize unknown real-space dynamics and implements a computational scattering forward model to relate real-space predictions to reciprocal-space observations. This method is shown to recover the dynamics of several computational model systems under various simulated conditions of measurement resolution and noise. Moreover, the trained model enables estimation of long-term dynamics well beyond the maximum observation time, which can be used to inform and refine experimental parameters in practice. Finally, we demonstrate an experimental proof-of-concept by applying our framework to recover the probe trajectory from a ptychographic scan. Our proposed framework bridges the wide existing gap between approximate models and complex data.

36 MATERIALS SCIENCE↗

Fast superconducting qubit control with subharmonic drives

Increasing the fidelity of single-qubit gates requires a combination of faster pulses and increased qubit coherence. However, these requirements can be contradictory. Additionally, increasing the drive power can heat the qubit’s environment and degrade coherence. In this work, we circumvent this issue and achieve rapid gates by pumping a transmon’s native Kerr at approximately one third of the qubit’s resonant frequency. The subharmonic Rabi rate of the process is proportional to applied drive amplitude cubed, allowing for rapid gates. In addition, we demonstrate that filtering can be used to protect the qubit’s coherence while performing rapid gates. Single qubit gates as short as 37.4 ns are demonstrated with fidelity of 99.91%. We present theoretical calculations indicating that drive induced multi-photon decay will not limit qubit lifetime; calculated power absorption also indicates that this technique could reduce cryostat heating for fast gates, a vital requirement for large-scale quantum computers.

Xia, Mingkang [Univ. of Pittsburgh, PA (United Sta↗

Reduced-basis method for few-body bound-state emulation

Recent advances in both theoretical and computational methods have enabled large-scale, precision calculations of the properties of atomic nuclei. With the growing complexity of modern nuclear theory, however, also comes the need for novel methods to perform systematic studies and quantify the uncertainties of models when confronted with experimental data. Here, this study presents an application of such an approach, the reduced basis method, to substantially lower computational costs by constructing a significantly smaller Hamiltonian subspace informed by previous solutions. Our method shows comparable efficiency and accuracy to other dimensionality reduction techniques on an artificial three-body bound system while providing a richer representation of physical information in its projection and training subspace. This methodological advancement can be applied in other contexts and has the potential to greatly improve our ability to systematically explore theoretical models and thus enhance our understanding of the fundamental properties of nuclear systems.

cluster models↗

Spin excitations and dynamics in 2D magnets: An overview of magnons and magnetic skyrmions

van der Waals magnetic materials open up exciting possibilities to investigate fundamental spin properties in low-dimensional systems and to build compact functional spintronic structures. This review focuses on the recent progress in two-dimensional (2D) magnets that explore beyond the homogenous magnetically ordered state, including magnons (spin waves), magnetic skyrmions, and complex magnetic domains. Properties of these spin and topology excitations in 2D magnets provide insights into spin-orbit interactions and other forms of coupling between electrons, phonons, and spin-dependent excitations. Such spin-based quasiparticles can also serve as information carriers for next-generation high-speed computing elements. Furthermore, we will first lay out the general theoretical basis of dynamical responses in magnetic systems, followed by detailed descriptions of experimental progress in magnons and spin textures (including magnetic domains and skyrmions). Discussion on the experimental techniques and future perspectives are also included.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nucleon valence quark distribution functions from Lattice QCD

We present results on the nucleon valence quark distribution extracted from Lattice QCD simulations, using a gauge ensemble of Nf=2+1 Wilson-Clover fermions with a pion mass of m?=350 MeV and lattice spacing of a=0.091 fm. We obtain reduced Ioffe Time Distributions (rITDs) by computing appropriate matrix elements on the lattice, and elaborate on the extraction of the desired quark distributions from the rITDs following the pseudo-PDF approach. A set of techniques are considered in order to ensure ground state dominance. Theoretical and experimental implications of our calculation are discussed.

Kallidonis, Christos↗

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗

Research in Dark Energy

This is the final report for the project aimed at developing techniques and understanding experimental data within the context of theoretical models, expressly in support of the Dark Energy Spectroscopic Instrument (DESI).

79 ASTRONOMY AND ASTROPHYSICS↗

Error statistics and scalability of quantum error mitigation formulas

Quantum computing promises advantages over classical computing in many problems. Nevertheless, noise in quantum devices prevents most quantum algorithms from achieving the quantum advantage. Quantum error mitigation provides a variety of protocols to handle such noise using minimal qubit resources. While some of those protocols have been implemented in experiments for a few qubits, it remains unclear whether error mitigation will be effective in quantum circuits with tens to hundreds of qubits. In this paper, we apply statistics principles to quantum error mitigation and analyse the scaling behaviour of its intrinsic error. We find that the error increases linearly O(ϵN) with the gate number N before mitigation and sublinearly O(ϵ'N γ ) after mitigation, where γ ≈ 0.5, ϵ is the error rate of a quantum gate, and ϵ' is a protocol-dependent factor. The $\sqrt{N}$ scaling is a consequence of the law of large numbers, and it indicates that error mitigation can suppress the error by a larger factor in larger circuits. We propose the importance Clifford sampling as a key technique for error mitigation in large circuits to obtain this result.

97 MATHEMATICS AND COMPUTING↗

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗

Unveiling Hidden Shake-Up Features in the Uranyl M 4 -Edge Spectrum

The M 4,5 -edge high energy resolution X-ray absorption near-edge structure (HR-XANES) spectra of actinyls offer valuable insights into the electronic structure and bonding properties of heavy-element complexes. To conduct a comprehensive spectral analysis, it is essential to employ computational methods that accurately account for relativistic effects and electron correlation. In this work, we utilize variational relativistic multireference configurational interaction methods to compute and analyze the X-ray M 4 -edge absorption spectrum of uranyl. By employing these advanced computational techniques, we achieve excellent agreement between the calculated spectral features and experimental observations. Moreover, the calculations unveil significant shake-up features, which arise from the intricate interplay between strongly correlated 3d core-electron and ligand excitations. This research provides important theoretical insights into the spectral characteristics of heavy-element complexes. Furthermore, it establishes the foundation for utilizing M 4,5 -edge spectroscopy as a means to investigate the chemical activities of such complexes. By leveraging this technique, we can gain a deeper understanding of the bonding behavior and reactivity of heavy-element compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The First Neptunium Dioxide Neutron Diffraction Experiment at HFIR

Here we present the first neutron diffraction experiment of a 50 mg sample of neptunium(IV) dioxide (NpO 2 ) powder at the High Flux Isotope Reactor (HFIR) using the wide-angle neutron diffractometer (WAND 2 ) HB-2C. We also provide a brief literature survey to pinpoint gaps in the crystal, electronic, and magnetic structure of NpO 2 to determine the optimum experiment to be pursued in the next HFIR user proposal cycle. Characterizing actinide materials with neutron diffraction, inelastic scattering, neutron imaging, polarized neutron studies, and small-angle neutron scattering techniques will increase our understanding of 5f electron behavior and elucidate atomic and magnetic structural properties. This strategy will help fill gaps in the literature and provide theoretical physicists with the data needed to build physical and chemical models to guide and validate computational methods. These methods may be used for predicting structural and chemical features of actinide-based materials. This proof-of-principle effort begins to uncover how to advance actinide science and technology at Oak Ridge National Laboratory (ORNL) using neutron capabilities. This effort has already sparked collaborative discussions between researchers outside the Radioisotope Science and Technology Division, including theoretical physicists from the Nuclear Nonproliferation Division and beamline scientists in the Neutron Scattering Division. We intend to continue pursuing neutron experiments and develop a scientific strategy focused on studying transuranium elements (e.g., Np, Pu, Am, and Cm) using the advanced neutron capabilities at ORNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Matter: Computation Meets Experiments

The past few years have seen tremendous advances in our understanding of strongly interacting quantum systems. By combining progress in theoretical concepts and methods with algorithmic advances, computational methods have shed new light on key open questions in the physics of quantum matter, both for materials with strong electronic correlations and for interacting quantum gases. At the same time, new materials and new techniques have greatly increased the range of experimental information available. This Aspen Winter Conference will bring leading theorists with expertise in a broad range of computational methods together with experimentalists to discuss the potential of new methods, the accomplishments of existing methods and opportunities for future experiment/ theory collaboration. The meeting will focus on a broad set of physics questions of current interest in the field, for which computational methods and experiments have brought or have the potential to bring new insights.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Structural and Dynamic Disorder, Not Ionic Trapping, Controls Charge Transport in Highly Doped Conducting Polymers

Doped organic semiconductors are critical to emerging device applications, including thermoelectrics, bioelectronics, and neuromorphic computing devices. It is commonly assumed that low conductivities in these materials result primarily from charge trapping by the Coulomb potentials of the dopant counter-ions. Here, we present a combined experimental and theoretical study rebutting this belief. Using a newly developed doping technique, we find the conductivity of several classes of high-mobility conjugated polymers to be strongly correlated with paracrystalline disorder but poorly correlated with ionic size, suggesting that Coulomb traps do not limit transport. A general model for interacting electrons in highly doped polymers is proposed and carefully parameterized against atomistic calculations, enabling the calculation of electrical conductivity within the framework of transient localisation theory. Theoretical calculations are in excellent agreement with experimental data, providing insights into the disordered-limited nature of charge transport and suggesting new strategies to further improve conductivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Localization of Nonblinking Point Sources Using Higher-Order-Mode Detection and Optical Heterodyning: Developing a Strategy for Extending the Scope of Molecular, Super-resolution Imaging

While the stochastic, “blinking” nature of fluorescent systems has enabled the super-resolution of their localization by the fitting of their point-spread functions (PSFs), this strategy cannot be exploited for similar resolution of “nonblinking” systems, such as those that might be encountered in a coherent Raman experiment. An alternative method for subdiffraction-limited imaging lies in the exploitation of optical heterodyning. For example, if a Gaussian PSF (a TEM00 mode) of a point emitter is displaced with respect to the origin of the optical system, photons in the higher-order TEM modes carry information about that displacement. Information concerning the displacement can be extracted from photons in these higher-order modes. These photons can be collected by optical heterodyning, which exploits the large gain in a detector’s response to an optical signal from an emitter coupled to a local oscillator, which is prepared in the TEM of interest, e.g., TEM10. We have generalized and developed the heterodyning technique to localize point emitters via the detection of higher-order spatial modes. We have developed a theoretical approach to find a practical estimation limit of the localization parameters using a realistic model that accounts for shot noise, background noise, and Gaussian noise. To demonstrate the applicability of the method, we designed experiments in which a laser is a surrogate for one and two point emitters. Using the Fisher information and its accompanying Cramér-Rao lower bound, we demonstrate super-resolution localization in these cases: we show that objects can be localized to roughly 2–3 orders of magnitude of their point-spread function’s size for a given optical system. Finally and most importantly, it is suggested that the results will ultimately be generalizable to multiple emitters and, most importantly, to “nonblinking” molecular systems, which will be essential for broadening the scope of super-resolution measurements beyond the limits of fluorescence-based techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enabling Hyper-Differential Sensitivity Analysis for Ill-Posed Inverse Problems

Inverse problems constrained by partial differential equations (PDEs) play a critical role in model development and calibration. In many applications, there are multiple uncertain parameters in a model that must be estimated. However, high dimensionality of the parameters and computational complexity of the PDE solves make such problems challenging. A common approach is to reduce the dimension by fixing some parameters (which we will call auxiliary parameters) to a best estimate and use techniques from PDE-constrained optimization to estimate the other parameters. In this article, hyper-differential sensitivity analysis (HDSA) is used to assess the sensitivity of the solution of the PDE-constrained optimization problem to changes in the auxiliary parameters. Foundational assumptions for HDSA require satisfaction of the optimality conditions which are not always practically feasible as a result of ill-posedness in the inverse problem. Here we introduce novel theoretical and computational approaches to justify and enable HDSA for ill-posed inverse problems by projecting the sensitivities on likelihood informed subspaces and defining a posteriori updates. Our proposed framework is demonstrated on a nonlinear multiphysics inverse problem motivated by estimation of spatially heterogeneous material properties in the presence of spatially distributed parametric modeling uncertainties.

97 MATHEMATICS AND COMPUTING↗

Understanding and Estimating Error Propagation in Neural Networks for Scientific Data Analysis

Neural networks are increasingly integrated into scientific discovery, where input data reduction and model quantization play a key role in accelerating inference. However, understanding and mitigating the impact of these techniques on output error is critical for ensuring reliable results, particularly in tasks demanding high numerical precision. This paper introduces a comprehensive framework for optimizing neural network inference in scientific computing by combining data reduction and weight quantization while maintaining error-controlled outcomes. We develop theoretical analyses to bound error propagation under these reductions and propose a framework that balances computational performance with error constraints. Evaluation on real-world learning-based combustion simulations and satellite image classification demonstrates that our derived error bounds accurately predict observed errors while enabling significant computational speedup under our framework. This work highlights the potential for further leveraging advancements in modern lossy compression algorithms and hardware accelerators that support lower-precision formats.

He, Weiming [New Jersey Institute of Technology]↗

SNS: A Solution-Based Nonlinear Subspace Method for Time-Dependent Model Order Reduction

Several reduced order models have been successfully developed for nonlinear dynamical systems. To achieve a considerable speed-up, a hyper-reduction step is needed to reduce the computational complexity due to nonlinear terms. Many hyper-reduction techniques require the construction of nonlinear term basis, which introduces a computationally expensive offline phase. A novel way of constructing nonlinear term basis within the hyper-reduction process is introduced. In contrast to the traditional hyper-reduction techniques where the collection of nonlinear term snapshots is required, the SNS method avoids collecting the nonlinear term snapshots. Instead, it uses the solution snapshots that are used for building a solution basis, which enables avoiding an extra data compression of nonlinear term snapshots. As a result, the SNS method provides a more efficient offline strategy than the traditional model order reduction techniques, such as the DEIM, GNAT, and ST-GNAT methods. The SNS method is theoretically justified by the conforming subspace condition and the subspace inclusion relation. It is useful for model order reduction of large-scale nonlinear dynamical problems to reduce the offline cost. It is especially useful for ST-GNAT that has shown promising results, such as a good accuracy with a considerable online speed-up for hyperbolic problems in a recent paper by Choi and Carlberg [SIAM J. Sci. Comput., 41 (2019), pp. A26--A58], because ST-GNAT involves an expensive offline cost related to collecting nonlinear term snapshots. Error analysis for the SNS method is presented. Numerical results support that the accuracy of the solution from the SNS method is comparable to the traditional methods and a considerable speed-up (i.e., a factor of two to a hundred) is achieved in the offline phase.

97 MATHEMATICS AND COMPUTING↗

On the modeling of $t\bar{t}W^\pm$ signatures at the LHC

We discuss our recent progress on improving the theoretical description of the pp \to t\bar{t}W^\pm p p → t t ‾ W ± process with particular focus on fiducial signatures that include top-quark decays. We employ a wide range of techniques from parton shower matching, narrow width approximation to full off-shell computations to assess the importance of the different modeling approaches.

Kraus, Manfred↗