Abstract
Phosphorene monolayers exhibit a range of advanced functional properties; however, their spin configurations and structural phase transitions remain unexplored in the context of enhancing performance in quantum computation and information. In this study, we employ a two-dimensional (2D) Ising model with nearest- and next-nearest-neighbor interactions, denoted as [Formula: see text] and [Formula: see text] respectively, alongside a Monte Carlo approach, to investigate the structural phase transitions and spin behavior of phosphorene. Using this model, we derive the temperature-dependent phase diagram of the phosphorene in terms of the interaction strength ratio ([Formula: see text]), revealing three distinct phases—checkerboard-ordered, glassy, and quasi-ordered—with a phase boundary at [Formula: see text]. Additionally, we analyze the effects of possible defects and strains in the phosphorene lattice induced by the substrate on the phase transitions. Our findings have potential applications in adiabatically switching states through gradual variations in coupling constants, offering new possibilities for adiabatic quantum computing (AQC) and π-ring arrays. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1038/s41598-025-21704-0.