24小时热门版块排行榜    

查看: 2482  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +9 我4大白菜 2026-08-20 9/450 2026-08-20 18:23 by 2549091707
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 言之凿凿,明天官网会通知放榜? +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见