Substrate-aware computational design of two-dimensional materials

Two-dimensional materials have attracted considerable attention due to their remarkable electronic, mechanical and optical properties, making them prime candidates for next-generation electronic and optoelectronic applications. Despite their widespread use in combination with substrates in practical applications, including the fabrication process and final device assembly, computational studies often neglect the effects of substrate interactions for simplicity. This study presents a novel method for predicting the atomic structure of 2D materials on arbitrary substrates by combining an evolutionary algorithm, a lattice-matching technique, an automated machine learning interatomic potentials training protocol, and the ab initio thermodynamics approach for predicting the possible conditions of experimental synthesis of the predicted 2D structures. Using the Mo-S system on a c-cut sapphire substrate as a case study, we reveal several new stable and metastable structures, including previously known 1H-MoS2 and newly found Pmma Mo3S2, P-1 Mo2S, P21m Mo5S3, and P4mm Mo4S, where the Mo4S structure is specifically stabilized by interaction with the substrate. Electronic band structure calculations of Mo3S2, Mo2S, Mo5S3, and Mo4S showed their metallic behavior, while phonon properties calculations indicated their dynamic stability and substrate-induced modulation of the phonon density of states. Finally, we use the ab initio thermodynamics approach to predict the synthesis conditions of the discovered structures in the parameter space of the commonly used CVD technique. These results provide insights into computational substrate engineering, allowing one to study the substrate effect on the thermodynamic and dynamical stability of 2D materials and to modulate their electronic and phonon properties for their future applications, as well as to provide guide maps for their experimental synthesis.