arxiv: 2604.26813 · v2 · submitted 2026-04-29 · 🪐 quant-ph

Recognition: no theorem link

Classical simulation of free-fermionic dynamics and quantum chemistry with magic input

Changhun Oh , Micha{\l} Oszmaniec , Oliver Reardon-Smith , Zolt\'an Zimbor\'as

Authors on Pith no claims yet

Pith reviewed 2026-05-13 07:19 UTC · model grok-4.3

classification 🪐 quant-ph

keywords free-fermionic dynamicsclassical simulationPfaffian polynomialnon-Gaussian statesquantum chemistrytransition amplitudespaired electronsWilson observables

0 comments

The pith

Block-product paired non-Gaussian fermionic states allow efficient classical approximation of transition amplitudes and overlaps under free dynamics via Pfaffian reduction.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that fermionic states with a specific block-product paired non-Gaussian structure admit efficient classical approximation of transition amplitudes, overlaps, and number correlators of any weight when evolved under free-fermionic dynamics. The key step is an algebraic compression that turns the full multiparticle interference pattern into the evaluation of one coefficient in a multivariate Pfaffian polynomial. A reader would care because these classical estimators reach the same additive error level as the statistical uncertainty of a quantum device that uses only a modest number of measurement shots, thereby supplying concrete benchmarks for experiments and for overlap subroutines in quantum chemistry.

Core claim

We prove that for block-product paired non-Gaussian fermionic states, essential quantum simulation primitives -- transition amplitudes, overlaps, and arbitrary-weight number correlators -- can be efficiently approximated to additive error under free-fermionic dynamics. This tractability stems from an algebraic reduction that compresses exponentially large multiparticle interference into a single coefficient of a multivariate Pfaffian polynomial. Because these classical estimators match the intrinsic O(1/sqrt(K)) statistical uncertainty of quantum hardware utilizing K measurement shots, they constitute a practical benchmark.

What carries the argument

Algebraic reduction that maps multiparticle interference to evaluation of a single coefficient in a multivariate Pfaffian polynomial, yielding Pfaffian-kernel estimators for the target quantities.

If this is right

Transition amplitudes between block-product paired states are classically approximable to additive error.
State overlaps for the same class of states are classically tractable.
Number correlators of arbitrary weight admit efficient classical estimators.
High-weight Wilson observables in noninteracting quenches of trapped-ion experiments can be classically benchmarked.
Overlap-based subroutines for antisymmetrized products of strongly orthogonal geminals in quantum chemistry become classically estimable.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The paired-electron scaffold itself is dequantized, so quantum advantage in these settings must rely on either unpaired magic or genuine interactions.
The same Pfaffian reduction may supply benchmarks for other paired fermionic platforms in condensed-matter simulation.
Classical estimators of this form could be combined with tensor-network methods to enlarge the simulable regime without increasing quantum resources.
Testing the reduction on small trapped-ion systems would immediately reveal whether real hardware noise already exceeds the additive-error threshold.

Load-bearing premise

The input states must be exactly block-product paired and non-Gaussian; any deviation from this paired structure, or any departure from purely free-fermionic dynamics, removes the Pfaffian reduction.

What would settle it

A measured transition amplitude or overlap for a concrete block-product paired non-Gaussian state evolved under free-fermionic dynamics that deviates from the corresponding Pfaffian-coefficient prediction by more than the claimed additive error, or that exceeds the O(1/sqrt(K)) quantum statistical limit.

Figures

Figures reproduced from arXiv: 2604.26813 by Changhun Oh, Micha{\l} Oszmaniec, Oliver Reardon-Smith, Zolt\'an Zimbor\'as.

**Figure 1.** Figure 1: FIG. 1 view at source ↗

**Figure 2.** Figure 2: FIG. 2. Fermionic sampling with view at source ↗

**Figure 3.** Figure 3: FIG. 3 view at source ↗

**Figure 4.** Figure 4: The detailed construction is given in App. E view at source ↗

**Figure 5.** Figure 5: summarizes the RMSE scaling with respect to the Monte Carlo sample count and the system size for several local uniformity profiles. To make the dependence on local uniformity more explicit, we also perform a direct sweep with uL = uR = u at fixed N and K, shown in view at source ↗

**Figure 6.** Figure 6: FIG. 6 view at source ↗

**Figure 7.** Figure 7: FIG. 7. Sample-complexity contours for real-time Hirsch view at source ↗

