SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Efficient Simulation of Open Quantum Systems on NISQ Trapped‐Ion Hardware

Abstract Simulating open quantum systems, which interact with external environments, presents significant challenges on noisy intermediate‐scale quantum (NISQ) devices due to limited qubit resources and noise. In this study, an efficient framework is proposed for simulating open quantum systems on NISQ hardware by leveraging a time‐perturbative Kraus operator representation of the system's dynamics. This approach avoids the computationally expensive Trotterization method and exploits the Lindblad master equation to represent time evolution in a compact form, particularly for systems satisfying specific commutation relations. The efficiency of this method is demonstrated by simulating quantum channels, such as the continuous‐time Pauli channel and damped harmonic oscillators, on NISQ trapped‐ion hardware, including IonQ Harmony and Quantinuum H1‐1. Additionally, hardware‐agnostic error mitigation techniques are introduced, including Pauli channel fitting and quantum depolarizing channel inversion, to enhance the fidelity of quantum simulations. These results show strong agreement between the simulations on real quantum hardware and exact solutions, highlighting the potential of Kraus‐based methods for scalable and accurate simulation of open quantum systems on NISQ devices. This framework opens pathways for simulating more complex systems under realistic conditions in the near term.

Burdine, Colin [Department of Electrical and Compu↗

DOE OSTI2025

Acceleration of the particle-in-cell code Osiris with graphics processing units

Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ~14 times faster and ~7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Nanopolysaccharide Builder: A User-Friendly Tool for Atomistic Models of Polysaccharide-Based Nanostructures

Here, we introduce Nanopolysaccharide Builder (NPB), a user-friendly software tool designed to construct polysaccharide nanostructures─mainly those based on cellulose, chitin, and chitosan─using experimental data or user-defined parameters. NPB enables the generation of cellulose and chitin allomorphs with customizable biochemical topologies and also facilitates the construction of large bundles that replicate nanostructures found in biological support systems, including plant cell walls and arthropod cuticles. The software outputs atomic Cartesian coordinates in Protein Data Bank (PDB) format and also provides atom connectivity files in PSF and PARM formats, ensuring seamless integration with major molecular dynamics (MD) engines such as NAMD, CHARMM, GROMACS, AMBER, OpenMM, and LAMMPS. Built on an interactive visualization framework, NPB features a graphical user interface (GUI) and supports both macOS and Linux operating systems. By enabling detailed atomic-scale studies of polysaccharide evolution in extracellular matrices and cell walls of algae, bacteria, fungi, and plants, NPB is poised to advance AI-guided research in sustainable chemical development and biomass utilization.

Wan, Zhangmin [Univ. of British Columbia, Vancouve↗

DOE OSTI2021

Suppressing Auger Recombination in Multiply Excited Colloidal Silicon Nanocrystals with Ligand-Induced Hole Traps

Nonradiative Auger recombination in multiply excited nanocrystals is a dominant efficiency loss pathway for nanocrystal-containing optoelectronic devices that rely on high-rate emission and absorption operating conditions. Overcoming Auger recombination in these quantum-confined systems is therefore a longstanding challenge to the synthetic nanocrystal as well as device manufacturing communities. Several successful strategies have been realized to reduce Auger recombination, but they rely on complex and time-consuming nanocrystal core/shell synthesis. Alternatively, controlling Auger rates by varying the nanocrystal–ligand-binding chemistry is a promising route to obtain functional tunability, which reduces the barrier to large-scale manufacturing. The covalent surface chemistry and the extremely long-lived photoexcited lifetimes of silicon nanocrystals (Si NCs) make them a unique system among colloidal semiconductor NCs to study the intersection between surface chemistry and photoexcited carrier dynamics. Here, we show that changing the functional group that binds a saturated dodecyl ligand to the surface of nonthermal plasma-synthesized Si NCs from alkyl to thiolate slows Auger recombination rates within multiply excited Si NCs. This reduction in Auger rate persists across Si NC sizes ranging from 3.5 to 8 nm in diameter, but the expected linear dependence of Auger rates on the NC volume is retained for both alkyl and alkylthiolate surface terminations. To understand the origin behind this elongation, we carry out steady-state and time-resolved photoluminescence measurements as well as time-resolved terahertz spectroscopy measurements. These measurements reveal that thiolate groups introduce mid-gap surface states, which, we argue, reduces the photoexcited electron–hole overlap and elongates Auger recombination times. These results highlight how a typically detrimental chemical species—mid-band gap NC surface states—can be beneficial under high-rate absorption/emission conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Scale Invariance of Hot Spot Formation in TATB High Explosives

Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Water-mediated ion transport in an anion exchange membrane

Water is a critical component in polyelectrolyte anion exchange membranes (AEMs). It plays a central role in ion transport in electrochemical systems. Gaining a better understanding of molecular transport and conductivity in AEMs has been challenged by the lack of a general methodology capable of capturing and connecting water dynamics, water structure, and ionic transport over time and length scales ranging from those associated with individual bond vibrations and molecular reorientations to those pertaining to macroscopic AEM performance. In this work, we use two-dimensional infrared spectroscopy and semiclassical simulations to examine how water molecules are arranged into successive solvation shells, and we explain how that structure influences the dynamics of bromide ion transport processes in polynorbornene-based materials. We find that the transition to the faster transport mechanism occurs when the reorientation of water molecules in the second solvation shell is fast, allowing a robust hydrogen bond network to form. Our findings provide molecular-level insights into AEMs with inherent transport of halide ions, and help pave the way towards a comprehensive understanding of hydroxide ion transport in AEMs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

A zeroth-order active-space frozen-orbital embedding scheme for multireference calculations

Multireference computations of large-scale chemical systems are typically limited by the computational cost of quantum chemistry methods. In this work, we develop a zeroth-order active space embedding theory [ASET(0)], a simple and automatic approach for embedding any multireference dynamical correlation method based on a frozen-orbital treatment of the environment. ASET(0) is combined with the second-order multireference driven similarity renormalization group and tested on several benchmark problems, including the excitation energy of 1-octene and bond-breaking in ethane and pentyldiazene. Finally, we apply ASET(0) to study the singlet–triplet gap of p-benzyne and 9,10-anthracyne diradicals adsorbed on a NaCl surface. Our results show that despite its simplicity, ASET(0) is a powerful and sufficiently accurate embedding scheme applicable when the coupling between the fragment and the environment is in the weak to medium regime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Dynamical Large Deviations of Two-Dimensional Kinetically Constrained Models Using a Neural-Network State Ansatz

We use a neural-network ansatz originally designed for the variational optimization of quantum systems to study dynamical large deviations in classical ones. We use recurrent neural networks to describe the large deviations of the dynamical activity of model glasses, kinetically constrained models in two dimensions. We present the first finite size-scaling analysis of the large-deviation functions of the two-dimensional Fredrickson-Andersen model, and explore the spatial structure of the high-activity sector of the South-or-East model. These results provide a new route to the study of dynamical large-deviation functions, and highlight the broad applicability of the neural-network state ansatz across domains in physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Improving Schrödinger Equation Implementations with Gray Code for Adiabatic Quantum Computers

We reformulate the continuous-space Schrödinger equation in terms of spin Hamiltonians. For the kinetic energy operator, the critical concept facilitating the reduction in model complexity is the idea of position encoding. A binary encoding of position produces a spin-1/2 Heisenberg-like model and yields exponential improvement in space complexity when compared to classical computing. Encoding with a binary reflected Gray code (BRGC), and a Hamming-distance-2 Gray code (H2GC) reduces the model complexity down to the 𝑋⁢𝑍 and transverse Ising model, respectively. For 𝐴 qubits BRGC yields 2 𝐴 positions and is reduced to its 2-local form with O⁡(𝐴) ancillary qubits. H2GC yields 2 𝐴/2+1 positions with O⁡(𝐴 2 ) three-local penalty terms. We also identify the bijective mapping between diagonal unitaries and the Walsh series, producing the mapping of any real potential to a series of 𝑘 -local Ising models through the fast Walsh transform. Finally, in a finite volume, we provide some numerical evidence to support the claim that the total time needed for adiabatic evolution is protected by the infrared cutoff of the system. As a result, initial state preparation from a free-field wave function to an interacting system is expected to exhibit polynomial time complexity with volume and constant scaling with respect to lattice discretization for all encodings. For H2GC, if the evolution starts with the transverse Hamiltonian due to hardware restrictions, then penalties are dynamically introduced such that the low-lying spectrum reproduces the energy levels of the Laplacian. The adiabatic evolution of the penalty Hamiltonian is therefore sensitive to the ultraviolet scale. It is expected to exhibit polynomial time complexity with lattice discretization, or exponential time complexity with respect to the number of qubits given a fixed volume.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

catalysis (heterogeneous)↗

DOE OSTI2022

Experimental Investigation of Polymer-Induced Fouling of Heater Tubes in the First-Ever Polymer Flood on Alaska North Slope Part II

A polymer flood pilot has been ongoing in the Schrader Bluff viscous oil reservoir at Milne Point on the Alaska North Slope. The results from the pilot are encouraging. However, a major concern of the operator is the influence of polymer on the production system after breakthrough, especially the fouling in heat exchangers. This work applies a multiexperimental approach to study the severity of polymer-induced fouling in both dynamic and static states of produced fluids to determine safe operating conditions. Dynamic scale loop (DSL) tests were conducted to study fouling due to polymer at different skin temperatures (165–350°F) in a dynamic state of fluid flow where the fluids’ flow mimics the residence time of fluids in the heat exchanger of the field pilot. Static deposit tests were also conducted at similar skin temperatures of 165–250°F using a novel experimental apparatus designed and built in-house. It was found that at higher skin temperatures of 250–350°F, tube blocking was observed in the DSL tests, whereas the tests at 165–200°F did not show any tube blocking, even in a more extended test period. The deposit test showed that the deposit rate generally increases with skin temperature, and the presence of polymer aggravates the fouling. The copper tube performs best when the skin temperatures are 165–200°F, while the stainless steel tube performs best at a skin temperature of 250°F. These experiments also manifested the influence of the cloudpoint of the solution as the deposit rate increased significantly when the skin temperature was higher than the solution cloudpoint. The study provides a source of practical guidance to the field operations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Use of Remote Sensing and In-Situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the ACME Model

The overachieving goal of the whole CMDV-MCS project is to improve understanding of warm season continental convection and to develop treatments of convection and microphysics capable of representing mesoscale convective systems (MCSs) features in large-scale models. Our tasks for this project contributing to the overachieving goal include: (1) Improve the ice nucleation formulation for MG2 and P3 cloud microphysics schemes; (2) Improve the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM; and (3) Test the performance of improved ice microphysics in E3SM with observation data. In this project, we have (1) Improved the ice nucleation parameterization for MG2 and P3 in E3SM by implementing two advanced empirical parameterizations with connection to aerosols. The two deterministic heterogeneous ice nucleation parameterizations (i.e., DeMott et al., 2015; Niemand et al., 2012) were merged with the MG2 and P3 cloud microphysics schemes in E3SM. Long-term simulations were conducted to examine the impacts of these new parameterizations on simulated cloud properties; (2) Improved the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM. We evaluated the double Gaussian PDF of vertical velocity simulated by the Cloud Layers Unified By Binormals (CLUBB) and the sub-column vertical velocity sampled from the Subgrid Importance Latin Hypercube Sampler (SILHS) in E3SM. We introduced the vertical velocity variance induced by topographic gravity waves for ice nucleation and droplet activation; and (3) Tested the performance of improved ice microphysics in E3SM with observation data. We tested the new treatments of ice nucleation in the single column model (SCM) mode for the stratiform mixed-phase clouds observed during 9-10 October 2004 in the DOE ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) and for the convective clouds observed on 20 May 2011 in the Midlatitude Continental Convective Clouds Experiment (MC3E). Modeled ice nucleating particles (INPs) concentrations were compared against observations collected around the globe.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

A Scalable Hardware-and-Human-in-the-Loop Grid-interactive Efficient Building Equipment Performance Dataset

This project developed a publicly available, high-fidelity dataset about the interactions among humans, homes, and heat pumps supporting grid interactive efficient buildings to balance demand on the grid with comfort for occupants. Laboratory measurements and simulations of the hardware capture the second-scale electric power dynamics of heat pumps providing grid services like load shifting and load shedding. Field measurements, behavior tracking, and qualitative surveys of people in their homes over multiple years—including experimentally adjusting the heating and cooling system to provide grid services—to capture the reciprocal effect of human behavior on grid services, and grid services on human comfort. Taken together, these data capture the complete Hardware and Human in the loop system for residential heat pumps, reducing large uncertainties in simulation for design, and models for control of heat pumps, and grid-interactive buildings.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

2.1.3.404 - WEC Array Power Management and Output Simulation Tool

An array of wave energy converter (WEC) devices has variations in power output due to the chaotic nature of the waves. Eliminating or mitigating the power fluctuations is important for reducing the integration impacts of WEC plants in both distribution and transmission grids, and in standalone isolated power systems. Reduced variability of WEC-generated power in combination with energy storage or power management control at each WEC and at the array level will help increasing hosting capacity of distribution feeders for this type of variable renewable generation, and minimization of electric losses. To reduce risks and risk perception, gain key stakeholder acceptance, and enable developers to design effective and compatible energy plants with arrays of WEC devices, it requires the use of modeling tools to simulate the resource environment, device dynamics, utility power system response, as well as the development of an interface for an array controller. The project will create a publicly accessible numerical modeling framework to empower the wave energy sector to design projects of various scales (kW-100s MW), which are optimized on a plant performance basis and are compatible with different power systems and wave conditions. The framework will integrate with WEC-Sim, a wave environment model (SWAN-FUNWAVE), as well as established relevant electrical analysis tools, to model the grid system and interconnection, and optimize power output and power management for the WEC array.

POWER TRANSMISSION AND DISTRIBUTION,TIDAL AND WAVE↗

DOE OSTI2023

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides (Final Report)

The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides. Final Report

36 MATERIALS SCIENCE↗

DOE OSTI2022

Large-Area Lasing in Nanoscale Complex Media: The Critical Role of Local Dielectric Environment

Controlling how electromagnetic waves interact with complex media is critical for applications in imaging and focusing. Such lightwave interactions with complex media can lead to dramatic optical effects like lasing. While much work in random lasing focus on understanding how gain and scattering co-operatively generate lasing, little work has focused on how to manipulate the lasing threshold without modifying the structural disorder. Here, a simple, mostly unexplored, strategy is demonstrated that employs atomic layer deposition (ALD) to tune the local near-field environment while preserving the underpinning disorder—controlling lasing in a nanoscale complex medium on a large scale (>cm 2 ). The nanoscale complex medium is a quasi-2D system of coupled zinc oxide nanospheres with overall thickness deep in the sub-wavelength regime (≈λ/4). Near-ultraviolet femtosecond spectroscopy probes the broadband response of the gain nanomaterial, details how ALD process fundamentally modifies the fast-picosecond and slow-nanosecond carrier dynamics, and informs on the relevant timescales critical for lasing. Full-field electromagnetic simulations provide critical insights about how near-field dielectric environment modifies the nanostructure's scattering cross-section, which ultimately results in enhanced lasing. Importantly, these results highlight a simple path to control how electromagnetic waves interact in a complex medium, a key step toward large-scale implementation of complex lasers.

36 MATERIALS SCIENCE↗

DOE OSTI2023

TOGETs Overview

The Transmission Optimization with Grid Enhancing Technologies (TOGETs) project will guide future research to fill in gaps on publicly available knowledge of Grid Enhancing Technologies (GETs) and conduct a full scale, multi-faceted field exercise. Ultimately the project aims to test, verify and validate deployed GETs in a transmission system.

24 POWER TRANSMISSION AND DISTRIBUTION↗