Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computing methods”

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 1,477 records · Page 82

GenAI4UQ: A software for forward and inverse uncertainty quantification using conditional generative AI

We introduce GenAI4UQ, a software package for forward and inverse uncertainty quantification in model calibration, parameter estimation, and ensemble forecasting. GenAI4UQ leverages a generative AI-based conditional modeling framework to address limitations of traditional inverse modeling techniques, such as Markov Chain Monte Carlo (MCMC) methods. By replacing computationally intensive iterative processes with a direct, learned mapping, GenAI4UQ enables efficient calibration of input parameters and generation of predictions directly from observations. The software supports rapid ensemble forecasting with robust uncertainty quantification while maintaining computational and storage efficiency. Built-in auto-tuning of hyperparameters simplifies model training, ensuring accessibility for users with varying expertise. Its versatile conditional generative framework is applicable across diverse scientific domains. While GenAI4UQ offers significant advantages in flexibility and efficiency, users should interpret its uncertainty estimates with caution in data-sparse scenarios, as the model may overestimate uncertainty—an effect common to all surrogate-based approaches including MCMC with surrogate models. Despite this, GenAI4UQ transforms inverse modeling by providing a fast, reliable, and user-friendly solution. It empowers researchers and practitioners to quickly estimate parameter distributions and generate model predictions for new observations, facilitating efficient decision-making and advancing the state of uncertainty quantification in computational modeling.

97 MATHEMATICS AND COMPUTING↗

Atomic-scale Imaging of PGM-free Catalyst Active Sites by 30 keV 4D-STEM

Platinum group metal (PGM)-free catalysts have attracted a large amount of attention due to their potential for enabling low-cost, commercially viable hydrogen fuel cells and electrolyzers [1]. Degradation of cell performance remains a significant challenge, however, currently limiting implementation of devices that utilize PGM-free materials [2]. Controlling degradation in these materials requires a better understanding of the atomic structure of active sites, thought to be FeN4 structures within a graphitic carbon lattice, which would enable more accurate prediction of the associated degradation pathways. The exact structural arrangement of active sites in these materials is still debated since a range of structures have been computationally predicted and methods to directly validate these models are still needed [1,3-5].While scanning transmission electron microscopy (STEM) has provided initial glimpses into the nature of the proposed active sites, detailed evaluation of the structure of these sites is challenging since they involve defects, edges, and nitrogen dopants in the graphitic lattice, which are susceptible to beam damage [6,7]. Conventional STEM imaging and spectroscopy methods exacerbate this problem with dose-inefficiency or nonideal contrast characteristics for imaging light and heavy elements simultaneously. A method that minimizes damage while generating dose-efficient, easily-interpretable contrast for both light and heavy elements is therefore needed to facilitate imaging of sensitive active site atomic structure.Here, we demonstrate direct atomic-scale imaging of PGM-free catalyst active sites by low-voltage four-dimensional (4D)-STEM. To accomplish this, we pair a 30 keV aberration-corrected probe, which minimizes knock-on damage, with a fast pixelated detector that has optimal performance at low beam energies [8]. We show how this setup enables relatively simple center-of-mass (CoM) techniques [9,10] to produce images with an increased signal-to-noise ratio (SNR) and light-element contrast over conventional imaging modes, allowing the entire active site structure to be imaged at the atomic scale. Moreover, we show how electron ptychography [11,12] enables images with further improved characteristics to be obtained, for example by minimizing residual aberrations, which provides a more accurate structural representation. The increased understanding that these low-voltage 4D-STEM techniques will provide about the atomic structure of PGM-free active sites and their associated degradation pathways will facilitate rational design of next-generation materials, promoting development of low-cost hydrogen fuel cells, electrolyzers, and other energy conversion devices [13]. References:[1] U Martinez et al., Adv. Mater. 31 (2019), p. 1806545.[2] Y Shao et al., Adv. Mater. 31 (2019), p. 1807615.[3] J Kneebone et al., J. Phys. Chem. C 121, (2017), p.16283.[4] T Mineva et al., ACS Catal. 9 (2019), p. 9359.[5] J Li et al., Nat. Catal. 4 (2021), p. 10.[6] T Susi et al., ACS Nano 6 (2012), p. 8837.[7] P Zelenay and DJ Myers, US Department of Energy Hydrogen and Fuel Cells Program 2017 Annual Merit Review and Peer Evaluation Meeting, Washington, DC (2017).[8] H Ryll et al., J. Inst. 11 (2016), p. P04006.[9] K Müller et al., Nat. Commun. 5 (2014), p. 5653.[10] I Lazic et al., Ultramicroscopy 160 (2016), p. 265.[11] H Yang et al., Nat. Commun. 7 (2016), p. 12532.[12] Y Jiang et al., Nature 559 (2018), p. 343. [13] Research sponsored by the Hydrogen and Fuel Cell Technologies Office, Office of Energy Efficiency and Renewable Energy, US Department of Energy (DOE). Research was conducted at the Center for Nanophase Materials Sciences, which is a DOE Office of Science User Facility.

