24小时热门版块排行榜    

查看: 3113  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 芦lty 2026-03-25 4/200 2026-03-25 23:25 by 芦lty
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 332求调剂 +6 032500 2026-03-25 6/300 2026-03-25 22:45 by 418490947
[考研] 考研调剂 +6 来好运来来来 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 材料学硕333求调剂 +4 北道巷 2026-03-24 4/200 2026-03-25 14:16 by mapenggao
[考研] 081700 调剂 267分 +10 迷人的哈哈 2026-03-23 10/500 2026-03-25 12:12 by userper
[考研] 289求调剂 +9 怀瑾握瑜l 2026-03-20 9/450 2026-03-25 11:02 by userper
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 11/550 2026-03-25 10:23 by userper
[考研] 0703化学求调剂 +6 奶油草莓. 2026-03-22 7/350 2026-03-25 10:00 by shangxh
[考研] 085601求调剂总分293英一数二 +3 钢铁大炮 2026-03-24 3/150 2026-03-24 22:03 by bingxueer79
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
信息提示
请填处理意见