Pith. Sign in

REVIEW

Fermion bag study of the two-dimensional Majorana-Hubbard model at finite hopping

Strongly Correlated Electrons · Submitted 22 Sep 2026 · 56 references

Pith's verdict: First finite-hopping QMC for the 2D Majorana-Hubbard model; the ED-validated algorithmic core and channel trends are real, but the gapless-semimetal claim rides on three-size power-law fits at fixed beta=L/4.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 24, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.26881 (version 1)
Pith IDQZOHEMGH, Pith's permanent record of this paper and its review
CitationEmilie Huffman, Mark Johnston. "Fermion bag study of the two-dimensional Majorana-Hubbard model at finite hopping." arXiv:2609.26881 (2026). https://pith.science/paper/QZOHEMGH
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Stable Pfaffian weights let quantum Monte Carlo reach finite hopping in the Majorana-Hubbard model, where data show a gapless Dirac semimetal with the charged dimerization doublet enhanced and the neutral channel suppressed as the…

The reading

The two-dimensional Majorana-Hubbard model places one Majorana fermion—a particle that is its own antiparticle—on each site of a $\pi$-flux square lattice with a plaquette interaction, making it one of the most minimal lattice models of interacting relativistic fermions. This paper claims to bring the model at nonzero hopping within reach of quantum Monte Carlo by expanding in continuous-time fermion bags and computing each configuration's weight as a stable Pfaffian that also fixes the sign; below $g\lesssim 0.7t$ the sign problem is fully absent for lattices up to 256 sites at the inverse temperatures studied. Finite hopping matters because the only earlier unbiased numerics sat at zero hopping, where the model reduces to a classical problem and the predicted semimetal-to-dimerization transition could not be tested. Within the sign-mild window the data are consistent with a gapless Dirac semimetal (a phase whose low-energy excitations are massless fermions with linear dispersion), and increasing $g$ strengthens the charged columnar-dimer doublet—the two-component order parameter that rotates under the emergent $U(1)$—while suppressing the neutral singlet channel, as the renormalization group predicts; the paper further reports that the emergent $U(1)$ symmetry is partially restored in the even bond-orientation component of the correlator while an odd, lattice-nematic anisotropy persists.

What carries the argument

The load-bearing machinery is the Hamiltonian fermion bag expansion with Pfaffian weights. Writing the Hamiltonian as $H=-\mu\sum_p D_p$ with $D_p=(1-ia\eta_{ij}\gamma_i\gamma_j)(1-ia\eta_{kl}\gamma_k\gamma_l)$ turns the partition function into a continuous-time sum over products of plaquette operators; the trace of each term is exactly $\operatorname{Tr}(1)\operatorname{Pf}(A)$, where $A$ is a $4k\times4k$ antisymmetric matrix built from the expansion order and the hopping factors. This Pfaffian gives the magnitude and the sign of each configuration stably, replacing the square-root determinant formula that is fast but leaves a sign ambiguity. The second piece is control of the sign problem by boundary conditions: the free hopping matrix has zero modes at $(\pm\pi/2,0)$, and with periodic boundary conditions these lie exactly on the lattice momenta whenever $4\mid L$, causing the sign to collapse, while antiperiodic boundary conditions shift the momenta by $\pi/L$ and keep the average sign flat. The physics comparison is carried by the three local operators that can gap the Dirac cone—the columnar dimer doublet $S_x,S_y$, which carries charge $2$ under the emergent $U(1)$, and the neutral singlet $M$—whose connected correlators at $r=L/2$ show the opposite trends as $g$ grows.

Extended reading notes

Core claim

