Key Laboratory of Artificial Structures and Quantum Control (Ministry of Education), Shenyang National Laboratory for Materials Science, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, People’s Republic of China
2.
Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, People’s Republic of China
3.
Institute of System Engineering, Beijing 100191, People’s Republic of China
ABC-stacked trilayer graphene has exhibited a variety of correlated phenomena owing to its relatively flat bands and gate-tunable bandgap. However, convenient methods are still lacking for identifying ABC graphene with nanometer-scale resolution. Here we demonstrate that the scanning near-field optical microscope working in ambient conditions can provide quick recognition of ABC trilayer graphene with no ambiguity and excellent resolution (20 nm). The recognition is based on the difference in their near-field infrared (IR) responses between the ABA and ABC trilayers. We show that in most frequencies, the response of the ABC trilayer is weaker than the ABA trilayer. However, near the graphene phonon frequency (1585 cm-1), ABC’s response increases dramatically when gated and exhibits a narrow and sharp Fano-shape resonant line, whereas the ABA trilayer is largely featherless. Consequently, the IR contrast between ABC and ABA becomes reversed and can even be striking (ABC/ABA 3) near the graphene phonon frequency. The observed near-field IR features can serve as a golden rule to quickly distinguish ABA and ABC trilayers with no ambiguity, which could largely advance the exploration of correlation physics in ABC-stacked trilayer graphene.
Future perspectivesFinding the accurate locations of ABC-stacked domains in trilayer graphene at nanometer scale is an important prerequisite for making ABC nano-devices to study their unique properties. The quick identification with nanoscale resolution proposed here using scanning near-field optical microscopy in ambient conditions will definitely reduce the barrier in fabricating ABC trilayer graphene nano-devices, and therefore greatly promote the exploration of novel correlated physics in ABC trilayers. Moreover, the observed near-field IR Fano resonance exhibits a very narrow peak at the graphene phonon mode. Considering that graphene phonon frequency is very sensitive to strains (1% change in lattice constant will induce 50 cm-1 change in graphene phonon frequency), we propose that the observed Fano resonance could also be used to detect local strains in ABC trilayer graphene with high sensitivity and high spatial resolution. As strains normally exist in low-dimensional materials and can significantly modify the materials’ intrinsic properties, it would be of great help to map out the strain distribution in ABC graphene to fully understand its novel correlated physics.
1.
Introduction
Trilayer graphene has two kinds of stacking orders: the Bernal ABA and the rhombohedral ABC. The electronic properties of ABC-stacked trilayer graphene are significantly different from those of the ABA-stacked trilayer due to its flat bands [1-5] and gate-tunable bandgap [4, 6, 7], leading to remarkable quantum and correlated phenomena [2, 3, 8-10], such as superconductivity, Mott insulator, correlated Chern insulator, and ferromagnetism, which have caught great research interests. The preparation of ABC graphene samples and especially finding the accurate locations of ABC domains is a prerequisite for exploring their unique correlated electronic properties. However, up to now, it is still not easy to find the exact locations of ABC-stacked domains and make ABC graphene devices, partially due to the relatively rare distribution of ABC graphene, and partially due to the lack of a convenient way to identify ABC graphene with nanometer-scale resolution. This has strongly hampered further extensive exploration of new correlated physics in ABC trilayer graphene.
The most widely used techniques to identify ABC graphene are traditional Raman spectroscopy [11-15] and infrared (IR) spectroscopy [16]. However, both of the above approaches are of low resolution due to the far-field optical diffraction limit. On the other hand, some techniques with high resolution, such as scanning tunneling microscopy and transmission electron microscopy, although can distinguish ABA and ABC graphene [17, 18], cannot serve as convenient tools due to either the requirement of ultrahigh vacuum or difficulties in making suspended samples. Some atomic force microscope (AFM) based techniques can distinguish ABA and ABC graphene at the nanometer scale, but their mechanism remains elusive [19, 20]. Recently, scanning near-field optical microscopy (SNOM) has acted as a powerful tool to study the optical properties of low-dimensional materials with high resolution (20 nm) and good efficiency [21-28]. A lot of nano-structures and phenomena can be directly detected in real space by the SNOM technique, such as plasmon polaritons in graphene and carbon nanotubes [21, 22, 26, 28, 29], phonon polaritons in hexagonal boron nitride (hBN) [30] and -MoO3 [31, 32], domain walls in few layers graphene [24, 29], coupling of polaritons between different nano-materials [33, 34], nanoscale local strains in bilayer graphene [35], hBN [36] and nano-scale semiconductors [37], etc. The SNOM has also been used to study ABC graphene domains in previous studies, as ABC and ABA graphene typically have different IR conductivities and therefore show different near-field IR contrasts [38]. However, it is still difficult to surely identify the ABC domains with the single SNOM technique. Typically, Raman spectroscopy or statistical analysis is used to further confirm the stacking orders.
It was reported that a perpendicular electric field can break the symmetry of rhombohedral-stacked few-layer graphene, and introduce a band gap as well as electric dipole moments [6, 39]. Consequently, the initially IR inactive phonons at 1585 cm-1 in bilayer graphene will become IR active through electron-phonon coupling [35], and a Fano system is formed when such electron-dressed phonons interact with IR photons, resulting in a strong IR response at the resonant frequency [35]. The properties of the ABC trilayer are similar to the bilayer graphene because both of them have inversion symmetry yet lack mirror symmetry [6]. In recent years, a far-field IR transmission spectroscopic study showed that the ABC trilayer also has obvious features near its optical phonon frequency [40], similar to the bilayer. However, the Bernal-stacked ABA trilayer appears largely featureless due to its mirror symmetry. The different IR responses provide possibilities to distinguish ABA and ABC trilayers.
Here we report quick identification of ABC trilayer graphene with nanoscale resolution using near-field IR nano-imaging. The identification is based on the fact that ABA and ABC graphene exhibit distinct near-field IR features: (1) in a wide range of excitation frequencies, ABC graphene responds weaker than ABA graphene due to ABC’s lower optical conductivity from electrons; and (2) near its optical phonon frequency (1585 cm-1), ABC graphene responses much stronger than ABA due to ABC’s higher phonon activity. A largely enhanced Fano-resonance peak was observed on ABC domains when subjected to a vertical electric field, while ABA appears more inert. These fingerprint features can be used to distinguish and identify ABC domains with no ambiguity. Moreover, the near-field IR spectroscopy can break the diffraction limit of light and reach an excellent spatial resolution of 20 nm, which is ideal for further device fabrication at the nanometer scale and can promote the exploration of ABC’s novel properties.
2.
Results
Figure 1(a) shows the atomic structures of ABA- and ABC-stacked trilayer graphene. The ABA trilayer has mirror symmetry while ABC does not, and ABA has four-fold degeneration at the charge neutral point while ABC has only two-fold, as figure 1(b) illustrates. The degeneration in ABC can be lifted by the external electrostatic field, and consequently a band-gap appears [6, 41, 42]. Nevertheless, ABA does not show such significant change under vertical fields owing to its mirror symmetry structure.
Figure
1.
The atomic structures of ABA- and ABC-stacked trilayer graphene and their identifications via scanning near-field optical microscopy (SNOM). (a) Atomic structure diagrams of ABA and ABC graphene. The red and blue balls represent A and B atoms, respectively. The bottom two pictures are the lateral views along the graphene armchair direction, as denoted by the black dashed lines. (b) The electronic bands of intrinsic ABA and ABC trilayers. (c) Illustration of the experiment setup. An infrared (IR) beam is illuminated onto a conductive AFM tip and the graphene sample, and the scattered light is collected to analyze the stacking orders of the trilayer graphene. The inset shows the electrostatic configuration. The tip and graphene keep grounded and the bottom voltage Vg is applied between the silicon substrate and graphene. (d)-(e) SNOM images of the same graphene region at Vg = -80 V, where the excitation frequency is 1480 cm-1 (far away from the graphene phonon frequency) in (d) and 1579.5 cm-1 (near the optical phonon frequency) in (e). The scale bar is 500 nm.
A homemade SNOM system built from a commercial AFM is used in our experiment. Figure 1(c) illustrates the near-field IR characterization of a trilayer graphene flake with different stacking domains. In short, a beam of quantum cascade laser (QCL) with changeable frequency is focused on a gold-coated AFM tip, and near-field signals return along the same path and are collected by a mercury cadmium telluride detector placed in the far field. Our graphene samples are mechanically exfoliated from bulk graphite onto the silicon substrate with a 285 nm-thick oxide layer. The thickness and stacking orders of the graphene samples are preliminarily characterized by optical microscopy, AFM, and Raman spectroscopy (corresponding data can be found in figures S1 and S2). In our SNOM measurements, a gate voltage is typically applied to tune the charge density of graphene, and more importantly to provide a vertical electrical field to activate the graphene phonon. The inset of figure 1(c) shows the electrostatic configuration of the back gate.
Figures 1(d) and (e) display two near-field IR images of the same graphene sample at two excitation frequencies of 1480 cm-1 and 1579.5 cm-1, respectively. Both images show two different kinds of domains in the sample, but the IR contrast of the two domain types is reversed at the two excitations. At 1480 cm-1 (figure 1(d)), the ABA is brighter than ABC. That is probably because ABA has higher optical conductivity than ABC, leading to its stronger near-field response. While at 1579.5 cm-1 (figure 1(e)) close to the graphene phonon, ABC becomes prominently brighter than ABA, due to the strong phonon response of ABC graphene (more details will be discussed later). These features indicate the capability of SNOM for identifying the stacking orders of trilayer graphene.
Next, we systematically investigated near-field IR responses of both ABA and ABC trilayers in a wide frequency range (about 940 cm-1-1630 cm-1). We found that in most frequencies, ABA domains respond strongly and display brighter than ABC graphene (figures 2(a)-(d) and (h)). This can be explained by the fact that the ABA trilayer possesses larger Drude conductivity than the ABC trilayer in the far infrared range, due to their different electronic bands. In the low energy range, ABA graphene contains a pair of linear bands with zero mass and a pair of parabolic bands (figure 1(b)) with a relatively small effective mass of 0.063me, whereas ABC graphene contains only a pair of cubic bands with a relatively large effective mass of 0.14me, where me is the free electron mass [43]. Considering that the Drude conductivity is inversely proportional to the effective mass [44], σ∝1/m, ABA’s Drude conductivity is therefore significantly greater than the ABC trilayers, leading to a brighter ABA contrast in most tested far-infrared frequencies.
Figure
2.
Near-field spectral response of ABA and ABC trilayers and the Fano resonance. (a)-(h) A series of near-field IR images taken at different excitation frequencies, from which obvious changes in the ABC domain can be seen. The measurements are performed at Vg = +80 V, and the scale bar is 500 nm. (i) The near field IR intensity of ABC trilayer (normalized using ABA graphene’s response) against photon frequency, behaving as a typical Fano-shape line with key Fano parameter q≈6.1. The inset is the normalized near-field IR intensity of ABC graphene in an extended range from 940 cm-1 to 1630 cm-1. The Fano fitting is carried out in a smaller range near the graphene phonon as denoted by the orange arrow.
However, the near-field IR intensity of ABC increases rapidly when incident light’s frequency is near the graphene intrinsic phonons, as shown in figures 2(e)-(g). This can be explained by the stronger IR activity of the ABC graphene phonon. Noted that the measurement is performed under a back gate of +80 V. The applied gate can provide a vertical electrical field that breaks the symmetry and endows electrical dipole moments into the ABC graphene phonon. This enables the initially IR inactive optical phonons of ABC to become IR active through electron-phonon coupling. The electron-dressed phonons can interact with IR photons directly and show significantly strong near-field signals. While for the ABA domains, the phonon always keeps silent to infrared light due to mirror symmetry.
We further plot the normalized near-field IR intensity of ABC as a function of excitation frequency in figure 2(i) (the inset shows more experimental data in an extended frequency range from 940 cm-1 to 1630 cm-1), which shows a prominent and sharp Fano peak at the graphene phonon frequency, the result of interference between dressed-phonon’s and electron’s transitions. When absorbing the energy from incident free IR photons, the transition of such electron-dressed phonon couples with continuous electron transitions, forming a typical Fano shape near-field signal line with a high slope [35]. We then fit data in figure 2(i) with the famous Fano formula [45]
A(E)=A0[qγ+E−E0]2(E−E0)2+γ2,
where E is the incident photon energy, E0 is the central energy, γ is the linewidth, A0 is a constant, and the dimensionless Fano parameter q describes the shape of Fano line, which determines whether electron’s or dressed-phonon’s transition plays the determining role. By fitting the experimental data, we got q≈6.1, a bit larger than the unit, meaning that the dressed phonon’s transitions are dominant, instead of the electron’s transitions. Noted that the peak is very narrow (2.5 cm-1) and asymmetric, which are the two merits of the Fano resonance (Fano fitting properties of more samples can be found in section 2 of supporting materials). The observed Fano resonance peak near the graphene phonon frequency can serve as a fingerprint to identify the ABC trilayer.
In the following, we study the influence of the gate voltage on the near-field IR signals of the trilayer graphene at two representative frequencies, 1595 cm-1 and 1583 cm-1. At 1595 cm-1 (out of the Fano resonance region), ABC always shows a weaker response than ABA at all gate voltages as shown in figures 3(a)-(c). Because at this frequency, phonon transitions almost do not contribute, and the near-field signals are mainly from the electronic contributions. As discussed above, ABC’s IR conductivity is lower than the ABA. Consequently, ABC is darker than ABA in the near-field IR images. In figure 3(g), we plot the near-field intensity of both ABA and ABC domains as a function of gate voltage. The SiO2 signal as the stable reference almost kept a constant, reflecting that the incident light and the electric system were stable during the whole measurement. Although the absolute intensities of ABA and ABC are different, they show the same tendency against the gate voltage Vg. The near-field intensity features a maximum at charge neutral, and decreases slightly under high gate voltage on both the electron and hole doping sides. This is probably because the inter-band transitions of electrons are suppressed when the Fermi level is lifted away from the charge-neutral state, as illustrated in figure S3.
Figure
3.
Gate dependences of the near-field IR responses of ABA and ABC trilayers at both non-resonant and resonant frequencies. (a)-(c) Near field IR images of a trilayer graphene sample at 1595 cm-1 (beyond resonant region) with effective bottom gate voltage Vg-VCN = 10 V, 20 V, 40 V, respectively. The locations of ABC and ABA are denoted in (a). (d)-(f) Near field IR intensity images of the same region in 1583 cm-1 (within Fano resonant region) at Vg-VCN = 10 V, 20 V, 40 V, respectively. Note that the contrast between ABC and ABA is reversed from (d) to (f), the result of gradually enhanced electron and phonon coupling. Scale bar: 500 nm. (g)-(h) The near-field IR responses of ABC and ABA trilayers under different effective gate Vg-VCN in 1595 cm-1 and 1583 cm-1, respectively. The SiO2 signal as the stable reference is also shown in the figures. (i) The ratio ABC/ABA under different effective gate voltages. The gray and pink dots represent data at 1595 cm-1 and 1583 cm-1, respectively.
Different gate dependence shows up for another excitation frequency of 1583 cm-1 (within the Fano resonance region) in figures 3(d)-(f). Although the ABC signal is weaker than the ABA near the charge neutral state (figure 3(d)), it increases quickly and shows the same contrast as the ABA at a mediate effective gate of Vg-VCN 20 V (figure 3(e)) and becomes much brighter at a higher effective gate voltage of Vg-VCN 40 V (figure 3(f)). The systematic changes of both stacking orders at the incident of 1583 cm-1 were plotted in figure 3(h). The ABA’s signal varies slightly with the gate voltage, showing the same gate dependence as that at 1595 cm-1. On the contrary, the ABC’s signal is very sensitive and largely enhanced at high gate voltages (especially the electron-charged side), which is the direct evidence of phonon-electron coupling [38]. More measurement results (see figure S4) show a similar gate dependence and confirm the phenomena observed in figure 3(h) is universal.
The near-field signal ratio of ABC to ABA is plotted in figure 3(i) as a function of effective gate voltage Vg-VCN at both 1595 cm-1 and 1583 cm-1. At 1595 cm-1, the ratio is roughly a constant (<1) for all gate voltages. At 1583 cm-1, the ratio is less than one only near the charge neutral state and increases to greater than one with gradually increased gate voltage. At a sufficient voltage, the ratio can raise to 2.77, reflecting a sharp contrast between the two stacking orders. This gate-dependent behavior could serve as a sensitive and reliable method to distinguish ABC and ABA domains.
To provide a full picture of ABC and ABA’s near-field response, we simultaneously tuned both the gate voltage and the incident light frequency, and measured the near-field signals of both ABC and ABA trilayers, the results of which are shown in figure 4. Figures 4(a) and (b) display the 2D maps for ABA and ABC’s near-field signals with excitation frequency varied in a range covered the graphene phonon frequency and gate voltage varied from -60 V to +60 V. The ABA’s signal is relatively flat in the whole mapping range, while ABC’s signal has two prominent peaks within a pretty small frequency range from 1582 cm-1 to 1585 cm-1 at high gate voltage (on both positive and negative sides). Note that the minor asymmetry of the near-field signals between the positive and negative gating sides probably originates from the different gate efficiencies for the electron and hole doping, which typically occurs when gating a device in the atmosphere due to the absorption of molecules. We further plotted the ratio ABC/ABA in figure 4(c), which indicates the contrast between ABC and ABA domains at certain measurement conditions. The most distinct contrast (ABC/ABA 2.95) happens in 1583 cm-1 at Vg = +60 V and could be even larger if the gate voltage is greater. Figure 4(c) could serve as a reference for readers using SNOM to identify ABA and ABC trilayer graphene.
Figure
4.
The near-field IR behaviors of ABA and ABC trilayers against both gate voltage and excitation frequency. (a)-(c) 2D maps of ABA, ABC signals, and the ratio ABC/ABA against both gate voltage and photon frequency. (a) and (b) are in the same colorbar range to make a comparison. In (c), prominent contrast between ABC and ABA can be observed within a small frequency range of 1582-1585 cm-1 at high Vg, where the ratio ABC/ABA of the near-field signal can reach 3. (d) The horizontal linecuts of (b), showing ABC trilayer’s behaviors against gate voltage Vg at a few representative excitations. Significant changes in the line shape can be observed when frequency varies. (e) The vertical linecuts of (b), displaying the ABC graphene’s near-field IR spectra under a few representative gate voltages. The red dashed lines denote the peak position, and the vertical blue dashed line denotes the approximate location of the graphene intrinsic phonon frequency, from which a gradually enhanced red shift of graphene phonon against increased gating can be seen.
To better present how the ABC trilayer’s near-field signal varies with incident frequency and gate voltage, we plotted in figures 4(d) and (e) a series of horizontal and vertical linecuts of figure 4(b). The horizontal linecuts in figure 4(d) show gate dependence of the IR signal at a few representative excitation frequencies, from which significant changes in the line shape can be seen. Specifically, when the excitation photon is away from the Fano resonance frequency, the near-field IR signal of ABC becomes slightly weaker when the Fermi level is gated far away from the Dirac point, since dominant electron transitions can be blocked at a high Fermi level. In contrast, when the excitation photon energy goes into the Fano resonance region, it continuously evolves to be a different shape with the intensity that becomes significantly larger when the Fermi level is pulled away from the Dirac point, and becomes the most violent in 1583 cm-1 for this sample due to the contribution from the dominant electron-dressed phonon transitions. The vertical linecuts of the ABC trilayer in figure 4(e) show the gradually enhanced intensity as well as the redshift of the Fano peak when bias becomes higher, in accordance with the previous far-field IR study [40]. On the other hand, the ABA behaves much more inertly against either gate voltage or excitation frequency (more details in section 4 of supporting materials).
3.
Discussions
The distinct near-field IR behaviors between the ABA and ABC trilayers presented above can be used to find the ABC graphene at nano-scale for studying the correlated physics. We compare our near-field method with the traditional Raman spectra in figure S7. The locations of ABC domains find by our SNOM method are in accordance with that judged by Raman, while the SNOM method has a much higher resolution which is only constrained by the AFM tip radius. Specifically, for those trilayer samples with at least one large domain of micrometer size, an alternative way to identify the ABC from ABA graphene at nanoscale resolution is to combine the SNOM mapping with a single-point Raman spectroscopy. One can first use the SNOM to map out the domains and classify them into two groups by their IR contrast, and then use the single-point Raman spectrum taken from the micrometer-sized domain to assign the ABA/ABC stacking order to all the domains based on the shape and linewidth of Raman peaks.
In our experiments, the gate voltage is applied to draw a full picture of the near-field response. However, to observe the Fano resonance and identify the ABC graphene, a gate voltage is normally not required due to the commonly existing initial doping, as illustrated in figure S8. Actually, the samples that are initially charged are in the majority when placed in the atmosphere. This enables to distinguish ABC graphene without the trouble to make a metal gate electrode.
4.
Conclusions
To summarize, we studied the near-field IR behaviors of both ABC and ABA trilayers in a wide frequency range covering the graphene optical phonon frequency under different gate voltages. The ABC trilayer displays a prominent Fano resonance feature especially when subjected to gate voltages, while the ABA trilayer is largely featureless. The ratio of the near-field IR signal ABC/ABA is typically smaller than one in most excitation frequencies, but it can be reversed near the graphene phonon frequency (1585 cm-1) due to the Fano resonance. Near this frequency range, the near-field signal of ABC can be much stronger than that of ABA (ABC/ABA3) under a high gate voltage. The big diversities in ABA and ABC’s behaviors against both gate voltage and incident laser frequency can serve as a sensitive and reliable method to distinguish different trilayer stacking domains with no ambiguity. Considering the high spatial resolution and the high measurement efficiency of SNOM working in ambient conditions, the reported near-field IR method provides quick identification of ABC-stacked graphene with nanometer-scale resolution, which will definitely lower the difficulties in making ABC trilayer devices.
5.
Methods
5.1
Sample preparations
The graphene samples used in this experiment were mechanically exfoliated from bulk graphite onto the silicon substrate with an oxide layer of about 285 nm. The graphene layer numbers and domains were first selected roughly by optical microscopy and SNOM system with a CO2 laser (10.6 m), and then confirmed by a homemade Raman spectroscopy system. The selected samples were then made to be the standard field effect transistors by Electron Beam Lithography System and Electron-beam Evaporation System, and the top metal gate is made of 3 nm Cr and 30 nm Au. Then the samples were annealed at 180 C with 30 sccm H2 plasma of 35 W power for one hour to remove the organic residues on the surface. Finally, we used AFM tips at lift mode to remove the remained contaminations on the graphene surface. The force of the AFM tip was set at a tiny value to avoid changing the domain shape.
5.2
Near-field infrared imaging
The main parts of the experiment were carried out with a homemade SNOM system built from a commercial AFM system and a frequency variable QCL was used as the light source. The near-field signals were extracted using a lock-in amplifier with the frequency at three times of the AFM tip’s oscillation frequency. The near-field signal of the third harmonic is purer than the first and second harmony signals. For this reason, it is more trustable, especially for the frequency range where the graphene optical conductivity is weaker.
Acknowledgments
This work is supported by the Open Research Fund of Songshan Lake Materials Laboratory (No. 2021SLABFK07), the National Key R&D Program of China (No. 2021YFA1202902), and the National Natural Science Foundation of China (No. 12074244). Z S acknowledges support from SJTU (21X010200846), and additional support from a Shanghai talent program.
Author contributions
Z S and P S conceived this project. P S, A D prepared the samples. P S, J C, B L, A D, S L performed the near-field infrared measurements. P S, J C, A D, Z Z, S M carried out Raman spectra. X Z carried out simulations for electronic bands of trilayer graphene. P S, X Z, J C, A D, B L, Z Z, S L, S M, B W, and Z S analyzed the data. P S and Z S wrote the paper with inputs from all authors. All authors discussed the results and commented on the manuscript.
Author to whom any correspondence should be addressed.
Wang H, Gao J H, Zhang F C 2013 Flat band electrons and interactions in rhombohedral trilayer graphene Phys. Rev. B87 155116 DOI: 10.1103/PhysRevB.87.155116
[2]
Shin J, Chittari B L, Jung J 2021 Stacking and gate-tunable topological flat bands, gaps, and anisotropic strip patterns in twisted trilayer graphene Phys. Rev. B104 045413 DOI: 10.1103/PhysRevB.104.045413
[3]
Ma Z, Li S, Zheng Y W, Xiao M M, Jiang H, Gao J H, Xie X C 2021 Topological flat bands in twisted trilayer graphene Sci. Bull.66 18-22 DOI: 10.1016/j.scib.2020.10.004
[4]
Chittari B L, Chen G R, Zhang Y B, Wang F, Jung J 2019 Gate-tunable topological flat bands in trilayer graphene boron-nitride moire superlattices Phys. Rev. Lett.122 016401 DOI: 10.1103/PhysRevLett.122.016401
[5]
Akhanda M S, Rezaei S E, Esfarjani K, Krylyuk S, Davydov V, Zebarjadi M 2021 Flat bands, electron interactions, and magnetic order in magic-angle mono-trilayer graphene Phys. Rev. Mater.5 084008 DOI: 10.1103/PhysRevMaterials.5.084008
[6]
Lui C H, Li Z Q, Mak K F, Cappelluti E, Heinz T F 2011 Observation of an electrically tunable band gap in trilayer graphene Nat. Phys.7 944-7 DOI: 10.1038/nphys2102
[7]
Trohalaki S 2011 Electrically tunable bandgap observed in ABC-trilayer graphene MRS Bull.36 951-2 DOI: 10.1557/mrs.2011.313
[8]
Zhang L Y, Zhang Y, Camacho J, Khodas M, Zaliznyak I 2011 The experimental observation of quantum Hall effect of l=3 chiral quasiparticles in trilayer graphene Nat. Phys.7 953-7 DOI: 10.1038/nphys2104
[9]
Yin L J, Yang L Z, Zhang L, Wu Q L, Fu X S, Tong L H, Yang G, Tian Y, Zhang L J, Qin Z H 2020 Imaging of nearly flat band induced atomic-scale negative differential conductivity in ABC-stacked trilayer graphene Phys. Rev. B102 241403 DOI: 10.1103/PhysRevB.102.241403
[10]
Li Y W, et al 2022 Observation of coexisting dirac bands and moire flat bands in magic-angle twisted trilayer graphene Adv. Mater.34 2205996 DOI: 10.1002/adma.202205996
[11]
Hao Y F, Wang Y Y, Wang L, Ni Z H, Wang Z Q, Wang R, Koo C K, Shen Z X, Thong J T L 2010 Probing layer number and stacking order of few-layer graphene by Raman spectroscopy Small6 195-200 DOI: 10.1002/smll.200901173
[12]
Cong C X, Yu T, Sato K, Shang J Z, Saito R, Dresselhaus G F, Dresselhaus M S 2011 Raman characterization of ABA- and ABC-stacked trilayer graphene Acs Nano5 8760-8 DOI: 10.1021/nn203472f
[13]
Lui C H, Li Z Q, Chen Z Y, Klimov P V, Brus L E, Heinz T F 2011 Imaging stacking order in few-layer graphene Nano Lett.11 164-9 DOI: 10.1021/nl1032827
[14]
Lui C H, Malard L M, Kim S, Lantz G, Laverge F E, Saito R, Heinz T F 2012 Observation of layer-breathing mode vibrations in few-layer graphene through combination Raman scattering Nano Lett.12 5539-44 DOI: 10.1021/nl302450s
[15]
Malard L M, Pimenta M A, Dresselhaus G, Dresselhaus M S 2009 Raman spectroscopy in graphene Phys. Rep.473 51-87 DOI: 10.1016/j.physrep.2009.02.003
[16]
Mak K F, Shan J, Heinz T F 2010 Electronic structure of few-layer graphene: experimental demonstration of strong dependence on stacking sequence Phys. Rev. Lett.104 176404 DOI: 10.1103/PhysRevLett.104.176404
[17]
Warner J H, Mukai M, Kirkland A I 2012 Atomic structure of ABC rhombohedral stacked trilayer graphene ACS Nano6 5680-6 DOI: 10.1021/nn3017926
[18]
Yin L J, Wang W X, Zhang Y, Ou Y Y, Zhang H T, Shen C Y, He L 2017 Observation of chirality transition of quasiparticles at stacking solitons in trilayer graphene Phys. Rev. B95 081402 DOI: 10.1103/PhysRevB.95.081402
[19]
Wu H J, et al 2021 Direct visualization and manipulation of stacking orders in few-layer graphene by dynamic atomic force microscopy J. Phys. Chem. Lett.12 7328-34 DOI: 10.1021/acs.jpclett.1c01579
[20]
Zhang S, et al 2022 Domino-like stacking order switching in twisted monolayer-multilayer graphene Nat. Mater.21 621 DOI: 10.1038/s41563-022-01232-2
[21]
Chen J N, et al 2012 Optical nano-imaging of gate-tunable graphene plasmons Nature487 77-81 DOI: 10.1038/nature11254
[22]
Fei Z, et al 2012 Gate-tuning of graphene plasmons revealed by infrared nano-imaging Nature487 82-85 DOI: 10.1038/nature11253
[23]
Lundeberg M B, Gao Y D, Woessner A, Tan C, Alonso-Gonzalez P, Watanabe K, Taniguchi T, Hone J, Hillenbrand R, Koppens F H L 2017 Thermoelectric detection and imaging of propagating graphene plasmons Nat. Mater.16 204-7 DOI: 10.1038/nmat4755
[24]
Jiang L L, Wang S, Shi Z W, Jin C H, Utama M I B, Zhao S H, Shen Y R, Gao H J, Zhang G Y, Wang F 2018 Manipulation of domain-wall solitons in bi- and trilayer graphene Nat. Nanotechnol.13 204 DOI: 10.1038/s41565-017-0042-6
[25]
Nikitin A Y, et al 2016 Real-space mapping of tailored sheet and edge plasmons in graphene nanoresonators Nat. Photon.10 239 DOI: 10.1038/nphoton.2016.44
[26]
Woessner A, et al 2015 Highly confined low-loss plasmons in graphene-boron nitride heterostructures Nat. Mater.14 421-5 DOI: 10.1038/nmat4169
[27]
Shi Z W, Hong X P, Bechtel H A, Zeng B, Martin M C, Watanabe K, Taniguchi T, Shen Y R, Wang F 2015 Observation of a Luttinger-liquid plasmon in metallic single-walled carbon nanotubes Nat. Photon.9 515-9 DOI: 10.1038/nphoton.2015.123
[28]
Zhao S H, et al 2018 Correlation of electron tunneling and plasmon propagation in a Luttinger liquid Phys. Rev. Lett.121 047702 DOI: 10.1103/PhysRevLett.121.047702
[29]
Ju L, et al 2015 Topological valley transport at bilayer graphene domain walls Nature520 650-U356 DOI: 10.1038/nature14364
[30]
Li P N, Dolado I, Alfaro-Mozaz F J, Casanova F, Hueso L E, Liu S, Edgar J H, Nikitin A Y, Velez S, Hillenbrand R 2018 Infrared hyperbolic metasurface based on nanostructured van der Waals materials Science359 892 DOI: 10.1126/science.aaq1704
[31]
Qu Y P, et al 2022 Tunable planar focusing based on hyperbolic phonon polaritons in alpha-MoO3Adv. Mater.34 2105590 DOI: 10.1002/adma.202105590
[32]
Wu Y J, et al 2020 Chemical switching of low-loss phonon polaritons in alpha-MoO3 by hydrogen intercalation Nat. Commun.11 2646 DOI: 10.1038/s41467-020-16459-3
[33]
Joshi T, et al 2017 Coupled one-dimensional plasmons and two-dimensional phonon polaritons in hybrid silver nanowire/silicon carbide structures Nano Lett.17 3662-7 DOI: 10.1021/acs.nanolett.7b00845
[34]
Zhang Y R, et al 2020 Tunable cherenkov radiation of phonon polaritons in silver nanowire/hexagonal boron nitride heterostructures Nano Lett.20 2770-7 DOI: 10.1021/acs.nanolett.0c00419
[35]
Du J, et al 2021 Fano resonance enabled infrared nano-imaging of local strain in bilayer graphene Chin. Phys. Lett.38 056301 DOI: 10.1088/0256-307X/38/5/056301
[36]
Lyu B S, et al 2019 Phonon polariton-assisted infrared nanoimaging of local strain in hexagonal boron nitride Nano Lett.19 1982-9 DOI: 10.1021/acs.nanolett.8b05166
[37]
Huber A J, Ziegler A, Kock T, Hillenbrand R 2009 Infrared nanoscopy of strained semiconductors Nat. Nanotechnol.4 153-7 DOI: 10.1038/nnano.2008.399
[38]
Qian J, Luan Y L, Kim M S, Ho K M, Shi Y, Wang C Z, Li Y, Fei Z 2021 Nonequilibrium phonon tuning and mapping in few-layer graphene with infrared nanoscopy Phys. Rev. B103 L201407 DOI: 10.1103/PhysRevB.103.L201407
[39]
Tang T T, et al 2010 A tunable phonon-exciton Fano system in bilayer graphene Nat. Nanotechnol.5 32-36 DOI: 10.1038/nnano.2009.334
[40]
Lui C H, Cappelluti E, Li Z Q, Heinz T F 2013 Tunable infrared phonon anomalies in trilayer graphene Phys. Rev. Lett.110 185504 DOI: 10.1103/PhysRevLett.110.185504
[41]
Aoki M, Amawashi H 2007 Dependence of band structures on stacking and field in layered graphene Solid State Commun.142 123-7 DOI: 10.1016/j.ssc.2007.02.013
[42]
Kumar S B, Guo J 2011 Multilayer graphene under vertical electric field Appl. Phys. Lett.98 222101 DOI: 10.1063/1.3595335
[43]
Mei H Y, Xu W, Wang C, Yuan H F, Zhang C, Ding L, Zhang J, Deng C, Wang Y F, Peeters F M 2018 Terahertz magneto-optical properties of bi- and tri-layer graphene J. Phys.: Condens. Matter30 175701 DOI: 10.1088/1361-648X/aab81d
[44]
Dressel M, Scheffler M 2006 Verifying the drude response Ann. Phys.15 535-44 DOI: 10.1002/andp.200651807-810
[45]
Fano U 1961 Effects of configuration interaction on intensities and phase shifts Phys. Rev.124 1866 DOI: 10.1103/PhysRev.124.1866
Jäckering, L., Wirth, K.G., Conrads, L. et al. Super-Resolution Imaging of Nanoscale Inhomogeneities in hBN-Covered and Encapsulated Few-Layer Graphene. Advanced Science, 2025.
DOI:10.1002/advs.202409039
2.
Luo, W., Fan, J., Kuzmenko, A.B. et al. Topological valley plasmons in twisted monolayer-bilayer graphene moiré superlattices. Physical Review B, 2024, 110(3): 035408.
DOI:10.1103/PhysRevB.110.035408
3.
Fox, C., Mao, Y., Zhang, X. et al. Stacking Order Engineering of Two-Dimensional Materials and Device Applications. Chemical Reviews, 2024, 124(4): 1862-1898.
DOI:10.1021/acs.chemrev.3c00618
DownLoad:
