skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Denoising of imaginary time response functions with Hankel projections
Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost—in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies. Published by the American Physical Society2024  more » « less
Award ID(s):
2310182
PAR ID:
10538870
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
Physical Review Research
Volume:
6
Issue:
3
ISSN:
2643-1564
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In the present work, a relativistic relation that connects the difference of interacting and noninteracting integrated two-particle correlation functions in finite volume to infinite volume scattering phase shift through an integral is derived. We show that the difference of integrated finite volume correlation functions converges rapidly to its infinite volume limit as the size of the periodic box is increased. The fast convergence of our proposed formalism is illustrated by analytic solutions of a contact interaction model, the perturbation theory calculation, and also the Monte Carlo simulation of a complex ϕ 4 lattice field theory model. Published by the American Physical Society2024 
    more » « less
  2. The formalism developed in the preceding papers that connects integrated correlation function of a trapped two-particle system to infinite volume scattering phase shift is further extended to coupled-channel systems in the present work. Using a trapped nonrelativistic two-channel system as an example, a new relation is derived that retains the same structure as in the single channel, and has explicit dependence on the phase shifts in both channels but not on the inelasticity. The relation is illustrated by a exactly solvable coupled-channel quantum mechanical model with contact interactions. It is further validated by path integral Monte Carlo simulation of a quasi-one-dimensional model that can admit general interaction potentials. In all cases, we found rapid convergence to the infinite volume limit as the trap size is increased, even at short times, making it potentially a good candidate to overcome signal-to-noise issues in Monte Carlo applications. Published by the American Physical Society2025 
    more » « less
  3. We design a quantum algorithm for ground state preparation in the early fault tolerant regime. As a Monte Carlo style quantum algorithm, our method features a Lindbladian where the target state is stationary. The construction of this Lindbladian is algorithmic and should not be seen as a specific approximation to some weakly coupled system-bath dynamics in nature. Our algorithm can be implemented using just one ancilla qubit and efficiently simulated on a quantum computer. It can prepare the ground state even when the initial state has zero overlap with the ground state, bypassing the most significant limitation of methods like quantum phase estimation. As a variant, we also propose a discrete-time algorithm, demonstrating even better efficiency and providing a near-optimal simulation cost depending on the desired evolution time and precision. Numerical simulations using Ising and Hubbard models demonstrate the efficacy and applicability of our method. Published by the American Physical Society2024 
    more » « less
  4. Midcircuit measurements (MCMs) are crucial ingredients in the development of fault-tolerant quantum computation. While there have been rapid experimental progresses in realizing MCMs, a systematic method for characterizing noisy MCMs is still under exploration. In this work, we develop a cycle benchmarking (CB)-type algorithm to characterize noisy MCMs. The key idea is to use a joint Fourier transform on the classical and quantum registers and then estimate parameters in the Fourier space, analogous to Pauli fidelities used in CB-type algorithms for characterizing the Pauli-noise channel of Clifford gates. Furthermore, we develop a theory of the noise learnability of MCMs, which determines what information can be learned about the noise model (in the presence of state preparation and terminating measurement noise) and what cannot, which shows that all learnable information can be learned using our algorithm. As an application, we show how to use the learned information to test the independence between measurement noise and state-preparation noise in an MCM. Finally, we conduct numerical simulations to illustrate the practical applicability of the algorithm. Similar to other CB-type algorithms, we expect the algorithm to provide a useful toolkit that is of experimental interest. Published by the American Physical Society2025 
    more » « less
  5. Dense spin ensembles in solids present a natural platform for studying quantum many-body dynamics. Multiple-pulse coherent control can be used to manipulate the magnetic dipolar interaction between the spins to engineer their dynamics. Here, we investigate the performance of a series of well-known pulse sequences that aim to suppress interspin dipolar couplings. We use a combination of numerical simulations and solid-state nuclear magnetic resonance experiments on adamantane to evaluate and compare sequence performance. We study the role of sequence parameters like interpulse delays and resonance offsets. Disagreements between experiments and theory are typically explained by the presence of control errors and experimental nonidealities. The simulations allow us to explore the influence of factors such as finite pulse widths, rotation errors, and phase transient errors. We also investigate the role of local disorder and establish that it is, perhaps unsurprisingly, a distinguishing factor in the decoupling efficiency of spectroscopic sequences (that preserve Hamiltonian terms proportional to S z ) and time-suspension sequences (which refocus all terms in the internal Hamiltonian). We discuss our findings in the context of previously known analytical results from average Hamiltonian theory. Finally, we explore the ability of time-suspension sequences to protect multispin correlations in the system. Published by the American Physical Society2025 
    more » « less