24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1351  |  回复: 15
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 xindushu 的 2 个金币 ,回帖就立即获得 2 个金币,每人有 1 次机会

xindushu

铜虫 (初入文坛)

[交流] 【求助】关于CASTEP结果文件里的能量

想知道结果文件里,开始的一段出现的每个原子的能量,如下,到底是什么意思呢?希望大家给点帮助啊~~

Pseudo atomic calculation performed for H 1s1

Converged in 15 iterations to a total energy of -12.1125 eV


Pseudo atomic calculation performed for O 2s2 2p4

Converged in 24 iterations to a total energy of -428.5316 eV


Pseudo atomic calculation performed for Al 3s2 3p1

Converged in 15 iterations to a total energy of -52.9882 eV
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水月慧心

金虫 (正式写手)

★ ★ ★
xindushu(金币+2):谢谢参与
ice_rain(金币+1):谢谢交流讨论 1-26 15:32
xindushu(金币+1): 1-27 08:40
对 H 1s1的价电子层进行赝势计算,在进行15次迭代后收敛,得到的能量值是-12.1125 eV,
对O 2s2 2p4的价电子层进行赝势计算,在进行24次迭代后收敛,得到的能量值是-428.5316 eV

对 Al 3s2 3p1的价电子层进行赝势计算,在进行15次迭代后收敛,得到的能量值是- -52.9882 eV
不知道对不对,还望高手指教更正,谢谢
2楼2010-01-26 09:35:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nyoxd

银虫 (小有名气)

★ ★ ★
xindushu(金币+2):谢谢参与
ice_rain(金币+1):谢谢交流 1-26 15:32
xindushu(金币+1): 1-27 08:40
本人觉得楼上说的有一定的道理,好像数据计算得不是很准确,H才-12.1125 eV
是否可以提高精度,采用较大基组计算
3楼2010-01-26 11:16:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyb4142

金虫 (正式写手)

★ ★ ★
xindushu(金币+2):谢谢参与
ice_rain(金币+1):xiexie 1-26 23:19
xindushu(金币+1): 1-27 08:40
2L正解。
3L为何说计算不准确呢,从哪来可以看出来?
就因为"H才-12.1125 eV“?这个貌似不能说明问题吧,精度在于使用的阶段能,收敛准则。。。
能量只有相对值,怎么看出12.1125不准确了?
4楼2010-01-26 21:57:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xindushu

铜虫 (初入文坛)

谢谢大家哦~~因为条件限制,我的计算精度提高的可能性很小。我算了一系列含上述原子的不同比例的体系,发现H、O、Al等都是这个能量~~那么这个能量大小有什么意义吗?和最后体系的如“BFGS: Final Enthalpy     = -1.93382236E+004 eV”有什么关系呢?
5楼2010-01-27 08:39:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyb4142

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★
ice_rain(金币+2):谢谢指导 1-27 19:24
aylayl08(金币+3):谢谢解释 1-27 19:25
xindushu(金币+3): 2010-01-28 08:19
a采用不同体系对于这个结果是不影响的。
这个结果给出的是,每种原子的能量,只要阶段能一样,不同体系计算结果不会有差别的(在精度之外的数字会有差别,这是数值计算的问题)
BFGS给出的是体系的总能量。
他们之间可以说没有直接的关系。
因为CASTEP头文件中的元素能量是指单原子的能力,和体系无关的。
6楼2010-01-27 19:09:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xindushu

铜虫 (初入文坛)

引用回帖:
Originally posted by wangyb4142 at 2010-01-27 19:09:12:
a采用不同体系对于这个结果是不影响的。
这个结果给出的是,每种原子的能量,只要阶段能一样,不同体系计算结果不会有差别的(在精度之外的数字会有差别,这是数值计算的问题)
BFGS给出的是体系的总能量。
他 ...

那么,体系的总能量,减去每种原子的n个能量,可以理解成结合能吗?
7楼2010-01-28 08:18:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jing15209791

金虫 (正式写手)

★ ★
xindushu(金币+2):谢谢参与
aylayl08(金币+1):谢谢提示 2010-01-28 14:12
提高精度,采用较大基组计算!好方法
8楼2010-01-28 09:32:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hgye

金虫 (小有名气)

xindushu(金币+2):谢谢参与
aylayl08(金币+1):谢谢提示 2010-01-28 19:06
xindushu(金币+1): 2010-01-30 18:54
这个数据好像是赝势文件里边提供的,与计算的体系及截断能都没有关系,但不知道它的作用,计算结合能需要重新计算单个原子的能量,不能直接用这个数据。
9楼2010-01-28 18:00:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyb4142

金虫 (正式写手)

★ ★
qasd(金币+2):谢谢 2010-01-30 17:38
xindushu(金币+3): 2010-01-30 18:54
引用回帖:
Originally posted by xindushu at 2010-01-28 08:18:57:



那么,体系的总能量,减去每种原子的n个能量,可以理解成结合能吗?

可以,这个在固体物理里面定义为内聚能(其实也是一种结合能,结合能不过是指一个过程的能量变化,可以根据不同初始-最终状态定义多个)
参见kittel的固体物理这个定义为cohesive energy ,有的固体物理书也称这个为结合能binding energy,其实只要指明你所用的公式就好
10楼2010-01-30 15:20:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xindushu 的主题更新
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复(可上传附件)
信息提示
请填处理意见