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 OSTI2017

A multilayer multiconfiguration time-dependent Hartree simulation of the reaction-coordinate spin-boson model employing an interaction picture

The multilayer multiconfiguration time-dependent Hartree method is applied in an interaction picture to simulate dynamics of the spin-boson model in the reaction-coordinate representation. The use of the interaction picture allows a more effective description of correlation effects, especially when the coupling strength between the reaction coordinate and the bath is very strong. Examples show that in most physical regimes the efficiency is improved significantly, in some cases up to several orders of magnitude. This opens up new avenues for studying quantum dynamical problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2017

Emptying Dirac valleys in bismuth using high magnetic fields

The Fermi surface of elemental bismuth consists of three small rotationally equivalent electron pockets, offering a valley degree of freedom to charge carriers. A relatively small magnetic field can confine electrons to their lowest Landau level. This is the quantum limit attained in other dilute metals upon application of sufficiently strong magnetic field. Here in this paper we report on the observation of another threshold magnetic field never encountered before in any other solid. Above this field, B empty , one or two valleys become totally empty. Drying up a Fermi sea by magnetic field in the Brillouin zone leads to a manyfold enhancement in electric conductance. We trace the origin of the large drop in magnetoresistance across B empty to transfer of carriers between valleys with highly anisotropic mobilities. The non-interacting picture of electrons with field-dependent mobility explains most results but the Coulomb interaction may play a role in shaping the fine details.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Communication: A difference density picture for the self-consistent field ansatz

We formulate self-consistent field (SCF) theory in terms of an interaction picture where the working variable is the difference density matrix between the true system and a corresponding superposition of atomic densities. As the difference density matrix directly represents the electronic deformations inherent in chemical bonding, this “difference self-consistent field (dSCF)” picture provides a number of significant conceptual and computational advantages. We show that this allows for a stable and efficient dSCF iterative procedure with wholly single-precision Coulomb and exchange matrix builds. We also show that the dSCF iterative procedure can be performed with aggressive screening of the pair space. These approximations are tested and found to be accurate for systems with up to 1860 atoms and >10 000 basis functions, providing for immediate overall speedups of up to 70% in the heavily optimized TERACHEM SCF implementation.

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

Resolution-of-identity stochastic time-dependent configuration interaction for dissipative electron dynamics in strong fields

In this contribution, we introduce a method for simulating dissipative, ultrafast many-electron dynamics in intense laser fields. The method is based on the norm-conserving stochastic unraveling of the dissipative Liouville-von Neumann equation in its Lindblad form. The N-electron wave functions sampling the density matrix are represented in the basis of singly excited configuration state functions. The interaction with an external laser field is treated variationally and the response of the electronic density is included to all orders in this basis. The coupling to an external environment is included via relaxation operators inducing transition between the configuration state functions. Single electron ionization is represented by irreversible transition operators from the ionizing states to an auxiliary continuum state. The method finds its efficiency in the representation of the operators in the interaction picture, where the resolution-of-identity is used to reduce the size of the Hamiltonian eigenstate basis. The zeroth-order eigenstates can be obtained either at the configuration interaction singles level or from a time-dependent density functional theory reference calculation. The latter offers an alternative to explicitly time-dependent density functional theory which has the advantage of remaining strictly valid for strong field excitations while improving the description of the correlation as compared to configuration interaction singles. The method is tested on a well-characterized toy system, the excitation of the low-lying charge transfer state in LiCN.

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

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗
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 OSTI2019

On the efficiency of different numerical methods for the calculation of intrapulse Raman scattering of optical solitons

In this paper we compare the performance of different numerical methods for the calculation of the asymptotic evolution of soliton self-frequency shift in the presence of intrapulse Raman scattering (IRS) in optical fibers. First we have calculated the order of global accuracy for the fundamental soliton and the second-order bound state of the unperturbed nonlinear Schrödinger equation for the following numerical methods: the simple split step (SS) method, the full SS method, the reduced SS method, the reduced SS method with fourth order Runge–Kutta (RK4), the Blow–Wood method, the fourth order Runge–Kutta in the interaction picture (RK4IP) method and the Agrawal SS method with one and two iterations. We have shown that the asymptotic evolution of soliton self-frequency shift in the presence of IRS in optical fiber can be best described by the Agrawal SS method (compared to the Blow–Wood method and the RK4IP method). The obtained numerical results for the soliton position and frequency are in agreement with the predictions for these parameters according to the perturbation theory. We have shown that in the presence of IRS the fundamental soliton quickly develops an oscillating tail on its left. The generated tail hardly influences the soliton propagation at large distances. At large distances the form of the soliton gradually becomes asymmetric. The increase of the frequency resolution leads to a notable increase of the maximum propagation distance of the numerical soliton under the influence of IRS.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Salt Solutions in Carbon Nanotubes: The Role of Cation-π Interactions

