Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Fizeau drag in graphene plasmonics


Dragging of light by moving media was predicted by Fresnel1 and verified by Fizeau’s celebrated experiments2 with flowing water. This momentous discovery is among the experimental cornerstones of Einstein’s special relativity theory and is well understood3,4 in the context of relativistic kinematics. By contrast, experiments on dragging photons by an electron flow in solids are riddled with inconsistencies and have so far eluded agreement with the theory5,6,7. Here we report on the electron flow dragging surface plasmon polaritons8,9 (SPPs): hybrid quasiparticles of infrared photons and electrons in graphene. The drag is visualized directly through infrared nano-imaging of propagating plasmonic waves in the presence of a high-density current. The polaritons in graphene shorten their wavelength when propagating against the drifting carriers. Unlike the Fizeau effect for light, the SPP drag by electrical currents defies explanation by simple kinematics and is linked to the nonlinear electrodynamics of Dirac electrons in graphene. The observed plasmonic Fizeau drag enables breaking of time-reversal symmetry and reciprocity10 at infrared frequencies without resorting to magnetic fields11,12 or chiral optical pumping13,14. The Fizeau drag also provides a tool with which to study interactions and nonequilibrium effects in electron liquids.

Access options

Rent or Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: Plasmonic Fizeau drag in graphene: theory and modelling.
Fig. 2: Experimental demonstration of the plasmonic Fizeau drag.
Fig. 3: Quantitative analysis of the plasmonic Fizeau drag.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Code availability

