Torsional and lateral eigenmode oscillations for atomic resolution imaging of hopg in air under ambient conditions


Torsional and lateral eigenmode oscillations for atomic resolution imaging of hopg in air under ambient conditions

Play all audios:


ABSTRACT Combined in-plane and out-of-plane multifrequency atomic force microscopy techniques have been demonstrated to be important tools to decipher spatial differences of sample surfaces


at the atomic scale. The analysis of physical properties perpendicular to the sample surface is routinely achieved from flexural cantilever oscillations, whereas the interpretation of


in-plane sample properties via force microscopy is still challenging. Besides the torsional oscillation, there is the additional option to exploit the lateral oscillation of the cantilever


for in-plane surface analysis. In this study, we used different multifrequency force microscopy approaches to attain better understanding of the interactions between a super-sharp tip and an


HOPG surface focusing on the discrimination between friction and shear forces. We found that the lateral eigenmode is suitable for the determination of the shear modulus whereas the


torsional eigenmode provides information on local friction forces between tip and sample. Based on the results, we propose that the full set of elastic constants of graphite can be


determined from combined in-plane and out-of-plane multifrequency atomic force microscopy if ultrasmall amplitudes and high force constants are used. SIMILAR CONTENT BEING VIEWED BY OTHERS


SIMULTANEOUS ACQUISITION OF CURRENT AND LATERAL FORCE SIGNALS DURING AFM FOR CHARACTERISING THE PIEZOELECTRIC AND TRIBOELECTRIC EFFECTS OF ZNO NANORODS Article Open access 03 February 2021


ATOMIC-FORCE-MICROSCOPY-BASED TIME-DOMAIN TWO-DIMENSIONAL INFRARED NANOSPECTROSCOPY Article 15 May 2024 OBSERVATION OF NANOSCALE OPTO-MECHANICAL MOLECULAR DAMPING AS THE ORIGIN OF


SPECTROSCOPIC CONTRAST IN PHOTO INDUCED FORCE MICROSCOPY Article Open access 10 November 2020 INTRODUCTION Understanding the in-plane nanomechanical behavior of graphitic surfaces at the


atomic scale under ambient conditions in air is of utmost importance to predict the long-term performance of graphene-based nanodevices such as van der Waals heterostructures1 or for DNA


sequencing2. Multifrequency atomic force microscopy (AFM)3 was shown to be an excellent tool for the quantification of forces at the nanoscale in out-of-plane as well as in in-plane


direction4,5,6. The out-of-plane force deconvolution and determination of elastic moduli from bimodal AFM spectroscopy data using the Sader method7 or the matrix method8 are well


established. In contrast, the determination of in-plane forces and shear moduli from spectroscopic data is still challenging9,10. One issue is the difficulty to excite the cantilever in the


in-plane direction. This problem can either be addressed by using qPlus sensors11,12,13 or by photothermal excitation14,15. The photothermal excitation technique is based on a


power-modulated laser focused at the base of the cantilever and a few micrometers off the cantilever length symmetry axis. This results in an in-plane oscillation if the excitation frequency


matches the resonance frequency of the desired eigenmode. Photothermal excitation facilitated atomic resolution imaging using the torsional-eigenmode oscillation in bimodal AFM16. A second


issue is the discrimination between torsional and lateral eigenmodes. In most publications the term “lateral” is used as synonym for both, the torsional and lateral eigenmode. This might


originate from the fact that in theory the lateral-eigenmode oscillation should not be observable using the beam detection methods as schematically explained by Ding et al_._17 and


illustrated in Fig. S1. However, they rationalized that the tip attached to the cantilever can induce coupling between lateral and torsional eigenmodes facilitating the detection of the


lateral resonance using a standard detection laser and segmented photodiodes. The authors demonstrated that both eigenmodes are useful for imaging and that the lateral oscillation might be


advantageous over the torsional one. The main issue in using the lateral eigenmode for imaging is the calibration of the inverse optical lever sensitivity (invOLS) to attain quantitative or


semi-quantitative data. There are a few approaches for the determination of the torsional invOLS18,19,20,21,22, whereas there is currently no procedure for the determination of the lateral


invOLS to the best of our knowledge. In this study, we compared the suitability of the torsional and the lateral eigenmodes for atomic resolution imaging of HOPG in air under ambient


conditions. The comparison is based on our recently published methodology named AMFlex2-FMTor1-FMFlex3 mode introduced in Ref.4 with a slightly modified setup (AMFlex2-FMLat1-FMFlex3). Both


setups exploit the second flexural eigenmode for the topographical feedback in amplitude modulation (AM) whereas the first torsional, first lateral and third flexural mode are


frequency-modulated (FM) and controlled by phase-locked-loop (PLL) electronics. Inspired by the work of colleagues in the field, we used higher flexural eigenmodes which were shown to be


beneficial for atomic resolution imaging due to the enhanced stiffness compared to the first flexural eigenmode23,24. Additionally, the combination of small free amplitudes and small


amplitude setpoints was used for imaging as it was shown by Santos et al_._25,26 and Lai et al_._27 that the method can lead to a very close proximity between tip and sample, being essential


for high resolution imaging. For some of the super-sharp tips we used, it was impossible to excite the second flexural and the first torsional eigenmode individually due to the close


proximity of their resonance frequencies. Interestingly, the coupling of the two eigenmodes facilitated atomic resolution imaging. A scheme of possible setups as well as the resulting


movement of the tip are schematically shown in Fig. 1. If there is coupling between the second flexural and the first torsional eigenmode, the torsional oscillation cannot be controlled


individually (open loop (OL)). From the second flexural and the first torsional phase shifts, dissipation in out-of-plane and in-plane direction can be determined, respectively.


Unfortunately, the AMFlex2-OLTor1 mode is not suitable for the deconvolution of in-plane or out-of-plane forces. Consequently, we extended the AMFlex2-OLTor1 mode, by two phase-locked loops


(PLL) for frequency modulation (FM) of the first lateral and the third flexural eigenmode. The resulting AMFlex2-OLTor1-FMLat1-FMFlex3 mode such as schematically depicted in Fig. 1


