Pith. sign in

REVIEW 6 minor 51 references

Vacancy-assisted superfluid drag

T0 review · 0 major / 6 minor · reviewed 2026-08-09 · deepseek-v4-flash

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

Pith's one-line read Single hole fixes superfluid drag at exactly 1−2/π

desk verdict Exact single-hole drag coefficient κ = 1 − 2/π in the hard-core two-component Bose-Hubbard model, backed by two independent numerics; the paper is solid and deserves serious refereeing. read the letter →

arxiv 2502.00542 v2 pith:IWLPWXHK submitted 2025-02-01 cond-mat.quant-gas cond-mat.str-el

Thomas G. Kiely , Chao Zhang , Erich J. Mueller This is my paper · ORCID
classification cond-mat.quant-gascond-mat.str-el
keywords superfluiddragAndreev-Bashkineffecttwo-componentBose-Hubbardmodelhard-corebosonskineticmagnetismspinpolaronquantumMonteCarlomatrixproductstates
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

This paper establishes an exact value for the superfluid drag coefficient in the two-component Bose-Hubbard model with infinitely strong repulsion on the square lattice. In the limit of vanishing hole density, a single mobile hole produces a strong dissipationless coupling between the currents of the two components, and the drag coefficient is exactly κ = 1 − 2/π. The authors derive this closed form by solving the single-hole Schrödinger equation with a variational ansatz that allows one flipped spin, and they verify it with infinite matrix product state calculations on cylinders and quantum Monte Carlo on finite tori. The result turns a strongly correlated many-body response into a solvable single-hole problem, and it predicts a measurable spin polaron around each hole.

What carries the argument

The central object is the single-hole variational wavefunction |ψ⟩ = (1/√N_s) ∑_j b_{j+} [ f_0 + ∑_s f_s b†_{j+s,−} b_{j+s,+} ] |ϕ⟩, which keeps exactly one spin flip relative to the fully polarized XY ferromagnet. Solving the resulting Schrödinger equations for the amplitudes f_0 and f_s reduces the drag coefficient to a Brillouin-zone integral, κ = ∫ d²p/(2π)² (cos 2p_x − 1)/(−4 + 2 cos p_x + 2 cos p_y), which evaluates to 1 − 2/π. The mechanism that carries the argument is correlated hopping: a flipped spin adjacent to the hole can hop over it, and this correlated process couples the density and spin channels, producing the drag.

What would settle it

A high-accuracy calculation of the energy change under a small counterflow twist at very low hole density on a two-dimensional square lattice would settle the claim: if the drag coefficient extrapolated to zero hole density disagrees with 1 − 2/π beyond finite-size error, the one-spin-flip ansatz is wrong. Alternatively, a quantum-gas-microscope measurement of the spin twist around a hole that does not show the predicted long-distance tail Υ_j ∝ (j_x q_x + j_y q_y)/|j|² would cast doubt on the mechanism.

Watch

Extended reading notes

Core claim

On its own terms, the paper claims that the superfluid density matrix of the hard-core two-component Bose-Hubbard model on a square lattice has eigenvalues ρ_+ = x_v and ρ_− = x_v/(π − 1) in the limit x_v → 0, giving a drag coefficient κ = (ρ_+ − ρ_−)/(ρ_+ + ρ_−) = 1 − 2/π ≈ 0.363. The argument is exact for a single hole: after a gauge transformation that imposes phase twists, the current-carrying ground state is an XY ferromagnet dressed by one flipped spin, and the correlated hopping of that spin with the hole depletes the spin-channel stiffness. The same calculation gives κ_1D = 1 on a chain, κ ≈ 0.23 on the triangular lattice, and κ ≈ 1/(z − 1) on hypercubic lattices, so the drag weakens as the coordination number grows. The authors verify the square-lattice value against tensor-network and Monte Carlo simulations at finite hole density, finding agreement in the extrapolated zero-doping limit.

Load-bearing premise

The derivation rests on the assumption that the wavefunction needs at most one flipped spin; if a second flipped spin changes the energy at the same order as the counterflow twist, the closed-form value of the drag coefficient would be different.

