作两个图

第一性原理新手...非常非常新

1楼2010-09-27 21:32:26

ylli

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1558.8
红花: 3
帖子: 267
在线: 41.6小时
虫号: 277021
注册: 2006-09-03
专业: 凝聚态物性 II ：电子结构

作两个图

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指点 2010-09-27 22:42:17

逐步增加Ecut,观测总能随着Ecut增加的变化,直到总能基本不变,此时的Ecut为合理的Ecut.

此外,还要测试k点选取,方法类似. 即:绘制随着k点增加,总能的改变.

» 本帖已获得的红花（最新10朵）

dashuai1133

2楼2010-09-27 21:48:05

xiaobai216

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 227
散金: 38
帖子: 43
在线: 26.4小时
虫号: 1062398
注册: 2010-07-22
性别: MM
专业: 半导体材料

引用回帖:

Originally posted by ylli at 2010-09-27 21:48:05:
逐步增加Ecut,观测总能随着Ecut增加的变化,直到总能基本不变,此时的Ecut为合理的Ecut.

此外,还要测试k点选取,方法类似. 即:绘制随着k点增加,总能的改变.

那就是要建立一个固定模型，固定住其他参数改变Ecut来进行结构优化得到体系总能？我这么做过，但得到的结果很费解，随着E的增大，Total Energy 出现了一个峰值，这是怎么回事呢

赞一下

第一性原理新手...非常非常新

3楼2010-09-27 22:20:12

xiaobai216

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 227
散金: 38
帖子: 43
在线: 26.4小时
虫号: 1062398
注册: 2010-07-22
性别: MM
专业: 半导体材料

引用回帖:

再麻烦问您下，收敛测试是不是在结构优化的过程中逐步改变E？

赞一下

第一性原理新手...非常非常新

4楼2010-09-27 22:24:51

ylli

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1558.8
红花: 3
帖子: 267
在线: 41.6小时
虫号: 277021
注册: 2006-09-03
专业: 凝聚态物性 II ：电子结构

在结构优化中改变Ecut

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+1):谢谢回复！ 2010-09-28 08:16:45

有峰值有可能,但是应该还会下来,也就是最后的能量是最低的,不然可能你的结构不是局域最小(第一性原理中的结构优化不是全局优化).

5楼2010-09-28 07:44:03

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

1ST强帖: 5
应助: 46 (小学生)
贵宾: 8.842
金币: 7440.4
散金: 9654
红花: 118
沙发: 23
帖子: 5306
在线: 1961.4小时
虫号: 537452
注册: 2008-04-01
性别: GG
专业: 理论和计算化学
管辖: 仿真模拟

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-11-09 13:50:20

全局优化=分子动力学优化+第一性原理优化？

6楼2010-09-28 08:17:19

limin2007

铁杆木虫 (著名写手)

应助: 4 (幼儿园)
金币: 9402.1
散金: 38
红花: 8
帖子: 1437
在线: 953.2小时
虫号: 454730
注册: 2007-11-10
专业: 凝聚态物性I:结构、力学和

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

我增加截断能，总能量一直还在增加，还没有出现最小的数值，是不是还需要一直怎加截断能计算？你说的绘制图形，是不是每一次的截断能和对应的计算能量表示出来，在orgi里面绘制？

赞一下(1人)

7楼2010-11-08 12:03:38

zzy870720z

荣誉版主 (文坛精英)

1ST强帖: 3
应助: 47 (小学生)
贵宾: 9.05
金币: 30914.3
散金: 5613
红花: 68
沙发: 99
帖子: 12592
在线: 23567.6小时
虫号: 745488
注册: 2009-04-10
性别: GG
专业: 凝聚态物性I:结构、力学和
管辖: 分子模拟

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ym23(金币+2):感谢指导 2010-11-08 16:33:57

引用回帖:

Originally posted by limin2007 at 2010-11-08 12:03:38:

我增加截断能，总能量一直还在增加，还没有出现最小的数值，是不是还需要一直怎加截断能计算？你说的绘制图形，是不是每一次的截断能和对应的计算能量表示出来，在orgi里面绘制？

不用试最小值，而是基本不变，比如小数点后面四五位在变化等等，当然由你自己的判别标准就可以
增加Ecut不会有最小值的，理论上肯定是Ecut越大能量应该最小，计算就需要考虑计算时间和计算精度的平衡了
直观上肯定是图形方便判断，E-Ecut的能量曲线，当Ecut增加到一定时，E就基本不再变化了就可以了
另外至于用什么画图，我的感觉是测试的话基本都差不多，如果文章中需要的话可以考虑用origin

博学、审问、慎思、明辨、笃学

8楼2010-11-08 14:45:01

limin2007

铁杆木虫 (著名写手)

应助: 4 (幼儿园)
金币: 9402.1
散金: 38
红花: 8
帖子: 1437
在线: 953.2小时
虫号: 454730
注册: 2007-11-10
专业: 凝聚态物性I:结构、力学和

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by zzy870720z at 2010-11-08 14:45:01:

不用试最小值，而是基本不变，比如小数点后面四五位在变化等等，当然由你自己的判别标准就可以
增加Ecut不会有最小值的，理论上肯定是Ecut越大能量应该最小，计算就需要考虑计算时间和计算精度的平衡了
直观 ...

谢谢你，我选择一个截断能计算后，能量收敛图形也趋于一条直线了，但是我增加截断能后在计算，castep文件里面的总能量还是减小，按照您的意思增加截断能总能量是不变才最好是这样吗？
那我的问题是不是还需改变截断能继续计算？

[ Last edited by limin2007 on 2010-11-8 at 15:44 ]

赞一下(1人)

9楼2010-11-08 15:28:16

zzy870720z

荣誉版主 (文坛精英)

1ST强帖: 3
应助: 47 (小学生)
贵宾: 9.05
金币: 30914.3
散金: 5613
红花: 68
沙发: 99
帖子: 12592
在线: 23567.6小时
虫号: 745488
注册: 2009-04-10
性别: GG
专业: 凝聚态物性I:结构、力学和
管辖: 分子模拟

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+1):谢谢交流 2010-11-09 00:14:14

引用回帖:

Originally posted by limin2007 at 2010-11-08 15:28:16:

谢谢你，我选择一个截断能计算后，能量收敛图形也趋于一条直线了，但是我增加截断能后在计算，castep文件里面的总能量增大了，按照您的意思增加截断能总能量是不变才最好是这样吗？
那我的问题是不是还需改变 ...

看最后的能量，最初来Etot vs Ecut 曲线在判断比较好
总能量应该是随着Ecut的增加减小才对