A Study on Graphene—metal Contact

The contact resistance between graphene and metal electrodes is crucial for the achievement of high-performance graphene devices. In this study, we review our recent study on the graphene–metal contact characteristics from the following viewpoints: (1) metal preparation method; (2) asymmetric conductance; (3) annealing effect; (4) interfaces impact.


Introduction
Graphene is a single atomic sheet of carbon atoms in a honeycomb lattice, where the carbon-carbon bonds in the plane are sp 2 hybridized.Ever since it was firstly isolated from micromechanical cleavage of bulky graphite flakes in 2004 [1], graphene has been widely investigated owing to its unique and attractive electrical, physical and chemical properties [1][2][3][4][5][6][7][8][9][10][11].Graphene is an energy gapless material, in which the conduction and valence bands meet at the Dirac point.Its linear electronic dispersion relation results in a zero effective mass of the carrier with a high Fermi velocity.This material thus exhibits an extraordinarily high carrier mobility of more than 200,000 cm 2 V −1 s −1 [12].

OPEN ACCESS
It is considered to be a potential channel material for field-effect transistors (FETs) application.In addition, graphene has been applied to various practical device applications, such as high frequency devices [9], gas sensors [5,13], flexible electronics [14,15], and photonics [16,17].In spite of the very high mobility in the graphene channel, however, the contact resistance between graphene and metal electrodes is crucial to achieving a high performance from the graphene [18][19][20], especially the high on-state current.It has been reported that the control of contact properties is more important than the intrinsic channel mobility, as otherwise the merit of high mobility from the graphene will be diminished significantly [18].In this paper, we provide a study on the characteristics of graphene-metal contact from the perspectives of: (1) metal preparation method; (2) asymmetric conductance; (3) annealing effect; and (4) interface impact.These findings offer insightful information to achieve high performance graphene devices via process optimization.

