24小时热门版块排行榜    

查看: 2969  |  回复: 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, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by 小范范1989 at 2015-01-15 09:56:40
谢谢枪神的回复,我打算按照sob的指点试试,激发态我用cis优化,并且算频率。但是对于另外的参数,比如说绝热激发能和垂直激发能,我是不是应该按照我的b3lyp/6-31g*的形式去找?但是,激发态6-31g*出不来,3-21g出 ...

那6-31g可以出得来吗?垂直激发能应该没问题,毕竟不涉及结构优化。绝热激发能是怎么个情况?我在理论上确实小白。。。
就像sob神说的,如果你单单只需要频率数据的话其精度的影响不会太大。
另外不好意思我看错了一点,你算的应该是三重态的基态,也就是T0,如果是T1的话应该用td(triplet)来算。T1的计算应该也有很多老帖讨论过的
最后cis的写法是替换td,不是b3lyp,如果我没记错的话
7楼2015-01-15 11:02:05
已阅   关注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的回帖
信息提示
请填处理意见