pith. sign in

arxiv: 2509.01637 · v2 · submitted 2025-09-01 · 🪐 quant-ph · cond-mat.quant-gas· cond-mat.str-el

Phase-Sensitive Measurements on a Fermi-Hubbard Quantum Processor

Pith reviewed 2026-05-18 19:27 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.quant-gascond-mat.str-el
keywords Fermi-Hubbard modelLoschmidt echoquantum simulationfermionic atomsoptical latticespectral propertiesimaginary time evolutionquantum quench
0
0 comments X

The pith

Hardware-efficient protocol measures complex Loschmidt echoes to reveal spectral properties of the Fermi-Hubbard model.

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

The paper develops a hardware-efficient way to measure complex Loschmidt echoes on a fermionic quantum processor implementing the Fermi-Hubbard model. It combines global quench dynamics with short imaginary time evolution realized through tailored pulse sequences from plaquette product states. This works for both half-filling and doped cases. The approach yields these echoes efficiently for large systems across a wide spectral range. Such measurements would enable access to spectral functions like the local density of states and open paths to finite-temperature studies in existing simulators.

Core claim

Numerical results show that complex Loschmidt echoes can be efficiently obtained for large many-body states over a broad spectral range using the proposed protocol. This allows one to measure spectral properties of the Fermi-Hubbard model, such as the local density of states.

What carries the argument

Architecture-tailored pulse sequences that realize short imaginary-time evolution starting from a product state of plaquettes, combined with global quench dynamics.

If this is right

  • Complex Loschmidt echoes become accessible for large many-body fermionic states.
  • Local density of states can be measured for the Fermi-Hubbard model at half-filling and finite doping.
  • The protocol supports investigation of finite-temperature properties in fermionic quantum simulators.

Where Pith is reading between the lines

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

  • This could connect to other quantum simulation platforms by adapting the pulse sequences for different architectures.
  • Extending the method might allow probing out-of-equilibrium dynamics in doped Hubbard systems.
  • Validation on small systems could guide improvements for larger-scale implementations.

Load-bearing premise

The architecture-tailored pulse sequences accurately realize the short imaginary-time evolution without introducing uncontrolled errors that invalidate the extracted complex Loschmidt echoes.

What would settle it

Running the protocol on a small system where exact Loschmidt echoes are known from classical computation and observing large discrepancies would falsify the accuracy claim.

Figures

