Search NASA⌕ Search

SEARCH · Search NASA

Results for “real-time lattice simulations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First-principles simulation of light-ion microscopy of graphene

The extreme sensitivity of 2D materials to defects and nanostructure requires precise imaging techniques to verify presence of desirable and absence of undesirable features in the atomic geometry. Helium-ion beams have emerged as a promising materials imaging tool, achieving up to 20 times higher resolution and 10 times larger depth-of-field than conventional or environmental scanning electron microscopes. Here, we offer first-principles theoretical insights to advance ion-beam imaging of atomically thin materials by performing real-time time-dependent density functional theory simulations of single impacts of 10–200 keV light ions in free-standing graphene. Here we predict that detecting electrons emitted from the back of the material (the side from which the ion exits) would result in up to three times higher signal and up to five times higher contrast images, making 2D materials especially compelling targets for ion-beam microscopy. This predicted superiority of exit-side emission likely arises from anisotropic kinetic emission. The charge induced in the graphene equilibrates on a sub-fs time scale, leading to only slight disturbances in the carbon lattice that are unlikely to damage the atomic structure for any of the beam parameters investigated here.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Microscopic force driving the photoinduced ultrafast phase transition: Time-dependent density functional theory simulations of IrTe 2

Photoinduced phase transitions can have complex and intriguing behaviors more than material ground-state dynamics. Understanding the underlying mechanism can help us to design new ways to manipulate the materials. A variety of mechanisms has been proposed to explain the photoinduced phase transitions of IrTe 2 , but a consensus has yet to be reached. Here, we study the photo-induced phase transitions of IrTe 2 by performing the real-time time-dependent density functional theory (rt-TDDFT) simulations in combination with the occupation constrained DFT method. We reveal that the microscopic driving force for the photo-induced phase transitions arises from the tendency to lower the energy levels occupied by the excited carriers, through the increase or decrease of the associated atomic pair distances, depending on whether the newly occupied states are antibonding or bonding states, respectively. The geometric constraints between different bonds represented by the Poisson ratio can bring together different tendencies from different atomic pairs, thus forming a complex intriguing dynamic picture depending on the intensity of the excitation. Furthermore, we also find that phonons don't play a primary role, but can assist the phase transition. These findings imply that one can control the structural phase transitions by selectively exciting photocarriers into designated atomic states using appropriate photon sources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electronically induced defect creation at semiconductor/oxide interface revealed by time-dependent density functional theory

Carrier induced defect creation at the semiconductor-oxide interface has been known as the origin of electronic device degradation for a long time, but how exactly the interface lattice can be damaged by carriers (especially low-energy ones) remains unclear. Here we carry out real-time time-dependent density functional theory simulations on concrete Si/SiO 2 interfaces to study the interaction between excited electrons and interface bonds. We show that the normal interface Si-H bonds are generally resistant to electrons due to the delocalized nature and high energy level of the Si-H antibonding states, and due to the high-energy barrier to break the Si-H bond. However, if an additional hydrogen atom exists by attaching to a nearby oxygen atom (forming a “Si-H···H-O” complex), the Si-H bond will be greatly weakened, including the reduction of energy barrier for bond breaking, and the lowering of the antibonding state energy level which favors electron injection. Together with the multiple vibrational excitation process, the corresponding Si-H bond can be broken much more easily. Thus we propose that the Si-H···H-O complex will be the center for defect creation and device degradation. Finally, we also explain why such a center might be relatively easy to form during the hydrogen annealing process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ab initio electron dynamics in high electric fields: Accurate prediction of velocity-field curves

Electron dynamics in external electric fields governs the behavior of solid-state electronic devices. First-principles calculations enable precise predictions of charge transport in low electric fields. However, studies of high-field electron dynamics remain elusive due to a lack of accurate and broadly applicable methods. In this work, we develop an efficient approach to solve the real-time Boltzmann transport equation with both the electric field term and ab initio electron-phonon collisions. These simulations provide field-dependent electronic distributions in the time domain, allowing us to investigate both transient and steady-state transport in electric fields ranging from low to high (> 10 kV/cm). The broad capabilities of our approach are shown by computing nonequilibrium electron occupations and velocity-field curves in Si, GaAs, and graphene, obtaining results in quantitative agreement with experiment. Our approach sheds light on microscopic details of transport in high electric fields, including the dominant scattering mechanisms and valley occupation dynamics. Our results demonstrate quantitatively accurate calculations of electron dynamics in low to high electric fields, with broad application to power and micro-electronics, optoelectronics, and sensing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Dynamic short-range correlation in photoinduced disorder phase transitions

