Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonequilibrium lattice models”

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.

49 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Learning stochastic dynamics and predicting emergent behavior using transformers

We show that a neural network originally designed for language processing can learn the dynamical rules of a stochastic system by observation of a single dynamical trajectory of the system, and can accurately predict its emergent behavior under conditions not observed during training. We consider a lattice model of active matter undergoing continuous-time Monte Carlo dynamics, simulated at a density at which its steady state comprises small, dispersed clusters. We train a neural network called a transformer on a single trajectory of the model. The transformer, which we show has the capacity to represent dynamical rules that are numerous and nonlocal, learns that the dynamics of this model consists of a small number of processes. Forward-propagated trajectories of the trained transformer, at densities not encountered during training, exhibit motility-induced phase separation and so predict the existence of a nonequilibrium phase transition. Transformers have the flexibility to learn dynamical rules from observation without explicit enumeration of rates or coarse-graining of configuration space, and so the procedure used here can be applied to a wide range of physical systems, including those with large and complex dynamical generators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Manipulating Growth and Propagation of Correlations in Dipolar Multilayers: From Pair Production to Bosonic Kitaev Models

Here, we study the nonequilibrium dynamics of dipoles confined in multiple stacked two-dimensional layers realizing a long-range interacting quantum spin 1/2 XXX model. We demonstrate that strong in-plane interactions can protect a manifold of collective layer dynamics. This then allows us to map the many-body spin dynamics to bosonic models. In a bilayer configuration we show how to engineer the paradigmatic two-mode squeezing Hamiltonian known from quantum optics, resulting in exponential production of entangled pairs and generation of metrologically useful entanglement from initially prepared product states. In multilayer configurations we engineer a bosonic variant of the Kitaev model displaying chiral propagation along the layer direction. Our study illustrates how the control over interactions, lattice geometry, and state preparation in interacting dipolar systems uniquely afforded by AMO platforms such as Rydberg and magnetic atoms, polar molecules, or trapped ions allows for the control over the temporal and spatial propagation of correlations for applications in quantum sensing and quantum simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The role of nonequilibrium LO phonons, Pauli exclusion, and intervalley pathways on the relaxation of hot carriers in InGaAs/InGaAsP multi-quantum-wells

Under continuous-wave laser excitation in a lattice-matched In 0.53 Ga 0.47 As/In 0.8 Ga 0.2 As 0.44 P 0.56 multi-quantum-well (MQW) structure, the carrier temperature extracted from photoluminescence rises faster for 405 nm compared with 980 nm excitation, as the injected carrier density increases. Ensemble Monte Carlo simulation of the carrier dynamics in the MQW system shows that this carrier temperature rise is dominated by nonequilibrium LO phonon effects, with the Pauli exclusion having a significant effect at high carrier densities. Further, we find a significant fraction of carriers reside in the satellite L-valleys for 405 nm excitation due to strong intervalley transfer, leading to a cooler steady-state electron temperature in the central valley compared with the case when intervalley transfer is excluded from the model. Good agreement between experiment and simulation has been shown, and detailed analysis has been presented. This study expands our knowledge of the dynamics of the hot carrier population in semiconductors, which can be applied to further limit energy loss in solar cells.

42 ENGINEERING↗
DOE OSTI2025

Unraveling the transformation pathway of the 𝛽 to 𝛾 phase transition in Ga 2 ⁢O 3 from atomistic simulations

