应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请教：截断能为什么是409这种奇怪的数字？为什么要交代未成对电子数Ns？
101/1返回列表
查看: 1220  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-09-21 22:34:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有2个 )

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★ ★ ★
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 the
value 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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
cenwanglai: 回帖置顶 2012-09-22 17:14:59
harder version potential 和standard potential,以及soft potential有什么关系呢？在standard, harder和soft版本的时候，依据是什么？
赞一下 回复此楼
10楼2012-09-22 17:02:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢交流 2012-09-22 10:53:56
cenwanglai: 金币+3, ★★★很有帮助, 谢谢回复。409是Cl的hard势里面的ENMAX. 2012-09-22 16:53:46
1.“409”，很有可能是赝势中的ENMAX或者ENMAX的1.3倍，也有可能是作者自己特意设置的这么个值。我觉得跟400或者450基本差别不大。

2.Ns是成单电子数目，关乎是否自旋极化、高低自旋态等。

个人理解，可以告知文献出处么？
赞一下 回复此楼
2楼2012-09-22 00:37:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +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朵）

cenwanglai
3楼2012-09-22 06:35:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

书万里

木虫 (著名写手)

历史控+原著控
★
fzx2008: 金币+1, 谢谢回帖！ 2012-09-22 11:51:13
引用回帖:
3楼: Originally posted by ChemiAndy at 2012-09-22 06:35:22
因为VASP截断半径设置的单位是Ry
而它用的是15Ry, 1Ry=27.21eV, 15Ry=409eV

不对吧，1Ry=13.6ev。1a.u=27.21ev吧。
不过极有可能是用的vasp，截断能是30Ry。
赞一下 回复此楼
喂人民服雾。
4楼2012-09-22 08:07:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主
楼上的Ry单位很有道理
赞一下 回复此楼
5楼2012-09-22 11:36:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
引用回帖:
2楼: Originally posted by fzx2008 at 2012-09-22 00:37:48
1.“409”，很有可能是赝势中的ENMAX或者ENMAX的1.3倍，也有可能是作者自己特意设置的这么个值。我觉得跟400或者450基本差别不大。

2.Ns是成单电子数目，关乎是否自旋极化、高低自旋态等。

个人理解，可以告知 ...

http://muchong.com/bbs/viewthread.php?tid=4986418&page=1#pid1

1楼附件中给出了文献。
赞一下 回复此楼
6楼2012-09-22 16:47:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
引用回帖:
4楼: Originally posted by 书万里 at 2012-09-22 08:07:17
不对吧，1Ry=13.6ev。1a.u=27.21ev吧。
不过极有可能是用的vasp，截断能是30Ry。...

贴图中已经交代了是用VASP算的。
赞一下 回复此楼
7楼2012-09-22 16:48:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
送鲜花一朵
引用回帖:
3楼: Originally posted by ChemiAndy at 2012-09-22 06:35:22
因为VASP截断半径设置的单位是Ry
而它用的是15Ry, 1Ry=27.21eV, 15Ry=409eV

谢谢回复。

邮件了作者，问题和回复已经帖出来了。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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