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 811 records · Page 45

Exact-Factorization-Based Surface Hopping without Velocity Adjustment

While surface hopping has emerged as a powerful method for simulating non-adiabatic dynamics in large molecules, the ad hoc nature of the necessary velocity adjustments and decoherence corrections in the algorithm somewhat reduces its reliability. Here we propose a new scheme that eliminates these aspects by combining the nuclear equation from the quantum-trajectory surface-hopping approach with the electronic equation derived from the exact-factorization approach. Furthermore, the resulting method, denoted QTSH-XF, yields a surface-hopping method on firmer ground than previous and is shown to successfully capture dynamics in Tully models and in a linear vibronic coupling model of the photoexcited uracil cation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Perpendicular crossing chains enable high mobility in a noncrystalline conjugated polymer

The nature of interchain π-system contacts, and their relationship to hole transport, are elucidated for the high-mobility, noncrystalline conjugated polymer C16-IDTBT by the application of scanning tunneling microscopy, molecular dynamics, and quantum chemical calculations. The microstructure is shown to favor an unusual packing motif in which paired chains cross-over one another at near-perpendicular angles. By linking to mesoscale microstructural features, revealed by coarse-grained molecular dynamics and previous studies, and performing simulations of charge transport, it is demonstrated that the high mobility of C16-IDTBT can be explained by the promotion of a highly interconnected transport network, stemming from the adoption of perpendicular contacts at the nanoscale, in combination with fast intrachain transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the equivalence of nonequilibrium and equilibrium measurements of slip in molecular dynamics simulations

Here, we show that nonequilibrium and equilibrium measurements of slip are consistent, provided the hydrodynamic wall location associated with the equilibrium measurement is properly taken into account. The latter is a strong function of the fluid state and wall-fluid interaction and cannot be neglected as it typically has been. Our results are based on an alternative approach for calculating the hydrodynamic wall location which alleviates most of the difficulties associated with its calculation via a Green-Kubo integral that appear to have contributed to its neglect. Extensive molecular dynamics simulations are used to validate our approach including a model for calculating the slip length that does not involve a Green-Kubo integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Floquet-engineered fast SNAP gates in weakly coupled circuit-QED systems

Superconducting cavities with high quality factors, coupled to a fixed-frequency transmon, provide a state-of-the-art platform for quantum information storage and manipulation. The commonly used selective number-dependent arbitrary phase ( SNAP ) gate faces significant challenges in ultrahigh-coherence cavities, where the weak dispersive shifts necessary for preserving high coherence typically result in prolonged gate times. Here, in this work, we propose a protocol to achieve high-fidelity SNAP gates that are orders of magnitude faster than the standard implementation, surpassing the speed limit set by the bare dispersive shift. We achieve this enhancement by dynamically amplifying the dispersive coupling via sideband interactions, followed by quantum optimal control on the Floquet-engineered system. We also present a unified perturbation theory that explains both the gate acceleration and the associated benign drive-induced decoherence, corroborated by Floquet-Markov simulations. These results pave the way for the experimental realization of high-fidelity, selective control of weakly coupled, high-coherence cavities, and expanding the scope of optimal control techniques to a broader class of Floquet quantum systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CFD Analysis of RLUOB Zone 1 HEPA Filter Plenum and Testing Manifolds

The purpose of this computational fluid dynamics (CFD) analysis is to ensure that Camfil Farr’s upstream and downstream injection and sampling manifolds can meet or exceed the requirements outlined in the ASME AG-1 1997 a(2000) Code on Nuclear Air and Gas Treatment for the testing of HEPA and adsorbent filters. This paper will present a numerical simulation of airflow in the Radiological Laboratory Utility Office Building (RLUOB) zone 1 HEPA filter plenum and testing manifolds using the commercial CFD software ANSYS FLUENT 2020R1. The CFD analysis focuses on the investigation of the air flow distribution and air-aerosol mixing uniformity. The evaluation was done for all steps of the modeling process: grid generation, physics setup, simulation, and post-processing. The mass flow rate in each section of the zone 1 injection and sampling manifolds is also reported.

42 ENGINEERING↗

