24小时热门版块排行榜    

查看: 1839  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡开发发

专家顾问 (著名写手)

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的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
9楼: Originally posted by Tonisam at 2014-12-18 09:29:18
补充一下,我做固定体积的几何优化是想得到E-V曲线,然后利用准谐德拜近似来算热力学近似。论坛上有人说,这个E-V曲线应该是采用加静水压力来改变,而不是直接改动晶格常数。我固体物理学的一般,也没仔细研究德拜 ...

可以这么做,不过对于Al这样的格子受到均匀并且不太大压力三个轴向的压力分量应该一样,形状应该还是立方的。如果是加压在特定的轴向,这种情况光靠恒定体积应该解决不了问题,需要进行静水压。

我没做过固体的力学性质方面的计算,仅供参考。
不一定挂在论坛,计算问题问题欢迎留言。
10楼2014-12-19 09:48:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

Tonisam

银虫 (小有名气)

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

zsl_321123

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
Tonisam(ljw4010代发): 金币+1, 谢谢交流! 2014-12-17 21:12:59
Tonisam: 金币+9, 有帮助 2014-12-22 17:27:17
你不优化晶格常数体积不就孤定了吗

[ 发自手机版 http://muchong.com/3g ]
宁愿骑着自行车笑,不愿坐着宝马哭。。。。。
6楼2014-12-17 19:21:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tonisam

银虫 (小有名气)

这也是个法子,但是总感觉没有固定体积的方法好。

[ 发自手机版 http://muchong.com/3g ]
7楼2014-12-18 09:23:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tonisam

银虫 (小有名气)

引用回帖:
5楼: Originally posted by 卡开发发 at 2014-12-17 14:03:34
Al晶格本身有对称性,可变化的自由度只有晶格常数,如果你固定了体积,还要保证原有的对称性,就没有可优化的自由度,程序算完能量自然只是输出单点能。正确的做法应该是先进行Make P1。

不过固定体积这种事情本 ...

太感谢了,我忘记考虑对称性了。

[ 发自手机版 http://muchong.com/3g ]
8楼2014-12-18 09:24:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tonisam

银虫 (小有名气)

引用回帖:
5楼: Originally posted by 卡开发发 at 2014-12-17 14:03:34
Al晶格本身有对称性,可变化的自由度只有晶格常数,如果你固定了体积,还要保证原有的对称性,就没有可优化的自由度,程序算完能量自然只是输出单点能。正确的做法应该是先进行Make P1。

不过固定体积这种事情本 ...

补充一下,我做固定体积的几何优化是想得到E-V曲线,然后利用准谐德拜近似来算热力学近似。论坛上有人说,这个E-V曲线应该是采用加静水压力来改变,而不是直接改动晶格常数。我固体物理学的一般,也没仔细研究德拜准谐近似的原理。也不确定这个方法是不是对的。

[ 发自手机版 http://muchong.com/3g ]
9楼2014-12-18 09:29:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Tonisam 的主题更新
信息提示
请填处理意见