We prove that scanning tunneling microscopy (STM) images of sub-surface Mn atoms in GaAs are formed by hybridization of the impurity state with intrinsic surface states. They cannot be interpreted in terms of bulk-impurity wavefunction imaging. High atomic resolution images obtained using a low-temperature apparatus are compared with advanced, parameter-free tight-binding simulations accounting for both the buckled (110) surface and vacuum electronic properties.PACS numbers: 71.15.Ap, 73.61.Eyi, The development of scanning tunneling microscopy has expended the applications of imaging to new areas of nanosciences and engineering such as atom or molecule identification, manipulation and nanostructuring [1,2,3]. Both spectroscopic and topological information is accessible, and quantum-size objects can be probed in real space with atomic resolution. In particular, there is growing interest in the STM images of quantum dots [4] and sub-surface impurity states in semiconductors [5,6,7,8]. Although the well-accepted Tersoff and Hamann's theory [9] simply relates the tunneling current to the local density of states (LDOS) at the tip position, the interpretation of specific experiments on localized states is still a matter of debate: on the one hand, the interaction between the surface and the quantum object under investigation must be examined; on the other hand, particularly in the case of deep bound-states, the carrier escape toward the extended band states may play a role [7]. In the last few years, acceptors in GaAs have attracted much attention because the shape of the images revealed strong and unexpected chemical signatures: from triangles in the case of shallow acceptors like carbon [7] and zinc [6], to asymmetric butterfly for the deep acceptor manganese [5]. So far, theoretical work has concentrated mostly on comparison of STM images with cross-sections of the bulk impurity wavefunction [5]. In this letter, we report on new experimental data obtained with a low-temperature apparatus and compare these results with advanced TB calculations. We first point that the STM image cannot reflect the LDOS of the impurity: indeed, empty-state STM images show only the group-III elements whereas the neutral acceptor wavefunction is distributed over the different chemical species. In fact, less than 40% of the bulk impurity LDOS shows up in the STM image. Then we prove that the image actually results from the hybridization between intrinsic surface states and the impurity state. A perturbative model is discussed and indicates that a quantitative relation between STM images and bulk wave function is not straightforward. Supercell calculations based on the sp 3 d 5 s * tight binding model and including the (110) surface as well as vacuum are performed and reproduce nicely the experimental images.The sample used in this study was grown by molecular beam epitaxy on a GaAs(001) substrate at 420 C. It consists of two 40 nm-thick GaAs layers doped with 2.1018 Mn cm −3 , embedded between 30 nm-thick GaAs conducting layers doped wi...