Figures reproduced from arXiv: 2509.01637 by Alberto R. Cavallar, Benjamin F. Schiffer, J. Ignacio Cirac, Luis Escalera-Moreno, Philipp M. Preiss, Timon Hilker, Titus Franz.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (b). There are four main contributions with an overlap squared of around 8%, which partially overlap in the figure, corresponding to all possible combinations of N´eel states on both plaquettes. The next contributions with around 2% population each correspond to one of sixteen combinations of a N´eel state on one plaquette and a striped (e.g., [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Ground state energy density of the Fermi–Hubbard [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Shot noise analysis of the broadened local density of states (LDOS) for a [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Instantaneous Fock state population throughout the [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Adiabatic preparation of doped plaquettes. [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
read the original abstract

Fermionic quantum processors are a promising platform for quantum simulation of correlated fermionic matter. In this work, we study a hardware-efficient protocol for measuring complex expectation values of the time-evolution operator, commonly referred to as Loschmidt echoes, with fermions in an optical superlattice. We analyze the algorithm for the Fermi-Hubbard model at half-filling as well as at finite doping. The method relies on global quench dynamics and short imaginary time evolution, the latter being realized by architecture-tailored pulse sequences starting from a product state of plaquettes. Our numerical results show that complex Loschmidt echoes can be efficiently obtained for large many-body states over a broad spectral range. This allows one to measure spectral properties of the Fermi-Hubbard model, such as the local density of states, and paves the way for the study of finite-temperature properties in current fermionic quantum simulators.

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.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a hardware-efficient protocol for measuring complex Loschmidt echoes in the Fermi-Hubbard model on a fermionic quantum processor realized in an optical superlattice. The approach combines global quench dynamics with short imaginary-time evolution, the latter implemented via architecture-tailored pulse sequences that begin from a product state of plaquettes. Numerical simulations are reported to demonstrate that complex Loschmidt echoes can be obtained efficiently for large many-body states over a broad spectral range, enabling extraction of spectral properties such as the local density of states and opening routes to finite-temperature studies.

Significance. If the protocol holds, the work would be significant for advancing phase-sensitive measurements in current fermionic quantum simulators, where accessing complex time-evolution operators has been challenging. The numerical demonstration of efficiency on large systems is a clear strength that supports potential scalability. Credit is due for the architecture-specific tailoring of the pulse sequences and the focus on practical implementation in existing hardware.

major comments (2)
  1. [Method description (pulse sequences)] The central claim depends on the accuracy of the architecture-tailored pulse sequences in realizing the required short imaginary-time evolution e^{-τH} starting from plaquette product states. No explicit fidelity benchmarks, error analysis, or comparison to the ideal propagator are provided for the interacting Fermi-Hubbard Hamiltonian at half-filling or finite doping, leaving open the possibility of uncontrolled Trotter-like or lattice-specific errors that would directly affect the extracted complex echoes.
  2. [Numerical results] Numerical results section: efficiency is asserted across a broad spectral range for large many-body states, yet the manuscript provides no quantitative error bars, explicit system sizes, or direct comparisons to exact diagonalization benchmarks. This makes it impossible to verify whether post-hoc parameter choices influence the central claim that complex Loschmidt echoes can be reliably obtained.
minor comments (2)
  1. [Abstract] Clarify in the abstract and introduction the precise range of interaction strengths U/t and doping levels for which the numerical demonstrations were performed.
  2. [Throughout] Ensure consistent notation for the imaginary-time parameter τ and the Loschmidt echo definition across equations and figures.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for providing constructive feedback. We address each of the major comments below and have made revisions to incorporate additional details as suggested.

read point-by-point responses
  1. Referee: The central claim depends on the accuracy of the architecture-tailored pulse sequences in realizing the required short imaginary-time evolution e^{-τH} starting from plaquette product states. No explicit fidelity benchmarks, error analysis, or comparison to the ideal propagator are provided for the interacting Fermi-Hubbard Hamiltonian at half-filling or finite doping, leaving open the possibility of uncontrolled Trotter-like or lattice-specific errors that would directly affect the extracted complex echoes.

    Authors: We appreciate the referee's concern regarding the validation of the pulse sequences. The sequences are specifically designed for the optical superlattice architecture to implement the short imaginary-time evolution starting from plaquette product states. To address this, we have added explicit fidelity benchmarks and error analysis in the revised manuscript for the Fermi-Hubbard model at half-filling and finite doping. These benchmarks compare the realized evolution to the ideal propagator and demonstrate that errors are controlled for the short evolution times used, supporting the accuracy of the complex Loschmidt echoes. revision: yes

  2. Referee: Numerical results section: efficiency is asserted across a broad spectral range for large many-body states, yet the manuscript provides no quantitative error bars, explicit system sizes, or direct comparisons to exact diagonalization benchmarks. This makes it impossible to verify whether post-hoc parameter choices influence the central claim that complex Loschmidt echoes can be reliably obtained.

    Authors: We thank the referee for this comment on the numerical results. Our simulations demonstrate the efficiency of the protocol for large many-body states across a broad spectral range. In response, we have revised the numerical results section to include quantitative error bars, explicit system sizes, and comparisons to exact diagonalization for small systems that benchmark the method. The parameter choices are determined by the requirements of the hardware-efficient protocol and are not post-hoc, as explained in the methods section. These revisions allow for better verification of the central claims. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation chain is self-contained and independent of fitted inputs or self-citation load-bearing steps.

full rationale

The paper describes a hardware-efficient protocol relying on global quench dynamics combined with short imaginary-time evolution implemented via architecture-tailored pulse sequences starting from plaquette product states. Numerical results are presented to show that complex Loschmidt echoes can be obtained efficiently over a broad spectral range for the Fermi-Hubbard model, enabling extraction of quantities such as the local density of states. No equation or step reduces by construction to a parameter fitted from the target data itself, nor does any central claim rest on a self-citation chain that is unverified or equivalent to the input. The protocol is framed as a direct measurement method whose validity is assessed through independent numerical simulation rather than tautological re-derivation of the same observables.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are stated. The protocol implicitly assumes that the Fermi-Hubbard Hamiltonian can be realized with sufficient fidelity in the optical superlattice and that short imaginary-time evolution can be approximated by the described pulses.

axioms (2)
  • domain assumption The optical superlattice can be tuned to realize the Fermi-Hubbard model at half-filling and finite doping with controllable parameters.
    Stated as the target system in the abstract.
  • domain assumption Short imaginary-time evolution can be implemented via architecture-tailored pulse sequences starting from plaquette product states.
    Central to the measurement protocol described in the abstract.

pith-pipeline@v0.9.0 · 5710 in / 1518 out tokens · 32377 ms · 2026-05-18T19:27:18.093566+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

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

  1. Programming long-range interactions in analog quantum simulators

    quant-ph 2026-04 unverdicted novelty 6.0

    A hybrid optimization strategy using classical pre-compilation, iterative extrapolation, and noise-aware quantum refinement achieves orders-of-magnitude gains in fidelity for state preparation in analog simulators wit...

  2. Programmable Fermionic Quantum Processors with Globally Controlled Lattices

    quant-ph 2026-04 unverdicted novelty 6.0

    Constructive protocols prove that arbitrary fermionic quantum operations are achievable with global controls in optical lattice Fermi-Hubbard systems.

Reference graph

Works this paper leans on

66 extracted references · 66 canonical work pages · cited by 2 Pith papers · 3 internal anchors

  1. [1]

    For the doped case, remove a doublon on every fourth plaquette

    Prepare the initial plaquette ground states: (a) Initialize the Fock state consisting of doublons |↑ ↓⟩on every second site. For the doped case, remove a doublon on every fourth plaquette. (b) Adiabatically prepare the plaquette ground state

  2. [2]

    2(f), each described by the pulse sequence in Fig

    Implement short ITE by applying the unitariesQ4 m=1 Vm, from Fig. 2(f), each described by the pulse sequence in Fig. 3(a). Employ corresponding pulses for the doped plaquettes

  3. [3]

    Time-evolve the full system for a duration τ in the analog mode

  4. [4]

    Measure in the Fock basis by taking a snapshot of the state and count the occurrences of the initial Fock state

    For the measuring step: Reverse the adiabatic ramps (see Appendix A 3) and uncompute the state preparation. Measure in the Fock basis by taking a snapshot of the state and count the occurrences of the initial Fock state

  5. [5]

    The maximum time corresponds to a minimally possible filter width (cf

    Repeat all previous steps for each time step τ. The maximum time corresponds to a minimally possible filter width (cf. Fig. 4). Step 4 is only included when estimating the amplitudes r(τ ± iϑ), but not for the r(τ)

  6. [6]

    Quantum Systems

    Finally, perform classical post-processing of the collected data: (a) Estimate the survival probabilities r(τ)2 and r(τ ± iϑ)2 from the sampled data. (b) Include normalization factors [cf. Eq. (7)], which can be computed numerically. (c) Estimate the phase ϕ(τ) of each Loschmidt echo from the amplitudes r(τ ± iϑ) following Eq. (3) and the ensuing numerica...

  7. [7]

    Adiabatic preparation of plaquette states We outline the preparation of the plaquette ground state, which is the first step of the experimental protocol. Starting from two isolated double wells, each double well in the |↑ ↓, 0⟩ Fock state, the ground state can be prepared by a two-step, piecewise-linear adiabatic path — in the spirit of earlier proposals ...

  8. [8]

    We perform a variational ground state search using DMRG [47] on matrix-product states in a snake-like geometry for different values of the maximum bond dimension χ

    Numerical plaquette ground state energies We include the data used for the estimation of the plaquette ground state energy density for different plaquette sizes. We perform a variational ground state search using DMRG [47] on matrix-product states in a snake-like geometry for different values of the maximum bond dimension χ. The simulations are implemente...

  9. [9]

    Projective measurement for adiabatically prepared states We discuss how to undo the adiabatic preparation, which is required for estimating the amplitude of the Loschmidt echo through projective Fock-basis measurements (cf. Sec. IV). We denote the (perfectly) adiabatic state preparation as the unitary U0→1, mapping the ground state |ψ0⟩ of an initial Hami...

  10. [10]

    IV we describe how the unitaryV implementing the (normalized) ITE exp( −HΛϑ), with HΛ =P m Hm can be decomposed into local unitaries Vm as defined in Eq

    Details for tiling the ITE unitaries In Sec. IV we describe how the unitaryV implementing the (normalized) ITE exp( −HΛϑ), with HΛ =P m Hm can be decomposed into local unitaries Vm as defined in Eq. (7). We now show why these unitaries can be defined acting on the same initial state |Ψ⟩ instead of acting each on the evolved state |Ψm−1⟩ := Vm−1 · · · V1|Ψ...

  11. [11]

    Details on the numerical pulse optimization In Sec. IV we describe how a pulse sequence modulating the hopping strengths of the Fermi–Hubbard 12 Hamiltonian on the subsystem of plaquettes Pµ and Pν can realize the (normalized) ITE unitary Vµν. Here, we describe how the unitary is obtained with an optimal control approach in the spirit of the GRAPE protoco...

  12. [12]

    F ock state contributions to doped ITE An ITE pulse sequence for 12.5%-doped plaquette states was presented in Sec. VII [cf. Fig. 5(e)]. We include here the Fock state populations throughout the evolution in Fig. 8. Unlike the half-filled case in Fig. 3(b), there do not seem to be clear dominant contributions from specific Fock states for the considered d...

  13. [13]

    Alternative for doped plaquette states We briefly discuss an alternative way to construct an initial plaquette state for 12.5% doping. Instead of combining three half-filled plaquettes A with a doped plaquette B that has two spins removed, we can consider two A plaquettes and a doped plaquette C missing a spin-down particle together with a plaquette D whe...

  14. [14]

    Sachdev, Quantum Phase Transitions , 2nd ed

    S. Sachdev, Quantum Phase Transitions , 2nd ed. (Cambridge University Press, 2011)

  15. [15]

    P. A. Lee, N. Nagaosa, and X.-G. Wen, Doping a mott in- sulator: Physics of high-temperature superconductivity, Rev. Mod. Phys. 78, 17 (2006)

  16. [16]

    Jiang and T

    H.-C. Jiang and T. P. Devereaux, Superconductivity in the doped hubbard model and its interplay with next- nearest hopping t′, Science 365, 1424–1428 (2019)

  17. [17]

    M. Qin, T. Sch¨ afer, S. Andergassen, P. Corboz, and E. Gull, The hubbard model: A computational perspective, Annual Review of Condensed Matter Physics 13, 275–302 (2022)

  18. [18]

    Stanisic, J

    S. Stanisic, J. L. Bosse, F. M. Gambetta, R. A. Santos, W. Mruczkiewicz, T. E. O’Brien, E. Ostby, and A. Montanaro, Observing ground-state properties of the fermi-hubbard model using a scalable algorithm on a quantum computer, Nature Communications 13 (2022)

  19. [19]

    Bourgund, T

    D. Bourgund, T. Chalopin, P. Bojovi´ c, H. Schl¨ omer, S. Wang, T. Franz, S. Hirthe, A. Bohrdt, F. Grusdt, I. Bloch, and T. A. Hilker, Formation of individual stripes in a mixed-dimensional cold-atom fermi–hubbard system, Nature 637, 57–62 (2025)

  20. [20]

    M. Xu, L. H. Kendrick, A. Kale, Y. Gang, C. Feng, S. Zhang, A. W. Young, M. Lebrat, and M. Greiner, A neutral-atom hubbard quantum simulator in the cryogenic regime, Nature 642, 909–915 (2025)

  21. [21]

    W. S. Bakr, Z. Ba, and M. L. Prichard, Mi- croscopy of ultracold fermions in optical lattices (2025), arXiv:2507.04042 [cond-mat.quant-gas]

  22. [22]

    L. W. Cheuk, M. A. Nichols, K. R. Lawrence, M. Okan, H. Zhang, E. Khatami, N. Trivedi, T. Paiva, M. Rigol, and M. W. Zwierlein, Observation of spatial charge and spin correlations in the 2d fermi-hubbard model, Science 353, 1260–1264 (2016)

  23. [23]

    Koepsell, J

    J. Koepsell, J. Vijayan, P. Sompet, F. Grusdt, T. A. Hilker, E. Demler, G. Salomon, I. Bloch, and C. Gross, Imaging magnetic polarons in the doped fermi–hubbard model, Nature 572, 358–362 (2019)

  24. [24]

    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, Nature 629, 317–322 (2024)

  25. [25]

    Holten, L

    M. Holten, L. Bayha, K. Subramanian, S. Brandstetter, C. Heintze, P. Lunt, P. M. Preiss, and S. Jochim, Observation of cooper pairs in a mesoscopic two- dimensional fermi gas, Nature 606, 287–291 (2022)

  26. [26]

    Havl´ ıˇ cek, M

    V. Havl´ ıˇ cek, M. Troyer, and J. D. Whitfield, Operator locality in the quantum simulation of fermionic models, Physical Review A 95 (2017). 14

  27. [27]

    A. M. Childs, Y. Su, M. C. Tran, N. Wiebe, and S. Zhu, Theory of Trotter Error with Commutator Scaling, Physical Review X 11, 011020 (2021)

  28. [28]

    A. Y. Kitaev, Quantum measurements and the abelian stabilizer problem (1995), arXiv:quant-ph/9511026 [quant-ph]

  29. [29]

    M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information: 10th Anniversary Edition (Cambridge University Press, 2010)

  30. [30]

    Lin and Y

    L. Lin and Y. Tong, Heisenberg-limited ground-state energy estimation for early fault-tolerant quantum computers, PRX Quantum 3 (2022)

  31. [31]

    N. S. Blunt, L. Caune, R. Izs´ ak, E. T. Campbell, and N. Holzmann, Statistical Phase Estimation and Error Mitigation on a Superconducting Quantum Processor, PRX Quantum 4, 040341 (2023)

  32. [32]

    S. Lu, M. C. Ba˜ nuls, and J. I. Cirac, Algorithms for Quantum Simulation at Finite Energies, PRX Quantum 2, 020321 (2021)

  33. [33]

    Cleve, A

    R. Cleve, A. Ekert, C. Macchiavello, and M. Mosca, Quantum algorithms revisited, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 454, 339 (1998)

  34. [34]

    T. E. O’Brien, S. Polla, N. C. Rubin, W. J. Huggins, S. McArdle, S. Boixo, J. R. McClean, and R. Babbush, Error mitigation via verified phase estimation, PRX Quantum 2 (2021)

  35. [35]

    H´ emery, K

    K. H´ emery, K. Ghanem, E. Crane, S. L. Campbell, J. M. Dreiling, C. Figgatt, C. Foltz, J. P. Gaebler, J. Johansen, M. Mills, S. A. Moses, J. M. Pino, A. Ransford, M. Rowe, P. Siegfried, et al. , Measuring the Loschmidt Amplitude for Finite-Energy Properties of the Fermi- Hubbard Model on an Ion-Trap Quantum Computer, PRX Quantum 5 (2024)

  36. [36]

    Y. Yang, A. Christianen, M. C. Ba˜ nuls, D. S. Wild, and J. I. Cirac, Phase-Sensitive Quantum Measurement without Controlled Operations, Phys. Rev. Lett. 132, 220601 (2024)

  37. [37]

    Clinton, T

    L. Clinton, T. S. Cubitt, R. Garcia-Patron, A. Montanaro, S. Stanisic, and M. Stroeks, Quantum phase estimation without controlled unitaries (2024), arXiv:2410.21517 [quant-ph]

  38. [38]

    B. F. Schiffer, D. S. Wild, N. Maskara, M. D. Lukin, and J. I. Cirac, Hardware-efficient quantum phase estimation via local control (2025), arXiv:2506.18765 [quant-ph]

  39. [39]

    Bluvstein, H

    D. Bluvstein, H. Levine, G. Semeghini, T. T. Wang, S. Ebadi, M. Kalinowski, A. Keesling, N. Maskara, H. Pichler, M. Greiner, V. Vuleti´ c, and M. D. Lukin, A quantum processor based on coherent transport of entangled atom arrays, Nature 604, 451–456 (2022)

  40. [40]

    D. K. Mark, H.-Y. Hu, J. Kwan, C. Kokail, S. Choi, and S. F. Yelin, Efficiently measuring d-wave pairing and beyond in quantum gas microscopes (2024), arXiv:2412.13186 [cond-mat.quant-gas]

  41. [41]

    & Gonz´ alez-Tudela, A

    C. Tabares, C. Kokail, P. Zoller, D. Gonz´ alez- Cuadra, and A. Gonz´ alez-Tudela, Programming optical- lattice fermi-hubbard quantum simulators (2025), arXiv:2502.05067 [quant-ph]

  42. [42]

    Probing quantum many-body dynamics using subsystem loschmidt echos,

    S. Karch, S. Bandyopadhyay, Z.-H. Sun, A. Impertro, S. Huh, I. P. Rodr´ ıguez, J. F. Wienand, W. Ketterle, M. Heyl, A. Polkovnikov, I. Bloch, and M. Aidelsburger, Probing quantum many-body dynamics using subsys- tem loschmidt echos (2025), arXiv:2501.16995 [cond- mat.quant-gas]

  43. [43]

    Gkritsis, D

    F. Gkritsis, D. Dux, J. Zhang, N. Jain, C. Gogolin, and P. M. Preiss, Simulating chemistry with fermionic optical superlattices, PRX Quantum 6 (2025)

  44. [44]

    High-fidelity collisional quantum gates with fermionic atoms

    P. Bojovi´ c, T. Hilker, S. Wang, J. Obermeyer, M. Barendregt, D. Tell, T. Chalopin, P. M. Preiss, I. Bloch, and T. Franz, High-fidelity collisional quantum gates with fermionic atoms (2025), arXiv:2506.14711 [cond-mat.quant-gas]

  45. [45]

    Trebst, U

    S. Trebst, U. Schollw¨ ock, M. Troyer, and P. Zoller, d- wave resonating valence bond states of fermionic atoms in optical lattices, Physical Review Letters 96 (2006)

  46. [46]

    Nascimb` ene, Y.-A

    S. Nascimb` ene, Y.-A. Chen, M. Atala, M. Aidelsburger, S. Trotzky, B. Paredes, and I. Bloch, Experimental Realization of Plaquette Resonating Valence-Bond States with Ultracold Atoms in Optical Superlattices, Physical Review Letters 108, 205301 (2012)

  47. [47]

    Jordan and E

    P. Jordan and E. Wigner, ¨Uber das Paulische ¨Aquivalenzverbot, Zeitschrift fur Physik 47, 631 (1928)

  48. [48]

    S. B. Bravyi and A. Y. Kitaev, Fermionic quantum computation, Annals of Physics 298, 210–226 (2002)

  49. [49]

    Verstraete and J

    F. Verstraete and J. I. Cirac, Mapping local hamiltonians of fermions to local hamiltonians of spins, Journal of Statistical Mechanics: Theory and Experiment 2005, P09012–P09012 (2005)

  50. [50]

    J. Yu, Y. Liu, S. Sugiura, T. V. Voorhis, and S. Zeytino˘ glu, Clifford circuit based heuristic optimization of fermion-to-qubit mappings (2025), arXiv:2502.11933 [quant-ph]

  51. [51]

    Corboz, Improved energy extrapolation with infinite projected entangled-pair states applied to the two- dimensional hubbard model, Physical Review B 93 (2016)

    P. Corboz, Improved energy extrapolation with infinite projected entangled-pair states applied to the two- dimensional hubbard model, Physical Review B 93 (2016)

  52. [52]

    Quantum Computation by Adiabatic Evolution

    E. Farhi, J. Goldstone, S. Gutmann, and M. Sipser, Quantum computation by adiabatic evolution (2000), arXiv:quant-ph/0001106 [quant-ph]

  53. [53]

    W. S. Bakr, J. I. Gillen, A. Peng, S. F¨ olling, and M. Greiner, A quantum gas microscope for detecting single atoms in a hubbard-regime optical lattice, Nature 462, 74–77 (2009)

  54. [54]

    M. F. Parsons, F. Huber, A. Mazurenko, C. S. Chiu, W. Setiawan, K. Wooley-Brown, S. Blatt, and M. Greiner, Site-resolved imaging of fermionic 6li in an optical lattice, Physical Review Letters 114 (2015)

  55. [55]

    L. W. Cheuk, M. A. Nichols, M. Okan, T. Gersdorf, V. V. Ramasesh, W. S. Bakr, T. Lompe, and M. W. Zwierlein, Quantum-gas microscope for fermionic atoms, Physical Review Letters 114 (2015)

  56. [56]

    Khaneja, T

    N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbr¨ uggen, and S. J. Glaser, Optimal control of coupled spin dynamics: design of NMR pulse sequences by gradient ascent algorithms, Journal of Magnetic Resonance 172, 296 (2005)

  57. [57]

    A. M. Childs and N. Wiebe, Hamiltonian Simula- tion Using Linear Combinations of Unitary Oper- ations, Quantum Information and Computation 12, 10.26421/QIC12.11-12 (2012)

  58. [58]

    There exist other choices to break the degeneracy, such as µk = δk,0t + δk,1t, but they might be more difficult to implement experimentally

  59. [59]

    Trivedi, A

    R. Trivedi, A. F. Rubio, and J. I. Cirac, Quantum advantage and stability to errors in analogue quantum simulators, Nature Communications 15, 6507 (2024)

  60. [60]

    Schollw¨ ock, The density-matrix renormalization group 15 in the age of matrix product states, Annals of physics 326, 96 (2011), publisher: Elsevier

    U. Schollw¨ ock, The density-matrix renormalization group 15 in the age of matrix product states, Annals of physics 326, 96 (2011), publisher: Elsevier

  61. [61]

    Fishman, S

    M. Fishman, S. White, and E. Stoudenmire, The itensor software library for tensor network calculations, SciPost Physics Codebases (2022)

  62. [62]

    B. F. Schiffer, D. v. Vreumingen, J. Tura, and S. Polla, Virtual mitigation of coherent non-adiabatic transitions by echo verification, Quantum 8, 1346 (2024)

  63. [63]

    R. H. Byrd, P. Lu, J. Nocedal, and C. Zhu, A Limited Memory Algorithm for Bound Constrained Optimization, SIAM Journal on Scientific Computing16, 1190 (1995)

  64. [64]

    C. Zhu, R. H. Byrd, P. Lu, and J. Nocedal, Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization, ACM Trans. Math. Softw. 23, 550–560 (1997)

  65. [65]

    Weinberg and M

    P. Weinberg and M. Bukov, QuSpin: a python package for dynamics and exact diagonalisation of quantum many body systems part I: spin chains, SciPost Physics 2 (2017)

  66. [66]

    Weinberg and M

    P. Weinberg and M. Bukov, QuSpin: a python package for dynamics and exact diagonalisation of quantum many body systems. Part II: bosons, fermions and higher spins, SciPost Physics 7 (2019)