A Symmetry-Enabled Direct Quantum Protocol for Many-Body Green's Functions

arXiv:2509.03288 · quant-ph · Submitted 2025-09-03 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: I'm Kai.

Mira: And I'm Mira, here with Lev.

Kai: We're looking at the paper: "A Symmetry-Enabled Direct Quantum Protocol for Many-Body Green's Functions".

Mira: Here's the gist of the paper. The gist The protocol introduces a symmetry-enabled direct quantum protocol, tailored quench spectroscopy (TQS), for computing both the real and imaginary parts of two-point time correlators without controlled unitaries or ancilla qubits in the correlator-estimation circuit<ref:2509.03288#pg12>.

Kai: This method exploits parity symmetry to simplify readout, providing a practical route to estimating symmetry-resolved dynamical correlators on near-term and early fault-tolerant quantum hardware<ref:2509.03288#pg18>.

Mira: How it works The protocol relies on exploiting parity symmetry, satisfied by a broad class of Hamiltonians in condensed matter physics and quantum chemistry, including the Fermi–Hubbard and Heisenberg models<ref:2509.03288#pg17>.

Kai: The core idea is to introduce a tailored quench spectroscopy scheme that recovers both the real and imaginary parts of two-point time correlators, from which Green’s functions can be reconstructed via efficient classical signal processing<ref:2509.03288#pg18>.

Title and authors: Kai: So we're looking at this paper today, "A Symmetry-Enabled Direct Quantum Protocol for Many-Body Green's Functions." It sounds pretty technical, but the core idea is that they figured out a way to calculate those many-body Green's functions without needing any extra qubits or complex unitaries in the circuit itself.

Mira: Yeah, it’s about using parity symmetry to get both the real and imaginary parts of two-point time correlators directly, which are essential for studying strongly correlated quantum systems. The authors are Changhao Yi and Cunlu Zhou.

Lev: From an engineering standpoint, that's huge because it means we can do this using only native time evolution and measurements already available on current hardware platforms without needing fancy ancilla qubits or complex setups to get those correlators.

Kai: Right, so the focus is on making this practical for near-term and even early fault-tolerant quantum hardware. The paper says they exploit parity symmetry which is satisfied by models like the Fermi–Hubbard and Heisenberg models.

Mira: Exactly, it's about taking that symmetry and using it to simplify how we read out the data, which allows them to reconstruct those Green’s functions using classical signal processing routines.

Lev: It avoids the high sample complexity issues you see with some other methods like the direct Krylov space method when you scale up the problem size.

The paper's summary: Kai: So, what does this protocol actually do? The paper introduces a tailored quench spectroscopy scheme that recovers both the real and imaginary parts of two-point time correlators from which Green’s functions can be reconstructed via efficient classical signal analysis.

Mira: It uses specific quench operators to get the real part of the time correlator, and they show that for a ground state initial condition, you don't even need the thermal state preparation step.

Lev: The complexity analysis shows that to hit an additive error epsilon with constant success probability, you need shots per time point scaling as O(epsilon-two). That’s a concrete number we can work with on simulators or hardware.

Kai: It also says they developed a tailored symmetric quantum Gibbs sampler to prepare parity-resolved thermal states, which are defined as rhoS for the symmetric and rhoA for the antisymmetric state.

Mira: Those states allow them to tackle finite-temperature extensions, which is something you usually have to simulate separately, and they define those thermal states using a tailored quantum Gibbs sampler based on dissipative state preparation.

Lev: For real hardware, that means we're talking about simulating the Davies generator LD or LD,S under a Lindbladian framework to get those symmetric and antisymmetric states accurately.

The paper's improvements: Kai: One key improvement they highlight is how they handle error mitigation. They show that symmetry-based symmetrization can quadratically suppress parity-breaking state-preparation errors, reducing the leading error contribution from O(mu) to O(mu two) under specific conditions <ref:2509.03288#pg1>.

Mira: That's a big deal because it means the state preparation errors don't just stay as linear; they get suppressed by a square factor. They achieve this by applying a symmetrization channel S(rho) equals rho plus P rhoP2, which removes coherences between opposite parity sectors.

Lev: That quadratic suppression is something we need when we think about running these on real hardware where state preparation errors are notoriously tricky to control. It reduces the leading error contribution from O(mu) to O(mu two) under specific conditions <ref:2509.03288#pg1>.

Kai: They also extend this approach naturally to out-of-time-order correlators, or OTOCs, which is interesting because it means the same symmetry primitive works for those too.

Mira: For an OTOC(A, B, t), the experiment reduces to estimating modified quench functions derived from the OTOC structure using tailored quench operators UIm and URe. This shows that this measurement primitive extends naturally to out-of-time-ordered correlators.