**Figure 8.** Figure 8: FIG. 8. Sample-complexity contours for extent-based view at source ↗

read the original abstract

Establishing the precise computational boundary between classically tractable fermionic systems and those capable of genuine quantum advantage is a central challenge in quantum simulation. While injecting non-Gaussian ``magic" inputs into free-fermion circuits is widely expected to generate intractable complexity, we identify a physically motivated intermediate regime. We prove that for block-product paired non-Gaussian fermionic states, essential quantum simulation primitives -- transition amplitudes, overlaps, and arbitrary-weight number correlators -- can be efficiently approximated to additive error under free-fermionic dynamics. This tractability stems from an algebraic reduction that compresses exponentially large multiparticle interference into a single coefficient of a multivariate Pfaffian polynomial. Because these classical estimators match the intrinsic $O(1/\sqrt{K})$ statistical uncertainty of quantum hardware utilizing $K$ measurement shots, they constitute a practical benchmark. Building on this foundation, we construct an additive-error estimator for high-weight Wilson observables in the noninteracting quench of recent trapped-ion experiments, providing a rigorous classical benchmark. Extending this to quantum chemistry, we demonstrate that core overlap-based subroutines for antisymmetrized products of strongly orthogonal geminals admit efficient additive-error Pfaffian-kernel estimators. Ultimately, these results sharpen the boundary of quantum advantage, establishing that the paired-electron scaffold is dequantized and clarifying where quantum resources are indispensable.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

They've identified a block-product paired regime of non-Gaussian fermionic states that stays classically tractable under free evolution by reducing key primitives to a single Pfaffian coefficient.

read the letter

The main point is that this paper carves out a concrete intermediate class of states—block-product paired non-Gaussian fermions—where injecting magic still leaves transition amplitudes, overlaps, and arbitrary-weight number correlators classically approximable to additive error after free-fermionic evolution. The reduction works by factoring the problem into a multivariate Pfaffian polynomial whose relevant coefficient can be pulled out efficiently because the block structure keeps the matrix size linear rather than exponential. That is the actual new piece relative to standard free-fermion and magic-state results. They then use it to build explicit additive-error estimators for Wilson observables in trapped-ion quenches and for overlap subroutines in antisymmetrized geminal quantum chemistry, which is useful because the error scale matches the 1/sqrt(K) shot noise on real hardware. Those are practical benchmarks rather than purely theoretical statements. The limitation is straightforward: the whole thing requires the states to be exactly block-product paired. Any deviation from that structure or any departure from purely free dynamics breaks the reduction, and the paper does not claim robustness beyond that. The algebraic steps look internally consistent from the stress-test description, but the additive-error bounds and any hidden polynomial factors would need checking in the full proofs. This is aimed at people mapping the boundary of quantum advantage in fermionic simulation and at quantum chemists who need reliable classical references for geminal overlaps. It deserves a serious referee because the claimed regime is new, the application to existing experiments is direct, and the core reduction is falsifiable algebra rather than fitting.

Referee Report

0 major / 2 minor

Summary. The manuscript proves that for block-product paired non-Gaussian fermionic states, transition amplitudes, overlaps, and arbitrary-weight number correlators can be efficiently approximated to additive error under free-fermionic dynamics by reducing the problem to extracting a single coefficient from a multivariate Pfaffian polynomial. It applies this to construct classical benchmarks for high-weight Wilson observables in trapped-ion experiments and for overlap-based subroutines in quantum chemistry using antisymmetrized products of strongly orthogonal geminals.

Significance. If the algebraic reduction holds, the result is significant because it dequantizes a physically motivated class of non-Gaussian states in free-fermion dynamics, supplying classical estimators whose additive error matches the intrinsic O(1/sqrt(K)) statistical uncertainty of K-shot quantum hardware. The work sharpens the boundary between classically tractable and quantum-advantaged fermionic simulation and supplies concrete benchmarks for both trapped-ion quenches and geminal-based quantum chemistry methods.

minor comments (2)

