| 查看: 1491 | 回复: 10 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
yangjun102银虫 (小有名气)
|
[交流]
【求助】siesta中如何将原子设置成离子？已有3人参与
|
||
|
比如说Na，我想把他设置成Na离子，让他带一个单位的正点，他的 block PAO.Basis如何设置，价电子在哪如何设置？ 请教各位高手了！ |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
yangjun102(金币+10):答难解疑，共同讨论！ 2010-04-21 20:53
|
这个我倒是没有试验过，但是我认为如果你的体系中的正负离子带电荷数目相当的话，你也设置一下NetCharge 0，虽然这个是默认的情况。ionic net charge这个地方设置不设置没有什么影响，既然你采用的就是离子赝势。即使你不设置，同样的这个也会自动加上的。你要是严格设置NetCharge 0，体系会认为每个原子都是中性的话，我觉得还有一个可以实验的办法就是你可以让体系带上非常少量的电荷，假如说NetCharge 0.000001，这样程序就会去检查你的这个电荷的来源，所以就会考虑你的离子赝势。这个你可以试验一下。 [ Last edited by zhangguangping on 2010-4-21 at 13:00 ] |
11楼2010-04-21 19:09:44
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
2楼2010-04-21 09:36:54
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★
gavinliu7390(金币+1):谢谢交流！ 2010-04-21 11:29
gavinliu7390(金币+1):谢谢交流！ 2010-04-21 11:29
|
楼主，你可以看一下siesta-3.0-b的手册的第24-25页，这个问题可能能够解决。当时我看手册的时候把这个地方忽略了。 ---------------------------- %block PAO.Basis # Define Basis set O 2 nodes 1.0 # Label, l_shells, type (opt), ionic_charge (opt) n=2 0 2 E 50.0 2.5 # n (opt if not using semicore levels),l,Nzeta,Softconf(opt 3.50 3.50 # rc(izeta=1,Nzeta)(Bohr) 0.95 1.00 # scaleFactor(izeta=1,Nzeta) (opt) 1 1 P 2 # l, Nzeta, PolOrb (opt), NzetaPol(opt) 3.50 # rc(izeta=1,Nzeta)(Bohr) H 1 # Label, l_shells, type (opt), ionic_charge (opt) 0 2 S 0.2 # l, Nzeta, Per-shell split norm parameter 5.00 0.00 # rc(izeta=1,Nzeta)(Bohr) %endblock PAO.Basis ------------------------------------------------------- ionic charge(is): Optional input. Net charge of species is. This is only used for basis set generation purposes. Default value: 0.0 (neutral atom). Note that if the pseudopotential was generated in an ionic configuration, and no charge is specied in PAO.Basis, the ionic charge setting will be that of pseudopotential generation. 看这个解释，好像直接可以在生成赝势的时候解决这个问题。 [ Last edited by zhangguangping on 2010-4-21 at 05:39 ] |
3楼2010-04-21 10:10:26
yangjun102
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 486.5
- 散金: 24
- 帖子: 248
- 在线: 351.2小时
- 虫号: 763027
- 注册: 2009-05-04
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2010-04-21 11:19:43