Time-resolved force microscopy using delay-time modulation method
Abstract
We developed a time-resolved force microscopy technique by integrating atomic force microscopy using a tuning-fork-type cantilever with the delay time modulation method for optical pump-probe light. We successfully measured the dynamics of surface recombination and diffusion of photoexcited carriers in bulk WSe2, which is challenging owing to the effect of tunneling current in time-resolved scanning tunneling microscopy. The obtained results were comprehensively explained with the model based on the dipole-dipole interaction induced by photo illumination.
Research and development aiming to develop new functionalities by rapidly controlling quantum processes in materials and devices is flourishing, becoming one of the critical challenges. In fundamental properties such as phase transitions, the fluctuations in the dynamics of non-equilibrium systems at the nanoscale define their characteristics. Efforts are also being made in semiconductors to control regions to be as small as nm. To achieve these objectives, it is essential to properly understand and utilize the elementary processes in non-equilibrium states. For instance, it is indispensable to elucidate the ultrafast dynamics of carriers and phase changes induced by external fields with high temporal and spatial resolution.
The scanning tunneling microscopy (STM) system can obtain various information in real space with atomic-level spatial resolution, but its time resolution is constrained by factors such as circuitry. On the other hand, with the advancement in quantum optical techniques, pump-probe methods using ultrashort pulse lasers have achieved time-resolved measurements using pulse widths in the femtosecond and attosecond domains. However, typically, there are limitations in spatial resolution. Therefore, over the past decade, efforts have been made to combine these cutting-edge technologies, enabling the development of measurement and control techniques for local spectroscopy with optical time resolution while observing atomic structures and electronic states with atomic-level spatial resolution. Several time-resolved microscopy methods have been realizedTerada ; Yoshida_NatNanotech ; Yoshida_ACSPhoton ; Mogi_NPJ ; Arashida_ACSPhoton ; Wang ; Plankl ; Jelic ; Yoshioka ; Abdo (see Note 1 in Supplementary Information for more details concerning time-resolved STM).
Time-resolved STM is a highly promising method, but owing to the use of tunneling current in STM, the samples are limited to those with conductivity. Therefore, as a complementary measurement method, the development of time-resolved atomic force microscopy (AFM), which is one of the probe microscopy methods and measures the force between the probe tip and the sample, has been attemptedHamers ; Jahng_APL ; Schumacher_PNAS ; Araki ; Schumacher_APL ; Kuroiwa . For example, interactions based on surface photovoltage (SPV)Hamers and dipole induced in probe tip by light (PiFM)Rajapaksa ; Yamanishi , have been used as a time-resolved force probeHamers ; Jahng_APL ; Schumacher_PNAS ; Araki . Initially, methods synchronizing the cantilever vibration with pulse light excitation were used; however, they had limited time resolutionAraki . Subsequently, picosecond-range time resolution was achieved, but the S/N ratio was not particularly highSchumacher_APL . Other innovations include techniques that use the difference in modulation frequency between pump light and pulse light to remove noise effectsJahng_APL and self-detection tuning-fork-type cantilever driven by a quartz oscillator to eliminate the effects of light irradiation used for force measurementJahng . However, these methods all relied on the intensity modulation of the excitation light, similar to the macroscopic optical pump-probe method.
Here, we demonstrate a time-resolved force microscopy technique developed by integrating a tuning-fork-type frequency modulation (FM)-AFM with the delay time modulation method. We successfully measured the dynamics of surface recombination and diffusion of photoexcited carriers in bulk WSe2, which are challenging in time-resolved STM owing to the effect of tunneling current.
Figure 1(a) shows an overview of the system developed. In this study, we used an optical system with the delay time modulation of two pulse lasers (532 nm, 45 ps pulse width, 500 kHz) that can be driven by an external voltage pulse generatorIwaya ; Mogi_APEX1 . The delay time modulation shown in Fig. 1(b) is a method by which the delay time between the pump light and pulse light periodically changes (in this case at 20 Hz) between the measurement time and the sufficiently separated time , and the corresponding signal is detected by lock-in detection. With the delay time modulation method, it is possible to not only eliminate the effects of thermal expansion but also extract only the force component that reflects the dynamics of the sample excited by the pump lightTerada ; Mogi_NPJ ; Iwaya ; Mogi_APEX1 ; Yoshida_APEX .
As shown in Fig. 1, a tuning-fork-type cantilever was used, eliminating the need for an optical system for probe vibration detection. This choice makes the method implementable even within complex devices inside a UHV chamber. The spring constant can be high, allowing for a reduced vibration amplitude. The probe tip was attached at an angle to allow for light irradiation directly from the direction of the sample surface. By choosing the probe tip type, both STM and AFM measurements are possible. A tungsten (W) probe, electropolished to a length of mm, was attached to the cantilever using silver epoxy paste (The frequency resonance Q-factor was confirmed in Fig. 1(c)). Laser light was introduced through an objective lens (N.A. ) positioned right at the top (below) of the probe tip, allowing for light irradiation with a spot diameter of several m without interference with the probe structure (Fig. 1(d)). With a 550 nm short-pass filter, the laser light and the light from the lamp for illumination were combined. The position of the spot was observed with a CCD camera placed on the objective lens. For time-resolved force measurement, the amplitude of the frequency shift responding to the delay time modulation was determined as a time-resolved signal using a lock-in amplifier.

