Local hydration structures at the solid–liquid interface around boundary edges on heterostructures are key to an atomic-level understanding of various physical, chemical and biological processes. Recently, we succeeded in visualising atomic-scale three-dimensional hydration structures by using ultra-low noise frequency-modulation atomic force microscopy. However, the time-consuming three-dimensional-map measurements on uneven heterogeneous surfaces have not been achieved due to experimental difficulties, to the best of our knowledge. Here, we report the local hydration structures formed on a heterogeneously charged phyllosilicate surface using a recently established fast and nondestructive acquisition protocol. We discover intermediate regions formed at step edges of the charged surface. By combining with molecular dynamics simulations, we reveal that the distinct structural hydrations are hard to observe in these regions, unlike the charged surface regions, possibly due to the depletion of ions at the edges. Our methodology and findings could be crucial for the exploration of further functionalities.

There are currently no experimental techniques that combine atomic-resolution imaging with elemental sensitivity and chemical fingerprinting on single molecules. The advent of using molecular-modified tips in noncontact atomic force microscopy (nc-AFM) has made it possible to image (planar) molecules with atomic resolution. However, the mechanisms responsible for elemental contrast with passivated tips are not fully understood. Here, we investigate elemental contrast by carrying out both nc-AFM and Kelvin probe force microscopy (KPFM) experiments on epitaxial monolayer hexagonal boron nitride (hBN) on Ir(111). The hBN overlayer is inert, and the in-plane bonds connecting nearest-neighbor boron and nitrogen atoms possess strong covalent character and a bond length of only ∼1.45 Å. Nevertheless, constant-height maps of both the frequency shift Δf and the local contact potential difference exhibit striking sublattice asymmetry. We match the different atomic sites with the observed contrast by comparison with nc-AFM image simulations based on the density functional theory optimized hBN/Ir(111) geometry, which yields detailed information on the origin of the atomic-scale contrast.