Ultrafast photoexcitation can induce a nonequilibrium dynamic with electron-lattice interaction, offering an effective way to study photoinduced phase transitions (PIPTs) in solids. The issue that atomic displacements after photoexcitation belong to a coherent change or disordered process has become controversial in the PIPT community. Using real-time, time-dependent density functional theory (rt-TDDFT) simulations, we were able to obtained both the coherent and the disordered PIPTs (dimer dissociation) in IrTe 2 with different electronic occupations. More importantly, we found that in the disordered phase transition there exists a local correlation between different dimers regarding their dissociation status. We define these Ir-Ir dimers directly connected by Te atoms, including intralayer and vertically across the layers, as a group (group I). Other Ir-Ir dimers separated by five Ir atoms from Ir-Ir dimers in group I are divided into another group (group II). The dimers in the same group will dissociate in a correlated fashion; they either all dissociate or all do not dissociate. On the other hand, the dimers in neighboring groups will have an anticorrelation: If the dimers in one group dissociate, the dimers in the neighboring group tend not to be dissociated, and vice versa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Steps toward quantum simulations of hadronization and energy loss in dense matter

A framework for simulating the real-time dynamics of composite particles in a simple model of dense matter that is amenable to quantum computers is developed. As a demonstration, we perform classical simulations of heavy-hadrons propagating through a dense medium in the Schwinger model. Measurements of the time-dependent energy and charge density are used to identify mechanisms responsible for energy loss and hadron production (hadronization). A study of entanglement dynamics highlights the importance of quantum coherence between the particles that make up the dense medium. Throughout this work, care is taken to isolate, and remove, phenomena that arise solely from a finite lattice spacing. It is found that signatures of entanglement are more sensitive to lattice artifacts than other observables. Toward quantum simulations, we present an efficient method and the corresponding quantum circuits for preparing ground states in the presence of heavy mesons. Finally, these circuits are used to estimate the resources required to simulate in-medium energy loss and hadronization in the Schwinger model using quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Real-time quantum calculations of phase shifts using wave packet time delays

We present a method to extract the phase shift of a scattering process using the real-time evolution in the early and intermediate stages of the collision in order to estimate the time delay of a wave packet. This procedure is convenient when using noisy quantum computers for which the asymptotic out-state behavior is unreachable. We demonstrate that the challenging Fourier transforms involved in the state preparation and measurements can be implemented in 1+1 dimensions with current trapped ion devices and IBM quantum computers. We compare quantum computations of the time delays obtained in the one-particle quantum mechanics limit and the scalable quantum field theory formulation with accurate numerical results. We discuss the finite volume effects in the Wigner formula connecting time delays to phase shifts. The results reported involve two- and four-qubit calculations, and we discuss the possibility of larger scale computations in the near future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Nuclear two point correlation functions on a quantum computer

The calculation of dynamic response functions is expected to be an early application benefiting from rapidly developing quantum hardware resources. The ability to calculate real-time quantities of strongly correlated quantum systems is one of the most exciting applications that can easily reach beyond the capabilities of traditional classical hardware. Response functions of fermionic systems at moderate momenta and energies corresponding roughly to the Fermi energy of the system are a potential early application because the relevant operators are nearly local, and the energies can be resolved in moderately short real time, reducing the spatial resolution and gate depth required. This is particularly the case in quasielastic electron and neutrino scattering from nuclei, a topic of great interest in the nuclear and particle physics communities and directly related to experiments designed to probe neutrino properties. In this paper we use current quantum hardware and error mitigation protocols to calculate response functions for a highly simplified nuclear model through calculations of a 2-point real time correlation function for a modified Fermi-Hubbard model in two dimensions with three distinguishable nucleons on four lattice sites.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Liouvillian dynamics of the open Schwinger model: String breaking and kinetic dissipation in a thermal medium