Impact of Metal Preparation
When measuring contact resistance, the contribution from the graphene sheet resistance should be minimized or eliminated.There are several methodologies to extract the graphene-metal contact resistance, including the transfer length [21,22], four-probe/two-probe [18,19,23], and residual resistance methods [20,24].We use the last two methods to extract the contact resistance.Ti has been used as metal electrodes in carbon-based devices because of its excellent adhesive capability in SiO 2 substrate or other insulators [25].In this section, we take Ti/graphene as an example to study the impact of different metal deposition methods on R c .
The total resistance between two contacts is the sum of the semiconductor resistance, the contact resistance, and the metal resistance, where the metal resistance can be neglected owing to good control of metal electrode growth.By subtracting the graphene resistance from the total resistance, one obtains the total contact resistance.The contact resistance for each contact is thus obtained as half of the total resistance [18,26,27]: where R c is the contact resistance between graphene and metal, which is assumed to the identical for the two contacts in back-gate graphene FETs; R total is the total resistance measured between source and drain; R g is the graphene resistance in channel, derived by R g = (V 2 -V 3 )/I, schematically shown in inset in Figure 1b.
Figure 1a shows the I -V characteristics of SLG (single-layer graphene)/Ti devices, where Ti is grown by electron-beam evaporation (EBM) and sputter processes, respectively.In our experiments, graphene is produced by mechanical exfoliation from bulk graphite.Subsequently, graphene is transferred onto degenerated silicon substrate with a thermally grown 90 nm SiO 2 , serving as a conventional back-gate dielectric.R total can be extracted from the linear I -V curves.The optical image of measured device is shown in inset in Figure 1(a).The linear I -V characteristic indicates that the Ti/SLG contact is ohmic.Note that there are also ohmic contacts between graphene and other metals in our samples, e.g., Al, Ag, Pd, Ni, Au.However, there is a Schottky junction between graphene nano-ribbon (GNR) and metals [28].The nonlinear I -V characteristic has been reported in GNR-Au, GNR-Al contacts theoretically and experimentally [29].This Schottky effect in GNR-metal contacts stems from the creation of energy band gap in GNR, where the energy band gap increases with reduction of the GRN width [8].The graphene resistance R g is measured at the force current from 1 µA to 1 mA, shown in Figure 1b.It was found that R g is almost the same at different force currents.By combining the results from two-probe (Figure 1a) and four-probe measurement (Figure 1b), R c is quantitatively addressed.The R c of EBM sample is observed to have a lower R c of 0.83 kΩ as compared to sputter one of 4.2 kΩ.Note that the poor on/off ratio is observed in our graphene FETs, owing to its gapless characteristic.Note that in the four-probe measurement, the back-gate electrode of graphene FETs is floating, and R g is almost a constant.While in residual resistance measurement, the R g is changeable according to the modulation of carrier density in graphene channel.The R c then can be determined as follows [20]: where V bg is back-gate voltage; R g is V bg dependent, owing to the modification of carrier density in graphene channel; The residual resistance R resi , which is the saturated resistance at high V bg , is equal to 2R c due to negligible R g at large negative V bg , indicating the total resistance is mainly from the contribution of contact resistance [20], schematically illustrated in Figure 2a,b.Using residual resistance measurement method, one may obtain that R c is 0.78 kΩ and 4.1 kΩ for EBM and sputter Ti/SLG devices, respectively, as shown in Figure 3.It was found that R c extracted from the two-probe/four-probe and residual resistance measurement methods is in excellent agreement with each other.The I ds − V ds characteristics at various back-gate voltages also exhibit the linear behavior, shown in the inset in Figure 2. The measured contact resistivity ρ c of SLG, DLG (double-layer graphene), MLG (multi-layer graphene)/Ti devices, prepared by EBM and sputter processes, is summarized in Figure 3.The standard error is calculated from eight samples for each group.For all the devices fabricated by the EBM process, the ρ c does not exhibit strong dependence on the number of graphene layers.It is consistent with the early work, where R c of SLG, DLG, and TLG FETs (also EBM samples) is insensitive to layer thickness [30], while for the devices prepared by sputter process, the ρ c exhibits layer dependence and increases with decreasing the number of layers.It is worth noting that the ρ c of Ti/SLG and Ti/DLG from the former is significantly smaller than that from latter.Note that there is a negligible difference in ρ c between the EBM and sputter samples.Thus it is proposed that the sputter Ti atoms would only effectively affect the top layers (up to two layers) of graphene in our case.It is believed that, for MLG, after the top layers create the vacancies, the metal can "penetrate" through the vacancies to contact with the bottom layer directly.Figure 4 shows the ρ c distribution of Ti/SLG devices prepared by sputter processes at various powers.It is also observed that the ρ c increases more than two orders of magnitude as sputter power enhances.It was reported that the ρ c can reach 10 9 Ω μm in case of sputtered Ti on SLG or MLG [18].It infers more carbon atoms are milled away when the sputter power increases.
The distribution of ρ c of Ti/SLG, DLG, MLG devices prepared by EBM and sputter processes, respectively.The solid symbols represent the average ρ c and the error bars come from standard error.The average ρ c of sputtered Ti/SLG and DLG devices is larger than EBM ones.Furthermore, the ρ c of the SLG device is larger than that of the DLG (MLG) devices for the sputter process.Reprinted with permission from [27].Copyright (2011) by the Electrochemical Society.

Figure 4.
The distribution of the ρ c of Ti/SLG devices prepared by sputter processes at different powers (power 1 < power 2).The ρ c increases more than two orders of magnitude when the sputtering power increases.
Figure 5 shows the Raman spectra of SLG, DLG, TLG (triple-layer graphene) with and without metal deposition, where they are prepared by (a) EBM; and (b) sputter processes, respectively.Raman spectroscopy has been used to physically probe the electronic structure of graphite and graphene without damaging the sample [31][32][33][34].Since our graphene sheet is prepared by mechanical exfoliation from bulky graphite, there is no D band in Raman spectra of graphene before metal deposition.There are only two main characteristic peaks of G and 2D bands before metal growth, where the peak position for G and 2D bands are around 1580 and 2670 cm −1 , respectively.For EBM samples, there is negligible D band, around at 1350 cm −1 , after Ti deposition by EBM process.It is worth noting that D bands are formed after receiving sputtered Ti for all SLG/DLG/TLG samples, significantly different from the ones after EBM process.The D band is caused by disordered or defected structure of graphene sheet [31,35,36].The intensity ratio of D band to G band (I D /I G ) is usually used to estimate the amount of defects in carbon materials.Accordingly, compared to sputtered Ti/graphene contact, a Ti/graphene contact prepared by EBM corresponds to very low defects or carbon vacancies in graphene.It has been reported that the formation energy of carbon vacancies in graphite is around 7.4 eV [37,38].Thus, the ion energy in the sputter process should be larger than that to create the carbon vacancies in graphene.Generally, it is believed that sputtered Ti atoms possess larger kinetic energy compared to the EBM case, and the energy could be transferred to the graphene layer, resulting in the removal of carbon atoms from the graphene lattice and creating of the carbon vacancies.In addition to the presence of the D band, the G and 2D band shifts are also observed in sputtered SLG, DLG, and MLG/Ti contacts.In the meantime, it is found that there is no noticeable band shift after Ti deposition onto graphene, where Ti is grown by EBM process, indicating no change in lattice constant in graphene and graphite, shown in Figure 5a.On the contrary, the 2D band and G band have an obvious red shift of ~30 cm −1 and ~20 cm −1 after Ti deposition by sputter process.G band shift caused by charge doping through phonon-electron coupling has been reported by Yan et al. [39].The 2D band is a second-order two-phonon process and exhibits a strong frequency dependence on the excitation laser energy [32].However, the additional band shift will be narrow at ~10 cm −1 in case of (c) graphene charge doping [39,40].Therefore, this band shift cannot be attributed only to the charge doping effect.Tensile stress is also reported to cause the red shift in graphene; the red shift expands with increasing tensile stress [41,42].As the 2D band originates from the two-phonon double-resonance process, it is closely related to the band structure of graphene layers.It is thus believed that the sputter process can lead to tensile strain in graphene underneath Ti owing to defects formation, thereby enlarging the lattice constant of graphene.The red shift of the 2D and G bands can thus be understood as the tensile strain weakening the bond and thus lowering the vibration frequency due to the elongation of C-C bonds [42].
A schematic model has been proposed to explain the difference in graphene-metal contacts between large and small grains in contact metal [22].In their model, the large grains and rough surface of contact metal play an important role in the contact area and small contact area results in large contact resistance.Then, contact resistance can be affected by the grain size and the uniformity of the contact metal films.For sputtered SLG and DLG/Ti devices, there is also the possibility of a smaller contact area, as parts of carbon atoms are milled away from the pristine graphene structure, accordingly increasing the contact resistance.It must be pointed out that the defected carbon or carbon vacancies not only occur in the graphene underneath the metal, but also at the adjunct region of graphene channel and graphene-metal contact.It was reported that the defected graphene will break the symmetry of regular hexagonal C-C bond structures, thus resulting in the intervalley electron scattering [43,44]; accordingly, the mobility will degrade compared to the pristine one.Note that there is possible carbide formation after Ti deposition onto the graphene, similar to surface-carbide formation in additional Ni deposition on graphene [45].However, the number of Ti-C bonds is very low and their influence on electron transport is negligible [46].It was also observed that the significant increase in D band sputtered Al/SLG, Al/DLG contacts, asshown in Figure 6.The I D /I G is larger than 2. The 2D band even disappears as the sample is subjected to high sputter power, indicating that the pristine structure of graphene has been significantly destroyed.Note that the reduced R c is observed in high vacuum deposition condition, suggesting deposition pressure also has a significant influence on the quality of R c [30].