The code used to analyse data are available from the corresponding author upon reasonable request.


  1. 1.

    Fresnel, A., Théorie de la lumière. Cinquième section: questions diverses d’optique [Theory of light. Fifth section: various questions of optics]. Letter to F. Arago, September 1818. In Oeuvers completes d’Augustin Fresnel Vol. 2, 627–636 (Imprimerie Impériale, 1868).

  2. 2.

    Fizeau, H. Sur les hypotheses relativesal, áther lumineux, et sur une expérience qui parait démontrer que le mouvement des corps change la vitesse avec laquelle la lumiere se propage dans leur intérieur [On the relative hypotheses, luminous ether, and on an experiment which seems to demonstrate that the movement of bodies changes the speed with which light propagates in their interior.]. CR Hebd. Acad. Sci. 33, 349–355 (1851).

    Google Scholar 

  3. 3.

    Landau, L. & Lifshitz, E. The Classical Theory of Fields Vol. 2 4th edn (Pergamon, 1975).

  4. 4.

    Lorentz, H. A. The Theory of Electrons and Its Applications to the Phenomena of Light and Radiant Heat (Teubner, 1916).

  5. 5.

    Unz, H. Relativistic magneto-ionic theory for drifting plasma in longitudinal direction. Phys. Rev. 146, 92–95 (1966).

    ADS  CAS  Google Scholar 

  6. 6.

    Moss, T., Burrell, G. & Hetherington, A. Measurement of Fresnel drag produced by electron motion in semiconductors. Proc. Roy. Soc. London Ser. A 308, 125–132 (1968).

    ADS  Google Scholar 

  7. 7.

    Almazov, L., Vas’ko, F. & Dykman, I. Influence of carrier drift on the propagation of electromagnetic wave in a solid-state plasma. JETP Lett. 16, 241–216 (1972).

    Google Scholar 

  8. 8.

    Basov, D. N., Fogler, M. M. & García de Abajo, F. J. Polaritons in van der Waals materials. Science 354, aag1992 (2016).

    PubMed  Google Scholar 

  9. 9.

    Low, T. et al. Polaritons in layered two-dimensional materials. Nat. Mater. 16, 182–194 (2017).

    ADS  CAS  PubMed  Google Scholar 

  10. 10.

    Caloz, C. et al. Electromagnetic nonreciprocity. Phys. Rev. Appl. 10, 047001 (2018).

    ADS  CAS  Google Scholar 

  11. 11.

    Jin, D. et al. Infrared topological plasmons in graphene. Phys. Rev. Lett. 118, 245301 (2017).

    ADS  PubMed  Google Scholar 

  12. 12.

    Tokura, Y. & Nagaosa, N. Nonreciprocal responses from non-centrosymmetric quantum materials. Nat. Commun. 9, 3740 (2018).

    ADS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Kumar, A. et al. Chiral plasmon in gapped Dirac systems. Phys. Rev. B 93, 041413 (2016).

    ADS  Google Scholar 

  14. 14.

    Song, J. C. W. & Rudner, M. S. Chiral plasmons without magnetic field. Proc. Natl Acad. Sci. USA 113, 4658–4663 (2016).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Jablan, M., Buljan, H. & Soljačić, M. Plasmonics in graphene at infrared frequencies. Phys. Rev. B 80, 245435 (2009).

    ADS  Google Scholar 

  16. 16.

    Woessner, A. et al. Highly confined low-loss plasmons in graphene–boron nitride heterostructures. Nat. Mater. 14, 421 (2015).

    ADS  CAS  PubMed  Google Scholar 

  17. 17.

    Dai, S. et al. Graphene on hexagonal boron nitride as a tunable hyperbolic metamaterial. Nat. Nanotechnol. 10, 682–686 (2015).

    ADS  CAS  PubMed  Google Scholar 

  18. 18.

    Ni, G. X. et al. Fundamental limits to graphene plasmonics. Nature 557, 530–533 (2018).

    ADS  CAS  PubMed  Google Scholar 

  19. 19.

    Son, S.-K. et al. Graphene hot-electron light bulb: incandescence from hBN-encapsulated graphene in air. 2D Mater. 5, 011006 (2017).

    Google Scholar 

  20. 20.

    Borgnia, D. S., Phan, T. V. & Levitov, L. S. Quasi-relativistic Doppler effect and non-reciprocal plasmons in graphene. Preprint at: (2015).

  21. 21.

    Sabbaghi, M., Lee, H.-W., Stauber, T. & Kim, K. S. Drift-induced modifications to the dynamical polarization of graphene. Phys. Rev. B 92, 195429 (2015).

    ADS  Google Scholar 

  22. 22.

    Duppen, B. V., Tomadin, A., Grigorenko, A. N. & Polini, M. Current-induced birefringent absorption and non-reciprocal plasmons in graphene. 2D Mater. 3, 015011 (2016).

    Google Scholar 

  23. 23.

    Bliokh, K., Rodríguez-Fortuño, F. J., Bekshaev, A., Kivshar, Y. & Nori, F. Electric-current-induced unidirectional propagation of surface plasmon-polaritons. Opt. Lett. 43, 963–966 (2018).

    ADS  CAS  PubMed  Google Scholar 

  24. 24.

    Sabbaghi, M., Lee, H.-W. & Stauber, T. Electro-optics of current-carrying graphene. Phys. Rev. B 98, 075424 (2018).

    ADS  CAS  Google Scholar 

  25. 25.

    W. Zhao et. al, Efficient Fizeau drag from Dirac electrons in monolayer graphene. Nature (2021).

  26. 26.

    Yu, Y.-J. et al. Tuning the graphene work function by electric field effect. Nano Lett. 9, 3430–3434 (2009).

    ADS  CAS  PubMed  Google Scholar 

  27. 27.

    Svintsov, D., Vyurkov, V., Ryzhii, V. & Otsuji, T. Hydrodynamic electron transport and nonlinear waves in graphene. Phys. Rev. B 88, 245444 (2013).

    ADS  Google Scholar 

  28. 28.

    Lucas, A. & Fong, K. C. Hydrodynamics of electrons in graphene. J. Phys. Condens. Matter 30, 053001 (2018).

    ADS  PubMed  Google Scholar 

  29. 29.

    Sun, Z., Basov, D. N. & Fogler, M. M. Universal linear and nonlinear electrodynamics of a Dirac fluid. Proc. Natl Acad. Sci. USA 115, 3285–3289 (2018).

    ADS  MathSciNet  CAS  PubMed  PubMed Central  MATH  Google Scholar 

  30. 30.

    Gao, H., Dong, Z. & Levitov, L. Plasmonic drag in a flowing Fermi liquid. Preprint at (2020).

  31. 31.

    Brower, D. L. et al. Fizeau interferometer for measurement of plasma electron current. Rev. Sci. Instrum. 75, 3399–3401 (2004).

    ADS  CAS  Google Scholar 

  32. 32.

    Tyson, R. E. et al. Far-infrared studies of the plasmon resonance of a drifting 2DEG. Superlattices Microstruct. 12, 371–374 (1992).

    ADS  CAS  Google Scholar 

  33. 33.

    Gramila, T. J., Eisenstein, J. P., MacDonald, A. H., Pfeiffer, L. N. & West, K. W. Mutual friction between parallel two-dimensional electron systems. Phys. Rev. Lett. 66, 1216–1219 (1991).

    ADS  CAS  PubMed  Google Scholar 

  34. 34.

    Gorbachev, R. V. et al. Strong Coulomb drag and broken symmetry in double-layer graphene. Nat. Phys. 8, 896–901 (2012).

    CAS  Google Scholar 

  35. 35.

    Gurevich, Y. G. & Mashkevich, O. L. The electron–phonon drag and transport phenomena in semiconductors. Phys. Rep. 181, 327–394 (1989).

    ADS  CAS  Google Scholar 

  36. 36.

    Narozhny, B. N. & Levchenko, A. Coulomb drag. Rev. Mod. Phys. 88, 025003 (2016).

    ADS  Google Scholar 

  37. 37.

    Lundeberg, M. B. et al. Tuning quantum nonlocal effects in graphene plasmonics. Science 357, 187–191 (2017).

    ADS  CAS  PubMed  Google Scholar 

  38. 38.

    Morgado, T. A. & Silveirinha, M. G. Nonlocal effects and enhanced nonreciprocity in current-driven graphene systems. Phys. Rev. B 102, 075102 (2020).

    ADS  CAS  Google Scholar 

  39. 39.

    Papaj, M. & Lewandowski, C. Plasmonic nonreciprocity driven by band hybridization in moiré materials. Phys. Rev. Lett. 125, 066801 (2020).

    ADS  CAS  PubMed  Google Scholar 

  40. 40.

    Bistritzer, R. & MacDonald, A. H. Moiré bands in twisted double-layer graphene. Proc. Natl Acad. Sci. USA 108, 12233–12237 (2011).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80–84 (2018).

    ADS  CAS  Google Scholar 

  42. 42.

    Dyakonov, M. & Shur, M. Shallow water analogy for a ballistic field effect transistor: new mechanism of plasma wave generation by dc current. Phys. Rev. Lett. 71, 2465–2468 (1993).

    ADS  CAS  PubMed  Google Scholar 

  43. 43.

    Gupta, R. & Ridley, B. Two-stream instability in two-dimensional degenerate systems. Phys. Rev. B 39, 6208 (1989).

    ADS  CAS  Google Scholar 

  44. 44.

    Aizin, G. R., Mikalopas, J. & Shur, M. Current-driven plasmonic boom instability in three-dimensional gated periodic ballistic nanostructures. Phys. Rev. B 93, 195315 (2016).

    ADS  Google Scholar 

  45. 45.

    Zolotovskii, I. O. et al. Plasmon-polariton distributed-feedback laser pumped by a fast drift current in graphene. Phys. Rev. A 97, 053828 (2018).

    ADS  CAS  Google Scholar 

  46. 46.

    Smetanin, I. V., Bouhelier, A. & Uskov, A. V. Coherent surface plasmon amplification through the dissipative instability of 2D direct current. Nanophotonics 8, 135–143 (2018).

    Google Scholar 

  47. 47.

    Zhang, Y., Small, J. P., Pontius, W. V. & Kim, P. Fabrication and electric-field-dependent transport measurements of mesoscopic graphite devices. Appl. Phys. Lett. 86, 073104 (2005).

    ADS  Google Scholar 

  48. 48.

    Liu, S. et al. Single crystal growth of millimeter-sized monoisotopic hexagonal boron nitride. Chem. Mater. 30, 6222–6225 (2018).

    CAS  Google Scholar 

  49. 49.

    Dean, C. R. et al. Boron nitride substrates for high-quality graphene electronics. Nat. Nanotechnol. 5, 722–726 (2010).

    ADS  CAS  PubMed  Google Scholar 

  50. 50.

    Purdie, D. G. et al. Cleaning interfaces in layered materials heterostructures. Nat. Commun. 9, 5387 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Wang, L. et al. One-dimensional electrical contact to a two-dimensional material. Science 342, 614–617 (2013).

    ADS  CAS  PubMed  Google Scholar 

  52. 52.

    Ben Shalom, M. et al. Quantum oscillations of the critical current and high-field superconducting proximity in ballistic graphene. Nat. Phys. 12, 318–322 (2016).

    Google Scholar 

  53. 53.

    Ocelic, N., Huber, A. & Hillenbrand, R. Pseudoheterodyne detection for background-free near-field spectroscopy. Appl. Phys. Lett. 89, 101124 (2006).

    ADS  Google Scholar 

  54. 54.

    Krishna Kumar, R. et al. Superballistic flow of viscous electron fluid through graphene constrictions. Nat. Phys. 13, 1182 (2017).

    CAS  Google Scholar 

  55. 55.

    Meeker, W. Q. & Escobar, L. A. Teaching about approximate confidence regions based on maximum likelihood estimation. Am. Stat. 49, 48–53 (1995).

    Google Scholar 

  56. 56.

    Jiang, B.-Y., Zhang, L. M., Castro Neto, A. H., Basov, D. N. & Fogler, M. M. Generalized spectral method for near-field optical microscopy. J. Appl. Phys. 119, 054305 (2016).

    ADS  Google Scholar 

