| 查看: 3240 | 回复: 16 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
[交流]
【求助】如何在siesta中设置ghost state原子的计算
|
||||||
|
用siesta计算体相金刚石中一个空位，形成能好像不准，单胞越大，形成能也越大，没有合理的解释。 师兄说，siesta用原子轨道基组，可能存在鬼态原子，需要在空位处设置一个假想的原子，我查了一下手册，没有完全搞清楚如何设置？ 不知道哪位大侠有做过？ |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta问题 | SIESTA | siesta基本概念 |
» 猜你喜欢
297，工科调剂?
已经有5人回复
-
接受任何调剂
已经有6人回复
-
22408 312求调剂
已经有26人回复
-
收到复试调剂但是去不了
已经有8人回复
-
260求调剂
已经有5人回复
-
急需调剂
已经有10人回复
-
申博/考博
已经有4人回复
-
化工学硕294分，求导师收留
已经有37人回复
-
求调剂
已经有11人回复
-
一志愿华中农业071010，320求调剂
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 为什么采用CASTEP计算CaTiO3的光学性质其中的静态介电常数与实验值相差很大
已经有14人回复
- 关于Discover动力学计算的一个问题
已经有26人回复
- 如何计算空位部分所带的电子
已经有5人回复
- 动力学计算过程：时间步长、步数一般设置为多少？
已经有10人回复
- siesta 计算电子占据数的问题
已经有3人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- siesta计算弱相互作用VDW
已经有17人回复
- 【求助】siesta进行ordern计算出错求助
已经有4人回复
- 【求助】siesta关于金属计算的参数设置
已经有5人回复
- 【求助】Ms计算磁性要不要设置每个原子的磁矩？？
已经有10人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
青岛大学/康复大学导师课题组2026年博士研究生
+3/353
-
坐标山东，诚征女友
+1/347
-
复旦大学 课题组 2026年9月批入学-博士招聘
+1/85
-
燕山大学亚稳材料全国重点实验室2026年博士研究生招生信息
+1/83
-
江西科技师范大学086生物与医药
+1/40
-
上海交通大学-化学化工学院 博士后招聘
+1/39
-
硕士招生
+1/38
-
爱情短剧
+1/28
-
26年二次电池储能方向招博士-中山大学四青人才课题组-已招到，此贴作废，谢谢！
+5/20
-
中科院合肥物质研究院招聘水系锌电博士后（长期有效）诚聘海内外优秀博士后
+1/17
-
岭南师范学院化学化工学院生物与医药制药工程专业需要大量调剂生
+1/17
-
南京林业大学-国家级青年人才团队 招2026级申请考核制博士 （有机合成化学方向）
+1/16
-
吉林大学李凯伟课题组招光纤传感器、触觉传感器方向博士生
+1/13
-
中国科学技术大学苏州高等研究院凤建岗课题组诚招博士后
+1/8
-
2026年华南师范大学光电科学与工程学院魏嵬课题组招收硕士调剂生化学专业
+1/7
-
基金祈福
+1/4
-
稀土发光分析｜核壳结构｜寿命拟合
+1/4
-
研究助理招聘 | 中山大学光电材料与技术全国重点实验室--集成信息光电子研究团队
+1/4
-
文献求助
+1/3
-
招07理学调剂生，剩余2个名额
+1/2
7楼2011-03-08 13:21:23
★ ★
zzy870720z(金币+2): 谢谢指教 2011-03-07 22:20:35
hw45888792(金币+4): 好厉害呀，好像回答过我好多问题啊~ 2011-03-08 10:18:18
zzy870720z(金币+2): 谢谢指教 2011-03-07 22:20:35
hw45888792(金币+4): 好厉害呀，好像回答过我好多问题啊~ 2011-03-08 10:18:18
|
参照手册上的ChemicalSpeciesLabel的说明，对于gost原子的原子序数用一个负值就行，然后最好给他起一个别的元素名称，赝势也相应的改变一下。这个时候他的原子质量就会是1.0d30 a.u.如果你计算的是金刚石，那么gost原子的赝势就用C的就行，和你就算中用的其他的C用一个赝势就行，只不过最好起一个别的名字，然后对赝势也重命名一下，将其元素序号用一个负值比如-6表示。 |
2楼2011-03-07 22:14:00
|
你好，这是我测试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
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+5): 谢谢zhang兄指教，说的太好了 2011-03-08 15:56:38
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+5): 谢谢zhang兄指教，说的太好了 2011-03-08 15:56:38
|
其实说句实话，我是真不厉害，很菜的。不是谦虚，我是自我表白。我在论坛里面比较活跃不是显摆我知道的多，而是恰是我知道少，但是身边有没有人用这个软件。所以只能依靠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板块，越来越红火，能成就一些新手。希望老虫多来哺养一下我们这些新虫。 末了，言辞不妥之处，望新老虫友见谅。  |
4楼2011-03-08 11:32:09
