A symmetry-oriented crystal structure prediction method for crystals with rigid bodies

We have developed an efficient crystal structure prediction (CSP) method for desired chemical compositions, specifically suited for compounds featuring recurring molecules or rigid bodies. We applied this method to two metal chalcogenides: Li₃PS₄ and Na₆Ge₂Se₆, treating PS₄ as a tetrahedral rigid body and Ge₂Se₆ as an ethane-like dimer rigid body. Initial trials not only identified the experimentally observed structures of these compounds but also uncovered several novel phases, including a new stannite-type Li₃PS₄ structure and a potential metastable structure for Na₆Ge₂Se₆ that exhibits significantly lower energy than the observed phase, as evaluated by density functional theory (DFT) calculations. We compared our results with those obtained using USPEX, a popular CSP package leveraging genetic algorithms. Both methods predicted the same lowest energy structures in both compounds. However, our method demonstrated better performance in predicting metastable structures. The method is implemented with Python code which is available at https://github.com/ColdSnaap/sgrcsp.git.