Download references


Research on the physics and imaging of the plasmonic Fizeau effect at Columbia was supported by the US Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), under award no. DE-SC0018426. D.A.B. acknowledges the support from MIT Pappalardo Fellowship. M.M.F. is supported by the Office of Naval Research under grant ONR-N000014-18-1-2722. Work in the P.J.-H. group was supported by AFOSR grant FA9550-16-1-0382, the Gordon and Betty Moore Foundation’s EPiQS Initiative through grant GBMF9643, and Fundacion Ramon Areces. The development of new nanofabrication and characterization techniques enabling this work has been supported by the US DOE Office of Science, BES, under award DE-SC0019300. The development of the universal cryogenic platform used for scanning probe measurements is supported as part of the Energy Frontier Research Center on Programmable Quantum Materials funded by the US Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), under award no. DE-SC0019443. The development of infrared nano-optics is supported by Vannevar Bush Faculty Fellowship to D.N.B., ONR-VB: N00014-19-1-2630. D.N.B. is Moore Investigator in Quantum Materials EPIQS #9455. This work also made use of the Materials Research Science and Engineering Center Shared Experimental Facilities supported by the National Science Foundation (NSF) (grant no. DMR-0819762). Support from the Materials Engineering and Processing programme of the National Science Foundation, award number CMMI 1538127 for hBN crystal growth is also greatly appreciated.