Editorial extensions

If this is right

  • In the dilute-hole limit, counterflow currents are strongly suppressed relative to mass currents, so moving one component inevitably drags the other along.
  • The spin-channel superfluid density at low doping is ρ_− = x_v/(π − 1), a sharp quantitative prediction that deviates from the mean-field expectation ρ_+ = x_v.
  • The drag coefficient sets the magnitude and sign of vortex–vortex interactions between the two species, so a drag of order 0.36 should produce readily observable vortex coupling.
  • The single-hole calculation is exact for any regular lattice, yielding testable predictions for chains, hypercubic lattices, and the triangular lattice.
  • The spin polaron around a hole, with twist decaying as Υ_j ∝ (j_x q_x + j_y q_y)/|j|² at long distances, is observable with quantum gas microscopes.

Reading between the lines

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

  • The single-hole solution could serve as the zeroth order of a controlled expansion in hole density, giving a systematic analytic route to finite-density drag beyond numerics.
  • An exact two-dimensional transport coefficient of this kind is rare, so the closed form provides a benchmark against which tensor-network and quantum Monte Carlo codes for strongly correlated lattice models can be tested.
  • Measuring the hole-spin correlator in a quantum gas microscope could provide a model-independent test of the drag mechanism, since the predicted polaron shape encodes the same physics as κ.
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

0 major / 6 minor

Summary. This paper studies superfluid drag in the two-component hard-core Bose-Hubbard model in the limit of infinite on-site repulsion. On the square lattice, the authors derive an exact expression for the drag coefficient in the vanishing-hole-density limit, κ = 1 − 2/π (Eqs. 21–23), by solving the single-hole Schrödinger equation in the presence of a phase twist to second order in the twist. They also obtain exact values for finite cylinders and other lattices. The analytic result is compared with two independent numerical methods: infinite cylinder iMPS (VUMPS) calculations and quantum Monte Carlo on finite tori, finding good agreement. The paper further characterizes the spin-hole correlations (spin polaron) and discusses experimental realization in cold atoms, Rydberg arrays, and transmon arrays.

Significance. The central analytic result is a rare exact, parameter-free prediction for a strongly correlated two-component lattice superfluid: the drag coefficient is fixed by the Schrödinger equation without fitting, and the same calculation yields explicit finite-cylinder values that match the numerical data. The paper's convergence of three methods—analytic perturbation theory, iMPS, and QMC—gives strong support to the claim. The result is falsifiable and has clear experimental relevance through the spin-hole correlator Υ_j. The work also provides a useful general picture: in the dilute-hole limit, the spin channel stiffness is strongly reduced by correlated-hopping processes, giving an order-unity drag coefficient.

minor comments (6)
  1. [Appendix A, Eq. (A3)] In the definition of Π^y_s, the second term should read Λ_{s_x, s_y − 1} rather than Λ_{s_x, s_y + 1}; the following equations use the antisymmetric combination, so this is a typographical error.
  2. [Appendix A, text after Eq. (A7)] The intermediate expression E = ϵ_0 − 2q^2t Π^x_x/(1 − Π^x_x) should have a plus sign in the denominator (1 + Π^x_x); the subsequent expression for E/t in Eq. (A8) is consistent with the plus sign, so the minus sign is a typo.
  3. [Sec. V, around Eq. (7)] The phrase 'fs,0 ∝ q is small' is misleading: the no-flip amplitude f0 is order one (f0 = 1 + O(q^2)), while the flip amplitudes f_s with s≠0 are O(q). Please reword to avoid confusion.
  4. [Sec. V, finite-cylinder results] The reported values for C = 3, 4, 5 should be written with parentheses to avoid ambiguity, e.g., κ = 2 − √(7/3) for C = 3; without parentheses, 2 − √7/3 could be read as 2 − (√7)/3.
  5. [Sec. VIII and Appendix B] The C→∞ extrapolation of ρ_−(C) rests on only three cylinder widths (C = 3, 4, 5) and an assumed 1/C^2 scaling; an estimate of the extrapolation uncertainty would be helpful, although the central analytic claim does not rely on this extrapolation.
  6. [Sec. V, generalization to other lattices] The statement that κ_nD ∼ 1/(z−1) on hypercubic lattices is given without derivation; a brief derivation or a more explicit statement of the asymptotic regime would improve the presentation.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the analytic κ=1−2/π is derived from the single-hole Schrödinger equation without fitting, and the numerical benchmarks are independent.