facilitates quantification of in-plane and out-of-plane forces, reconstructed from lateral- and third-eigenmode flexural frequency-shift data. In this work, we analyzed the imaging


capability of coupled flexural/torsional oscillations as well as the suitability of the first lateral eigenmode for atomic resolution imaging and provide an approach for the calibration of


the lateral eigenmode optical lever sensitivity. In addition, we calculated in-plane forces from torsional and lateral frequency-shift data and analyzed the origin of the forces regarding


friction or shear. Analyzing the tip-trajectories, effective torsional and lateral amplitudes as well as the maximum indentation of the tip into the surface (with and without compression)


were estimated. The plausibility of the determined lateral forces was confirmed by calculating the lateral displacement using the literature value of the shear modulus of HOPG. RESULTS AND


DISCUSSION IMAGING CAPABILITY OF THE AMFLEX2-OLTOR1-FMLAT1 MODE FOR ATOMIC RESOLUTION From Fig. 2, we get an overview of the capability of different imaging channels for atomic resolution


imaging in air under ambient conditions depending on the _z_-sensor position. Increasing _z_-sensor position (from left to right) causes a reduction of the average tip-sample distance.


Figure 2a–e show schemes of the expected tip-trajectory corresponding to the AMFlex2-OLTor1-FMLat1 mode with a lateral-eigenmode amplitude setpoint of _A__lat,_1 = 863 pm. The shown


tip-trajectories were calculated by plotting the out-of-plane deflection _z_(_t_) _vs._ the in-plane deflection _x_(_t_) of the cantilever using linear combinations of cosine functions such


as shown by Benaglia et al_._28 for the _z_-component and assuming that _Φ__tor,_1 = 90° (see Eqs. (S1, S2) and Fig. S2 of the Supplementary Information for details). The lateral-eigenmode


invOLS was calibrated by imaging a wrinkle of a graphene layer on an HOPG sample with different lateral amplitude setpoints while oscillating perpendicular to the wrinkle and comparing the


images with the results of a Savitzky-Golay-filtered reference image. A detailed description of the technique can be found in Section 4 (Calibration of lateral-oscillation-eigenmode


sensitivity) of the Supplementary Information. From Fig. 2a–e the second flexural-eigenmode-amplitude setpoint is reduced from 700 to 454 pm, causing a reduction of the torsional amplitude


due to the coupling of both modes (Fig. 2, second row). The corresponding approximate _z_-sensor position can be taken from the third row of Fig. 2. From spectroscopic experiments (local


amplitude distance curves recording all available observables), we can visualize the dependence of the second flexural amplitude (Fig. 2f), the first torsional amplitude (Fig. 2l) and the


first lateral frequency shift (Fig. 2r) on the _z_-sensor position. The crossed circles mark the positions where the AFM images of Fig. 2 where taken. The height images are shown in Fig. 


2g–k, the “error images” of the first torsional amplitude in Fig. 2m–q and the first lateral frequency-shift images in Fig. 2s–w. “Error images” of the first torsional amplitude here means


the deviation of the actual torsional amplitude relative to the torsional-amplitude setpoint corresponding to the second flexural-amplitude setpoint due to the coupling of both eigenmodes.


Comparing the amplitude _vs. z_-position curves in Fig. 2f,l, it becomes evident that although the second flexural and the first torsional oscillations are coupled, they show different


dependencies on the tip-sample distance. These differences are consistent with results from uncoupled flexural–torsional AFM studies in open-loop configuration16. The lateral frequency shift


_vs. z_-sensor position curve in Fig. 2r shows an increasing repulsive interaction with decreasing tip-sample distance. Atomic resolution imaging was feasible in all three channels as e.g.


visible in Fig. 2k,q,w at a _z_-sensor position of approximately 3.7 nm. Interestingly, the lateral frequency-shift images show atomically resolved structures at setpoint ratios up to 93% of


the second flexural-eigenmode amplitude as observable in Fig. 2s, which corresponds to a _z_-sensor value of approximately 2.2 nm. Compared to imaging with “uncoupled” cantilevers, where


atomic resolution imaging was impossible to achieve using setpoint ratios larger than 15%, it is reasonable to believe that performing measurements with “coupled” cantilevers implies more


gentle conditions, preserving the integrity of both, tip and sample. Note, that the torsional amplitude ratio was already at 63% which fits to our observation that the stiffness of the


second flexural eigenmode was higher compared to the first torsional eigenmode stiffness for the type of cantilevers used in this study (see “Materials and methods” for details). At this


position, however, neither in the height image (Fig. 2g) nor in the torsional-amplitude error image (Fig. 2m) atomic contrast could be resolved. Lowering the second


flexural-eigenmode-amplitude setpoint, atomic contrast becomes more and more apparent in the height images (Fig. 2m–q). The growing corrugation amplitude (measured vertical distance between


atomic and hollow side) is assumed to originate from an increase in the mean cantilever deflection16,29,30. Strikingly, the height image at a _z_-sensor position of 3.7 nm (Fig. 2k) shows a


step which matches the height of individual graphene layers. This step is neither observable in the Δ_A__tor_,1 image (Fig. 2q) nor in the Δ_f__lat_,1 image at the same _z_-sensor position.


Consequently, the topographic feedback seems to function extremely precise. Another interesting observation is, that the atomic structure at both, the lower and the upper plateau of the


graphene step, exhibit the same arrangement of carbon atoms. This leads us to the conclusion that we imaged a monoatomic step edge covered by a small number of graphene layers which is in


line with results shown by Abooalizadeh et al_._31. From this observation the question arises why we were not able to resolve the step at smaller _z_-sensor positions. To answer the


question, we exemplarily drew two cross sectional profiles through the height images shown in Fig. 2k,j (see Fig. S3 in the Supplementary Information). We assume that the visibility of the


step edge depends on whether the average tip-sample force (represented by the mean deflection, averaged over several oscillation cycles) is predominantly attractive or repulsive. On the


right side of Fig. S3, we schematically sketched the alleged interaction between tip and sample on a covered step edge. In the predominantly attractive regime, we assume that several


graphene layers cover the step such as depicted in the top right part of Fig. S3. As discussed in our recent work and by others4,32, it is reasonable to believe that the topmost carbon