Conclusion: Kai: So, wrapping up on this paper, the main point is they’ve created a way to compute many-body Green's functions using only native time evolution and straightforward measurements available on current hardware platforms by exploiting parity symmetry.

Mira: They show how you can get both the real and imaginary parts of two-point time correlators directly, which lets you reconstruct the Green's functions using classical signal processing routines like MUSIC.

Lev: The method is robust because it uses time-independent Hamiltonian evolution, avoiding the simulation of time-dependent dynamics, which reduces experimental complexity. But they also flag a limitation: they need to resolve the frequency content of the measured timedomain correlator using methods like MUltiple SIgnal Classification for spectral reconstruction.

Kai: And they say that their error analysis shows that with O(epsilon-two) scaling for shots per time point, and the quadratic suppression of state preparation errors, this is a practical route to estimating symmetry-resolved dynamical correlators on near-term hardware.

Mira: It opens up several exciting directions, including the measurement of higher-order and nonlinear response functions and exploring analogous ideas in non-Hermitian or effective open-system settings.

Lev: For real hardware, this is a promising first milestone because it's compatible with both digital platforms and analog simulators, making it a good starting point for experimental demonstration on small Hubbard or spin-chain instances.

Changhao Yi, Cunlu Zhou

Department of Physics, Shanghai University · Department of Computer Science, Universitພ de Sherbrooke · Institut quantique, Universitພ de Sherbrooke

quant-ph

Submitted: 2025-09-03

Updated: 2026-10-19

Comments: The new version includes added references (including the dataset), updated acknowledgments, and general text polishing

Code: https://github.com/CYI1995/TQS

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 78/100

The gist: The gist The protocol introduces a symmetry-enabled direct quantum protocol, tailored quench spectroscopy (TQS), for computing both the real and imaginary parts of two-point time correlators without

Key concepts

Parity Symmetry: This is a fundamental symmetry in many physical systems where a specific operator P leaves the Hamiltonian unchanged. The protocol exploits this property to simplify measurements, allowing researchers to extract information about different states (parity sectors) without needing extra qubits or complex operations.
Tailored Quench Spectroscopy (TQS): A specialized technique applied during a quantum quench—a sudden change in the system's Hamiltonian. TQS is designed specifically to recover both the real and imaginary components of two-point time correlators, which are essential for reconstructing the full Green's function of the system.
Parity-Resolved Thermal States: When dealing with finite temperatures, this refers to preparing thermal states that are separated into even (symmetric) and odd (antisymmetric) parity sectors. This is achieved using a specialized quantum Gibbs sampler that simulates the Davies generator under a Lindbladian framework.

Summary

The gist The protocol introduces a symmetry-enabled direct quantum protocol, tailored quench spectroscopy (TQS), for computing both the real and imaginary parts of two-point time correlators without controlled unitaries or ancilla qubits in the correlator-estimation circuit<ref:2509.03288#pg12>. This method exploits parity symmetry to simplify readout, providing a practical route to estimating symmetry-resolved dynamical correlators on near-term and early fault-tolerant quantum hardware<ref:2509.03288#pg18>.

How it works

The protocol relies on exploiting parity symmetry, satisfied by a broad class of Hamiltonians in condensed matter physics and quantum chemistry, including the Fermi–Hubbard and Heisenberg models<ref:2509.03288#pg17>. The core idea is to introduce a tailored quench spectroscopy scheme that recovers both the real and imaginary parts of two-point time correlators, from which Green’s functions can be reconstructed via efficient classical signal processing<ref:2509.03288#pg18>.

The protocol is formalized by considering Hamiltonians that commute with a parity operator P, while the observables of interest A and B anti-commute with P<ref:2509.03288#pg18>. By designing tailored quench operators, the authors extend the standard QS method to access bosonic correlation functions, enabling direct evaluation of the real part of time correlators<ref:2509.03288#pg18>. Specifically, they introduce the quench operator URe:= P + A√2<ref:2509.03288#pg18>.

Symmetry and State Preparation

The protocol requires that the initial state be parity-symmetric, [ρ, P] = 0<ref:2509.03288#pg18>. For finite-temperature Green’s functions, the authors introduce a tailored symmetric quantum Gibbs sampler for preparing symmetric and antisymmetric thermal states<ref:2509.03288#pg18>. These states are defined as rhoS = ΠSρβΠS and rhoA = ΠAρβΠA, which decompose the thermal state into parity sectors<ref:2509.03288#pg18>.

Reconstruction and Advantages

