24小时热门版块排行榜    

查看: 1645  |  回复: 12

大瑶Alice

新虫 (小有名气)

引用回帖:
10楼: Originally posted by 锐利的碎片 at 2011-12-14 08:43:07:
金属可以只动半个电子,参考elk的基组文件 ,还要修改一下e0
'Cu+'                                      : spsymb
'copper with pseudo core-hole'             : spname
  -29.0000                       ...

前一段时间一直在研究手册来着。。。可是还是不知道您说的elk基组文件是什么文件?还有您说的修改eo,应该对应在哪个文件里呢?您的意思是就只把core内的电子占有数减少了1吗?别的还用改什么吗?新手见笑了,感激指教~
11楼2012-03-06 08:59:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大瑶Alice

新虫 (小有名气)

引用回帖:
4楼: Originally posted by stractor at 2011-12-04 20:56:47:
没使用过Win2K,但听过一个老外用Win2K计算TEM的报告,他指定了原子内壳层的电子数。这么说来,你的问题Win2K可以做到。

我计算了两个W原子的情况,把其中一个W原子的.inc中的第一行也就是1s轨道的occupy number由两个改成了一个,.inst中没有做相应的修改,.inm中-1,结果发现差了好多,能量从原来的-64664.50361174变成了-59524.62760543,这个原因是什么呢?
12楼2012-03-15 13:10:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大瑶Alice

新虫 (小有名气)

引用回帖:
2楼: Originally posted by nright at 2011-12-03 18:59:29:
没有做过这方面的计算。
不过,也许你可以试试更改输入文件case.inst文件,该文件给的是电子在各层的分布情况,你可以在特定电子层减掉一个电子,再计算。这样做以后,原子就带电了,也许还需要考虑更多的东西。 ...

我分别对含有两个W原子的晶胞其中的一个W原子做了如下改动:
1.W2.inc中第一行occupation数-1,;W2.inm中 -1 YES。没有对W2.inst做任何修改,这样的能量计算结果比原来的大了大约50000Ry。
2.W2.in2中的NE-1;.inm -1 YES。这样的总能量比原来小了1Ry
3.W2.inst中最后一行的occupation数-1;.inm -1 YES。这样的计算结果差不多与原来相同。
我觉得2和3的情况是一样的,只是精度的问题。我的k取得是1,shift,没有设定收敛条件,只是run_lapw。不知道我这样算精度有没有问题。也不知道1中不知道要不要改.inst 。至于能量的这种结果,道理还不是很清楚。
13楼2012-03-20 16:12:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 大瑶Alice 的主题更新
信息提示
请填处理意见