A-B Transition in Superfluid $^3$He and Cosmological Phase Transitions

First-order phase transitions in the very early universe are a prediction of many extensions of the Standard Model of particle physics and could provide the departure from equilibrium needed for a dynamical explanation of the baryon asymmetry of the Universe. They could also produce gravitational waves of a frequency observable by future space-based detectors such as the Laser Interferometer Space Antenna. All calculations of the gravitational wave power spectrum rely on a relativistic version of the classical nucleation theory of Cahn-Hilliard and Langer, due to Coleman and Linde. The high purity and precise control of pressure and temperature achievable in the laboratory made the first-order A to B transition of superfluid $^3$He ideal for test of classical nucleation theory. As Leggett and others have noted, the theory fails dramatically. The lifetime of the metastable A phase is measurable, typically of order minutes to hours, far faster than classical nucleation theory predicts. If the nucleation of B phase from the supercooled A phase is due to a new, rapid intrinsic mechanism that would have implications for first-order cosmological phase transitions as well as predictions for gravitational wave production in the early universe. Here we discuss studies of the A-B phase transition dynamics in $^3$He, both experimental and theoretical, and show how the computational technology for cosmological phase transition can be used to simulate the dynamics of the A-B transition, support the experimental investigations of the A-B transition in the QUEST-DMC collaboration with the goal of identifying and quantifying the mechanism(s) responsible for nucleation of stable phases in ultra-pure metastable quantum phases.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum fluctuations drive nonmonotonic correlations in a qubit lattice

Abstract Fluctuations may induce the degradation of order by overcoming ordering interactions, consequently leading to an increase of entropy. This is particularly evident in magnetic systems characterized by nontrivial, constrained disorder, where thermal or quantum fluctuations can yield counterintuitive forms of ordering. Using the proven efficiency of quantum annealers as programmable spin system simulators, we present a study based on entropy postulates and experiments on a platform of programmable superconducting qubits to show that a low level of uncertainty can promote ordering in a system impacted by both thermal and quantum fluctuations. A set of experiments is proposed on a lattice of interacting qubits arranged in a triangular geometry with precisely controlled disorder, effective temperature, and quantum fluctuations. Our results demonstrate the creation of ordered ferrimagnetic and layered anisotropic disordered phases, displaying characteristics akin to the elegant order-by-disorder phenomenon. Extensive experimental evidence is provided for the role of quantum fluctuations in lowering the total energy of the system by increasing entropy and defect clustering. Our thorough and comprehensive application of an intentionally introduced noise on a quantum platform provides insight into the dynamics of defects and fluctuations in quantum devices, which may help to reduce the cost associated with quantum processing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High-speed x-ray phase contrast imaging and digital image correlation analysis of microscale shock response of an additively manufactured energetic material simulant

The performance of energetic materials subjected to dynamic loading significantly depends on their micro- and meso-scale structural morphology. The geometric versatility offered by additive manufacturing opens new pathways to tailor the performance of these materials. Additively manufactured energetic materials (AMEMs) have a wide range of structural characteristics with a hierarchy of length scales and process-inherent heterogeneities, which are hitherto difficult to precisely control. It is important to understand how these features affect AMEMs’ response under dynamic/shock loading. Therefore, temporally and spatially resolved measurements of both macroscopic behavior and micro- and meso-level processes influencing macroscopic behavior are required. In this paper, we analyze the shock compression response of an AMEM simulant loaded under several impact conditions and orientations. Furthermore, x-ray phase contrast imaging (PCI) is used to track features across the observed shock front and determine the linear shock velocity vs particle velocity equation of state, as well as to quantify the interior deformation fields via digital image correlation (DIC) analyses. Photon Doppler velocimetry is simultaneously used to measure the particle velocities of the specimens, which are consistent with those obtained from x-ray PCI. The DIC analyses provide an assessment of the average strain fields inside the material, showing that the average axial strain depends on the loading intensity and reaches as high as 0.23 for impact velocities up to 1.5 km/s. The overall results demonstrate the utility of x-ray PCI for probing “in-material” equation of state and interior strains associated with dynamic shock compression behavior of the AMEM simulant.

