应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】如何在siesta中设置ghost state原子的计算
51/1返回列表
查看: 2937  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hw45888792

金虫 (正式写手)

[交流] 【求助】如何在siesta中设置ghost state原子的计算

用siesta计算体相金刚石中一个空位，形成能好像不准，单胞越大，形成能也越大，没有合理的解释。
师兄说，siesta用原子轨道基组，可能存在鬼态原子，需要在空位处设置一个假想的原子，我查了一下手册，没有完全搞清楚如何设置？
不知道哪位大侠有做过？
回复此楼

» 收录本帖的淘帖专辑推荐

siesta问题 SIESTA siesta基本概念

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-07 21:27:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hw45888792

金虫 (正式写手)
引用回帖:
Originally posted by zhangguangping at 2011-03-08 11:39:57:
其实，我之前从来没有做过ghost state的计算，并且对ghost state的知识和用处不了解。仅仅是我之前看手册的时候记得手册中提到了ghost atom的处理。但是具体怎么处理，我是没往心里记的。昨天晚上临回宿舍前，上 ...

这里我还有一个疑问，我的师兄说，由于轨道基组的原因，在用siesta计算从体相超单胞一个空位处拿掉的单个C原子能量的时候，需要在相同超单胞里面的相同位置放置一个C原子，然后其他地方全部放鬼原子，从而得到这个C原子的里面。
    还有一个方法，就在这个超单胞里面直接放一个C原子，直接计算。
    平面波方法不存在这个问题，VASP，可以通过模拟最小bulk单胞计算单个C原子的能量。

    我测试了一下，第一个方法得到的单原子能量最低，但是目前没有得到其它相关数据的支持。
赞一下 回复此楼
8楼2011-03-08 13:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

zhangguangping

木虫 (著名写手)
★ ★
zzy870720z(金币+2): 谢谢指教 2011-03-07 22:20:35
hw45888792(金币+4): 好厉害呀，好像回答过我好多问题啊~ 2011-03-08 10:18:18
引用回帖:
Originally posted by hw45888792 at 2011-03-07 14:27:16:
用siesta计算体相金刚石中一个空位，形成能好像不准，单胞越大，形成能也越大，没有合理的解释。
师兄说，siesta用原子轨道基组，可能存在鬼态原子，需要在空位处设置一个假想的原子，我查了一下手册，没有完全搞 ...

参照手册上的ChemicalSpeciesLabel的说明，对于gost原子的原子序数用一个负值就行，然后最好给他起一个别的元素名称，赝势也相应的改变一下。这个时候他的原子质量就会是1.0d30 a.u.如果你计算的是金刚石，那么gost原子的赝势就用C的就行，和你就算中用的其他的C用一个赝势就行，只不过最好起一个别的名字，然后对赝势也重命名一下，将其元素序号用一个负值比如-6表示。
赞一下(1人) 回复此楼
2楼2011-03-07 22:14:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hw45888792

金虫 (正式写手)
引用回帖:
Originally posted by zhangguangping at 2011-03-07 22:14:00:
参照手册上的ChemicalSpeciesLabel的说明，对于gost原子的原子序数用一个负值就行，然后最好给他起一个别的元素名称，赝势也相应的改变一下。这个时候他的原子质量就会是1.0d30 a.u.如果你计算的是金刚石，那么 ...

你好，这是我测试NV在金刚石2*2*2超单胞C62N中的测试输入数据：

计算鬼态为Cg，在搭建结构的时候，我在空位处加上一个C原子，写成输入文件以后，将其改成Cg，标号为3。其中，C标号为1，N标号为2。赝势方面，我将C.psf直接复制成Cg.psf，然后将Cg.psf里面的元素名称改成Cg，不知是否可行？

SystemName      C62N_ghost
SystemLabel     C62N_ghost
NumberOfSpecies         3

%block ChemicalSpeciesLabel
  1   6   C
  2   7   N
  3  -6   Cg
%endblock ChemicalSpeciesLabel


如果是常规计算，就不用鬼态，输入文件如下：
SystemName      C62N
SystemLabel     C62N
NumberOfSpecies         2

