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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » siesta计算弱相互作用VDW
52/1返回列表
查看: 3659  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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 ]
回复此楼

» 猜你喜欢
宠辱不惊，看庭前花开花落；去留无意，望天空云卷云舒。
1楼 2011-04-20 15:02:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

franch

版主 (著名写手)
引用回帖:
Originally posted by 漂泊四方 at 2011-05-30 22:45:22:
你可以问下小燕同学的，她也用siesta做过一些vdW作用力的测试计算，可能有部分结果。

我以前用siesta也做过vdW的计算，vdW赝势也自己做过，不过测试的效果不是太好。问题主要是基组的问题，vdW作用能比较小 ...

我用MS 5.5测试过vdw，，效果也不行。。。。    感觉那些vdw修正都不具有普遍性。
赞一下 回复此楼
遇弱则弱
13楼2011-05-31 09:42:24
已阅   回复此楼   关注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的赝势。
赞一下(1人) 回复此楼
弘德明志博学笃行
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的以前的邮件的。
赞一下(1人) 回复此楼
弘德明志博学笃行
4楼2011-04-20 21:36:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 考博 +4 Eternitf 2026-04-29 4/200 2026-05-02 16:10 by yejenny
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 16:06 by zahslv263y
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[教师之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:36 by zahslv263y
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:21 by zahslv263y
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:06 by zahslv263y
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:21 by zahslv263y
[教师之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:06 by zahslv263y
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 3/150 2026-05-02 10:21 by ovidyohpyz
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-02 10:21 by ovidyohpyz
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-02 10:06 by ovidyohpyz
[教师之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 3/150 2026-05-02 09:36 by dz7y1p1i98
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 5/250 2026-05-02 09:21 by dz7y1p1i98
[考博] 药化及相关博士的申请 5+3 中药学渣 2026-04-29 3/150 2026-05-01 18:14 by babero
[论文投稿] 一篇MDPI论文改变了学习工作和生活 +3 helsinglee 2026-04-29 4/200 2026-05-01 11:45 by bobvan
[有机交流] 一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇 50+4 zhquanbing 2026-04-28 4/200 2026-04-30 22:45 by zyqchem
[文学芳草园] 太白金星有点烦 +3 yexuqing 2026-04-27 3/150 2026-04-30 21:56 by myrtle
[考博] 中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业 +3 liuxh329 2026-04-29 3/150 2026-04-30 15:42 by 顺利毕业zs
[基金申请] 收到国自然专家邀请后几年才会有本子送过来评 +3 barley5 2026-04-30 3/150 2026-04-30 14:58 by aspect3000
[论文投稿] 有没有快的中文核心比较快录用的，纳米材料光催化 5+3 2914439618 2026-04-27 4/200 2026-04-29 08:59 by 北京莱茵润色
信息提示
关闭
请填处理意见
关闭