24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3183  |  回复: 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的回帖

枪下游魂

木虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+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的回帖
查看全部 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的回帖

阿Q~~

至尊木虫 (文坛精英)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小范范1989: 金币+10, 谢谢妹子 2015-01-15 09:56:55
路过的看了一下,帮顶~~~~~~~~~~~~~~~~~~
自强不息,厚德载物;独立精神,自由思想。
4楼2015-01-15 09:31:18
已阅   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +7 关忆北. 2026-04-14 7/350 2026-04-16 01:30 by L0alice
[考研] 初试324 中药学 一志愿天中医 求调剂 +3 李的Lucky 2026-04-10 3/150 2026-04-15 23:54 by 陈皮皮
[考研] 273求调剂 +6 白居不易. 2026-04-09 8/400 2026-04-15 22:02 by wooluyong
[考研] 材料相关专业344求调剂双非工科学校或课题组 +21 hualkop 2026-04-12 23/1150 2026-04-15 22:02 by hualkop
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 化学070300 求调剂 +23 哈哈哈^_^ 2026-04-12 23/1150 2026-04-14 16:30 by zhouxiaoyu
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 290求调剂 +18 柯淮然 2026-04-12 20/1000 2026-04-13 12:56 by cyh—315
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 085400 328分 求调剂 +10 喂你一个大橙子 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 346,工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 0858求调剂 5+5 Gky09300550, 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 一志愿中科大070300化学,314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
信息提示
请填处理意见