24小时热门版块排行榜    

查看: 3109  |  回复: 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, 给个红包,谢谢回帖
引用回帖:
8楼: Originally posted by 小范范1989 at 2015-01-15 12:05:19
谢谢神枪指点,我还有几个问题:
1:绝热激发能,就是s1的稳态能量和s0的稳态能量之差,我当时问过,这个帖子:http://ccc.keinsci.com/forum.php?mod=viewthread&tid=365&extra=
2:那我就采用sob的提议 ...

叫小弟“神”是在担当不起。。。
1、如果是这个意思的话那激发态要优化的
3、这个嘛这两个貌似都可以,这个问题以前也有过讨论,到底用哪种输入,由于我三重态做得不多,所以也无法确切回答你,你可以按照sob教你的方法。
    不知道氧气的基态算不算T0?
4、应该可以
9楼2015-01-15 13:01:27
已阅   关注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的回帖

阿Q~~

至尊木虫 (文坛精英)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小范范1989: 金币+10, 谢谢妹子 2015-01-15 09:56:55
路过的看了一下,帮顶~~~~~~~~~~~~~~~~~~
自强不息,厚德载物;独立精神,自由思想。
4楼2015-01-15 09:31:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +5 niko- 2026-03-22 5/250 2026-03-23 22:01 by fuyu_
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 287求调剂 +8 晨昏线与星海 2026-03-19 9/450 2026-03-22 17:01 by i_cooler
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
信息提示
请填处理意见