应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 高斯09优化激发态并且计算激发态频率的讨论与求教。
101/1返回列表
查看: 2966  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小范范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.
1楼 2015-01-14 17:18:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

小范范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的回帖

小范范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的回帖

小范范1989

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

对，还有一事相求：
我的优化t1的表头，采用的是0  3的方式，
优化的时候也是  0  3，类似上面写的，
这个应该没问题哈，第一次计算t1，还请您鉴定一下。
ps：结果很顺利的出来了，没有虚频，耗时不是很长。
赞一下 回复此楼
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
6楼2015-01-15 09:59:02
已阅   关注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的回帖

小范范1989

木虫 (著名写手)
引用回帖:
7楼: Originally posted by 枪下游魂 at 2015-01-15 11:02:05
那6-31g可以出得来吗？垂直激发能应该没问题，毕竟不涉及结构优化。绝热激发能是怎么个情况？我在理论上确实小白。。。
就像sob神说的，如果你单单只需要频率数据的话其精度的影响不会太大。
另外不好意思我看错 ...

谢谢神枪指点，我还有几个问题：
1：绝热激发能，就是s1的稳态能量和s0的稳态能量之差，我当时问过，这个帖子：http://ccc.keinsci.com/forum.php ... =365&extra=
2：那我就采用sob的提议，CIS计算频率。
3：难道以前我都理解错了？我以为我的这个方式就是优化T1,
sob的回复：又是三重态又是root=1，直接用0 3计算就行了，不需要TDDFT，又快又不会出现这些问题
我当时在这个帖子问过：
http://ccc.keinsci.com/forum.php ... =369&extra=
再就是，三重态基态，是T0？T0没见过啊，我以前都认为是t1，还请枪神指点。
4：我昨晚上看了看，是不是这样写  #p RCIS(root=1,nstates=40)/6-31g* opt freq
赞一下 回复此楼
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
8楼2015-01-15 12:05:19
已阅   关注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的回帖

小范范1989

木虫 (著名写手)
引用回帖:
9楼: Originally posted by 枪下游魂 at 2015-01-15 13:01:27
叫小弟“神”是在担当不起。。。
1、如果是这个意思的话那激发态要优化的
3、这个嘛这两个貌似都可以，这个问题以前也有过讨论，到底用哪种输入，由于我三重态做得不多，所以也无法确切回答你，你可以按照sob教你 ...

你就是很厉害，很早之前就关注你了，很多帖子你回复的都很好，而且很到位。你在我心里就是大神一样的人物。
哈哈哈，谢谢你的帮助，希望以后常交流。谢谢你
赞一下 回复此楼
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
10楼2015-01-15 15:02:50
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小范范1989 的主题更新
101/1返回列表
信息提示
关闭
请填处理意见
关闭