Imaging flat band electron hydrodynamics in biased bilayer graphene
Pith reviewed 2026-05-21 11:38 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [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.
- [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
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
-
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
-
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
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
free parameters (1)
- electron-electron scattering length
axioms (1)
- domain assumption The sample geometry is sensitive to both laminar and vortical flow components.
Forward citations
Cited by 3 Pith papers
-
Cryogenic shock exfoliation for ultrahigh mobility rhombohedral graphite nanoelectronics
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.
-
Characterizing electronic scattering rates with transport in multiterminal devices
A five-terminal geometry diagnoses ballistic-hydrodynamic-Ohmic crossovers and extracts momentum-relaxing and conserving scattering rates from current partition in electron liquids.
-
Characterizing electronic scattering rates with transport in multiterminal devices
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
-
[1]
Gurzhi, R. N. Hydrodynamic effects in solids at low tem- perature.Soviet Physics Uspekhi11, 255 (1968)
work page 1968
-
[2]
Bandurin, D. A.et al.Negative local resistance caused by viscous electron backflow in graphene.Science351, 1055–1058 (2016)
work page 2016
-
[3]
Krishna Kumar, R.et al.Superballistic flow of viscous electron fluid through graphene constrictions.Nature Physics13, 1182–1185 (2017)
work page 2017
-
[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)
work page 2018
-
[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)
work page 2019
-
[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)
work page 2019
-
[7]
Ku, M. J. H.et al.Imaging viscous flow of the dirac fluid in graphene.Nature583, 537–541 (2020)
work page 2020
-
[8]
Kim, M.et al.Control of electron-electron interaction in graphene by proximity screening.Nature Communi- cations11, 2339 (2020)
work page 2020
-
[9]
Jenkins, A.et al.Imaging the breakdown of ohmic trans- port in graphene.Phys. Rev. Lett.129, 087701 (2022)
work page 2022
-
[10]
Kumar, C.et al.Imaging hydrodynamic electrons flowing without landauer–sharvin resistance.Nature609, 276– 281 (2022)
work page 2022
-
[11]
Huang, W.et al.Electronic poiseuille flow in hexagonal boron nitride encapsulated graphene field effect transis- tors.Phys. Rev. Res.5, 023075 (2023)
work page 2023
-
[12]
Krebs, Z. J.et al.Imaging the breaking of electrostatic dams in graphene for ballistic and viscous fluids.Science 379, 671–676 (2023)
work page 2023
-
[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)
work page 2024
-
[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]
Koshino, M. & McCann, E. Trigonal warping and berry’s phasenπin abc-stacked multilayer graphene.Phys. Rev. B80, 165409 (2009)
work page 2009
-
[16]
Nano Letters10, 1046–1049 (2010)
Finkler, A.et al.Self-aligned nanoscale squid on a tip. Nano Letters10, 1046–1049 (2010)
work page 2010
-
[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)
work page 2013
-
[18]
Aharon-Steinberg, A.et al.Direct observation of vortices in an electron fluid.Nature607, 74–80 (2022)
work page 2022
- [19]
-
[20]
Fritz, L. & Scaffidi, T. Hydrodynamic electronic trans- port.Annual Review of Condensed Matter Physics15, 17–44 (2024)
work page 2024
-
[21]
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)
work page 2015
-
[22]
Levitov, L. & Falkovich, G. Electron viscosity, current vortices and negative nonlocal resistance in graphene. Nature Physics12, 672–676 (2016)
work page 2016
-
[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)
work page 2005
-
[24]
Jung, J. & MacDonald, A. H. Accurate tight-binding models for theπbands of bilayer graphene.Phys. Rev. B89, 035405 (2014)
work page 2014
-
[25]
Zhou, H.et al.Isospin magnetism and spin-polarized superconductivity in bernal bilayer graphene.Science 375, 774–778 (2022)
work page 2022
-
[26]
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)
work page 2022
-
[27]
Seiler, A. M.et al.Quantum cascade of correlated phases in trigonally warped bilayer graphene.Nature608, 298– 302 (2022)
work page 2022
-
[28]
Zhang, Y.et al.Enhanced superconductivity in spin– orbit proximitized bilayer graphene.Nature613, 268– 273 (2023)
work page 2023
-
[29]
Li, C.et al.Tunable superconductivity in electron- and hole-doped bernal bilayer graphene.Nature631, 300– 306 (2024)
work page 2024
-
[30]
Holleis, L.et al.Nematicity and orbital depairing in superconducting bernal bilayer graphene.Nature Physics 21, 444–450 (2025). 8
work page 2025
-
[31]
Zhang, Y.et al.Twist-programmable superconductivity in spin–orbit-coupled bilayer graphene.Nature641, 625– 631 (2025)
work page 2025
-
[32]
Holleis, L.et al.Fluctuating magnetism and pomer- anchuk effect in multilayer graphene.Nature640, 355– 360 (2025)
work page 2025
-
[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)
work page 2024
-
[34]
Andreev, A. V., Kivelson, S. A. & Spivak, B. Hydro- dynamic description of transport in strongly correlated electron systems.Phys. Rev. Lett.106, 256804 (2011)
work page 2011
-
[35]
Lucas, A. & Hartnoll, S. A. Kinetic theory of trans- port for inhomogeneous electron fluids.Phys. Rev. B 97, 045105 (2018)
work page 2018
-
[36]
Colvin, R. V., Legvold, S. & Spedding, F. H. Electrical resistivity of the heavy rare-earth metals.Phys. Rev. 120, 741–745 (1960)
work page 1960
-
[37]
Fisher, M. E. & Langer, J. S. Resistive anomalies at magnetic critical points.Phys. Rev. Lett.20, 665–668 (1968)
work page 1968
-
[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)
work page 2015
-
[39]
Andreev, A. F. & Kosevich, Y. A. Kinetic phenomena in semiquantum liquids.Zh. Eksp. Teor. Fiz77, 2518–2523 (1979)
work page 1979
-
[40]
Meric, I.et al.Current saturation in zero-bandgap, top- gated graphene field-effect transistors.Nature Nanotech- nology3, 654–659 (2008)
work page 2008
-
[41]
Massicotte, M., Soavi, G., Principi, A. & Tielrooij, K.- J. Hot carriers in graphene – fundamentals and applica- tions.Nanoscale13, 8376–8411 (2021)
work page 2021
-
[42]
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]
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]
Cook, C. Q. & Lucas, A. Electron hydrodynamics with a polygonal fermi surface.Phys. Rev. B99, 235148 (2019)
work page 2019
-
[45]
Varnavides, G., Jermyn, A. S., Anikeeva, P., Felser, C. & Narang, P. Electron hydrodynamics in anisotropic mate- rials.Nature Communications11, 4710 (2020)
work page 2020
-
[46]
Redekop, E.et al.Direct magnetic imaging of fractional chern insulators in twisted mote2.Nature635, 584–589 (2024)
work page 2024
-
[47]
Huber, M.et al.Dc squid series array amplifiers with 120 mhz bandwidth.IEEE Transactions on Applied Su- perconductivity(2001)
work page 2001
-
[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)
work page 2004
-
[49]
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)
work page 2017
-
[50]
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)
work page 1954
-
[51]
Mendl, C. B. & Lucas, A. Dyakonov-shur instability across the ballistic-to-hydrodynamic crossover.Applied Physics Letters112, 124101 (2018)
work page 2018
-
[52]
Ledwith, P., Guo, H., Shytov, A. & Levitov, L. To- mographic dynamics and scale-dependent viscosity in 2d electron systems.Physical Review Letters123, 116601 (2019)
work page 2019
-
[53]
Farrell, J. H. & Lucas, A. Simple devices that distin- guish viscous, ballistic, and diffusive transport.to appear (2026)
work page 2026
-
[54]
Farrell, J. H. FermiHarmonics: Linearized 2d fermi-liquid boltzmann transport on unstructured meshes.Zenodo (2026)
work page 2026
-
[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)
work page 2022
-
[56]
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...
work page 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.