3D printing↗

Towards a classification of bifurcations in Vlasov equations

Here, we propose a classification of bifurcations of Vlasov equations, based on the strength of the resonance between the unstable mode and the continuous spectrum on the imaginary axis. We then identify and characterize a new type of generic bifurcation where this resonance is weak, but the unstable mode couples with a stable mode and a Casimir invariant of the system to form a size-3 Jordan block. We derive a three-dimensional reduced noncanonical Hamiltonian system describing this bifurcation. Comparison of the reduced dynamics with direct numerical simulations on a test case gives excellent agreement. We finally discuss the relevance of this bifurcation to specific physical situations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mechanical Resonant Sensing of Spin Texture Dynamics in a 2D Antiferromagnet

Abstract The coupling between the spin degrees of freedom and macroscopic mechanical motions, including striction, shearing, and rotation, has attracted wide interest with applications in actuation, transduction, and information processing. Experiments so far have established the mechanical responses to the long‐range ordered or isolated single spin states. However, it remains elusive whether mechanical motions can couple to a different type of magnetic structure, the non‐collinear spin textures, which exhibit nanoscale spatial variations of spin (domain walls, skyrmions,etc.) and are promising candidates to realize high‐speed computing devices. Here, collective spin texture dynamics is detected with nanoelectromechanical resonators fabricated from 2D antiferromagnetic (AFM) MnPS 3 with 10 −9 strain sensitivity. By examining radio frequency mechanical oscillations under magnetic fields, new magnetic transitions are identified with sharp dips in resonant frequency. They are attributed to collective AFM domain wall motions as supported by the analytical modeling of magnetostriction and large‐scale spin‐dynamics simulations. Additionally, an abnormally large modulation in the mechanical nonlinearity at the transition field infers a fluid‐like response due to ultrafast domain motion. The work establishes a strong coupling between spin texture and mechanical dynamics, laying the foundation for electromechanical manipulation of spin texture and developing quantum hybrid devices.

Chemistry↗

Microgap breakdown with floating metal rod perturbations

Here, we report the characterization of microgap breakdown with perturbations from a metal rod floating between anode and cathode electrodes. The effects of the metal rod on the electric field distribution and the field enhancement factor are evaluated by numerical simulation and the conformal mapping method, and they indicate that the field emission regime is not reached. The breakdown voltages in the Townsend discharge regime are determined based on the voltage–current characteristics, which are obtained from two-dimensional fluid simulations. It is found that the breakdown characteristics can be significantly modulated by the floating metal rod, and the breakdown curve (breakdown voltage vs the net gap distance) is no longer U-shaped, which deviates from the conventional Paschen's law. The underlying physical mechanisms are related to the electric field enhancement, curved breakdown path, and nonuniform ion flux caused by the electric shielding effect. The results provide insights into breakdown characteristics in microscale discharges, which may promote conventional investigation of simplified clean gaps toward more complex conditions (e.g., with floating microparticles) in miniaturized plasma devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Statistical Properties of Undulator Radiation

Two experiments were carried out to study the statistical properties of undulator radiation in the Integrable Optics Test Accelerator (IOTA) storage ring at Fermilab. The first experiment studied the turn-to-turn fluctuations in the power of the radiation generated by an electron bunch. The magnitude of these fluctuations depends on the 6D phase-space distribution of the electron bunch. In IOTA, we demonstrated that this effect can be used to measure some electron bunch parameters, small transverse emittances in particular. In the second experiment, a single electron was stored in the ring, emitting a photon only once per several hundred turns. In this regime, any classical interference-related collective effects were eliminated, and the quantum fluctuations could be studied in detail to search for possible deviations from the expected Poissonian photon statistics. In addition, the photocount arrival times were used to track the longitudinal motion of a single electron and to compare it with simulations. This allowed us to determine several dynamical parameters of the storage ring such as the rf cavity phase jitter and the dependence of the synchrotron motion period on amplitude.

43 PARTICLE ACCELERATORS↗

Enhancement of Photoemission on p-Type GaAs Using Surface Acoustic Waves