layers are lifted by the attractive tip-sample interactions. Consequently, the tip will not sense the step edge. If imaging takes place in the predominantly repulsive regime, we assume that


the graphene cover layers are tightly stretched over the step edge as sketched in the bottom right of Fig. S3. From the imaging results shown in Fig. 2 we can conclude that the coupled


motion of the cantilever (AMFlex2-OLTor1) seems to promote atomic resolution at higher setpoint ratios (_A__flex,_2/_A_0_(flex,_2_)_) compared to uncoupled ones. This can, on the one hand,


be attributed to the enhanced dynamic stiffness of the coupled flexural/torsional oscillation of the cantilever and, on the other hand, to the oval shape of the tip trajectory. This might


lead to a reduced influence of the interactions far away from the center of the oscillation, because the tip will have the closest distance to the sample and therefore the highest


interaction at the currently measured local position. Nevertheless, the largest drawback of imaging with coupled cantilever modes is that the relation between the flexural and the torsional


amplitudes cannot be set individually. Additionally, it needs to be considered that the combined in-plane oscillation of torsional and lateral components needs to be analyzed in detail in


order to get a better understanding regarding the interaction with the sample. An idea for future studies would be tailoring cantilevers e.g. by a controlled introduction of holes at


different positions of the cantilever such as shown by Eslami et al_._33, in order to analyze the influence of different resonance-frequency ratios. COMPARISON OF TORSIONAL AND LATERAL


FREQUENCY-SHIFT IMAGES AS A FUNCTION OF THE AMPLITUDE SETPOINTS In order to gain deeper insights into the different oscillation behaviors of the torsional and the lateral eigenmodes, we


analyzed several 5 × 5 nm2 frequency-shift images taken in the AMFlex2-FMLat1-FMFlex3 mode and the AMFlex2-FMTor1-FMFlex3 mode at different lateral- and torsional-eigenmode setpoints. The


results are shown in Fig. 3. In Fig. 3a,e and Fig. 3d,h lateral and torsional frequency-shift images (2.5 × 2.5 nm2 zoom-in) taken at _A__lat,_1 = 518 pm (a), _A__lat,_1 = 3451 pm (e) and


_A__tor,_1 = 158 pm (d), _A__tor,_1 = 628 pm (h) are shown exemplarily. The corresponding histograms of the lateral and the torsional frequency-shift images can be seen in Fig. 3b,c,


respectively. The open circles in Fig. 3f,g show the lateral and torsional frequency shifts at maximum counts and the dashed lines represent the full width at half maximum (FWHM) as a


function of the lateral- and torsional-amplitude setpoints, respectively. If we compare the frequency-shift images taken in the AMFlex2-FMLat1-FMFlex3 mode in Fig. 3a,e with the images taken


in the AMFlex2-FMTor1-FMFlex3 mode in Fig. 3d,h we observe that atomic contrast was achieved in all images but the range of the frequency-shift values strongly differs. While the lateral


and torsional frequency-shift images taken at the higher-amplitude setpoints (_A__lat,_1 = 3451 pm and _A__tor,_1 = 628 pm) shown in Fig. 3e,h show similar frequency-shift ranges, the


frequency-shift images at smaller-amplitude setpoints (_A__lat,_1 = 518 pm and _A__tor,_1 = 158 pm) in Fig. 3a,e clearly differ regarding their frequency-shift ranges and the distribution of


the values. This becomes even more evident by looking at the histograms corresponding to frequency-shift images taken at five different lateral- (Fig. 3b) and torsional-amplitude (Fig. 3c)


setpoints. In general, the full width at half maximum (FWHM) of the peaks increases with decreasing amplitude setpoints (lateral/torsional). The histograms of the lateral frequency-shift


images remain symmetric around a center value of approximately 120 Hz, whereas the histograms of the torsional frequency-shift images become asymmetric for torsional-amplitude setpoints


smaller than 314 pm. As a result, the maxima of the histograms are shifted to more repulsive values. Moreover, the third eigenmode flexural frequency shift is only slightly influenced by the


lateral- or torsional-eigenmode amplitude which corroborates the assumption of the independency of the in-plane (lateral/torsional) and the out-of-plane (flexural) eigenmodes (Fig. 3f,g


(top graphs)). Interestingly, we observe that although the torsional- and lateral-amplitude setpoints are in most cases substantially larger than the interatomic spacing of the carbon atoms,


we can still obtain atomic contrast in the frequency-shift images. To further investigate this phenomenon, we calculated the tip-trajectory resulting from the combined lateral-flexural


(Fig. 4a,b) or torsional-flexural (Fig. 4c,d) oscillations of the cantilever according to equations (S1, S2) for two different lateral and torsional amplitudes, respectively. It needs to be


mentioned that in Fig. 4 only an extract of the full oscillation trajectory is shown. While imaging one position, almost every point in the area spanned by _A__flex_,2 + _A__flex_,3 and


_A__tor,_1 or _A__lat,_1, respectively, is reached such as plotted in the Supplementary Information (Fig. S6). From Fig. 4 it becomes obvious that the interaction between tip and sample


strongly depends on the curvature of the tip. For the situations shown in Fig. 4a,b,d we interpret, that the minimum distance from the lower turning point of the tip-trajectory to the


intersection of the tip-trajectory and the horizontal line (purple dotted) defined by the interatomic spacing between two atoms (here 246 pm, due to the tip oscillating perpendicular to the


carbon bonds) determines the indent (without distortion of the carbon atoms) of the tip into the carbon ring (vertical purple line). This indent without distortion of carbon atoms


corresponds to the corrugation amplitude (_A__C-H_). Thus, we defined the effective in-plane amplitudes as half of the horizontal distance between two carbon atoms (here: _A__lat,_1_ (eff)_ 


= _A__tor,_1 _(eff)_ = 123 pm because the tip oscillates perpendicular to the carbon bonds) for the situations shown in Fig. 4a,b,d. From the situation shown in Fig. 4c we found that if the


tip-trajectory shows a higher curvature, the effective amplitude is no longer equivalent to half of the distance between two carbon bonds. This is a result of the maximum possible