Understanding the dynamics of bound state formation is one of the fundamental questions in confining quantum field theories such as Quantum Chromodynamics (QCD). One hadronization mechanism that has garnered significant attention is the breaking of a string initially connecting a fermion and an antifermion. Deepening our understanding of real-time string-breaking dynamics with simpler, lower dimensional models like the Schwinger model can improve our understanding of the hadronization process in QCD and other confining systems found in condensed matter and statistical systems. Here in this paper, we consider the string-breaking dynamics within the Schwinger model and investigate its modification inside a thermal medium, treating the Schwinger model as an open quantum system coupled to a thermal environment. Within the regime of weak coupling between the system and environment, the real-time evolution of the system can be described by a Lindblad evolution equation. We analyze the Liouvillian gaps of this Lindblad equation and the time dependence of the system’s von Neumann entropy. We observe that the late-time relaxation rate decreases as the environment correlation length increases. Moreover, when the environment correlation length is infinite, the system exhibits two steady states, one in each of the sectors with definite charge-conjugation-parity (CP) quantum numbers. For parameter regimes where an initial string breaks in vacuum, we observe a delay of the string breaking in the medium, due to kinetic dissipation effects. Conversely, in regimes where an initial string remains intact in vacuum time evolution, we observe string breaking (melting) in the thermal medium. We further discuss how the Liouvillian dynamics of the open Schwinger model can be simulated on quantum computers and provide an estimate of the associated Trotter errors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Quantum real-time evolution of entanglement and hadronization in jet production: Lessons from the massive Schwinger model

The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantum real-time evolution using tensor renormalization group methods

We introduce an approach for approximate real-time evolution of quantum systems using tensor renormalization group (TRG) methods originally developed for imaginary time. We use higher-order TRG to generate a coarse-grained time evolution operator for a 1+1⁢D transverse Ising model with a longitudinal field. We show that the standard tensor norm used for the singular value decomposition-based truncation is degenerate and propose an alternate method to discriminate. We show that it is effective and efficient in evolving Gaussian wave packets for one and two particles in the disordered phase, while ordered phase behavior is more challenging to capture. We compare our algorithm with local simulators for universal quantum computers and discuss possible benchmarking in the near future.

lattice gauge theory↗
DOE OSTI2023

Emergent Glassy Behavior in a Kagome Rydberg Atom Array

Here, we present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal emergent glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glassy region is demarcated using the Edwards-Anderson order parameter, and its phase transitions to the two proximate valence bond solids—as well as the crossover towards a trivial paramagnetic phase—are identified. We demonstrate the intrinsically slow (imaginary) time dynamics deep inside the glassy phase and discuss experimental considerations for detecting such a quantum disordered phase with numerous nearly degenerate local minima. Our proposal paves a new route to the study of real-time glassy phenomena and highlights the potential for quantum simulation of a distinct phase of quantum matter beyond solids and liquids in current-generation Rydberg platforms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Random forest prediction of crystal structure from electron diffraction patterns incorporating multiple scattering

Diffraction is the most common method to solve for unknown or partially known crystal structures. However, it remains a challenge to determine the crystal structure of a new material that may have nanoscale size or heterogeneities. Here, in this study, we train an architecture of hierarchical random forest models capable of predicting the crystal system, space group, and lattice parameters from one or more unknown two-dimensional electron diffraction patterns. Our initial model correctly identifies the crystal system of a simulated electron diffraction pattern from a 20-nm-thick specimen of arbitrary orientation 67% of the time. We achieve a topline accuracy of 79% when aggregating predictions from ten patterns of the same material but different zone axes. The space group and lattice predictions range from 70% to 90% accuracy and median errors of 0.01-0.5Å, respectively, for cubic, hexagonal, trigonal, and tetragonal crystal systems while being less reliable on orthorhombic and monoclinic systems. We apply this architecture to a four-dimensional scanning transmission electron microscopy scan of gold nanoparticles, where it accurately predicts the crystal structure and lattice constants. These random forest models can be used to significantly accelerate the analysis of electron diffraction patterns, particularly in the case of unknown crystal structures. Additionally, due to the speed of inference, these models could be integrated into live transmission electron microscopy experiments, allowing real-Time labeling of a specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2026