Zachman, Michael↗

Predictive Synthesis

Current solid state synthesis intrinsically involves a multidimensional space which is challenging to parametrize and predict. The diversity of extended structures comes from the diversity of basic properties of elements of the Periodic system which may exhibit a variety of bonding modes. The fundamental challenges of prediction of the preparative outcome are further complicated by practical synthetic issues. Current development in computational and experimental methods calls for collaborative efforts to make solid state synthesis more predictable. Overall, this Perspective discusses several steps in this direction, including integration of predictions of synthetic conditions with new structure predictions, widespread in situ studies to obtain a panoramic view of the reaction mechanism, and the creation of a synthetic database to properly document all synthetic efforts, including the unsuccessful ones.

36 MATERIALS SCIENCE↗

Zero-Field Splitting Calculations by Multiconfiguration Pair-Density Functional Theory

Zero-field splitting (ZFS) is a fundamental molecular property that is especially relevant for single-molecule magnets (SMMs), electron paramagnetic resonance spectra, and quantum computing. Developing a method that can accurately predict ZFS parameters can be very powerful for designing new SMMs. One of the challenges is to include external correlation in an inherently multiconfigurational open-shell species for the accurate prediction of magnetic properties. Previously available methods depend on expensive multireference perturbation theory calculations to include external correlation. In this paper, we present spin-orbit-inclusive multiconfiguration and multistate pair-density functional theory (MC-PDFT) calculations of ZFSs; these calculations have a cost comparable to complete-active-space self-consistent field (CASSCF) theory, but they include correlation external to the active space. We found that combining a multistate formulation of MC-PDFT, namely, compressed-state multistate pair-density functional theory, with orbitals optimized by weighted-state-averaged CASSCF, yields reasonably accurate ZFS results

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hetero-site Double Core Ionization Energies with Sub-electronvolt Accuracy from Delta-Coupled-Cluster Calculations

Benchmark scalar-relativistic delta-coupled-cluster calculations of hetero-site double core ionization energies of small molecules containing second-row elements are reported. The present study has focused on the high-spin triplet components of two-site double core-ionized states, which are single reference in character and consistent with the use of standard coupledcluster methods. Contributions to computed double core ionization energies from electron-correlation and basis-set effects as well as corrections to the core-valence separation approximation have been analyzed. Furthermore, on the basis of systematic convergence of computational results with respect to these effects, delta-coupled-cluster calculations have been shown to be capable of providing accurate double core ionization energies with remaining errors estimated to be below 0.3 eV, and thus they are recommended for use to facilitate experimental studies of two-site double core-ionized states that are involved in X-ray pump/X-ray probe studies of electronic and molecular dynamics following inner shell ionization or excitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanism of O 2 Activation and Cysteine Oxidation by the Unusual Mononuclear Cu(I) Active Site of the Formylglycine-Generating Enzyme

The formylglycine-generating enzyme (FGE) catalyzes the selective oxidation of a peptidyl-cysteine to form formylglycine, a critical cotranslational modification for type I sulfatase activation and a useful bioconjugation handle. Previous studies have shown that the substrate peptidyl-cysteine binds to the linear bis-thiolate Cu(I) site of FGE to form a trigonal planar tris-thiolate Cu(I) structure that activates O 2 for the oxidation of the C β –H of the cysteine substrate via an unknown mechanism. Here, we employed a combination of stopped-flow kinetic, spectroscopic (UV–vis absorption, XAS, and EPR), and computational (DFT/TD-DFT calculations) methods to observe and characterize the key intermediates in this reaction for FGE from Streptomyces coelicolor. Our results define the reaction coordinate of FGE, which involves H-atom abstraction from the C β –H bond of the cysteine substrate by a reactive Cu(II)–O 2 •– species to form the now experimentally observed Cu(I)–OOH intermediate bound to a peptidyl-thioaldehyde, which proceeds to oxidize one of the protein-derived cysteine residues to a sulfenate that is end-on O-coordinated to Cu(I). These results provide fundamental insights into how the unusual mononuclear Cu(I) site of FGE activates O 2 for cysteine oxidation and stores oxidizing equivalents during catalysis by employing a Cu(I)–sulfenate intermediate with an end-on O-coordination that is unprecedented in biology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing Surface Plasmon Dynamics in Periodic Nanostructures through Ultrafast Electron Microscopy

