REVIEW 3 major objections 4 minor 57 references
Special Issue: Commemorating the 110th Anniversary of TANG Au-chin's Birthday Calculation of the Green's function on near-term quantum computers via Cartan decomposition
T0 review · 3 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read Constant-depth circuits compute Green's functions for any time
desk verdict Genuine but incremental: analytic gradients for Cartan decomposition plus a Green's-function pipeline; the two-site Hubbard test is solid, but the TFIM results and the fixed-depth claim are unverified. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The Cartan decomposition of the Hamiltonian's Lie algebra, a split g=k⊕m satisfying [k,k]⊂k, [m,m]⊂k, [k,m]⊂m, together with the KHK theorem: for H∈m there exists a fixed unitary K∈e^k and a Cartan subalgebra element h∈h such that H=KhK†. This turns time evolution into U(t)=K e^{-i h t} K†, and because h is Abelian, e^{-i h t} is a product of commuting Pauli rotations, giving fixed-depth circuits for arbitrary t. The paper finds K by minimizing the Killing form f(K)=Tr(K v K† H) over factorized unitaries K=∏ e^{i θ_i k_i}, with analytical gradients derived for efficient BFGS optimization; the gradients come from the factorized form, allowing each exponential to act on v and H as a Pauli rota
What would settle it
Run the Cartan-BFGS procedure on a transverse-field Ising chain with six or more sites, compute the residual norm ‖K†HK−h‖ after optimization, and compare the resulting U(t) against exact diagonalization for a long time t. If the residual is not at machine precision, or if the Green's function deviates from the exact result beyond rounding, the fixed-depth claim fails for that Hamiltonian.
Extended reading notes
Core claim
The central claim is that for a Hamiltonian H whose Pauli-string Lie algebra admits a Cartan decomposition with H in the noncompact part m, the real-time evolution operator is exactly U(t)=K_0 e^{-i h t} K_0†, with K_0 time-independent and h a Cartan subalgebra element. Because h is Abelian, e^{-i h t} is a product of commuting rotations around single Pauli strings, so its circuit depth is constant in t. The retarded Green's function is then obtained from ground-state overlaps of the form ⟨Ψ| c_a K_0 e^{-i h t} K_0† c_b† |Ψ⟩, evaluated by Hadamard tests, plus the conjugate term. Numerical experiments on the two-site Fermi-Hubbard model show the fixed-depth circuits reproduce the exact Green'
Load-bearing premise
The algorithm's correctness rests on the classical optimizer always reaching a global minimum of f(K)=Tr(KvK†H) so that K†HK exactly lands in the Cartan subalgebra; the paper demonstrates this on its benchmark Hamiltonians but provides no guarantee, and if the optimizer stalls the time-evolution operator is not e^{-iHt} at all.
Editorial extensions
If this is right
- For any Hamiltonian admitting the Cartan decomposition, the quantum circuit needed to evolve a state to arbitrary time has constant depth, eliminating Trotter error accumulation and variational drift over long times.
- The spectral function A_k(ω) is recovered by Fourier transforming the real-time Green's function, so the method gives direct access to Hubbard bands and quasiparticle features without summing over excited states.
- The construction can be reused for any initial state: once K_0 and h are fixed classically, each Green's function entry is just a Hadamard-test overlap on the same short circuit.
- The demonstrated crossover from a correlated metal at U=3 to a Mott insulator at U=6 in the two-site Fermi-Hubbard model shows the fixed-depth circuits preserve enough dynamical information to distinguish phases.
- The method separates the expensive classical search for K_0 from the quantum sampling, which is the natural division for near-term hardware.
Reading between the lines
- Beyond the paper's examples, the same fixed-depth identity applies to any retarded or advanced correlation function of the form ⟨Ψ| O_1 e^{-iHt} O_2 |Ψ⟩, since the time-evolution operator is the only part that needs U(t).
- If the classical optimization of K can be made locality-aware, the approach could move past the exponential scaling the authors note for the Lie algebra dimension; a natural test is a one-dimensional chain with local interactions where a block-structured K may suffice.
- The method's advantage is likely to be most visible on devices where two-qubit gate errors dominate: once K is compiled, the same circuit is reused for every time step, so calibration overhead is amortized over the whole spectral function.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a hybrid quantum–classical algorithm for computing retarded Green's functions on near-term quantum hardware. The key idea is to use a Cartan decomposition of the Hamiltonian Lie algebra: after a classical optimization step, the authors seek a unitary K and an Abelian subalgebra element h such that H = K h K†. Time evolution then takes the fixed-depth form U(t) = K e^{-i h t} K†, independent of t. The Green's function is assembled from Hadamard-test overlaps of time-evolved excited states, and spectral functions are obtained by Fourier transformation. The method is demonstrated on the two-site Fermi–Hubbard model with exact benchmarks in Figs. 1–3, and on transverse-field Ising chains of two, four, and six sites in Fig. 4. The conclusion explicitly acknowledges that the dimension of K scales exponentially and limits the method to small systems.
Significance. The central algebraic reduction is valuable and largely correct: Eq. (23) reduces the retarded Green's function to two time-evolution overlaps, Eq. (37) is the correct Jordan–Wigner two-site Hubbard Hamiltonian, and the basis h in Eq. (39) is pairwise commuting, as required for a Cartan subalgebra. For the Hubbard model, the exact comparisons in Figs. 1–3 give credible evidence that the Cartan decomposition succeeded for the parameters used. If the fixed-depth identity Eq. (19) can be reliably achieved, the algorithm offers a genuine advantage over Trotter-based simulation for small NISQ-era systems, and the analytical gradients in Eqs. (14)–(18) are a useful technical contribution. The paper is less convincing for the transverse-field Ising model, where no exact benchmark or optimization residual is provided, and the physical discussion of quantum criticality in Fig. 4 rests on unvalidated curves.
major comments (3)
- [One-dimensional Spin Chains / Fig. 4] The TFIM demonstration does not validate the central claim. Unlike Figs. 1–3, Fig. 4 contains only Cartan-derived spectral functions: no exact (diagonalization or Trotter) reference curves, no time-domain overlap fidelity, and no residual measure ||K†HK − h|| or BFGS success statistics are reported. Since Eq. (19) is exact only if the optimization of Eq. (6) actually finds K satisfying Eq. (7), the TFIM results cannot be distinguished from an optimization failure. Please add exact benchmarks for at least one system size and report residual norms or success statistics for all Cartan decompositions used in the paper.
- [Method, Eqs. (6)–(9)] The statement after Eq. (8) that the exponential form K = e^{Σ ia_i k_i} and the product form K = ∏ e^{ia_i k_i} are 'equivalent' is not justified, and it is load-bearing. A local minimum of f(K) over the fixed-order product submanifold Eq. (9) is only a constrained stationary point; the tangent space of that submanifold is not obviously all of k at every point, so the critical-point condition that would imply K†HK ∈ h need not hold. BFGS has no convergence guarantee for this nonconvex problem. The paper should either prove the equivalence for the relevant Lie algebras or report residual norms and multiple-random-restart success rates; otherwise the identity U(t) = K e^{-iht} K† is an assumption rather than a derived result.
- [Abstract / Conclusion] The abstract claims 'an efficient algorithm... requires only fixed-depth quantum circuits for arbitrarily long time simulations.' The conclusion correctly states that the dimension of K scales exponentially and limits the method to roughly 10 qubits. This is a significant qualification: the fixed-depth property holds only after a successful classical optimization whose cost is exponential and whose success is, for TFIM, unverified. Please state this limitation prominently in the abstract or introduction, and separate the quantum circuit depth from the classical preprocessing cost in the claims.
minor comments (4)
- [Eq. (38)] The last basis element of k is written as Z0Y1Z2X0, which repeats Z0 and is inconsistent with the pattern of the other elements; it should almost certainly be Z0Y1Z2X3.
- [Eq. (24)] The Jordan–Wigner expression is written as Q_a^† ⊗ Z_{a−1} ⊗ ... ⊗ Z_1; the ordering of the tensor factors should be specified explicitly to avoid ambiguity about which qubit is the least significant.
- [One-dimensional Spin Chains] The notation k is used both for momentum labels in the Hubbard section and for the system-size index in Eq. (45), where 'different k means different sizes.' This is confusing; use N or n for the chain length.
- [Throughout] There are several typographical and grammatical issues: 'homomrophism', 'Cartan decomposition has been emerged', and 'We use [PiPj to represent' should be corrected.
Circularity Check
No significant circularity: the Cartan decomposition is an exact Hamiltonian identity, and the Green's-function predictions are benchmarked against independent exact results.
full rationale
The paper's central construction is not circular. The parameters θ in K are obtained by minimizing f(K)=Tr(KvK†H) (Eq. 6), a classical optimization whose target is the exact identity H=K0 h K0† (Eq. 7). This is a Hamiltonian decomposition, not a fit to Green's-function data; the Green's function is then computed from the resulting U(t)=K0 e^{-iht} K0† (Eq. 19). No quantity entering the Cartan decomposition is defined in terms of the target Green's function or spectral function, so the prediction is not forced by construction. The numerical benchmarks in Figs. 1–3 compare against independently computed exact results, providing external validation. The only self-citation (ref. [12]) appears in a broad list of prior quantum algorithms and is not load-bearing; no uniqueness claim or ansatz is imported from the authors' own prior work. Concerns about BFGS convergence or the absence of exact TFIM comparisons are correctness/robustness issues, not circularity.
Assumptions & free parameters
free parameters (3)
- eta (spectral broadening) =
0.2
- Cartan parameters theta_i of K =
not reported (determined by BFGS)
- gamma (coefficients of generic v) =
unspecified transcendental number
assumptions (4)
- standard math KHK theorem: for any m in m there exists K in e^k and h in h with m = K h K†.
- domain assumption The Hamiltonian H belongs to m for the chosen Cartan decomposition g = k ⊕ m.
- domain assumption BFGS minimization of the nonconvex f(K) reaches a global (not merely local) critical point with K†HK ∈ h.
- domain assumption The exact ground state is available (via ADAPT-VQE) and Hadamard tests give the required overlaps.
Cite this review
Pith. "Pith review of Special Issue: Commemorating the 110th Anniversary of TANG Au-chin's Birthday Calculation of the Green's function on near-term quantum computers via Cartan decomposition." pith.science (2026). https://pith.science/paper/26YNDIIY
@misc{pith2026250909248,
author = {Pith},
title = {Pith review of: Special Issue: Commemorating the 110th Anniversary of TANG Au-chin's Birthday Calculation of the Green's function on near-term quantum computers via Cartan decomposition},
year = {2026},
howpublished = {\url{https://pith.science/paper/26YNDIIY}},
note = {Machine review of arXiv:2509.09248}
}
read the original abstract
Accurate computation of the Green's function is crucial for connecting experimental observations to the underlying quantum states. A major challenge in evaluating the Green's function in the time domain lies in the efficient simulation of quantum state evolution under a given Hamiltonian-a task that becomes exponentially complex for strongly correlated systems on classical computers. Quantum computing provides a promising pathway to overcome this barrier by enabling efficient simulation of quantum dynamics. However, for near-term quantum devices with limited coherence times and fidelity, the deep quantum circuits required to implement time-evolution operators present a significant challenge for practical applications. In this work, we introduce an efficient algorithm for computing Green's functions via Cartan decomposition, which requires only fixed-depth quantum circuits for arbitrarily long time simulations. Additionally, analytical gradients are formulated to accelerate the Cartan decomposition by leveraging a unitary transformation in a factorized form. The new algorithm is applied to simulate long-time Green's functions for the Fermi-Hubbard and transverse-field Ising models, extracting the spectral functions through Fourier transformation.
Figures
Reference graph
Works this paper leans on
-
[1]
R. P. Feynman,Int. J. Theor. Phys., 1982,21, 467–488
1982
-
[2]
J. Du, N. Xu, X. Peng, P. Wang, S. Wu and D. Lu,Phys. Rev. Lett., 2010,104, 030502
2010
-
[3]
Peruzzo, J
A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik and J. L. O’ Brien, Nat. Commun., 2014,5, 4213
2014
-
[4]
P. J. J. O’ Malley, R. Babbush, I. D. Kivlichan, J. Romero, J. R. McClean, R. Barends, J. Kelly, P. Roushan, A. Tranter, N. Ding, B. Campbell, Y. Chen, Z. Chen, B. Chiaro, A. Dunsworth, A. G. Fowler, E. Jef- frey, E. Lucero, A. Megrant, J. Y. Mutus, M. Neeley, C. Neill, C. Quintana, D. Sank, A. Vainsencher, J. Wen- ner, T. C. White, P. V. Coveney, P. J. L...
2016
-
[5]
Kandala, A
A. Kandala, A. Mezzacapo, K. Temme, M. Takita, M. Brink, J. M. Chow and J. M. Gambetta,Nature, 2017, 549, 242
2017
-
[6]
Hempel, C
C. Hempel, C. Maier, J. Romero, J. McClean, T. Monz, H. Shen, P. Jurcevic, B. P. Lanyon, P. Love, R. Babbush, A. Aspuru-Guzik, R. Blatt and C. F. Roos,Phys. Rev. X, 2018,8, 031022
2018
-
[7]
Arute, K
F. Arute, K. Arya, R. Babbush, D. Bacon, J. C. Bardin, R. Barends, S. Boixo, M. Broughton, B. B. Buck- ley, D. A. Buell, B. Burkett, N. Bushnell, Y. Chen, Z. Chen, B. Chiaro, R. Collins, W. Courtney, S. Demura, A. Dunsworth, E. Farhi, A. Fowler, B. Foxen, C. Gid- ney, M. Giustina, R. Graff, S. Habegger, M. P. Harri- gan, A. Ho, S. Hong, T. Huang, W. J. Hu...
2020
-
[8]
Y. Cao, J. Romero, J. P. Olson, M. Degroote, P. D. John- son, M. Kieferov´ a, I. D. Kivlichan, T. Menke, B. Per- opadre, N. P. D. Sawaya, S. Sim, L. Veis and A. Aspuru- Guzik,Chem. Rev., 2019,119, 10856–10915
2019
Show all 57 references
-
[9]
McArdle, S
S. McArdle, S. Endo, A. Aspuru-Guzik, S. C. Benjamin and X. Yuan,Rev. Mod. Phys., 2020,92, 015003
2020
-
[10]
Zhong, H
H.-S. Zhong, H. Wang, Y.-H. Deng, M.-C. Chen, L.-C. Peng, Y.-H. Luo, J. Qin, D. Wu, X. Ding, Y. Hu, P. Hu, X.-Y. Yang, W.-J. Zhang, H. Li, Y. Li, X. Jiang, L. Gan, G. Yang, L. You, Z. Wang, L. Li, N.-L. Liu, C.-Y. Lu and J.-W. Pan,Science, 2020,370, 1460–1463
2020
-
[11]
Cerezo, R
M. Cerezo, R. B. Andrew Arrasmith, S. C. Benjamin, S. Endo, K. Fujii, J. R. McClean, K. Mitarai, X. Yuan, L. Cincio and P. J. Coles,Nat. Rev. Phys., 2021,3, 625– 644
2021
-
[12]
J. Liu, Y. Fan, Z. Li and J. Yang,Chem. Soc. Rev., 2022, 51, 3263–3279
2022
-
[13]
Shang, L
H. Shang, L. Shen, Y. Fan, Z. Xu, C. Guo, J. Liu, W. Zhou, H. Ma, R. Lin, Y. Yang, F. Li, Z. Wang, 9 Y. Zhang and Z. Li, Proceedings of the International Conference on High Performance Computing, Network- ing, Storage and Analysis, 2022
2022
-
[14]
W. J. Huggins, B. A. O’Gorman, N. C. Rubin, D. R. Reichman, R. Babbush and J. Lee,Nature, 2022,603, 416–420
2022
-
[15]
Robledo-Moreno, M
J. Robledo-Moreno, M. Motta, H. Haas, A. Javadi- Abhari, P. Jurcevic, W. Kirby, S. Martiel, K. Sharma, S. Sharma, T. Shirakawa, I. Sitdikov, R.-Y. Sun, K. J. Sung, M. Takita, M. C. Tran, S. Yunoki and A. Mezza- capo,Sci. Adv., 2025,11, eadu9991
2025
-
[16]
A. L. Fetter and J. D. Walecka,Quantum theory of many- particle systems, Courier Corporation, 2012
2012
-
[17]
A. R. Williams, P. J. Feibelman and N. D. Lang,Phys. Rev. B, 1982,26, 5433
1982
-
[18]
Kotliar, S
G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. Parcollet and C. A. Marianetti,Rev. Mod. Phys., 2006,78, 865–951
2006
-
[19]
McClain, J
J. McClain, J. Lischner, T. Watson, D. A. Matthews, E. Ronca, S. G. Louie, T. C. Berkelbach and G. K.-L. Chan,Phys. Rev. B, 2016,93, 235139
2016
-
[20]
K. M. Nakanishi, K. Mitarai and K. Fujii,Phys. Rev. Res., 2019,1, 033062
2019
-
[21]
S. Endo, I. Kurata and Y. O. Nakagawa,Phys. Rev. Res., 2020,2, 033281
2020
-
[22]
R. M. Parrish, E. G. Hohenstein, P. L. McMahon and T. J. Mart ´ ınez,Phys. Rev. Lett., 2019,122, 230401
2019
-
[23]
Rizzo, F
J. Rizzo, F. Libbi, F. Tacchino, P. J. Ollitrault, N. Marzari and I. Tavernelli,Phys. Rev. Res., 2022,4, 043011
2022
-
[24]
P. J. Ollitrault, A. Kandala, C.-F. Chen, P. K. Barkout- sos, A. Mezzacapo, M. Pistoia, S. Sheldon, S. Woerner, J. M. Gambetta and I. Tavernelli,Phys. Rev. Research, 2020,2, 043140
2020
-
[25]
T. E. Baker,Phys. Rev. A, 2021,103, 032404
2021
-
[26]
Jamet, A
F. Jamet, A. Agarwal, C. Lupo, D. E. Browne, C. Weber and I. Rungger,Krylov variational quantum algorithm for first principles materials simulations, 2021,https: //arxiv.org/abs/2105.13298
2021 arXiv
-
[27]
Cai, W.-H
X. Cai, W.-H. Fang, H. Fan and Z. Li,Phys. Rev. Res., 2020,2, 033324
2020
-
[28]
Huang, X
K. Huang, X. Cai, H. Li, Z.-Y. Ge, R. Hou, H. Li, T. Liu, Y. Shi, C. Chen, D. Zheng, K. Xu, Z.-B. Liu, Z. Li, H. Fan and W.-H. Fang,J. Chem. Theory Comput., 2022,13, 9114–9121
2022
-
[29]
H. Chen, M. Nusspickel, J. Tilly and G. H. Booth,Phys. Rev. A, 2021,104, 032405
2021
-
[30]
K¨ okc¨ u, H
E. K¨ okc¨ u, H. Labib, J. Freericks and A. F. Kemper,Nat. Commun., 2024,15, 3881
2024
-
[31]
D. W. Berry, G. Ahokas, R. Cleve and B. C. Sanders, Commun. Math Phys., 2007,270, 359–371
2007
-
[32]
Roggero and J
A. Roggero and J. Carlson,Phys. Rev. C, 2019,100, 034610
2019
-
[33]
Kosugi and Y.-i
T. Kosugi and Y.-i. Matsushita,Phys. Rev. A, 2020,101, 012330
2020
-
[34]
D. W. Berry, A. M. Childs, R. Cleve, R. Kothari and R. D. Somma,Phys. Rev. Lett., 2015,114, 090502
2015
-
[35]
Preskill,Quantum, 2018,2, 79
J. Preskill,Quantum, 2018,2, 79
2018
-
[36]
Li and S
Y. Li and S. C. Benjamin,Phys. Rev. X, 2017,7, 021050
2017
-
[37]
S. Endo, J. Sun, Y. Li, S. C. Benjamin and X. Yuan, Phys. Rev. Lett., 2020,125, 010501
2020
-
[38]
Libbi, J
F. Libbi, J. Rizzo, F. Tacchino, N. Marzari and I. Taver- nelli,Phys. Rev. Res., 2022,4, 043038
2022
-
[39]
Sakurai, W
R. Sakurai, W. Mizukami and H. Shinaoka,Phys. Rev. Res., 2022,4, 023219
2022
-
[40]
Y.-X. Yao, N. Gomes, F. Zhang, C.-Z. Wang, K.-M. Ho, T. Iadecola and P. P. Orth,PRX Quantum, 2021,2, 030307
2021
-
[41]
Gyawali and M
G. Gyawali and M. J. Lawler,Phys. Rev. A, 2022,105, 012413
2022
-
[42]
Mootz, T
M. Mootz, T. Iadecola and Y.-X. Yao,J. Chem. Theory Comput., 2024,20, 8689–8710
2024
-
[43]
K¨ okc¨ u, T
E. K¨ okc¨ u, T. Steckmann, Y. Wang, J. K. Freericks, E. F. Dumitrescu and A. F. Kemper,Phys. Rev. Lett., 2022, 129, 070501
2022
-
[44]
Steckmann, T
T. Steckmann, T. Keen, E. K¨ okc¨ u, A. F. Kemper, E. F. Dumitrescu and Y. Wang,Phys. Rev. Res., 2023,5, 023198
2023
-
[45]
H. R. Grimsley, S. E. Economou, E. Barnes and N. J. Mayhall,Nat. Commun., 2019,10, 1–9
2019
-
[46]
B. C. Hall, inQuantum Theory for Mathematicians, Springer, 2013, pp. 333–366
2013
-
[47]
Khaneja and S
N. Khaneja and S. J. Glaser,Chem. Phys., 2001,267, 11–23
2001
-
[48]
Jordan and E
P. Jordan and E. Wigner,Z. Phys., 1928,47, 631–651
1928
-
[49]
Aharonov, V
D. Aharonov, V. Jones and Z. Landau, Proceedings of the thirty-eighth annual ACM symposium on Theory of computing, 2006, pp. 427–436
2006
-
[50]
M. C. Gutzwiller,Phys. Rev. Lett., 1963,10, 159
1963
-
[51]
Kanamori,Prog
J. Kanamori,Prog. Theor. Phys., 1963,30, 275–289
1963
-
[52]
R. S. (London),Proceedings of the Royal Society, 1869, vol. 17
-
[53]
J. R. McClean, N. C. Rubin, K. J. Sung, I. D. Kivlichan, X. Bonet-Monroig, Y. Cao, C. Dai, E. S. Fried, C. Gid- ney, B. Gimbyet al.,Quantum Sci. Technol., 2020,5, 034014
2020
-
[54]
N. F. Mott,Rev. Mod. Phys., 1968,40, 677
1968
-
[55]
Meinders, H
M. Meinders, H. Eskes and G. Sawatzky,Phys. Rev. B, 1993,48, 3916
1993
-
[56]
Radicevic,arXiv preprint arXiv:1809.07757, 2018
D. Radicevic,arXiv preprint arXiv:1809.07757, 2018
2018 arXiv
-
[57]
B. T. McDonough, C. Yin, A. Lucas and C. Zhang,arXiv preprint arXiv:2502.02652, 2025
2025
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.