24小时热门版块排行榜    

查看: 3350  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

hw45888792

金虫 (正式写手)

[求助] siesta计算弱相互作用VDW

siesta trunk版本可以计算两个体系之间的弱相互作用VDW,泛函为

XC.functional VDW
XC.authors DRSLL

    不知道有人用过没?我测试发现,好像需要自己重新做赝势?

read_vps: Pseudopotential generation method:
read_vps: ATM3      Troullier-Martins
Total valence charge:    1.00000

xc_check: Exchange-correlation functional:
xc_check: VDW Dion-Rydberg-Schroeder-Langreth-Lundqvist
xc_check: WARNING: Pseudopotential generated with Ceperley-Alder functional
V l=0 = -2*Zval/r beyond r=  1.2343
V l=1 = -2*Zval/r beyond r=  1.2189
V l=2 = -2*Zval/r beyond r=  1.2189
All V_l potentials equal beyond r=  1.2343
This should be close to max(r_c) in ps generation
All pots = -2*Zval/r beyond r=  1.2343

VLOCAL1: 99.0% of the norm of Vloc inside     28.493 Ry
VLOCAL1: 99.9% of the norm of Vloc inside     64.935 Ry
atom: Maximum radius for 4*pi*r*r*local-pseudopot. charge    1.45251
atom: Maximum radius for r*vlocal+2*Zval:    1.21892
GHOST: No ghost state for L =  0
GHOST: WARNING: Ghost state for L =  1
GHOST: No ghost state for L =  2
KBgen: WARNING: Ghost states have been detected
KBgen: WARNING: Some parameter should be changed in the
KBgen: WARNING: pseudopotential generation procedure.
Stopping Program from Node:    0
p0_20173:  p4_error: : 1

[ Last edited by hw45888792 on 2011-4-20 at 15:03 ]
回复此楼
宠辱不惊,看庭前花开花落;去留无意,望天空云卷云舒。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漂泊四方

金虫 (小有名气)

★ ★
franch(金币+2): 鼓励交流 2011-05-31 09:40:18
引用回帖:
Originally posted by hw45888792 at 2011-04-20 15:02:38:
siesta trunk版本可以计算两个体系之间的弱相互作用VDW,泛函为

XC.functional VDW
XC.authors DRSLL

    不知道有人用过没?我测试发现,好像需要自己重新做赝势?

read_vps: Pseudopotential generat ...

你可以问下小燕同学的,她也用siesta做过一些vdW作用力的测试计算,可能有部分结果。

我以前用siesta也做过vdW的计算,vdW赝势也自己做过,不过测试的效果不是太好。问题主要是基组的问题,vdW作用能比较小,而BSSE的影响非常大。
我与我纠缠久
12楼2011-05-30 22:45:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

zhangguangping

木虫 (著名写手)

【答案】应助回帖


zzy870720z(金币+1): 鼓励交流 2011-04-20 21:39:21
hw45888792(金币+10): 2011-05-05 21:45:05
引用回帖:
Originally posted by hw45888792 at 2011-04-20 08:02:38:
siesta trunk版本可以计算两个体系之间的弱相互作用VDW,泛函为

XC.functional VDW
XC.authors DRSLL

    不知道有人用过没?我测试发现,好像需要自己重新做赝势?

read_vps: Pseudopotential generat ...

杨老师的学生胡伟师兄?
我对这个问题不在行。但是好像我听说过一点。
如果是需要改变参数的话,那么你去这个版本的软件包内的atom中看一下,应该这个版本中的atom能够产生VDW的赝势。
弘德明志博学笃行
2楼2011-04-20 16:53:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hw45888792

金虫 (正式写手)

引用回帖:
Originally posted by zhangguangping at 2011-04-20 16:53:30:
杨老师的学生胡伟师兄?
我对这个问题不在行。但是好像我听说过一点。
如果是需要改变参数的话,那么你去这个版本的软件包内的atom中看一下,应该这个版本中的atom能够产生VDW的赝势。

恩,我是胡伟,你是?

这个泛函好像不能计算单个水分子,两个水分子就没有问题了,静态计算VDW比PBE总能还要低~

但是优化起来,收敛性好差,不知道是否还需要设置些什么,我也是第一个玩这个,本来打算用VASP做,但是体系有点大,两三百原子,太慢了~

siesta-L里面有人讨论过这个问题,可惜我一直打不开,不爽~
宠辱不惊,看庭前花开花落;去留无意,望天空云卷云舒。
3楼2011-04-20 21:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖


zzy870720z(金币+1): 鼓励交流 2011-04-20 21:39:31
引用回帖:
Originally posted by hw45888792 at 2011-04-20 14:24:31:
恩,我是胡伟,你是?

这个泛函好像不能计算单个水分子,两个水分子就没有问题了,静态计算VDW比PBE总能还要低~

但是优化起来,收敛性好差,不知道是否还需要设置些什么,我也是第一个玩这个,本来打算用 ...

我是胡伟的师弟。你如果打不开那个网站,可以搜到了之后,然后用网页代理打开。这个方法准灵。http://www.7daili.com/
我就是这样对付SIESTA-L的以前的邮件的。
弘德明志博学笃行
4楼2011-04-20 21:36:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见