Abstract
Electron-electron interactions play a critical role in many condensed matter phenomena, and it is tempting to find a way to control them by changing the interactions’ strength. One possible approach is to place a studied system in proximity of a metal, which induces additional screening and hence suppresses electron interactions. Here, using devices with atomically-thin gate dielectrics and atomically-flat metallic gates, we measure the electron-electron scattering length in graphene and report qualitative deviations from the standard behavior. The changes induced by screening become important only at gate dielectric thicknesses of a few nm, much smaller than a typical separation between electrons. Our theoretical analysis agrees well with the scattering rates extracted from measurements of electron viscosity in monolayer graphene and of umklapp electron-electron scattering in graphene superlattices. The results provide a guidance for future attempts to achieve proximity screening of many-body phenomena in two-dimensional systems.
Similar content being viewed by others
Introduction
Electrostatic screening by conducting gates has previously been employed to suppress charge inhomogeneity in graphene1,2,3, alter its plasmon spectra4,5, and renormalize electronic spectra of monolayer semiconductors6,7. Elementary electrostatics tells us that the electron charge \(e\) placed at the distance \(d\) from a bulk conductor leads to a dipole potential evolving as \(2ed^2/r^3\) at large in-plane distances \(r \gg d\), which is much weaker than the original, unscreened Coulomb potential, \(e/r\). Accordingly, a metallic gate placed sufficiently close to another electronic system can alter its electron–electron (e–e) interactions. They can be parametrized by the e–e scattering length \(\ell _{{\mathrm{ee}}}\). From the above considerations, one can infer that what matters most is the ratio \(d/D\), where \(D\approx1/\sqrt n\) is the average separation between electrons and \(n\) is the carrier concentration. For a two-dimensional (2D) electron system with typical \(n = 10^{12}\,{\mathrm{cm}}^{ - 2}\), \(D \approx 10\,{\mathrm{nm}}\) and, therefore, the inferred gate separation \(d \approx D\) is relatively easy to achieve experimentally. However, as shown below, the naïve expectations fail because of a small numerical factor \(\delta\) such that e–e interactions for massless Dirac fermions are altered only if \(d \le \delta D \approx0.03\,\varepsilon D\), where \(\varepsilon\) is the gate dielectric’s permittivity. For typical gate dielectrics with \(\varepsilon \,{<}\, 5\), the required separation falls into a \(1\,{\mathrm{nm}}\) range. For massive charge carriers such as those in bilayer graphene and 2D semiconductors, even smaller (atomic scale) \(d\) are necessary for efficient screening (Methods). It seems impossible to realize the required small \(d\) because of, for example, inevitable surface roughness of the conducting and insulating films used for gating and electrical leakage through dielectrics of nanometer thickness.
In this communication, we achieve the extremely challenging conditions for proximity-gate screening by using van der Waals heterostructures with atomically thin dielectric layers and atomically flat gates. Measurements of viscous electron flow in graphene and umklapp scattering rates in graphene superlattices provide two independent but complementary ways to quantify the effect of proximity screening and its dependence on \(d\), \(n\), and temperature (\(T\)).
Results
Experimental devices and measurement setup
Our devices were graphene monolayers encapsulated between hexagonal boron nitride (hBN) crystals whereas graphite monocrystals served as a bottom gate (Fig. 1). These heterostructures were fabricated using the standard dry-transfer procedures1 described in Methods. Multiterminal Hall bar devices with several point contacts and closely placed voltage probes (Fig. 1a) were then defined by electron-beam lithography and plasma etching. An extra metal gate was deposited on top of the heterostructures, which allowed us to vary \(n\) without applying voltages to the bottom screening gate. This was particularly important for our case of ultra-thin dielectrics to avoid their accidental breakdown and electrical leakage. The minimum thickness \(d\) for the gate dielectric (Fig. 1b) was limited to four hBN layers (i.e. ∼\(1.3\,{\mathrm{nm}}\)) because thinner crystals exhibited notable electron tunneling8. The devices typically had low-\(T\) mobility \(\mu\) of about \(10^6\,{\mathrm{cm}}^2\,{\mathrm{V}}^{ - 1}\,{\mathrm{s}}^{ - 1}\) and highly reproducible characteristics such that, at finite \(T\), their longitudinal resistivity \(\rho\) was practically independent of \(d\) (Supplementary Note 1; Supplementary Fig. 1). This ensured that the reported behavior of \(\ell _{{\mathrm{ee}}}\) was due to changes in \(d\) rather than transport characteristics. Because graphite is a semimetal with a relatively low carrier density of ∼\(10^{19}\,{\mathrm{cm}}^{ - 3}\), we also crosschecked that our conclusions were independent of the gate material by using screening gates made from layered metals such as Bi2Sr2CaCu2O8+x and TaS2 (Methods; Supplementary Fig. 2).
To demonstrate that e–e interactions can be tuned by proximity-gate screening, a reliable diagnostic tool is essential. Many quantum transport characteristics are known to be affected by the strength of e–e interactions. For example, the phase breaking length depends on it and can be measured in quantum interference experiments9 (other possibilities are discussed in ref. 10). In principle, it should be possible to use such mesoscopic physics tools to probe e–e interactions in graphene but, because of its ballistic transport at micrometer-scale distances, the approach is not easy to implement in practice and its results could be difficult to interpret. On the other hand, recent experiments have shown that graphene at finite \(T\) and away from the charge neutrality point (NP) exhibits pronounced hydrodynamic effects11,12,13,14,15, which allowed measurements of the kinematic electron viscosity \(\nu _0\). The extracted values of \(\ell_{{\mathrm{ee}}}\), obtained by inverting the relation \(\nu _0 = v_{\mathrm{F}}\ell_{{\mathrm{ee}}}/4\) derived in refs. 16,17,18 (\(v_{\mathrm{F}} \approx 10^6\,{\mathrm{ms}}^{ - 1}\) is the Fermi velocity), were in quantitative agreement with theory (because viscosity stems from e–e interactions, the proportionality is not surprising and dimensionally correct). The viscosity measurements can be carried out using three complementary approaches: vicinity resistance12,13, point-contact geometry14,18, and the viscous Hall effect15. Below we use all three to show that \(\ell _{{\mathrm{ee}}}\) changes with \(d\). In a complementary approach, we demonstrate that umklapp e–e scattering in graphene superlattices19 is also affected by proximity-gate screening.
Enhanced electron viscosity
First, we demonstrate the screening effect qualitatively. Figure 1c shows that the vicinity resistance \(R_{\mathrm{V}}\) is notably affected if a thin gate dielectric is employed. Vicinity measurements are discussed in detail in ref. 12 but, briefly, an electric current is injected through a narrow contact into a wide graphene channel. The negative voltage drop arising locally from a viscous electron flow is detected using a vicinity contact at a short distance \(L\) from the current-injecting contact (Fig. 1a). One can see from Fig. 1c that, as \(T\) increases, \(R_{\mathrm{V}}\) first decreases and then becomes negative. This indicates a transition from the ballistic transport regime (positive \(R_{\mathrm{V}}\)) into a regime where ballistics is strongly affected by e–e scattering13. The minimum in \(R_{\mathrm{V}}(T)\) corresponds to the condition \(\ell _{{\mathrm{ee}}} \approx L\) and indicates an onset of hydrodynamic behavior13. As \(\ell _{{\mathrm{ee}}}\) decreases further with increasing \(T\), \(R_{\mathrm{V}}\) becomes less negative and eventually positive, being dominated by currents caused by electron–phonon scattering12,13. The dependences \(R_{\mathrm{V}}(T)\) shown in Fig. 1c were measured for two similar devices at the same \(L\). One had \(d \approx 300\,{\mathrm{nm}}\) (conventional Si back gate) whereas the other was made using four-layer hBN as the gate dielectric. Despite the similar behavior of \(R_{\mathrm{V}}(T)\), the curve for \(d \approx 1.3\,{\mathrm{nm}}\) is clearly shifted to higher \(T\). The shift direction indicates that the nearby gate causes an increase in \(\ell _{{\mathrm{ee}}}\), which is equivalent to a reduction in electron temperature by ∼\(30\,{\mathrm{K}}\). Note that for \(T\) above \(100\,{\mathrm{K}}\) where the hydrodynamic regime develops, electron transport in high-quality graphene is universal and insensitive to experimental details.
Similar phenomenology was observed in the point-contact geometry (Fig. 1a). Again, the \(T\)-dependence of the point-contact resistance \(R_{{\mathrm{PC}}}\) exhibits a clear minimum due to a viscous flow14. The shift to higher \(T\) for the device with a proximity gate (Fig. 1d) indicates an increase in \(\ell _{{\mathrm{ee}}}\) for a given \(T\) (also, see Supplementary Note 3). Such influence of the proximity gating was consistently observed in all our experiments. The \(R_{\mathrm{V}}\) and \(R_{{\mathrm{PC}}}\) dependences could also be used to extract \(\ell _{{\mathrm{ee}}}(T)\) following the recipe reported in refs. 12,14. Unfortunately, we found that, for atomically thin gate dielectrics, detailed behavior of \(R_{\mathrm{V}}\) and, to some extent, \(R_{{\mathrm{PC}}}\) notably varied between different devices with nominally the same \(d\). Those variations can be traced back to the fact that \(R_{\mathrm{V}}\) is sensitive to current injector’s geometry12 whereas a viscous contribution to \(R_{{\mathrm{PC}}}\) becomes smaller for close-gate devices as compared to those with thicker gate dielectrics.
In contrast to the vicinity and point-contact measurements, the viscous Hall effect15 was found to be very robust, yielding quantitatively same results for different devices with same \(d\). Accordingly, for quantitative analysis of how \(\ell _{{\mathrm{ee}}}\) depended on \(d\), we focused on the latter measurements. The Hall viscosity experiments utilize the already discussed vicinity geometry (Fig. 1a) but a non-quantizing magnetic field \(B\) is applied perpendicular to graphene15. The field leads to an asymmetry in the potential created by the viscous flow around the injection contact (insets of Fig. 1e, f). The viscous contribution asymmetric in \(B\) is called the viscous Hall resistance \(R_{\mathrm{A}}\) and given by15,17
where \(\xi \left( x \right)\) is a dimensionless function17, \(\tau\) is the transport scattering time, \(B_0 = E_{\mathrm{F}}/(8|e|\nu _0)\) is a characteristic magnetic field, and \(E_{\mathrm{F}}\) is the Fermi energy. Because \(|\xi \left( x \right)|\) is a monotonically decreasing function of its argument for \(x \,{> }\, 0\), \(|R_{\mathrm{A}}|\) increases with increasing \(\ell _{{\mathrm{ee}}}\) and, accordingly, devices with weaker e–e scattering should exhibit larger \(|R_{\mathrm{A}}|\).
To illustrate the effect of proximity-gate screening on Hall viscosity, Fig. 1e, f plot \(R_{\mathrm{A}}(B)\) for two representative devices with \(d \approx 1.7\) and \(300\,{\mathrm{nm}}\). The curves are taken under exactly the same conditions for several same\(n\). As the two devices exhibited close \(\rho\) and \(\tau\) (Supplementary Fig. 1), the profound difference between Fig. 1e and f can only be attributed to different screening. The device with the thin dielectric exhibited much larger Hall viscosity than the reference device, and the effect was most pronounced at low \(n\). This behavior proves again that the proximity screening suppresses e–e scattering, in agreement with the conclusions reached from the vicinity and point-contact measurements.
For the known transport characteristics (\(\rho\) and \(\tau\)), Eq. (1) allows us to convert \(R_{\mathrm{A}}\) into \(\ell _{{\mathrm{ee}}}\), as described in detail in ref. 15. Figure 2a shows examples of \(\ell _{{\mathrm{ee}}}(T)\) found for close-gate and reference devices. At all \(T\), the screened device displays \(\ell _{{\mathrm{ee}}}\) approximately twice longer than that in the standard device of the same electronic quality. This agrees well with many-body theory (solid curves in Fig. 2; Supplementary Fig. 3). Importantly, the proximity-gate screening qualitatively changes the dependence \(\ell _{{\mathrm{ee}}}(n)\) so that, away from the NP, \(\ell _{{\mathrm{ee}}}\) decreases with increasing \(n\) (Fig. 2b). This contrasts with monotonically increasing \(\ell _{{\mathrm{ee}}}(n)\) for the reference devices, which was also reported previously14,15. Figure 2c summarizes our results by showing \(\ell _{{\mathrm{ee}}}\) measured for more than 10 different devices at characteristic \(n\) and \(T\) where viscous effects become most pronounced in graphene. Despite the experimental scatter, Fig. 2c clearly shows that \(\ell _{{\mathrm{ee}}}\) can be altered appreciably by using thin gate dielectrics, if \(d\) is smaller than a few nm.
To explain the observed dependences of \(\ell _{{\mathrm{ee}}}\) on \(n\) and \(d\), we carried out numerical calculations in the random phase approximation for the dynamically screened interactions10,21,22. The gate was modeled as a perfect conductor, and small departures from this model caused by a finite carrier density were estimated in Supplementary Note 4. The results are shown by the solid curves in Fig. 2. No fitting parameters were used, except for multiplying all the theoretical curves by the same small factor of \(1.3\) (its non-Fermi-liquid origins are discussed in Supplementary Note 4). However, to gain better insight about the observed behavior, we also derived the following analytical expression:
where \(k_{\mathrm{F}} = \sqrt {\pi n}\) and \(q_{{\mathrm{TF}}} = 4\alpha _{{\mathrm{ee}}}k_{\mathrm{F}}\) are the Fermi and Thomas–Fermi wavenumbers, respectively. Here, \(\alpha _{{\mathrm{ee}}} \equiv e^2/({\it{\epsilon }}\hbar v_{\mathrm{F}}) \approx 2.2/\varepsilon\) is graphene’s coupling constant and \(k_{\mathrm{B}}\) is the Boltzmann constant (Supplementary Note 4 discusses the case of generally anisotropic \(\varepsilon\)). The expression is accurate in the Fermi-liquid regime (\(k_{\mathrm{B}}T \ll E_{\mathrm{F}}\)), where it matches our numerical results (Supplementary Note 4). The last term in Eq. (2) appears due to the gate presence, and the key parameter describing its screening effect is \(dq_{{\mathrm{TF}}}\). In the far-gate regime, \(d \gg 1/q_{{\mathrm{TF}}}\), Eq. (2) reduces to the standard unscreened expression22. In the opposite limit, \(d \ll 1/q_{{\mathrm{TF}}}\), e–e scattering is strongly reduced due to screening, and \(\ell _{{\mathrm{ee}}}\) increases with decreasing both \(d\) and \(n\), as \(1/d^2\) and approximately \(1/\sqrt n\), respectively, in agreement with our experiment (Fig. 2). The latter dependence is opposite to the unscreened case, where \(\ell _{{\mathrm{ee}}}\) increases as \(\sqrt n\), in agreement with the results of Fig. 2b. The crossover between the far- and close-gate regimes occurs at a critical distance \(d_{\mathrm{c}}\) such that \(d_{\mathrm{c}} \approx 1/2q_{{\mathrm{TF}}} = 1/(8\alpha _{{\mathrm{ee}}}k_{\mathrm{F}})\), which translates into the previously introduced parameter \(\delta \approx0.03\varepsilon\). For hBN with \(\varepsilon \approx 3.5\) and at typical \(n = 10^{12}\,{\mathrm{cm}}^{ - 2}\), we obtain \(d_{\mathrm{c}} \approx 1.1\,{\mathrm{nm}}\), which explains why the gate screening becomes noticeable only for our smallest \(d\) (Fig. 2c). Further information about our theoretical analysis is provided in Supplementary Note 4.
Suppression of e–e umklapp scattering
To check how robust our conclusions are, we have also examined the effect of gate-induced screening on umklapp e–e scattering19 that dominates resistivity \(\rho\) of graphene-on-hBN superlattices at elevated \(T\). We made several superlattice devices with the moiré periodicity \(\lambda \approx 15\,{\mathrm{nm}}\), as confirmed by the periodicity of Brown–Zak oscillations23 and the appearance of secondary NPs24,25,26,27 at the expected \(n\) (Fig. 3a). One of the devices was the standard Hall bar with \(d = 300\,{\mathrm{nm}}\), like those reported previously19. The other two were same in design but had a bottom graphite gate placed at short \(d\), as in the above viscosity experiments. Figure 3 shows typical \(\rho (n,T)\) measured for these graphene superlattices. For \(d = 300\,{\mathrm{nm}}\), the observed behavior was same as reported previously, and the \(T\)-dependent part (\(\Delta \rho\)) of graphene superlattice’s resistivity could be described quantitatively by umklapp e–e scattering19. It is responsible for the rapid increase of \(\Delta \rho \propto T^2\) (Fig. 3b). The proximity-gate screening notably suppressed \(\Delta \rho (T)\), by a factor \(> 2\) for \(d \approx 1.3\,{\mathrm{nm}}\). Our theoretical analysis (Supplementary Note 5) shows that \({\mathrm{\Delta }}\rho\) for the close-gate devices should exhibit the same \(T\) dependence (\({\propto} T^2\)) but with a reduced absolute value. The umklapp e–e scattering length, \(\ell _{{\mathrm{ee}}}^{\mathrm{U}}\), is governed by distinctive processes with a momentum transfer of \({\sim} \hbar g\) where \(g = \frac{{4\pi }}{{\sqrt 3 }}\lambda ^{ - 1}\) is the superlattice reciprocal vector. As shown in Supplementary Note 5, proximity screening for \(\ell _{{\mathrm{ee}}}^{\mathrm{U}}\) becomes important if \(d \,{<}\, 0.1\lambda\), which means that few-nm-thick gate dielectrics are essential to observe the screening effect, similarly to the case of hydrodynamic transport. It is convenient to quantify the changes in umklapp scattering by the dimensionless ratio, \(\Delta \rho \left( \infty \right)/\Delta \rho \left( d \right)\approx\ell _{{\mathrm{ee}}}^{\mathrm{U}}(d)/\ell _{{\mathrm{ee}}}^{\mathrm{U}}(\infty )\). The results are plotted in the inset of Fig. 3b and show good agreement with theory (for details, see Supplementary Note 5 and Supplementary Fig. 4).
Discussion
Electron–electron scattering in monolayer graphene at finite \(n\) can strongly be suppressed if a metallic gate is placed at \(d\) of about \(1\,{\mathrm{nm}}\). This “close-gate” regime has become accessible due to the use of van der Waals assembly that allows atomically sharp interfaces and ultra-thin dielectrics. It is tempting to exploit the outlined strategy to assess interaction phenomena near the NP where low \(n\) allow the condition \(d \ll 1/\sqrt n\) to be satisfied easier but interpretation of some previous observations had proven difficult. Other interesting candidates for proximity screening are exotic phenomena driven by strong correlations (e.g., various many-body phases in twisted bilayer graphene28,29,30) and, especially, interaction effects governed by lengths longer than \(\ell _{{\mathrm{ee}}}\). The experimental challenge to reach the close-gate regime can partially be mitigated by using high-\(\varepsilon\) dielectrics.
Methods
Device fabrication
Our heterostructures were assembled using stamps made from polypropylene carbonate (PPC) as a sacrificial polymer placed on polydimethylsiloxane (PDMS). Such polymer stamps were used to pick up exfoliated thin crystals in the following sequence: top hBN (typically thicker than \(30\,{\mathrm{nm}}\)), monolayer graphene, and thin bottom hBN. The latter served as a gate dielectric in the final device configuration (Fig. 1b), and its thickness was determined by atomic force microscopy. The resulting hBN/graphene/hBN stack was then released onto relatively small graphite crystals with thickness of 3–10 nm, which were prepared in advance on an oxidized Si wafer. The stack was large enough to extend outside the bottom graphite region, which allowed us to make quasi-one-dimensional contacts to graphene31 without electrically contacting the graphite gate. The metallic contacts were defined by electron-beam lithography. We first used a mixture of CHF3 and O2 to plasma-etch hBN/graphene and expose the required contact regions. This was followed by deposition of \(2\,{\mathrm{nm}}\) Cr/\(60\,{\mathrm{nm}}\) Au to make Ohmic contacts to graphene. A gold top gate was then fabricated using another round of electron-beam lithography and, also, served as an etching mask for the final etching step to define the Hall bar geometry.
The devices with other metallic gates (Bi2Sr2CaCu2O8+x and TaS2) required fabrication in an oxygen- and moisture-free atmosphere of a glovebox32 to avoid deterioration of the metal surfaces. Even using glovebox encapsulation, we observed a notable reduction in graphene’s quality for the above gate materials, presumably because of electrical charges at the exposed surfaces (for small \(d\), typical \(\mu\) became \(< 10^5\,{\mathrm{cm}}^2\,{\mathrm{V}}^{ - 1}{\mathrm{s}}^{ - 1}\) and charge inhomogeneity near the NP considerably increased). Accordingly, reliable measurements of \(\ell _{{\mathrm{ee}}}\) in this case were only possible at high \(n \,{\gtrsim}\, 2.0 \times 10^{12}\,{\mathrm{cm}}^{ - 2}\) (Supplementary Note 2). We also note that encapsulated graphene devices with the conventional gates made by metal deposition on top of a thin gate dielectric (\(d \,{<}\, 2\,{\mathrm{nm}}\)) exhibited extremely low \(\mu\) of only \({\sim} 10^4\,{\mathrm{cm}}^2\,{\mathrm{V}}^{ - 1}{\mathrm{s}}^{ - 1}.\) Such poor electronic quality made it impossible to carry out the \(\ell _{{\mathrm{ee}}}\) measurements described in the main text.
Electrical measurements
The devices were measured in a variable temperature insert that allowed stable \(T\) between \(2\) and \(300\,{\mathrm{K}}\). The standard lock-in amplifier techniques were employed using excitation currents of typically 0.1–1 μA at a frequency of \(30.5\,{\mathrm{Hz}}\). For measurements of Hall viscosity, we used the same vicinity geometry as shown in the schematic of Fig. 1a. The distance between injector and detector contacts was usually between \(0.5\) and \(1.5\,\upmu {\mathrm{m}}\). The viscous Hall resistance was determined as an antisymmetric-in-\(B\) component of the vicinity resistance in fields below \(\pm 30\,{\mathrm{mT}}\). For the point-contact measurements, we employed the quasi-four-probe geometry by driving the current through the wide contacts (on the left and right in Fig. 1a) and using the leads next to the studied constrictions as voltage probes.
Proximity screening for systems with the parabolic spectrum
The close-gate condition depends on the density of states at the Fermi energy of the material one wants to control. We have studied graphene not only because of its electronic quality but also because of the low-density of states provided by its Dirac spectrum. For a 2D system with the conventional parabolic spectrum, the close-gate condition is much more difficult to achieve. In the latter case, a conducting gate can provide efficient screening of e–e interactions only for distances \(d\) below \(\sim\varepsilon \,m_{\mathrm{e}}a_{\mathrm{B}}/(2\,N_{\mathrm{f}}\,m_{{\mathrm{eff}}})\), where \(a_{\mathrm{B}} \approx 0.5\) Å is the Bohr radius, \(m_{\mathrm{e}}\) and \(m_{{\mathrm{eff}}}\) are the free-electron and effective masses, respectively, and \(N_{\mathrm{f}}\) is the number of spin/valley flavors. Here, \(\varepsilon = \varepsilon _ \bot\) is the perpendicular component of the dielectric permittivity of a gate dielectric. For bilayer graphene33,34 with \(N_{\mathrm{f}} = 4\), \(m_{{\mathrm{eff}}} \ge 0.03\,m_{\mathrm{e}}\) and using hBN as a dielectric (\(\varepsilon _ \bot \approx 3.5\)), the close-gate condition requires d < 7 Å, which is essentially out of experimental reach.
Data availability
The data that support our findings are available upon reasonable request from M.K.
Code availability
Only commercial software was used in this study.
Change history
11 June 2020
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
References
Yankowitz, M. et al. van der Waals heterostructures combining graphene and hexagonal boron nitride. Nat. Rev. Phys. 1, 112–125 (2019).
Ponomarenko, L. A. et al. Tunable metal-insulator transition in double-layer graphene heterostructures. Nat. Phys. 7, 958–961 (2011).
Lu, C.-P. et al. Local, global, and nonlinear screening in twisted double-layer graphene. Proc. Natl. Acad. Sci. USA 113, 6623–6628 (2016).
Lundeberg, M. B. et al. Tuning quantum nonlocal effects in graphene plasmonics. Science 357, 187–191 (2017).
Alcaraz Iranzo, D. et al. Probing the ultimate plasmon confinement limits with a van der Waals heterostructure. Science 360, 291–295 (2018).
Qiu, Z. et al. Giant gate-tunable bandgap renormalization and excitonic effects in a 2D semiconductor. Sci. Adv. 5, eaaw2347 (2019).
Waldecker, L. et al. Rigid band shifts in two-dimensional semiconductors through external dielectric screening. Phys. Rev. Lett. 123, 206403 (2019).
Britnell, L. et al. Electron tunneling through ultrathin boron nitride crystalline barriers. Nano Lett. 12, 1707–1710 (2012).
Imry, Y. Introduction to Mesoscopic Physics (Oxford University Press, Oxford, 1997).
Giuliani, G. F. & Vignale, G. Quantum Theory of the Electron Liquid (Cambridge University Press, Cambridge, 2005).
Lucas, A. & Fong, K. C. Hydrodynamics of electrons in graphene. J. Phys. Condens. Matter 30, 053001 (2018).
Bandurin, D. A. et al. Negative local resistance caused by viscous electron backflow in graphene. Science 351, 1055–1058 (2016).
Bandurin, D. A. et al. Fluidity onset in graphene. Nat. Commun. 9, 4533 (2018).
Krishna Kumar, R. et al. Superballistic flow of viscous electron fluid through graphene constrictions. Nat. Phys. 13, 1182–1185 (2017).
Berdyugin, A. I. et al. Measuring Hall viscosity of graphene’s electron fluid. Science 364, 162–165 (2019).
Principi, A., Vignale, G., Carrega, M. & Polini, M. Bulk and shear viscosities of the two-dimensional electron liquid in a doped graphene sheet. Phys. Rev. B 93, 125410 (2016).
Pellegrino, F. M. D., Torre, I. & Polini, M. Nonlocal transport and the Hall viscosity of two-dimensional hydrodynamic electron liquids. Phys. Rev. B 96, 195401 (2017).
Guo, H., Ilseven, E., Falkovich, G. & Levitov, L. S. Higher-than-ballistic conduction of viscous electron flows. Proc. Natl Acad. Sci. USA 114, 3068–3073 (2017).
Wallbank, J. R. et al. Excess resistivity in graphene superlattices caused by umklapp electron–electron scattering. Nat. Phys. 15, 32–36 (2019).
Torre, I., Tomadin, A., Geim, A. K. & Polini, M. Nonlocal transport and the hydrodynamic shear viscosity in graphene. Phys. Rev. B 92, 165433 (2015).
Li, Q. & Das Sarma, S. Finite temperature inelastic mean free path and quasiparticle lifetime in graphene. Phys. Rev. B 87, 085406 (2013).
Polini, M. & Vignale, G. in No-nonsense Physicist: An Overview of Gabriele Giuliani’s Work and Life (eds. Polini, M. et al.) 107–124 (Scuola Normale Superiore, 2016).
Krishna Kumar, R. et al. High-temperature quantum oscillations caused by recurring Bloch states in graphene superlattices. Science 357, 181–184 (2017).
Yankowitz, M. et al. Emergence of superlattice Dirac points in graphene on hexagonal boron nitride. Nat. Phys. 8, 382–386 (2012).
Ponomarenko, L. A. et al. Cloning of Dirac fermions in graphene superlattices. Nature 497, 594 (2013).
Dean, C. R. et al. Hofstadter’s butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598–602 (2013).
Hunt, B. et al. Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure. Science 340, 1427–1430 (2013).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018).
Pizarro, J. M., Rösner, M., Thomale, R., Valentí, R. & Wehling, T. O. Internal screening and dielectric engineering in magic-angle twisted bilayer graphene. Phys. Rev. B 100, 161102 (2019).
Goodwin, Z. A. H. et al. Critical role of device geometry for the phase diagram of twisted bilayer graphene. Phys. Rev. B 101, 165110 (2020).
Wang, L. et al. One-dimensional electrical contact to a two-dimensional material. Science 342, 614–617 (2013).
Cao, Y. et al. Quality heterostructures from two-dimensional crystals unstable in air by their assembly in inert atmosphere. Nano Lett. 15, 4914–4921 (2015).
Hwang, E. H. & Das Sarma, S. Screening, Kohn anomaly, Friedel oscillation, and RKKY interaction in bilayer graphene. Phys. Rev. Lett. 101, 156802 (2008).
Borghi, G., Polini, M., Asgari, R. & MacDonald, A. H. Dynamical response functions and collective modes of bilayer graphene. Phys. Rev. B 80, 241402(R) (2009).
Acknowledgements
This work was supported by the European Research Council, Graphene Flagship, EPSRC Grand Challenges (EP/N010345/1), the Royal Society and Lloyd’s Register Foundation. A.I.B. and W.K. were supported by Graphene NowNANO Doctoral Training Programme.
Author information
Authors and Affiliations
Contributions
M.K. and S.G.X. carried out the project and analyzed the experimental data with the help of A.I.B., A.P., M.P., and A.K.G. A.K.G. suggested and supervised the project. S.G.X., N.X., P.K., M.H., and R.V.G. fabricated devices. M.K., A.I.B., W.K., and R.K.K. performed electrical measurements. A.P., S.S., V.I.F., and M.P. provided theoretical support. K.W. and T.T. supplied hBN crystals. M.K., S.G.X., A.I.B., A.P., S.S., I.V.G., V.I.F., M.P., and A.K.G. wrote the manuscript. All the authors contributed to discussions.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Kin Chung Fong and the other, anonymous, reviewer(s) 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.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kim, M., Xu, S.G., Berdyugin, A.I. et al. Control of electron-electron interaction in graphene by proximity screening. Nat Commun 11, 2339 (2020). https://doi.org/10.1038/s41467-020-15829-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-15829-1
This article is cited by
-
Charge transport and hydrodynamics in materials
Nature Reviews Materials (2023)
-
Para-hydrodynamics from weak surface scattering in ultraclean thin flakes
Nature Communications (2023)
-
Ultra-broadband photoconductivity in twisted graphene heterostructures with large responsivity
Nature Photonics (2023)
-
Giant magnetoresistance of Dirac plasma in high-mobility graphene
Nature (2023)
-
Direct observation of vortices in an electron fluid
Nature (2022)
Comments
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.