Asymmetry Conductance
Because graphene is an energy gapless material, it is very similar to the metal-metal contact when graphene is brought into contact with metal.With metal-metal contact, a small redistribution of electrons can screen the potential difference due to the large carrier density, where the potential difference originates from the work function difference in two metals.The screening length λ is expressed by ) ( 41 is the DOS at the E F in graphene.The λ is less than 1 nm for metal-metal contact.However, for graphene-metal contact, DOS is very small compared to metal, especially at the Dirac point, and λ is much larger than that of metal.In addition, graphene itself is sensitive to external perturbations owing to its all-surface and zero volume nature [48].Charge transfer between graphene and metal may dope the underlying graphene into n-type or p-type, depending on work function differences.Therefore, the metal doping effect must be considered in transfer characteristics of graphene FETs.In Figure 2, it is shown that p-type and n-type branches show a moderate asymmetry and the conductance in the n-type branch is always smaller than that in the p branch.It is consistent with the previous report [18].This asymmetric conductance behavior was first investigated by comparing devices with invasive electrodes crossing the whole graphene channel width and external electrodes connected to the side of the graphene channel [23].The graphene FETs with invasive electrode show the asymmetry of conductance in n-type and p-type branches; however, the asymmetry almost never occurs in devices with external electrodes [18,23].Giovannetti and Khomyakov et al. [49,50] applied density functional theory (DFT) to study how graphene is doped by various metals, including Al, Ag, Cu, Au, Pt.The calculated Fermi energy shift (∆E F ) with respect to the conical point can increase by decreasing the distance between graphene and metal.Xia et al., by using scanning photocurrent microscopy (SPCM), illustrated that the charge doping occurs not only underneath the metal, but also extends hundreds of nanometers into adjacent regions in the graphene channel [51,52].They separated a graphene FET into three segments: (1) the metal-controlled graphene; (2) the transition region affected by both the metal and the back-gate; and (3) the bulk graphene region controlled only by the back-gate.With the modulation of back-gate voltage, the width of transition region (segment 2) is accordingly modified.In their experiments, the detected photocurrent comes from the p-n junction in the graphene sheet, since the photocurrent is proportional to the potential gradient at the excitation position.It exhibits the different polarities along the graphene channel.Lee et al. [53] observed the photocurrent in graphene FETs with the help of SPCM, resulting from p-n junction formation in the graphene sheet.It was also reported that chemical doping can introduce conductance asymmetry.It is caused by a combination of the neutrality point misalignment at the electrode-channel interface and the non-constant DOS of the graphene electrodes [54].Previously, we implemented scanning Kelvin probe microscopy (SKPM) to experimentally address the work function difference between graphene channel and metal electrodes [55].The graphene underneath the metal is assumed to be pinned, and the constant potential   can be described as follows [52]: where  is the Plank constant divided by  2 ; flat bg V is the flat band voltage; ν F is Fermi velocity of 0.8 × 10 6 m/s; α is 2.24 × 10 11 cm −2 v −1 , which is estimated from a simple capacitor model [56].The results show that the effective constant potential   for Ti/SLG, Al/SLG, Ag/SLG Pd/SLG is 0.1 eV, −0.22 eV, 0.08 eV, 0.45 eV, respectively.For Ti/SLG, our calculated   is consistent with Xia's report of ~0.12 eV [51] or 0.1 eV [52], Figure 7a-c shows the R total − V bg plots of graphene FETs with different metal electrodes Ag, Pd, Al, respectively.They exhibit an asymmetric behavior in hole and electron transportation branches, including the graphene FETs with Ti electrode, shown in Figure 2.
The numerical simulation has shown that the Fermi level of graphene underneath the metal is shifted and may result in asymmetric transfer characteristics [57], which is similar to the aforementioned.The odd resistance is used to characterize the asymmetry in graphene FETs with various metal electrodes.
The R odd is defined as follows [23,52]: where ∆V bg = V bg − V Dirac ; The R odd is positive for Ti, Ag, Au, Pd, and negative for a Al/SLG FET, indicating hole branch dominating transportation for the Ti, Ag, Au, Pd cases, while electron branch dominates the transportation for Al.It is observed that the Dirac point of measured graphene locates at the region of positive voltage, indicating that our fabrication process introduces a p-type doping in the graphene channel.For Ag (Ti, Au, Pd)/SLG FETs, the p-p-p (p-n-p) junction forms when ).The excess resistance in the electron transportation region for Ag (Ti, Au, Pd)/SLG FETs is the result of the p-n-p junction formation, whereas the n-p-n junction is believed to be the cause of high resistance in the hole transportation region for Al/SLG FETs.All R total vs. (V bg − V Dirac ) curves show the asymmetry in hole and electron transportation regions.R odd is positive for Ag, Ti, Au, Pd-graphene devices, and negative for Al/SLG FETs.

