Chemical modifications and stability of phosphorene with impurities: A first principles study

We perform a systematic first-principles study of phosphorene in the presence of typical monovalent (hydrogen, fluorine) and divalent (oxygen) impurities. The results of our modeling suggest a decomposition of phosphorene into weakly bonded one-dimensional (1D) chains upon single- and double-side hydrogenation and fluorination. In spite of a sizable quasiparticle band gap (2.29 eV), fully hydrogenated phosphorene found to be dynamically unstable. In contrast, full fluorination of phosphorene gives rise to a stable structure, being an indirect gap semiconductor with the band gap of 2.27 eV. We also show that fluorination of phosphorene from the gas phase is significantly more likely than hydrogenation due to the relatively low energy barrier for the dissociative adsorption of F2 (0.19 eV) compared to H2 (2.54 eV). At low concentrations, monovalent impurities tend to form regular atomic rows phosphorene, though such patterns do not seem to be easily achievable due to high migration barriers (1.09 and 2.81 eV for H2 and F2, respectively). Oxidation of phosphorene is shown to be a qualitatively different process. Particularly, we observe instability of phosphorene upon oxidation, leading to the formation of disordered amorphous-like structures at high concentrations of impurities.

💡 Research Summary

The paper presents a comprehensive first‑principles investigation of how typical monovalent (hydrogen, fluorine) and divalent (oxygen) impurities affect the structural, electronic, and dynamical properties of phosphorene. Using density‑functional theory (DFT) with the PBE functional and projector‑augmented wave potentials, the authors model single‑ and double‑sided adsorption of H, F, and O atoms on a phosphorene monolayer. They further refine the electronic band gaps with GW quasiparticle corrections and assess dynamical stability through phonon calculations performed with the Phonopy package. Transition‑state searches via the nudged elastic band (NEB) method provide activation energies for dissociative adsorption of H₂ and F₂ molecules and for surface diffusion of the adsorbed H and F atoms.

The key findings can be grouped into three impurity classes. First, both hydrogenation and fluorination break the original P–P network, leading to a reconstruction into weakly bonded one‑dimensional chains. Fully hydrogenated phosphorene exhibits a sizable quasiparticle gap of 2.29 eV, but its phonon spectrum contains imaginary modes, indicating a dynamical instability that would cause spontaneous decomposition at low temperature. In contrast, fully fluorinated phosphorene retains an indirect gap of 2.27 eV and shows no imaginary phonons, confirming that the fluorinated structure is mechanically robust. The difference stems from the stronger P–F covalency and the ability of fluorine to saturate dangling bonds without over‑straining the lattice.

Second, the thermodynamics and kinetics of impurity incorporation differ dramatically between H₂ and F₂. The calculated dissociative adsorption barrier for F₂ is only 0.19 eV, making fluorination from the gas phase highly favorable under mild conditions. By comparison, H₂ dissociation requires 2.54 eV, implying that hydrogenation would need either elevated temperatures or catalytic assistance. At low impurity concentrations, both H and F tend to align along regular atomic rows on the phosphorene surface, a pattern that could in principle enable anisotropic functionalization. However, the migration barriers for H and F atoms (1.09 eV and 2.81 eV, respectively) are sufficiently high to prevent facile diffusion, suggesting that achieving ordered rows experimentally would be challenging.

Third, oxidation follows a qualitatively distinct pathway. Oxygen atoms readily form P–O bonds, which disrupt the underlying P–P framework. At high oxygen coverages the phosphorene sheet collapses into a disordered, amorphous‑like network lacking long‑range order. This structural degradation is accompanied by a rapid closure of the electronic gap, potentially rendering the material metallic. The authors therefore conclude that oxygen is a destructive impurity for phosphorene, and its exposure must be minimized or mitigated by protective coatings.

From a practical standpoint, the study identifies full fluorination as the most promising route to chemically modify phosphorene while preserving its semiconducting character and mechanical integrity. Hydrogenation, despite producing a comparable band gap, is ruled out due to its intrinsic dynamical instability. Oxidation, on the other hand, is identified as a failure mode that must be avoided in device fabrication and handling. The authors propose that controlled fluorination could be employed to tune carrier concentrations, introduce anisotropic transport pathways, or engineer heterostructures with other two‑dimensional materials. Overall, the work provides a detailed theoretical framework for impurity‑driven engineering of phosphorene and offers clear guidelines for experimentalists seeking to exploit this emerging 2D semiconductor in electronic and optoelectronic applications.