Second harmonic generation (SHG) from three-dimensional topological insulators originates from both surface and bulk, which does not allow probing of surface states unless the measurement can separate the two contributions. In this study, we used combined measurements of transmitted and reflected SHG from epitaxially grown Bi 2 Se 3 thin films of different thickness on BaF 2 , and a bulk Bi 2 Se 3 crystal, to deduce surface and bulk nonlinear susceptibilities of Bi 2 Se 3 separately. We found that the surface * E-mail: yrshen@berkeley.edu, wtliu@fudan.edu.cn 2 contribution to SHG was comparable to that from the bulk of the crystal, but becomes dominant in ultrathin films. In the latter case, contributions from both air/Bi 2 Se 3 and Bi 2 Se 3 /BaF 2 interfaces were significant, and exhibited a strong out-of-plane polar ordering. The bulk contribution came mainly from the space charge region (SCR), which was formed by Se vacancies aggregated at the air/Bi 2 Se 3 interface; its magnitude can provide an estimate on the field strength in the SCR. Clarification of surface and bulk contributions to SHG can help nonlinear optical techniques be used as a versatile in situ probe for topological insulators.3