Search NASA⌕ Search

SEARCH · Search NASA

Results for “interaction picture”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Piecewise interaction picture density matrix quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

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

Hybridized Methods for Quantum Simulation in the Interaction Picture

Conventional methods of quantum simulation involve trade-offs that limit their applicability to specific contexts where their use is optimal. In particular, the interaction picture simulation has been found to provide substantial asymptotic advantages for some Hamiltonians, but incurs prohibitive constant factors and is incompatible with methods like qubitization. We provide a framework that allows different simulation methods to be hybridized and thereby improve performance for interaction picture simulations over known algorithms. These approaches show asymptotic improvements over the individual methods that comprise them and further make interaction picture simulation methods practical in the near term. Physical applications of these hybridized methods yield a gate complexity scaling as log 2 ⁡ Λ in the electric cutoff Λ for the Schwinger Model and independent of the electron density for collective neutrino oscillations, outperforming the scaling for all current algorithms with these parameters. For the general problem of Hamiltonian simulation subject to dynamical constraints, these methods yield a query complexity independent of the penalty parameter λ used to impose an energy cost on time-evolution into an unphysical subspace.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Efficient simulation of open quantum systems coupled to a reservoir through multiple channels

It is challenging to simulate open quantum systems that are connected to a reservoir through multiple channels. For example, vibrations may induce fluctuations in both energy gaps and electronic couplings, which represent two independent channels of system–bath couplings. Systems of this kind are ubiquitous in the processes of excited state radiationless decay. Combined with density matrix renormalization group (DMRG) and matrix product states (MPS) methods, we develop an interaction-picture chain mapping strategy for vibrational reservoirs to simulate the dynamics of these open systems, resulting in time-dependent spatially local system–bath couplings in the chain-mapped Hamiltonian. This transformation causes the entanglement generated by the system–bath interactions to be restricted within a narrow frequency window of vibrational modes, enabling efficient DMRG/MPS dynamical simulations. We demonstrate the utility of this approach by simulating singlet fission dynamics using a generalized spin-boson Hamiltonian with both diagonal and off-diagonal system–bath couplings. This approach generalizes an earlier interaction-picture chain mapping scheme, allowing for efficient and exact simulation of systems with multi-channel system–bath couplings using matrix product states, which may further our understanding of nonlocal exciton–phonon couplings in exciton transport and the non-Condon effect in energy and electron transfer.

Chemistry↗
DOE OSTI2024

The Gauge Picture of Quantum Dynamics

Although local Hamiltonians exhibit local time dynamics, this locality is not explicit in the Schrödinger picture in the sense that the wavefunction amplitudes do not obey a local equation of motion. We show that geometric locality can be achieved explicitly in the equations of motion by "gauging" the global unitary invariance of quantum mechanics into a local gauge invariance. That is, expectation values $\langle \psi|A|\psi \rangle$ are invariant under a global unitary transformation acting on the wavefunction $|\psi\rangle \to U |\psi\rangle$ and operators $A \to U A U^\dagger$, and we show that it is possible to gauge this global invariance into a local gauge invariance. To do this, we replace the wavefunction with a collection of local wavefunctions $|\psi_J\rangle$, one for each patch of space $J$. The collection of spatial patches is chosen to cover the space; e.g. we could choose the patches to be single qubits or nearest-neighbor sites on a lattice. Local wavefunctions associated with neighboring pairs of spatial patches $I$ and $J$ are related to each other by dynamical unitary transformations $U_{IJ}$. The local wavefunctions are local in the sense that their dynamics are local. That is, the equations of motion for the local wavefunctions $|\psi_J\rangle$ and connections $U_{IJ}$ are explicitly local in space and only depend on nearby Hamiltonian terms. (The local wavefunctions are many-body wavefunctions and have the same Hilbert space dimension as the usual wavefunction.) We call this picture of quantum dynamics the gauge picture since it exhibits a local gauge invariance. The local dynamics of a single spatial patch is related to the interaction picture, where the interaction Hamiltonian consists of only nearby Hamiltonian terms. We can also generalize the explicit locality to include locality in local charge and energy densities.

Physics↗
DOE OSTI2026

Discrete Superconvergence Analysis for Quantum Magnus Algorithms of Unbounded Hamiltonian Simulation

