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
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
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
- 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
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.