The success of aberration correctors has led to a remarkable improvement in the spatial resolution in scanning transmission electron microscopy (STEM) to better than half an angstrom [1], which is sufficient for determining atomic structure in most materials. However, this high resolution applies only in the lateral two dimensions. The resolution (probe elongation) along the optical axis typically remains several nanometers, hindering the development of three-dimensional atomic-scale imaging in electron microscopy. There are several approaches towards achieving high-resolution 3D imaging and spectroscopy, including electron tomography, hollow-cone illumination STEM and confocal STEM. Though electron tomography has achieved excellent 3D atom-by-atom imaging in some cases [3], it is not well-suited to typical plate-shape thin specimens. Since the successful development of higher-order aberration correctors [2], it has become possible to increase the illumination angle (α) up to 70 mrad. Therefore, optical depth sectioning with large-angle illumination (LAI) STEM [3] may allow us to solve three-dimensional materials problems. However, two technical difficulties must be overcome to perform atomic-scale optical depth sectioning via LAI-STEM [4]. One difficulty is electron dose. For larger-angle illumination, the electron probe is very finely focused in all three dimensions and we may need to perform observations under relatively low electron dose (shot noise) to avoid excessive specimen damage. The other difficulty is chromatic aberration, where the electron probe is elongated along z-direction by the energy spread of the electrons. In this study, we explore single dopant visibility using LAI-STEM imaging under realistic dose and energy spread conditions through full dynamical image simulations. We also show application of LAI-STEM imaging to thickness measurement and surface imaging [5]. Figure 1 shows simulated ADF-STEM images of Ce-doped w-AlN, for the conditions of 300 kV electrons, α = 60 mrad, C c = 1 mm, and ΔE = 0.3 eV. No significant difference is visible between the images with and without chromatic aberration. Figure 2 shows, for a sample 7.4 nm thick, a cross-sectional-defocusview along the line X-X' for electron doses of 25×2 m e -/pix (m = 0, 1, … ,7). Though the chromatic aberration gives a slight elongation along the z-direction, the bright Z-contrast is well localized within ±2 unit cells along z-direction, comparable with the length of point spread function. Therefore, we may efficiently remove the channeling effect by using large-angle illumination. Z-contrast reduction by chromatic aberration is relatively small, enabling focal-series imaging at the low-dose condition of 200 e -/pix. However, for multiple dopants (Fig. 2(c) and (d)), the bright Z-contrast is delocalized along the zdirection, making it difficult not only to determine the dopant depths but even to reliably estimate the number of dopants. Therefore, it will be necessary to develop a new chromatic aberration corrector or better electr...