As the sample, we used WSe2, one of the layered materials (Fig. 1(e)). In previous study using time-resolved STM, photoexcited carrier dynamics, including photoinduced SPV, recombination, and diffusion, were investigatedYoshida_APEX ; Mogi_APEX1 . Considering that the results of our time-resolved AFM can be compared with the findings of the previous study, we conducted experiments on multilayered WSe2.
First, to investigate the band alignment and photoexcited carriers near the surface of the bulk WSe2 sample, we conducted force spectroscopy (- curves, : sample bias)Arai under both dark and illuminated conditions (Fig. 2(a)). We selected three representative conditions: (i) dark state, (ii) weak excitation: only illuminating with the continuous LED lamp used for microscopy (28 W; spot diameter, 200 m; wavelength, nm; namely, excitation density of pW/), and (iii) strong excitation: in addition to (ii), we illuminated with a pulse laser (1 mW; spot diameter, m; central wavelength, 532 nm; implying excitation density of ). The averaged results from 20 curves obtained for each condition are shown in Fig. 2(a). The feedback was conducted in the repulsive regime.
By fitting the obtained - curves with a quadratic function for each, we derived the contact potential difference (), under the three conditions: (i) V, (ii) V, and iii) V, indicating the existence of SPV. SPV is the difference in with and without illumination, producing the force between the probe tip and the sample. The quadratic curvature became pronounced as the light intensity increased, and at the peak tended to increase. The curvature depends on the distance between the probe tip and the sampleArai . These results suggest that, owing to the attractive interaction produced upon illumination, the probe tip-sample distance decreased to counteract this attraction owing to the feedback ( V and Hz), increasing the repulsive force at the peak.
The - curve obtained by simultaneously capturing the current between the probe tip and the sample while measuring the - curve is shown in Fig. 2(b). In the (i) dark and (ii) weak excitation states, almost no current flowed. In contrast, in the (iii) strong excitation state, we detected a current on the order of several hundred fA. This result corresponds to the observation on , where the distance between the probe tip and the sample decreased, thinning the tunneling barrier. The minor current observed within the bandgap region was not noted in previously recorded results of time-resolved STMYoshida_APEX . Namely, for the present sample, it is considered that defects or impurities, which can induce states within the gap, exist at the surfaceAddou .