Understanding the structure of aqueous electrolytes at interfaces is essential for predicting and optimizing device performance for a wide variety of emerging energy and environmental technologies. In this work, we investigate the structure of two com- mon salt solutions, NaCl and KCl, at a hydrophobic interface within narrow carbon nanotubes (CNTs). Using a combination of first-principles and classical molecular dynamics simulations in conjunction with molecular orbital analysis, we find that the solvation structure of the cations in the CNTs can deviate substantially from the conventional weakly interacting hydrophobic picture. Instead, interactions between solvated ions and π-orbitals of the CNTs are found to play a critically important role. Specifically, the ion solvation structure is ultimately determined by a complex interplay between cation- interactions and the intrinsic flexibility of the solvation shell. In the case of K + , these effects result in an unusually strong propensity to partially desolvate and reside closer to the carbon wall than both Na + and Cl - , in sharp contrast to the known ion ordering at the water-vapor interface.

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

Partition-DFT on the water dimer

As is well known, the ground-state symmetry group of the water dimer switches from its equilibrium C s -character to C 2h -character as the distance between the two oxygen atoms of the dimer decreases below R O–O ~2.5 Å. For a range of R O–O between 1 and 5 Å, and for both symmetries, we apply Partition Density Functional Theory (PDFT) to find the unique monomer densities that sum to the correct dimer densities while minimizing the sum of the monomer energies. We calculate the work involved in deforming the isolated monomer densities and find that it is slightly larger for the Cs geometry for all R O–O . Finally, we discuss how the PDFT densities and the corresponding partition potentials support the orbital-interaction picture of hydrogen-bond formation.

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

Versatile laser-free trapped-ion entangling gates

We present a general theory for laser-free entangling gates with trapped-ion hyperfine qubits, using either static or oscillating magnetic-field gradients combined with a pair of uniform microwave fields symmetrically detuned about the qubit frequency. By transforming into a 'bichromatic' interaction picture, we show that either ${\hat{\sigma }}_{\phi }\otimes {\hat{\sigma }}_{\phi }$ or ${\hat{\sigma }}_{z}\otimes {\hat{\sigma }}_{z}$ geometric phase gates can be performed. The gate basis is determined by selecting the microwave detuning. The driving parameters can be tuned to provide intrinsic dynamical decoupling from qubit frequency fluctuations. The ${\hat{\sigma }}_{z}\otimes {\hat{\sigma }}_{z}$ gates can be implemented in a novel manner which eases experimental constraints. Here, we present numerical simulations of gate fidelities assuming realistic parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

Optimal control of coupled quantum systems based on the first-order Magnus expansion: Application to multiple dipole-dipole-coupled molecular rotors

This paper presents a method for performing approximate optimal control simulations for quantum systems with multiple coupled degrees of freedom. In this work, the time evolution is simulated using the first-order Magnus expansion in the interaction picture, where the couplings between different degrees of freedom are treated as the perturbation. A numerical implementation procedure is presented that leverages upon pairwise couplings and the separability of the zeroth-order time evolution operator to achieve a reduced computational cost, which is analyzed with respect to the number of degrees of freedom. The formulation is compatible with gradient-free methods to optimize the control field, and a stochastic hill climbing algorithm is adopted for this purpose. As illustrations, optimal control simulations are performed for systems of two and three dipole-dipole-coupled molecular rotors under the influence of a control field. For the two-rotor system, the field is optimized to achieve either orientation or entanglement objectives. For the three-rotor system, the field is optimized either to orient all three rotors in the same direction or to orient one rotor in a particular direction while the other two rotors point in the opposite direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Proton decays in 16 Ne and 18 Mg and isospin-symmetry breaking in carbon isotopes and isotones

We report that proton-rich nuclei possess unique properties in the nuclear chart. Due to the presence of both continuum coupling and Coulomb interaction, phenomena such as halos, Thomas-Ehrman shift, and proton emissions can occur. Relevant experimental data are difficult to obtain, so that theoretical calculations are needed to understand nuclei at drip lines and to guide experimentalists. In particular, the 16 Ne and 18 Mg isotopes are supposed to be one-proton and/or two-proton emitting nuclei, but associated experimental data are either incomplete or even unavailable. Consequently, we performed Gamow shell model calculations of carbon isotones bearing A = 15 - 18 . Isospin-symmetry breaking occurring in carbon isotones and isotopes is also discussed. It is hereby shown that the mixed effects of continuum coupling and Coulomb interaction at drip lines generate complex patterns in isospin multiplets. Added to that, it is possible to determine the one-proton and two-proton widths of 16 Ne and 18 Mg . Obtained decay patterns are in agreement with those obtained in previous experimental and theoretical works. Moreover, to our knowledge, this is the first theoretical calculation of binding energy and partial decay widths of 18 Mg in a configuration interaction picture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