pith. sign in

arxiv: 2603.11175 · v2 · pith:QV5I4K3Wnew · submitted 2026-03-11 · ❄️ cond-mat.mes-hall · cond-mat.str-el

Imaging flat band electron hydrodynamics in biased bilayer graphene

Pith reviewed 2026-05-21 11:38 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall cond-mat.str-el
keywords electron hydrodynamicsbilayer grapheneflat bandcurrent imagingSQUID microscopyBoltzmann transportdisplacement field tuning
0
0 comments X

The pith

In dual-gated bilayer graphene, hydrodynamic electron flow is strongest in the flat band regime with scattering lengths down to 50 nanometers.

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

This paper establishes that electron hydrodynamics in bilayer graphene can be made more robust by biasing the system into a flat band, where the effective mass is large. By mapping local current densities with a scanning superconducting sensor in a geometry that distinguishes different flow types, the authors map out ballistic, hydrodynamic, and diffusive regimes as functions of density and displacement field. The key finding is that the hydrodynamic regime reaches its peak strength when the band is flat, with the electron-electron scattering length shrinking to roughly the Fermi wavelength of 50 nm. This shortening allows hydrodynamic phenomena to persist at much smaller device scales than in conventional graphene.

Core claim

Imaging local current flow in biased bilayer graphene identifies three transport regimes across the carrier density and displacement field phase space, with the strongest hydrodynamic transport occurring in the flat band regime where the electron-electron scattering length is comparable to the Fermi wavelength of ~50 nm.

What carries the argument

Scanning superconducting magnetic sensor for local current imaging in a sample geometry sensitive to laminar and vortical flow, analyzed with a unified Boltzmann transport model.

If this is right

  • High current densities produce nonlinear flow patterns in the hydrodynamic regime.
  • The tunable flat band allows hydrodynamic transport to be observed at length scales around 50 nm.
  • Enhanced hydrodynamics via large effective mass supports development of compact hydrodynamic electronic devices.

Where Pith is reading between the lines

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

  • Flat band tuning could similarly enhance hydrodynamics in other rhombohedral multilayer graphenes.
  • Short scattering lengths may allow observation of hydrodynamic effects in nanoscale constrictions or vortices.
  • Nonlinear hydrodynamic flow could be harnessed for novel current rectification or amplification mechanisms.

Load-bearing premise

The sample geometry and the unified Boltzmann transport model allow clean separation of hydrodynamic flow from ballistic and diffusive contributions at the probed length scales.

What would settle it

Direct measurement showing that the electron-electron scattering length in the flat band remains much larger than 50 nm would falsify the enhanced hydrodynamics claim.

Figures

Figures reproduced from arXiv: 2603.11175 by Aidan Keough, Andrea F. Young, Andrew Lucas, Ania C. Bleszynski Jayich, Canxun Zhang, David Gong, Evgeny Redekop, Hari Stoyanov, Jack H. Farrell, Kenji Watanabe, Ludwig Holleis, Martin E. Huber, SungHoon Kim, Takashi Taniguchi, Youngjoon Choi.

Figure 1
Figure 1. Figure 1: ). Electron transport in the channel distinguishes hydrodynamic flow, characterized by a Poiseuille velocity profile with a higher current density in the center, from ballistic or diffusive behavior with flatter or even concave profiles[1, 6, 7, 11, 38]. In contrast, current flow in the lat￾eral chambers probes the local violation of Ohm’s law— stationary vortices or whirlpools can emerge in either the bal… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [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_p006_4.png] view at source ↗
read the original abstract

