Towards accurate ionic relaxation algorithms for two-dimensional chalcogenide van der Waals materials - A first-principles study

• Published in: Physica. E, Low-dimensional systems & nanostructures Vol. 140; p. 115223
• Main Authors: el-Attar, Mahmoud, Allam, Nageh K.
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01.06.2022
• Subjects: Convergence; DFT; Interlayer spacing; Ionic relaxation; Polytypes; Transition metal dichalcogenides; Van der waals forces; Ionic relaxation; Van der waals forces; DFT; Interlayer spacing; Polytypes; Transition metal dichalcogenides; Convergence
• ISSN: 1386-9477; 1873-1759
• DOI: 10.1016/j.physe.2022.115223

Recently, two-dimensional (2D) layered materials, known as van der Waals (vdW) materials, have attracted remarkable attention due to their unique electrochemical and physical properties. Transition metal dichalcogenides (TMDs) are 2D materials, which commonly exist in three polymorphs (1T, 2H, and 3 R) in a honeycomb structure (hexagonal lattice). From a structure relaxation point of view, the 1T phase has less structural symmetry along the z-direction, rendering it more challenging than the 2H counterpart. Herein, the importance of enabling vdW forces as well as the associated challenges during the simulation of 1T polymorph are tackled. The two most widely used ionic relaxation methods (conjugate gradient and RMM-DIIS) were used for six TMDs (TaS2, TaSe2, MoS2, MoSe2, WS2, WSe2), and the results were compared to those reported in the literature. Moreover, a newly designed advanced conversion algorithm (TPSCA) is developed by incorporating shell scripts, DFT electronic relaxation, and Python coding. The developed algorithm is tested for 1T-TaS2 to confirm the obtained results and get better insights into the dependence of the energy on the different combinations of those parameters. [Display omitted] •A three-parameters simultaneous convergence algorithm is presented.•The algorithm enabled the consideration of vdW forces during DFT calculations.•The minimum energy convergence point can be directly identified using the algorithm.•A good agreement with the CG model and the reported experimental values.•Algorithm enables the investigation of 2D materials upon considering the vdW forces.