Author information




D.A.B. and D.N.B. conceived and supervised the project. Y.D., L.X., S.Z., D.A.B. and D.N.B. designed the experiments. I.Y.P. and D.A.B. fabricated the devices. S.L. and J.H.E. provided the isotopic hBN crystals. Y.D., L.X., A.S.M. and R.P. performed the experimental measurements. Y.D., L.X., R.J., F.L.R., D.A.B. and D.N.B. analysed the experimental data. Z.S., M.M.F., A.J.M. and L.S.L. developed the theoretical analysis of the experimental data with input from P.J.-H., H.G. and Z.D. Y.D., L.X., Z.S., M.M.F., D.A.B. and D.N.B. co-wrote the manuscript with input from all co-authors.

Corresponding authors

Correspondence to D. A. Bandurin or D. N. Basov.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature thanks Joel Cox, Jiahua Duan and Hugen Yan for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 Gate-dependent transport of a typical device.

Two-terminal resistance R2pt of a typical device as a function of the back-gate voltage Vg at T = 170 K. Inset shows the transport measurement configuration, where a source meter (Keithley 2450) was used to source gate voltage, and a lock-in amplifier (SR 830) was used to measure the resistance of the entire device.

Extended Data Fig. 2 Experimental configuration for measuring current-gating effects.

Voltage is applied across the source/drain electrode and SPP imaging is performed close to the drain. The black streamlines symbolize the d.c. current. The voltmeter measures the electrostatic potential of the SPP launcher as a function of the biasing current.

Extended Data Fig. 3 Current–voltage characteristics of a typical device.

The black symbols represent the source–drain voltage while sourcing current through the device. The blue symbols represent the simultaneously measured voltage on the SPP launcher, representing the current-gating effect induced by the biasing current. Inset shows a magnified view of the low-voltage region.

Extended Data Fig. 4 Characteristic real-space nano-infrared and topography images.

