Abstract
Stacking order plays a crucial role in determining the crystal symmetry and has significant impacts on electronic, optical, magnetic, and topological properties. Electron-phonon coupling, which is central to a wide range of intriguing quantum phenomena, is expected to be intricately connected with stacking order. Understanding the stacking order-dependent electron-phonon coupling is essential for understanding peculiar physical phenomena associated with electron-phonon coupling, such as superconductivity and charge density waves. In this study, we investigate the effect of stacking order on electron-infrared phonon coupling in graphene trilayers. By using gate-tunable Raman spectroscopy and excitation frequency-dependent near-field infrared nanoscopy, we show that rhombohedral ABC-stacked trilayer graphene has a significant electron-infrared phonon coupling strength. Our findings provide novel insights into the superconductivity and other fundamental physical properties of rhombohedral ABC-stacked trilayer graphene, and can enable nondestructive and high-throughput imaging of trilayer graphene stacking order using Raman scattering.
Similar content being viewed by others
Introduction
Stacking order is a unique structural degree of freedom presented in two-dimensional (2D) layered materials. It plays a significant role in governing symmetry breaking and a wide range of fascinating electronic, optical, magnetic, and topological phenomena. For instance, 3R-stacked transition metal dichalcogenides (TMDs) lack spatial inversion symmetry and exhibit intriguing valley physics and second-order nonlinear responses such as electric-dipole-allowed second harmonic generation (SHG). In contrast, inversion symmetry is restored in 2H-stacked TMD bilayers, resulting in nil valley and even-order nonlinear responses1,2. ABA-stacked (Bernal) and ABC-stacked (rhombohedral) trilayer graphene have distinct electronic and optical properties3,4,5,6,7,8,9,10. ABC-stacked trilayer graphene exhibits an electric-field-tunable band gap and features van Hove singularities at or near the band edge, where the density of states diverges. It can also exhibit tunable Mott insulator11, superconductor12,13, and ferromagnetic14,15 behavior, while these phenomena have not been observed in ABA-stacked trilayer graphene16,17.
Electron–phonon coupling is a fundamental interaction between elementary excitations that plays a significant role in a variety of physical phenomena and quantum phase transitions in condensed matter physics. For example, electron–phonon coupling sets the intrinsic limit of electron mobility and leads to ultra-low thermal conductivity18,19,20,21. It also lays a foundation for charge density waves22,23, electron hydrodynamics24,25,26, superfluidity27,28, and superconductivity29,30,31. As both the electronic band structures and phonon dispersions depend strongly on stacking geometry, a stacking order-governed electron–phonon coupling is expected in principle. Understanding the dependence of electron–phonon coupling on stacking order is fundamentally important for both comprehending a rich variety of peculiar physical phenomena and engineering novel device applications.
In this study, we report the observation of stacking order-governed electron–infrared phonon coupling in trilayer graphene. Specifically, we demonstrate that ABC-stacked trilayer graphene has strong electron–infrared phonon coupling strength. Our findings are supported by gate-tunable Raman spectroscopy and excitation frequency-dependent near-field infrared spectroscopy. The observed giant electron–infrared phonon coupling in ABC-stacked trilayer graphene sheds new light on its peculiar quantum properties, such as strong correlation, superconductivity, and ferromagnetism. Moreover, the distinct electron–infrared phonon coupling between ABC-stacked and ABA-stacked trilayer graphene offers a nondestructive and high-throughput imaging approach based on Raman scattering for identifying the stacking order of trilayer graphene.
The stacking order of trilayer graphene
Trilayer graphene possesses two distinct stacking geometries that have different symmetries and electronic properties. The ABC trilayer graphene is centrosymmetric and semiconducting, while the ABA trilayer graphene is non-centrosymmetric and semimetallic (Fig. 1a and b). This unique characteristic of trilayer graphene makes it an exciting platform for studying stacking order-driven electron–phonon coupling. Both ABC and ABA trilayer graphene can be directly obtained by mechanical exfoliation from bulk crystals as they belong to structurally metastable/stable carbon allotropes. A representative trilayer graphene region with an area >5000 μm2 (highlighted by the dashed black line) on 285 nm-SiO2/Si substrates can be seen in the white-light microscopic image (Fig. 1c). However, having the same optical contrast, the ABC and ABA trilayer graphene cannot be distinguished from the optical image. Scanning near-field optical microscope (SNOM), on the other hand, can sensitively probe the electronic band structure and optical conductivity of samples, offering an effective technique to resolve the stacking order of graphene32,33,34,35. Figure 1d shows a typical SNOM image taken from a trilayer graphene with an excitation frequency of 940 cm−1. From it, the ABC and ABA trilayer graphene can be clearly distinguished as dark and bright regions, respectively.
Strong electron–infrared phonon coupling in the ABC trilayer graphene
Due to the interlayer coupling and stacking geometry, high-frequency optical phonon modes in graphene couple with different valence and conduction bands, resulting in different electron–phonon interaction strengths. The ABA and ABC trilayer graphene have different point group symmetries D3h and D3d, respectively, corresponding to different mode symmetries: 2E’ + E” for ABA trilayer graphene, and 2Eg + Eu for the ABC trilayer graphene, as shown in Fig. 2a. The E’ mode is both Raman and infrared active, the E” and Eg modes are Raman active, the antisymmetric vibrational mode Eu is infrared active, and can be simultaneously Raman active through inversion symmetry breaking36,37,38,39. Figure 2b shows typical Raman spectra including the G mode and 2D mode of ABC and ABA trilayer graphene at heavy hole doping concentrations (ABC at P++-red line, ABA at P++-blue line) and charge neutral points (ABC at CNP-orange line, ABA at CNP-black line), excited by 532 nm (2.33 eV) laser at a temperature of 10 K. At heavy hole doping concentrations, the ABA trilayer graphene shows a single Raman peak at ∼1588 cm−1, assigned as Raman phonon G mode. In striking contrast, ABC trilayer graphene clearly shows an additional low wave number peak at ∼1572 cm−1, which is absent in ABA trilayer graphene. This characteristic Raman peak is also confirmed by measurements under 633 nm-laser excitation and at 300 K (refer to Fig. S1). Besides, the 2D modes arising from a double-resonant electronic process of ABC and ABA trilayer graphene at P++ have almost identical line shapes. At charge neutral points, the ABC and ABA trilayer graphene shows both a single Raman peak at ∼1580 and ∼1582 cm−1, respectively. The unique low wave number peak of ABC trilayer graphene at CNP disappears. Moreover, the 2D modes of ABC and ABA trilayer graphene at CNP have different line shapes. The ABC trilayer graphene exhibits a more asymmetric and wider shape than the ABA trilayer graphene, with the 2D mode of ABC trilayer graphene including a sharp peak on the left and a flat shoulder on the right, consistent with previous observations5,6. The new Raman peak at ∼1572 cm−1 may stem from an infrared active phonon mode. To confirm this hypothesis, we carried out far-field infrared spectroscopy measurements, and the corresponding spectra are shown in Fig. 2c. Remarkably, the infrared absorption peak matches the low wavenumber Raman peak, suggesting that this low wavenumber Raman peak is the infrared active phonon mode.
As mentioned above, this infrared active phonon mode only appears in the ABC trilayer graphene with a broken inversion symmetry. To further quantitatively access the electron–infrared phonon coupling strength, we thus performed Raman spectra of both ABC and ABA trilayer graphene at different gate voltages (\({V}_{{{{{{\rm{g}}}}}}}\)) in Figs. S3–S4. Then we estimated the carrier density n and Fermi level \({E}_{{{{{{\rm{F}}}}}}}\) from \(n={C}_{{{{{{\rm{g}}}}}}}({V}_{{{{{{\rm{g}}}}}}}-{V}_{{{{{{{\rm{cnp}}}}}}}})/e\) and \({E}_{{{{{{\rm{F}}}}}}}=-{{{{{\mathrm{sgn}}}}}}(n){\hslash}{v}_{{{{{{\rm{F}}}}}}}\sqrt{\left(\pi |n|\right)}\), where \({C}_{{{{{{\rm{g}}}}}}}=115\,{{{{{\rm{{aF}}}}}}}/{\mu{{{{{\rm{m}}}}}}}^{2}\), e, \({V}_{{{{{{{\rm{cnp}}}}}}}}\) and \({v}_{{{{{{\rm{F}}}}}}}\) are the gate capacitance, the electron charge, the gate voltage corresponding to the CNP and the Fermi velocity40,41. Notably, the frequencies \(\omega\) and full-width at half-maximum (FWHM) \(\varGamma\) of phonon G mode of both ABA and ABC trilayer graphene are strongly \(n\) dependent. Figure 2d presents the phonon energy extracted by the Lorentz fitting as a function of \(n\) and \({E}_{{{{{\rm{{F}}}}}}}\). With increasing the hole density, the mode of ABA trilayer graphene and high wavenumber mode of ABC trilayer graphene harden, indicating both are symmetric Raman G mode. In marked contrast, the low wavenumber component of ABC trilayer graphene softens with doping density, confirming the antisymmetric infrared active nature. In addition, \(\varGamma\) of phonon G mode of both ABA and ABC trilayer graphene sharply decreases with increasing \(n\) (Fig. 2e).
The observation of infrared active phonon mode in centrosymmetric ABC trilayer graphene by Raman spectroscopy is quite surprising since the parity is a conserved quantity, and thus Raman and infrared transitions are mutually exclusive. We firstly carried out circular polarization and linear polarization Raman measurements (Fig. S9), revealing that the phonon G mode and infrared active phonon mode have the same circular polarization responses; while the linear polarization responses XX and XY are almost identical, indicating that the G peak splitting we observed is not from the stress42,43 and boundary44,45,46. Besides, Fig. S10 shows that the D peak at 1350 cm−1 is almost invisible, indicating that the defects in the sample are rare, and the split G peak is not from the defects47,48. We thus assign this mode to the infrared active phonon mode caused by the dielectric environment doping, which breaks the inversion symmetry of ABC trilayer graphene. It is worth noting that the infrared active phonon disappears with hole doping below \(\left|n\right|=\left|-10\times {10}^{12}\right|\,{{{\rm {cm}}}}^{-2}\), and can only be clearly seen with hole doping above \(\left|n\right|=\left|-15\times {10}^{12}\right|\,{{{\rm {cm}}}}^{-2}\) (Figs. 2d and S4). In principle, for intrinsic graphene (\({V}_{{{\rm {CNP}}}}=0\,{\rm {V}}\)) on 285-nm thick \({{{\rm {SiO}}}}_{2}\) substrates, a gate voltage of \({V}_{{\rm {g}}}=-210\,{\rm {V}}\) needs to be applied to tune the carrier density up to \(n=-15\times {10}^{12}\,{{{\rm {cm}}}}^{-2}\). Such high gate voltage is not practical as which far exceeds the breakdown limit of the \({{{\rm {SiO}}}}_{2}\) gate dielectric. Thus, we measured two strong hole doping devices and one device that detected the CNP to normalize the carrier density \(n\) and Fermi level \({E}_{{\rm {F}}}\) based on the transport curve and the gate voltage corresponding to the same frequencies \(\omega\) of phonon G mode in Figs. S2–S4.
We quantitatively calculated the electron–phonon coupling strength through the frequencies \(\omega\) and full width at half maximum (FWHM) \(\varGamma\) of phonon G mode. The change of G phonon energy is described by \({\hslash}\omega-{\hslash}{\omega}_{{{{{{{\rm{CNP}}}}}}}}=\lambda\left\{\left|{E}_{{{{{{\rm{F}}}}}}}\right|+\frac{{{\hslash}}\omega}{4}{{{{\mathrm{ln}}}}}\left|\frac{2\left|{E}_{{{{{{\rm{F}}}}}}}\right|-{{\hslash}}\omega}{2\left|{E}_{{{{{\rm{{F}}}}}}}\right|+{{\hslash}}\omega}\right|\right\}\), where \({\omega}_{{{{{{{\rm{CNP}}}}}}}}\) is \(\omega\) at CNP, \(\lambda={A}_{{{{{{{\rm{uc}}}}}}}}{D}^{2}/2\pi{{\hslash}}\omega M{v}_{{{{{{\rm{F}}}}}}}^{2}\). \({A}_{{{{{{{\rm{uc}}}}}}}}\) is the area of the graphene unit cell, \({M}\) is the carbon atom mass, \({v}_{{{{{{{\rm{F}}}}}}}}\) is the Fermi velocity and D is the electron–phonon coupling strength40,41. The fits of the linear segments in Fig. 2d (dashed black lines) give electron–phonon coupling strength \({D}_{{{{{{{\rm{ABA}}}}}}}}\approx 7.0\,{{{{{\rm{{eV}}}}}}}/\mathring{{{{{\rm{A}}}}}}\) for ABA trilayer graphene and \({D}_{{{{{{{\rm{ABC}}}}}}\; {{{{{\rm{G}}}}}}}({E}_{{{{{{\rm{g}}}}}}})}\approx 8.6\,{{{{{{\rm{eV}}}}}}}/\mathring{{{{{\rm{A}}}}}}\) for ABC trilayer graphene \({E}_{{{{{{\rm{g}}}}}}}\) phonon mode. In addition, the electron–phonon coupling strength can also be calculated from \(\Delta \varGamma=\frac{{A}_{{{{{{{\rm{uc}}}}}}}}{D}^{2}}{8M{v}_{{{{{{\rm{F}}}}}}}^{2}}\), \(\Delta \varGamma={\varGamma }_{{{{{{{\rm{CNP}}}}}}}}-{\varGamma }_{0}\), where Γ0 is the residual linewidth from processes that are not related to Landau damping40,41. Figure 2e shows \(\varDelta \varGamma \approx 4\,{{{{{\rm{{{cm}}}}}}}}^{-1}\) for both ABA and ABC trilayer graphene \({E}_{{{{{\rm{{g}}}}}}}\) phonon mode, giving \({D}_{{{{{{{\rm{ABA}}}}}}}}\approx 8.0\,{{{{{{\rm{eV}}}}}}}/\mathring{{{{{\rm{A}}}}}}\) and \({D}_{{{{{{{\rm{ABC}}}}}}\; {{{{{\rm{G}}}}}}}({E}_{{{{{{\rm{g}}}}}}})}\approx 8.0\,{{{{{\rm{{eV}}}}}}}/\mathring{{{{{\rm{A}}}}}}\), consistent with the electron–phonon coupling strength calculated by frequencies \(\omega\); while \(\varDelta \varGamma \approx 7{\,{{{{{\rm{{cm}}}}}}}}^{-1}\) for ABC trilayer graphene \({E}_{{{{{\rm{{u}}}}}}}\) infrared active phonon mode gives \({D}_{{{{{{{{{\rm{ABC}}}}}}\; {{{{{\rm{G}}}}}}}}}({E}_{{{{{{\rm{u}}}}}}})}\approx 10.6\,{{{{{\rm{{eV}}}}}}}/\mathring{{{{{\rm{A}}}}}}\). Note that we did not see the infrared phonon in the Raman spectroscopy of ABA trilayer graphene, so we cannot directly compare the coupling strength of electron–infrared phonon between ABA and ABC trilayer graphene. But we can see that the electron–phonon coupling strength \({D}_{{{{{{\rm{{ABA}}}}}}}}\approx 8.0\,{{{{{{\rm{eV}}}}}}}/\mathring{{{{{\rm{A}}}}}}\) and \({D}_{{{{{{{\rm{ABC}}}}}}\; {{{{{\rm{G}}}}}}}({E}_{{{{{{\rm{g}}}}}}})}\approx 8.0\,{{{{{\rm{{eV}}}}}}}/\mathring{{{{{\rm{A}}}}}}\) are almost the same but smaller than \({D}_{{{{{{{\rm{ABC}}}}}}\; {{{{{\rm{G}}}}}}}({E}_{{{{{{\rm{u}}}}}}})}\approx 10.6\,{{{{{\rm{{eV}}}}}}}/\mathring{{{{{\rm{A}}}}}}\). In addition, the \({D}_{{{{{{{\rm{ABC}}}}}}\; {{{{{\rm{G}}}}}}}({E}_{u})}\approx 10.6\,{{{{{\rm{{eV}}}}}}}/\mathring{A}\) is also stronger than the electron–phonon coupling strength 6.4\({{{{{{\rm{eV}}}}}}}/\mathring{{A}}\) in bilayer graphene38.
Notably, recent studies have revealed superconductivity in ABC trilayer graphene12,13 but not in ABA trilayer graphene16,17. Although superconductivity has also been observed in bilayer graphene in the latest research, it needs to be induced by either parallel magnetic fields49 or WSe2 proximity-induced spin–orbit coupling50. The stronger electron–infrared phonon coupling in ABC trilayer graphene may provide new insights into the superconductivity, as electron–phonon coupling may play a significant role.
Resonance with the infrared active phonon in ABC trilayer graphene
We performed SNOM measurements using lasers of different frequencies to confirm the strong electron–infrared phonon coupling in ABC trilayer graphene. As shown in Fig. 3a, when the laser frequency is not at the phonon resonance frequency, the SNOM experiment reveals a low scattering amplitude intensity in ABC trilayer graphene compared to ABA trilayer graphene, indicating weak conductivity. However, when the excitation laser frequency is resonant with the infrared active phonon, the near-field signal should be dominated by the infrared phonon response if the electron–infrared phonon coupling is strong. Figure 3b shows the SNOM image excited by laser radiation resonant with the infrared active phonon at 1572 cm−1 (also see Fig. S7). We performed variable frequency spectroscopy measurements (i.e., point spectroscopy measurements) on different ABA and ABC trilayer graphene regions (labeled as 1–6 in Fig. 3b). All region signal intensities were normalized using the scattering intensity of the SiO2 substrate. We conclude that ABC has a stronger response intensity than ABA at near-infrared phonons frequencies in Fig. 3c. Electron–phonon interactions primarily arise from the coupling between phonon resonances and interband transitions, leading to the transfer of interband resonances to phonon resonances. These results suggest that the electron–infrared phonon coupling strength in ABC trilayer graphene is significantly stronger than that in ABA trilayer graphene.
Identifying the stacking order based on the electron–infrared phonon coupling
Previously, techniques for identifying the stacking order of graphene include SNOM32,33,34,35, SHG51, scanning Kelvin probe microscopy (SKPM)52, and conductive atomic force microscopy (CFM)53,54. Thanks to the activated infrared phonon in ABC trilayer graphene, Raman spectroscopy provides an alternative and high-throughput technique for imaging the stacking order of trilayer graphene. Figure 4a shows an optical image of mechanically exfoliated trilayer graphene. Figure 4b presents the SNOM image taken from the trilayer graphene within the dotted square area shown in Fig. 4a. ABC and ABA trilayer graphene can be clearly identified by the dark and bright regions, respectively. Figure 4c shows the Raman intensity mapping of infrared active phonon mode taken from the same region as in Fig. 4b. Obviously, the Raman mapping matches well with the SNOM mapping in terms of distinguishing the ABC and ABA regions. The Raman mapping provides a fast, simple, non-destructive, and high-throughput technology for identifying the stacking order of multilayer graphene.
In summary, our research has shown that the electron–infrared phonon coupling strength in ABC trilayer graphene is significantly stronger than in ABA trilayer graphene and AB bilayer graphene. This discovery provides new perspectives on the understanding of superconductivity and other physical properties of ABC trilayer graphene. Additionally, the differences in electron–infrared phonon coupling between ABC and ABA graphene can be utilized to develop a novel, non-destructive, and high-throughput Raman scattering imaging technique for identifying the stacking order of multilayer graphene. It also provides a new perspective for exploring and searching for new superconducting systems in graphene, such as rhombohedral multilayer graphene and twisted multilayer graphene.
Methods
Device fabrication
The few-layer graphene flakes were mechanically exfoliated by tape from the graphite bulk crystal and then transferred to a substrate of 285 nm-SiO2/Si (P++). In order to remove adsorbents on the substrate surface and ensure cleanliness of the interface, the SiO2/Si (P++) substrate was ultrasonically cleaned for 2 h in acetone and deionized water, and then subjected to oxygen plasma before the process of mechanical exfoliation. Then, the transferred graphene samples were subjected to high-temperature (450 °C) annealing in an Ar/H2 atmosphere for 12 h to ensure a clean surface of the graphene samples. The electrode contact was carried out by coating silver glue or evaporating Cr/Au by electron-beam metal evaporation on the edge of the graphene sample. Directly applied gate voltage through Si (P++) substrate and used 285 nm-SiO2 as the dielectric layer. The hexagonal boron nitride h-BN encapsulated graphene device is prepared through dry transfer, and a standard electron-beam lithography process and electron-beam metal evaporation.
Raman measurement
Raman spectra are acquired using a micro-Raman spectrometer (Horiba LabRAM HR Evolution) in a confocal backscattering geometry. A solid-state laser at 532/633 nm is focused onto the samples along the z direction by a ×50/×100 objective. The backscattered signal is collected by the same objective and dispersed by a 600/1800-groove mm−1 grating. The laser power during Raman measurement is kept below 100 μW in order to avoid sample damage and excessive heating. The integration time is 10–60 s. Low-temperature (10 K) measurements were performed using a closed-cycle optical cryostat. Electrical measurements were performed in a probe station and under dark conditions with semiconductor parameter analyzers.
Near- and far-field infrared measurement
A scattering scanning near-field optical microscope (Neaspec) equipped with a wavelength-tunable quantum cascade laser (890–2000 cm−1) was used to image optical near fields. A metalized cantilever atomic force microscope tip served as a scattering near-field probe. The tip was illuminated with monochromatic p-polarized infrared light from a quantum cascade laser. We utilized traditional Pt metal-coated probes for the experiment. The tapping frequency of the probe was adjusted to approximately 270 kHz, while the amplitude was set to around 100 nm. Far- field infrared measurements were performed by FTIR microscopy (Thermo Fisher Nicolet iN10).
Data availability
The data that support the findings of this study are available from the corresponding authors on a request.
References
Suzuki, R. et al. Valley-dependent spin polarization in bulk MoS2 with broken inversion symmetry. Nat. Nanotechnol. 9, 611–7 (2014).
Xu, X. et al. Towards compact phase-matched and wave guided nonlinear optics in atomically layered semiconductors. Nat. Photonics 16, 698–706 (2022).
Lui, C. et al. Observation of an electrically tunable band gap in trilayer grapheme. Nat. Phys. 7, 944–947 (2011).
Zou, K. et al. Transport studies of dual-gated ABC and ABA trilayer graphene: band gap opening and band structure tuning in very large perpendicular electric fields. Nano Lett. 13, 369–73 (2013).
Lui, C. et al. Imaging stacking order in few-layer graphene. Nano Lett. 11, 164–9 (2011).
Cong, C. et al. Raman characterization of ABA- and ABC-stacked trilayer graphene. ACS Nano 5, 8760–8768 (2011).
Lui, C. H. et al. Tunable infrared phonon anomalies in trilayer graphene. Phys. Rev. Lett. 110, 185504 (2013).
Qian, J. et al. Nonequilibrium phonon tuning and mapping in few-layer graphene with infrared nanoscopy. Phys. Rev. B 103, L201407 (2021).
Luan, Y. et al. Imaging stacking-dependent surface plasmon polaritons in trilayer graphene. Phys. Rev. Appl. 18, 024052 (2022).
Yang, J. et al. Spectroscopy signatures of electron correlations in a trilayer graphene/hBN moiré superlattice. Science 375, 6586–1299 (2022).
Chen, G. et al. Evidence of a gate-tunable Mott insulator in a trilayer graphene moiré superlattice. Nat. Phys. 15, 237–241 (2019).
Zhou, H. et al. Superconductivity in rhombohedral trilayer graphene. Nature 598, 434–438 (2021).
Chen, G. et al. Signatures of tunable superconductivity in a trilayer graphene moiré superlattice. Nature 572, 215–219 (2019).
Zhou, H. et al. Half and quarter metals in rhombohedral trilayer graphene. Nature 598, 429–433 (2021).
Chen, G. et al. Tunable correlated Chern insulator and ferromagnetism in a moiré superlattice. Nature 579, 56–61 (2020).
Zibrov, A. A. et al. Emergent Dirac gullies and gully-symmetry-breaking quantum Hall states in ABA trilayer graphene. Phys. Rev. Lett. 121, 167601 (2018).
Winterer, F. et al. Spontaneous gully-polarized quantum Hall states in ABA trilayer graphene. Nano Lett. 22, 3317–3322 (2022).
Li, C. W. et al. Orbitally driven giant phonon anharmonicity in SnSe. Nat. Phys. 11, 1063–1069 (2015).
Dean, C. R. et al. Boron nitride substrates for high-quality graphene electronics. Nat. Nanotechnol. 5, 722–726 (2010).
Wang, L. et al. One-dimensional electrical contact to a two-dimensional material. Science 342, 614–617 (2013).
Sadeghi, M. M. et al. Tunable electron–flexural phonon interaction in graphene heterostructures. Nature 617, 282–286 (2023).
Luo, H. et al. Electronic nature of charge density wave and electron-phonon coupling in kagome superconductor KV3Sb5. Nat.comm 13, 1 (2022).
Quan, J. et al. Phonon renormalization in reconstructed MoS2 moiré superlattices. Nat. Mater 20, 1100 (2021).
Bandurin, D. A. et al. Negative local resistance caused by viscous electron backflow in graphene. Science 351, 1055–1058 (2016).
Crossno, J. et al. Observation of the Dirac fluid and the breakdown of the Wiedemann–Franz law in graphene. Science 351, 1058–1061 (2016).
Krishna Kumar, R. et al. Superballistic flow of viscous electron fluid through graphene constrictions. Nat. Phys. 13, 1182–1185 (2017).
Liu, X. et al. Quantum Hall drag of exciton condensate in graphene. Nat. Phys. 13, 746–750 (2017).
Li, J. I. A. et al. Excitonic superfluid phase in double bilayer graphene. Nat. Phys. 13, 751–755 (2017).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018).
Gadelha, A. C. et al. Localization of lattice dynamics in low-angle twisted bilayer graphene. Nature 590, 405 (2021).
Drozdov, A. P. et al. Conventional superconductivity at 203 kelvin at high pressures in the sulfur hydride system. Nature 525, 73 (2015).
Sunku, S. S. et al. Photonic crystals for nano-light in moiré graphene superlattices. Science 362, 1153–1156 (2018).
Ju, L. et al. Topological valley transport at bilayer grapheme domain walls. Nature 520, 650 (2015).
Jiang, L. et al. Manipulation of domain-wall solitons in bi-and trilayer graphene. Nat. Nanotechnol. 13, 204–208 (2018).
Li, H. et al. Global control of stacking-order phase transition by doping and electric field in few-layer graphene. Nano Lett. 20, 3106–3112 (2020).
Yan, J. A. et al. Phonon dispersions and vibrational properties of monolayer, bilayer, and trilayer graphene: density-functional perturbation theory. Phys. Rev. B 77, 125401 (2008).
Ando, T. et al. Field effects on optical phonons in bilayer graphene. J. Phys. Soc. Jpn. 78, 034709 (2009).
Malard, L. M. et al. Observation of distinct electron–phonon couplings in gated bilayer graphene. Phys. Rev. Lett. 101, 257401 (2008).
Yan, J. et al. Optical phonon mixing in bilayer graphene with a broken inversion symmetry. Phys. Rev. B 80, 241417 (2009).
Yan, J. et al. Electric field effect tuning of electron–phonon coupling in graphene. Phys. Rev. Lett. 98, 166802 (2007).
Yan, J. et al. Raman scattering and tunable electron-phonon coupling in single layer graphene. Solid State Commun. 143, 39–43 (2007).
Mohiuddin, T. et al. Uniaxial strain in graphene by Raman spectroscopy: G peak splitting, Grüneisen parameters, and sample orientation. Phys. Rev. B 79, 205433 (2009).
Frank, O. et al. Raman 2D-band splitting in graphene: theory and experiment. ACS Nano 5, 2231–2239 (2011).
Yang, R. et al. Observation of Raman G-peak split for graphene nanoribbons with hydrogen-terminated zigzag edges. Nano Lett. 11, 4083–4088 (2011).
Cançado, L. G. et al. Anisotropy of the Raman spectra of nanographite ribbons. Phys. Rev. Lett. 93, 047403 (2004).
Ren, W. et al. Edge phonon state of mono- and few-layer graphene nanoribbons observed by surface and interference co-enhanced Raman spectroscopy. Phys. Rev. B 81, 035412 (2010).
Eckmann, A. et al. Probing the nature of defects in graphene by Raman spectroscopy. Nano Lett. 12, 3925–3930 (2012).
Venezuela, P. et al. Theory of double-resonant Raman spectra in graphene: Intensity and line shape of defect-induced and two-phonon bands. Phys. Rev. B 84, 035433 (2011).
Zhou, H. et al. Isospin magnetism and spin-polarized superconductivity in Bernal bilayer graphene. Science 375, 774–778 (2022).
Zhang, Y. et al. Enhanced superconductivity in spin–orbit proximitized bilayer graphene. Nature 613, 268–273 (2023).
Shan, Y. et al. Stacking symmetry governed second harmonic generation in graphene trilayers. Sci. Adv. 4, eaat0074 (2018).
Yu, J. et al. Imaging graphene Moiré superlattices via scanning Kelvin probe microscopy. Nano Lett. 21, 3280–3286 (2021).
Zhang, S. et al. Domino-like stacking order switching in twisted monolayer–multilayer graphene. Nat. Mater. 21, 621–626 (2022).
Li, Y. et al. Symmetry breaking and anomalous conductivity in a double moiré superlattice. Nano Lett. 22, 6215–6222 (2022).
Acknowledgements
This work is supported by the National Key Research and Development Program of China (Grant Nos. 2023YFA1407000, 2021YFA1202900, 2020YFA0309600), National Science Foundation of China (NSFC, Grant Nos. 12274447, 61888102, 51925203, 52102160, 51972074), the Strategic Priority Research Program of CAS (Grant No. XDB0470101) and the Key-Area Research and Development Program of Guangdong Province (Grant No. 2020B0101340001).
Author information
Authors and Affiliations
Contributions
G.Z. supervised the project. X.Z. and X.G. designed the experiments. X.Z. fabricated the devices, performed the Raman spectroscopy and transport measurements with the help of Z.H. X.G., X.Y., and Q.D. carried out far-field infrared spectroscopy and SNOM variable frequency measurements. A.D. and Z.S. performed SNOM characterization measurements. X.Z. analyzed the data and prepared the figures. X.Z., X.G., L.D., and G.Z. wrote the paper. L.L., F.W., Y.Y., J.Z., Y.P., L.L., Y.Z., X.L., J.Z., J.D., W.Y., and D.S. were involved in discussions of this work.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Tiancheng Song, Liuyan Zhao, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zan, X., Guo, X., Deng, A. et al. Electron/infrared-phonon coupling in ABC trilayer graphene. Nat Commun 15, 1888 (2024). https://doi.org/10.1038/s41467-024-46129-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-46129-7
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.