Annealing Effect on Contact Junction
In order to improve the performance of graphene-metal contact, the forming gas annealing (FGA) were performed to examine the R c of the sputtered Ti/SLG junction.Figure 8a shows the transfer characteristics of a SLG field-effect transistor before and after forming gas annealing (FGA), whereby the process is undergoing with N 2 :H 2 mixing at 425 °C for half hours.After FGA treatment, the R c could not be reduced as expected.The Raman spectra of the Ti/SLG contact are shown in Figure 8b.From Figure 8b, it is observed that the presence of the D band during the sputter process does not disappear, but rather expands.Concurrently, the FWHM of the G band also expands with the D band.The broadening of the D peak is correlated to the distribution of cluster with different orders and dimensions [58], thus suggesting that there are C-H sp 3 bonds [59] or carbon defects in the graphene layer underneath the metal after FGA treatment, which is very similar to amorphous carbon [58].From Figure 8a again, it is interesting to observe that the Dirac point is shifted to the right direction, or positive polarity.The right shift of V Dirac , and the ∆V Dirac is more than 50 V, indicating a strong hole doping effect after FGA.The doped charge density, which is induced by FGA treatment, can be estimated by the following equation [60,61], where C 0 is the capacitance between channel and back-gate per unit area; ∆V Dirac is the change in the voltage at the Dirac point, and q is the unit charge.For a 90-nm-thick SiO 2 , C 0 = ε 0 ε/d = 38.3nF/cm 2 , where ε 0 , ε is the permeability of vacuum and relative permittivity of the gate dielectric (3.9 for SiO 2 ), respectively.Accordingly, the FGA treatment-induced hole concentration is around ~1 × 10 13 cm −2 .To further confirm this, Raman measurement in graphene channel is performed, as well.The absence of the D band in the channel also infers that the D band only stems from the high energetic sputter metal in our sputter sample.It was found that the G band had a right shift and the decrease is observed in its FWHM.The hole doping can be reflected in the stiffening and sharpening of the G band [39,62], as shown in Figure 8c.Raman examination of the blue shift reveals a decrease in full width at half maximum (FWHM) in the G band on the SLG channel, which indicates that the hole doping increases after FGA, leading to the Dirac point shifting toward high positive voltage [63].A low temperature rapid thermal annealing (RTA) is proposed to eliminate process-induced resistance [64], while RTA at high temperature will degrade the performance of the device.

Interface Impact
There are two interfaces in a graphene-metal contact system: graphen/subtrate and graphene/metal interfaces.Graphene always physically lies on the surface of SiO 2 /Si substrate, or another oxide/Si substrate, thus the first interface must necessarily be taken into account.Previously, the properties of graphene between SiO 2 , Al 2 O 3 , and HfO 2 are studied by using Raman and atomic force microscopy (AFM) [65].In addition to strong adhesion between graphene and dielectric, a higher hole concentration also occurs as graphene lies in a higher dielectric constant oxide.Figure 9a,b present the schematics illustration of SiO 2 surface treated by HF dipping and the re-oxidation process, respectively [66].The difference between HF dipping and the re-oxidation process is that the SiO 2 surface is terminated by OH bonds for the former, and the latter is an O-terminated SiO 2 surface.The former is hydrophilic and dopes graphene into the p-type, while the latter is hydrophobic and dopes graphene into the n-type, which is verified in the Raman G band position [66,67].It is also reported by Casiraghi et al. [68] that the G band position on substrates are within the range of fluctuation (1580-1588 cm −1 ) by unintentional electron or hole doping effect for more than 40 graphene samples on SiO 2 /Si substrate.Watanabe's observation shows that R c is independent of work function of contact metal such as Ti, Ag, Co, Cr, Fe, Ni, and Pd [22].Robinson et al. [69] also examined the effect of work function difference on ρ c with Al, Ti, Cu, Pd, Ni, Pt, where the work function difference with graphene varies from −0.2 eV to 1.2 eV.There is little difference in ρ c as the difference in metal/graphene work function varies, indicating that the device fabrication process heavily influences one's ability to "dope" graphene [69].This is the possible reason for contact resistance results widely varying among reported experiments, even for the same metal.The second interface is also important for understanding graphene-metal contacts.Antonio et al. [70] ascribed double dip in I ds − V bg to charge transfer between the graphene and the metal electrodes.Nouchi et al. [71] have also reported that anomalously distorted transportation originated from the partially formed oxide layer in Ni/graphene contact.The previous work assumed that carrier density of graphene underneath metal, described as segment 1 [52], is pinned [51,52,55].It is independent from back-gate voltage, which was believed to modulate the graphene channel, rather than the graphene underneath metal.However, Raman and transfer measurements show that carrier density of graphene underneath metal can be modulated, indicating graphene contacting metal is still graphene because of a weak interaction [72].In their experiments, there exists residual photo-resistor between (PR) graphene and metal.Chen [73], Xia [22] and Song's [74] observations also demonstrated that contact resistance is V bg dependent, and also quantum capacitance is observed [74].To eliminate the side effect of residual PR, Hsu et al. [75] intentionally inserted an Al sacrificial layer during lithography.Graphene surface roughness underneath the ohmic contacts is reduced; they were able to achieve the improved contact resistance to 200-500 Ω μm.Very recently, Nagashio et al. [76] developed a photo-resistor-free method to fabricate graphene FETs with Ni electrode.The modulation of the G band of Ni/graphene contact with and without residual PR is compared.The position of G band, corresponding to the doping effect, can be modulated up to 5 cm −1 for Ni/graphene with residual PR, while the shift of the G band is limited to 1-2 cm −1 in PR-free samples.It indicates that the graphene underneath metal even without residual PR is still graphene, but the modulation is restricted, due to strong interaction from π − d coupling between graphene and Ni [50,77].In the case of Ni/graphene with PR-free process, R c is not improved compared to the case of PR.Theoretical calculation by Matsuda et al. [78] shows R c will reduce in an "end-contacted" graphene-metal interface, where carbon pπ orbitals, as well as pσ orbitals, play important roles in cohesion.Therefore, a sandwich contact structure (metal-graphene-metal) [79,80] and graphitic contact [76] has been suggested to reduce R c .