The full complex correlator can be reconstructed from a linear combination of six experimentally accessible quench functions, simplifying to only two when the observables anti-commute with P<ref:2509.03288#pg18>. Theorem 1 states that for all t ∈ R, Im[C(A, B, t)] = Qρ,Be−iHtUIm and Re[C(A, B, t)] = p Qρ,Be−iHtURe when the initial state has definite parity p<ref:2509.03288#pg18>. This approach requires only time-independent Hamiltonian evolution, avoiding the simulation of time-dependent dynamics and reducing experimental complexity<ref:2509.03288#pg18>.

Complexity and Error Analysis

The complexity analysis shows that to achieve additive error ε with constant success probability, the number of shots per time point scales as O(ε−2)<ref:2509.03288#pg18>. Reconstructing the Green’s function G(z) requires resolving the frequency content of the measured timedomain correlator using classical signal processing methods such as MUltiple SIgnal Classification (MUSIC)<ref:2509.03288#pg18>. The total number of circuit executions scales as O(M ε−2) = O(∆f−1δt−1ε−2), up to logarithmic factors and constants depending on the noise level<ref:2509.03288#pg18>.

The error analysis demonstrates that symmetry-based symmetrization can quadratically suppress parity-breaking state-preparation errors, reducing the leading error contribution from O(µ) to O(µ2) under specific conditions<ref:2509.03288#pg18>. Numerical simulations on Fermi–Hubbard and XXZ instances support the robustness of the approach under Trotterization and imperfect state preparation<ref:2509.03288#pg18>.

Generalization to OTOCs

The approach can be extended to a class of n-point correlators, including out-of-time-order correlators (OTOCs)<ref:2509.03288#pg18>. For an OTOC(A, B, t), the experiment reduces to estimating the modified quench functions Qρβ,Ae iHtBe−iHtUIm/Re<ref:2509.03288#pg18>. This demonstrates that the same symmetry-based measurement primitive extends naturally to out-of-time-ordered correlators<ref:2509.03288#pg18>.

The protocol is compatible with both digital platforms and analog simulators, making it a promising first milestone for experimental demonstration in small Hubbard or spin-chain instances<ref:2509.03288#pg18>. This work opens several new exciting directions, including the measurement of higher-order and nonlinear response functions[48] and exploring analogous ideas in non-Hermitian or effective open-system settings<ref:2509.03288#pg18>. The paper concludes by benchmarking the protocol under realistic noise and comparing against other available approaches to clarify its practical advantage<ref:2509.03288#pg18>.

Improvements for AI systems

No-Ancilla Real/Imaginary Reconstruction: The improved AI system can compute the real and imaginary parts of two-point time correlators using only native time evolution and straightforward measurements available on current hardware platforms, specifically by implementing the tailored quench spectroscopy (TQS) protocol described in Theorem 1.
Finite-Temperature State Preparation: The system can prepare parity-resolved thermal states, such as the symmetric state ρS and antisymmetric state ρA, using a tailored symmetric quantum Gibbs sampler based on dissipative state preparation, which is achievable by simulating the Davies generator LD or LD,S under a Lindbladian framework.
OTOC Estimation: The protocol can be extended to estimating out-of-time-ordered correlators (OTOCs) by introducing tailored quench operators UIm and URe, allowing for the estimation of modified quench functions derived from the OTOC structure, as shown in Section III C.
Error Mitigation via Symmetrization: The system can suppress state-preparation errors quadratically by applying a symmetrization channel S(ρ) = ρ + P ρP2, which removes coherences between opposite parity sectors, thereby reducing the leading error contribution from O(µ) to O(µ 2).
Spectral Reconstruction via MUSIC: The system can reconstruct the Green’s function G(z) by using classical signal processing methods like MUltiple SIgnal Classification (MUSIC) on measured time-domain correlators, where the noise-subspace correlation function R(ω) provides a direct and visually stable diagnostic of the recovered spectral support.

Abstract

We present a symmetry-enabled direct quantum protocol for computing many-body Green's functions, a central tool for studying strongly correlated quantum systems. Our protocol relies only on native time evolution and straightforward measurements available on current hardware platforms. By exploiting parity symmetry -- satisfied by a broad class of Hamiltonians in condensed matter physics and quantum chemistry, including the Fermi--Hubbard and Heisenberg models -- we introduce a tailored quench spectroscopy scheme that recovers both the real and imaginary parts of two-point time correlators, from which Green's functions can be reconstructed via efficient classical signal processing. We further develop a tailored symmetric quantum Gibbs sampler that prepares parity-resolved (symmetric and antisymmetric) thermal states, enabling finite-temperature extensions within the same framework. Finally, we show that the same symmetry-based measurement primitive extends naturally to out-of-time-ordered correlators (OTOCs). Our results provide a practical route to estimating symmetry-resolved dynamical correlators on near-term and early fault-tolerant quantum hardware.

Sources

A circuit-differentiation framework for Green's functions on quantum computers