| 查看: 1587 | 回复: 9 | ||
cenwanglai荣誉版主 (知名作家)
老和山猥琐派九段
|
[求助]
请教：截断能为什么是409这种奇怪的数字？为什么要交代未成对电子数Ns？
|
|
|
看看这个交代计算的细节。ENCUT = 409 eV是如何得到的？ Ns的设置是如何操作的，有什么意思？  hehe.jpg[ Last edited by cenwanglai on 2012-9-22 at 16:46 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : jp207342z[1]_Cl_JPCC.pdf
2012-09-22 16:46:28, 2.93 M
» 猜你喜欢
能否退出参与的面上项目解除限项
已经有15人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
-
科研孤儿太难了
已经有19人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp_DFTP+phonopy 出现负频且有一个声学支斜率为负值
已经有19人回复
- 石墨烯原胞截断能测试，出现问题，请教，谢谢！
已经有5人回复
- CASTEP计算石墨烯如何设置能得到零带隙的结果？
已经有9人回复
- 参数测试的问题
已经有6人回复
- VASP收敛测试
已经有12人回复
- 【分享】常用函数公式及技巧【已搜索无重复】
已经有129人回复
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
★ ★ ★
cenwanglai: 回帖置顶 2012-09-22 16:49:38
fzx2008: 金币+3, 谢谢分享 2012-09-23 16:30:23
cenwanglai: 回帖置顶 2012-09-22 16:49:38
fzx2008: 金币+3, 谢谢分享 2012-09-23 16:30:23
|
给作者发邮件问了，问题和回复如下： 1. Why the ENCUT = 409 eV is used? How dit you get the value 409? ENCUT=409eV corresponds to the cutoff used for Cl in the "hard" pseudopotential (Cl_h). If you looked for the keyword "ENMAX" into the POTCAR file corresponding to the hard Cl pseudopotential, you'll get: ENMAX = 409.136 To avoid comparing calculations obtained using different basis sets (e.i., different cutoff for the plane waves expension), we explicitly set ENCUT=409 in all the calculations including those not involving Cl (for instance the oxygen molecule in the gas phase). > 2. The number of unpaired electrons (Ns) was considered. How can you assign the Ns and fix it during geo optimaztion? Could you please show me the corresponding TAG in INCAR or other ways to make it? The number of unpaired electrons is fixed by the keyword "NUPDOWN". For instance "NUPDOWN=1" in the INCAR file means Ns=1 in the paper. The value of "NUPDOWN" appears in the "OSZICAR" and "out" files after each electronic minimization. Look at the number given after "mag" (for instance "mag= 0.0000" means "NUPDOWN = 0"  . If you don't set thevalue of "NUPDOWN" in the INCAR file, VASP will compute the magnetization automatically. The problem with this approach is that the total magnetic moment may vary from one step to another during the geometry optimization. Personally, I prefer to fix NUPDOWN during the calculation and try different value of NUPDOWN. This ensure that I get the lowest "spin state" for a given geometry. Furthermore, this approach ensure the fulfilment of the "spin conservation" (see S. Chretien and H. Metiu J. Chem. Phys. 129(7), 074705 (2008)). |
8楼2012-09-22 16:49:15
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
10楼2012-09-22 17:02:09
fzx2008
荣誉版主 (著名写手)
- 1ST强帖: 9
- 应助: 272 (大学生)
- 贵宾: 1.045
- 金币: 21814.1
- 散金: 710
- 红花: 33
- 沙发: 1
- 帖子: 1808
- 在线: 1601小时
- 虫号: 621341
- 注册: 2008-10-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
2楼2012-09-22 00:37:48
ChemiAndy
木虫 (正式写手)
- 1ST强帖: 4
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.8小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢交流 2012-09-22 10:54:02
cenwanglai: 金币+3, ★★★很有帮助, 感谢回复~ 2012-09-22 17:02:42
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢交流 2012-09-22 10:54:02
cenwanglai: 金币+3, ★★★很有帮助, 感谢回复~ 2012-09-22 17:02:42
|
因为VASP截断半径设置的单位是Ry 而它用的是15Ry, 1Ry=27.21eV, 15Ry=409eV |
» 本帖已获得的红花（最新10朵）
3楼2012-09-22 06:35:22
书万里
木虫 (著名写手)
历史控+原著控
- 应助: 70 (初中生)
- 金币: 3295.6
- 散金: 742
- 红花: 8
- 沙发: 4
- 帖子: 1519
- 在线: 466.8小时
- 虫号: 1244505
- 注册: 2011-03-25
- 性别: GG
- 专业: 中国通史
4楼2012-09-22 08:07:17
fzx2008
荣誉版主 (著名写手)
- 1ST强帖: 9
- 应助: 272 (大学生)
- 贵宾: 1.045
- 金币: 21814.1
- 散金: 710
- 红花: 33
- 沙发: 1
- 帖子: 1808
- 在线: 1601小时
- 虫号: 621341
- 注册: 2008-10-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
5楼2012-09-22 11:36:55
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
6楼2012-09-22 16:47:22
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
7楼2012-09-22 16:48:17
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
送鲜花一朵 |
谢谢回复。 邮件了作者，问题和回复已经帖出来了。409是Cl的hard势的ENMAX. VASP手册上说： For some elements several PAW versions exist. The standard version has generally no extension. An extension _h implies that the potential is harder than the standard potential and hence requires a larger energy cutoff. The enxtension _s means that the potential is softer than the standard version. |
9楼2012-09-22 16:59:00
20