24小时热门版块排行榜    

查看: 241  |  回复: 2
当前主题已经存档。

journalemu

[交流] 【求助】Gaussian如何打印出双电子积分?

如题,请求高手给予指点,如蒙您指点,不胜感激。多谢。
比如,就 hf 计算一下 H2分子的能量,选择某个基组,基组个数为10吧,我想要全部的双电子积分数据,就是计算的时候设置不考虑对成性。总共:1540个双电子积分吧。多谢高手帮助。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wolfzhong

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
journalemu(金币+10,VIP+0):非常感谢,我尝试了一下,通过不了,嘿嘿,一直报错,不知道是为啥,不知道高手是否能给尝试一下?我不知道哪里错了。谢谢了。 1-14 09:09
journalemu(金币+10,VIP+0):多谢了。 1-14 10:19
From 百度


gaussian 默认是不输出各种电子积分的。如果想要得到这些积分,需要如下的工作。

在你的routes中加入#P ExtraLinks=L316   iop(3/33=6,5/33=3)   symm=noint noraff

ExtraLinks=L316 表示要执行links 316,意思是pints 2-electron integrals

iop(3/33=6,5/33=3) 表示控制积分打印,3/33=6,表示打印双电子积分和单电子积分;其中双电子积分的格式为debug格式。5/33=3表示控制输出,在每一次叠代后输出结果。

symm=noint 表示在计算积分的时候不使用对称性。如果使用对称性,可以减少积分的数量。

noraff 也可以写成 integeral(noraff),表示双电子输出使用regular integral format。

例子:#p hf/sto-3g extralinks(L316) IOP(3/33=6) scf=conventional Integral(noraff)
2楼2010-01-14 00:33:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wolfzhong

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
journalemu(金币+20,VIP+0):我再尝试一下看看,多谢了。 1-14 10:17
我做了个H2的 HF/sto-3g的结果可以的哦
3楼2010-01-14 09:46:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 journalemu 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见