REVIEW 4 major objections 4 minor 29 references
This paper gives a quantum algorithm that ε-approximates modular flow in O~(κ²|t| log(κ²/ε)) queries, with a matching lower bound.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
A QSVT-based algorithm implements modular flow of an operator with respect to a density matrix in O~(κ²|t|log(κ²/ε)) queries to a block encoding of the state.
T0 review reviewed 2026-08-05 challenge →
load-bearing objection Modular flow via QSVT is a good idea, but the central lemma has a bogus error bound and the advertised complexity doesn't follow. the 4 major comments →
A quantum algorithm for modular flow
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
The paper's central claim is that modular flow is efficiently implementable from a block encoding of ρ. Theorem 1 states that for a block-encoded ρ with smallest nonzero eigenvalue 1/κ, one can construct Õ(t) with ||Õ(t)-ρ^{-it}Oρ^{it}||<ε using O~(κ²|t| log(κ²/ε)) queries to the encoding. The proof constructs an explicit polynomial approximating (1/2 log(1/2κ)) log(x²) on the spectrum outside a gap of width 1/κ, uses the quantum singular value transform to turn it into an approximate modular Hamiltonian, and feeds that into Hamiltonian simulation to produce ρ^{it}. A reduction from von Neumann entropy estimation to modular flow then gives a query lower bound, so the algorithm's scaling cann
What carries the argument
The modular Hamiltonian polynomial: a degree-2N Chebyshev partial sum f_N(x) = -log 2 + Σ_{n=1}^N T_{2n}(x)(-1)^{n-1}/n, rescaled by 1/log(1/2κ) and multiplied by a rectangular polynomial that suppresses |x|≤1/2κ, making it bounded by 1 and of definite parity. This polynomial is admissible for the quantum singular value transform (QSVT), which applies it to the singular values of a block-encoded ρ, producing an ε-approximation to -log ρ; Hamiltonian simulation then applies ρ^{it}. The degree N=O(κ² log(κ²/ε)) of the log polynomial is the dominant cost.
Load-bearing premise
The complexity guarantee is measured only against an oracle that block-encodes ρ; if no such encoding can be made—likely for holographic states—the algorithm's cost bounds do not carry over.
What would settle it
Take a 2×2 state ρ = diag(1/κ, 1−1/κ), classically compute ρ^{-i}Oρ^i for a fixed traceless operator O, then run the algorithm with a block-encoding circuit and count the queries. If the query count exceeds C κ²|t| log(κ²/ε) for a fixed constant C, or the implemented operator differs from the exact flow by more than ε in spectral norm, Theorem 1 is false.
If this is right
- For any block-encoded ρ with smallest nonzero eigenvalue 1/κ, the modular flow of any operator can be ε-approximated at modular time t using O~(κ²|t| log(κ²/ε)) queries, making modular flow a concrete computational primitive rather than a formal object.
- The reduction from von Neumann entropy estimation to modular flow implies any modular-flow algorithm must scale polynomially with the dimension n, so the κ² dependence in Theorem 1 cannot be improved to log n.
- In a topologically ordered gapped ground state, running modular flow across a tripartite disk and measuring the entanglement entropy of one region reveals the chiral central charge c_- through the slope πc_-/3 of the entropy as a function of modular time.
- Given oracle access to a block encoding of the reduced state and to operator insertions, modular flow plus quantum phase estimation can evaluate the correlation function proposed for probing behind the horizon in the holographic bulk-observer scenario.
- When modular flow is applied to a purification |ψ⟩, the κ-dependence can be tamed: choosing κ=O(d²/δ²) gives error at most δ with cost linear in |t| and polynomial in the dimension d.
Where Pith is reading between the lines
- I infer that the Chebyshev expansion of log(x²) is the true workhorse: the same technique should apply to other operator functions with logarithmic singularities, giving κ²-type costs whenever the spectrum has a gap.
- The paper's caveat that holographic states may lack efficient block encodings suggests the practical bottleneck is state preparation, not flow simulation; a testable next step is finding circuit families for which ρ^{it} is easy even though ρ itself is hard to block-encode.
- If the purification construction is combined with existing ground-state preparation oracles for gapped systems, modular flow could become a practical probe of topological order at modest system sizes, since the cost no longer singles out tiny eigenvalues.
- The polylog(n) entropy estimate for fixed κ indicates that the hardness of von Neumann entropy estimation resides in resolving tiny eigenvalues; for states with known condition number, entropy estimation may be much easier than the worst-case n² lower bound.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a QSVT-based algorithm for implementing modular flow ρ^{-it}Oρ^{it} for a block-encoded density matrix ρ. The central construction approximates -log ρ by a truncated Chebyshev series, then feeds the resulting approximate modular Hamiltonian into a QSVT Hamiltonian-simulation routine. The authors claim a query complexity of O~(κ²|t| log(κ²/ε)), where 1/κ is the smallest nonzero eigenvalue of ρ. They also discuss applications to chiral central charge extraction, to the Jafferis–Lamprou/Gao–Lamprou holographic observer proposal, and to von Neumann entropy estimation, and they claim a query-complexity lower bound for modular flow based on the known entropy-estimation lower bound. The paper is clearly written and situates the algorithm within standard QSVT tools, but the main technical lemma and the lower-bound argument have serious defects.
Significance. If correct, the paper would provide the first explicit provable algorithm for modular flow and would connect a well-studied physics quantity to the QSVT toolbox. The applications, especially to chiral central charge and to holographic bulk-observer probes, are potentially interesting. The paper also has the virtue of being constructive and explicit. However, the central polynomial approximation lemma is demonstrably incorrect, and the claimed lower bound is not a lower bound. The significance is therefore conditional: the conceptual idea is natural, but the results as stated are not supported.
major comments (4)
- [§3.1, Lemma 2] The error bound in Lemma 2 is false. The proof identifies the tail of ∑ T_{2n}(x)(-1)^{n-1}/n with the Maclaurin tail of log(1+z) at z=cos(2 arccos x), but the actual tail at x=cos θ is ∑_{n>N} (-1)^{n-1} cos(2nθ)/n = Re ∑_{n>N} (-1)^{n-1} e^{2inθ}/n, which is not a power series in cos(2θ). At x=1 the tail is the alternating harmonic tail ∑_{n>N} (-1)^{n-1}/n ~ 1/(2N), not exponentially small. At x=1/κ, the tail is -∑_{n>N} cos(2n/κ)/n, and Abel summation with |∑_{n=1}^M cos(2n/κ)| ≤ 1/sin(1/κ) gives |tail| ≤ 2κ/(N+1). Thus for N=O(κ² log(κ²/ε)) the error at x=1/κ is ~1/(κ log(κ²/ε)), which is far larger than ε for large κ and small ε. Consequently Eq. (12), the degree bound Eq. (13), and the query complexity of Theorem 1 are unsupported.
- [§4.3, 'query complexity lower bound'] This section does not prove a lower bound on modular flow. The argument says that the von Neumann entropy lower bound Ω(n²/ε) forces any modular flow implementation to scale polynomially in n. But a reduction from entropy estimation to modular flow only shows that the particular entropy-estimation algorithm built on modular flow must respect the entropy lower bound; it does not constrain the complexity of an arbitrary modular-flow oracle. The inequality in Eq. (47) compares an upper bound on the proposed entropy algorithm with a known lower bound; that is a consistency check, not a lower bound on modular flow. The sentence 'There can be no modular flow implementation scaling like log n' is a non sequitur. A genuine lower bound would require showing that any algorithm implementing modular flow with q queries can be converted into an entropy estimator whose complexity violates the known bo
- [§4.3, Eq. (39)–(41)] The QPE-based entropy estimation omits the phase-precision cost. The entropy is recovered from the mean of phases θ_n = π log λ_n / log κ. To achieve additive error ε in S(ρ), the phases must be estimated to precision O(ε/log κ), and QPE then requires O(log κ/ε) controlled-U applications per sample (or a comparable number of ancilla qubits). This factor is absent from the sample count in Eq. (40) and from the total query count in Eq. (41). In addition, Eq. (33) defines QPE(ρ,U) as ∑ e^{iθ_k}⟨k|ρ|k⟩ = Tr(ρU), while Eq. (38) uses QPE(V)_ρ = ∑ θ_n λ_n; these are different quantities, and QPE does not directly output Tr(ρU) for a general U. The section therefore does not provide a valid estimate of von Neumann entropy with the stated complexity.
- [§4.2, Eq. (33)–(36)] The claimed computation of the correlation function W(s,t) by QPE is not justified. For a mixed state ρ and a unitary U that does not commute with ρ, QPE on ρ gives a distribution over U-eigenstates with weights ⟨k|ρ|k⟩ in the eigenbasis of U; it does not return Tr(ρU). The operators ilde U = ρ^{-is}U_1ρ^{is}U_2 are not generally diagonal in the eigenbasis of ρ, so the identity in Eq. (35) does not follow. A controlled-Hadamard test on a purification would be a standard alternative, but no such procedure is described. Thus the holographic application, as written, is not established.
minor comments (4)
- [§2.1] The definition of a block encoding is imprecise: 'acts on H_R⊗H_A as |0⟩⟨0|⊗O' should be 'there exists a unitary U whose top-left block equals O/α', with α≥∥O∥.
- [Eq. (33) and Eq. (38)] The notation QPE(ρ,U) is used inconsistently: Eq. (33) writes it as a trace with e^{iθ_k}, while Eq. (38) writes it as a phase average with θ_n. These are different quantities and should be distinguished.
- [Appendix A and Conclusions] Typos: 'Cheybshev' should be 'Chebyshev'; 'langauge' should be 'language'.
- [Appendix B and §3.2] The signs around log(1/2κ) are confusing. The text says P^{MH} uses 1/(2 log(1/2κ)), which is negative for κ>1/2, while Appendix B says f^{log}_{κ',ε} is bounded by log(2κ). Please normalize consistently to avoid ambiguity.
Circularity Check
No significant circularity: the central claim derives from QSVT and an external entropy lower bound; self-citations appear only in applications.
full rationale
The paper's central derivation is self-contained. The modular flow algorithm approximates -log ρ by constructing a Chebyshev polynomial for log(x^2) (Lemmas 1-2) and then applies QSVT-based Hamiltonian simulation. No parameter is fitted to the modular flow itself, and the target quantity ρ^{-it}Oρ^{it} is not used to define the logarithmic approximation. The query-complexity bound follows from the polynomial degree and the known Hamiltonian simulation scaling; it is not obtained by assuming the conclusion. The lower-bound argument in Section 4.3 reduces von Neumann entropy estimation to modular flow and then invokes an external lower bound from [22]; it does not assume the modular-flow query complexity, so the reduction is a legitimate complexity-theoretic implication rather than a circular step. Self-citations [9,10,13] occur only in the applications discussion (chiral central charge, modular commutator) and are not load-bearing for Theorem 1 or Theorem 2. The skeptical concern about Lemma 2's error bound is a mathematical-correctness issue, not a circularity issue; even if the bound were wrong, that would invalidate the proof but not constitute the derivation reducing to its own input. Overall, no circular step can be exhibited from the manuscript's own equations.
Axiom & Free-Parameter Ledger
free parameters (2)
- rectangular function parameter ε' =
min(2ε/(5 log κ), log(1/2κ)/N(κ,ε))
- truncation point κ in Section 3.4 =
(d/(|t|ε))^{2/3}
axioms (4)
- standard math QSVT can implement any polynomial of the singular values of a block-encoded operator with degree equal to the number of encoding queries.
- standard math Hamiltonian simulation via QSVT has query complexity Θ(α|t| + log(1/ε)/log(e + log(1/ε)/(α|t|))) for a Hamiltonian with norm bound α.
- domain assumption The algorithm has oracle access to a block encoding of ρ.
- domain assumption The known query lower bound Ω(n² log(1/δ)/ε) for von Neumann entropy estimation holds in the same oracle model.
Cite this review
Pith. "Pith review of A quantum algorithm for modular flow." pith.science (2026). https://pith.science/paper/7VXZ7D6P
@misc{pith2026250816826,
author = {Pith},
title = {Pith review of: A quantum algorithm for modular flow},
year = {2026},
howpublished = {\url{https://pith.science/paper/7VXZ7D6P}},
note = {Machine review of arXiv:2508.16826}
}
read the original abstract
Entanglement is a defining property of quantum systems. For a subsystem of a larger quantum system, one can formally define an operator known as the modular Hamiltonian, which is closely linked to the entanglement properties of that subsystem, and a corresponding operator flow called the modular flow. Algorithms for estimating the von Neumann entropy, the best-known entanglement measure, are well-established, but no equivalent procedures have been previously described for the modular flow. In this work, we briefly review the quantum singular value transform (QSVT) framework for developing quantum algorithms, and then discuss the implementation of modular flow within this framework. We conclude by describing select applications of our modular flow algorithm, such as extracting the chiral central charge of a topologically ordered system and simulating the experience of the bulk observer in holography. We also prove a query complexity lower bound for modular flow, which shows that our method cannot be improved further substantially.
Figures
Reference graph
Works this paper leans on
-
[1]
Topological entanglement entropy.Physical review letters, 96(11):110404, 2006
Alexei Kitaev and John Preskill. Topological entanglement entropy.Physical review letters, 96(11):110404, 2006
work page 2006
-
[2]
Michael Levin and Xiao-Gang Wen. Detecting topological order in a ground state wave function.Physical review letters, 96(11):110405, 2006
work page 2006
-
[3]
The large-n limit of superconformal field theories and supergravity
Juan Maldacena. The large-n limit of superconformal field theories and supergravity. International journal of theoretical physics, 38(4):1113–1133, 1999
work page 1999
-
[4]
Shinsei Ryu and Tadashi Takayanagi. Holographic Derivation of Entanglement En- tropy from AdS/CFT.Physical Review Letters, 96(18):181602, May 2006. arXiv:hep- th/0603001
-
[5]
A finite entanglement entropy and the c-theorem.Physics Letters B, 600(1-2):142–150, 2004
H Casini and M Huerta. A finite entanglement entropy and the c-theorem.Physics Letters B, 600(1-2):142–150, 2004
work page 2004
-
[6]
H Casini and Marina Huerta. Renormalization group running of the entanglement entropy of a circle.Physical Review D—Particles, Fields, Gravitation, and Cosmology, 85(12):125016, 2012
work page 2012
-
[7]
Fusion rules from entanglement.Annals of Physics, 418:168164, 2020
Bowen Shi, Kohtaro Kato, and Isaac H Kim. Fusion rules from entanglement.Annals of Physics, 418:168164, 2020
work page 2020
-
[8]
Uni- versal lower bound on topological entanglement entropy
Isaac H Kim, Michael Levin, Ting-Chun Lin, Daniel Ranard, and Bowen Shi. Uni- versal lower bound on topological entanglement entropy. Physical Review Letters, 131(16):166601, 2023
work page 2023
-
[9]
Chiral central charge from a single bulk wave function.Physical Review Letters, 128(17):176402, 2022
Isaac H Kim, Bowen Shi, Kohtaro Kato, and Victor V Albert. Chiral central charge from a single bulk wave function.Physical Review Letters, 128(17):176402, 2022
work page 2022
-
[10]
Kim, Bowen Shi, Kohtaro Kato, and Victor V
Isaac H. Kim, Bowen Shi, Kohtaro Kato, and Victor V. Albert. Modular commutator in gapped quantum many-body systems.Physical Review B, 106(7):075147, August
-
[11]
Strict area law entangle- ment versus chirality, 2024
Xiang Li, Ting-Chun Lin, John McGreevy, and Bowen Shi. Strict area law entangle- ment versus chirality, 2024
work page 2024
-
[12]
Ruihua Fan. From entanglement generated dynamics to the gravitational anomaly and chiral central charge.Physical Review Letters, 129(26):260403, 2022
work page 2022
-
[13]
Modular commutators in conformal field theory.Physical Review Letters, 129(26):260402, 2022
Yijian Zou, Bowen Shi, Jonathan Sorce, Ian T Lim, and Isaac H Kim. Modular commutators in conformal field theory.Physical Review Letters, 129(26):260402, 2022
work page 2022
-
[14]
Bulk locality from modular flow, April 2018
Thomas Faulkner and Aitor Lewkowycz. Bulk locality from modular flow, April 2018. arXiv:1704.05464
Pith/arXiv arXiv 2018
-
[15]
Jafferis, Aitor Lewkowycz, Juan Maldacena, and S
Daniel L. Jafferis, Aitor Lewkowycz, Juan Maldacena, and S. Josephine Suh. Relative entropy equals bulk relative entropy, December 2015. arXiv:1512.06431
Pith/arXiv arXiv 2015
-
[16]
Inside the Hologram: Reconstructing the bulk observer’s experience, August 2021
Daniel Louis Jafferis and Lampros Lamprou. Inside the Hologram: Reconstructing the bulk observer’s experience, August 2021. arXiv:2009.04476 [hep-th]
Pith/arXiv arXiv 2021
-
[17]
On black hole in- terior reconstruction, singularities and the emergence of time, November 2022
Jan de Boer, Daniel Louis Jafferis, and Lampros Lamprou. On black hole in- terior reconstruction, singularities and the emergence of time, November 2022. arXiv:2211.16512 [hep-th]
Pith/arXiv arXiv 2022
-
[18]
Realizing topologically ordered states on a quantum processor.Science, 374(6572):1237–1241, 2021
KJ Satzinger, Y-J Liu, A Smith, C Knapp, M Newman, C Jones, Zhaoshi Chen, Céline Quintana, Xiaohua Mi, A Dunsworth, et al. Realizing topologically ordered states on a quantum processor.Science, 374(6572):1237–1241, 2021
work page 2021
-
[19]
Ruben Verresen, Mikhail D Lukin, and Ashvin Vishwanath. Prediction of toric code topological order from rydberg blockade.Physical Review X, 11(3):031005, 2021
work page 2021
-
[20]
András Gilyén, Yuan Su, Guang Hao Low, and Nathan Wiebe. Quantum singu- lar value transformation and beyond: exponential improvements for quantum matrix 16 arithmetics. InProceedings of the 51st Annual ACM SIGACT Symposium on Theory of Computing, pages 193–204, June 2019. arXiv:1806.01838 [quant-ph]
Pith/arXiv arXiv 2019
-
[21]
Seeing behind black hole horizons in SYK.Jour- nal of High Energy Physics, 2022(6):143, June 2022
Ping Gao and Lampros Lamprou. Seeing behind black hole horizons in SYK.Jour- nal of High Energy Physics, 2022(6):143, June 2022. arXiv:2111.14010 [cond-mat, physics:gr-qc, physics:hep-th, physics:quant-ph]
Pith/arXiv arXiv 2022
-
[22]
Sublinear quantum algorithms for estimating von Neumann entropy, November 2021
Tom Gur, Min-Hsiu Hsieh, and Sathyawageeswar Subramanian. Sublinear quantum algorithms for estimating von Neumann entropy, November 2021. arXiv:2111.11139 [quant-ph]
Pith/arXiv arXiv 2021
-
[23]
John M. Martyn, Zane M. Rossi, Andrew K. Tan, and Isaac L. Chuang. A Grand Unification of Quantum Algorithms. PRX Quantum, 2(4):040203, December 2021. arXiv: 2105.02859
Pith/arXiv arXiv 2021
-
[24]
Guang Hao Low and Isaac L. Chuang. Optimal Hamiltonian Simulation by Quan- tum Signal Processing. Physical Review Letters, 118(1):010501, January 2017. arXiv:1606.02685 [quant-ph]
Pith/arXiv arXiv 2017
-
[25]
Guang Hao Low and Isaac L. Chuang. Hamiltonian Simulation by Qubitization. Quantum, 3:163, July 2019. arXiv:1610.06546 [quant-ph]
Pith/arXiv arXiv 2019
-
[26]
Hubeny, Mukund Rangamani, and Tadashi Takayanagi
Veronika E. Hubeny, Mukund Rangamani, and Tadashi Takayanagi. A Covari- ant Holographic Entanglement Entropy Proposal.Journal of High Energy Physics, 2007(07):062–062, July 2007. arXiv:0705.0016 [hep-th]
Pith/arXiv arXiv 2007
-
[27]
Jayadev Acharya, Ibrahim Issa, Nirmal V. Shende, and Aaron B. Wagner. Measuring Quantum Entropy.IEEE Journal on Selected Areas in Information Theory, 1(2):454– 468, August 2020. arXiv:1711.00814 [quant-ph]
Pith/arXiv arXiv 2020
-
[28]
Daniel Jafferis, Alexander Zlokapa, Joseph D. Lykken, David K. Kolchmeyer, Saman- tha I. Davis, Nikolai Lauk, Hartmut Neven, and Maria Spiropulu. Traversable worm- hole dynamics on a quantum processor.Nature, 612(7938):51–55, December 2022. 17
work page 2022
-
[2022]
Publisher: American Physical Society
This paper was first reviewed by deepseek-v4-flash on August 5, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.