Hydrodynamic electron transport arises when carrier kinetics are dominated by interelectron collisions rather than the relaxation of momentum out of the electron system. In recent years, signatures of electron hydrodynamics have been reported in graphene devices owing to the low disorder and weak electron-phonon coupling. However, these experiments have been performed in regimes where the carrier mass is light, and the electron-electron collision length--though smaller than corresponding lengths for phonon or impurity scattering--remains large in absolute terms, typically several hundred nanometers. This restricts hydrodynamic transport phenomena to large length scales, limiting miniaturization of devices based on hydrodynamic flow. The advent of dual-gated rhombohedral graphene multilayers introduces a new route toward enhanced hydrodynamic behavior via their large--and tunable--effective mass. Here, we employ a scanning superconducting magnetic sensor to image local current flow in dual-gated bilayer graphene. Exploiting a sample geometry sensitive to both laminar and vortical flow, we identify three distinct transport regimes--ballistic, hydrodynamic, and diffusive--across the full phase space spanned by carrier density and displacement field. The strongest hydrodynamic transport is observed in the flat band regime, where fitting our results to a unified Boltzmann transport model reveals the electron-electron scattering length to be comparable to the Fermi wavelength of ~50 nm. High-current measurements, meanwhile, reveal striking nonlinearities in the flow pattern. Our results pave the way for miniaturized electronic devices based on linear and nonlinear electron hydrodynamics.

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 reports scanning SQUID imaging of local current flow in dual-gated bilayer graphene devices. By mapping transport across carrier density and displacement field, the authors identify three regimes—ballistic, hydrodynamic, and diffusive—using a geometry sensitive to both laminar and vortical flow. The strongest hydrodynamic signatures appear in the flat-band regime, where fitting local current images to a unified Boltzmann transport model yields an electron-electron scattering length l_ee comparable to the Fermi wavelength (~50 nm). High-current data additionally show nonlinear flow patterns.

Significance. If the model extraction is robust, the result establishes that flat-band systems support hydrodynamic transport at length scales an order of magnitude shorter than in conventional graphene, directly addressing the miniaturization limit noted in prior work. The direct visualization of flow patterns and the use of displacement-field tuning to access the flat band constitute clear experimental advances. The unified Boltzmann framework that simultaneously describes all three regimes is a methodological strength.

major comments (2)
  1. [Boltzmann model fitting] § on Boltzmann model fitting and parameter extraction: the central claim that l_ee reaches ~50 nm rests on the assumption that the dual-gated geometry and unified model cleanly isolate the hydrodynamic contribution from ballistic, diffusive, phonon, and disorder channels at the probed ~50 nm scales. Please supply the explicit functional form of the model (including how momentum-relaxing lengths enter) and a sensitivity analysis or covariance matrix showing that the fitted l_ee is not degenerate with other scattering lengths.
  2. [Flat-band regime results] Results section on flat-band regime: the large effective mass narrows the hydrodynamic window; demonstrate quantitatively (via additional simulations or temperature-dependent controls) that the observed current patterns cannot be reproduced by a purely ballistic or diffusive model with realistic edge scattering or residual disorder at the device length scales.
minor comments (2)
  1. [Figure captions] Figure captions should explicitly state the criteria (e.g., vorticity threshold or current-profile shape) used to assign each data point to ballistic, hydrodynamic, or diffusive regime.
  2. [Notation consistency] Notation for the displacement field D and carrier density n should be consistent between the abstract, main text, and figure axes.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major point below and have revised the manuscript to incorporate additional details and analyses as requested.