[Abstract] Abstract: the phrase 'multivariate Pfaffian polynomial' is used without a preceding definition or reference; a short parenthetical or footnote would improve immediate readability.
[Abstract] The error analysis for the additive approximation is asserted to be polynomial-time but the explicit dependence on the number of blocks and the Pfaffian degree is not summarized in the abstract; adding a one-sentence complexity statement would strengthen the claim.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of our manuscript, accurate summary of the algebraic reduction to Pfaffian polynomials, and recommendation for minor revision. The significance statement correctly highlights the dequantization of block-product paired states and the practical benchmarking value for trapped-ion and geminal-based quantum chemistry applications. No specific major comments were provided in the report, so we have no point-by-point rebuttals to address.

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The paper derives classical tractability of amplitudes, overlaps, and number correlators for block-product paired non-Gaussian fermionic states under free-fermionic dynamics via an algebraic reduction to extraction of a single coefficient from a multivariate Pfaffian polynomial. This follows directly from the block-product factorization (keeping matrix size linear in blocks) combined with the linear action of free-fermionic evolution on Majorana operators, both of which are standard algebraic facts independent of the claimed result. No load-bearing step reduces by construction to a fitted parameter, self-definition, or unverified self-citation chain; the additive-error bounds and polynomial-time claims are obtained from the Pfaffian properties themselves without circular renaming or smuggling of ansatzes. The derivation is therefore self-contained against external algebraic benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on algebraic identities of multivariate Pfaffians and the structural assumption of block-product pairing; no free parameters or new entities are introduced.

axioms (1)

standard math Multivariate Pfaffian polynomials admit efficient coefficient extraction that compresses multiparticle interference for paired states
Invoked to justify the classical estimator for amplitudes and correlators

pith-pipeline@v0.9.0 · 5550 in / 1246 out tokens · 64578 ms · 2026-05-13T07:19:57.718766+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Two Layers, No Swaps: Biplanar SPOQC Architecture Improves Runtime of Fermi-Hubbard Simulation
quant-ph 2026-05 unverdicted novelty 6.0

The biplanar architecture maps Fermi-Hubbard spin sectors to two planes, eliminating swaps and cutting each Trotter step depth to 4t_synth + 90 logical timesteps versus 6t_synth + 354 in single-plane methods, yielding...

Reference graph

Works this paper leans on

58 extracted references · 58 canonical work pages · cited by 1 Pith paper

[1]

The transition amplitude expands as a sum of Slater determinants: A= 2 −N/2 X i∈{0,1}N det USx,R(i) ,(18) whereS x ⊂[4N] is the set of occupied output modes in |x⟩

or|0011⟩(i t = 1), and letR(i)⊂[4N] denote the resulting set of occupied input modes (exactly two per block). The transition amplitude expands as a sum of Slater determinants: A= 2 −N/2 X i∈{0,1}N det USx,R(i) ,(18) whereS x ⊂[4N] is the set of occupied output modes in |x⟩. To circumvent the exponential cost of evaluating 2 N determinants, we reorganize t...

work page
[2]

Equivalently, defining ˜A≡2 −N/2 ˜S, we obtain| ˜A−A| ≤ ϵ/ √ 2N with probability at least 1−δ

Since|Z(b)| ≤1, applying Hoeffding’s inequality implies that it suffices to takeK=O(ϵ −2 logδ −1) to achieve additive errorϵwith failure probability at mostδ. Equivalently, defining ˜A≡2 −N/2 ˜S, we obtain| ˜A−A| ≤ ϵ/ √ 2N with probability at least 1−δ. Theorem 1(Additive-error transition amplitude for block-magic input).Let ˆUbe a passive FLO unitary on ...

work page
[3]

Exact” denotes the exact untrotterized free-fermion evolution, “FLO

Consider the transition amplitudeA= ⟨Ψ| ˆU|Ψ⟩. There exists a classical randomized algorithm that outputs an estimator ˜AsatisfyingPr h | ˜A−A|> ϵ i ≤ δ. The algorithm requires a sample complexityK= O(ϵ−2 logδ −1), with each sample processed inO(N 3) time. C. Extension to block-product APSG states While the|Ψ 4⟩⊗N block-magic state illustrates the core mi...

work page 2024
[4]

Arute, K

F. Arute, K. Arya, R. Babbush, D. Bacon, J. C. Bardin, R. Barends, A. Bengtsson, S. Boixo, M. Broughton, B. B. Buckley,et al., Observation of separated dynamics of charge and spin in the fermi-hubbard model, arXiv preprint arXiv:2010.07965 (2020)

work page arXiv 2010
[5]

A. J. Daley, I. Bloch, C. Kokail, S. Flannigan, N. Pearson, M. Troyer, and P. Zoller, Practical quantum advantage in quantum simulation, Nature607, 667 (2022)

work page 2022
[6]

