Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum dynamics simulation”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects

As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.

force field↗
DOE OSTI2024

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗
DOE OSTI2026

Interplay between Ultrafast Electronic and Librational Dynamics in Liquid Nitrobenzene Probed with Two-Color Four-Wave Mixing

Here, we present an experimental and theoretical study of the interplay between ultrafast electron dynamics and librational dynamics in liquid nitrobenzene. A femtosecond ultraviolet pulse and two femtosecond near-infrared pulses interact with nitrobenzene molecules, generating a four-wave mixing nonlinear signal measured in the Optical Kerr Effect geometry. The signal is measured to be nonzero only at negative time delays, corresponding to the near-infrared pulses arriving before the ultraviolet pulse. We perform time-dependent Quantum Master Equation calculations with classical libration to simulate the experiment. The simulations support the conclusion that the near-infrared pulses launch librational motion while creating electronic coherences resulting in a libration-modulated electronic nonlinear response. The analysis of the phase-matched four-wave mixing signals suggests a nonparametric process leaving the molecules in an excited electronic state, providing new insight into ultrafast nonlinear optical interactions in liquids and advancing toward probing ultrafast electronic coherences in complex molecular liquids.

Shivaram, Niranjan [Lawrence Berkeley National Lab↗
DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Institutional Computing: Annual Progress Report, Precision early universe simulations to constrain nuclear reactions and beyond standard model physics

The current one-year allocation (w20 qburst) is the first year in a two year allocation of the same name. The project has supported two publications in peer-review journals, another manuscript in preparation for publication, a white paper for the Astro2020 Decadal Survey, and a number of talks (see below). Code development continues on implementing a neutrino quantum kinetic equation solver into BURST, now called QBURST. While implementing coherent and forward-scattering processes into the dynamical neutrino evolution equations, multiple numerical issues have appeared. We have remedied a number of these issues through extensive debugging and testing, which include streamlining interpolation and implementing new ODE integration techniques. Further testing is required to solve the problem at a precision of better than one part in 10 6 . We expect that we will understand the numerical difficulties and implement solutions by end of year 2021.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

ADoSPSiQS (Automated Detection of Symmetry-Protected Subspaces in Quantum Simulations ) [SWR-23-55]

Symmetry-protected subspaces in quantum systems are subspaces demarcated by a conserved quantity of that quantum system. We introduce two classical algorithms, which efficiently compute and elucidate features of these subspaces. The first algorithm explores the entire symmetry-protected subspace of an initial state in time complexity linear to the size of the subspace by closing local basis state-to-basis state transitions. The second algorithm determines, with bounded error, if a given measurement outcome of a dynamically-generated state is within the symmetry-protected subspace of the state in which the dynamical system is initialized. This software is paired with the paper of the same name at https://arxiv.org/pdf/2302.08586.pdf.

Rotello, Caleb↗
DOE OSTI2025

Dominant balance-based adaptive mesh refinement for incompressible fluid flows

This work introduces a novel adaptive mesh refinement (AMR) method that utilizes dominant balance analysis (DBA) for efficient and accurate grid adaptation in computational fluid dynamics (CFD) simulations. The proposed method leverages a Gaussian mixture model (GMM) to classify grid cells into active and passive regions based on the dominant physical interactions within the equation space. By modeling truncation error probabilistically from discretized terms, the method identifies regions of high interaction where numerical accuracy is most sensitive to resolution. Unlike traditional AMR strategies, this approach does not rely on heuristic-based sensors or user-defined thresholds, providing a fully automated and problem-independent framework for AMR. Applied to the incompressible Navier-Stokes equations for steady and unsteady flow past a cylinder, the DBA-based AMR method achieves comparable accuracy to high-resolution grids while reducing computational costs by up to 70 %. The validation highlights the method’s effectiveness in capturing complex flow features while minimizing grid cells, directing computational resources toward regions with the most critical dynamics. This modular and scalable strategy is adaptable to a wide range of applications, presenting a promising tool for efficient high-fidelity simulations in CFD and other multiphysics domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

High-throughput spin-bath characterization of spin defects in semiconductors

Detailed knowledge of the local environments of spin defects in semiconductors, such as nitrogenvacancy (NV) centers in diamond or divacancies in silicon carbide, is crucial for optimizing control and entanglement protocols in quantum sensing and information applications. However, at present a direct experimental characterization of individual defect environments is not scalable, as conventional spin-bath measurements are time consuming and difficult to automate. Achieving high-throughput characterization requires short experiments to probe the spin bath. However, with fewer and noisier measurements, the inverse problem of recovering spin-bath properties from measured data becomes ill posed, with multiple spin baths having a high likelihood of yielding the same data. In this work, we present a set of computational tools to resolve the ill-posed inverse problem of recovering the atomic positions and hyperfine couplings of random nuclei surrounding spin defects from sparse, noisy experimental coherence data, which can be obtained in hours. Here, we use a trans-dimensional Bayesian approach that incorporates ab initio data to yield full posterior distributions over nuclear spin environments, enabling robust recovery from limited data. We also provide practical tools and guidelines to determine the limits of detectability for hyperfine couplings under specific dynamical decoupling sequences and sampling conditions. In addition, we demonstrate how the tools developed here, in combination with ab initio simulations of spin baths, can guide the design of efficient experimental protocols for application-specific high-throughput screening. To showcase the utility of our approach, we apply it to design fast dynamical decoupling experiments to characterize the spin baths often individual NV centers in diamond. While the primary focus is on accelerating spin-bath characterization of spin defects, this Bayesian approach also lays the foundation for digital-twin studies of spin defects, where a virtual model of the spin-defect system evolves in real time with ongoing experimental measurements. Together, the set of tools we designed and applied paves the way for scalable deployment of spin defects in semiconductors for quantum sensing and information applications.

Bayesian methods↗
DOE OSTI2022

Transport model comparison studies of intermediate-energy heavy-ion collisions

Transport models are the main method to obtain physics information on the nuclear equation of state and in-medium properties of particles from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. To this end, calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, which test separately selected ingredients of a transport code, and more realistic calculations of heavy-ion collisions. Over the years, six studies have been performed within this project. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. After a compact description of the underlying transport approaches, we review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known under the names of Boltzmann–Uehling–Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. Furthermore, our evaluation studies often indicate improved strategies in performing transport simulations and thus can provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗
DOE OSTI2021

On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity

The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solving the Hele–Shaw flow using the Harrow–Hassidim–Lloyd algorithm on superconducting devices: A study of efficiency and challenges

The development of quantum processors for practical fluid flow problems is a promising yet distant goal. Recent advances in quantum linear solvers have highlighted their potential for classical fluid dynamics. In this study, we evaluate the Harrow–Hassidim–Lloyd (HHL) quantum linear systems algorithm (QLSA) for solving the idealized Hele–Shaw flow. Our focus is on the accuracy and computational cost of the HHL solver, which we find to be sensitive to the condition number, scaling exponentially with problem size. This emphasizes the need for preconditioning to enhance the practical use of QLSAs in fluid flow applications. Moreover, we perform shots-based simulations on quantum simulators and test the HHL solver on superconducting quantum devices, where noise, large circuit depths, and gate errors limit performance. Error suppression and mitigation techniques improve accuracy, suggesting that such fluid flow problems can benchmark noise mitigation efforts. Finally, our findings provide a foundation for future, more complex application of QLSAs in fluid flow simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Room-temperature skyrmion lattice in a layered magnet (Fe 0.5 Co 0.5 ) 5 GeTe 2

Novel magnetic ground states have been stabilized in two-dimensional (2D) magnets such as skyrmions, with the potential next-generation information technology. Here, we report the experimental observation of a Néel-type skyrmion lattice at room temperature in a single-phase, layered 2D magnet, specifically a 50% Co-doped Fe 5 GeTe 2 (FCGT) system. The thickness-dependent magnetic domain size follows Kittel's law. The static spin textures and spin dynamics in FCGT nanoflakes were studied by Lorentz electron microscopy, variable-temperature magnetic force microscopy, micromagnetic simulations, and magnetotransport measurements. Current-induced skyrmion lattice motion was observed at room temperature, with a threshold current density, j th = 1 × 10 6 A/cm 2 . This discovery of a skyrmion lattice at room temperature in a noncentrosymmetric material opens the way for layered device applications and provides an ideal platform for studies of topological and quantum effects in 2D.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

I ntera C hem : Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics

InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree–Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.

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

Realistic Ion Dynamics through Charge Renormalization in Nonaqueous Electrolytes

While many practically important electrolytes contain lithium ions, interactions of these ions are particularly difficult to probe experimentally because of their small X-ray and neutron scattering cross sections and large neutron absorption cross sections. Molecular dynamics (MD) is a powerful tool for understanding the properties of nonaqueous electrolyte solutions from the atomic level, but the accuracy of this computational method crucially depends on the physics built into the classical force field. In this study, we demonstrate that several force fields for lithium bistriflimide (LiTFSI) in acetonitrile yield a solution structure that is consistent with the neutron scattering experiments, yet these models produce dramatically different ion dynamics in solution. Such glaring discrepancies indicate that inadequate representation of long-range interactions leads to excessive ionic association and ion-pair clustering. We show that reasonable agreement with the experimental observations can be achieved by renormalization of the ion charges using a “titration” method suggested herewith. This simple modification produces realistic concentration dependencies for ionic diffusion and conductivity in <2 M solutions, without loss in quality for simulation of the structure.

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

Reaction nanoscopy of ion emission from sub-wavelength propanediol droplets

Abstract Droplets provide unique opportunities for the investigation of laser-induced surface chemistry. Chemical reactions on the surface of charged droplets are ubiquitous in nature and can provide critical insight into more efficient processes for industrial chemical production. Here, we demonstrate the application of the reaction nanoscopy technique to strong-field ionized nanodroplets of propanediol (PDO). The technique’s sensitivity to the near-field around the droplet allows for the in-situ characterization of the average droplet size and charge. The use of ultrashort laser pulses enables control of the amount of surface charge by the laser intensity. Moreover, we demonstrate the surface chemical sensitivity of reaction nanoscopy by comparing droplets of the isomers 1,2-PDO and 1,3-PDO in their ion emission and fragmentation channels. Referencing the ion yields to gas-phase data, we find an enhanced production of methyl cations from droplets of the 1,2-PDO isomer. Density functional theory simulations support that this enhancement is due to the alignment of 1,2-PDO molecules on the surface. The results pave the way towards spatio-temporal observations of charge dynamics and surface reactions on droplets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