Pith. sign in

REVIEW 3 major objections 6 minor 5 cited by

Majorana string simulation of nonequilibrium dynamics in two-dimensional lattice fermion systems

T0 review · 3 major / 6 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper reports a Majorana-string propagation algorithm with a Trotter-consistent truncation that reaches accurate real-time Fermi-Hubbard dynamics in two dimensions, matching a cold-atom experiment on magnetic polaron formation and exte

desk verdict The method is promising and the numerics are largely clean, but the 'preserves Trotter accuracy' claim is not established: the per-layer truncation argument does not control cross-layer accumulation. read the letter →

arxiv 2511.02809 v2 pith:NHGGWSVV submitted 2025-11-04 cond-mat.quant-gas cond-mat.str-elquant-ph

Matteo D'Anna , Jannes Nys , Juan Carrasquilla This is my paper
classification cond-mat.quant-gascond-mat.str-elquant-ph
keywords MajoranastringsFermi-Hubbardmodelreal-timedynamicsHeisenbergpictureTrottererrortruncationmagneticpolaroncold-atomquantumsimulator
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper introduces Majorana propagation, a Heisenberg-picture algorithm that expands a fermionic observable as a sum of Majorana strings and updates that expansion gate by gate under a Trotterized time evolution. For quadratic Hamiltonians the expansion stays inside a fixed low-weight sector, so two-dimensional free-fermion dynamics is reproduced exactly. For interacting Fermi-Hubbard models the algorithm truncates strings by small coefficients and by the number of unpaired Majoranas, using a cap chosen to stay consistent with the Trotter error; on a 100-site one-dimensional quench it remains accurate to timescales beyond the tensor-network benchmarks used, and on a 7×7 lattice at strong coupling it yields hole dynamics in qualitative agreement with a cold-atom experiment. The payoff is a classically computable, systematically improvable route to short- and intermediate-time fermionic dynamics in two dimensions.

What carries the argument