Defect spinel 𝛾−Ga 2 ⁢O 3 is the least stable polymorph of Ga 2 ⁢O 3 , so its frequent appearance as a structural defect within or on the surface of monoclinic 𝛽−Ga 2 ⁢O 3 remains a mystery. Through first-principles calculations, we explore potential pathways for the phase transition from 𝛽−Ga 2⁢ O 3 to 𝛾−Ga 2 ⁢O 3 , and examine two key driving forces: tensile strain and Ga deficiency. When configurational entropy contributions to phase energies are included, the 𝛾 phase becomes energetically competitive with the 𝛽 phase, with the free energy difference between these phases diminishing even further under Ga-deficient conditions. Notably, a stability crossover occurs at room temperature at high vacancy concentrations ([V$^{3−}_{Ga}$]>3%) . A simple model 𝛽 → 𝛾 transformation pathway is identified, comprising two primary reactions, that enables the formation of the 𝛾 phase via simultaneous migration of Ga atoms from tetrahedral lattice sites to octahedral interstitial positions. The transformation barriers are prohibitively large in pristine Ga 2 ⁢O 3 , but can be substantially reduced by: (1) the presence of Ga vacancies, (2) elongational strains along the crystallographic 𝑎-axis, and (3) when volumetric relaxations are possible during transformation. These results elucidate prior experimental observations, where 𝛾−Ga 2⁢ O 3 is seen on damaged surfaces or in highly 𝑛-type 𝛽−Ga 2⁢ O 3 environments, which support Ga deficiency and mechanical strain. The insights into the driving forces and mechanisms of 𝛾−Ga 2⁢ O 3 formation enhance understanding of how localized strain and nonequilibrium defect concentrations may facilitate its formation from the 𝛽 phase.

Defects↗
DOE OSTI2020

Driven-dissipative quantum mechanics on a lattice: Simulating a fermionic reservoir on a quantum computer

The driven-dissipative many-body problem remains one of the most challenging unsolved problems in quantum mechanics. The advent of quantum computers may provide a unique platform for efficiently simulating such driven-dissipative systems. But, there are many choices for how one can engineer the reservoir. One can simply employ ancilla qubits to act as a reservoir and then digitally simulate them via algorithmic cooling. A more attractive approach, which allows one to simulate an infinite reservoir, is to integrate out the bath degrees of freedom and describe the driven-dissipative system via a master equation, that can also be simulated on a quantum computer. In this work, we consider the particular case of noninteracting electrons on a lattice driven by an electric field and coupled to a fermionic thermostat. Then, we provide two different quantum circuits: the first one reconstructs the full dynamics of the system using Trotter steps, while the second one dissipatively prepares the final nonequilibrium steady state in a single step. We run both circuits on the IBM quantum experience. For circuit (i), we achieved up to five Trotter steps. When partial resets become available on quantum computers, we expect that the maximum simulation time can be significantly increased. Lastly, the methods developed here suggest generalizations that can be applied to simulating interacting driven-dissipative systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Theoretical Description of Pump-Probe Experiments in Charge-Density-Wave Materials out to Long Times

We describe coupled nonequilibrium electron-phonon systems semiclassically—Ehrenfest dynamics for the phonons and quantum mechanics for the electrons—using a classical Monte Carlo approach that determines the nonequilibrium response to a large pump field. The semiclassical approach is expected to be accurate, because the phonons are excited to average energies much higher than the phonon frequency, eliminating the need for a quantum description. The numerical efficiency of this method allows us to perform a self-consistent time evolution out to very long times (tens of picoseconds), enabling us to model pump-probe experiments of a charge-density-wave (CDW) material. Our system is a half-filled, one-dimensional (1D) Holstein chain that exhibits CDW ordering due to a Peierls transition. The chain is subjected to a time-dependent electromagnetic pump field that excites it out of equilibrium, and then a second probe pulse is applied after a time delay. By evolving the system to long times, we capture the complete process of lattice excitation and subsequent relaxation to a new equilibrium, due to an exchange of energy between the electrons and the lattice, leading to lattice relaxation at finite temperatures. We employ an indirect (impulsive) driving mechanism of the lattice by the pump pulse due to the direct driving of the electrons. We identify two driving regimes, where the pump can either cause small perturbations or completely invert the initial CDW order. Our work successfully describes the ringing of the amplitude mode in CDW systems that has long been seen in experiment but never successfully explained by microscopic theory. We also describe the fluence-dependent crossover that inverts the CDW order parameter and changes the phonon dynamics. Finally, we illustrate how this method can examine a number of different types of experiments including photoemission, x-ray diffraction, and two-dimensional (2D) spectroscopy. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

Magnetism in the two-dimensional dipolar XY model

