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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 高斯高手请进
101/1返回列表
查看: 809  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

x628

金虫 (小有名气)

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

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

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

经验

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-07-23 10:23:49
已阅   回复此楼   关注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 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