%block ChemicalSpeciesLabel
  1   6   C
  2   7   N
%endblock ChemicalSpeciesLabel

不知道是否正确？

计算过程没有任何问题，对比结果我也大概写在下面：

加入鬼态：
siesta: Program's energy decomposition (eV):
siesta: Eions   =     16680.576136
siesta: Ena     =      3653.202863
siesta: Ekin    =      7163.588024
siesta: Enl     =      -820.400928
siesta: DEna    =      -182.458272
siesta: DUscf   =        15.575306
siesta: DUext   =         0.000000
siesta: Exc     =     -3054.613637
siesta: eta*DQ  =         0.000000
siesta: Emadel  =         0.000000
siesta: Ekinion =         0.000000
siesta: Eharris =     -9905.679557
siesta: Etot    =     -9905.682781
siesta: FreeEng =     -9905.682783

siesta: Final energy (eV):
siesta:       Kinetic =    7163.588024
siesta:       Hartree =     894.608032
siesta:    Ext. field =       0.000000
siesta:   Exch.-corr. =   -3054.613637
siesta:  Ion-electron =   -8393.758885
siesta:       Ion-ion =   -6515.506314
siesta:       Ekinion =       0.000000
siesta:         Total =   -9905.682781


不加鬼态
siesta: Program's energy decomposition (eV):
siesta: Eions   =     16680.576136
siesta: Ena     =      3653.202863
siesta: Ekin    =      7164.385562
siesta: Enl     =      -820.554764
siesta: DEna    =      -182.943454
siesta: DUscf   =        15.589130
siesta: DUext   =         0.000000
siesta: Exc     =     -3054.694636
siesta: eta*DQ  =         0.000000
siesta: Emadel  =         0.000000
siesta: Ekinion =         0.000000
siesta: Eharris =     -9905.587400
siesta: Etot    =     -9905.591435
siesta: FreeEng =     -9905.591437

siesta: Final energy (eV):
siesta:       Kinetic =    7164.385562
siesta:       Hartree =     894.861106
siesta:    Ext. field =       0.000000
siesta:   Exch.-corr. =   -3054.694636
siesta:  Ion-electron =   -8394.637153
siesta:       Ion-ion =   -6515.506314
siesta:       Ekinion =       0.000000
siesta:         Total =   -9905.591435


计算结果显示，加入鬼态后，总能量略低0.1eV左右，而且鬼态原子附近有大约0.2单自旋电荷分布，其它C原子是2个电荷分布。
赞一下 回复此楼
3楼2011-03-08 11:11:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+5): 谢谢zhang兄指教，说的太好了 2011-03-08 15:56:38
引用回帖:
Originally posted by hw45888792 at 2011-03-07 14:27:16:
用siesta计算体相金刚石中一个空位，形成能好像不准，单胞越大，形成能也越大，没有合理的解释。
师兄说，siesta用原子轨道基组，可能存在鬼态原子，需要在空位处设置一个假想的原子，我查了一下手册，没有完全搞 ...

其实说句实话，我是真不厉害，很菜的。不是谦虚，我是自我表白。我在论坛里面比较活跃不是显摆我知道的多，而是恰是我知道少，但是身边有没有人用这个软件。所以只能依靠emuch论坛的这个版块。通过讨论可以相互交流一下，哪怕我的是错误的，我也可以说一下，这样高手看见了，可以来指正我，这也是一个自我的进步。另外咱们这个版块不是很活跃。貌似高手现在都不大来了，或者来了之后看到问题的幼稚就嗤之一笑，然后不懒的回复了。希望老虫还是多来指点指点我们这些幼虫们。

