| 查看: 3179 | 回复: 16 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
[交流]
【求助】如何在siesta中设置ghost state原子的计算
|
||||||
|
用siesta计算体相金刚石中一个空位，形成能好像不准，单胞越大，形成能也越大，没有合理的解释。 师兄说，siesta用原子轨道基组，可能存在鬼态原子，需要在空位处设置一个假想的原子，我查了一下手册，没有完全搞清楚如何设置？ 不知道哪位大侠有做过？ |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta问题 | SIESTA | siesta基本概念 |
» 猜你喜欢
281求调剂（0805）
已经有8人回复
-
环境领域全国重点实验室招收博士1-2名
已经有3人回复
-
材料专硕306英一数二
已经有10人回复
-
301求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有7人回复
-
302求调剂
已经有6人回复
-
26博士申请
已经有3人回复
-
268求调剂
已经有3人回复
-
311求调剂
已经有10人回复
-
被我言中：新模板不强调格式了，假专家开始管格式了
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 为什么采用CASTEP计算CaTiO3的光学性质其中的静态介电常数与实验值相差很大
已经有14人回复
- 关于Discover动力学计算的一个问题
已经有26人回复
- 如何计算空位部分所带的电子
已经有5人回复
- 动力学计算过程：时间步长、步数一般设置为多少？
已经有10人回复
- siesta 计算电子占据数的问题
已经有3人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- siesta计算弱相互作用VDW
已经有17人回复
- 【求助】siesta进行ordern计算出错求助
已经有4人回复
- 【求助】siesta关于金属计算的参数设置
已经有5人回复
- 【求助】Ms计算磁性要不要设置每个原子的磁矩？？
已经有10人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
南京医科大学生殖医学与子代健康国重实验室-董飞宏课题组-学术硕士招生
+1/97
-
【代人发帖勿站内信】【成都征男友】89年未婚博士漂亮女孩
+1/71
-
上海交大化院功能大分子团队招2026年硕士、联培、博士后
+1/38
-
上海交通大学电催化方向博后招聘-有锂电、高分子、燃料电池背景者优先
+1/32
-
2026年工科硕士调剂-上海大学全国重点实验室团队-材料数据挖掘方向-研究生3-5人
+1/30
-
海南师范大学孙元元老师招生博士，名额1~2个，2026年9月份入学（光电/光热催化方向）
+1/18
-
医学327分专硕调剂
+1/16
-
【2026博士招生/博后招聘】北京航空航天大学潘彪课题组——AI芯片设计方向
+1/11
-
南京林业大学有机功能分子实验室博士生招生启事
+1/8
-
哈尔滨理工大学材料与化学工程学院 纳米功能材料与光电催化团队 招收硕士研究生
+1/7
-
上海科技大学物质学院招收有机化学方向2026级工程博士
+1/7
-
聊城大学化院青椒招收化学相关专业调剂学生
+1/7
-
福州大学国家级人才汤育欣教授课题组2026年招收”申请考核制“储能电池方向博士研究生
+1/6
-
“机械、材料”招生
+1/6
-
北京高校副校长团队招收机械类，环境类学硕和专硕
+1/6
-
发paper
+1/5
-
2026考研的同学，分数不够？想调剂上岸的同学，来看
+1/5
-
福建师范大学环境微生物技术课题组接受2026级硕士调剂
+1/2
-
上海交通大学化学化工学院张智涛课题组诚聘博士后
+1/2
-
北京大学环境学院张子帅研究员招收科研助理
+1/1
★ ★ ★ ★ ★ ★
小木虫(金币+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
★ ★
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(金币+2): 继续奖励 2011-03-08 15:56:51
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 继续奖励 2011-03-08 15:56:51
|
其实，我之前从来没有做过ghost state的计算，并且对ghost state的知识和用处不了解。仅仅是我之前看手册的时候记得手册中提到了ghost atom的处理。但是具体怎么处理，我是没往心里记的。昨天晚上临回宿舍前，上了一下小木虫，发现了你的这个问题。怕你急于解答，所以我就临时的翻看了一下手册，并且到网上搜索了一下ghos state在siesta中的讨论而已。具体的加入ghost atom之后计算结果应该是什么样子是预期，我不清楚，这个应该你是高手。 需要指明的是，C_g.psf里面的内容不需要改的。表头部分的，程序是不读的，只是给使用者一个参数介绍，让使用者知道这个赝势的类型和参数而已。 |
5楼2011-03-08 11:39:57