The paper's central claim is that the finite-hopping Majorana-Hubbard model on the $\pi$-flux square lattice can be simulated with quantum Monte Carlo without a severe sign problem, and that the accessible data realize the weak-coupling side of the predicted phase diagram. Each continuous-time configuration's trace is expressed as $\operatorname{Tr}(1)\operatorname{Pf}(A)$, with an antisymmetric matrix $A$ of dimension $4k$ set by the number of inserted plaquettes, which is stable and fixes the sign; antiperiodic boundary conditions move the free-theory zero modes off the lattice momenta, and the average sign then remains one up to $g\simeq 0.6t$ and mild up to $g\simeq 0.9t$. The algorithm is validated on $L=4$ systems against exact diagonalization for the same correlators. On lattices $L=8,12,16$ at $\beta=L/4$, the charged-dimer and neutral-mass correlators at separation $r=L/2$ follow power laws in $L$ up to $g=0.9t$, which the paper reads as a gapless Dirac semimetal with no condensed mass in this window. The charged doublet is enhanced relative to the free theory (to $3.9\pm0.3$ times the free value at $L=16$, $g=0.7t$) while the neutral channel is suppressed (to $0.40$ of free), with the two exponents moving in opposite directions, matching the renormalization-group prediction for $g>0$. The paper also claims that the approach to the emergent $U(1)$ symmetry is only partial: the off-diagonal, orientation-even component of the $\Phi$ doublet correlator is driven toward zero by increasing $g$, while the orientation-odd anisotropy component acquires a negative shift that persists with distance, leaving a lattice-nematic $C_4$-odd signal.

Load-bearing premise

The conclusions assume that three lattice sizes (8, 12, and 16 sites across), each simulated for a time interval equal to one quarter of its linear size, already behave like the zero-temperature infinite lattice when a power law is fitted to their farthest-separation correlators, and no check of longer simulated times is reported.

Editorial extensions

If this is right

  • Because the sign is under control in the window $g\lesssim 0.9t$, larger lattices and reweighting into the sign-problematic region become feasible; the paper puts the predicted transition near $g\simeq 0.9t$ at the edge of what it can currently reach.
  • If the central claim holds, the sign-mild region realizes a gapless Dirac semimetal with no condensed mass, so the columnar dimerization transition must lie at or beyond the edge of this window.
  • The opposite trends of the charged doublet and neutral singlet identify the $U(1)$-charged columnar dimer as the leading instability for $g>0$, in line with the renormalization-group flow toward the supersymmetric side.
  • The partial restoration of $U(1)$—even component restoring, odd component not—means any lattice calculation approaching the critical point must contend with a persistent $C_4$ anisotropy, which may be dangerously irrelevant or a separate nematic tendency.

Reading between the lines

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

  • Beyond the paper's explicit conclusions, the persistent $C_4$-odd anisotropy suggests a concrete way to separate the two proposed interpretations: measure whether a nematic order parameter shows a finite-size crossing near $g\simeq 0.9t$ and, if so, whether it crosses at a different coupling from the dimerization signal.
  • The zero-mode mechanism the paper identifies may transfer to other interacting Majorana models: shifting accidental zero modes off the lattice momenta by boundary conditions could create sign-mild simulation windows wherever free Majorana cones sit exactly at allowed momenta.
  • If the gapless-semimetal interpretation survives larger systems, this model becomes a minimal lattice realization of 2+1-dimensional interacting Dirac physics with a single two-component Dirac fermion, and the measured channel asymmetry can be used to test higher-order continuum renormalization-group predictions.
Original abstract

The two-dimensional Majorana-Hubbard model on the $\pi$-flux lattice involves the most local interaction possible for a lattice with one relativistic Majorana fermion per site. It is predicted to host a Majorana semimetal, dimerized phases, and quantum-critical points with an emergent $U(1)$ symmetry. Nonperturbative lattice calculations for this model remain a major open challenge, thus far limited to zero hopping. We present quantum Monte Carlo simulations of the model at finite hopping using Hamiltonian fermion bags, which are based on a continuous-time expansion with Boltzmann weights that are Pfaffians of a matrix with dimension set by the expansion order, and which give the sign of the weight. Sign problem severity is related to the location of zero modes in the free hopping matrix, and for couplings $g\lesssim 0.7$ and $\beta=L/4$ it is fully removed for lattices up to at least $256$ sites. The interaction strength $g>0$ enhances the $U(1)$-charged dimerization doublet and suppresses the neutral channel, as the renormalization group predicts. Increasing the interactions also restores $U(1)$ symmetry in the component of the dimerization correlator that is even under bond orientation exchange, but not its odd component. There is a persistent anisotropy in the odd component with the symmetry of a lattice nematic. The data is consistent with a gapless semimetal throughout the sign-problem-minimal window.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-24: 1 request, 7,086 input tokens.