Conclusions
We discuss the characteristics of graphene-metal contact from the viewpoints of metal preparation, asymmetric transportation, annealing effect, and interface impact.These findings provide insightful information into process optimization not only for graphene devices, but also for other layer-material devices such as MoS 2 , WS 2 .

Figure 1 .
Figure 1.(a) The I 2p -V 2p plots of Ti/SLG devices at for EBM and sputter processes, respectively.The optical image of measured device is shown in inset of Figure 1a; (b) The graphene resistance R g (R 4p ) between two probes as a function of force current.The schematic for four-probe measurement is shown in inset of Figure 1b.R c can thus be quantitatively obtained by combining the results from two-probe and four-probe measurements.Reprinted with permission from [27].Copyright (2011) by the Electrochemical Society.

Figure 2 .
Figure 2. (a) The schematic of residual resistance measurement; (b) The equivalent circuit for graphene FETs.R g can be modulated by back-gate voltage; (c) R − V bg plots of the Ti/SLG device prepared by EBM and sputter processes, respectively.The residual resistance R resi equals to 2R c at high negative voltage.The I ds − V ds plots as a function of V bg are shown in the inset in Figure 2. R c is larger for the sputtered Ti/SLG device than EBM one.

Figure 5 .
Figure 5. Raman spectra of 1-3 layers graphene with and without Ti (~9 nm) deposition by (a) EBM; and (b) sputter process, respectively.The noticeable D band in sputtered Ti/grapheme junction, which is related to the defects (carbon vacancies) in graphene, is marked by a red dash box in Figure 5b; (c) The optical image of Raman sample.1, 2, 3, and 4 in Figure 5c represent SLG, DLG, TLG, and MLG (multi-layer graphene), respectively.Reprinted with permission from [27].Copyright (2011) by the Electrochemical Society.

Figure 6 .
Figure 6.The Raman spectra of SLG, DLG before and after sputtered Al.The optical image of measured Al/SLG, Al/DLG contacts is presented in the inset.

Figure 7 .
Figure 7. R odd − V bg plots of SLG FETs with various metal electrodes (a) Ag; (b) Pd; (c) Al, respectively.Insets show the Fermi level (red dash) with respect to the Dirac point at the hole and electron transportation regions; (d) R odd vs. (V bg − V Dirac ) plots with various electrodes.All R total vs. (V bg − V Dirac ) curves show the asymmetry in hole and electron transportation regions.R odd is positive for Ag, Ti, Au, Pd-graphene devices, and negative for Al/SLG FETs.

Figure 8 .
Figure 8.(a) R total of sputtered Ti/SLG device as a function of V bg before and after forming gas annealing (FGA).The FGA does not improve R c of Ti/SLG; (b) Raman spectra of Ti/SLG junction before and after FGA, the D band does not disappear; (c) Raman spectra of the graphene device channel before and after FGA.

Figure 9 .
Figure 9.The schematics illustration of SiO 2 surface treated by (a) HP dipping; and (b) re-oxidation process.After HF solution dipping, the SiO 2 surface will mainly be terminated by OH bonds and some hydrocarbons (e.g., −CH 3 ).While for the reoxidation-treated sample, it is an O-terminated SiO 2 surface.