Surface plasmon polaritons (SPPs) can be manipulated to localize and guide light in subwavelength distances, enabling them to find applications in a wide range of areas, from sensing to quantum computing. Among several methods of SPP excitation, periodic arrays of nano- and microstructures are of particular interest, as they enable engineering SPP properties through structural parameters. Here, in this study, using the photon-induced near-field electron microscopy (PINEM) technique, we investigated the mode formation, coupling, interference, and decay of SPPs in square and hexagonal arrays of circular nanoholes under both visible and near-infrared excitation. Polarization-resolved analysis revealed the key factors governing SPP localization and interference patterns, showing that the periodicity and symmetry of the array primarily determine the SPP interference patterns and their orientation, while pump polarization mainly modulates their intensity. Time-resolved PINEM measurements demonstrated the spatial dependence of the SPP temporal characteristics. In addition, cathodoluminescence (CL) spectroscopy was employed to examine the intrinsic plasmonic characteristics of the structure. Finite difference time domain (FDTD) simulations showed strong agreement with both PINEM and CL measurements on the spatial and spectral behavior of SPPs. Understanding the spatiotemporal dynamics of SPPs on nanostructures beyond the diffraction limit is crucial for optimizing plasmonic structures for advanced photonic and quantum technologies.

Plasmonics↗

NLML: A Deep Neural Network Emulator for the Exact Nonlinear Interactions in a Wind Wave Model

Nonlinear wave interactions describe the resonant energy transfer between wave components, playing a fundamental role in the evolution of ocean wave spectra. Nonlinear wave interactions significantly influence wave growth and development, making them essential for accurate wave modeling. However, resolving the full six-dimensional Boltzmann integral of the exact nonlinear wave interactions (Webb-Resio-Tracy method, WRT) is computationally expensive, limiting its application in real-time operational wave forecasting and for research purposes. Current approximations, such as the Discrete Interaction Approximation (DIA), prioritize computational speed over accuracy, resulting in significant errors in wave mean parameters. Here, we introduce NLML, a machine learning (ML) emulator designed to approximate the exact nonlinear wave interactions within WAVEWATCH III (WW3), with the goal of achieving the accuracy of WRT while maintaining the stability and computational speed of DIA. By leveraging GPU capabilities such as half precision inference, we achieved substantial speedups, up to 136x mathematical equation faster than the WRT and only a modest 1.04x mathematical equation slowdown relative to DIA, while achieving 2x mathematical equation the accuracy of DIA in global wave spectral energy and mean wave parameters, with up to 7x mathematical equation higher accuracy in some regions. Unlike previous ML approaches, NLML maintained inherent stability throughout model integration in a standalone, year-long WW3 simulation, without requiring additional constraints. Our new ML parameterization bridges the gap between accuracy and efficiency, offering a promising alternative for improving wave modeling in operational settings and research purposes.

16 TIDAL AND WAVE POWER↗

Enabling reactive microscopy with MicroMator

Microscopy image analysis has recently made enormous progress both in terms of accuracy and speed thanks to machine learning methods and improved computational resources. This greatly facilitates the online adaptation of microscopy experimental plans using real-time information of the observed systems and their environments. Applications in which reactiveness is needed are multifarious. Here we report MicroMator, an open and flexible software for defining and driving reactive microscopy experiments. It provides a Python software environment and an extensible set of modules that greatly facilitate the definition of events with triggers and effects interacting with the experiment. We provide a pedagogic example performing dynamic adaptation of fluorescence illumination on bacteria, and demonstrate MicroMator’s potential via two challenging case studies in yeast to single-cell control and single-cell recombination, both requiring real-time tracking and light targeting at the single-cell level.

47 OTHER INSTRUMENTATION↗

Machine learning for deep elastic strain engineering of semiconductor electronic band structure and effective mass