corrugation amplitude which was reported by Kawai et al_._ to equal 85 pm in the repulsive regime32. This value was additionally verified by performing dynamic spectroscopy experiments such


as shown in Fig. S5 of the Supplementary Information. Consequently, the effective amplitude is still determined by the intersection between the tip-trajectory and the horizontal line (purple


dotted) defined by the interatomic spacing between two atoms, but it becomes smaller than 123 pm due to the limited corrugation amplitude. If we calculate the effective torsional amplitude


from the tip-trajectory in Fig. 4c, assuming a corrugation amplitude of 85 pm, it would be approximately 83 pm, which is in line with the 79 pm determined in our recent work4. We can


conclude that imaging the atomic structure of HOPG in the frequency-shift channel is also possible with in-plane oscillations larger than half of the interatomic spacing, but the resolution


in the height images is strongly reduced due to the small corrugation amplitudes. Additionally, it needs to be considered that the frequency-shift images reflect the tip-sample interaction


averaged over a few carbon hexagons. This needs to be kept in mind if forces are determined from the frequency-shift images. Nevertheless, due to the formation of plateaus in the torsional


and lateral frequency shift values for larger in-plane oscillation amplitudes (see Fig. 3f,g (bottom graph)), we propose that conclusions about different interaction mechanisms, _i.e._


friction or shear, as a function of the different in-plane oscillation amplitudes can be drawn. We point out that the most obvious way of determining the corrugation amplitude would be the


analysis of the topography images. Unfortunately, the topography images are strongly influenced by the mean cantilever deflection, so that the corrugation amplitude is not directly


accessible4. CALCULATION OF LATERAL AND TORSIONAL FORCES From Fig. 3, we additionally observed that the full width at half maximum becomes much broader for the smaller lateral-amplitude


setpoints compared to the torsional-amplitude setpoints. In order to identify the origin of the differences between torsional and lateral frequency-shift behavior, we calculated the forces


from the images taken at _A__lat,_1 = 518 pm and _A__tor,_1 = 158 pm (_A__flex,_2 = 110 pm and _A__flex,_3 = 145 pm), respectively. The results are presented in Fig. 5. Schemes of the


interaction between tip and sample are shown in Fig. 5a,d. The lateral and torsional frequency-shift images are shown in Fig. 5b,e. In Fig. 5c,f forces calculated from the frequency-shift


images in Fig. 5b,e by using the Fourier method4,35 are shown. Obviously, both, the torsional and the lateral force image in Fig. 5c,f reproduce the hexagonal structure of graphene. It needs


to be mentioned that this is only the case if the forces are calculated by using the effective torsional or lateral amplitude. Interestingly, forces determined from the lateral


frequency-shift image in Fig. 5b are around 50 times larger compared to the forces (Fig. 5f) calculated from the torsional frequency-shift images in Fig. 5e. We propose that this large


difference in forces originates partially from the difference in force constants but also/mainly from different types of interaction mechanisms between tip and sample, i.e. friction and


shear. While Weymouth et al_._11,36 observed “non-contact friction” during imaging of an H-terminated Si(100) surface in dynamic lateral force microscopy using a qPlus sensor13, we assume


that the lateral frequency shift determined on HOPG at 518 pm lateral amplitude mainly originates from a local shear due to relatively weak van-der-Waals interactions between the carbon


layers. We assume this to happen for two reasons: First, the angle under which the tip touches the surface at the respective local position is less steep for the shown lateral oscillation


compared to the torsional one. Second, the lateral stiffness of the cantilever is higher compared to the torsional stiffness which might lead to shear for a certain range of lateral


oscillation amplitudes whereas torsional amplitudes of the same magnitude result in frictional interaction between tip and sample. In order to verify our assumption about shearing


interaction for the constellation shown in Fig. 5a–c we calculated the lateral displacement _Δx_ of the topmost carbon layer relative to the second one if we assume a shear modulus of 4.18 


GPa in x-z-direction37. In a simple approximation, the shear modulus _G__xz_ can be calculated by $${G}_{xz}=\frac{{F}_{lat}\cdot d}{\Delta x\cdot {A}_{c}},$$ (1) where _d_ is the distance


between two carbon layers (334 pm) and _A__c_ is the contact area between tip and sample which can be calculated according to Hertz from the tip radius _R_ and the indentation _δ_38


$${A}_{c}=\pi R\delta .$$ (2) In order to determine the indentation _δ_ which is caused by the tip oscillating over a carbon atom, we used Hertz contact mechanics model38


$$F=\frac{4}{3}{E}^{*}\sqrt{R{\delta }^{3}},$$ (3) where _E_* is the effective Young’s modulus which can be considered here as equal to the Young’s modulus _E_zz of graphite due to the high


stiffness of the diamond-like-carbon tip. The maximum value of _F__flex,_3 was calculated in our resent work to be approximately 520 pN for the AMFlex2-FMTor1-FMFlex3 mode4. Due to the


similarity of the third-eigenmode flexural frequency-shift values we chose this value as representative for all the other experiments. By rearranging Eq. (3) and inserting _E_zz = 36.5 GPa


as reported by Blakslee et al_._39 we determined an indentation of approximately 50 pm. Using _δ_ = 50 pm and rearranging Eq. (1), we calculated a local displacement of approximately _Δx_ = 


157 pm between the first and the second graphene layer. This implies a shift by the length of approximately one carbon–carbon bond (142 pm) which seems to be reasonable in distance. We


corroborated our findings by performing additional friction force microscopy images, however, using a comparably large vertical deflection setpoint (43 nm ≙ 1100 nN) in order to ensure


permanent contact between tip and sample. By means of the analysis of cross sections drawn through lateral deflection trace and retrace images, the actual lateral tip-position was


calculated. The results are shown in Fig. S5 of the Supplementary Information. As expected, the typical stick–slip movement of the tip imaging an HOPG surface was observed40, however, the


stiction was stretched out over a region of approximately 150 pm, which can be interpreted as a temporary shift of the topmost carbon layer relative to the second one. For the tip-trajectory


of the torsional oscillation shown in Fig. 5d, we observe a very similar oscillatory behavior, however, much more compact in the in-plane direction, resulting from the smaller torsional


