| 查看: 1613 | 回复: 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
» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp_DFTP+phonopy 出现负频且有一个声学支斜率为负值
已经有19人回复
- 石墨烯原胞截断能测试，出现问题，请教，谢谢！
已经有5人回复
- CASTEP计算石墨烯如何设置能得到零带隙的结果？
已经有9人回复
- 参数测试的问题
已经有6人回复
- VASP收敛测试
已经有12人回复
- 【分享】常用函数公式及技巧【已搜索无重复】
已经有129人回复
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
- 专业: 理论和计算化学
- 管辖: 仿真模拟
★ ★ ★
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
- 专业: 理论和计算化学
- 管辖: 仿真模拟
送鲜花一朵 |
谢谢回复。 邮件了作者，问题和回复已经帖出来了。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
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
5