24小时热门版块排行榜    

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

小范范1989

木虫 (著名写手)

[交流] 高斯09优化激发态并且计算激发态频率的讨论与求教。

先说一下我的背景:点上计算资源有限,50多个人用一个很老的集群,有些老师感觉太满,都自己带电脑,用单机跑分子。
然后,我要计算的分子大约含有原子110个,曾经,我计算一个分子(大约80个原子),激发态的频率计算了半年,出来结果还有一个虚频,我当时的心,就和现在只穿一个背心逛街感觉是一样的。
我的计算的最终目的是:我需要的是基态的频率,还有激发态的频率。

其次:说一下我遇到的问题
    1:基态的优化和频率计算都是在b3lyp/6-31g*的基础上完成的,还算很顺利,没有虚频。耗时也不是很久。
    2:优化激发态,一直出错,采用的是b3lyp/6-31g*.我最近好多帖子都是问这个问题,相信很多虫友已经见过了,基本都能看到我的阴魂。我在想,激发态优化,我都优化2个月,如果计算激发态的频率,是不是要等我博6毕业才能出来?所以我等不起。
    3:看了sob的帖子,说激发态的优化,可以采用pbe0,或者是m062x等,我这就不过多的累述了。给sob一个邮件,sob的回复是:
      至于频率计算,如果你在TDDFT下进行,由于Gaussian没有TDDFT解析二阶导数,这是几乎做不动的,耗时我估计是600~1200次激发能计算的耗时。必须要计算频率的话,只能是在CIS下优化并计算频率,因为CIS下有解析二阶导数,速度和TDDFT计算频率耗时完全不在一个数量级。频率对计算级别的敏感度远小于激发能,因此CIS下的频率的准确度应当没问题。
    这是的想问的主要问题:sob建议我激发态的优化和频率的计算使用CIS是吧?我第一感觉,是不是CIS不更耗时?再就是,我激发态用cis,基态用b3lyp,然后拿这两个频率去进行下一步的计算,(清华的momoap计算杜欣斯基转动矩阵)。能靠谱吗?
   还有一事就是,CIS计算,就是把我的b3lyp换成CIS就行了是吧?其他的输入文件的地方不用改?
谢谢指点,谢谢。
回复此楼

» 收录本帖的淘帖专辑推荐

经验

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

阿Q~~

至尊木虫 (文坛精英)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小范范1989: 金币+10, 谢谢妹子 2015-01-15 09:56:55
路过的看了一下,帮顶~~~~~~~~~~~~~~~~~~
自强不息,厚德载物;独立精神,自由思想。
4楼2015-01-15 09:31:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

小范范1989

木虫 (著名写手)

同时,我还有一个问题:
我优化出来了T1的结构的方式是:
#p b3lyp/6-31* opt scf(maxcycle=500)

test

0 3

计算频率我的方式是:
#p b3lyp/6-31* freq

test

0 3

都正常结束了,而且也没有虚频,不知道这样的结果可靠不?
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
2楼2015-01-15 00:00:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
小范范1989: 金币+10 2015-01-15 09:56:45
基态和激发态可以用不同的理论水平。CIS会比较不耗时,但是精确度不行,所以sob建议你激发态能量取DFT的,但是频率可以用CIS的。
如果你T1的频率是在T1的驻点结构上算的话就没有问题了。
3楼2015-01-15 08:53:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

引用回帖:
3楼: Originally posted by 枪下游魂 at 2015-01-15 08:53:31
基态和激发态可以用不同的理论水平。CIS会比较不耗时,但是精确度不行,所以sob建议你激发态能量取DFT的,但是频率可以用CIS的。
如果你T1的频率是在T1的驻点结构上算的话就没有问题了。

谢谢枪神的回复,我打算按照sob的指点试试,激发态我用cis优化,并且算频率。但是对于另外的参数,比如说绝热激发能和垂直激发能,我是不是应该按照我的b3lyp/6-31g*的形式去找?但是,激发态6-31g*出不来,3-21g出来,能否用基态是b3lyp/6-31g*,激发态优化是3-21g,来找绝热激发能?
可是我感觉,我整个的计算的流程,基态和激发态优化的方式不一样,基泰和激发态频率的方式不一样,在进行下一步,借助清华momap去算辐射跃迁,无辐射跃迁速率,这样一系列的省略下来,不知道最后的结果和实验差不多少。一直担心这个。
谢谢神枪的指点,后面的如果运到其他我不能单独解决的问题,还希望您指点一二,thanks
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
5楼2015-01-15 09:56:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 复试调剂 +21 积极向上; 2026-04-10 23/1150 2026-04-15 12:50 by 西北望—风沙
[考研] 367求调剂 +11 hffQAQ 2026-04-09 11/550 2026-04-14 17:48 by lhj2009
[考研] 药学305求调剂 +10 玛卡巴卡boom 2026-04-10 10/500 2026-04-14 15:55 by zs92450
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 332求调剂 +15 蕉蕉123 2026-04-10 15/750 2026-04-13 23:12 by pies112
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 一志愿京区985,085401电子信息,本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 初试分332,一志愿报考西北工业大学, +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
信息提示
请填处理意见