amplitude. As a consequence, the angle at which the tip approaches and finally touches the surface is much steeper and therefore a negligible shifting of the first carbon layer relative to


the second is assumed to take place. We suppose that the torsional forces involved in this mechanism can be interpreted as strongly localized frictional forces which is on the one hand in


accordance with the magnitude of frictional forces on graphene determined by Lee et al_._41. On the other hand, the interpretation is in compliance with the observation of direction


dependent torsional forces observed in our recent work4. Additionally, it needs to be highlighted, that the flexural/torsional oscillation of the cantilever simulated here is symmetric in


out-of-plane direction but not in in-plane direction such as schematically shown in Fig. S6 of the Supplementary Information. This can explain the shift between topography and


frequency-shift images such as shown in our recent work4. For the flexural/lateral oscillation the tip-trajectory is overall symmetric in both, in-plane and out-of-plane direction (see Fig.


S6). Here, no shift between the height and the frequency-shift images could be observed. However, it needs to be mentioned that the resolution in the height images was strongly reduced if


imaging was accomplished with the lateral instead of the torsional eigenmode. The effect can be attributed to the shear interaction between tip and sample. If the lateral or the torsional


amplitudes become larger in size, we observe averaged friction along the surface which might explain the formation of the plateaus in the frequency shifts in Fig. 3f,g (bottom graph). Here,


we also observed that the flexural frequency-shift values and therefore the flexural forces are only slightly influenced by the in-plane oscillation amplitudes. In theory, it should


additionally be possible to determine the in-plane Young’s moduli _E__xx_ = _E__yy_ and _E__xy_ by using higher torsional or lateral eigenmodes of the cantilever. Unfortunately, it is


nontrivial to control these higher eigenmodes and determine their invOLS as well as the corresponding force constants. Nevertheless, this approach could be considered as an alternative


method for the determination of the Young’s modulus of graphene instead of the frequently used techniques where graphene is spanned over membranes or holes and _E__xx_ is calculated from


nanoindentation experiments42,43. Therefore, we propose that the simultaneous determination of strongly localized in-plane and out-of-plane sample properties such as shear moduli and Young’s


moduli using multifrequency AFM is feasible if very small amplitudes and adequate force constants are used for the analysis. The main advantage of determining elastic properties at the


atomic scale is the possibility to locally assess the influence of atomic defects in the material, which were shown to strongly influence the mechanical stability of graphene samples on a


larger scale43. Serving the large scientific interest to analyze the origin of friction anisotropy on graphene and graphite samples44,45,46,47, we are convinced that our presented


multifrequency AFM method can substantially contribute to gain deeper insights into in-plane sample properties by comparing the mechanical properties on the atomic scale with that on the


nano-/micrometer scale. The ability to compare material properties at different length scales in air under ambient conditions is of particular interest for the analysis of adsorbate


formation46,47 or ripples44,45 on a graphitic surface which strongly influences the performance of graphene-based nanodevices. Additionally, the availability of torsional and lateral


eigenmodes extends the spectrum of oscillation modes, facilitating the simultaneous acquisition of different in-plane sample properties, using one and the same cantilever. CONCLUSION In


summary, we demonstrated that performing multifrequency AFM with in-plane and out-of-plane components of the tip-motion has additional potential for atomic resolution imaging of HOPG


surfaces. The coupling of the second flexural and the first torsional eigenmode facilitates imaging with atomic contrast at higher amplitude-setpoint ratios, although the quantitative


interpretation of the results is challenging if the torsional and the lateral oscillation amplitude are simultaneously excited. In order to estimate the inverse optical lever sensitivity of


the lateral cantilever eigenmode, a calibration procedure was proposed which is based on imaging a graphene wrinkle with different lateral amplitude setpoints while oscillating perpendicular


to the wrinkle and comparing the images with the results of a Savitzky-Golay-filtered reference image. We compared the results from imaging in the AMFlex2-FMLat1-FMFlex3 and the


AMFlex2-FMTor1-FMFlex3 mode. Analyzing the tip-trajectories and the in-plane forces resulting from both modes with the respective amplitudes, we proposed that imaging with small lateral


amplitudes can result in localized shear of the first graphite layer with respect to the second one with a relative displacement of 157 pm. The applied out-of-plane forces caused an


indentation of approximately 50 pm of the top graphene layer. Both values are in a reasonable range to propose that the AMFlex2-FMLat1-FMFlex3 method can be used for the simultaneous


determination of _E__xz_ = _E__yz_ and _G__xz_ = _G__yz_. Imaging with the AMFlex2-FMTor1-FMFlex3 in contrast, promotes rather frictional forces which we mainly attributed to the smaller


force constant and the curved tip-trajectory compared to the lateral eigenmode. Based on our findings, we suggest that the use of higher in-plane cantilever eigenmodes or generally stiffer


cantilevers can provide a strategy to determine the shear moduli on HOPG and graphene also with torsional eigenmodes. The method carries great potential for future assessment of graphene in


nanodevices where local differences in mechanical properties e.g. induced by defects or adsorbates, play a major role. To this end, the in-plane oscillation amplitudes need to be in the


range of the interatomic spacings, which is still challenging to achieve while imaging in air under ambient conditions due to the small signal-to-noise ratio. With the presented method we


aim to promote the analysis of friction anisotropy observed on graphene by the comparison of atomic and nano-/micrometer resolution images at different relative orientations between the


hexagonal carbon lattice and the in-plane cantilever oscillation using torsional and/or lateral eigenmodes. This can help to investigate the origin of fundamental friction mechanisms


underlaying graphene-based systems. Additionally, we expect that the presented method might not only be limited to the analysis of stiff samples but also can be applied for a broad range of


materials, including soft matter, facilitating a complete in- and out-of-plane sample surface characterization. MATERIALS AND METHODS SAMPLE The HOPG sample (grade 2, mosaic spread angle:


0.8 ± 0.2°) was purchased from SPI Supplies (Structure Probe, Inc., West Chester, PA, USA). Cleaving was performed prior to the AFM experiments with adhesive tape in air under ambient


conditions to expose a fresh clean surface. CANTILEVERS Super-sharp cantilevers of the type HiResC15/Cr-Au purchased from Mikromasch (Innovative Solutions Bulgaria Ltd., Sofia, Bulgaria)