read point-by-point responses
  1. Referee: [Boltzmann model fitting] § on Boltzmann model fitting and parameter extraction: the central claim that l_ee reaches ~50 nm rests on the assumption that the dual-gated geometry and unified model cleanly isolate the hydrodynamic contribution from ballistic, diffusive, phonon, and disorder channels at the probed ~50 nm scales. Please supply the explicit functional form of the model (including how momentum-relaxing lengths enter) and a sensitivity analysis or covariance matrix showing that the fitted l_ee is not degenerate with other scattering lengths.

    Authors: We agree that explicit details on the model are important for robustness. The unified Boltzmann transport model solves the Boltzmann equation with a momentum-conserving electron-electron collision integral characterized by scattering length l_ee together with momentum-relaxing terms for impurity scattering (l_imp), phonon scattering (l_ph), and boundary scattering. The functional form is implemented numerically via a relaxation-time approximation augmented by the conserving collision operator, with the local current obtained by integrating the nonequilibrium distribution over the device geometry. In the revised manuscript we have added the explicit equations to the Methods section. We have also included a sensitivity analysis and the covariance matrix of the fits in the Supplementary Information; these confirm that l_ee remains well constrained and is not strongly degenerate with the other lengths at the ~50 nm scales probed by the SQUID sensor. revision: yes

  2. Referee: [Flat-band regime results] Results section on flat-band regime: the large effective mass narrows the hydrodynamic window; demonstrate quantitatively (via additional simulations or temperature-dependent controls) that the observed current patterns cannot be reproduced by a purely ballistic or diffusive model with realistic edge scattering or residual disorder at the device length scales.

    Authors: We appreciate the referee highlighting the effect of the large effective mass on the hydrodynamic window. To address this quantitatively, we have carried out additional numerical simulations of purely ballistic and diffusive transport using realistic edge scattering (specular-diffuse mixing) and disorder strengths consistent with the measured mobility and device dimensions. These simulations, now presented in the revised Supplementary Information, fail to reproduce the observed vortical flow patterns and the spatial current distributions measured in the flat-band regime. Our existing temperature-dependent data further support this distinction, as the hydrodynamic signatures persist in a manner inconsistent with temperature-independent ballistic or diffusive models. revision: yes

Circularity Check

0 steps flagged

Experimental imaging and Boltzmann model fitting yield l_ee ~ λ_F without derivation reducing to inputs by construction

full rationale

The paper's results derive from direct scanning SQUID imaging of local current in dual-gated bilayer graphene across density and displacement field, identifying ballistic/hydrodynamic/diffusive regimes via comparison to a standard unified Boltzmann transport model. The reported l_ee ≈ 50 nm in the flat-band regime is an output of fitting observed flow patterns to that model rather than a self-definitional or fitted-input-called-prediction step; the model is an independent theoretical framework with stated assumptions about scattering lengths, and the imaging data provide external falsifiability. No load-bearing self-citation chain, uniqueness theorem, or ansatz smuggling is required for the central claim, which retains independent experimental content.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the unified Boltzmann transport model for separating transport regimes and on the assumption that the scanning sensor accurately maps local current without significant back-action or calibration offsets.

free parameters (1)
  • electron-electron scattering length
    Extracted by fitting the observed flow patterns to the Boltzmann model; its value is not independently measured.
axioms (1)
  • domain assumption The sample geometry is sensitive to both laminar and vortical flow components.
    Invoked to distinguish hydrodynamic from ballistic signatures.