Abstract The controlled introduction of elastic strains is an appealing strategy for modulating the physical properties of semiconductor materials. With the recent discovery of large elastic deformation in nanoscale specimens as diverse as silicon and diamond, employing this strategy to improve device performance necessitates first-principles computations of the fundamental electronic band structure and target figures-of-merit, through the design of an optimal straining pathway. Such simulations, however, call for approaches that combine deep learning algorithms and physics of deformation with band structure calculations to custom-design electronic and optical properties. Motivated by this challenge, we present here details of a machine learning framework involving convolutional neural networks to represent the topology and curvature of band structures in k -space. These calculations enable us to identify ways in which the physical properties can be altered through “deep” elastic strain engineering up to a large fraction of the ideal strain. Algorithms capable of active learning and informed by the underlying physics were presented here for predicting the bandgap and the band structure. By training a surrogate model with ab initio computational data, our method can identify the most efficient strain energy pathway to realize physical property changes. The power of this method is further demonstrated with results from the prediction of strain states that influence the effective electron mass. We illustrate the applications of the method with specific results for diamonds, although the general deep learning technique presented here is potentially useful for optimizing the physical properties of a wide variety of semiconductor materials.

Tsymbalov, Evgenii↗

Structural evolution of two-dimensional silicates using a “bond-switching” algorithm

Silicates are the most abundant materials in the earth's crust. In recent years, two-dimensional (2D) versions of them grown on metal supports (known as bilayer silicates) have allowed their study in detail down to the atomic scale. These structures are self-containing. Furthermore, they are not covalently bound to the metal support but interact with it through van der Waals forces. Like their three-dimensional counterparts, the 2D-silicates can form both crystalline and vitreous structures. Furthermore, the interconversion between vitreous to crystalline structures has been experimentally observed at the nanoscale. While theoretical work has been carried out to try to understand these transformations, a limitation for ab initio methods, and even molecular dynamics methods, is the computational cost of studying large systems and long timescales. In this work, we present a simple and computationally inexpensive approach, that can be used to represent the evolution of bilayer silicates using a bond-switching algorithm. This approach allows reaching equilibrium ring size distributions as a function of a parameter that can be related to the ratio between temperature and the energy required for the bond-switching event. The ring size distributions are compared to experimental data available in the literature.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

An Entropy-Maximization Approach to Automated Training Set Generation for Interatomic Potentials

Machine learning-based interatomic potentials are currently garnering a lot of attention as they strive to achieve the accuracy of electronic structure methods at the computational cost of empirical potentials. Given their generic functional forms, the transferability of these potentials is highly dependent on the quality of the training set, the generation of which can be highly labor-intensive. Good training sets should at once contain a very diverse set of configurations while avoiding redundancies that incur cost without providing benefits. We formalize these requirements in a local entropy-maximization framework and propose an automated sampling scheme to sample from this objective function. We show that this approach generates much more diverse training sets than unbiased sampling and is competitive with hand-crafted training sets.

74 ATOMIC AND MOLECULAR PHYSICS↗

Synergistic Effect of Diatomic Mo-B Site Confined in Graphene-Like C 2 N Enables Electrocatalytic Nitrogen Reduction via Novel Mechanism

Structural modulation of active site with atomic-level precision is of great importance to meet the activity and selectivity challenges that electrocatalysts are commonly facing. In this work, we have designed a metal (M) - nonmetal (NM) diatomic site embedded in graphene-like C 2 N (denoted as Mo-B@C 2 N), where electrocatalytic N 2 reduction reaction (eNRR) was thoroughly explored using density functional theory combined with computational hydrogen electrode method. Compared to M-M diatomic sites, Mo-B site can generate pronounced synergistic effect that led to eNRR proceeding via a novel quasi-dissociative reaction mechanism that has not been reported relative to the conventional enzymatic, consecutive, distal, and alternating associative mechanism. This newly uncovered mechanism in which N-N bond scission takes place immediately after the first proton-coupled electron transfer (PCET) step (i.e., *NH-*N + H + + e – ® *NH 2 *N) has demonstrated much advantage in PCET process over the four conventional mechanism in terms of thermodynamic barrier, except that the adsorption of side-on *N 2 seemed thermodynamically unfavorable (DG ads = 0.61eV). Our results have revealed that the activation of inert NºN triple bond is dominated by π*-backdonation mechanism as a consequence of charge transfers from both B and Mo site, and unexpectedly, from the substrate C 2 N itself as well. Moreover, the hybrid Mo-B diatomic site demonstrated superior performance over either Mo-Mo or B-B site for driving eNRR. Furthermore, our study could provide insight into the delicate relationships among atomic site, substrate and electrocatalytic performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strong-Weak Duality via Jordan-Wigner Transformation: Using Fermionic Methods for Strongly Correlated su(2) Spin Systems