were used. Although of the same type, the cantilevers showed slightly varying resonance frequencies resulting in “coupled” and “uncoupled” eigenmodes. For the AFM images shown in Fig. 2 we


used a cantilever with “coupled” second flexural and first torsional eigenmodes as it depicted a close proximity of both resonance frequencies. For the AFM images shown in Figs. 3 and 5


cantilevers with “uncoupled” eigenmodes were used. The resonance frequencies, quality factors and force constants of the used cantilevers are listed in Table 1. Details for the determination


of the force constants are given in Section 7 of the Supplementary Information. The tip radius of the cantilevers was _R_ = 1 nm according to the manufacturer's data sheet. In Section 


8 of the Supplementary Information we list the ratios between the flexural, torsional and lateral resonance frequencies determined for several cantilevers of the type HiResC15/Cr-Au for


further considerations. ENVIRONMENTAL CONDITIONS The AFM-lab was equipped with a controlled ventilation system, which provided stable environmental conditions also inside the AFM chamber.


The temperature and the relative humidity in the AFM chamber were tracked by sensors and remained stable during imaging (relative humidity: 22 ± 2%, temperature: 26 ± 2 °C). AFM SETUP A


Cypher S atomic force microscope (Asylum Research, Oxford Instruments, Santa Barbara, CA, USA) equipped with a built-in blueDrive photothermal excitation setup for dynamic AFM modes was used


for all experiments. The blueDrive laser was focused at the fixed end of the cantilever and the lateral position was optimized by finding the spot that led to the maximum possible in-plane


amplitude. A scheme for the approximate positioning of the laser spots on the cantilever is shown in Fig. S1e of the Supplementary Information. The respective frequency shifts (depending on


the method used) were tracked with additional phase-locked loops (HF2PLL, Zurich Instruments, Zurich, Switzerland). The drive amplitudes of the frequency-modulated eigenmodes were adjusted


in order to maintain constant amplitudes by two PID controllers implemented in the same instrument. DATA PROCESSING The topography and the torsional-amplitude error images were first-order


flattened to remove any tilt from the images using the Igor Pro v6.36 software (WaveMetrics Inc., Lake Oswego, OR, USA). In order to remove noise from the small-scale images a 3 × 2 Gauss


filter was applied to all AFM images except for the ones used for the determination of the lateral invOLS (Fig. S4). For the calculation of the torsional and the lateral force images in Fig.


 5 we wrote a Matlab code (MATLAB R2018a, MathWorks Inc., Natick, MA, USA) for the implementation of the Fourier method based on the script provided by Seeholzer et al_._4,35. The calculated


force images were smoothed with a Savitzky–Golay filter over nine points with a first-order polynomial. DATA AVAILABILITY All data will be made available from the authors upon reasonable


request. REFERENCES * Wang, J., Zhang, D. & Zhou, B. Achieving an Ohmic contact in graphene-based van der Waals heterostructures by intrinsic defects and the inner polarized electric


field of Janus AlGaSSe. _New J. Chem._ 45, 21178–21187. https://doi.org/10.1039/D1NJ03861C (2021). Article  CAS  Google Scholar  * Heerema, S. J. & Dekker, C. Graphene nanodevices for


DNA sequencing. _Nat. Nanotechnol._ 11, 127–136. https://doi.org/10.1038/nnano.2015.307 (2016). Article  ADS  CAS  PubMed  Google Scholar  * Garcia, R. & Herruzo, E. T. The emergence of


multifrequency force microscopy. _Nat. Nanotechnol._ 7, 217–226. https://doi.org/10.1038/nnano.2012.38 (2012). Article  ADS  CAS  PubMed  Google Scholar  * Eichhorn, A. L. & Dietz, C.


Simultaneous deconvolution of in-plane and out-of-plane forces of HOPG at the atomic scale under ambient conditions by multifrequency atomic force microscopy. _Adv. Mater. Interfaces_ 8,


2101288. https://doi.org/10.1002/admi.202101288 (2021). Article  CAS  Google Scholar  * Uhlig, M. R., Martin-Jimenez, D. & Garcia, R. Atomic-scale mapping of hydrophobic layers on


graphene and few-layer MoS2 and WSe2 in water. _Nat. Commun._ 10, 2606. https://doi.org/10.1038/s41467-019-10740-w (2019). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  *


Baykara, M. Z. _et al._ Exploring atomic-scale lateral forces in the attractive regime: A case study on graphite (0001). _Nanotechnology_ 23, 405703.


https://doi.org/10.1088/0957-4484/23/40/405703 (2012). Article  ADS  CAS  PubMed  Google Scholar  * Sader, J. E. & Jarvis, S. P. Accurate formulas for interaction force and energy in


frequency modulation force spectroscopy. _Appl. Phys. Lett._ 84, 1801–1803. https://doi.org/10.1063/1.1667267 (2004). Article  ADS  CAS  Google Scholar  * Giessibl, F. J. A direct method to


calculate tip–sample forces from frequency shifts in frequency-modulation atomic force microscopy. _Appl. Phys. Lett._ 78, 123–125. https://doi.org/10.1063/1.1335546 (2000). Article  ADS 


Google Scholar  * Kawai, S. _et al._ Ultrasensitive detection of lateral atomic-scale interactions on graphite (0001) via bimodal dynamic force measurements. _Phys. Rev. B_ 81, 085420.


https://doi.org/10.1103/PhysRevB.81.085420 (2010). Article  ADS  CAS  Google Scholar  * Dietz, C. Sensing in-plane nanomechanical surface and sub-surface properties of polymers: Local shear


stress as function of the indentation depth. _Nanoscale_ 10, 460–468. https://doi.org/10.1039/c7nr07147g (2018). Article  CAS  Google Scholar  * Weymouth, A. J. _et al._ Atomic structure


affects the directional dependence of friction. _Phys. Rev. Lett._ 111, 126103. https://doi.org/10.1103/PhysRevLett.111.126103 (2013). Article  ADS  CAS  PubMed  Google Scholar  *