What the abstract itself states, as probabilities: claims first 0.21 · states limitation 0.81 · human participants 0.01 · controlled comparison 0.06 · randomized assignment 0.01 · states quantitative result 0.96 · compares with prior results 0.57 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

The two-dimensional Majorana-Hubbard model on the $\pi$-flux lattice involves the most local interaction possible for a lattice with one relativistic Majorana fermion per site. It is predicted to host a Majorana semimetal, dimerized phases, and quantum-critical points with an emergent $U(1)$ symmetry. Nonperturbative lattice calculations for this model remain a major open challenge, thus far limited to zero hopping. We present quantum Monte Carlo simulations of the model at finite hopping using Hamiltonian fermion bags, which are based on a continuous-time expansion with Boltzmann weights that are Pfaffians of a matrix with dimension set by the expansion order, and which give the sign of the weight. Sign problem severity is related to the location of zero modes in the free hopping matrix, and for couplings $g\lesssim 0.7$ and $\beta=L/4$ it is fully removed for lattices up to at least $256$ sites. The interaction strength $g>0$ enhances the $U(1)$-charged dimerization doublet and suppresses the neutral channel, as the renormalization group predicts. Increasing the interactions also restores $U(1)$ symmetry in the component of the dimerization correlator that is even under bond orientation exchange, but not its odd component. There is a persistent anisotropy in the odd component with the symmetry of a lattice nematic. The data is consistent with a gapless semimetal throughout the sign-problem-minimal window.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The two-dimensional Majorana-Hubbard model on the $\pi$-flux lattice involves the most local interaction possible for a lattice with one relativistic Majorana fermion per site. background 1.00 0.06 0.02 asserted 1.00
s1 It is predicted to host a Majorana semimetal, dimerized phases, and quantum-critical points with an emergent $U(1)$ symmetry. background 1.00 0.09 0.09 hedged possibility 0.55
s2 Nonperturbative lattice calculations for this model remain a major open challenge, thus far limited to zero hopping. gap 0.98 0.04 0.02 asserted 1.00
s3 We present quantum Monte Carlo simulations of the model at finite hopping using Hamiltonian fermion bags, which are based on a continuous-time expansion with Boltzmann weights that are Pfaffians of a matrix with dimension set by the expansion order, and which give the sign of the weight. method 0.74 0.96 0.07 asserted 1.00
s4 Sign problem severity is related to the location of zero modes in the free hopping matrix, and for couplings $g\lesssim 0.7$ and $\beta=L/4$ it is fully removed for lattices up to at least $256$ sites. result 0.90 0.87 0.11 asserted 0.99
s5 The interaction strength $g>0$ enhances the $U(1)$-charged dimerization doublet and suppresses the neutral channel, as the renormalization group predicts. result 0.80 0.69 0.18 asserted 1.00
s6 Increasing the interactions also restores $U(1)$ symmetry in the component of the dimerization correlator that is even under bond orientation exchange, but not its odd component. result 0.99 0.89 0.11 asserted 1.00
s7 There is a persistent anisotropy in the odd component with the symmetry of a lattice nematic. result 0.93 0.88 0.08 asserted 1.00
s8 The data is consistent with a gapless semimetal throughout the sign-problem-minimal window. result 0.69 0.91 0.49 hedged possibility 1.00

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 2. FIG. 2. (a) Channel ratio
Figure 3. FIG. 3. Connected correlators of the charged doublet
Figure 4. FIG. 4. With Φ =

Showing 3 of the paper's 4 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.26881 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.