The Jordan-Wigner transformation establishes a duality between $su(2)$ and fermionic algebras. We present qualitative arguments and numerical evidence that when mapping spins to fermions, the transformation makes strong correlation weaker, as demonstrated by the Hartree-Fock approximation to the transformed Hamiltonian. This result can be rationalized in terms of rank reduction of spin shift terms when transformed to fermions. Conversely, the mapping of fermions to qubits makes strong correlation stronger, complicating its solution when one uses qubit-based correlators. The presence of string operators poses challenges to the implementation of quantum chemistry methods on classical computers, but these can be dealt with using established techniques of low computational cost. Here, our proof of principle results for XXZ and J1-J2 Heisenberg (in 1D and 2D) indicate that the JW transformed fermionic Hamiltonian has reduced complexity in key regions of their phase diagrams, and provides a better starting point for addressing challenging spin problems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Optical and transport properties of plasma mixtures from ab initio molecular dynamics

Predicting the charged particle transport properties of warm dense matter/hot dense plasma mixtures is a challenge for analytical models. High accuracy ab initio methods are more computationally expensive, but can provide critical insight by explicitly simulating mixtures. In this work, we investigate the transport properties and optical response of warm dense carbon–hydrogen mixtures at varying concentrations under either conserved electronic pressure or mass density at a constant temperature. We compare options for mixing the calculated pure species properties to estimate the results of the mixtures. We find that a combination of the Drude model with the Matthiessen's rule works well for DC electron transport and low-frequency optical response. This breaks down at higher frequencies, where a volumetric mix of pure-species AC conductivities works better.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tutorial: Defects in topological semimetals

Three-dimensional topological semimetals are a class of electronic materials in which their bulk and surface states contain linear band touching nodes near the Fermi level. Like semiconductors, their properties will be affected by point and extended defects in their crystal structures, although the extent to which defects and disorders influence topological semimetals may differ in key ways due to their unique electronic structures. In this Tutorial, we provide an overview of the defects in topological semimetals, covering both computational and experimental methods for exploring defect-property relationships. We also include a discussion on open questions that still need to be explored further.

36 MATERIALS SCIENCE↗

Accurate noncovalent interactions in atomistic systems via quantum Drude oscillators

Accurately modeling polarization and van der Waals (vdW) interactions in atomistic systems typically requires high-level quantum-mechanical methods that are computationally expensive, hence limited in applicability. To address this challenge, efficient yet physically grounded models are needed—ones that not only enable accurate predictions but also provide insight into how noncovalent interactions scale in complex molecular and material systems. This review highlights the quantum Drude oscillator (QDO) model, a physically motivated and computationally efficient framework that captures the essential features of electronic response, including polarization and dispersion forces, across a wide range of chemical and material systems. We discuss how the QDO model quantitatively reproduces the polarization response of many-electron atoms and how key components of noncovalent interactions—exchange-repulsion, polarization, and dispersion—emerge naturally in QDO dimers. Furthermore, the model provides predictive scaling laws that elucidate trends in polarizability and dispersion across the periodic table and in molecular assemblies. By uniting interpretability, accuracy, and efficiency, the QDO model offers a versatile approach for modeling noncovalent interactions in systems ranging from isolated molecules to complex condensed phases and nanostructured materials.

Khabibrakhmanov, Almaz [Univ. of Luxembourg, Luxem↗

Geometry dependence of the fluctuation intensity in gyrokinetic turbulence

The findings of an investigation into the properties of the three dimensional (3D) saturated fluctuation intensity of the electric potential in gyrokinetic turbulence simulations is presented. Scans in flux surface elongation and Shafranov shift are used to isolate the tokamak geometric dependencies. The potential intensity required in order to compute exact fluxes by a quasilinear method is determined using linear eigenmodes computed with the gyrokinetic code. A model of this non-linear intensity is constructed using the linear eigenmode properties and the geometry shape functions obtained from the 3D intensity spectrum. The model computes the poloidal wavenumber spectrum of the electron and ion energy fluxes with unprecedented accuracy. In conclusion, new insights are gained into the way zonal flow mixing saturates ion-scale turbulence by controlling the radial wavenumber width of the turbulence spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