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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamics of collision-induced energy transfer

Collisions are of fundamental importance in many relevant physico-chemical phenomena in the gas phase and at the gas-surface interface. Collision-induced energy transfer not only offers a sensitive probe of interatomic interactions, but also impacts an array of processes ranging from transport to reactivity. This perspective highlights some recent advances on collision-induced energy transfer dynamics. In the gas phase, cold collisions revealed remarkable quantum effects and stereodynamics of bimolecular inelastic scattering, which motivated development of new state-to-state quantum scattering methods. For surface processes, experiments of hydrogen atom scattering from various surfaces combined with first-principles molecular dynamics simulations shed valuable light on the relative importance of adiabatic energy transfer to surface phonons and nonadiabatic energy transfer to surface electrons. We outline further challenges in these fields and expect more fruitful interplay between experiment and theory.

Cold collision↗
DOE OSTI2026

Single-Molecule Electron Transport in Peptoids

Peptoids are structural analogs of peptides in which side chains are appended to the backbone nitrogen rather than the α-carbon. The sequence-defined modularity of peptoids enables precise control over structure−function relationships, enabling applications in energy storage and biomedical materials. Despite recent progress, the role of sequence and conformation on electron transport in peptoid molecules is not fully understood. Here, we synthesize a library of peptoid oligomers and characterize their molecular electronic properties using the scanning tunneling microscope-break junction (STM-BJ) technique. Our results show well-defined electron transport behavior for peptoid sequences containing aromatic side groups lacking hydrogen bonds (H-bonds) and without chemical substitutions at the N−C α position. This behavior fundamentally differs from electron transport in peptides, where H-bond interactions give rise to higher conductance states. All-atom molecular dynamics (MD) simulations are used to understand the conformational heterogeneity of peptoids, and molecular conformations obtained from MD simulations are used in quantum mechanical calculations based on the nonequilibrium Green’s function−density functional theory (NEGF-DFT) formalism. In all cases, computational results are in reasonable qualitative agreement with experiments. Our work demonstrates that the conductance behavior of peptoids depends on monomer identity, including side-chain aromaticity and substitution at the N−C α position. Overall, this work provides new insights into the structure−function relationships governing electron transport in peptoid-based materials and establishes design rules for peptoid-based molecular junctions.

Charge transport↗
DOE OSTI2021

Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator

The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3–4 orders of magnitude fewer Hamiltonian evaluations and 2–3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2–3/1–2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ~1.3×10 9 parameters to obtain accurate dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Quantum Decision Maker Theory and Simulation

A quantum device simulating the human decision making process is introduced. It consists of quantum recurrent nets generating stochastic processes which represent the motor dynamics, and of classical neural nets describing the evolution of probabilities of these processes which represent the mental dynamics.

Quantum↗
DOE OSTI2021

Surface chemistry models for GaAs epitaxial growth and hydride cracking using reacting flow simulations