F. Alam, J. L. Bosse, I. ˇCepait˙ e, A. Chapman, L. Clinton, M. Crichigno, E. Crosson, T. Cubitt, C. Derby, O. Dowinton,et al., Fermionic dynamics on a trapped- ion quantum computer beyond exact classical simulation, arXiv preprint arXiv:2510.26300 (2025)

work page arXiv 2025
[7]

Arute, K

F. Arute, K. Arya, R. Babbush, D. Bacon, J. C. Bardin, R. Barends, S. Boixo, M. Broughton, B. B. Buckley, et al., Hartree-fock on a superconducting qubit quantum computer, Science369, 1084 (2020)

work page 2020
[8]

Oszmaniec, N

M. Oszmaniec, N. Dangniam, M. E. Morales, and Z. Zimbor´ as, Fermion sampling: a robust quantum computational advantage scheme using fermionic linear optics and magic input states, PRX Quantum3, 020328 (2022)

work page 2022
[9]

B. M. Terhal and D. P. DiVincenzo, Classical simulation of noninteracting-fermion quantum circuits, Physical Review A65, 032325 (2002)

work page 2002
[10]

Jozsa and A

R. Jozsa and A. Miyake, Matchgates and classical simulation of quantum circuits, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences464, 3089 (2008)

work page 2008
[11]

L. G. Valiant, Quantum circuits that can be simulated classically in polynomial time, SIAM Journal on Computing31, 1229 (2002)

work page 2002
[12]

Fermionic Linear Optics and Matchgates

E. Knill, Fermionic linear optics and matchgates, arXiv preprint quant-ph/0108033 (2001)

work page Pith review arXiv 2001
[13]

Oszmaniec and Z

M. Oszmaniec and Z. Zimbor´ as, Universal extensions of restricted classes of quantum operations, arXiv preprint arXiv:1705.11188 (2017)

work page arXiv 2017
[14]

Hebenstreit, R

M. Hebenstreit, R. Jozsa, B. Kraus, and S. Strelchuk, Computational power of matchgates with supplementary resources, Physical Review A102, 052604 (2020)

work page 2020
[15]

Dias and R

B. Dias and R. Koenig, Classical simulation of non- gaussian fermionic circuits, Quantum8, 1350 (2024)

work page 2024
[16]

D. A. Ivanov, Computational complexity of exterior products and multi-particle amplitudes of non- interacting fermions in entangled states, arXiv preprint arXiv:1603.02724 (2016). 15

work page arXiv 2016
[17]

Aaronson and A

S. Aaronson and A. Arkhipov, The computational complexity of linear optics, inProceedings of the forty- third annual ACM symposium on Theory of computing (2011) pp. 333–342

work page 2011
[18]

D. J. Brod and E. F. Galvao, Extending matchgates into universal quantum computation, Physical Review A—Atomic, Molecular, and Optical Physics84, 022310 (2011)

work page 2011
[19]

McArdle, S

S. McArdle, S. Endo, A. Aspuru-Guzik, S. C. Benjamin, and X. Yuan, Quantum computational chemistry, Reviews of Modern Physics92, 015003 (2020)

work page 2020
[20]

W. J. Huggins, B. A. O’Gorman, N. C. Rubin, D. R. Reichman, R. Babbush, and J. Lee, Unbiasing fermionic quantum monte carlo with a quantum computer, Nature 603, 416 (2022)

work page 2022
[21]

P. J. Ollitrault, C. L. Cortes, J. F. Gonthier, R. M. Parrish, D. Rocca, G.-L. Anselmetti, M. Degroote, N. Moll, R. Santagati, and M. Streif, Enhancing initial state overlap through orbital optimization for faster molecular electronic ground-state energy estimation, Physical Review Letters133, 250601 (2024)

work page 2024
[22]

U. Baek, D. Hait, J. Shee, O. Leimkuhler, W. J. Huggins, T. F. Stetina, M. Head-Gordon, and K. B. Whaley, Say no to optimization: A nonorthogonal quantum eigensolver, PRX Quantum4, 030307 (2023)

work page 2023
[23]

A. L. Fetter and J. D. Walecka,Quantum theory of many- particle systems(Courier Corporation, 2012)

work page 2012
[24]

D. E. Knuth, Overlapping pfaffians, arXiv preprint math/9503234 (1995)

work page arXiv 1995
[25]

Wimmer, Algorithm 923: Efficient numerical computation of the pfaffian for dense and banded skew- symmetric matrices, ACM Transactions on Mathematical Software (TOMS)38, 1 (2012)

