24小时热门版块排行榜    

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

xiaowandouer

木虫 (正式写手)

[交流] 【求助】8核12GB内存的机器,MP2能算多大体系?(以基函数为准)已有10人参与

如题,恳请赐教。目前一个含S、O、C、H的体系,40个原子,单点貌似算不下去,老是这样就over了:

The electronic state of the initial guess is 2-A.
Initial guess = 0.0000 = 0.0000 = 0.5000 = 0.7500 S= 0.5000
Requested convergence on RMS density matrix=1.00D-08 within 500 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
回复此楼

» 猜你喜欢

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

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

mscic

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
[quote]Originally posted by zhangmt at 2010-12-28 17:35:57:
不会吧。600多个基函,是个小体系啊。


600多个basis 还小么?而且是用MP2方法……

要不用组合方法算单点吧
16楼2010-12-29 14:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mscic

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by xiaowandouer at 2010-12-30 10:32:13:



对,是DFT优化好的结构。想着用mp2优化不知道跑到猴年马月,就打算在DFT结果的基础上算个单点了事,不曾想这个也是问题多多。加个P试试,到时候还得请您多指点

为什么用MP2算单点啊?G4不行么?
21楼2010-12-30 14:03:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mscic

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by hairan at 2010-12-30 17:45:48:



G4可能会有些困难,这么大的体系
光是CCSD(T)/6-311G(d)、MP4/6-311G(2df)和MP2/6-311+G(3df,p)这三步能把楼主等崩溃,虽然都是单点

[ Last edited by hairan on 2010-12-30 at  ...

不用程序自带的关键字G4……
24楼2010-12-30 19:11:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mscic

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by xiaowandouer at 2010-12-30 16:01:06:
为什么用MP2算单点啊?G4不行么?

汗,木用过G4,这个能考虑扩散相互作用么?对范德华相互作用的描述怎么样? [/quote]

不知道,你说的这个作用应该在做结构优化的时候就考虑吧?如果B3LYP方法不能很好的描述这个问题,得到的结构应该不可靠吧。
25楼2010-12-30 19:14:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mscic

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by hairan at 2010-12-30 22:49:59:


我一般不用G4,我算的体系都很小,CCSD(T)/aug-cc-pVTZ单点都能做得动,就不用G4费事的组合了。
另外,你是经常用G4但不知道G4是怎么算的?

不需要直接用CCSD(T)/aug-cc-pVTZ算单点,不然就跟你说的一样,不用那么费事了。

%chk=h2.chk
# BMK/6-31+G(2df,p) Opt

Hydrogen molecule G4(MP2)-6X calculation

0 1
H 0.000000 0.000000 0.371049
H 0.000000 0.000000 -0.371049

--Link1--
%chk=h2.chk
# Geom=AllCheck Guess=Read BMK/6-31+G(2df,p) Freq

--Link1--
%chk=h2.chk
# Geom=AllCheck Guess=Read CCSD(T,FrzG4)/GTBas1

--Link1--
%chk=h2.chk
# Geom=AllCheck Guess=Read MP2(FrzG4)/GTMP2LargeXP

--Link1--
%chk=h2.chk
# Geom=AllCheck Guess=Read HF/GFHFB3

--Link1--
%chk=h2.chk
# Geom=AllCheck Guess=Read HF/GFHFB4

这个例子中的那些基组是怎么回事?

[ Last edited by mscic on 2010-12-31 at 11:04 ]
30楼2010-12-31 10:59:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowandouer 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见