||
关注:
1) 为提高结构搜索的效率及成功概率,哪些关键参数需尽可能准确设置
2) 使用过程中遇到的问题及解决策略
摘录:
叶 2:08:01 PM
CALYPSO4.0 程序预测和优化一切正常(优化过程中保留了OUTCAR_*),但pso_sor_*文件中没有焓值,请问是怎么回事?
王2:11:29 PM
文件是空的吗?
叶 2:12:36 PM
不是空的,OUTCAR正常优化了哈
王 2:12:52 PM
文件里什么内容
叶 2:14:55 PM
这是手动提取的焓值:
压力 体积 能量 OUTCAR_*
0.09 82.07 -18.23328041 OUTCAR_62
0.11 72.88 -16.76458672 OUTCAR_20
0.15 116.14 -25.91471775 OUTCAR_42
0.15 134.12 -30.72169013 OUTCAR_52
0.17 112.53 -25.79816780 OUTCAR_38
0.17 68.66 -18.24874537 OUTCAR_26
0.19 117.67 -23.66091375 OUTCAR_40
0.22 140.92 -28.54269814 OUTCAR_32
0.41 111.08 -26.45216629 OUTCAR_16
0.69 88.98 -12.94127388 OUTCAR_48
叶 2:15:14 PM
进行的是常压预测和优化计算
叶 2:18:14 PM
有点诡异。难道是焓值的源程序有问题?这是从getenth.py程序中摘录的部分内容
def enthalpy(caldir):
if pressure(caldir):
e = outEnergy(caldir)
else:
e = 610612509
# h = hardness()
return e
王2:20:44 PM
常压PSTRESS=0
王2:20:55 PM
INCAR是这么设置的吗?
叶2:21:57 PM
grep PST INCAR_*
INCAR_1:PSTRESS = 0.0
INCAR_2:PSTRESS = 0.0
INCAR_3:PSTRESS = 0.0
王 2:22:20 PM
你把PSTRESS=0.00001
王 2:22:32 PM
就可以了
叶 2:23:08 PM
啊,为什么
王 2:23:30 PM
程序实际提取的是enthalpy
叶 2:23:45 PM
这个小细节很难注意到哈
王2:23:54 PM
但当把PSTESS=0的时候没有enthalpy,只有能量
王2:24:30 PM
当把这个设的很小的时候也不会有什么影响
叶
嗯,我试试,能改动这个程序,避免出现这种情况吗?
def enthalpy(caldir):
if pressure(caldir):
e = outEnergy(caldir)
else:
e = 610612509
# h = hardness()
return e
叶 2:28:23 PM
4.0版CALYPSO是不是不需要编译哈
叶 2:28:33 PM
我直接用的下载的可执行文件
王2:32:00 PM
是的
探索铈在沉积过程中可能存在的相变规律;
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 22:22
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社