24小时热门版块排行榜    

查看: 1031  |  回复: 9

x628

金虫 (小有名气)

[交流] 高斯高手请进 已有1人参与

wb97xd默认omega0.2,输入文件加什么关键词能把omega定义为0, 0.1, 0,15?
回复此楼

» 收录本帖的淘帖专辑推荐

经验

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
zhou2009: 金币+5 2014-07-24 10:32:24
参看IOP指令集中IOp(3/108)和IOp(3/109)两项,这两项是调整omega值的,前者针对Ex部分,后者针对Ec部分。查阅了蔡政达、Head-Gordon的wb97xd原文献,range separated的部分是在Ex的部分,并且omega长程短程的值相等,如取0.15,我想应该是加上IOp(3/108=0150001500)。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2014-07-24 00:23:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
zhou2009: 金币+2 2014-07-24 10:32:56
额,忘了还有HF项,IOp(3/107=0150001500)。
泛函的Ex部分由LR-HF、SR-HF和SR-B97构成。
不一定挂在论坛,计算问题问题欢迎留言。
3楼2014-07-24 00:35:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x628

金虫 (小有名气)

太好了,谢谢您呀!!
4楼2014-07-24 10:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x628

金虫 (小有名气)

呃,加下这个关键词后,高斯秒错,l502。
5楼2014-07-24 11:46:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by x628 at 2014-07-24 11:46:07
呃,加下这个关键词后,高斯秒错,l502。

maybe我理解错了,IOp(3/107)不加没问题,加上就出错。
我试了一下,只考虑IOp(3/108=0150001500)就没问题。
不一定挂在论坛,计算问题问题欢迎留言。
6楼2014-07-24 11:53:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x628

金虫 (小有名气)

Symmetry not used in FoFCou.
FMM levels:  10  Number of levels for PrismC:   9
nThT2=3 NYI in PrismS.
Error termination via Lnk1e in /scratch/g09/l502.exe at Fri Jul 25 11:01:50 2014.
Job cpu time:  0 days  0 hours  4 minutes 40.3 seconds.
7楼2014-07-24 11:54:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x628

金虫 (小有名气)

是的,不用107就能算了。另外omega是不是只是单指long?那0.15是不是应该IOp(3/108=0000001500)?
8楼2014-07-24 12:16:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
8楼: Originally posted by x628 at 2014-07-24 12:16:12
是的,不用107就能算了。另外omega是不是只是单指long?那0.15是不是应该IOp(3/108=0000001500)?

B97的部分只有短程部分。
LR-HF、SR-HF和SR-B97。
其形式为,a*LR-HF(omega=0.15)+b*SR-HF(omega=0.15)+c*LR-B97(omega=0.15)+d*SR-B97(omega=0.15)=ExwB97

这里c已经=0了,对于这项的omega=多少无所谓,我认为可以不用调,因为无论长程omega=多少,B97长程的term系数为0,这样就无关紧要了。

原文中HF的term 和B97的term对应的omega并没有采用不同的标志,长程和短程也没用不同的omega,所以设定为同一个值就可以了。

omega指的是分开长程和短程的一个参数,不是只是指的long的部分。
1/r12=erf(r*omega)/r12+erfc(r*omega)/r12.

理解的对不对就看各位前辈斧正了。
不一定挂在论坛,计算问题问题欢迎留言。
9楼2014-07-24 12:58:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x628

金虫 (小有名气)

谢谢您
10楼2014-07-24 21:07:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 x628 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +20 木水思豆 2026-08-25 23/1150 2026-08-26 06:35 by echo8914667
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 PbwVEMK5haaE 2026-08-25 5/250 2026-08-26 05:09 by cNXvBfCpiZOM
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Qj98d49dHFb9 2026-08-25 5/250 2026-08-26 04:57 by cNXvBfCpiZOM
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 今天务委会开完了,明天出结果吗 +18 angus9576 2026-08-25 22/1100 2026-08-26 00:41 by merchancy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见