We demonstrate that photoemission properties of p-type GaAs can be altered by surface acoustic waves (SAWs) generated on the GaAs surface due to dynamical piezoelectric fields of SAWs. Multiphysics simulations indicate that charge-carrier recombination is greatly reduced, and electron effective lifetime in p-doped GaAs may increase by a factor of 10× to 20×. It implies a significant increase, by a factor of 2× to 3×, of quantum efficiency (QE) for GaAs photoemission applications, like GaAs photocathodes. Conditions of different SAW wavelengths, swept SAW intensities, and varied incident photon energies were investigated. Essential steps in SAW device fabrication on a GaAs substrate are demonstrated, including deposition of an additional layer of ZnO for piezoelectric effect enhancement, measurements of current–voltage (I–V) characteristics of the SAW device, and ability to survive high-temperature annealing. Results obtained and reported in this study provide the potential and basis for future studies on building SAW-enhanced photocathodes, as well as other GaAs photoelectric applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Chemical Simulations of CO 2 and N 2 Capture in Reline, a Prototypical Deep Eutectic Solvent

Deep eutectic solvents such as reline are an emerging class of low-cost, environmentally friendly solvents with tunable properties that are potentially applicable for the capture and separation of CO 2 . Experimental measurements showed that a reline-based membrane contactor can capture and separate CO 2 via physisorption through a dissolution process with 96.7% purity from a mixed gas containing CO 2 and N 2 (50:50% molar ratio). Here, we examine the nature of the interaction of CO 2 and N 2 with reline employing quantum chemical methods. We focus on explaining the mechanism by which CO 2 and N 2 bind to reline and the reason for the high selectivity for absorption of CO 2 compared to N 2 . We analyze the dynamics, energetics, and binding motifs for CO 2 and N 2 in reline employing density functional theory, density functional tight binding, and ab initio molecular dynamics. We also investigate the effect of reline on the vibrational spectra of CO 2 and reline. Our simulations indicate that the selective capture of CO 2 from the mixture of CO 2 and N 2 is due to the interplay between attractive electrostatic and charge polarization forces with opposing entropic effects, which shift the energetic balance and make the N 2 absorption unfavorable in reline.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

3D Cell-centered hydrodynamics with subscale closure model and multi-material remap

We extend a higher-order finite volume cell-centered hydrodynamic (CCH) formulation to include an interface-aware subscale closure model and a multi-material remap for simulating 3D compressible hydrodynamic problems within an arbitrary Lagrangian-Eulerian (ALE) framework. This CCH formulation involves a multidirectional approximate Riemann solution using quadratic polynomial reconstructions of the stress tensor and the velocity. At the subscale level, we determine pair-wise material interactions by solving a distinct approximate Riemann problem at the common interface, using the volume of fluids (VOF) method to find the interface. Material interactions are constrained to ensure smooth pressure equilibration among materials. The accuracy and robustness of the ALE method is demonstrated by simulating a suite of 3D Cartesian multi-material problems covering both gas and solid dynamics, where each test case has two or more materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy-efficient pathway for selectively exciting solute molecules to high vibrational states via solvent vibration-polariton pumping

Abstract Selectively exciting target molecules to high vibrational states is inefficient in the liquid phase, which restricts the use of IR pumping to catalyze ground-state chemical reactions. Here, we demonstrate that this inefficiency can sometimes be solved by confining the liquid to an optical cavity under vibrational strong coupling conditions. For a liquid solution of 13 CO 2 solute in a 12 CO 2 solvent, cavity molecular dynamics simulations show that exciting a polariton (hybrid light-matter state) of the solvent with an intense laser pulse, under suitable resonant conditions, may lead to a very strong (>3 quanta) and ultrafast (<1 ps) excitation of the solute, even though the solvent ends up being barely excited. By contrast, outside a cavity the same input pulse fluence can excite the solute by only half a vibrational quantum and the selectivity of excitation is low. Our finding is robust under different cavity volumes, which may lead to observable cavity enhancement on IR photochemical reactions in Fabry–Pérot cavities.

74 ATOMIC AND MOLECULAR PHYSICS↗