From the results above, the band structure of the sample and probe tip, as considered in the dark state, is shown in Fig. 2(c). Bulk WSe2 exhibits -type behavior internally, with an electron affinity of eV and a band gap of eVYoshida_APEX ; Smyth . Near the surface, theoretical calculations and results of photoelectron spectroscopy have shown the presence of inherent surface band bendingLiu . Furthermore, Se vacancy defects on the surface have also been reported to generate donor-like statesKim . Consequently, as shown in Fig. 2(d), a very large SPV is generated, which is the origin of the variation shown in Fig. 2(a).
Next, we present the results of time-resolved measurements conducted at three different light intensities in Figs. 3(a)-3(c). The measurements were conducted under stable feedback conditions in the repulsive regime ( Hz, corresponding tip-sample force nN in this setup), with an amplitude of 4 nm, the of approximately 0 V, and the of less than 10 fA. The mid light intensity corresponds to that used for the V condition in Fig. 2(a). We successfully obtained a time-resolved signal with a very high S/N ratio as shown in Fig. 3(a). Fitting the data with a two-component exponential function () yielded ns and ns (Fig. 3(b)).
A schematic model explaining the observed results is shown in Figs. 3(d)-3(g). Previously, when measured by time-resolved STM, both fast and slow components were identifiedYoshida_APEX . The slow component on the 100 ns scale was considered to correspond to carrier diffusion in the bulk, because the carrier lifetime inside the crystal is extremely long on the s scale. On the other hand, the fast component was considered to correspond to the process of extracting electrons that accumulated on the surface as the tunneling current (not shown here). Indeed, for a current of 500 pA, a lifetime of ns was observed. As the current decreases (i.e., increasing the distance between the probe and the sample, thus extracting fewer electrons), the lifetime increases (250 pA leads to ns), providing a good explanation for the experimental findings.
The slow component in this study roughly matches the results obtained from STM measurements, and this component can be attributed to the diffusion process within the bulk (Fig. 3(g)). The fast component is related to the relaxation of photoexcited carriers on the surface (initial relaxation of SPV, Fig. 3(f)), similar to the case with time-resolved STM. However, the fast component in this experiment is about ns, slightly longer than that in STM measurements. We conducted experiments using the AFM system as a probe to exploit the force interactions under conditions where almost no tunneling current flows ( V, fA). Therefore, the relaxation of about ns after photoexcitation can be attributed to the recombination or capture of surface electrons via the gap states suggested by the - curve measurements. This process is considered the reason why the lifetime appears slightly longer than that in the STM experiment.
For the data obtained at the three excitation intensities, the results of fitting analysis using a two-component exponential function are shown in Figs. 3(b) and 3(c). These results are derived by averaging the fitting outcomes on both the positive and negative sides of the delay time. As shown in Fig. 3(b), the time constant did not change significantly with intensity, but the signal strength shown in Fig. 3(c) changed. Concerning the signal strength, when the light intensity was reduced (0.6 mW), the intensity ratio of the fast component to the slow one became particularly large ( (0.6 mW), (1.1 mW), and (2.3 mW)). This trend also matches the results of time-resolved STM measurement. In the case of STM, under the weak excitation condition, a small amount of photoexcited carriers were easily pulled out from the surface as tunneling current; therefore, a less diffusion component from the bulk was observed. On the other hand, as the excitation became stronger, a slow component became evident owing to the bulk-side carriers diffusing in the direction perpendicular to the surface after the faster recombination processes of photoexcited carriers via tunneling. The observation that the time constant did not change with light intensity in this study also supports the idea that the recombination processes via the surface states determine the decay rate. The time-resolved signals obtained in this study are consistent with the results obtained by time-resolved STM, indicating the successful measurement of the dynamics of photoexcited carriers using the system developed.

Owing to the use of delay time modulation, the typical forces that the probe experiences when light is irradiated onto a semiconductor sample are the force induced by SPV (Fig. 4(a)) and the image dipole force discussed for PiFM (Fig. 4(b)). The SPV-induced force is a coulomb force acting between the probe tip and the sample due to the photoexcited carriers that accumulate under the probe after photoexcitation (electrons in this case, Fig. 2(d)). Under the experimental conditions where the distance between the probe tip and the sample does not change with light irradiation (slow feedback condition), the change in SPV due to light irradiation corresponds to the change in , producing the force between the probe tip and the sample. As long as SPV exists under the probe tip, the force due to it will persist, and its time integral becomes the measured force signal.
On the other hand, according to the discussion on PiFM, the gradient and scattering forces operate owing to the interaction with the optical fieldJahng_PRB . When a dipole is induced on the probe tip side by the optical field, an image dipole is induced below the probe tip by assuming a semi-infinite sample shape as the experiments in Fig. 3(a)-(c). As a result, an attractive force operates between these dipolesJahng ; Jahng_PRB . Here, and are the complex effective polarizabilities of the tip and sample, respectively, within the point dipole approximation. is the incident light electric field vectors. When using pulsed light, a dipole is momentarily induced within the pulse width, and this force operates while plasmons exist on the probe tip side. The lifetime of plasmons in metals is very short on the fs scale compared with that within the laser pulse width (45 ps) used in this studyBrongersma , so it can be approximated that the force which reflects instantaneous only operates while the laser pulse is being applied.
Figures 4(c) and 4(d) show how the SPV-induced or image dipole force is detected in time-resolved measurement, respectively. Considering the ultrafast time domain, as shown in Figs. 3(d) to 3(e), there is a certain delay from when the carriers are excited until they accumulate on the surface, leading to corresponding change in the generating attractive interactions (rising time shown in Figs. 4(c) and 4(d), see Note 2 in Supplementary Information for more details). When the delay time is sufficiently long, the responses to the pump and probe light excitations are induced equally. Let us consider the case of irradiating the probe pulse while photoexcited carriers by the pump pulse are present by reducing . For simplicity, assuming the flat band condition (where the amount of photoexcited carriers on the surface is saturated, therefore SPV and also saturate for each light irradiation), the change induced by the probe light will increase to the same peak height as that caused by the pump light.
For the case of the force induced by SPV, as previously mentioned, since the force induced by SPV appears as a time integral (green area in Fig. 4(c)), the shorter the delay time, the weaker the measured attractive force becomes. This trend is similar to the signal detection principle of time-resolved STM when optical absorption bleaching is presentTerada ; Iwaya ; Mogi_APEX1 .
On the other hand, for the image dipole force shown in Fig. 4(b), when contemplating the force transient response in the high-speed regime, the force is only measured during light irradiation that is sufficiently shorter than the rising time of due to the carrier accumulation on the surface, as shown by the orange and red bars in Fig. 4(d). Furthermore, when using delay time modulation, since the value at longer delay times serves as a reference, the change in the magnitude of the force induced by the pump light indicated by the red bars in Fig. 4(d) is probed at , . Therefore, a larger attractive force is detected for shorter delay times. That is, the observed signal with the delay time modulation shows opposite polarity in case of the force derived from SPV shown in Fig. 4(c) and image dipole shown in Fig. 4(d).