Martin-Jimenez, D. _et al._ Chemical bond imaging using torsional and flexural higher eigenmodes of qPlus sensors. _Nanoscale_ 14, 5329–5339. https://doi.org/10.1039/D2NR01062C (2022).


Article  CAS  PubMed  Google Scholar  * Giessibl, F. J. The qPlus sensor, a powerful core for the atomic force microscope. _Rev. Sci. Instrum._ 90, 011101. https://doi.org/10.1063/1.5052264


(2019). Article  ADS  CAS  PubMed  Google Scholar  * Pini, V., Tiribilli, B., Gambi, C. M. C. & Vassalli, M. Dynamical characterization of vibrating AFM cantilevers forced by


photothermal excitation. _Phys. Rev. B_ 81, 054302. https://doi.org/10.1103/PhysRevB.81.054302 (2010). Article  ADS  CAS  Google Scholar  * Kiracofe, D., Kobayashi, K., Labuda, A., Raman, A.


& Yamada, H. High efficiency laser photothermal excitation of microcantilever vibrations in air and liquids. _Rev. Sci. Instrum._ 82, 013702. https://doi.org/10.1063/1.3518965 (2011).


Article  ADS  CAS  PubMed  Google Scholar  * Meier, T., Eslami, B. & Solares, S. D. Multifrequency force microscopy using flexural and torsional modes by photothermal excitation in


liquid: Atomic resolution imaging of calcite (10(1)over-bar4). _Nanotechnology_ 27, 085702. https://doi.org/10.1088/0957-4484/27/8/085702 (2016). Article  ADS  CAS  PubMed  Google Scholar  *


Ding, R.-F., Yang, C.-W., Huang, K.-Y. & Hwang, I.-S. High-sensitivity imaging with lateral resonance mode atomic force microscopy. _Nanoscale_ 8, 18421–18427.


https://doi.org/10.1039/C6NR04151E (2016). Article  CAS  PubMed  Google Scholar  * Mullin, N. & Hobbs, J. K. A non-contact, thermal noise based method for the calibration of lateral


deflection sensitivity in atomic force microscopy. _Rev. Sci. Instrum._ 85, 113703. https://doi.org/10.1063/1.4901221 (2014). Article  CAS  PubMed  Google Scholar  * Ogletree, D. F.,


Carpick, R. W. & Salmeron, M. Calibration of frictional forces in atomic force microscopy. _Rev. Sci. Instrum._ 67, 3298–3306. https://doi.org/10.1063/1.1147411 (1996). Article  ADS  CAS


  Google Scholar  * Varenberg, M., Etsion, I. & Halperin, G. An improved wedge calibration method for lateral force in atomic force microscopy. _Rev. Sci. Instrum._ 74, 3362–3367.


https://doi.org/10.1063/1.1584082 (2003). Article  ADS  CAS  Google Scholar  * Zhang, G., Li, P., Wei, D., Hu, K. & Qiu, X. Lateral force calibration for atomic force microscope


cantilevers using a suspended nanowire. _Nanotechnology_ 31, 475703. https://doi.org/10.1088/1361-6528/abae32 (2020). Article  ADS  CAS  PubMed  Google Scholar  * Munz, M. Force calibration


in lateral force microscopy: A review of the experimental methods. _J. Phys. D Appl. Phys._ 43, 063001. https://doi.org/10.1088/0022-3727/43/6/063001 (2010). Article  ADS  CAS  Google


Scholar  * Kawai, S., Kitamura, S.-I., Kobayashi, D., Meguro, S. & Kawakatsu, H. An ultrasmall amplitude operation of dynamic force microscopy with second flexural mode. _Appl. Phys.


Lett._ 86, 193107. https://doi.org/10.1063/1.1923200 (2005). Article  ADS  CAS  Google Scholar  * Korolkov, V. V. _et al._ Supramolecular networks stabilise and functionalise black


phosphorus. _Nat. Commun._ 8, 1385. https://doi.org/10.1038/s41467-017-01797-6 (2017). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Santos, S., Olukan, T. A., Lai, C.-Y.


& Chiesa, M. Hydration dynamics and the future of small-amplitude AFM imaging in air. _Molecules_ 26, 7083. https://doi.org/10.3390/molecules26237083 (2021). Article  CAS  PubMed  PubMed


Central  Google Scholar  * Santos, S. _et al._ Stability, resolution, and ultra-low wear amplitude modulation atomic force microscopy of DNA: Small amplitude small set-point imaging. _Appl.


Phys. Lett._ 103, 063702. https://doi.org/10.1063/1.4817906 (2013). Article  ADS  CAS  Google Scholar  * Lai, C. Y., Santos, S. & Chiesa, M. Systematic multidimensional quantification


of nanoscale systems from bimodal atomic force microscopy data. _ACS Nano_ 10, 6265–6272. https://doi.org/10.1021/acsnano.6b02455 (2016). Article  CAS  PubMed  Google Scholar  * Benaglia,


S., Amo, C. A. & Garcia, R. Fast, quantitative and high resolution mapping of viscoelastic properties with bimodal AFM. _Nanoscale_ 11, 15289–15297. https://doi.org/10.1039/c9nr04396a


(2019). Article  CAS  PubMed  Google Scholar  * Kawai, S. _et al._ Time-averaged cantilever deflection in dynamic force spectroscopy. _Phys. Rev. B_ 80, 085422.


https://doi.org/10.1103/PhysRevB.80.085422 (2009). Article  ADS  CAS  Google Scholar  * Yalcin, S. E., Legg, B. A., Yeşilbaş, M., Malvankar, N. S. & Boily, J.-F. Direct observation of


anisotropic growth of water films on minerals driven by defects and surface tension. _Sci. Adv._ 6, eaaz9708. https://doi.org/10.1126/sciadv.aaz9708 (2020). Article  ADS  CAS  PubMed  PubMed


Central  Google Scholar  * Abooalizadeh, Z., Sudak, L. J. & Egberts, P. Nanoscale spatial mapping of mechanical properties through dynamic atomic force microscopy. _Beilstein J.


Nanotechnol._ 10, 1332–1347. https://doi.org/10.3762/bjnano.10.132 (2019). Article  CAS  PubMed  PubMed Central  Google Scholar  * Kawai, S. & Kawakatsu, H. Surface-relaxation-induced


