Recognition: unknown
Gauge-covariant projected entangled paired states for interacting systems in a magnetic field
Pith reviewed 2026-05-08 04:10 UTC · model grok-4.3
The pith
Projected entangled-pair states with virtual flux tensors simulate interacting systems in uniform magnetic fields while keeping all physical observables translation invariant by construction.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By introducing a pattern of virtual flux tensors into the PEPS, the wavefunction can be made such that all physical observables are translation invariant by construction, possibly within an enlarged unit cell if symmetry breaking occurs. The usual contraction and optimization methods for translation-invariant PEPS apply directly, with the flux per plaquette acting as a continuous parameter. This approach works independently of the gauge choice for the vector potential and avoids considering extended magnetic unit cells.
What carries the argument
Virtual flux tensors arranged in a pattern within the PEPS, which encode the magnetic field effects virtually to enforce physical translation invariance.
If this is right
- The method enables direct use of standard PEPS algorithms for magnetic systems.
- Physical quantities can be computed without dependence on the vector potential gauge.
- Enlarged unit cells are only needed if the target state itself breaks translation symmetry.
- Interacting many-body systems in magnetic fields become simulable with continuous flux parameter.
Where Pith is reading between the lines
- This framework might extend to time-dependent simulations or excited states in magnetic fields.
- It could facilitate comparisons between different lattice models of quantum Hall effects.
- Applications to systems with inhomogeneous fields or other gauge structures may follow from similar virtual tensor patterns.
Load-bearing premise
A suitable pattern of virtual flux tensors exists for any target ground state, allowing the PEPS to produce translation-invariant observables while permitting unmodified use of standard contraction and optimization algorithms.
What would settle it
A calculation on a specific model such as the Hofstadter-Hubbard model where the ground state energy or correlation functions computed with this PEPS deviate from known exact or alternative numerical results when the virtual flux pattern is altered or removed.
Figures
read the original abstract
The Hamiltonian for a system of itinerant particles on a two-dimensional lattice in a uniform magnetic field reduces the translational symmetry to a magnetic translation group, because of the need to choose a particular gauge for the vector potential. Nonetheless, in many situations all physical observables of the ground state remain entirely translation invariant. In this work, we introduce a projected entangled-pair state (PEPS) wavefunction with a pattern of virtual flux tensors, for which all physical expectation values are translation invariant by construction, possibly within an enlarged unit cell reflecting any symmetry breaking in the target state. Moreover, we show that the usual contraction and optimization methods for translation-invariant PEPS can be used, with the magnetic flux per plaquette only entering as a continuous parameter in the tensor network contractions. Therefore, our approach provides a method for simulating an interacting many-body system in a uniform magnetic field independently of the gauge choice for the vector potential and bypassing the need to consider extended magnetic unit cells.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a projected entangled-pair state (PEPS) ansatz incorporating a fixed pattern of virtual flux tensors. This construction is designed so that all physical observables remain translation-invariant by construction (possibly on an enlarged cell reflecting any spontaneous symmetry breaking), while the magnetic flux per plaquette enters the tensor network only as a continuous scalar parameter. The authors claim that standard translation-invariant PEPS contraction and optimization algorithms can be applied without modification, yielding a gauge-independent representation that avoids the need to enlarge the unit cell to a magnetic supercell.
Significance. If the virtual-flux pattern can be shown to exist for arbitrary interacting ground states and to preserve unmodified contraction/optimization routines, the approach would constitute a meaningful technical advance for tensor-network studies of lattice models in uniform magnetic fields (e.g., Hofstadter-Hubbard or fractional quantum Hall systems). It would remove a long-standing practical obstacle—gauge dependence and the associated proliferation of magnetic unit cells—while retaining the computational advantages of translation-invariant PEPS.
major comments (3)
- [construction section / abstract] The central claim (abstract and the construction section) that a universal, state-independent pattern of virtual flux tensors exists for any target ground state such that (i) physical observables are translation-invariant by construction, (ii) the representation is gauge-independent, and (iii) standard TI-PEPS contraction and optimization routines remain unmodified is load-bearing but not demonstrated. No general algorithmic procedure or existence proof is supplied; the pattern is introduced as an ansatz whose suitability for a given interacting state is assumed rather than constructed or verified.
- [contraction section] § on contraction properties: the assertion that the flux per plaquette enters the network solely as a continuous parameter inside otherwise unchanged TI-PEPS routines is not accompanied by explicit tensor definitions, bond-dimension scaling, or numerical benchmarks on interacting models. Without these, it is impossible to confirm that the virtual fluxes do not alter the contraction graph or introduce additional computational cost.
- [discussion / numerical results] The weakest assumption—that a suitable virtual-flux pattern always exists without state-dependent modifications to the tensor network or its contraction graph—is not tested against known counter-examples (e.g., states with complex anyonic statistics or strong gauge-dependent correlations). A concrete counter-example or a proof of existence would be required to substantiate the bypass of extended magnetic unit cells.
minor comments (2)
- [tensor definitions] Notation for the virtual flux tensors is introduced without an explicit index diagram or bond-dimension table; adding one would improve readability.
- [introduction] The manuscript would benefit from a short comparison table contrasting the new ansatz with conventional magnetic-unit-cell PEPS in terms of bond dimension, gauge freedom, and contraction cost.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive comments. We address each major comment point by point below. Revisions have been made to clarify the ansatz nature of the construction, provide explicit tensor details, and expand the discussion of limitations.
read point-by-point responses
-
Referee: [construction section / abstract] The central claim (abstract and the construction section) that a universal, state-independent pattern of virtual flux tensors exists for any target ground state such that (i) physical observables are translation-invariant by construction, (ii) the representation is gauge-independent, and (iii) standard TI-PEPS contraction and optimization routines remain unmodified is load-bearing but not demonstrated. No general algorithmic procedure or existence proof is supplied; the pattern is introduced as an ansatz whose suitability for a given interacting state is assumed rather than constructed or verified.
Authors: We thank the referee for this observation. The virtual flux pattern is an ansatz, but one that is explicitly state-independent: the flux tensors are fixed solely by the magnetic flux per plaquette and the lattice geometry to enforce gauge covariance of the physical wavefunction. This ensures translation invariance of all physical observables by construction, without reference to the target state. We have revised the construction section to include a step-by-step algorithmic procedure for inserting the flux tensors and clarified that the pattern does not depend on the ground state. A general existence proof for arbitrary states lies outside the scope of this work, which focuses on practical gauge-independent simulations for interacting lattice models. revision: yes
-
Referee: [contraction section] § on contraction properties: the assertion that the flux per plaquette enters the network solely as a continuous parameter inside otherwise unchanged TI-PEPS routines is not accompanied by explicit tensor definitions, bond-dimension scaling, or numerical benchmarks on interacting models. Without these, it is impossible to confirm that the virtual fluxes do not alter the contraction graph or introduce additional computational cost.
Authors: We agree that these details strengthen the claim. The revised manuscript now includes explicit definitions of the virtual flux tensors and their placement within the PEPS network. The bond dimension is unchanged, as the fluxes are encoded via phase factors on existing virtual indices. The contraction graph and optimization routines remain identical to standard TI-PEPS, with the flux entering only as a continuous scalar parameter. A scaling analysis confirming no extra cost has been added. Numerical benchmarks on interacting models are deferred to future work, as the present manuscript centers on the formal construction. revision: yes
-
Referee: [discussion / numerical results] The weakest assumption—that a suitable virtual-flux pattern always exists without state-dependent modifications to the tensor network or its contraction graph—is not tested against known counter-examples (e.g., states with complex anyonic statistics or strong gauge-dependent correlations). A concrete counter-example or a proof of existence would be required to substantiate the bypass of extended magnetic unit cells.
Authors: We acknowledge this limitation. The pattern is constructed to be independent of the state and requires no modifications to the tensor network or contraction graph for any flux value. We have expanded the discussion section to address potential counter-examples, such as states with anyonic statistics, noting that our ansatz applies directly to models with translation-invariant observables. While we do not provide a universal proof or exhaustive counter-example search, the construction demonstrably bypasses magnetic supercells for the Hofstadter-type systems considered. Further exploration of exotic cases is left for future study. revision: partial
Circularity Check
No significant circularity; construction is self-contained by design
full rationale
The manuscript introduces a specific PEPS ansatz incorporating a fixed pattern of virtual flux tensors chosen so that physical observables are translation-invariant by construction (possibly on an enlarged cell) and the flux per plaquette appears only as a scalar parameter in otherwise standard TI-PEPS contractions. This is a definitional choice of wavefunction form rather than a derivation in which a claimed result is obtained by fitting to itself or by a self-referential loop. No load-bearing step reduces an output to an input by construction, no self-citation is invoked to establish uniqueness or existence, and the flux is treated as an external physical input. The central claim therefore stands as an independent proposal whose validity rests on whether the ansatz can represent the target states, not on internal circularity.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Standard PEPS contraction and optimization algorithms remain valid when the virtual flux pattern is inserted and the physical flux per plaquette is supplied as a continuous parameter.
invented entities (1)
-
virtual flux tensors
no independent evidence
Reference graph
Works this paper leans on
-
[1]
and the fermionic [23, 24] versions of this model, but the need for accommodating the magnetic unit cell puts a stringent limitation of the values for the flux. Physically, however, we expect that implementing these large unit cells is overkill: While the Peierls substi- tution introduces position-dependent phases in the hop- ping amplitudes, this does no...
-
[2]
A for the 4 explicit entries up to second-order perturbation theory)
when applied to the product state (see App. A for the 4 explicit entries up to second-order perturbation theory). As a result, the flux tensors in Eq. (16) will appear on the virtual level of the resulting PEPS, with the pattern matching the gauge choice in the Hamiltonian. When choosing the Landau gaugeA=−Bxe y as in Fig. 1, we find precisely the pattern...
-
[3]
J. M. Luttinger, The effect of a magnetic field on elec- trons in a periodic potential, Phys. Rev.84, 814 (1951)
1951
-
[4]
Kohn, Theory of Bloch electrons in a magnetic field: The effective Hamiltonian, Phys
W. Kohn, Theory of Bloch electrons in a magnetic field: The effective Hamiltonian, Phys. Rev.115, 1460 (1959)
1959
-
[5]
Peierls, Zur theorie des diamagnetismus von leitungse- lektronen, Zeitschrift f¨ ur Physik80, 763 (1933)
R. Peierls, Zur theorie des diamagnetismus von leitungse- lektronen, Zeitschrift f¨ ur Physik80, 763 (1933)
1933
-
[6]
P. G. Harper, Single Band Motion of Conduction Elec- trons in a Uniform Magnetic Field, Proceedings of the Physical Society. Section A68, 874 (1955)
1955
-
[7]
M. Ya. Azbel’, Energy Spectrum of a Conduction Elec- tron in a Magnetic Field, Journal of Experimental and Theoretical Physics19, 634 (1964)
1964
-
[8]
D. R. Hofstadter, Energy levels and wave functions of Bloch electrons in rational and irrational magnetic fields, Phys. Rev. B14, 2239 (1976)
1976
-
[9]
Brown, Bloch electrons in a uniform magnetic field, Phys
E. Brown, Bloch electrons in a uniform magnetic field, Phys. Rev.133, A1038 (1964)
1964
-
[10]
Zak, Magnetic translation group, Phys
J. Zak, Magnetic translation group, Phys. Rev.134, A1602 (1964)
1964
-
[11]
A. S. Sørensen, E. Demler, and M. D. Lukin, Fractional 6 quantum Hall states of atoms in optical lattices, Phys. Rev. Lett.94, 086803 (2005)
2005
-
[12]
M¨ oller and N
G. M¨ oller and N. R. Cooper, Composite fermion theory for bosonic quantum hall states on lattices, Phys. Rev. Lett.103, 105303 (2009)
2009
-
[13]
Sterdyniak, N
A. Sterdyniak, N. Regnault, and G. M¨ oller, Particle en- tanglement spectra for quantum hall states on lattices, Phys. Rev. B86, 165314 (2012)
2012
-
[14]
A. M. L¨ auchli, Z. Liu, E. J. Bergholtz, and R. Moessner, Hierarchy of fractional Chern insulators and competing compressible states, Phys. Rev. Lett.111, 126802 (2013)
2013
-
[15]
M¨ oller and N
G. M¨ oller and N. R. Cooper, Fractional Chern insulators in Harper-Hofstadter bands with higher Chern number, Phys. Rev. Lett.115, 126401 (2015)
2015
-
[16]
J. I. Cirac, D. P´ erez-Garc´ ıa, N. Schuch, and F. Ver- straete, Matrix product states and projected entangled pair states: Concepts, symmetries, theorems, Rev. Mod. Phys.93, 045003 (2021)
2021
-
[17]
Xiang,Density Matrix and Tensor Network Renor- malization(Cambridge University Press, 2023)
T. Xiang,Density Matrix and Tensor Network Renor- malization(Cambridge University Press, 2023)
2023
-
[18]
Gerster, M
M. Gerster, M. Rizzi, P. Silvi, M. Dalmonte, and S. Mon- tangero, Fractional quantum Hall effect in the interacting Hofstadter model via tensor networks, Phys. Rev. B96, 195123 (2017)
2017
-
[19]
D. L. Kovrizhin, Density matrix renormalization group for bosonic quantum hall effect, Phys. Rev. B81, 125130 (2010)
2010
-
[20]
W. Zhu, S. S. Gong, F. D. M. Haldane, and D. N. Sheng, Topological characterization of the non-abelian Moore- Read state using density-matrix renormalization group, Phys. Rev. B92, 165106 (2015)
2015
-
[21]
A. G. Grushin, J. Motruk, M. P. Zaletel, and F. Poll- mann, Characterization and stability of a fermionicν= 1/3 fractional Chern insulator, Phys. Rev. B91, 035136 (2015)
2015
-
[22]
Motruk, A
J. Motruk, A. G. Grushin, F. de Juan, and F. Pollmann, Interaction-driven phases in the half-filled honeycomb lattice: An infinite density matrix renormalization group study, Phys. Rev. B92, 085147 (2015)
2015
-
[23]
Renormalization algorithms for Quantum-Many Body Systems in two and higher dimensions
F. Verstraete and J. I. Cirac, Renormalization algorithms for quantum-many body systems in two and higher di- mensions, arXiv:cond-mat/0407066 (2004)
work page internal anchor Pith review arXiv 2004
-
[24]
E. L. Weerda and M. Rizzi, Fractional quantum Hall states with variational projected entangled-pair states: A study of the bosonic Harper-Hofstadter model, Phys. Rev. B109, L241117 (2024)
2024
- [25]
- [26]
- [27]
-
[28]
Vanderstraeten, M
L. Vanderstraeten, M. Mari¨ en, J. Haegeman, N. Schuch, J. Vidal, and F. Verstraete, Bridging perturbative expan- sions with tensor networks, Phys. Rev. Lett.119, 070401 (2017)
2017
-
[29]
Haegeman and F
J. Haegeman and F. Verstraete, Diagonalizing transfer matrices and matrix product operators: A medley of ex- act and computational methods, Annual Review of Con- densed Matter Physics8, 355 (2017)
2017
-
[30]
M. T. Fishman, L. Vanderstraeten, V. Zauner-Stauber, J. Haegeman, and F. Verstraete, Faster methods for con- tracting infinite two-dimensional tensor networks, Phys. Rev. B98, 235148 (2018)
2018
-
[31]
Vanderstraeten, L
L. Vanderstraeten, L. Burgelman, B. Ponsioen, M. Van Damme, B. Vanhecke, P. Corboz, J. Haegeman, and F. Verstraete, Variational methods for contracting projected entangled-pair states, Phys. Rev. B105, 195140 (2022)
2022
-
[32]
Vanderstraeten, J
L. Vanderstraeten, J. Haegeman, P. Corboz, and F. Ver- straete, Gradient methods for variational optimization of projected entangled-pair states, Phys. Rev. B94, 155123 (2016)
2016
-
[33]
Liao, J.-G
H.-J. Liao, J.-G. Liu, L. Wang, and T. Xiang, Differ- entiable programming tensor networks, Phys. Rev. X9, 031041 (2019)
2019
-
[34]
Capogrosso-Sansone, S
B. Capogrosso-Sansone, S. G. S¨ oyler, N. Prokof’ev, and B. Svistunov, Monte Carlo study of the two-dimensional Bose-Hubbard model, Phys. Rev. A77, 015602 (2008)
2008
-
[35]
R. O. Umucalilar and M. O. Oktel, Phase boundary of the boson mott insulator in a rotating optical lattice, Phys. Rev. A76, 055601 (2007)
2007
-
[36]
S. S. Natu, E. J. Mueller, and S. Das Sarma, Competing ground states of strongly correlated bosons in the Harper- Hofstadter-Mott model, Phys. Rev. A93, 063610 (2016)
2016
-
[37]
Song and S.-J
Y.-F. Song and S.-J. Yang, Quantum phases for bosons in a square lattice with artificial magnetic flux, Journal of Physics B: Atomic, Molecular and Optical Physics52, 118001 (2019)
2019
-
[38]
D. S. Goldbaum and E. J. Mueller, Vortex lattices of bosons in deep rotating optical lattices, Phys. Rev. A 77, 033629 (2008)
2008
-
[39]
Powell, R
S. Powell, R. Barnett, R. Sensarma, and S. Das Sarma, Bogoliubov theory of interacting bosons on a lattice in a synthetic magnetic field, Phys. Rev. A83, 013612 (2011)
2011
-
[40]
M¨ oller and N
G. M¨ oller and N. R. Cooper, Condensed ground states of frustrated Bose-Hubbard models, Phys. Rev. A82, 063625 (2010)
2010
-
[41]
Corboz, Variational optimization with infinite pro- jected entangled-pair states, Phys
P. Corboz, Variational optimization with infinite pro- jected entangled-pair states, Phys. Rev. B94, 035133 (2016)
2016
-
[42]
Zauner-Stauber, L
V. Zauner-Stauber, L. Vanderstraeten, M. T. Fishman, F. Verstraete, and J. Haegeman, Variational optimization algorithms for uniform matrix product states, Phys. Rev. B97, 045145 (2018)
2018
-
[43]
Jiang and Y
S. Jiang and Y. Ran, Symmetric tensor networks and practical simulation algorithms to sharply identify classes of quantum phases distinguishable by short-range physics, Phys. Rev. B92, 104414 (2015)
2015
-
[44]
Mambrini, R
M. Mambrini, R. Or´ us, and D. Poilblanc, Systematic con- struction of spin liquids on the square lattice from tensor networks with SU(2) symmetry, Phys. Rev. B94, 205124 (2016)
2016
-
[45]
Hackenbroich, A
A. Hackenbroich, A. Sterdyniak, and N. Schuch, Inter- play of SU(2), point group, and translational symmetry for projected entangled pair states: Application to a chi- ral spin liquid, Phys. Rev. B98, 085151 (2018)
2018
-
[46]
Hasik and P
J. Hasik and P. Corboz, Incommensurate order with translationally invariant projected entangled-pair states: Spiral states and quantum spin liquid on the anisotropic triangular lattice, Phys. Rev. Lett.133, 176502 (2024)
2024
-
[47]
Mortier, L
Q. Mortier, L. Devos, L. Burgelman, B. Vanhecke, N. Bultinck, F. Verstraete, J. Haegeman, and L. Van- derstraeten, Fermionic tensor network methods, SciPost Phys.18, 012 (2025). 7
2025
-
[48]
T. B. Wahl, H.-H. Tu, N. Schuch, and J. I. Cirac, Pro- jected entangled-pair states can describe chiral topologi- cal states, Phys. Rev. Lett.111, 236805 (2013)
2013
-
[49]
Dubail and N
J. Dubail and N. Read, Tensor network trial states for chiral topological phases in two dimensions and a no- go theorem in any dimension, Phys. Rev. B92, 205307 (2015)
2015
-
[50]
S. Yang, T. B. Wahl, H.-H. Tu, N. Schuch, and J. I. Cirac, Chiral projected entangled-pair state with topological or- der, Phys. Rev. Lett.114, 106803 (2015)
2015
-
[51]
Poilblanc, J
D. Poilblanc, J. I. Cirac, and N. Schuch, Chiral topo- logical spin liquids with projected entangled pair states, Phys. Rev. B91, 224431 (2015)
2015
-
[52]
J.-Y. Chen, L. Vanderstraeten, S. Capponi, and D. Poil- blanc, Non-abelian chiral spin liquid in a quantum anti- ferromagnet revealed by an iPEPS study, Phys. Rev. B 98, 184409 (2018)
2018
-
[53]
Hasik, M
J. Hasik, M. Van Damme, D. Poilblanc, and L. Vander- straeten, Simulating chiral spin liquids with projected entangled-pair states, Phys. Rev. Lett.129, 177201 (2022). Appendix A: Perturbative expansion A general procedure for finding the most natural struc- ture for a tensor network representation in a given gapped quantum phase is by perturbation theory...
2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.