Hydride vapor phase epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultrafast growth rates (~500 µm/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is necessary to understand the kinetics of the growth process and its coupling with transport phenomena, so as to enable fast and uniform epitaxial growth. In this work, we derive a kinetic model using experimental data and integrate it into a computational fluid dynamics simulation of an HVPE growth reactor. We also modify an existing hydride cracking model that we validate against numerical simulations and experimental data. Here, we show that the developed growth model and the improved cracking model are able to reproduce experimental growth measurements of GaAs in an existing HVPE system.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

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

Deuteration removes quantum dipolar defects from KDP crystals

Abstract Dielectric properties of the hydrogen-bonded ferroelectric crystal KH 2 PO 4 (KDP) differ significantly from those of KD 2 PO 4 (DKDP). It is well established that deuteration affects the interplay of hydrogen-bond switches and heavy ion displacements that underlie the emergence of macroscopic polarization, but a detailed microscopic model is missing. We show that all-atom path integral molecular dynamics simulations can predict the isotope effects, revealing the microscopic mechanism that differentiates KDP and DKDP. Proton tunneling generates phosphate configurations that do not contribute to the polarization. At low temperatures, these quantum dipolar defects are substantial in KDP but negligible in DKDP. These intrinsic defects explain why KDP has lower spontaneous polarization and transition entropy than DKDP. The prominent role of quantum fluctuations in KDP is related to the unusual strength of the hydrogen bonds and should be equally important in other crystals of the KDP family, which exhibit similar isotope effects.

Yang, Bingjia (ORCID:0000000340749553)↗
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Attosecond XUV probing of vibronic quantum superpositions in Br 2 +

Ultrafast laser excitation can create coherent superpositions of electronic states in molecules and trigger ultrafast flow of electron density on a few femtosecond timescale. While recent attosecond experiments have addressed real-time observation of these primary photochemical processes, the underlying roles of simultaneous nuclear motions and how they modify and disturb the valence electronic motion remain uncertain. In this work, we investigate coherent electronic-vibrational dynamics induced among multiple vibronic levels of ionic bromine (Br 2 + ), including both spin-orbit (X 2 Π 3/2,g – X 2 Π 1/2,g ) and valence (A 2 Π 3/2,u – 4 Σ 3/2,u ) electronic superpositions, using attosecond transient absorption spectroscopy. Decay, revival, and apparent frequency shifts of electronic coherences are measured via characteristic quantum beats on the Br-3d core-level absorption signals. Quantum-mechanical simulations attribute the observed electronic decoherence to broadened phase distributions of nuclear wave packets on anharmonic potentials. Molecular vibronic structure is further revealed to be imprinted as discrete progressions in electronic beat frequencies. These results provide a future basis to interpret complex charge-migration dynamics in polyatomic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Non-Kasha fluorescence of pyrene emerges from a dynamic equilibrium between excited states

Pyrene fluorescence after a high-energy electronic excitation exhibits a prominent band shoulder not present after excitation at low energies. The standard assignment of this shoulder as a non-Kasha emission from the second-excited state (S 2 ) has been recently questioned. To elucidate this issue, we simulated the fluorescence of pyrene using two different theoretical approaches based on vertical convolution and nonadiabatic dynamics with nuclear ensembles. To conduct the necessary nonadiabatic dynamics simulations with high-lying electronic states and deal with fluorescence timescales of about 100 ns of this large molecule, we developed new computational protocols. The results from both approaches confirm that the band shoulder is, in fact, due to S 2 emission. We show that the non-Kasha behavior is a dynamic-equilibrium effect not caused by a metastable S 2 minimum. However, it requires considerable vibrational energy, which can only be achieved in collisionless regimes after transitions into highly excited states. This strict condition explains why the S 2 emission was not observed in some experiments.

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

Testing strong-field QED with the avalanche precursor

A two-beam high-power laser facility is essential for the study of one of the most captivating phenomena predicted by strong-field quantum electrodynamics (QED) and yet unobserved experimentally: the avalanche-type cascade. In such a cascade, the energy of intense laser light can be efficiently transformed into high-energy radiation and electron–positron pairs. The future 50-petawatt-scale laser facility NSF OPAL will provide unique opportunities for studying such strong-field QED effects, as it is designed to deliver two ultra-intense, tightly focused laser pulses onto the interaction point. In this work, we investigate the potential of such a facility for studying elementary particle and plasma dynamics deeply in the quantum radiation-dominated (RD) regime, and the generation of QED avalanches. With three-dimensional particle-in-cell simulations, we demonstrate that QED avalanche precursors can be reliably triggered under realistic laser parameters and layout (namely, focusing f/2, tilted optical axes, and non-ideal co-pointing) with the anticipated capabilities of NSF OPAL. We demonstrate that seed electrons can be efficiently injected into the laser focus by using targets of three types: a gas of heavy atoms, an overcritical plasma, and a thin foil. A strong positron and high-energy photon signal is generated in all cases. The cascade properties can be identified from the final particle distributions, which have a clear directional pattern. At increasing laser field intensity, such distributions provide signatures of the transition, first, to the RD interaction regime, and then to a QED avalanche. Our findings can also be used for designing related future experiments.

Lasers↗
DOE OSTI2021

Machine learning the dynamics of quantum kicked rotor

Highlights: • Quantum kicked rotor. • Quantum phase transition. • Machine learning. • LSTM network. • Convolutional neural network. Using the multilayer convolutional neural network (CNN), we can detect the quantum phases in random electron systems, and phase diagrams of two and higher dimensional Anderson transitions and quantum percolations as well as disordered topological systems have been obtained. Here, instead of using CNN to analyze the wave functions, we analyze the dynamics of wave packets via long short-term memory network (LSTM). We adopt the quasi-periodic quantum kicked rotors, which simulate the three and four dimensional Anderson transitions. By supervised training, we let LSTM extract the features of the time series of wave packet displacements in localized and delocalized phases. We then simulate the wave packets in unknown phases and let LSTM classify the time series to localized and delocalized phases. We compare the phase diagrams obtained by LSTM and those obtained by CNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

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

Generative deep-learning reveals collective variables of Fermionic systems

Complex processes of fermionic systems ranging from protein folding to nuclear fission often follow a low-dimensional reaction path parametrized in terms of a few collective variables. In nuclear theory, variables related to the shape of the nuclear density in a mean-field picture are key to describing the large amplitude collective motion of the neutrons and protons. Exploring the adiabatic energy landscape spanned by these degrees of freedom reveals the possible reaction channels while simulating the dynamics in this reduced space yields their respective probabilities. Unfortunately, this theoretical framework breaks down whenever the systems encounters a quantum phase transition with respect to the collective variables. Here, in this study, we introduce a novel generative deep-learning algorithm designed to build reaction paths that ensure that the many-fermion wave function stays differentiable with respect to the collective variables. This approach is applicable to any fermionic system described by a coherent state. We use the case of potential energy curves in the 16 O nucleus within the Hartree-Fock theory to illustrate its main features.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Hidden quantum criticality and entanglement in quench dynamics

Entanglement exhibits universal behavior near the ground-state critical point where correlations are long ranged and the thermodynamic entropy is vanishing. On the other hand, a quantum quench imparts extensive energy and results in a build up of entropy, hence no critical behavior is expected at long times. In this work, we present a new paradigm in the quench dynamics of integrable spin chains which exhibit a ground-state order-disorder phase transition at a critical line. Specifically, we consider a quench along the critical line which displays a volume-law behavior of the entropy and exponentially decaying correlations; however, we show that quantum criticality is hidden in higher-order correlations and becomes manifest via measures such as the mutual information and logarithmic negativity. Furthermore, we showcase the scale invariance of the Rényi mutual information between disjoint regions as further evidence for genuine critical behavior. We attribute the emergent quantum criticality to the soft mode not getting excited in spite of the quench. Moreover, the results presented here are universal to models whose low-energy or long-wavelength dynamics are well described by a free-fermionic field theory. Our results are amenable to an experimental realization on different quantum simulator platforms, particularly the Rydberg simulators. Published by the American Physical Society 2024

Paul, Sanku↗
DOE OSTI2024

Extraction of the electron excess temperature in terahertz quantum cascade lasers from laser characteristics

Abstract We propose a method to extract the upper laser level’s (ULL’s) excess electronic temperature from the analysis of the maximum light output power ( P max ) and current dynamic range Δ J d = ( J max − J th ) of terahertz quantum cascade lasers (THz QCLs). We validated this method, both through simulation and experiment, by applying it on THz QCLs supporting a clean three-level system. Detailed knowledge of electronic excess temperatures is of utmost importance in order to achieve high temperature performance of THz QCLs. Our method is simple and can be easily implemented, meaning an extraction of the excess electron temperature can be achieved without intensive experimental effort. This knowledge should pave the way toward improvement of the temperature performance of THz QCLs beyond the state-of-the-art.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