| 查看: 665 | 回复: 3 | |||
| 当前主题已经存档。 | |||
xirainbow木虫 (正式写手)
|
[交流]
【原创】What is nelec
|
||
|
以下内容来自我的博客：http://blog.sina.com.cn/s/blog_5f15ead20100cw4o.html Parameter nelec is explained in PWscf as " number of electron in the unit cell". It seems very easy to be understood. But for a Carbon atom, nelec is 6 or 4!?. I use US PP and tried 6, then scf will give very large error. It seems nelec is concerned with the PP. Furthermore, I found a example in "User’s Guide for Quantum-ESPRESSO (version 3.2)" "Example of input file (Benzene Molecule): &control title = ’ Benzene Molecule ’, calculation = ’cp’, pseudo_dir=’/scratch/acv0/benzene/’, outdir=’/scratch/ac ..... &system nat = 12, ntyp = 2, nbnd = 15, nelec = 30, ecutwfc = 40.0, ATOMIC_SPECIES C 12.0d0 c_blyp_gia.pp H 1.00d0 h.ps ATOMIC_POSITIONS (bohr) C 2.6 0.0 0.0 C 1.3 -1.3 0.0 C -1.3 -1.3 0.0 C -2.6 0.0 0.0 C -1.3 1.3 0.0 C 1.3 1.3 0.0 H 4.4 0.0 0.0 H 2.2 -2.2 0.0 H -2.2 -2.2 0.0 H -4.4 0.0 0.0 H -2.2 2.2 0.0 H 2.2 2.2 0.0" Therefore, I believe that the nelec for carbon atom should be the valence electron number, not all electron number. In the US PP file, I find the following content: " SLA PZ NOGX NOGC PZ Exchange-Correlation functional 4.00000000000 Z valence -10.68218164300 Total energy" |
回复此楼» 猜你喜欢
生物学308分求调剂（一志愿华东师大）
已经有12人回复
-
调剂求助（生物与医药）
已经有3人回复
-
一志愿吉大化学327求调剂
已经有7人回复
-
一志愿生物与医药，296分，求调剂
已经有13人回复
-
327求调剂
已经有5人回复
-
初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文
已经有8人回复
-
327考研调剂推荐
已经有6人回复
-
328求调剂
已经有3人回复
-
304求调剂
已经有9人回复
-
化工求调剂！
已经有16人回复
doudouding
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2306.2
- 红花: 1
- 帖子: 182
- 在线: 47.2小时
- 虫号: 64507
- 注册: 2005-04-16
- 专业: 凝聚态物性I:结构、力学和
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+2,VIP+0):欢迎讨论 10-24 22:23
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+2,VIP+0):欢迎讨论 10-24 22:23
| 和nelec相关联的就是这个nbnd ，说明书里说“number of electronic states (bands) to be calculated.” default里说是“for an insulator, nbnd = number of valence bands (nbnd=nelec/2, see below for nelec); for a metal, 20% more (minimum 4 more)”。没有地方能找到这个东西是做什么用的，什么时候需要设定这个值，学这个软件真是有点窝火，说明书写的太简单，好像是写给作者自己看的，只是给自己有个提示就可以了。论坛里也很少有人能像楼主一样说出来自己的体会，估计是大家都没什么谱，怕自己蒙的不对。 |
2楼2009-10-24 15:10:28
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
★ ★ ★
aylayl08(金币+3,VIP+0):版里貌似学PWscf的不是很多，呵呵，谢谢分享你的经验 10-25 09:36
aylayl08(金币+3,VIP+0):版里貌似学PWscf的不是很多，呵呵，谢谢分享你的经验 10-25 09:36
|
当初自己在安装PWscf的时候也是一个人在那里瞎折腾。因为以前连Linux系统都没接触过，而且身边没有任何可以求助的朋友，所以遇到了N多问题，及其痛苦的折腾了1一个月，终于把Linux系统装好，然后又把PWscf编译成功。期间遇到的问题都是通过google上查找到的方法解决的。我非常感谢那些网上的素不相识的朋友些出的经验总结，没有他们的帖子，我什么都做不了。因此，我觉得我应该吧自己关于计算的一些心得写出来，希望可以给他人以参考：） 这个软件的说明确实写的有点简单，需要有足够的理论基础，才能理解其真正的含义：） 至于nbnd的在计算中究竟代表什么，你可以参看以下http://blog.sina.com.cn/s/blog_5f15ead20100e1l5.html 这只是我个人的体会，还请多多指正；） |
3楼2009-10-24 22:38:32
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+4):感谢讨论 2010-03-15 08:41
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+4):感谢讨论 2010-03-15 08:41
| 个人觉得这个已经说的很明白了，金属和绝缘体的区别，这个需要你的固体物理的知识。你不理解我想是你没有把你学过的知识和软件结合起来而已。nbnd当然就是能带数目，每个能带可以占据2个电子，也就是对于理想晶体而言，尤其是零温体系，价电子数目刚好占据满价带，且导带与价带之间有很宽的带宽，电子一般不会跃迁上去。而对于金属，由于导带和价带没有间隔，电子可以在导带上迁移，因此产生电流，仅仅限制电子在Fermi能级以下是不对的。当然，这里又涉及到费米能级的定义，这又是另外的问题了。 |
4楼2010-03-14 21:57:47