full rationale

The paper's central analytic result is self-contained. The ansatz Eq. (7) restricts the wavefunction to zero or one flipped spin relative to the fully polarized XY ferromagnet, but this is controlled by the perturbation structure: the spin-flip terms in Eq. (6) are linear in q and change the spin-flip number by ±1, so the O(q^2) energy shift is exactly given by one-flip intermediate states; multi-flip amplitudes enter only at O(q^4). Equations (12)–(13) are the complete O(q^2) Schrödinger equations, and Appendix A solves them without fitting, yielding the integral κ = ∫ (cos 2p_x − 1)/(−4 + 2 cos p_x + 2 cos p_y) = 1 − 2/π. The identification κ = Π^x_x in Eq. (A8) is algebraic comparison of the q^2 coefficient with the definition of κ; it is not a fit. The finite-density iMPS extrapolation uses ρ_−(C) = ρ_−∞ + α/C^2 with C = 3, 4, 5, but this α is a fitting parameter only for the finite-density curve, not for the x_v→0 limit, which is separately verified by the analytic stars and by QMC agreement. Self-citations [32,36] describe the VUMPS methodology and a gauge-invariance warning; they are not used to establish the value of κ or to rule out alternatives. Limitations (small-C extrapolation, QMC small-x_v) are acknowledged but do not reduce the central claim to its inputs.

Assumptions & free parameters 1 free parameters · 4 assumptions · 1 invented entities

The exact single-hole result uses no fitted constants; it is a parameter-free derivation from the hard-core Hamiltonian. The trusted inputs are Nagaoka and kinetic magnetism, the standard stiffness definition, and a one-spin-flip truncation that is justified only to second order in q. The finite-density numeric extrapolation introduces one fitted coefficient (alpha in the 1/C^2 form). No new microscopic entities are added; the spin polaron is a derived correlation with a falsifiable spatial profile.

free parameters (1)
  • Extrapolation coefficient alpha in rho_minus(C) = rho_minus_inf + alpha/C^2 = not quoted in paper
    Chosen to extrapolate cylinder iMPS data to C -> infinity in Figs. 1 and 2. It is a numerical fitting parameter and is not used in the exact single-hole formula for kappa.
assumptions (4)
  • domain assumption For a single hole at t/U=0, the ground state at zero twist is the fully polarized XY ferromagnet |psi0> = (1/sqrt(Ns)) sum_j b_{j+} |phi>.
    Invoked in Sec. V and relies on Nagaoka and kinetic magnetism [16,17]. The entire perturbative expansion is built around this state.
  • domain assumption To second order in the phase twist q, only single-spin-flip configurations are needed; states with more flipped spins are suppressed by higher powers of q.
    Used in the ansatz Eq. (7) and Appendix A. Standard perturbation counting, but the paper does not give a rigorous bound.
  • standard math The superfluid density tensor is extracted from the constrained minimization of the gauge-transformed Hamiltonian over translationally invariant states (Sec. IV).
    This is the standard definition of the stiffness matrix and is stated in Eq. (5).
  • domain assumption The t/U -> 0 limit replaces interactions by a hard-core constraint and neglects the superexchange J = 4 t^2/U (Sec. III).
    Model definition; the paper argues this limit is physically relevant for t/U << 1.
invented entities (1)
  • Spin polaron (a hole dressed by a cloud of flipped spins) independent evidence
    purpose: Explains the reduced spin-channel stiffness and the spatial dependence of the spin twist near a hole.
    The polaron is derived from the wavefunction, not postulated as a new degree of freedom. It predicts the spin-hole correlator Upsilon_j ~ (j_x q_x + j_y q_y)/|j|^2, which the paper argues is measurable in quantum gas microscopes.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Vacancy-assisted superfluid drag." pith.science (2026). https://pith.science/paper/IWLPWXHK