M. Wimmer, Algorithm 923: Efficient numerical computation of the pfaffian for dense and banded skew- symmetric matrices, ACM Transactions on Mathematical Software (TOMS)38, 1 (2012)

work page 2012
[26]

Gurvits, On the complexity of mixed discriminants and related problems, inInternational Symposium on Mathematical Foundations of Computer Science (Springer, 2005) pp

L. Gurvits, On the complexity of mixed discriminants and related problems, inInternational Symposium on Mathematical Foundations of Computer Science (Springer, 2005) pp. 447–458

work page 2005
[27]

Aaronson and T

S. Aaronson and T. Hance, Generalizing and derandomizing gurvits’s approximation algorithm for the permanent, arXiv preprint arXiv:1212.0025 (2012)

work page arXiv 2012
[28]

C. Oh, Y. Lim, Y. Wong, B. Fefferman, and L. Jiang, Quantum-inspired classical algorithms for molecular vibronic spectra, Nature Physics20, 225 (2024)

work page 2024
[29]

Ikai, On the theory of pfaffians based on exponential maps in exterior algebras, Linear algebra and its applications434, 1094 (2011)

H. Ikai, On the theory of pfaffians based on exponential maps in exterior algebras, Linear algebra and its applications434, 1094 (2011)

work page 2011
[30]

Reardon-Smith, M

O. Reardon-Smith, M. Oszmaniec, and K. Korzekwa, Improved simulation of quantum circuits dominated by free fermionic operations, Quantum8, 1549 (2024)

work page 2024
[31]

Cudby and S

J. Cudby and S. Strelchuk, Gaussian decomposition of magic states for matchgate computations (2024), arXiv:2307.12654 [quant-ph]

work page arXiv 2024
[32]

Reardon-Smith, The fermionic linear optical extent is multiplicative for 4 qubit parity eigenstates (2024), arXiv:2407.20934 [quant-ph]

O. Reardon-Smith, The fermionic linear optical extent is multiplicative for 4 qubit parity eigenstates (2024), arXiv:2407.20934 [quant-ph]

work page arXiv 2024
[33]

Bravyi, D

S. Bravyi, D. Browne, P. Calpin, E. Campbell, D. Gosset, and M. Howard, Simulation of quantum circuits by low- rank stabilizer decompositions, Quantum3, 181 (2019)

work page 2019
[34]

J. R. McClean, M. E. Kimchi-Schwartz, J. Carter, and W. A. De Jong, Hybrid quantum-classical hierarchy for mitigation of decoherence and determination of excited states, Physical Review A95, 042308 (2017)

work page 2017
[35]

J. F. Stanton and R. J. Bartlett, The equation of motion coupled-cluster method. a systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties, The Journal of chemical physics98, 7029 (1993)

work page 1993
[36]

E. J. Sundstrom and M. Head-Gordon, Non-orthogonal configuration interaction for the calculation of multielectron excited states, The Journal of chemical physics140(2014)

work page 2014
[37]

W. J. Huggins, J. Lee, U. Baek, B. O’Gorman, and K. B. Whaley, A non-orthogonal variational quantum eigensolver, New Journal of Physics22, 073009 (2020)

work page 2020
[38]

Miller, J

A. Miller, J. Favre, Z. Holmes, ¨O. Salehi, R. Chakraborty, A. Nyk¨ anen, Z. Zimbor´ as, A. Glos, and G. Garc´ ıa- P´ erez, Simulation of fermionic circuits using majorana propagation, arXiv preprint arXiv:2503.18939 (2025)

work page arXiv 2025
[39]

D’Anna, J

M. D’Anna, J. Nys, and J. Carrasquilla, Majorana string simulation of nonequilibrium dynamics in two- dimensional lattice fermion systems, arXiv preprint arXiv:2511.02809 (2025)

work page arXiv 2025
[40]

Seron, L

B. Seron, L. Novo, A. Arkhipov, and N. J. Cerf, Efficient validation of boson sampling from binned photon-number distributions, Quantum8, 1479 (2024)

work page 2024
[41]

Lim and C

Y. Lim and C. Oh, Efficient classical algorithms for linear optical circuits, arXiv preprint arXiv:2502.12882 (2025)

work page arXiv 2025
[42]

Bak´ o, Z

B. Bak´ o, Z. Kolarovszki, and Z. Zimbor´ as, Fermionic born machines: Classical training of quantum generative models based on fermion sampling, arXiv preprint arXiv:2511.13844 (2025)

work page arXiv 2025
[43]