另外你的这个问题实际上可以通过多认真的读一读手册是可以得到解决的。使用siesta工作的是手边要么有一本打印的手册。要么在电脑上打开一份电子的手册。有什么问题，直接搜索一下相关的关键字。比如你的问题：直接在手册里面用ctrl+F搜索ghost,立马会定位到ChemicalSpeciesLabel位置，再搜，会定位到AtomicMass位置，接下来在PS.KBprojectors会发现对ghost设置的一些条件，这个需要检查一下产生的PS.KBprojectors对ghost atom是不是只有一个，如果是的话，那就是程序给你处理好了。如果不是的话，计算出来预料不到结果，需要考虑这个问题。PAO.Basis在这儿，解释了程序对ghost atom的基矢的处理办法。
到此为止，应该对siesta处理ghost atom有一个了解。如果不是量化高手的话，可能对手册上的解释还是知道如何下手，这个时候需要google了。在google里面输入ghost atom in siesta,一搜，google大多数时候不会辜负你的。
http://www.google.com.hk/search? ... i=&aql=&oq=
这个头两条就是siesta前辈们讨论关于怎么设置ghost。
很懊恼的是打不开，这个时候我来告诉你，网页代理来帮助：http://www.7daili.com/，首先打开这个网页代理，然后打开你要打开但是打不的那个网页，从浏览器上把这个网址复制到你的网页代理里面，然后让它给你打开。不要嫌弃这个太慢，耐心点，能打开总比打不开好多了。灯打开之后，看到里面讨论的时候，欣喜不由从心底产生。在之内的链接可以直接打开了（其实还是网页代理给你打开的）。不管你从这些讨论中得没得到你要的结果，但是起码，你学习了，搜索了，看到之前也有人讨论这个问题了。这就够了。如果问题解决了，皆大欢喜，如果还是没有解决，在走下一条通往罗马的大道。总之，这个过程是一个需要的过程，我认为也是发现问题，解决问题的作为学生的过程。可能你会说，作为学生，时间短，论文压力大。没有时间花费在这个上面，直接问问知道就行。前提是如果有知道在身边什么都好说。如果没有呢？难道咱们的问题不不解决了？向小木虫丢出来一个求助，不一定有人立马回复你。这个期间还得自己去走上面的过程，只不过小木虫的讨论给上面的过程多了一个希望而已。

这些只是我的感想而已。我也是初学者。来到小木虫，首先，要搜索一下关于你的这个话题的帖子。最起码，来到小木虫定居之后，要把之前讨论过的有关你用的这个软件的所有帖子全部看一遍，知道别人都讨论了些什么。结果怎么样，哪怕自己用不到，说不定以后用的时候，就需要回来再次检索。平时没事，一定要天天登录一下小木虫看一下自己用的这个程序有什么新的问题，别人是怎么解决的。这样对该软件在小木虫的历史就有了了解。对历史不懂得人，会走很多弯路的。比如最近yashumimi发的那个帖子transista的错误，其实，别人之前早就讨论过了。只是她没有关心，以为自己不会碰到这个问题，压根就没有看这个帖子。但是她就应为自己对这个帖子的默然，也像mail-list求助了，也发帖子了。我也看到一段时间没有人回复。估计期间她也费了好大劲，做了不少计算，并且花了不少机时。这个都是对历史不解或者上小木虫之后没有首先检索惹的祸。

最后，小木虫不是一个有问题就往上丢，没问题的时候压根不来的地方，而是一个同行之间相互了解相互交流的平台，窗口。就像一些学术会议一样。只不过小木虫上天天再开“学术会议”，“作报告的人”都是一些资源的提供者，更或者对我们这些板块，“作报告的人”是一些遇到问题的人，报告的内容就是问题本身而已。

当自己不是“作报告的人”时候要常来当听报告的。因为这个问题被别人讨论之后，再重复的提出来，可能会使得一些高手失去回答的耐心。不但对自己没有好处，然高手失去了对咱们新一代新手的信心，对别的新手来说，是一个莫大的损失。

我也看到有些虫友都工作了，都当教授了，还是常常忙里偷闲来答复一下新虫的问题。我看了，确实有些感动。

最后，祝愿小木虫SIESTA板块，越来越红火，能成就一些新手。希望老虫多来哺养一下我们这些新虫。


末了，言辞不妥之处，望新老虫友见谅。
赞一下(5人) 回复此楼
4楼2011-03-08 11:32:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