Here, motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [C. Chen et al., Nature (London) 616, 691 (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U⁡(1) symmetric setting. Our simulations perform extensive thermometry in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N≈0.5, close to the one at the critical temperature, S c /N=0.585⁢(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the nonequilibrium dynamics of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simulating cosmological evolution by quantum quench of an atomic Bose-Einstein condensate

In cosmological evolution, it is the homogeneous scalar field (inflaton) that drives the universe to expand isotropically and to generate standard model particles. However, to simulate cosmology, atomic gas research has focused on the dynamics of Bose-Einstein condensates (BEC) with continuously applied forces. In this paper we argue that a complementary approach needs also to be pursued; we thus consider the analog BEC experiments in a nondriven and naturally closed atomic system. We implement this using a BEC in an optical lattice which, after a quench, freely transitions from an unstable to a stable state. Furthermore, this dynamical evolution displays the counterpart “preheating,” “reheating,” and “thermalization” phases of cosmology. Importantly, our studies of these analog processes yield tractable analytic models. Additionally, of great utility to the cold atom community, such understanding elucidates the dynamics of nonadiabatic condensate preparation. Indeed, the dynamical processes discussed here are generic and in future cold atom reequilibration experiments it will be important to observe both the preheating stage, corresponding to a fragmented condensate, and the reheating stage, corresponding to a particle cloud.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Tunable fragile topology in Floquet systems

We extend the notion of fragile topology to periodically driven systems. Here we demonstrate driving-induced fragile topology in two different models, namely, the Floquet honeycomb model and the Floquet π-flux square-lattice model. In both cases, we discover a rich phase diagram that includes Floquet fragile topological phases protected by crystalline rotation or mirror symmetries, Floquet Chern insulators, and trivial atomic phases, with distinct boundary features. Remarkably, the transitions between different phases can be feasibly achieved by simply tuning the driving amplitudes, which is a unique feature of driving-enabled topological phenomena. Moreover, corner-localized fractional charges are identified as a “smoking-gun” signal of fragile topology in our systems. Our work paves the way for studying and realizing fragile topology in Floquet systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Complex paths around the sign problem

The Monte Carlo evaluation of path integrals is one of a few general purpose methods to approach strongly coupled systems. It is used in all branches of physics, from QCD and nuclear physics to the correlated electron systems. However, many systems of great importance (dense matter inside neutron stars, the repulsive Hubbard model away from half filling, and dynamical and nonequilibrium observables) are not amenable to the Monte Carlo method as it currently stands due to the so-called sign problem. Here, a new set of ideas recently developed to tackle the sign problem based on the complexification of field space and the Picard-Lefshetz theory accompanying it is reviewed. The mathematical ideas underpinning this approach, as well as the algorithms developed thus far, are described together with nontrivial examples where the method has already been proved successful. Directions of future work, including the burgeoning use of machine learning techniques, are delineated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Probing configurational disorder in ZnGeN 2 using cluster-based Monte Carlo

ZnGeN 2 is sought as a semiconductor with comparable lattice constant to GaN and tunable band gap for integration in optoelectronic devices. Configurational disorder on the cation sublattice of ZnGeN 2 can strongly modify the electronic structure compared to the ordered material, and both ordered and disordered forms of ZnGeN 2 are candidates for light-emitting diodes and other emitters. The nonisovalent character of the disordered species (Zn 2+ and Ge 4+ ) subjects the cation ordering to strong short-range order effects. To model these effects, we use Monte Carlo (MC) simulations utilizing a cluster expansion to approximate formation enthalpy. Representative disordered configurations in 1024-atom supercells are relaxed in density functional theory calculations. From the MC structures, we extract a short-range order parameter (the N-cation coordination motif), and two long-range order parameters (Bragg-Williams and stretching parameters), and examine their correlations. We perform a thermodynamic integration to determine the mixing entropy and free energy. ZnGeN 2 exhibits a first-order phase transition with pronounced discontinuities in enthalpy and entropy, as well as in the structural order parameters. We discuss the relationship between the effective temperature used in the MC simulation and the growth temperatures in experiment in relation to the crossover from the nonequilibrium to the equilibrium growth regime. This work expands on current models of site disorder in ZnGeN 2 and provides atomic structure models with a systematic variation of the degree of cation disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