24小时热门版块排行榜    

CyRhmU.jpeg
查看: 792  |  回复: 8
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

csfn

荣誉版主 (知名作家)

优秀超版优秀版主优秀区长

[交流] 【讨论】关于量化软件中原子的能量定义的问题

最近在做报告的时候,有人提问,为什么用dmol算出来原子的能量会那么大,原子的能量是怎么定义的?dmol算出来,比如Ti原子的能量有-863.3838261 Ha,也就是23492.6739eV,这个能量是相当大的。量子力学里面,原子的能量应该是包含原子核的能量和电子的能量,考虑近似原子核几乎可以认为静止,能量就只要考虑运动的电子,电子包含动能项,势能项,以及考虑自旋的能量,那么这个考虑电子所带来的能量就不会很大,但是dmol算出来的单个原子的能量就很大~~  从dmol对总能定义来看,他是考虑了无相互作用粒子的动能,库仑作用产生的静电能以及多体贡献的交换关联能。这个能量为什么也会这么大了?又没有存在质量亏损?
在vasp里面能量是怎么定义的了?或者您现在用的或者熟悉的软件里面原子的能量是怎么定义的了?
请大家指点下~~


[ Last edited by csfn on 2008-1-10 at 09:30 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gkd_0000

新虫 (正式写手)

VASP里可以自己设定U值吗?castep好像不可以
成功地证实种种大胆的设想和揣测,这难道不正说明我们毕生的工作是无可非议和值得自豪的吗?--柯南道尔
7楼2007-12-22 14:33:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

acridine

木虫 (正式写手)

★ ★
zzgyb(金币+2,VIP+0):谢谢您的参与,欢迎您常来计算模拟版。
大概说一下
之前也没太注意过这个问题,看了楼主的帖子自己试了一下
用DMOL算的
首先计算了一个H氢原子,得到的能量大概为0.49 Ha.,而我们知道0.5 Ha 正好是氢原子核外电子轨道能量,因此Dmol算的应该是体系电子的总能量
之后有试了Ti原子,能量确实很大,但仔细观察结果文件会发现由于采用的是全电子计算,因此内层的电子的能量也考虑进去了,因此总能量会很低
如果选择赝势计算得到的总能量就会较小
2楼2007-12-08 12:04:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carlon

金虫 (小有名气)

★ ★
zzgyb(金币+2,VIP+0):谢谢您的参与,欢迎您常来计算模拟版。
能量=动能+势能+交换相关能
其中交换相关能主要用LDA和GGA等各种泛函(如PBE、PW91等)近似估计。
势能则包含电子-核、电子-电子的相互作用
dmol中计算的能量很大,是因为dmol用的是全电子的方法,
也就是体系中所有电子的能量都包含进去,所以体系的能量很大
而CASTEP中则是赝势的方法,不考虑内层电子的能量,所以计算的能量绝对值小很多。
dmol也可以不用全电子的方法计算,比如选择effective core potential处理core的话能量的绝对值应该会减小。
3楼2007-12-08 13:16:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)

谢谢楼上的讲解!!
我也想掺和一下
同样是采用赝势方法,而且渊源颇深的两个软件castep也vasp,怎么对同一个体系计算的总能也不一样呢?虽然只讨论总能没有多大意义,但这个区别究竟来自于何处?
望高人指点!!
4楼2007-12-08 19:46:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见