giant corrugation on graphite (0001). _Phys. Rev. B_ 79, 115440. https://doi.org/10.1103/PhysRevB.79.115440 (2009). Article  ADS  CAS  Google Scholar  * Eslami, B. & Damircheli, M.


Biharmonic versus bimodal AFM: Numerical and experimental study on soft matter. _J. Appl. Phys._ 126, 095301. https://doi.org/10.1063/1.5116794 (2019). Article  ADS  CAS  Google Scholar  *


Zhang, Z., Ryu, S., Ahn, Y. & Jang, J. Molecular features of hydration layers probed by atomic force microscopy. _Phys. Chem. Chem. Phys._ 20, 30492–30501.


https://doi.org/10.1039/C8CP06126B (2018). Article  CAS  PubMed  Google Scholar  * Seeholzer, T., Gretz, O., Giessibl, F. J. & Weymouth, A. J. A Fourier method for estimating potential


energy and lateral forces from frequency-modulation lateral force microscopy data. _New J. Phys._ 21, 083007. https://doi.org/10.1088/1367-2630/ab3353 (2019). Article  ADS  CAS  Google


Scholar  * Weymouth, A. J. Non-contact lateral force microscopy. _J. Phys. Condens. Matter_ 29, 323001. https://doi.org/10.1088/1361-648x/aa7984 (2017). Article  CAS  PubMed  Google Scholar


  * Michel, K. H. & Verberck, B. Theory of the elastic constants of graphite and graphene. _Phys. Status Solidi B Basic Solid State Phys._ 245, 2177–2180.


https://doi.org/10.1002/pssb.200879604 (2008). Article  ADS  CAS  Google Scholar  * Hertz, H. Ueber die Berührung fester elastischer Körper. _Journal für die reine und angewandte Mathematik_


156–171, 1882. https://doi.org/10.1515/crll.1882.92.156 (1882). Article  MATH  Google Scholar  * Blakslee, O. L., Proctor, D. G., Seldin, E. J., Spence, G. B. & Weng, T. Elastic


constants of compression-annealed pyrolytic graphite. _J. Appl. Phys._ 41, 3373–3382. https://doi.org/10.1063/1.1659428 (1970). Article  ADS  CAS  Google Scholar  * Holscher, H., Schwarz, U.


D., Zworner, O. & Wiesendanger, R. Consequences of the stick-slip movement for the scanning force microscopy imaging of graphite. _Phys. Rev. B_ 57, 2477–2481.


https://doi.org/10.1103/PhysRevB.57.2477 (1998). Article  ADS  CAS  Google Scholar  * Lee, H., Lee, N., Seo, Y., Eom, J. & Lee, S. Comparison of frictional forces on graphene and


graphite. _Nanotechnology_ 20, 325701. https://doi.org/10.1088/0957-4484/20/32/325701 (2009). Article  CAS  PubMed  Google Scholar  * Frank, I. W., Tanenbaum, D. M., Van der Zande, A. M.


& McEuen, P. L. Mechanical properties of suspended graphene sheets. _J. Vac. Sci. Technol. B_ 25, 2558–2561. https://doi.org/10.1116/1.2789446 (2007). Article  CAS  Google Scholar  *


Zandiatashbar, A. _et al._ Effect of defects on the intrinsic strength and stiffness of graphene. _Nat. Commun._ 5, 3186. https://doi.org/10.1038/ncomms4186 (2014). Article  ADS  CAS  PubMed


  Google Scholar  * Choi, J. S. _et al._ Correlation between micrometer-scale ripple alignment and atomic-scale crystallographic orientation of monolayer graphene. _Sci. Rep._ 4, 7263.


https://doi.org/10.1038/srep07263 (2014). Article  CAS  PubMed  PubMed Central  Google Scholar  * Choi Jin, S. _et al._ Friction anisotropy-driven domain imaging on exfoliated monolayer


graphene. _Science_ 333, 607–610. https://doi.org/10.1126/science.1207110 (2011). Article  ADS  CAS  PubMed  Google Scholar  * Gallagher, P. _et al._ Optical imaging and spectroscopic


characterization of self-assembled environmental adsorbates on graphene. _Nano Lett._ 18, 2603–2608. https://doi.org/10.1021/acs.nanolett.8b00348 (2018). Article  ADS  CAS  PubMed  Google


Scholar  * Gallagher, P. _et al._ Switchable friction enabled by nanoscale self-assembly on graphene. _Nat. Commun._ 7, 1–7. https://doi.org/10.1038/ncomms10745 (2016). Article  CAS  Google


Scholar  Download references ACKNOWLEDGEMENTS The authors thank the Deutsche Forschungsgemeinschaft (Project number 407750697) for financial support. FUNDING Open Access funding enabled and


organized by Projekt DEAL. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Physics of Surfaces, Institute of Materials Science, Technische Universität Darmstadt, Alarich-Weiss-Str. 2, 64287,


Darmstadt, Germany Anna L. Eichhorn & Christian Dietz Authors * Anna L. Eichhorn View author publications You can also search for this author inPubMed Google Scholar * Christian Dietz


View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS A.L.E. extended the AFM measurement method, performed the experiments, conducted data


analysis and wrote the manuscript. A.L.E. and C.D. designed the experiment. C.D. proposed the basic idea for the method, discussed the results, and revised the manuscript. CORRESPONDING


AUTHOR Correspondence to Christian Dietz. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION PUBLISHER'S NOTE Springer Nature


remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION. RIGHTS AND PERMISSIONS OPEN ACCESS


This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as


long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third


party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the


article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the


copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Eichhorn, A.L., Dietz, C.


Torsional and lateral eigenmode oscillations for atomic resolution imaging of HOPG in air under ambient conditions. _Sci Rep_ 12, 8981 (2022). https://doi.org/10.1038/s41598-022-13065-9


Download citation * Received: 02 February 2022 * Accepted: 12 May 2022 * Published: 28 May 2022 * DOI: https://doi.org/10.1038/s41598-022-13065-9 SHARE THIS ARTICLE Anyone you share the


following link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer


Nature SharedIt content-sharing initiative