输入/输出文件依赖于进行结构弛豫的外部代码和作业提交的类型。 我们的原理中一个重要的技术要素是多步结构弛豫策略,该方法有坚实的理论基础。为了提供依据 能量的正确的结构排列,最后一步的结构和能量必须是高质量的。大多数新产生的结构都距局部最 小值很远(例如:所包含的键不是太短就是太长)而它们的高质量弛豫的代价极其昂贵。如果第一 步的弛豫使用较为粗略的计算条件,只有最后一步的弛豫才需要高质量的计算,则这种高花费就可 以避免。第一步的结构弛豫可以使用代价较低的方式(例如:原子间的势使用GULP)。你可以改变 计算条件(基本设置,k点采样网格,赝势或者PAW势)或者近似的程度(原子间势、LDA 或GGA)或者结构弛豫的代码(GULP、DMACRYS、SIESTA、VASP、CP2K、QE 、CASTEP 或ATK)对于每个候选结构在结构弛豫期间。更进一步,为了保证结构弛豫的稳健性,也需要多步 策略——如果初始力和应力太大,全弛豫可能会导致结构的“爆炸”。因此,我们强烈建议初始优 化原子的位置时保证晶格参数为常数,或者优化晶胞的形状和原子位置时保证单位晶胞体积是常数。 当体积固定进行优化时,你就不必担心在平面波计算中的Pulay应力。因此,使用小一点的基本设置 也可以;然而,对于常压变晶胞弛豫,你需要一个高质量的基本设置。对于结构弛豫,你可以从一个 小一点的k点设置开始,但是不要忘记在最后一步结构弛豫中适当增加,以获得精确的能量。