从一篇文献中得到了不用脚本的一种办法：在煤分子上下添加多层石墨，石墨要尽可能靠近煤分子，然后对所有石墨原子固定笛卡尔坐标，跑几何优化。完成后，移动上下层石墨使其靠近煤分子，再跑几何优化，如此反复至合适的结构。当然，这个只能在上下平面模拟压实。