Here, the latter mechanism corresponds well to the experimental results. To clarify this point, we measured responsive to the delay time modulation and the time variation of the phase (Fig. S1). The attractive force was stronger when the delay was shorter, indicating that we captured the time-resolved signal through the image dipole force. By utilizing the dipole that forms momentarily when the probe light is irradiated, we successfully measured the carrier dynamics of the sample excited by the pump light, such as recombination through intragap states and diffusion. Using a tuning fork and setting the measurement conditions to the repulsive region with a small amplitude of 4 nm, it is considered that the measurement of the near-field forces by PiFM became straightforward.
In conclusion, we developed a time-resolved force microscopy technique by integrating AFM using a tuning fork-type cantilever with the delay time modulation method for optical pump-probe light. We successfully measured the dynamics of surface recombination and diffusion of photoexcited carriers in bulk WSe2, which is challenging owing to the effect of tunneling current in time-resolved STM. By combining with multiprobe systems and using even shorter-pulse laser lightMogi_NPJ ; Mogi_APEX2 , we anticipate new developments such as time-resolved measurements of light-induced ultrafast non-equilibrium dynamics, including measurements in the operation conditions of nanostructures and device structures.
Acknowledgements.
We acknowledge the financial supports of a Grant-in-Aid for Scientific Research (20H00341, 22H00289, 22K14597, 23H00264,) from Japan Society for the Promotion of Science, Japan Core Research Evolutional Science and Technology (CREST) (JPMJCR1875), and JST PRESTO (JPMJPR22AA).Data Availability Statement
The authors declare no competing interests.
Supporting Information
The Supporting Information is available at DOI 10.35848/1882-0786/ad0c04.
References
- (1) Y. Terada, S. Yoshida, O. Takeuchi and H. Shigekawa, Nat. Photonics 4, 869 (2010).
- (2) S. Yoshida, Y. Aizawa, Z.-H. Wang, R. Oshima, Y. Mera, E. Matsuyama, H. Oigawa, O. Takeuchi and H. Shigekawa, Nat. Nanotechnol. 9, 588 (2014).
- (3) S. Yoshida, Y. Arashida, H. Hirori, T. Tachizaki, A. Taninaka, H. Ueno, O. Takeuchi and H. Shigekawa, ACS Photonics 8, 315 (2021).
- (4) H. Mogi, Y. Arashida, R. Kikuchi, R. Mizuno, J. Wakabayashi, N. Wada, Y. Miyata, A. Taninaka, S. Yoshida, O. Takeuchi and H. Shigekawa, npj 2D Mater. Appl. 6, 72 (2022).
- (5) Y. Arashida, H. Mogi, M. Ishikawa, I. Igarashi, A. Hatanaka, N. Umeda, J. Peng, S. Yoshida, O. Takeuchi and H. Shigekawa, ACS Photonics 9, 3156 (2022).
- (6) L. Wang, Y. Xia and W. Ho, Science 376, 401 (2022).
- (7) M. Plankl, P. E. Faria Junior, F. Mooshammer, T. Siday, M. Zizlsperger, F. Sandner, F. Schiegl, S. Maier, M. A. Huber, M. Gmitra, J. Fabian, J. L. Boland, T. L. Cocker and R. Huber, Nat. Photonics 15, 594 (2021).
- (8) V. Jelic, K. Iwaszczuk, P. H. Nguyen, C. Rathje, G. J. Hornig, H. M. Sharum, J. R. Hoffman, M. R. Freeman and F. A. Hegmann, Nat. Phys. 13, 591 (2017).
- (9) K. Yoshioka, I. Katayama, Y. Minami, M. Kitajima, S. Yoshida, H. Shigekawa and J. Takeda, Nat. Photonics 10, 762 (2016).
- (10) M. Abdo, S. Sheng, S. Rolf-Pissarczyk, L. Arnhold, J. A. J. Burgess, M. Isobe, L. Malavolti and S. Loth, ACS Photonics 8, 702 (2021).
- (11) R. J. Hamers and D. G. Cahill, Appl. Phys. Lett. 57, 2031 (1990).
- (12) J. Jahng, J. Brocious, D. A. Fishman, S. Yampolsky, D. Nowak, F. Huang, V. A. Apkarian, H. K. Wickramasinghe and E. O. Potma, Appl. Phys. Lett. 106, 083113 (2015).
- (13) Z. Schumacher, R. Rejali, R. Pachlatko, A. Spielhofer, P. Nagler, Y. Miyahara, D. G. Cooke and P. Grütter, Proc. Natl. Acad. Sci. U. S. A. 117, 19773 (2020).
- (14) K. Araki, Y. Ie, Y. Aso, H. Ohoyama and T. Matsumoto, Commun. Phys. 2, 10 (2019).
- (15) Z. Schumacher, A. Spielhofer, Y. Miyahara and P. Grutter, Appl. Phys. Lett. 110, 053111 (2017).
- (16) T. Kuroiwa and T. Takahashi, Jpn. J. Appl. Phys. 61, SL1004 (2022).
- (17) I. Rajapaksa, K. Uenal and H. K. Wickramasinghe, Appl. Phys. Lett. 97, 3 (2010).
- (18) J. Yamanishi, H. Yamane, Y. Naitoh, Y. J. Li, N. Yokoshi, T. Kameyama, S. Koyama, T. Torimoto, H. Ishihara and Y. Sugawara, Nat. Commun. 12, 3865 (2021).
- (19) J. Jahng, H. Kwon and E. Lee, Sensors 19, 1530 (2019).
- (20) K. Iwaya, M. Yokota, H. Hanada, H. Mogi, S. Yoshida, O. Takeuchi, Y. Miyatake and H. Shigekawa, Sci. Rep. 13, 818 (2023).
- (21) H. Mogi, Z. Wang, R. Kikuchi, C. Hyun Yoon, S. Yoshida, O. Takeuchi and H. Shigekawa, Appl. Phys. Express 12, 025005 (2019).
- (22) S. Yoshida, Y. Terada, M. Yokota, O. Takeuchi, Y. Mera and H. Shigekawa, Appl. Phys. Express 6, 016601 (2013).
- (23) T. Arai and M. Tomitori, Phys. Rev. Lett. 93, 256101 (2004).
- (24) R. Addou and R. M. Wallace, ACS Appl. Mater. Interfaces 8, 26400 (2016).
- (25) C. M. Smyth, R. Addou, S. McDonnell, C. L. Hinkle and R. M. Wallace, 2D Mater. 4, 025084 (2017).
- (26) R.-Y. Liu, K. Ozawa, N. Terashima, Y. Natsui, B. Feng, S. Ito, W.-C. Chen, C.-M. Cheng, S. Yamamoto, H. Kato, T.-C. Chiang and I. Matsuda, Appl. Phys. Lett. 112, 211603 (2018).
- (27) J. Kim, H. Park, S. Yoo, Y. Im, K. Kang and J. Kim, Adv. Mater. Interfaces 8, 2100718 (2021).
- (28) J. Jahng, J. Brocious, D. A. Fishman, F. Huang, X. Li, V. A. Tamma, H. K. Wickramasinghe and E. O. Potma, Phys. Rev. B 90, 155417 (2014).
- (29) M. L. Brongersma, N. J. Halas and P. Nordlander, Nat. Nanotechnol. 10, 25 (2015).
- (30) H. Mogi, Z. Wang, T. Bamba, Y. Takaguchi, T. Endo, S. Yoshida, A. Taninaka, H. Oigawa, Y. Miyata, O. Takeuchi and H. Shigekawa, Appl. Phys. Express 12, 045002 (2019).