Motivated by various applications, unbounded Hamiltonian simulation has recently garnered great attention. Quantum Magnus algorithms, designed to achieve commutator scaling for time-dependent Hamiltonian simulation, have been found to be particularly efficient for such applications. When applied to unbounded Hamiltonian simulation in the interaction picture, they exhibit an unexpected superconvergence phenomenon. However, existing proofs are limited to the spatially continuous setting and do not extend to discrete spatial discretizations. Here, in this work, we provide the first superconvergence estimate in the fully discrete setting with a finite number of spatial discretization points N, and show that it holds with an error constant uniform in N. The proof is based on the two-parameter symbol class, which, to our knowledge, is applied for the first time in algorithm analysis. The key idea is to establish a semiclassical framework by identifying two parameters through the discretization number and the time step size rescaled by the operator norm, such that the semiclassical uniformity guarantees the uniformity of both. This approach may have broader applications in numerical analysis beyond the specific context of this work.

Borns-Weil, Yonah [University of California, Berke↗
DOE OSTI2025

Quantum Hamiltonian algorithms for maximum independent sets

ABSTRACT We compare two quantum Hamiltonian algorithms that address the maximum independent set problem: one based on the emergent non-Abelian gauge matrix in adiabatic evolution of an energetically isolated manifold of states; the other based on designed application of single-qubit operations. We demonstrate that they are mathematically equivalent in the sense that one is the other’s interaction picture. Despite their mathematical equivalence, our numerical simulations show significant differences between them in performance, which is explained analytically. Intriguingly, this equivalence unveils that the PXP model, recently prominent in quantum dynamics research, can be viewed as quantum diffusion over the median graph of all independent sets governed by the non-Abelian gauge matrix.

Science & Technology - Other Topics↗
DOE OSTI2022

Time-dependent Hamiltonian Simulation of Highly Oscillatory Dynamics and Superconvergence for Schrödinger Equation

We propose a simple quantum algorithm for simulating highly oscillatory quantum dynamics, which does not require complicated quantum control logic for handling time-ordering operators. To our knowledge, this is the first quantum algorithm that is both insensitive to the rapid changes of the time-dependent Hamiltonian and exhibits commutator scaling. Our method can be used for efficient Hamiltonian simulation in the interaction picture. In particular, we demonstrate that for the simulation of the Schrödinger equation, our method exhibits superconvergence and achieves a surprising second order convergence rate, of which the proof rests on a careful application of pseudo-differential calculus. Numerical results verify the effectiveness and the superconvergence property of our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Spin dynamics in the skyrmion-host lacunar spinel GaV 4 S 8

Here in the lacunar spinel GaV 4 S 8 , the interplay of spin, charge, and orbital degrees of freedom produces a rich phase diagram that includes an unusual Néel-type skyrmion phase composed of molecular spins. To provide insight into the interactions underlying this complex phase diagram, we study the spin excitations in GaV 4 S 8 through inelastic neutron scattering measurements on polycrystalline and single crystal samples. Using linear spin wave theory, we describe the spin wave excitations using a model where V 4 clusters decorate an fcc lattice. The effective cluster model includes a ferromagnetic interaction and a weaker antisymmetric Dzyaloshinskii-Moriya interaction between the neighboring molecular spins. Our work clarifies the spin interactions in GaV 4 S 8 and supports the picture of interacting molecular clusters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Statistical description of coalescing magnetic islands via magnetic reconnection

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.

Physics↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Computational projects with the Landau–Zener problem in the quantum mechanics classroom

The Landau–Zener problem, where a minimum energy separation is passed with constant rate in a two-state quantum-mechanical system, is an excellent model quantum system for a computational project. It requires a low-level computational effort, but has a number of complex numerical and algorithmic issues that can be resolved through dedicated work. It can be used to teach computational concepts, such as accuracy, discretization, and extrapolation, and it reinforces quantum concepts of time-evolution via a time-ordered product and of extrapolation to infinite time via time-dependent perturbation theory. Additionally, we discuss the concept of compression algorithms, which are employed in many advanced quantum computing strategies, and easy to illustrate with the Landau–Zener problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Spectral Signatures of a Negative Polaron in a Doped Polymer Semiconductor: Energy Levels and Hubbard U Interactions

The modern picture of negative charge carriers on conjugated polymers invokes the formation of a singly occupied (spin-up/spin-down) level within the polymer gap and a corresponding unoccupied level above the polymer conduction band edge. The energy splitting between these sublevels is related to on-site Coulomb interactions between electrons, commonly termed Hubbard U. However, spectral evidence for both sublevels and experimental access to the U value is still missing. Here, in this work, we provide evidence by n-doping the polymer P(NDI 2 OD-T 2 ) with [RhCp*Cp] 2 , [N-DMBI] 2 , and cesium. Changes in the electronic structure after doping are studied with ultraviolet photoelectron and low-energy inverse photoemission spectroscopies (UPS, LEIPES). UPS data show an additional density of states (DOS) in the former empty polymer gap while LEIPES data show an additional DOS above the conduction band edge. These DOS are assigned to the singly occupied and unoccupied sublevels, allowing determination of a U value of ~1 eV.

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

Ligand-modified nanoparticle surfaces influence CO electroreduction selectivity

Improving the kinetics and selectivity of CO 2 /CO electroreduction to valuable multi-carbon products is a challenge for science and is a requirement for practical relevance. Here we develop a thiol-modified surface ligand strategy that promotes electrochemical CO-to-acetate. We explore a picture wherein nucleophilic interaction between the lone pairs of sulfur and the empty orbitals of reaction intermediates contributes to making the acetate pathway more energetically accessible. Density functional theory calculations and Raman spectroscopy suggest a mechanism where the nucleophilic interaction increases the sp 2 hybridization of CO (ad) , facilitating the rate-determining step, CO* to (CHO)*. We find that the ligands stabilize the (HOOC–CH 2 )* intermediate, a key intermediate in the acetate pathway. In-situ Raman spectroscopy shows shifts in C–O, Cu–C, and C–S vibrational frequencies that agree with a picture of surface ligand-intermediate interactions. A Faradaic efficiency of 70% is obtained on optimized thiol-capped Cu catalysts, with onset potentials 100 mV lower than in the case of reference Cu catalysts.

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

Universal transport in periodically driven systems without long-lived quasiparticles

An intriguing regime of universal charge transport at high entropy density has been proposed for periodically driven interacting one-dimensional systems with Bloch bands separated by a large single-particle band gap. For weak interactions, a simple picture based on well-defined Floquet quasiparticles suggests that the system should host a quasisteady state current that depends only on the populations of the system's Floquet-Bloch bands and their associated quasienergy winding numbers. Here we show that such topological transport persists into the strongly interacting regime where the single-particle lifetime becomes shorter than the drive period. Analytically, we show that the value of the current is insensitive to interaction-induced band renormalizations and lifetime broadening when certain conditions are met by the system's nonequilibrium distribution function. We show that these conditions correspond to a quasisteady state. We support these predictions through numerical simulation of a system of strongly interacting fermions in a periodically modulated chain of Sachdev-Ye-Kitaev dots. Our paper establishes universal transport at high entropy density as a robust far from equilibrium topological phenomenon, which can be readily realized with cold atoms in optical lattices. Published by the American Physical Society 2024

Esin, Iliya (ORCID:0000000329590617)↗
DOE OSTI2025

Probing the influence of ion-pairing on ligand-field excited-state dynamics

Exploration of the photophysical and photochemical properties of transition metal complexes has driven ground-breaking advancements in solar energy conversion technologies, including photoredox catalysis. While significant research has been devoted to understanding excited state properties of second- and third-row transition metal complexes, earth-abundant first-row metal complexes have received comparatively little attention in this context until very recently. In particular, the role of ion-pairing – which has been identified as a potentially significant factor for Ir(III)-based photosensitizers – has not been examined with regard to its influence on the ligand-field excited states that dominate much of first-row photophysics. A key challenge in studying ion-pair interactions lies in quantifying the extent and nature of ion-pairing, particularly in non-aqueous media where the vast majority of photophysical studies are performed. Cobalt(III) polypyridyl complexes provide an attractive platform to address such questions due to their demonstrated potential for applications in photoredox catalysis involving ligand-field excited states. In the present study, we prepared a cobalt(III) polypyridyl complex, [Co(4,4′-OMebpy) 3 ](BAr F 4 ) 3 (where 4,4′-OMebpy is 4,4′-dimethoxy-2,2′-bipyridine and BAr F 4 is tetrakis(3,5-bis(trifluoromethyl)-phenyl)borate) to probe ion-pairing in non-aqueous solutions. Specifically, analysis of data acquired from both variable-temperature diffusion ordered spectroscopy (DOSY) NMR and 1-D rotating-frame nuclear Overhauser effect (ROE) experiments allowed us to identify and differentiate between solvent-separated ion pairs in high-dielectric media and contact ion pairs in a low-dielectric solvent. Time-resolved absorption spectroscopy was then used to measure ground-state recovery dynamics under these varying conditions of ion-pairing, the results of which revealed an increase in excited-state lifetime for contact ion-pairs that we suggest arises from a reduction in outer-sphere reorganization energy relative to conditions which favored solvent-separated ion pairs. We believe this study demonstrates that one can leverage broadly available NMR-based methods to understand ion-pairing in non-aqueous solutions, which in turn can provide a microscopic picture of intermolecular interactions that can impact the photophysical properties of transition metal-based chromophores.

Ghosh, Atanu [Michigan State University, East Lans↗
DOE OSTI2023

Thermodynamics of Minerals Stable near Earth's Surface

Project “Thermodynamics of Minerals Stable Near the Earth’s Surface” uses specialized calorimetric techniques to obtain new thermochemical data for solid phases of interest in the Earth’s “critical zone”, the near surface regions of the planet most directly affecting and affected by human activity. Major activity was in three areas: manganese oxide materials, aluminum clusters which are precursors to aluminum oxyhydroxide precipitates, and zeolites. All these classes of minerals are potential carriers of heavy metals and radionuclides in the environment, both within their crystal structures and as adsorbed layers or colloids. Their thermodynamic stability, as well as the sequence of stability of competing phases, governs the aqueous solubility of contaminants (and nutrients) in soil, groundwater, and sediments. The thermodynamics of these phases is strongly affected by particle size and crystallinity, since they often occur as very fine particles which do not coarsen in the near-surface environment. Many of these materials exhibit extensive polymorphism, multicomponent solid solution, and variable hydration states. Thus, one must define and systematize important thermodynamic variables, since one cannot make measurements on each and every possible phase, particle size, and composition. A major finding of the work on manganese oxides was labile water-induced oxidation-reduction reactions. These will be studied further by a combination of thermodynamic and synchrotron-based spectroscopic methods, the latter in collaboration with Peggy O’Day at UC Merced. The adsorption of other species, including simple organics, on such redox-active surfaces will also be investigated by a combination of thermochemical and spectroscopic studies. Studies of the energetics of water adsorption on mineral surfaces will be expanded to emphasize the adsorption of CO2 and small organic molecules, relevant to soil, groundwater, and the subsurface environment. Gas adsorption calorimetry, immersion calorimetry, and solution calorimetry will be combined with appropriate analytical and spectroscopic techniques to obtain a comprehensive picture of the interactions and their energetics. The evolution of aqueous aluminum clusters to form aluminum oxyhydroxide precipitates will be studied in detail by a combination of calorimetric and structural methods applied to aqueous samples aged at room temperature. The energetics of new classes of clusters will be explored. The project addresses several common themes relevant to the broad mission of the Department of Energy. Thermochemical data for minerals stable near the Earth’s surface, though still incomplete, are essential for understanding both natural and anthropogenic processes that affect our very existence on this planet. Such data contribute to the fundamental underpinning for mitigation and management strategies for the extraction and use of energy. Fundamental thermodynamic and mechanistic understanding is required to model processes ranging from heavy metal pollution to nuclear waste disposal to climate change.

58 GEOSCIENCES↗
DOE OSTI2024

Filtering higher-order datasets

Many complex systems often contain interactions between more than two nodes, known as higher-order interactions, which can change the structure of these systems in significant ways. Researchers often assume that all interactions paint a consistent picture of a higher-order dataset's structure. In contrast, the connection patterns of individuals or entities in empirical systems are often stratified by interaction size. Ignoring this fact can aggregate connection patterns that exist only at certain scales of interaction. To isolate these scale-dependent patterns, we present an approach for analyzing higher-order datasets by filtering interactions by their size. We apply this framework to several empirical datasets from three domains to demonstrate that data practitioners can gain valuable information from this approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

TMD factorisation for diffractive jets in photon-nucleus interactions

Abstract Using the colour dipole picture and the colour glass condensate effective theory, we study the diffractive production of two or three jets via coherent photon-nucleus interactions at high energy. We consider the hard regime where the photon virtuality and/or the transverse momenta of the produced jets are much larger than the saturation momentumQ s of the nuclear target. We show that, despite this hardness, the leading-twist contributions are controlled by relatively large parton configurations, with transverse sizesR~ 1/Q s , which undergo strong scattering and probe gluon saturation. We demonstrate that these leading-twist contributions admit transverse-momentum dependent (TMD) factorisation, in terms of quark and gluon diffractive TMD distribution functions, for which we obtain explicit expressions from first principles. We go beyond our previous work by evaluating the contributions involving the quark diffractive distributions and by establishing that their DGLAP evolution emerges via controlled calculations within the colour dipole picture. We find the same expression for the quark diffractive TMD in two different processes (semi-inclusive diffraction and the diffractive production of quark-gluon dijets), thus demonstrating its universality.

Physics↗