@misc{pith2026250200542,
  author       = {Pith},
  title        = {Pith review of: Vacancy-assisted superfluid drag},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IWLPWXHK}},
  note         = {Machine review of arXiv:2502.00542}
}
read the original abstract

We study superfluid drag in the two-component Bose-Hubbard model with infinitely strong repulsive interactions. In this system, all transport is mediated by the motion of empty sites, or ``holes", and it is hard to move one component without moving the other. We demonstrate, with a combination of analytic and numeric techniques, that the motion of holes leads to strong dissipationless coupling between currents in the two components. This behavior is attributable to polaronic correlations that emerge in the presence of spin currents, which can be observed in experiments. We derive a closed-form expression for the coupling on various lattices in arbitrary spatial dimensions, which we verify through numerical simulations on two dimensional lattices.

Figures

Figures reproduced from arXiv: 2502.00542 by the authors.

Figure 1
Figure 1. shows the superfluid fractions ρ±, normalized by the mean-field superfluid density ρ+,MF = xν(1−xν), as a func￾tion of total hole density for iMPS calculations on cylinders of circumference C = 3, 4, 5 (blue). Stars of the associated color show the perturbative results at infinitesimal doping, which are in close agreement with the iMPS calculations on every cylindrical geometry shown. Black stars show the perturba￾t… view at source ↗
Figure 2
Figure 2. FIG. 2. Superfluid drag coe [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 4
Figure 4. shows κ for square L×L lattices with L = 10, 16, and 20 at a inverse temperature β = L. As described in Sec. VII, the drag coefficient is calculated from winding number statis￾tics. In these calculations we find that κ grows with L. This is opposite to the C dependence of the zero temperature iMPS calculations shown in [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figures from the paper (1 more)
Figure 5
Figure 5. Figure 5: FIG. 5. Superfluid drag coe [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

51 extracted references · 40 canonical work pages

  1. [1]

    A. F. Andreev and E. P. Bashkin, Soviet Journal of Experimen- tal and Theoretical Physics 42, 164 (1975)

    work page 1975

  2. [2]

    D. V . Fil and S. I. Shevchenko, Phys. Rev. A72, 013616 (2005)

    work page 2005

  3. [3]

    D. V . Fil and S. I. Shevchenko, Low Temperature Physics 30, 770–777 (2004)

    work page 2004

  4. [4]

    Nespolo, G

    J. Nespolo, G. E. Astrakharchik, and A. Recati, New Journal of 9 0 0.25 0.50 0.1 0.2 0.3 0 0.05 0.1 0.15 0.2 FIG. 5. Superfluid drag coe fficientκ as a function of hole density xν from quantum Monte Carlo (QMC) simulations on a square lattice with linear system size L = 20. Data are shown for inverse temper- atures β = L/2 (green), β = L (blue), and β = 2...

    work page 2017

  5. [5]

    Karle, N

    V . Karle, N. Defenu, and T. Enss, Phys. Rev. A 99, 063627 (2019)

    work page 2019

  6. [6]

    Romito, C

    D. Romito, C. Lobo, and A. Recati, Phys. Rev. Res. 3, 023196 (2021)

    work page 2021

  7. [7]

    V . M. Kaurov, A. B. Kuklov, and A. E. Meyerovich, Phys. Rev. Lett. 95, 090403 (2005)

    work page 2005

  8. [8]

    Sellin and E

    K. Sellin and E. Babaev, Phys. Rev. B 97, 094517 (2018)

    work page 2018

Show all 51 references
  1. [9]

    Kuklov, N

    A. Kuklov, N. Prokof’ev, and B. Svistunov, Phys. Rev. Lett.92, 030403 (2004)

    2004

  2. [10]

    Kuklov, N

    A. Kuklov, N. Prokof’ev, and B. Svistunov, Phys. Rev. Lett.92, 050402 (2004)

    2004

  3. [11]

    Isacsson, M.-C

    A. Isacsson, M.-C. Cha, K. Sengupta, and S. M. Girvin, Phys. Rev. B 72, 184507 (2005)

    2005

  4. [12]

    A. Hu, L. Mathey, I. Danshita, E. Tiesinga, C. J. Williams, and C. W. Clark, Phys. Rev. A80, 023619 (2009)

    2009

  5. [13]

    A. B. Kuklov and B. V . Svistunov, Phys. Rev. Lett.90, 100401 (2003)

    2003

  6. [14]

    Altman, W

    E. Altman, W. Hofstetter, E. Demler, and M. D. Lukin, New Journal of Physics 5, 113 (2003)

    2003

  7. [15]

    Hubener, M

    A. Hubener, M. Snoek, and W. Hofstetter, Phys. Rev. B 80, 245109 (2009)

    2009

  8. [16]

    Nagaoka, Phys

    Y . Nagaoka, Phys. Rev.147, 392 (1966)

    1966

  9. [17]

    Eisenberg and E

    E. Eisenberg and E. H. Lieb, Phys. Rev. Lett. 89, 220403 (2002)

    2002

  10. [18]

    Parisi, G

    L. Parisi, G. E. Astrakharchik, and S. Giorgini, Phys. Rev. Lett. 121, 025302 (2018)

    2018

  11. [19]

    E. K. Dahl, E. Babaev, and A. Sudbø, Phys. Rev. Lett. 101, 255301 (2008)

    2008

  12. [20]

    M. F. Parsons, A. Mazurenko, C. S. Chiu, G. Ji, D. Greif, and M. Greiner, Science 353, 1253–1256 (2016)

    2016

  13. [21]

    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, Science 353, 1260–1264 (2016)

    2016

  14. [22]

    M. Boll, T. A. Hilker, G. Salomon, A. Omran, J. Nespolo, L. Pollet, I. Bloch, and C. Gross, Science 353, 1257–1260 (2016)

    2016

  15. [23]

    Koepsell, J

    J. Koepsell, J. Vijayan, P. Sompet, F. Grusdt, T. A. Hilker, E. Demler, G. Salomon, I. Bloch, and C. Gross, Nature 572, 358–362 (2019)

    2019

  16. [24]

    Koepsell, D

    J. Koepsell, D. Bourgund, P. Sompet, S. Hirthe, A. Bohrdt, Y . Wang, F. Grusdt, E. Demler, G. Salomon, C. Gross, and I. Bloch, Science 374, 82–86 (2021)

    2021

  17. [25]

    Bohrdt, Y

    A. Bohrdt, Y . Wang, J. Koepsell, M. K´anasz-Nagy, E. Demler, and F. Grusdt, Phys. Rev. Lett.126, 026401 (2021)

    2021

  18. [26]

    Linder and A

    J. Linder and A. Sudbø, Phys. Rev. A 79, 063610 (2009)

    2009

  19. [27]

    P. P. Hofer, C. Bruder, and V . M. Stojanovi´c, Phys. Rev. A 86, 033627 (2012)

    2012

  20. [28]

    Hartman, E

    S. Hartman, E. Erlandsen, and A. Sudbø, Phys. Rev. B 98, 024512 (2018)

    2018

  21. [29]

    V . E. Colussi, F. Cale ffi, C. Menotti, and A. Recati, SciPost Phys. 12, 111 (2022)

    2022

  22. [30]

    Contessi, D

    D. Contessi, D. Romito, M. Rizzi, and A. Recati, Phys. Rev. Res. 3, L022017 (2021)

    2021

  23. [31]

    Gr ´emaud and G

    B. Gr ´emaud and G. G. Batrouni, Phys. Rev. Lett. 127, 025301 (2021)

    2021

  24. [32]

    T. G. Kiely and E. J. Mueller, Phys. Rev. B105, 134502 (2022)

    2022

  25. [33]

    Schollw ¨ock, Annals of Physics 326, 96–192 (2011)

    U. Schollw ¨ock, Annals of Physics 326, 96–192 (2011)

    2011

  26. [34]

    Liang and H

    S. Liang and H. Pang, Phys. Rev. B 49, 9214 (1994)

    1994

  27. [35]

    Zauner-Stauber, L

    V . Zauner-Stauber, L. Vanderstraeten, M. T. Fishman, F. Ver- straete, and J. Haegeman, Phys. Rev. B 97, 045145 (2018)

    2018

  28. [36]

    T. G. Kiely and E. J. Mueller, Phys. Rev. B106, 235126 (2022)

    2022

  29. [37]

    Pollmann, S

    F. Pollmann, S. Mukerjee, A. M. Turner, and J. E. Moore, Phys. Rev. Lett. 102, 255701 (2009)

    2009

  30. [38]

    Prokof’ev, B

    N. Prokof’ev, B. Svistunov, and I. Tupitsyn, Journal of Experi- mental and Theoretical Physics 87, 310 (1998)

    1998

  31. [39]

    Prokof’ev, B

    N. Prokof’ev, B. Svistunov, and I. Tupitsyn, Physics Letters A 238, 253 (1998)

    1998

  32. [40]

    Capogrosso-Sansone, G

    B. Capogrosso-Sansone, G. S ¨oyler, N. V . Prokof’ev, and B. V . Svistunov, Phys. Rev. A81, 053622 (2010)

    2010

  33. [41]

    D. M. Ceperley and E. L. Pollock, Phys. Rev. B 39, 2084 (1989)

    1989

  34. [42]

    F. m. c. Cr ´epin, N. Laflorencie, G. Roux, and P. Simon, Phys. Rev. B 84, 054517 (2011)

    2011

  35. [43]

    M. Qiao, G. Emperauger, C. Chen, L. Homeier, S. Hollerith, G. Bornet, R. Martin, B. G ´ely, L. Klein, D. Barredo, S. Geier, N.-C. Chiu, F. Grusdt, A. Bohrdt, T. Lahaye, and A. Browaeys, Realization of a doped quantum antiferromagnet with dipolar tunnelings in a rydberg tweezer...

    2025 arXiv

  36. [44]

    Gadway, D

    B. Gadway, D. Pertot, R. Reimann, and D. Schneble, Phys. Rev. Lett. 105, 045303 (2010)

    2010

  37. [45]

    Carlini and S

    F. Carlini and S. Stringari, Phys. Rev. A 104, 023301 (2021)

    2021

  38. [46]

    Zheng, A

    Y .-G. Zheng, A. Luo, Y .-C. Shen, M.-G. He, Z.-H. Zhu, Y . Liu, W.-Y . Zhang, H. Sun, Y . Deng, Z.-S. Yuan, and J.-W. Pan, Na- ture Physics 10.1038/s41567-024-02732-5 (2025)

    2025 doi

  39. [47]

    Sekino, H

    Y . Sekino, H. Tajima, and S. Uchino, Phys. Rev. Res.5, 023058 (2023)

    2023

  40. [48]

    Pradhan, R

    N. Pradhan, R. Kanamoto, M. Bhattacharya, and P. K. Mishra, Signature of andreev-bashkin superfluid drag from cavity op- tomechanics (2024), arXiv:2410.21015 [cond-mat.quant-gas]

    2024 arXiv

  41. [49]

    Homeier, T

    L. Homeier, T. J. Harris, T. Blatz, S. Geier, S. Hollerith, U. Schollw¨ock, F. Grusdt, and A. Bohrdt, Phys. Rev. Lett. 132, 230401 (2024)

    2024

  42. [50]

    Yanay, J

    Y . Yanay, J. Braum ¨uller, S. Gustavsson, W. D. Oliver, and C. Tahan, npj Quantum Information 6, 10.1038 /s41534-020- 0269-1 (2020)

    2020

  43. [51]

    T. J. Harris, U. Schollw ¨ock, A. Bohrdt, and F. Grusdt, Kinetic magnetism and stripe order in the doped afm bosonic t− J model (2024), arXiv:2410.00904 [cond-mat.quant-gas]

    2024 arXiv

Pith tools

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