The central object is the Majorana string µ(v), a Hermitian product of Majorana operators indexed by a binary vector v, whose 'unpaired count' w_s(v) counts modes where exactly one of the pair γ_i, γ'_i appears. The carrying mechanism is the splitting rule for conjugation by a Trotter rotation exp(−iθµ(v')/2): a string either is unchanged or splits into cosine and sine terms µ(v) and µ(v+v'), so each gate branches the expansion. The truncation then caps w_s at S after each layer while allowing temporary growth to S+p inside the layer—S+2 for the second-order schemes used here—and drops strings with |λ|<ε. For quadratic Hamiltonians, Appendix D shows the weight w(v) is conserved under hopping

What would settle it

Using the paper's own 3×3 benchmark (U/t=8, checkerboard initial state with central hole, exact-diagonalization reference), run Majorana propagation with S=4 and ε=10⁻⁶, comparing an in-layer cap of S+2 against S+4—the paper's own 2p bound for a second-order scheme. If the two hole-probability curves agree to within the Trotter error up to τ/t≈2, the consistency claim is supported; if they diverge, the discarded strings were recombining into relevant low-unpaired strings and the truncation is not Trotter-consistent.

Watch

Extended reading notes

Core claim

The central claim is that controlled truncation in a Majorana-string basis—removing strings with coefficient below ε and strings whose unpaired-Majorana count w_s exceeds a cutoff S—can be made 'Trotter-consistent,' so that the discarded strings do not change the observable at the accuracy of the underlying second-order Trotter decomposition, provided the in-layer cap is S+2 for the second-order scheme used. On this basis the paper reports converged local densities in a 100-site 1D Fermi-Hubbard quench for S≥6 up to τ/t≈4–6, beyond the stable range of the MPS benchmarks it compares with; convergence to exact diagonalization on 3×3 lattices; and qualitative agreement with a cold-atom experime

Load-bearing premise

The scheme's accuracy rests on the Trotter-consistency claim that Majorana strings whose unpaired count temporarily exceeds the in-layer cap cannot recombine into allowed strings within the same Trotter step; if that bound is too tight, the truncation silently lowers the accuracy below the advertised order.

Editorial extensions

If this is right

  • For quadratic Hamiltonians, local observables evolve exactly within a fixed low-weight Majorana-string sector, making two-dimensional free-fermion scattering tractable without truncation, a regime where tensor-network methods struggle.
  • For interacting Fermi-Hubbard models, convergence in S and ε yields reliable local observables at intermediate times (1D: τ/t≈4–6 at 100 sites) beyond the stable range of the MPS benchmarks used in the paper.
  • The scheme accepts variational initial states (e.g., DMRG ground states) by evaluating ⟨ψ|µ(v)|ψ⟩ as an MPO, extending the method beyond simple Fock-state initial conditions.
  • On the 7×7 magnetic-polaron setup at U/t=8.72, hole dynamics agree qualitatively with a cold-atom experiment, while Lieb-Robinson commutators show observable spreading across the full lattice by τ/t<1, quantifying finite-size effects in such comparisons.
  • Coefficient truncation ε directly controls the tails of out-of-time-order correlation profiles, giving a tunable trade-off between computational cost and light-cone fidelity.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper leaves implicit that its Trotter-consistent in-layer cap could be imported into Pauli propagation for qubit circuits, upgrading ad hoc weight truncations into principled cuts that preserve product-formula accuracy.
  • The overlap distributions in Appendix C suggest a natural extension: design truncation rules that use the actual initial-state overlap weight rather than only w_s, which could extend Majorana propagation to highly entangled initial states.
  • The 19×19 four-hole commutator data imply that the embedded 7×7 regions in the experiment begin to interact by τ/t≈0.4, so quantitative experimental comparisons may require simulating the full surrounding region or treating region observables as non-commuting.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper introduces a Heisenberg-picture algorithm, Majorana propagation (MP), for real-time dynamics of lattice fermions. Observables are expanded in a Majorana-string basis and evolved through Trotterized gates via exact splitting rules. To control operator growth, two truncations are used: a coefficient cutoff ε and a cutoff S on the number of unpaired Majoranas per string. The authors claim that the latter is 'Trotter-consistent', i.e. it preserves the Trotter accuracy of the underlying product formula. The method is shown to be exact for quadratic Hamiltonians by a weight-conservation argument, and is benchmarked on free-fermion scattering, 1D Fermi-Hubbard quenches against MPS, 2D 3×3 ED convergence tests, and a 7×7 polaron quench against the experiment of Ref. [10]. The central numerical claim is that controlled truncation reaches reliable timescales comparable to or beyond variational tensor-network and experimental results.

Significance. If the central claims hold, this is a useful non-variational simulator for two-dimensional fermionic dynamics at short to intermediate times. The exact Gaussian limit is a genuine strength: Appendix D proves weight conservation and the 3×3 exact-diagonalization convergence studies (Figs. 8, 9) are clean and support the practical usefulness of the cutoffs. The paper also releases a Julia implementation, which is a concrete contribution. The comparison against the 7×7 experiment and the Lieb-Robinson diagnostics are valuable. However, the abstract and Section I make an explicit theoretical claim that the truncation scheme 'preserves Trotter accuracy'; as argued below, this claim is not established by the presented per-layer argument and needs either a proof or a significant softening plus numerical evidence. The empirical benchmarks are promising but do not by themselves justify the Trotter-consistency statement.

major comments (3)
  1. [§II.C, Eq. (19) and the three bullets after Eq. (22)] The 'Trotter-consistent' truncation argument is per-layer, but the algorithm applies the final cap S after every Trotter layer. A string with ws>S that is discarded after layer m can, in later layers, be reduced by hopping gates to ws≤S and ultimately to ws=0, thereby contributing to the final Fock-state expectation. The Campbell expansion (19) only controls contributions generated inside one layer from retained parent strings; it does not bound the influence of a discarded string over subsequent layers. This is the load-bearing gap for the abstract claim. The authors should either provide a multi-layer error bound (e.g. showing that any discarded string can only affect low-ws components at a specified order in δτ) or report a numerical δτ-scaling test at fixed S and ε that demonstrates global second-order convergence of the truncated algorithm. The reader's specific factor-of-two arithm
  2. [§II.C, Eq. (22); §III.B, Eq. (25)] The temporary cap S'=S+p is derived from a single Campbell expansion for the whole Trotter layer, but the code applies a gate-by-gate product of local hoppings (with the third factor applied in reverse order). Within one half-layer, a string can undergo many intermediate hopping gates and temporarily exceed S+p before returning to ws≤S at the end of the layer; such paths are not captured by the p-commutator bound. The 'temporary growth' rule therefore needs to be justified for the actual sequential gate implementation, not for an idealized layer-level expansion. Otherwise the cap may discard strings that the per-layer argument intended to keep, or keep strings that cannot be justified in either direction.
  3. [§II.C and Figs. 2(a), 5(a)] The coefficient truncation ε is enforced 'after each gate application' and is not covered by the Trotter-consistency argument for the S cutoff. The numerical results show visible ε-dependence at late times, especially in Fig. 5(a) for S=4,6 at τ/t≳1. If the phrase 'preserves Trotter accuracy' is meant to include ε truncation, an error estimate is needed; if not, the claim should be restricted to the S truncation and ε should be described as a separate heuristic control parameter. As written, the abstract conflates the two truncations.
minor comments (6)
  1. [Abstract / §III.C] Typo: 'projected entanged pair states' should be 'entangled'. Please correct.
  2. [Fig. 3] The labels 'ε=10□9' and 'ε=10□5' appear to be broken LaTeX; the exponent symbols are missing.
  3. [§III.B, Fig. 2(b)] The right axis showing 'maximum bond dimension' is confusingly placed; clarify whether the dashed grey line refers to χmax over all times or the final bond dimension.
  4. [§II.C] The temporary cap S' is introduced verbally but not defined with an equation. A boxed definition, e.g. S'=S+p for a p-th order layer, would help reproducibility.
  5. [§III.B, Eq. (25)] The statement 'gates in the third factor are applied in reverse order' is important but not standard terminology; briefly explain why reverse order is used (symmetric Trotterization) so readers can reproduce the exact gate sequence.
  6. [Appendix D] The proof that any weight-w Majorana string stays in the weight-w subspace for Gaussian dynamics is clear. It would be helpful to state explicitly that this holds for the unpaired count ws as well, or explain why weight conservation is the relevant invariant for the free-fermion benchmarks.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: algorithm parameters are convergence controls and benchmarks are external.

full rationale

The paper's derivation chain is self-contained at the level that matters for circularity. The Majorana-string splitting rules follow from the anticommutation relations and the phase calculus in the appendices; the Gaussian weight-conservation result is proven directly for quadratic Hamiltonians. The truncation cutoffs S and epsilon are convergence/truncation controls, not parameters fitted to the target data: the text explicitly optimizes S until observables converge and checks convergence against exact diagonalization on 3x3 lattices, MPS/fPEPS calculations, and the independent experimental data of Ref. [10]. No 'prediction' is obtained by renaming a fitted input. The self-citations present (e.g., Ref. [4] on neural-network variational methods, Ref. [20] on fermion-to-qudit mappings, Ref. [48] NetKet, Ref. [49] code release) are contextual or implementational and are not load-bearing for the central claim. The manuscript also openly acknowledges limitations, such as leaving dedicated truncation rules for generic ground states to future work and noting finite-size effects in the experimental comparison; these are correctness or scope caveats, not circular reasoning. The skeptic's concern about the Trotter-consistent truncation bound is a numerical-analysis gap rather than a circularity, since the truncation is not defined in terms of the benchmark results and the claimed accuracy is tested by convergence and external comparison. Therefore no circular step can be exhibited from the paper's own equations.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The paper introduces no new physical entities. The main intellectual burden is carried by two hand-set truncation parameters, S and ε, plus the heuristic that truncation error can be made comparable to Trotter error. The 'Trotter-consistent' rule is the least externally justified item in this ledger.

free parameters (3)
  • S (unpaired-Majorana cutoff) = 2–10 depending on simulation
    Chosen by hand; controls the truncation of strings with many unpaired Majoranas. The central interacting results use S=4–8.
  • ε (coefficient truncation threshold) = 10^-4 to 10^-9 depending on simulation
    Chosen by hand; strings with |λ_v|<ε are discarded. Affects accuracy and computational cost, especially at longer times.
  • δτ (Trotter time step) = 0.01 for Gaussian dynamics, 0.12 for 1D interacting; not explicitly stated for 2D
    Discretization parameter controlling Trotter error; chosen by hand.
assumptions (6)
  • standard math Majorana operators satisfy anticommutation relations and Majorana strings form a closed multiplication algebra with the phase factor of Eq. (7).
    Used throughout §II.A and Eq. (13); this is the algebraic foundation of the method.
  • standard math Trotter-Suzuki decomposition obeys U(δτ)=e^{-iHδτ}+O((δτ)^{p+1}), enabling the nested-commutator expansion of Campbell's identity.
    Invoked in §II.C to define 'Trotter-consistent' truncation; standard numerical-analysis background.
  • standard math The largest change in unpaired Majorana count from k nested commutators with hopping terms is at most 2k, with Δw_s∈{+2,0,-2}.
    Derived from weight-2 hopping strings; used to set the temporary cutoff S'. The factor-of-two bookkeeping is where the paper's argument appears inconsistent.
  • ad hoc to paper Strings with large unpaired-Majorana count contribute negligibly to the overlap with Fock-near initial states, so truncating them is safe.
    Motivates the S cutoff in §II.C; no a priori error bound is given, and Fig. 7 shows that generic DMRG ground states develop non-negligible higher-w_s overlaps.
  • ad hoc to paper Within a Trotter layer, temporarily relaxing the cap to S'=S+(2p)/2 preserves all contributions that can recombine to w_s≤S at order (δτ)^p.
    This is the 'Trotter-consistent cutting' premise in §II.C. Its arithmetic appears to permit only half the temporary growth that the paper's own bound allows.
  • domain assumption The relevant observables are spatially local, parity-even Majorana strings, and initial states are Fock-like checkerboard or DMRG-approximated states.
    The method is restricted to local observables in this paper, and the truncation rules are designed for initial states close to pure Fock states.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Majorana string simulation of nonequilibrium dynamics in two-dimensional lattice fermion systems." pith.science (2026). https://pith.science/paper/NHGGWSVV

@misc{pith2026251102809,
  author       = {Pith},
  title        = {Pith review of: Majorana string simulation of nonequilibrium dynamics in two-dimensional lattice fermion systems},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NHGGWSVV}},
  note         = {Machine review of arXiv:2511.02809}
}
read the original abstract

