24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1845  |  回复: 33

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的回帖

xiaowandouer

木虫 (正式写手)

为什么用MP2算单点啊?G4不行么? [/quote]


汗,木用过G4,这个能考虑扩散相互作用么?对范德华相互作用的描述怎么样?
22楼2010-12-30 16:01:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by mscic at 2010-12-30 14:03:50:


为什么用MP2算单点啊?G4不行么?


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 17:49 ]
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
23楼2010-12-30 17:45:48
已阅   回复此楼   关注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的回帖

hzw83

铁虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by xiaowandouer at 2010-12-28 09:44:00:
如题,恳请赐教。目前一个含S、O、C、H的体系,40个原子,单点貌似算不下去,老是这样就over了:

The electronic state of the initial guess is 2-A.
Initial guess = 0.0000 = 0.000 ...

跑你的体系应该没问题,我跑过200个原子的,在6-31+g*的基组下,在DFT下基组加到6-311++g**没问题
26楼2010-12-30 19:51:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

引用回帖:
Originally posted by mscic at 2010-12-30 19:14:31:



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

不知道,你说的这个作用应该在做结构优化的时候就考虑吧?如果B3LYP方法不能很好的描述这个问题,得到的结构应该不可 ... [/quote]


你说的没错,但是G03中的泛函貌似都没有考虑范德华相互作用,我的体系用MP2优化的话又不现实,只好做个单点。起码在现有构型下,用MP2得到的能量会更准确些

G09中有些泛函能够考虑范德华相互作用,但是其可靠性又无法验证。。。。。
27楼2010-12-30 21:23:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

引用回帖:
Originally posted by hzw83 at 2010-12-30 19:51:29:

跑你的体系应该没问题,我跑过200个原子的,在6-31+g*的基组下,在DFT下基组加到6-311++g**没问题

谢谢哈,题正在算,看看这次能不能跑下去
28楼2010-12-30 21:27:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫


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


不用程序自带的关键字G4……

我一般不用G4,我算的体系都很小,CCSD(T)/aug-cc-pVTZ单点都能做得动,就不用G4费事的组合了。
另外,你是经常用G4但不知道G4是怎么算的?
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
29楼2010-12-30 22:49:59
已阅   回复此楼   关注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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见