pathSQE : an automated workflow for single-crystal inelastic neutron scattering data processing and analysis

Inelastic neutron scattering (INS) experiments utilizing modern time-of-flight spectrometers enable the comprehensive mapping of the energy (E)- and momentum (Q)-resolved dynamical structure factor of single crystals, probing both the lattice and magnetic excitations. Yet, the large size and complexity of four-dimensional INS data are challenging current analysis workflows, often resulting in an underutilization of the measured information. To help address this issue, this paper introduces new software interfaced with the Mantid framework, pathSQE, designed to streamline the processing, analysis and interpretation of 4D single-crystal INS data. By automating key tasks such as 1D/2D slicing, symmetrization, Brillouin zone folding, data visualization, prioritization and filtering, and comparisons with simulations, pathSQE facilitates and accelerates INS data analysis workflows. Here, this paper outlines the features and implementation and provides several illustrations of the use of pathSQE on data collected on single crystals using direct-geometry time-of-flight spectrometers at the Spallation Neutron Source, including Ge, FeSi, MnO and SnS single-crystal measurements on the ARCS, HYSPEC and CNCS neutron spectrometers. Beyond streamlining post-experiment data processing, pathSQE establishes an automated and modular processing pipeline that could support future real-time experiment steering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The seeds and homogeneous nucleation of photoinduced nonthermal melting in semiconductors due to self-amplified local dynamic instability

Laser-induced nonthermal melting in semiconductors has been studied over the past four decades, but the underlying mechanism is still under debate. Here, by using an advanced real-time time-dependent density functional theory simulation, we reveal that the photoexcitation-induced ultrafast nonthermal melting in silicon occurs via homogeneous nucleation with random seeds originating from a self-amplified local dynamic instability. Because of this local dynamic instability, any initial small random thermal displacements of atoms can be amplified by a charge transfer of photoexcited carriers, which, in turn, creates a local self-trapping center for the excited carriers and yields the random nucleation seeds. Because a sufficient amount of photoexcited hot carriers must be cooled down to band edges before participating in the self-amplification of local lattice distortions, the time needed for hot carrier cooling is the response for the longer melting time scales at shorter laser wavelengths. This finding provides fresh insights into photoinduced ultrafast nonthermal melting.

42 ENGINEERING↗
NASA NTRS1989

The formation of topological defects in phase transitions

It was argued, and fought through numerical work that the results of non-dynamical Monte Carlo computer simulations cannot be applied to describe the formation of topological defects when the correlation length at the Ginzburg temperature is significantly smaller than the horizon size. To test the current hypothesis that infinite strings at formation are essentially described by Brownian walks of size the correlation length at the Ginzburg temperature, fields at the Ginzburg temperature were equilibrated. Infinite structure do not exist in equilibrium for reasonable definitions of the Ginzburg temperature, and horizons must be included in a proper treatment. A phase transition, from small-scale to large-scale string or domain wall structure, is found to occur very close to the Ginzburg temperature, in agreement with recent work. The formation process of domain walls and global strings were investigated through the breaking of initially ordered states. To mimic conditions in the early Universe, cooling times are chosen so that horizons exist in the sample volume when topological structure formation occurs. The classical fields are evolved in real-time by the numerical solution of Langevin equations of motion on a three dimensional spatial lattice. The results indicate that it is possible for most of the string energy to be in small loops, rather than in long strings, at formation.

Hodges, Hardy M.↗
DOE OSTI2020

Quantum Simulation of Field Theories Without State Preparation

We propose an algorithm for computing real-time observables using a quantum processor while avoiding the need to prepare the full quantum state. This reduction in quantum resources is achieved by classically sampling configurations in imaginary-time using standard lattice field theory. These configurations can then be passed to quantum processor for time-evolution. This method encounters a signal-to-noise problem which we characterize, and we demonstrate the application of standard lattice QCD methods to mitigate it.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