The study of real-time dynamics of fermions remains one of the last frontiers beyond the reach of classical simulations and is key to our understanding of quantum behavior in chemistry and materials, with implications for quantum technology. Here we introduce a Heisenberg-picture algorithm that propagates observables expressed in a Majorana-string basis using a truncation scheme that preserves Trotter accuracy and aims at maintaining computational efficiency. The framework is exact for quadratic Hamiltonians -- remaining restricted to a fixed low-weight sector determined by the physical observable -- admits variational initial states, and can be extended to interacting regimes via systematically controlled truncations. We benchmark our approach on one- and two-dimensional Fermi-Hubbard quenches, comparing against tensor network methods (MPS and fPEPS) and recent experimental data. The method achieves high accuracy on timescales comparable to state-of-the-art variational techniques and experiments, demonstrating that controlled Majorana-string truncation is a practical tool for simulating two-dimensional fermionic dynamics.

Figures

Figures reproduced from arXiv: 2511.02809 by the authors.

Figure 1
Figure 1. FIG. 1. Scattering of two spinless fermions on a 12 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. a) Expectation value of the density [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Space-time map of the OTOC [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Expectation values of the doublon occupation probability [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. a) Comparison of [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Expansion of [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Distribution of the overlaps against initial states prepared with a) [PITH_FULL_IMAGE:figures/full_fig_p011_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Dynamics of a hole in an anti-ferromagnetic background subject to [PITH_FULL_IMAGE:figures/full_fig_p012_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. The hole probability of a lattice site in the corner of the lattice, for the same dynamics as Fig. [PITH_FULL_IMAGE:figures/full_fig_p012_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10. Expansion of [PITH_FULL_IMAGE:figures/full_fig_p012_10.png]

Discussion (0). Sign in to comment.

Forward citations

Cited by 5 Pith papers

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

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

    quant-ph 2026-04 unverdicted novelty 7.0 of 10

    Block-product paired non-Gaussian fermionic states allow efficient classical additive-error approximation of transition amplitudes, overlaps, and high-weight correlators under free-fermionic dynamics using multivariat...

  2. Thermalization Dynamics in the Two-Dimensional Hubbard Model with Neural-Network Quantum States

    cond-mat.str-el 2026-06 unverdicted novelty 6.0 of 10

    Real-time dynamics in the 2D Hubbard model show thermalization of double occupancy below a critical U_c but clear breakdown of thermalization above it.

  3. Thermalization Dynamics in the Two-Dimensional Hubbard Model with Neural-Network Quantum States

    cond-mat.str-el 2026-06 conditional novelty 6.0 of 10

    In the 2D Hubbard model, the long-time double occupancy after a ramp quench matches the canonical thermal value for U≤3 but deviates above U≈3–4, suggesting thermalization breakdown.

  4. Combining non-parametric quantum states and MERA tensor networks for ground-state optimization

    quant-ph 2026-05 unverdicted novelty 6.0 of 10

    A hybrid method uses fixed quantum annealing states as boundary resources for classical MERA tensor networks to improve ground-state approximations without deeper quantum circuits.

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

    quant-ph 2026-04 unverdicted novelty 6.0 of 10

    Paired non-Gaussian fermionic states under free-fermionic dynamics admit efficient classical additive-error approximations for amplitudes, overlaps, and high-weight correlators via reduction to multivariate Pfaffian c...

Reference graph

Works this paper leans on

55 extracted references · 5 linked inside Pith · cited by 3 Pith papers

  1. [10]

    J. Nys, G. Pescia, A. Sinibaldi, and G. Carleo, Ab-initio variational wave functions for the time-dependent many-electron schr¨ odinger equation, Nature communications15, 9404 (2024)

    2024

  2. [1]

    (12) has|λ (m) v |< ε, for a fixedε, we remove µ(v) from the linear combination

    A first truncation is coefficient truncation: when- ever a stringµ(v) in the linear combination in Eq. (12) has|λ (m) v |< ε, for a fixedε, we remove µ(v) from the linear combination

  3. [2]

    The second truncation is based on the overlap be- tween Majorana strings and initial Fock states (16): we introduce a cutoffS∈N + on the number of un- paired Majoranasw s(v) in Eq. (4). To motivate the second rule, we observe that only Ma- jorana stringsµ(v) with fully paired Majoranas have a non-vanishing overlap with Fock states |n1 · · ·nN ⟩:= f † N nN...

  4. [3]

    As we propagate further towards the initial state, it becomes increasingly unlikely that stringsvwith a large ws(v) will contribute significantly to the final expectation value. A simple truncation rule based on estimations of the potential overlap with the initial state is not always available for a general initial stateρ(possibly beyond a pure Fock stat...

  5. [4]

    closeness condition

    Multiplicative factors We give the explicit form of the prefactorζ=ζ(v, v ′)∈ {±1,±i}appering in the “closeness condition” (7). We have [35] ζ(v, v′) = (−1)vT ωLv′+f(v,v ′)ivT ωv ′ ≡(−1) g(v,v ′)ivT ωv ′ (A7) for f(v, v′) = vT ωLv v′T ωLv′ +v T ωv ′ vT ωLv+v ′T ωLv′ + 1 , (A8) where all operations in (A7), (A8) are again to be understood as mod 2

  6. [5]

    It is therefore worth investigating how to write operator expressions, e.g

    Commutation relations with Majorana binary vectors Majorana strings are uniquely described in terms of their binary vectorv∈ {0,1} 2N , see (3). It is therefore worth investigating how to write operator expressions, e.g. the commutation relations, in terms of operations on binary vectors. In particular, we are interested in expressions to evaluatev T ωu, ...

  7. [6]

    With the definition of the Majorana strings (3), the expectation value ofµ(v) wrt to |n1 · · ·nN ⟩is computed as ⟨n1 · · ·nN |µ(v)|n 1 · · ·nN ⟩= = (i)vT ωLv⟨0|f n1 1

    Overlap with F ock basis states Consider the Fock basis states |n1 · · ·nN ⟩= f † N nN · · · f † 1 n1 |0⟩ for|0⟩the fermionic vacuum. With the definition of the Majorana strings (3), the expectation value ofµ(v) wrt to |n1 · · ·nN ⟩is computed as ⟨n1 · · ·nN |µ(v)|n 1 · · ·nN ⟩= = (i)vT ωLv⟨0|f n1 1 . . . fnN N γv1 1 (γ′ 1)v2 · · ·γv2N−1 N (γ′ N )v2N f † ...

  8. [7]

    Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics326, 96 (2011)

    U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics326, 96 (2011)

    2011

Show all 55 references
  1. [8]

    Georges, G

    A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions, Reviews of Modern Physics68, 13 (1996)

    1996

  2. [9]

    Troyer and U.-J

    M. Troyer and U.-J. Wiese, Computational complexity and fundamental limitations to fermionic quantum monte carlo simulations, Physical Review Letters94, 170201 (2005)

    2005

  3. [11]

    Fauseweh, Quantum many-body simulations on digital quantum computers: State-of-the-art and future challenges, Nature Communications15, 2123 (2024)

    B. Fauseweh, Quantum many-body simulations on digital quantum computers: State-of-the-art and future challenges, Nature Communications15, 2123 (2024)

    2024

  4. [12]

    Preskill, Quantum computing and the entanglement frontier (2012), arXiv:1203.5813 [quant-ph]

    J. Preskill, Quantum computing and the entanglement frontier (2012), arXiv:1203.5813 [quant-ph]

    2012 arXiv

  5. [13]

    Hirthe, T

    S. Hirthe, T. Chalopin, D. Bourgund, P. Bojovi´ c, A. Bohrdt, E. Demler, F. Grusdt, I. Bloch, and T. A. Hilker, Magnetically mediated hole pairing in fermionic ladders of ultracold atoms, Nature613, 463 (2023)

    2023

  6. [14]

    Hartke, B

    T. Hartke, B. Oreg, C. Turnbaugh, N. Jia, and M. Zwierlein, Direct observation of nonlocal fermion pairing in an attractive fermi–hubbard gas, Science381, 82 (2023)

    2023

  7. [15]

    Mazurenko, C

    A. Mazurenko, C. S. Chiu, G. Ji, M. F. Parsons, M. Kan´ asz-Nagy, R. Schmidt, F. Grusdt, E. Demler, D. Greif, and M. Greiner, A cold-atom fermi–hubbard antiferromagnet, Nature545, 462 (2017)

    2017

  8. [16]

    G. Ji, M. Xu, L. H. Kendrick, C. S. Chiu, J. C. Br¨ uggenj¨ urgen, D. Greif, A. Bohrdt, F. Grusdt, E. Demler, M. Lebrat, and M. Greiner, Coupling a mobile hole to an antiferromagnetic spin background: Transient dynamics of a magnetic polaron, Physical Review X11, 021022 (2021)

    2021

  9. [17]

    Lebrat, M

    M. Lebrat, M. Xu, L. H. Kendrick, A. Kale, Y. Gang, P. Seetharaman, I. Morera, E. Khatami, E. Demler, and M. Greiner, Observation of nagaoka polarons in a fermi–hubbard quantum simulator, Nature629, 317 (2024)

    2024

  10. [18]

    M. Xu, L. H. Kendrick, A. Kale, Y. Gang, A. W. Young, M. Lebrat, M. Greiner, C. Feng, and S. Zhang, A neutral-atom hubbard quantum simulator in the cryogenic regime, Nature 10.1038/s41586-025-09112-w (2025)

    2025 doi

  11. [19]

    Jordan and E

    P. Jordan and E. Wigner, ¨Uber das paulische ¨ aquivalenzverbot, Zeitschrift f¨ ur Physik47, 631 (1928)

    1928

  12. [20]

    S. B. Bravyi and A. Y. Kitaev, Fermionic quantum computation, Annals of Physics298, 210 (2002), arXiv:quant- ph/0003137

    2002

  13. [21]

    J. T. Seeley, M. J. Richard, and P. J. Love, The bravyi–kitaev transformation for quantum computation of electronic structure, The Journal of Chemical Physics137, 224109 (2012)

    2012

  14. [22]

    Barends, L

    R. Barends, L. Lamata, J. Kelly, L. Garc ´ ıa- ´Alvarez, A. G. Fowler,et al., Digital quantum simulation of fermionic and bosonic models in a superconducting circuit, Nature Communications6, 7654 (2015). 16

    2015

  15. [23]

    Wecker, M

    D. Wecker, M. B. Hastings, N. Wiebe, B. K. Clark, C. Nayak, and M. Troyer, Solving strongly correlated electron models on a quantum computer, Physical Review A92, 062318 (2015)

    2015

  16. [24]

    E. T. Campbell, Random compiler for fast hamiltonian simulation, Physical Review Letters123, 070503 (2019)

    2019

  17. [25]

    C. Cade, L. Mineh, A. Montanaro, and S. Stanisic, Strategies for solving the fermi-hubbard model on near-term quantum computers, Physical Review B102, 235122 (2020)

    2020

  18. [26]

    Carobene, A

    R. Carobene, A. Giachero, S. Barison, and J. Nys, Local fermion-to-qudit mappings: A practical recipe for four-level systems, Physical Review A112, 032619 (2025)

    2025

  19. [27]

    A. H. Karamlou, I. T. Rosen, S. E. Muschinske, C. N. Barrett, A. Di Paolo, L. Ding, P. M. Harrington, M. Hays, R. Das, D. K. Kim, B. M. Niedzielski, M. Schuldt, K. Serniak, M. E. Schwartz, J. L. Yoder, S. Gustavsson, Y. Yanay, J. A. Grover, and W. D. Oliver, Probing entangleme...

    2024

  20. [28]

    W. Chen, S. Zhang, J. Zhang, X. Su, Y. Lu, K. Zhang, M. Qiao, Y. Li, J.-N. Zhang, and K. Kim, Error-mitigated quantum simulation of interacting fermions with trapped ions, npj Quantum Information9, 10.1038/s41534-023-00784-8 (2023)

    2023 doi

  21. [29]

    Srinivasan, A

    D. Srinivasan, A. Beyer, D. Zhu, S. Churchill, K. Mehta, S. K. Sridhar, K. Chakrabarti, D. W. Steuerman, N. Chopra, and A. Dutt, Trapped-ion quantum simulation of the fermi–hubbard model as a lattice gauge theory using hardware-aware native gates (2024), arXiv preprint, arXiv:...

    2024

  22. [30]

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

    2002

  23. [31]

    Bravyi, Lagrangian representation for fermionic linear optics, Quantum Information & Computation5, 216 (2005), arXiv:quant-ph/0404180

    S. Bravyi, Lagrangian representation for fermionic linear optics, Quantum Information & Computation5, 216 (2005), arXiv:quant-ph/0404180

    2005 arXiv

  24. [32]

    S. B. Bravyi and A. Y. Kitaev, Fermionic quantum computation, Annals of Physics298, 210 (2002)

    2002

  25. [33]

    P. Rall, D. Liang, J. Cook, and W. Kretschmer, Simulation of qubit quantum circuits via pauli propagation, Physical Review A99, 062337 (2019)

    2019

  26. [34]

    Beguˇ si´ c, J

    T. Beguˇ si´ c, J. Gray, and G. K.-L. Chan, Fast and converged classical simulations of evidence for the utility of quantum computing before fault tolerance, Science Advances10, eadk4321 (2024), https://www.science.org/doi/pdf/10.1126/sciadv.adk4321

    2024 doi

  27. [35]

    M. S. Rudolph, T. Jones, Y. Teng, A. Angrisani, and Z. Holmes, Pauli propagation: A computational framework for simulating quantum systems (2025), arXiv:2505.21606 [quant-ph]

    2025 arXiv

  28. [36]

    Y. Shao, F. Wei, S. Cheng, and Z. Liu, Simulating noisy variational quantum algorithms: A polynomial approach, Physical Review Letters133, 120603 (2024)

    2024

  29. [37]

    Fontana, M

    E. Fontana, M. S. Rudolph, R. Duncan, I. Rungger, and C. C ˆ ırstoiu, Classical simulations of noisy variational quantum circuits, npj Quantum Information11, 84 (2025)

    2025

  30. [38]

    M. S. Rudolph, E. Fontana, Z. Holmes, and L. Cincio, Classical surrogate simulation of quantum systems with lowesa (2023), arXiv:2308.09109 [quant-ph]

    2023 arXiv

  31. [39]

    Beguˇ si´ c and G

    T. Beguˇ si´ c and G. K.-L. Chan, Real-time operator evolution in two and three dimensions via sparse pauli dynamics, PRX Quantum6, 020302 (2025)

    2025

  32. [40]

    Miller, Z

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

    2025

  33. [41]

    Bettaque and B

    V. Bettaque and B. Swingle, Structure of the majorana clifford group (2025)

    2025

  34. [42]

    B. C. Hall,Lie Groups, Lie Algebras, and Representations(Graduate Texts in Mathematics, 2015)

    2015

  35. [43]

    Z. Dai, Y. Wu, T. Wang, and M. P. Zaletel, Fermionic isometric tensor network states in two dimensions, Phys. Rev. Lett. 134, 026502 (2025)

    2025

  36. [44]

    M. P. Zaletel, R. S. Mong, C. Karrasch, J. E. Moore, and F. Pollmann, Time-evolving a matrix product state with long-ranged interactions, Physical Review B91, 165112 (2015)

    2015

  37. [45]

    Haegeman, C

    J. Haegeman, C. Lubich, I. Oseledets, B. Vandereycken, and F. Verstraete, Unifying time evolution and optimization with matrix product states, Physical Review B94, 165116 (2016)

    2016

  38. [46]

    Surace and L

    J. Surace and L. Tagliacozzo, Fermionic Gaussian states: an introduction to numerical approaches, SciPost Phys. Lect. Notes , 54 (2022)

    2022

  39. [47]

    Krumnow, L

    C. Krumnow, L. Veis, ¨O. Legeza, and J. Eisert, Fermionic orbital optimization in tensor network states, Physical review letters117, 210402 (2016)

    2016

  40. [48]

    Krumnow, J

    C. Krumnow, J. Eisert, and ¨O. Legeza, Towards overcoming the entanglement barrier when simulating long-time evolution, arXiv preprint arXiv:1904.11999 (2019)

    1904 arXiv

  41. [49]

    E. H. Lieb and D. W. Robinson, The finite group velocity of quantum spin systems, Communications in mathematical physics28, 251 (1972)

    1972

  42. [50]

    (Anthony) Chen, A

    C.-F. (Anthony) Chen, A. Lucas, and C. Yin, Speed limits and locality in many-body quantum dynamics, Reports on Progress in Physics86, 116001 (2023)

    2023

  43. [51]

    Schollw¨ ock, The density-matrix renormalization group, Reviews of Modern Physics77, 259 (2005)

    U. Schollw¨ ock, The density-matrix renormalization group, Reviews of Modern Physics77, 259 (2005)

    2005

  44. [52]

    Fishman, S

    M. Fishman, S. R. White, and E. M. Stoudenmire, The ITensor Software Library for Tensor Network Calculations, SciPost Phys. Codebases , 4 (2022)

    2022

  45. [53]

    Fishman, S

    M. Fishman, S. R. White, and E. M. Stoudenmire, Codebase release 0.3 for ITensor, SciPost Phys. Codebases , 4 (2022)

    2022

  46. [54]

    Vicentini, D

    F. Vicentini, D. Hofmann, A. Szab´ o, D. Wu, C. Roth, C. Giuliani, G. Pescia, J. Nys, V. Vargas-Calder´ on, N. Astrakhantsev, et al., Netket 3: Machine learning toolbox for many-body quantum systems, SciPost Physics Codebases , 007 (2022)

    2022

  47. [55]

    D’Anna and M

    M. D’Anna and M. S. Rudolph, MajoranaPropagation.jl (2025)

    2025

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.