1) Operando KPFM: KPFM measurements were carried out in heterodyne frequency modulation (FM) mode with an external lock-in amplifier (Zurich Instruments HF2LI-MOD) to measure contact potential difference (CPD) between SFM tip and sample. More details about the KPFM working principle and interpretation is provided in Supplementary Note
We first performed KPFM measurements on the LLZO or the Li3PO4 surface at position x without any external potential applied to the Li-CE. x is the distance to the Li-CE. Upon applying a constant external current or potential to the symmetric Li|LLZO|Li cell, we measured a CPD change at a specific position x on LLZO surface. It is defined as:
Here, is the contact potential difference in the OCV state. is the contact potential difference with an applied potential . is work function change of the LLZO surface at position x, which changes with a material’s composition change.
2) Tr-EFM: Tr-EFM measurements were performed with the external lock-in amplifier (Zurich Instruments HF2LI-MOD) as well. Details are provided in Supplementary Note
Next, we brought the tip into contact with the surface and applied a DC bias voltage of −3 V for 1.5 s to induce ion displacements in LLZO. Then, we grounded the tip for 1.5 s to allow ion relaxation, after which we retracted the tip and repeated this sequence at the next position. During this sequence, we recorded the vibration amplitude of the cantilever at its second resonance frequency (ω) to track changes of the electrostatic force . As lithium ions are the main mobile charges in LLZO, we attribute changes of the electrostatic force to lithium-ion displacement. In solid ionic conductors, the ionic transport follows a stretched-exponential time response due to the electric field between the sample and tip60 (link)–63 (link):
Here, is the total amplitude change at frequency ω, is the amplitude change, before ionic relaxation due to ultrafast vibrational and electronic polarization44 (link),64 (link), is the amplitude change until the system reaches a saturation state due to ionic relaxation, τ is a time constant and β is a stretch exponent representing ion diffusion properties44 (link). For simplicity, we set β to 1 to fit the vs. time curve. Differences in ion diffusivity ( ) can be measured by fitting Eq. (
The comparison of tr-EFM results on LLZO and Au shows that tr-EFM can effectively track ion diffusivity in ionic conductors (Supplementary Fig.