The physics of adhesion of one-dimensional nano structures such as nanotubes, nano wires, and biopolymers on different substrates is of great interest for the study of biological adhesion and the development of nano electronics and nano mechanics. In this paper, we present force spectroscopy experiments of individual single wall carbon nanotube loops using a home-made interferometric atomic force microscope. Characteristic force plateaus during the peeling process allow the quantitative measurement of the adhesion energy per unit length on various substrates: graphite, mica, platinum, gold, and silicon. Moreover, using a time-frequency analysis of the deflection of the cantilever, we estimate the dynamic stiffness of the contact, providing more information on the nanotube configurations and its intrinsic mechanical properties.
I. INTRODUCTION
Since their discovery,1 carbon nanotubes (CNTs) have attracted the interest of scientists for their unique electrical,2 thermal,3 and mechanical properties4 and are foreseen as a major material in a huge range of applications in the next few decades, from materials reinforcement5 to components of nanoscale electronics6 and mechanics, for instance, nanoswitches,7 motors,8 actuators,9–11 etc. They are also widely used as a bench system to study fundamental physical phenomena on the mesoscopic scale and represent as such an archetype of nano-objects. Whatever exceptional their intrinsic properties are, their use in any application is linked to their interactions with their environment, mainly through the weak adhesive Van der Waals (VdW) force. The study of the VdW interaction between nanotubes and the rest of the world can thus help to understand the physics of polymer nanocomposites,5,10,12 setae adhesion of geckos,13 protein filament adhesion in mussels,14 nanotube-tipped atomic force microscopy (AFM) probes,15 nanoscale sensors,16 gecko-foot-mimetic dry adhesives,17 etc.
Up to now, the adhesive properties of nanotubes have been mostly probed by various smart but indirect measurements. Among those, Hertel and coworkers18,19 imaged by AFM the shape of crossed nanotubes adsorbed on a silicon substrate. The profile of the top nanotube results from a balance of the deformation energy with the surface energy lost in this configuration, allowing the estimation of an adhesion energy when mechanical properties of the nanotubes are assumed. Kis and coworkers20 performed a direct measurement using an AFM tip to pull the inner core of a telescopic multi wall CNT. Their experiment demonstrated a friction free interaction between the concentric layers and provides an estimation of the adhesion for this very specific geometry and material.
However, in these experiments, one usually accesses either intrinsic properties of the nanotube or its interaction with its environment using hypotheses on the other properties. Direct measurement of adhesive interactions can provide quantitative values of several properties in one single test and are thus an appealing method. For example, measurements of adhesion energy and Young's modulus of graphene layers on SiO2 have been conducted in single bulge test experiments.21 For nanotubes, a peeling test is a potentially powerful technique to characterize the adhesion properties of carbon nanotubes or nano wires with various substrates. A few experiments have been conducted,22–27 along with theoretical/numerical modeling.28,29 Quantitative measurements however have not been easily achieved, and the experimental data analysis relies on complex comparison with numerical simulations. We recently proposed a simpler protocol to perform peeling test,30 which allow direct and quantitative characterization of the adhesion properties of nanotubes. A single wall CNT (SWCNT) is anchored to the tip of an AFM cantilever and is simply pushed almost perpendicularly to a flat surface. The VdW interaction causes part of the nanotube to adhere to the surface when the induced bending is strong enough, and the analysis of the force curve leads to quantitative information both on the adhesion process and on the nanotube itself.
In this paper, we present an extension of this protocol to CNT loops on various materials, leading to quantitative values of the energy of interaction of the nanotube with substrates of graphite, mica, platinum, gold, and silicon. We first present the samples (CNT and substrates) and the acquired data during approach-retract cycles: force, dynamic stiffness, and compression. We then give the analyzing framework and discuss the nanotube configurations during peeling. We finally conclude, giving quantitative values of the adhesion energy of a SWCNT on the five different substrates.
II. EXPERIMENTS
A. Samples preparation
The nanotubes are grown directly32 at the tip apex of AFM probes by Chemical Vapor Deposition (CVD): the bare silicon cantilevers are fully dipped into the catalyst solution, then gently dried in a nitrogen flux before being placed in the furnace. CNTs grow everywhere on the cantilever, and around 1 every 3 cantilevers present a CNT at the tip. The parameters are tuned to grow long SWCNT (a few micrometers), leading to a high probability of having nanotube loops (two anchoring points on the tip). A few Scanning Electron Microscopy (SEM) images of our samples right after the growth are presented in Figure 1. Note however that nanotubes are often detected during AFM measurements, while they are not visible on the SEM images (the opposite is also true): the first contacts between the nanotube and the substrate have a strong influence on the nanotube shape. They can lift some CNTs from the tip, for example, or break their link to the tip. Only samples with reproducible behavior during numerous force cycles are considered in the present work.
(Top) Scanning electron micrographs of several SWCNTs grown directly on AFM tips. Long nanotubes tend to form loops on the tip, with a typical diameter around 1 μm. (Bottom) When the nanotube loop is pushed against a flat surface, part of the nanotube adheres to the surface due to Van der Waals interactions. The radius of curvature Ra at the last point of contact on the CNT with the substrate is fixed by an equilibrium between the adhesion of the part in contact and the bending of the free standing part of the nanotube. From the measurement of the AFM cantilever deflection d (using differential interferometry31) and sample position zs, the force F acting on the nanotube and its compression zc can be recorded.
(Top) Scanning electron micrographs of several SWCNTs grown directly on AFM tips. Long nanotubes tend to form loops on the tip, with a typical diameter around 1 μm. (Bottom) When the nanotube loop is pushed against a flat surface, part of the nanotube adheres to the surface due to Van der Waals interactions. The radius of curvature Ra at the last point of contact on the CNT with the substrate is fixed by an equilibrium between the adhesion of the part in contact and the bending of the free standing part of the nanotube. From the measurement of the AFM cantilever deflection d (using differential interferometry31) and sample position zs, the force F acting on the nanotube and its compression zc can be recorded.
For the substrates, we start with graphite and mica, since a fresh layer is always easy to be cleaved before the test. The three other substrates we choose to investigate the relative adhesion energies are surfaces of platinum (Pt), gold (Au), and silicon (Si). For the peeling tests, very flat surfaces are necessary in order to get rid of the impact from the morphology of the substrate. We choose to use chips of commercial AFM cantilevers as substrates, bare (thus Si surface) or with a coating of Pt or Au. Before the experiments, these 3 substrates undergo a cleaning process in an ultrasonic bath of ethanol for 10 min, isopropanol for 10 min, and then in a plasma cleaner (medium power) for 10 min to avoid any contamination from the environment. They are finally kept in a clean and dry container for 6 h to eliminate any potential surface impact of the cleaning process.
B. Force curves
In the experiments, the nanotube loop is pushed against a flat sample, as shown in the schematic diagram of Figure 1. The translation of the substrate is performed with a piezo translation platform (Physik Instrumente—PI P527.3) operated in closed loop, featuring an accuracy of 0.3 nm rms. We measure the deflection d of the AFM cantilever with a home-made highly sensitive quadrature phase differential interferometer, which detects the optical path difference between the sensing beam (focused on the cantilever tip) and the reference beam (focused on the static base—see sketch in Figure 1).31,33,34 The deflection d and the sample vertical position zs are simultaneously recorded with high resolution acquisition cards (National Instruments—NI-PXI-4462) at 200 kHz.
With both zs and d being calibrated, using a proper definition of the origins, we can compute at any time the compression of the nanotube
where accounts for the inclination of the cantilever with the substrate. We can also compute the vertical force acting on the nanotube
with ks the static stiffness of the AFM cantilever (calibrated from its thermal noise35). Using compression instead of sample position allows us to take into account the compliance of the cantilever, thus to focus on the nanotube properties only in the force versus compression curves. An example of such a force curve is plotted in Figure 2 for a substrate of graphite. Several examples of force curves for different nanotubes are also provided in the supplementary material.36 A strong hysteresis can be noted during the approach-retract cycle, indicating that the nanotube is changing its configuration during compression. The force is mostly attractive (except at the end of the approach), hinting at adhesion as the main interaction process between the nanotube and the substrate. Most experiments are performed in air, but no difference in the force curves can be observed when experiments are performed in a dry atmosphere of nitrogen or in vacuum at 10−2 mbar: capillary forces are not involved in our experiments. Finally, we note long force plateaus during the retraction: this is the signature of a peeling process.30
Force F of a nanotube as a function of its compression zc on a graphite substrate. A strong hysteresis, due to the adhesion, can be noticed between approach (red, top curve) and retraction (blue, bottom curve). Force plateaus, characteristic of a peeling mechanism, are observed during retraction.
Force F of a nanotube as a function of its compression zc on a graphite substrate. A strong hysteresis, due to the adhesion, can be noticed between approach (red, top curve) and retraction (blue, bottom curve). Force plateaus, characteristic of a peeling mechanism, are observed during retraction.
Let us summarize the analysis presented in Ref. 30 to give grounds for this last claim. We denote by Ea the energy of adhesion per unit length of the nanotube on the surface. As soon as a small length of the nanotube is in contact with the substrate, the system tends to minimize its energy by maximizing the contact length. However, this process increases the bending of the free standing part of the nanotube and the associated curvature energy. The contact length is thus a balance between adhesion and bending, which leads to a constant radius of curvature Ra at the contact point. If the free standing part of the nanotube is long enough compared to Ra, its shape does not change much when it is being peeled from the surface. The vertical displacement δz needed to peel a small length δl from the surface is thus the same length: . In a quasi static displacement, the work produced by the pulling force is thus equivalent to the energy released , leading to (the factor 2 here accounts for the loop geometry of the CNT: we peel the same length δl for the two strands). Peeling a nanotube loop from a surface thus results in a flat force-compression curve, and the value of the force plateau gives directly the value of the adhesive energy per unit length. For the longest plateau of this nanotube on graphite, we read, for example, .
For every CNT and substrate, we perform at least 50 approach-retract cycles: during the first contacts of the nanotube, its configuration usually changes significantly before it reaches a stationary behavior. Some nanotubes can be lost during the process, others never reach a stable operation state. From the grown CNTs, we select only those having a clear peeling process signature: long force plateaus, similar force curves on various substrates, reproducible, and stable in time. Most indeed have many defects resulting in a too complex interpretation of the data (see first figure of supplementary material for an example36). For clarity, we present here the results corresponding to a single nanotube on various substrates. Its force signature is the clearest of our tests, and we observe similar behavior for other samples. No SEM image of this nanotube in its useful configuration is available. It is not visible in the SEM image performed right after the CNT growth. And as for all our nanotubes, it has eventually been lost after extensive force cycling on the various substrates, so no SEM image could be performed a posteriori. The force curves however are qualitatively similar to other CNT loops shown in the supplementary material.36
C. Dynamic stiffness
Following Refs. 30 and 37, during the approach-retract cycles, we also analyse the dynamic properties of the nanotube-substrate contact, by following the evolution of the thermal noise spectrum during the peeling process. Indeed, when the nanotube is in contact with the sample, the contact stiffness sums with the cantilever spring constant and leads to a higher resonance frequency of the system. The fluctuations of the deflection, driven by the random thermal noise excitation of the oscillator, allow us to track this frequency shift during contact, as illustrated in the inset of Figure 3: the power spectrum density (PSD) of the deflection Sd presents the characteristic Lorentzian shape of the thermal noise of a simple harmonic oscillator, peaked at frequencies ranging from 12 kHz (free cantilever, no contact) to 26 kHz (peeling configuration around s).
Time evolution of: (a) the sample position zs, (b) the deflection d, (c) the resonance frequency , and (d) the dynamic stiffness of the nanotube in contact during a slow approach-retract cycle (note the dual vertical scale of this last plot). The PSD Sd of the deflection, computed in a 0.1 s window around s (green, peeling configuration, plateau A), s (yellow, peeling configuration, plateau B), and s (red, no contact), is shown in the inset. The shift of the resonance frequency of the oscillator (cantilever + nanotube in contact) can be used to compute the dynamic stiffness of the contact (see Figure 4). A time-frequency analysis allows one to track the maxima of the PSD computed in every 5 ms window, and then plot this frequency shift as a function of time in (c).
Time evolution of: (a) the sample position zs, (b) the deflection d, (c) the resonance frequency , and (d) the dynamic stiffness of the nanotube in contact during a slow approach-retract cycle (note the dual vertical scale of this last plot). The PSD Sd of the deflection, computed in a 0.1 s window around s (green, peeling configuration, plateau A), s (yellow, peeling configuration, plateau B), and s (red, no contact), is shown in the inset. The shift of the resonance frequency of the oscillator (cantilever + nanotube in contact) can be used to compute the dynamic stiffness of the contact (see Figure 4). A time-frequency analysis allows one to track the maxima of the PSD computed in every 5 ms window, and then plot this frequency shift as a function of time in (c).
From this frequency shift, one can recover the equivalent stiffness of the CNT in contact with the substrate
where f0 and are, respectively, the resonance frequency out of and in contact, and k0 is the free cantilever dynamic stiffness. In this formula, we approximate the cantilever by a simple spring k0, however when the frequency shift is significant as in this measurement, a complete Euler-Bernoulli description of the cantilever is better suited. It leads to the following relation between the dynamic contact stiffness and the resonance frequency of the first mode of the cantilever:37
with
where is the first spatial eigenvalue of a clamped-free Euler-Bernoulli mechanical beam. In Figure 4, we plot both relations 3 and 4: the simple spring approximation is valid till , but the full relation is better suited as soon as the contact stiffness is larger than the cantilever one.
Dynamic stiffness kCNT deduced from the resonance frequency fCNT of the oscillator. When the dynamic stiffness is not too large compared to the static stiffness ks of the free cantilever, a simple harmonic oscillator (SHO) model with a spring is sufficient. For larger values of , a full Euler-Bernoulli model of the cantilever is required: a hard contact produces a finite frequency shift ( when ).
Dynamic stiffness kCNT deduced from the resonance frequency fCNT of the oscillator. When the dynamic stiffness is not too large compared to the static stiffness ks of the free cantilever, a simple harmonic oscillator (SHO) model with a spring is sufficient. For larger values of , a full Euler-Bernoulli model of the cantilever is required: a hard contact produces a finite frequency shift ( when ).
To have a continuous information about the contact stiffness at any compression, we perform a time-frequency analysis of the deflection. We plot in Figures 3(a) and 3(b) the time trace of the sample position zs and the deflection d, sampled at 200 kHz. We slice the latter signal in 5 ms time windows and compute in each window the PSD of the deflection around the time t. Due to the short 5 ms time window, the frequency resolution is only 200 Hz, and the spectra are too noisy to be easily fitted. However, the maximum of the PSD can be measured with the following estimator:
where is an adequate frequency interval centered on (self adapting procedure). The sampling rate for is thus . At such rate, one measures every 0.6 nm of sample position during retraction. We plot the result of this procedure in Figure 3(c). From Eqs. (4) and (5), one finally computes the dynamic stiffness of the nanotube in contact, and plot it in Figure 3(d). At the maximum compression of the nanotube, one sees from this analysis that we reach a hard contact between the AFM tip and the substrate: reaches huge values (several N/m) compared to the cantilever stiffness ( N/m) when .
D. Mean interaction curves for the five substrates
For every substrate, we perform at least 10 very slow approach-retract cycles and compute for each the force and dynamic stiffness versus compression curves during the cycles. We then average these curves for each substrate and plot them in Figure 5. We focus here on the retractions only, which allow to probe the properties of the nanotube in a peeling configuration. The same nanotube is used for all the substrates.
Mean force F and dynamic stiffness kCNT as a function of compression zc measured during retraction on five different substrates: graphite, gold, mica, platinum, and silicon. Each curve is the average of a minimum of 10 approach-retract cycles. Two ranges of compression corresponding to force plateaus are defined: plateau A for nm and plateau B for nm. The zero of compression zs is defined by the loss of contact between the nanotube and the substrate. Missing data in the dynamic stiffness curve for mica correspond to thermal noise spectra that are too noisy to extract the resonance frequency and thus .
Mean force F and dynamic stiffness kCNT as a function of compression zc measured during retraction on five different substrates: graphite, gold, mica, platinum, and silicon. Each curve is the average of a minimum of 10 approach-retract cycles. Two ranges of compression corresponding to force plateaus are defined: plateau A for nm and plateau B for nm. The zero of compression zs is defined by the loss of contact between the nanotube and the substrate. Missing data in the dynamic stiffness curve for mica correspond to thermal noise spectra that are too noisy to extract the resonance frequency and thus .
The curves for the five different substrates have the same generic features: a steep force versus compression dependence and a huge dynamic stiffness at large compression, corresponding to a hard contact between the tip and the surface, then two distinct force plateaus (labelled A and B). The overall behavior is quite different in these two ranges of compression. For the first plateau (A), both F and are rather flat. The average value of increases with the adhesion force on the substrate. On plateau B, the force is also rather flat, but it almost does not depend on the sample (with the exception of graphite). The stiffness presents a significant dependence on compression, with the highest values just before loosing contact, but is rather substrate independent. The shift between the 2 plateaus occurs at slightly different compressions for the five samples (the zero of compression zc is defined by the last contact point during retraction).
III. ADHESION ENERGY AND MECHANICAL PROPERTIES OF THE NANOTUBE
To understand the observed behavior, we use the framework of Ref. 30, corresponding to a straight nanotube initially almost perpendicular to the surface. In such case, a preliminary buckling is necessary to reach the adhered state. In the current work, we consider a nanotube loop that is inherently tangent to the surface before any contact. Adhesion processes are thus immediately relevant. The first contact leads directly to the adhesion of the nanotube. The loss of contact also starts directly with the CNT tangent to the surface, with no intermediate situation where only the tip of a straight nanotube is in contact. The behavior of the loop in contact is however equivalent to having two straight nanotubes sticking in parallel: the two free standing strands of the nanotube are independent since the part sticking to the substrate decouples them. The results of Ref. 30 are thus directly applicable, except that the measured force and contact stiffness should be divided by 2 for each strand of the nanotube!
Each part of the CNT in contact with the sample is thus described as an elastic line, incompressible along its axis. The shape of the line is given by a balance between adhesion and curvature, leading to the radius of curvature Ra at the last contact point before the free standing part
where Ea is the energy of adhesion per unit length, E the nanotube Young's modulus, and its quadratic moment ( its diameter, the wall thickness—0.34 nm for a SWCNT11). As long as the free standing part of the nanotube is much larger than Ra, if we neglect horizontal components of the interaction, the force of interaction with the substrate should be constant and equal to ,30 hence
for a CNT loop.
The static stiffness, defined as is thus zero on the force plateau. However, as clearly illustrated in Figure 3 or 5, the dynamic stiffness measured at the resonance frequency of the oscillator (cantilever + CNT in contact) is not zero: at a few tens of kHz, adhesion has no time to be switched on and off by the small thermal fluctuations, and the sticking part of the CNT can be considered as rigidly clamped to the substrate. In such a case, using the same hypotheses (long nanotube, negligible horizontal forces, and 2 strands), we compute for the dynamic stiffness30
The dynamic stiffness should therefore present a plateau in correspondence to the force one, with values scaling as .
In this analysis framework, we suppose for simplicity that the quadratic moment I of the nanotube is that of an ideal SWCNT. Our samples may differ from this model, as they may be made of bundles of a few single wall nanotubes, carry a significant amount of amorphous carbon,38 or be few-wall nanotubes. The inferred diameter should therefore be considered carefully, but nevertheless provides a reasonable order of magnitude for the probed nanotube.
To test those ideas, we compute for the five substrates the joint histograms of and F during all the retractions and plot them in Figure 6. As clearly illustrated in this figure, the CNT adopts two different configurations during peeling, corresponding to plateaus A and B. We also report in the bottom plot of Figure 6 the mean values of and F on the 2 plateaus defined in Figure 5 for each substrate. The data for plateau A are consistent with our model of peeling a CNT in the direction normal to the surface: constant force F, constant dynamic stiffness , both depending on the substrate material, but all sharing the same nanotube diameter .
To understand the configurations of the CNT during retraction, we plot the joint histograms of and F, for the five substrates. The relation expected between and F for adhered SWCNT loops of various diameters nm is superposed on the graphs. The nanotube explores mainly 2 different configurations: one at low force-low stiffness (plateau B) and one at higher force and stiffness (plateau A). The mean value and standard deviation of and F on the two plateaus are plotted in the bottom graph (plain circles for A, empty circles for B). The configuration of plateau A agrees well with the model, with an estimated diameter nm for all the substrates. The configuration of plateau B would correspond to nm, but show no substrate dependence and is attributed to adhesion on the AFM tip.
To understand the configurations of the CNT during retraction, we plot the joint histograms of and F, for the five substrates. The relation expected between and F for adhered SWCNT loops of various diameters nm is superposed on the graphs. The nanotube explores mainly 2 different configurations: one at low force-low stiffness (plateau B) and one at higher force and stiffness (plateau A). The mean value and standard deviation of and F on the two plateaus are plotted in the bottom graph (plain circles for A, empty circles for B). The configuration of plateau A agrees well with the model, with an estimated diameter nm for all the substrates. The configuration of plateau B would correspond to nm, but show no substrate dependence and is attributed to adhesion on the AFM tip.
The behavior on plateau B is not well described by the model: except for graphite which presents 2 different configurations during this range of compression, both F and depend on compression but not on the substrate. This can be the signature of an adhesion process of the nanotube loop on the silicon tip. Indeed, the two nanotube strands are not perfectly at the tip apex (or the hard contact we reach during the cycle would certainly break their clamping), so they can present some adhesion/peeling mechanism. In such case, the signature should be quite independent of the substrate on the other side of the nanotube. Moreover, the peeling of the nanotube from the AFM tip is not normal to the surface of the silicon tip, so our model may not be suited: the force component parallel to the surface is not negligible any more, so the force and dynamic stiffness can depend on compression. The case of the graphite substrate, implying larger forces, apparently leads to different nanotube configurations in the range of compression corresponding to plateau B.
We finally focus of the range of compression corresponding to plateau A only. The observations of the force and dynamic stiffness plateaus are well described by a peeling configuration, corresponding to a nanotube loop with . The corresponding values of the adhesion energies of the five substrates are reported in Table I: graphite has the strongest interaction with the CNT at nJ/m, while silicon has the lowest at nJ/m. Those numbers are in the same order of magnitude as others reported in the literature,18–24,39–42 though we notice a factor of 3 with our previous measurements on graphite and mica.30 However, the nanotubes of those 2 sets of experiments were grown at a different time in different laboratories and may have some different adhesion properties due to different chiralities, amount of amorphous carbon around the SWCNT,38 number of bundled SWCNT, etc. The ratio between the energy of adhesion is anyway preserved (nanotubes stick twice more on graphite than on mica), so the relative values of the substrates presented in Table I are of broad utility.
Adhesion energies of a SWCNT on various substrates. The values of Ea may differ for different nanotube parameters (diameter, chirality, purity, and additional amorphous carbon layer), but not their relative values. The uncertainty corresponds to the standard deviation in the range of compression corresponding to plateau A for all forces curves on each substrate.
Substrate . | Ea . |
---|---|
Graphite | 2.96 ± 0.33 |
Gold | 2.06 ± 0.35 |
Mica | 1.66 ± 0.12 |
Platinum | 1.39 ± 0.13 |
Silicon | 1.24 ± 0.11 |
Substrate . | Ea . |
---|---|
Graphite | 2.96 ± 0.33 |
Gold | 2.06 ± 0.35 |
Mica | 1.66 ± 0.12 |
Platinum | 1.39 ± 0.13 |
Silicon | 1.24 ± 0.11 |
IV. CONCLUSION
In this article, we present some experiments where a SWCNT loop, grown directly on an AFM tip, is pushed against various substrates. The adhesion force and dynamic stiffness of the nanotube in contact with the sample are recorded as a function of its compression. The experimental data are analyzed in the framework of an elastic loop sticking to a flat surface. During retraction, we observe force plateaus characteristic of a mechanism of peeling. The cross-information between force and dynamic stiffness helps to understand the configuration of the nanotube during the peeling process: some of the plateaus are attributed to adhesion on the substrate, while others are hinting at adhesion on the AFM tip. Quantitative values are derived for the diameter of the nanotube and its energy of adhesion per unit length on various substrates: graphite, gold, mica, platinum, and silicon.
Our experiments illustrate how CNT loops are useful nano-objects to probe peeling processes at the nanoscale, leading to a quantitative measurement of their Van der Walls interaction with flat substrates. This work provides an interesting insight into the physical mechanism of adhesion and should be helpful in the design of nanotube-based nanomechanical devices: if adhesion is sought (for clamping purposes for example), graphite is a good candidate, whereas low interaction can be obtained using a silicon substrate. Other materials of technological relevance could be characterized using our simple protocol.
ACKNOWLEDGMENTS
We thank F. Vittoz and F. Ropars for technical support, L. Champougny, M. Geitner, A. Petrosyan, J. P. Aimé, and Z. Sun for stimulating discussions. This work has been supported by the ANR project HiResAFM (ANR-11-JS04-012-01) of the Agence Nationale de la Recherche in France. Finally, T. J. Li thanks the Chinese Scholar Council (CSC) for financial support. The authors acknowledge the Plateforme Nanofils et Nanotubes Lyonnaise of the University Lyon 1.