This contribution explores the systematic substitution of phosphorene monoflakes (Mfs) and biflakes (Bfs) with aluminum, silicon, and sulfur. These systems were investigated using density functional theory employing the TPSS exchange-correlation functional and complete active space self-consistent field (CASSCF) calculations. Al and Si substitution produces significant structural changes in both Mfs and Bfs compared to S-substituted and pristine systems. However, in Mfs, all heteroatoms generate a decrease in band gap and the ionization potentials (IP), and an increase in electron affinity (EA) in comparison with pristine phosphorene. Al doping improves the hole mobility in the phosphorene monoflake, while Si and S substitutions exhibit a similar behavior on EAs and reorganization energies. For Bfs, the presence of Si-Si and Al-P interlaminar interactions causes structural changes and higher binding energies for Si-Bfs and Al-Bfs. Regarding the electronic properties of Bfs, substitution with Si does not produce significant variations in the band gap. Nevertheless, it conduces the formation of hole transport materials, which does not occur in Si-Mfs. The same is observed for Al systems, whereas no correlation was identified between the doping level and reorganization energies for S complexes. The substitution with Al and S leads to an opposite behavior of the band gap and IP values, while the EA variation is similar. In summary, the nature of heteroatom and the doping degree can modify the semiconductor character and electronic properties of phosphorene mono- and biflakes, whose trends are closely related to the atomic properties considered. Overall, these computational calculations provide significant insights into the study of doped phosphorene materials.
Keywords: Aluminum; DFT; Doped phosphorene nanoflakes; Silicon; Sulfur.