M. S. Rudolph, S. Lerch, S. Thanasilp, O. Kiss, O. Shaya, S. Vallecorsa, M. Grossi, and Z. Holmes, Trainability barriers and opportunities in quantum generative modeling, npj Quantum Information10, 116 (2024)

work page 2024
[44]

Recio-Armengol, S

E. Recio-Armengol, S. Ahmed, and J. Bowles, Train on classical, deploy on quantum: scaling generative quantum machine learning to a thousand qubits, arXiv preprint arXiv:2503.02934 (2025)

work page arXiv 2025
[45]

Tecmer and K

P. Tecmer and K. Boguslawski, Geminal-based electronic structure methods in quantum chemistry. toward a geminal model chemistry, Physical Chemistry Chemical Physics24, 23026 (2022)

work page 2022
[46]

P. R. Surj´ an, ´A. Szabados, P. Jeszenszki, and T. Zoboki, Strongly orthogonal geminals: size-extensive and variational reference states, Journal of Mathematical Chemistry50, 534 (2012)

work page 2012
[47]

L. Zhao, J. Goings, K. Shin, W. Kyoung, J. I. Fuks, J.-K. Kevin Rhee, Y. M. Rhee, K. Wright, J. Nguyen, J. Kim,et al., Orbital-optimized pair-correlated electron simulations on trapped-ion quantum computers, npj Quantum Information9, 60 (2023)

work page 2023
[48]

L. Zhao, Q. Wang, J. J. Goings, K. Shin, W. Kyoung, S. Noh, Y. M. Rhee, and K. Kim, Enhancing the electron pair approximation with measurements on trapped-ion quantum computers, npj Quantum Information10, 76 (2024)

work page 2024
[49]

Helgaker, P

T. Helgaker, P. Jorgensen, and J. Olsen,Molecular electronic-structure theory(John Wiley & Sons, 2013)

work page 2013
[50]

D. A. Mazziotti, Two-electron reduced density matrix as the basic variable in many-electron quantum chemistry and physics, Chemical Reviews112, 244 (2012)

work page 2012
[51]

Motta, J

M. Motta, J. Shee, S. Zhang, and G. K.-L. Chan, Efficient ab initio auxiliary-field quantum monte carlo calculations in gaussian bases via low-rank tensor decomposition, 16 Journal of chemical theory and computation15, 3510 (2019)

work page 2019
[52]

Leimkuhler and K

O. Leimkuhler and K. B. Whaley, Exponential quantum speedups in quantum chemistry with linear depth, arXiv preprint arXiv:2503.21041 (2025)

work page arXiv 2025
[53]

L. G. Valiant, The complexity of computing the permanent, Theoretical computer science8, 189 (1979)

work page 1979
[54]

Ikenmeyer and M

C. Ikenmeyer and M. Walter, Hyperpfaffians and geometric complexity theory, arXiv preprint arXiv:1912.09389 (2019)

work page arXiv 1912
[55]

J. E. Hirsch, Discrete hubbard-stratonovich transformation for fermion lattice models, Physical Review B28, 4059 (1983)

work page 1983
[56]

Dias and R

B. Dias and R. Koenig, Classical simulation of non- gaussian fermionic circuits, Quantum8, 1350 (2024). Appendix A: Proof of the Pfaffian-Wedge Identity In this Appendix, we provide an explicit derivation of the fundamental algebraic identity relating the Pfaffian of a skew-symmetric matrixA∈C 2N×2N to theN-th exterior power of its associated two- form fo...

work page 2024
[57]

Because the HS transformation requires complex coupling parameters (λ=a+ibwitha= Reλ), the one-body propagators ˆG(σ) are inherently non-unitary

A conservative worst-case envelope for the Pfaffian estimator The fundamental obstruction to pushing this classical simulation to arbitrary times is not algorithmic correctness, but variance. Because the HS transformation requires complex coupling parameters (λ=a+ibwitha= Reλ), the one-body propagators ˆG(σ) are inherently non-unitary. SinceV(σ ℓ) is diag...

work page
[58]

Extent-based simulation of interacting dynamics An alternative approach to simulating interacting dynamics is to generalize the extent-based simulation algorithms by combining them with our mixed-Pfaffian estimators. Assume we are interested in computing a quantity of the form ⟨Φ| ˆU|Ψ⟩,(F10) 21 1 2 3 4 5 6 Interaction Strength ( W ) 0.5 1.0 1.5 2.0 2.5 E...

work page