a, Near-field image taken in the vicinity of a gold launcher at T = 170 K and Vg = 50 V. Gold-launched λp fringes and tip-launched λp/2 fringes near the graphene edge are clearly visible. Dashed rectangles mark the regions magnified in c and d. b, AFM topography image taken simultaneously with the near-field image in a. The graphene region is uniform with minimal topographic variations. c, d, Magnified near-field images near graphene edges, showing the tip-launched and edge-reflected SPP fringes with λp/2 periodicity.

Extended Data Fig. 5 Real-space nano-infrared and topography line profiles under current.

ad, Representative results for simultaneously taken topography (a, b) and near-field (c, d) data as a function of current density at T = 170 K and Vg = −50 V. a, AFM topography collected in the vicinity of a gold launcher on the left of the field of view. The 2D plot is assembled from AFM line profiles measured at different current densities while scanning along the same line in real space. Red and black arrows and dashed lines indicate positions where the averaged line profiles in b are acquired. b, Averaged line profiles of AFM topography for current densities of ±0.75 mA μm−1, whose topography signals are essentially the same. One of the line profiles is shifted vertically for clarity. c, Near-field data, taken simultaneously with a. A standard one-dimensional Fourier filter was applied here to reduce noise. d, Averaged line profiles of the near-field signal in c for the same current densities as the topography data in b. A Fizeau shift is clearly visible.

Extended Data Fig. 6 Uncertainty analysis for fitted SPP wavelength.

a, Distribution of the least-squares estimate of the SPP wavelength in equation (7) generated by Monte Carlo simulation. b, Examples of typical simulated SPP line profiles used for analysis in a. c, I: Dependence of variance (bright and dark red lines) and bias (blue line) of wavelength estimate on spatial resolution (SR). Bright red line corresponds to zero signal noise (σy/A = 0). Dark red and blue lines correspond to σy/A = 10%; II: Dependence of bias and variance in wavelength estimate on pixel size in units of wavelength. The pixel size has minimal effect on σλ as long as one samples above the Nyquist rate, as indicated by the vertical green dashed line. d, I: Strong dependence of error in wavelength estimate on signal noise σy. The variance of the wavelength estimate (bright red, dark red and green lines) will increase roughly linearly with σy until about 25%, and the bias (bright and dark blue lines) is less than 1 nm for SR = 20 nm; II: Dependence of error in wavelength estimate on SPP propagation length 1/q2. Both the variance (bright and dark red lines) and the bias (blue line) of the wavelength estimate improve with 1/q2, even more so when there is positioning noise σx (dark red line). e, Assessing the statistical significance of the Fizeau shift using an F test. Solid red and purple lines represent the dependence of F statistics on the sample wavelength-shift standard deviation sΔλ. Purple line assumes α = 0 and red line assumes α is finite (see text). Cyan shaded region corresponds to F statistics that reject the null hypothesis of no Fizeau shift (F > Fcrit = 2.65). Vertical dashed line corresponds to wavelength-shift standard deviation sΔλ estimated from data in Fig. 3a, b.

Extended Data Fig. 7 Additional datasets revealing plasmonic Fizeau drag in graphene.

Other representative data when scanning along the same line at different d.c. currents (averaged ±25 μA μm−1 for each profile) for different gate voltages, temperatures and devices. Within a set of polariton line profiles, the first polariton fringes are aligned to enable better visual inspection of Fizeau shifts. Line profiles are shifted vertically for clarity. Within each panel, the fitted line profiles of the smallest and largest current densities are shown in the lower panel for visual comparison. The images of the devices are near-field scattering amplitude measured at 170 K without gating. a, Device 1, T = 170 K, Vg = −47 V; b, Device 1, T = 170 K, Vg = +47 V; c, Device 1, T = 170 K, Vg = −60 V; d, Device 1, T = 60 K, Vg = +60 V; e, Device 2, T = 170 K, Vg = +50 V; f, Device 3, T = 60 K, Vg = +60 V.

Supplementary information

Supplementary Information

This file contains a list of the notations, details regarding the theory of Fizeau plasmon drag in graphene (this section includes Supplementary Figures 1-4), Supplementary Table 1 and the Supplementary References.

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Dong, Y., Xiong, L., Phinney, I.Y. et al. Fizeau drag in graphene plasmonics. Nature 594, 513–516 (2021).

Download citation


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links