|
|
在考察某种元素含量对金属熔体性质的影响时，采用了以下计算流程：
1）预先跑A模型在1673K下，通过GFN1-xTB采用NPT系综预平衡了16ps，初前几步SCF不收敛外，后续SCF、压力、温度、晶胞尺寸等均收敛；
2）然后采用A模型的restart文件通过删除同种元素原子数目与对应的速度建立6个模型，分别进行10ps续跑，获得构型；
3）目前6个模型中4个完成计算，SCF、压力、温度、晶胞尺寸等均收敛；
4）上述计算收敛后的构型另采用DZVP-MOLOPT-SR-GTH重新跑AIMD
说明：之所以采用GFN1-xTB，除了为了验证方法可行外，为了减少packmol建的模型中不合理接触
异常：模型6启动计算后，起始60余步SCF收敛正常，温度1700K左右；之后连续不收敛，温度基本在1800K左右（该模型的输入文件、ener文件和out文件见附件1，任务尚在运行中，故out并非正常结束）
6-GFN-xTB-NPT-t2.zip
(1.07 MB, 下载次数 Times of downloads: 0)
问题1：上述计算流程是否可行
问题2：模型6连续不收敛情况，但温度尚未发生巨变，是否继续让其运行，等待SCF再次收敛
问题3：可否抽取模型6运行过程中60余步的结构，重新新建AIMD任务进行续跑
问题4：AIMD过程中温度变化范围如果在3~5倍✖设定温度，是否可接受？如果温度超过3~5倍✖设定温度，或者骤升至10000K以上，是否应该马上暂停任务？
|
|
发表于 Post on 
收藏 Add to favorites