pith-pipeline@v0.9.0 · 5852 in / 1210 out tokens · 31915 ms · 2026-05-21T11:38:00.307877+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 3 Pith papers

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

  1. Cryogenic shock exfoliation for ultrahigh mobility rhombohedral graphite nanoelectronics

    cond-mat.mes-hall 2026-04 unverdicted novelty 8.0

    Cryogenic shock exfoliation yields large rhombohedral graphene devices over 1300 square micrometers with 90% fabrication yield, mean free path exceeding 200 micrometers, and signatures of electron hydrodynamics.

  2. Characterizing electronic scattering rates with transport in multiterminal devices

    cond-mat.mes-hall 2026-05 unverdicted novelty 6.0

    A five-terminal geometry diagnoses ballistic-hydrodynamic-Ohmic crossovers and extracts momentum-relaxing and conserving scattering rates from current partition in electron liquids.

  3. Characterizing electronic scattering rates with transport in multiterminal devices

    cond-mat.mes-hall 2026-05 unverdicted novelty 6.0

    Current partition in a five-terminal geometry diagnoses ballistic-hydrodynamic-Ohmic crossovers and extracts momentum-relaxing and conserving scattering rates in 2D electron systems.

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages · cited by 2 Pith papers

  1. [1]

    Gurzhi, R. N. Hydrodynamic effects in solids at low tem- perature.Soviet Physics Uspekhi11, 255 (1968)

  2. [2]

    A.et al.Negative local resistance caused by viscous electron backflow in graphene.Science351, 1055–1058 (2016)

    Bandurin, D. A.et al.Negative local resistance caused by viscous electron backflow in graphene.Science351, 1055–1058 (2016)

  3. [3]

    Krishna Kumar, R.et al.Superballistic flow of viscous electron fluid through graphene constrictions.Nature Physics13, 1182–1185 (2017)

  4. [4]

    A.et al.Fluidity onset in graphene.Nature Communications9, 4533 (2018)

    Bandurin, D. A.et al.Fluidity onset in graphene.Nature Communications9, 4533 (2018)

  5. [5]

    I.et al.Measuring hall viscosity of graphene’s electron fluid.Science364, 162–165 (2019)

    Berdyugin, A. I.et al.Measuring hall viscosity of graphene’s electron fluid.Science364, 162–165 (2019)

  6. [6]

    A.et al.Visualizing poiseuille flow of hydro- dynamic electrons.Nature576, 75–79 (2019)

    Sulpizio, J. A.et al.Visualizing poiseuille flow of hydro- dynamic electrons.Nature576, 75–79 (2019)

  7. [7]

    Ku, M. J. H.et al.Imaging viscous flow of the dirac fluid in graphene.Nature583, 537–541 (2020)

  8. [8]

    Kim, M.et al.Control of electron-electron interaction in graphene by proximity screening.Nature Communi- cations11, 2339 (2020)

  9. [9]

    Jenkins, A.et al.Imaging the breakdown of ohmic trans- port in graphene.Phys. Rev. Lett.129, 087701 (2022)

  10. [10]

    Kumar, C.et al.Imaging hydrodynamic electrons flowing without landauer–sharvin resistance.Nature609, 276– 281 (2022)

  11. [11]

    Huang, W.et al.Electronic poiseuille flow in hexagonal boron nitride encapsulated graphene field effect transis- tors.Phys. Rev. Res.5, 023075 (2023)

  12. [12]

    J.et al.Imaging the breaking of electrostatic dams in graphene for ballistic and viscous fluids.Science 379, 671–676 (2023)

    Krebs, Z. J.et al.Imaging the breaking of electrostatic dams in graphene for ballistic and viscous fluids.Science 379, 671–676 (2023)

  13. [13]

    L.et al.Observation of current whirlpools in graphene at room temperature.Science384, 465–469 (2024)

    Palm, M. L.et al.Observation of current whirlpools in graphene at room temperature.Science384, 465–469 (2024)

  14. [14]

    URLhttps://arxiv.org/abs/2406.13799

    Talanov, A.et al.Observation of electronic viscous dis- sipation in graphene magneto-thermal transport (2024). URLhttps://arxiv.org/abs/2406.13799. 2406.13799

  15. [15]

    & McCann, E

    Koshino, M. & McCann, E. Trigonal warping and berry’s phasenπin abc-stacked multilayer graphene.Phys. Rev. B80, 165409 (2009)

  16. [16]

    Nano Letters10, 1046–1049 (2010)

    Finkler, A.et al.Self-aligned nanoscale squid on a tip. Nano Letters10, 1046–1049 (2010)

  17. [17]

    Nature Nanotechnology8, 639–644 (2013)

    Vasyukov, D.et al.A scanning superconducting quantum interference device with single electron spin sensitivity. Nature Nanotechnology8, 639–644 (2013)

  18. [18]

    Aharon-Steinberg, A.et al.Direct observation of vortices in an electron fluid.Nature607, 74–80 (2022)

  19. [19]

    & Fong, K

    Lucas, A. & Fong, K. C. Hydrodynamics of electrons in graphene.Journal of Physics: Condensed Matter30, 053001 (2018)

  20. [20]

    & Scaffidi, T

    Fritz, L. & Scaffidi, T. Hydrodynamic electronic trans- port.Annual Review of Condensed Matter Physics15, 17–44 (2024)

  21. [21]

    N., Gornyi, I

    Narozhny, B. N., Gornyi, I. V., Titov, M., Sch¨ utt, M. & Mirlin, A. D. Hydrodynamics in graphene: Linear- response transport.Phys. Rev. B91, 035414 (2015)

  22. [22]

    & Falkovich, G

    Levitov, L. & Falkovich, G. Electron viscosity, current vortices and negative nonlocal resistance in graphene. Nature Physics12, 672–676 (2016)

  23. [23]

    & Vignale, G.Quantum Theory of the Elec- tron Liquid(Cambridge University Press, 2005)

    Giuliani, G. & Vignale, G.Quantum Theory of the Elec- tron Liquid(Cambridge University Press, 2005)

  24. [24]

    & MacDonald, A

    Jung, J. & MacDonald, A. H. Accurate tight-binding models for theπbands of bilayer graphene.Phys. Rev. B89, 035405 (2014)

  25. [25]

    Zhou, H.et al.Isospin magnetism and spin-polarized superconductivity in bernal bilayer graphene.Science 375, 774–778 (2022)

  26. [26]

    C.et al.Cascade of isospin phase transi- tions in bernal-stacked bilayer graphene at zero magnetic field.Nature Physics18, 771–775 (2022)

    de la Barrera, S. C.et al.Cascade of isospin phase transi- tions in bernal-stacked bilayer graphene at zero magnetic field.Nature Physics18, 771–775 (2022)

  27. [27]

    M.et al.Quantum cascade of correlated phases in trigonally warped bilayer graphene.Nature608, 298– 302 (2022)

    Seiler, A. M.et al.Quantum cascade of correlated phases in trigonally warped bilayer graphene.Nature608, 298– 302 (2022)

  28. [28]

    Zhang, Y.et al.Enhanced superconductivity in spin– orbit proximitized bilayer graphene.Nature613, 268– 273 (2023)

  29. [29]

    Li, C.et al.Tunable superconductivity in electron- and hole-doped bernal bilayer graphene.Nature631, 300– 306 (2024)

  30. [30]

    Holleis, L.et al.Nematicity and orbital depairing in superconducting bernal bilayer graphene.Nature Physics 21, 444–450 (2025). 8

  31. [31]

    Zhang, Y.et al.Twist-programmable superconductivity in spin–orbit-coupled bilayer graphene.Nature641, 625– 631 (2025)

  32. [32]

    Holleis, L.et al.Fluctuating magnetism and pomer- anchuk effect in multilayer graphene.Nature640, 355– 360 (2025)

  33. [33]

    M.et al.Interaction-driven quasi-insulating ground states of gapped electron-doped bilayer graphene

    Seiler, A. M.et al.Interaction-driven quasi-insulating ground states of gapped electron-doped bilayer graphene. Phys. Rev. Lett.133, 066301 (2024)

  34. [34]

    V., Kivelson, S

    Andreev, A. V., Kivelson, S. A. & Spivak, B. Hydro- dynamic description of transport in strongly correlated electron systems.Phys. Rev. Lett.106, 256804 (2011)

  35. [35]

    & Hartnoll, S

    Lucas, A. & Hartnoll, S. A. Kinetic theory of trans- port for inhomogeneous electron fluids.Phys. Rev. B 97, 045105 (2018)

  36. [36]

    V., Legvold, S

    Colvin, R. V., Legvold, S. & Spedding, F. H. Electrical resistivity of the heavy rare-earth metals.Phys. Rev. 120, 741–745 (1960)

  37. [37]

    Fisher, M. E. & Langer, J. S. Resistive anomalies at magnetic critical points.Phys. Rev. Lett.20, 665–668 (1968)

  38. [38]

    Torre, I., Tomadin, A., Geim, A. K. & Polini, M. Non- local transport and the hydrodynamic shear viscosity in graphene.Phys. Rev. B92, 165433 (2015)

  39. [39]

    Andreev, A. F. & Kosevich, Y. A. Kinetic phenomena in semiquantum liquids.Zh. Eksp. Teor. Fiz77, 2518–2523 (1979)

  40. [40]

    Meric, I.et al.Current saturation in zero-bandgap, top- gated graphene field-effect transistors.Nature Nanotech- nology3, 654–659 (2008)

  41. [41]

    & Tielrooij, K.- J

    Massicotte, M., Soavi, G., Principi, A. & Tielrooij, K.- J. Hot carriers in graphene – fundamentals and applica- tions.Nanoscale13, 8376–8411 (2021)

  42. [42]

    K.et al.Revealing electron-lattice de- coupling by peltier thermometry and nanoscale thermal imaging in graphene (2025)

    Srivastav, S. K.et al.Revealing electron-lattice de- coupling by peltier thermometry and nanoscale thermal imaging in graphene (2025). URLhttps://arxiv.org/ abs/2506.21523. 2506.21523

  43. [43]

    URLhttps: //arxiv.org/abs/2008.04862

    Geurs, J.et al.Rectification by hydrodynamic flow in an encapsulated graphene tesla valve (2020). URLhttps: //arxiv.org/abs/2008.04862. 2008.04862

  44. [44]

    Cook, C. Q. & Lucas, A. Electron hydrodynamics with a polygonal fermi surface.Phys. Rev. B99, 235148 (2019)

  45. [45]

    S., Anikeeva, P., Felser, C

    Varnavides, G., Jermyn, A. S., Anikeeva, P., Felser, C. & Narang, P. Electron hydrodynamics in anisotropic mate- rials.Nature Communications11, 4710 (2020)

  46. [46]

    Redekop, E.et al.Direct magnetic imaging of fractional chern insulators in twisted mote2.Nature635, 584–589 (2024)

  47. [47]

    Huber, M.et al.Dc squid series array amplifiers with 120 mhz bandwidth.IEEE Transactions on Applied Su- perconductivity(2001)

  48. [48]

    Feldmann, D. M. Resolution of two-dimensional cur- rents in superconductors from a two-dimensional mag- netic field measurement by the method of regularization. Phys. Rev. B69, 144515 (2004)

  49. [49]

    Y., Levin, E

    Meltzer, A. Y., Levin, E. & Zeldov, E. Direct recon- struction of two-dimensional currents in thin films from magnetic-field measurements.Phys. Rev. Appl.8, 064030 (2017)

  50. [50]

    L., Gross, E

    Bhatnagar, P. L., Gross, E. P. & Krook, M. A model for collision processes in gases. i. small amplitude processes in charged and neutral one-component systems.Physical Review94, 511–525 (1954)

  51. [51]

    Mendl, C. B. & Lucas, A. Dyakonov-shur instability across the ballistic-to-hydrodynamic crossover.Applied Physics Letters112, 124101 (2018)

  52. [52]

    & Levitov, L

    Ledwith, P., Guo, H., Shytov, A. & Levitov, L. To- mographic dynamics and scale-dependent viscosity in 2d electron systems.Physical Review Letters123, 116601 (2019)

  53. [53]

    Farrell, J. H. & Lucas, A. Simple devices that distin- guish viscous, ballistic, and diffusive transport.to appear (2026)

  54. [54]

    Farrell, J. H. FermiHarmonics: Linearized 2d fermi-liquid boltzmann transport on unstructured meshes.Zenodo (2026)

  55. [55]

    Proceedings of the JuliaCon Conferences1, 77 (2022)

    Ranocha, H.et al.Adaptive numerical simulations with Trixi.jl: A case study of Julia for scientific computing. Proceedings of the JuliaCon Conferences1, 77 (2022)

  56. [56]

    mock reconstructions

    Schlottke-Lakemper, M., Winters, A. R., Ranocha, H. & Gassner, G. J. A purely hyperbolic discontinu- ous Galerkin approach for self-gravitating gas dynamics. Journal of Computational Physics442, 110467 (2021). METHODS Sample fabrication Dual-gated R2G was fabricated through van der Waals stacking of bilayer graphene, hexagonal boron nitride (hBN) dielectr...