Black phosphorus (BP) has recently emerged as a promising semiconducting twodimensional material. However, its viability is threatened by its instability in ambient conditions, and by the significant decrease of its band gap in multilayers. We show that one could solve all the aforementioned problems by interfacing BP with hexagonal boron nitride (hBN). To this end, we simulate large, rotated hBN/BP interfaces using linear-scaling density functional theory (DFT). We predict that hBN-encapsulation preserves the main electronic properties of the BP monolayer, while hBN spacers can be used to counteract the band gap reduction in stacked BP. Finally, we propose a model for a tunneling field effect transistor (TFET) based on hBN-spaced BP bilayers.Such BP TFETs would sustain both low-power and fast-switching operations, including negative differential resistance behaviour with peak-to-valley ratios of the same order of magnitude as those encountered in transition metal dichalcogenide TFETs.1 Keywords 2D heterostructures; tunneling transistor; linear-scaling DFT; electric fields Since its inception, the two-dimensional (2D) material community has switched focus several times between different materials, aiming to satisfy different technological requirements. For instance, the lack of a band gap in graphene caused an increased interest in semiconducting transition metal dichalcogenides (TMDCs). While their many qualities have enabled their remarkable success in semiconductor electronics, 1,2 spintronics, 3 and optoelectronics, 2,4 the associated flaws of TMDCs are also of concern. Most notably, they only exhibit a direct band gap when in monolayer form, 5 the carrier mobilities are far smaller than those encountered in graphene, and their large band gap implies that the infrared (IR) spectrum cannot be harnessed by TMDC optoelectronics.In this context, the spotlight has recently shifted to layered black phosphorus (BP), 6,7 a material which rectifies the most important shortcomings of TMDCs. The band gap of BP is direct even in stacked forms, and its value changes significantly with the number of layers: from 0.3 eV in bulk to 1.5-2.0 eV for the monolayer. 8,9 This range of values allows BP to