24小时热门版块排行榜    

查看: 1840  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Tonisam

银虫 (小有名气)

[求助] CASTEP固定体积优化遇到问题请教已有2人参与

本人最近在学习如何用CASTEP进行固定晶胞体积的优化,参照论坛中前人提出的修改.cell文件进行固定体积几何优化的方法,查看最后的任务结果中发现,直接跳过了几何优化的过程,仅仅进行能量计算完成之后就结束了。
以金属Al为例,具体操作过程如下:
1.导入金属Al单胞
2.设置几何优化计算参数(勾选-optimize cell)
3.点击files并保存
4.找到.cell文件后,用记事本打开,并在最后一行中加入 fix_vol : true
5. 直接运行.param文件
最后计算的结果是:
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -5.64703957E+001  0.00000000E+000                         3.24  <-- SCF
      1  -5.64715447E+001  7.78315099E+000   1.14900107E-003       3.42  <-- SCF
      2  -5.64715466E+001  7.78315066E+000   1.89855493E-006       3.56  <-- SCF
      3  -5.64715546E+001  7.78314770E+000   8.05779705E-006       3.68  <-- SCF
      4  -5.64715591E+001  7.78311413E+000   4.43893732E-006       3.81  <-- SCF
      5  -5.64715588E+001  7.78311411E+000  -2.55095320E-007       3.92  <-- SCF
      6  -5.64715588E+001  7.78311403E+000   3.67090164E-008       4.04  <-- SCF
------------------------------------------------------------------------ <-- SCF

Final energy, E             =  -56.46858492319     eV
Final free energy (E-TS)    =  -56.47155884207     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -56.47007188263     eV

For future reference: finite basis dEtot/dlog(Ecut) =      -0.080753eV
Total energy corrected for finite basis set =     -56.471572 eV

WARNING - there is nothing to optimise - skipping relaxation


根据warning的提示是,没有进行几何优化,跳过了弛豫过程就结束了。这不就相当于做了个能量计算么?有么有好心人,大神们帮忙分析一下原因。
ps:我用的是materials studio 7.0的版本,难道是版本问题不成???
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tonisam

银虫 (小有名气)

此贴别沉啊,求救!
4楼2014-12-17 11:27:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

tsglss

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
Tonisam: 金币+20, ★★★很有帮助, 可不可以说的详细一点,这样我把后面 所有的金币都给你 2014-12-16 11:20:20
ljw4010: 金币+1, 谢谢交流! 2014-12-17 21:12:31
此方法是错的!建议不要使用。

[ 发自手机版 http://muchong.com/3g ]
做好自己
2楼2014-12-15 22:17:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tonisam

银虫 (小有名气)

引用回帖:
2楼: Originally posted by tsglss at 2014-12-15 22:17:05
此方法是错的!建议不要使用。

那正确的方法是不是加压之后进行几何优化?
有人说现在castep可以加负压了
3楼2014-12-16 11:19:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
ljw4010: 金币+2, 谢谢交流! 2014-12-17 21:12:44
Al晶格本身有对称性,可变化的自由度只有晶格常数,如果你固定了体积,还要保证原有的对称性,就没有可优化的自由度,程序算完能量自然只是输出单点能。正确的做法应该是先进行Make P1。

不过固定体积这种事情本身就没什么物理意义,我猜测这个方法可能本身就是为了计算表面体系的a,b用的,即便这样,使用起来很可能力会很难收敛。
不一定挂在论坛,计算问题问题欢迎留言。
5楼2014-12-17 14:03:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见