“…Most strikingly, intrinsic even-order THz harmonics, which are only generated in systems that break inversion symmetry, have never been demonstrated in an equilibrium material, having been observed only in superconductors with a net propagating supercurrent [46,47] or in carefully engineered devices [48]. The lack of THz second harmonic generation (SHG) in the studies of the bismuth chalcogenides [40,44,45] is of particular note, as SHG-and even-order HG in general-originating from the surface has been a well-established feature of the nonlinear optical response outside of the THz regime [49][50][51][52][53][54][55]. Furthermore, as the bismuth chalcogenides are prototypical topological insulators [56][57][58][59] with a centrosymmetric bulk and inversion symmetry-breaking surfaces, the second-order optical response of THz-SHG offers a pathway to measuring the properties of the topological surface state while intrinsically avoiding the properties of the bulk band, without resorting to doping [45].…”