24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1799  |  回复: 14

Captain-Jack

木虫 (著名写手)

[求助] 用高斯计算拉曼强度的时候为何总是有虚频?

我用高斯计算某一个激发态下的拉曼强度,使用的分子坐标是经过优化和频率计算(没有虚频)的坐标,但是在拉曼强度的计算中计算的频率却总是有虚频!一下是输入文件:
# b3lyp-td(singlet,root=2)/6-311++g(d,2p) Force(EnOnly) Test scrf=(pcm,solvent=acetonitrile)

3-OPD

0 1
N                 -1.16167000   -1.13379100    0.00000000
C                 -1.12257400    0.20356100    0.00000000
N                  0.00000000    1.04896200    0.00000000
C                  1.25010700    0.30531700    0.00000000
C                  1.23521200   -1.05167500    0.00000000
C                 -0.01273500   -1.74897400    0.00000000
O                 -0.13514300    2.26224700    0.00000000
H                 -2.02556300    0.65801300    0.00000000
H                  2.16353700    0.88185300    0.00000000
H                  2.15648100   -1.61953300    0.00000000
H                 -0.05911400   -2.82848100    0.00000000

36.64

--Link1--
# b3lyp/6-311++g(d,2p)   Freq(HPModes) Test   scrf=(pcm,solvent=acetonitrile)

3-OPD

0 1                                                
N                 -1.16167000   -1.13379100    0.00000000
C                 -1.12257400    0.20356100    0.00000000
N                  0.00000000    1.04896200    0.00000000
C                  1.25010700    0.30531700    0.00000000
C                  1.23521200   -1.05167500    0.00000000
C                 -0.01273500   -1.74897400    0.00000000
O                 -0.13514300    2.26224700    0.00000000
H                 -2.02556300    0.65801300    0.00000000
H                  2.16353700    0.88185300    0.00000000
H                  2.15648100   -1.61953300    0.00000000
H                 -0.05911400   -2.82848100    0.00000000

36.64
可是输出文件中的频率第一个就是虚频,求高手解释!!本人很急!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
cenwanglai: 金币+1, 谢谢回复~ 2012-11-07 09:41:59
你这结构是优化好的吗????????????
2楼2012-11-05 10:08:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Captain-Jack

木虫 (著名写手)

引用回帖:
2楼: Originally posted by lishijunzong at 2012-11-05 10:08:53
你这结构是优化好的吗????????????

结构绝对是优化好的,因为这个结构就是我优化并且计算频率之后没有虚频的结构。大侠,这到底是怎么回事。我其他的物质很好,就这个物质很奇怪。
3楼2012-11-06 08:15:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuzhubing

金虫 (初入文坛)

引用回帖:
3楼: Originally posted by Captain-Jack at 2012-11-06 08:15:01
结构绝对是优化好的,因为这个结构就是我优化并且计算频率之后没有虚频的结构。大侠,这到底是怎么回事。我其他的物质很好,就这个物质很奇怪。...

我这两天也在算这个,指令和你的一模一样,只是基组大小不同。我算的没有虚频啊。用优化好没有虚频的结构。
4楼2012-11-06 10:54:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)

引用回帖:
3楼: Originally posted by Captain-Jack at 2012-11-06 08:15:01
结构绝对是优化好的,因为这个结构就是我优化并且计算频率之后没有虚频的结构。大侠,这到底是怎么回事。我其他的物质很好,就这个物质很奇怪。...

你调整一下就结构,在优化一下。可能是局部最小点

» 本帖已获得的红花(最新10朵)

5楼2012-11-06 14:32:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Captain-Jack

木虫 (著名写手)

送鲜花一朵
引用回帖:
5楼: Originally posted by lishijunzong at 2012-11-06 14:32:14
你调整一下就结构,在优化一下。可能是局部最小点...

如何调整结构?你的意思是把优化好的没有虚频的坐标输出来再一次优化吗?
6楼2012-11-08 09:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Captain-Jack

木虫 (著名写手)

引用回帖:
4楼: Originally posted by xuzhubing at 2012-11-06 10:54:23
我这两天也在算这个,指令和你的一模一样,只是基组大小不同。我算的没有虚频啊。用优化好没有虚频的结构。...

兵哥,你用的是什么基组?我要是把这个基组更改了就和之前所有的数据对不上了。
7楼2012-11-08 09:23:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
Captain-Jack: 金币+3, 谢谢指点,我回去再看看。 2012-11-09 18:36:45
gmy1990: 金币+1 2012-11-09 19:26:59
引用回帖:
6楼: Originally posted by Captain-Jack at 2012-11-08 09:22:09
如何调整结构?你的意思是把优化好的没有虚频的坐标输出来再一次优化吗?...

把优化好的结构加上opt=tight优化,看是否是最稳定状态;查看分子中有没有氢键等作用,尽量是能量最低
8楼2012-11-08 09:35:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gujing136

禁虫 (文坛精英)

【答案】应助回帖

那肯定是结构有问题啊··

» 本帖已获得的红花(最新10朵)

叫老猫的兔子
9楼2012-11-08 09:55:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuzhubing

金虫 (初入文坛)

引用回帖:
7楼: Originally posted by Captain-Jack at 2012-11-08 09:23:32
兵哥,你用的是什么基组?我要是把这个基组更改了就和之前所有的数据对不上了。...

6-31++g(d,p),何师兄......做一下势能面扫描,看看还有没有能量更低的点。

» 本帖已获得的红花(最新10朵)

10楼2012-11-09 11:42:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Captain-Jack 的主题更新
信息提示
请填处理意见