| 查看: 3126 | 回复: 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就行了是吧？其他的输入文件的地方不用改？ 谢谢指点，谢谢。 |
回复此楼» 收录本帖的淘帖专辑推荐
 经验 |
» 猜你喜欢
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有188人回复
-
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
已经有11人回复
-
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
-
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有17人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
070300总分298求调剂
已经有2人回复
-
一志愿厦门大学化学学硕307求调剂
已经有10人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 这是高斯09计算圆锥交叉点的算例，求解关键字
已经有11人回复
- 关于激发态T1的构型优化方法哪种比较好？
已经有6人回复
- 高斯分子优化及频率计算
已经有4人回复
- 高斯算频率的问题
已经有5人回复
- 高斯计算-激发态优化
已经有13人回复
- Gaussian TDDFT 激发态 频率计算 输入文件
已经有8人回复
- 关于荧光的高斯计算
已经有13人回复
- 利用TDDFT在优化基础上做激发态性质计算，得到的激发态能级是怎么来的。
已经有5人回复
- 求助有关高斯计算激发态能量准确度问题
已经有6人回复
- 求教：高斯如何计算单重态（激发态）过渡态
已经有7人回复
- 急求用高斯09怎么计算电子激发态
已经有7人回复
- td做激发态优化出错
已经有19人回复
- 高斯计算激发态并行问题
已经有6人回复
- 用高斯09计算的tddft方法优化激发态出现问题如何解决？
已经有15人回复
- 【求助】关于激发态优化及发射光谱的计算
已经有36人回复
- 【求助】TDDFT方法优化激发态问题
已经有5人回复
- 【求助】关于第一激发态跃迁能？？
已经有10人回复
- 【求助】用TDDFT优化出来激发态之后的布居分析
已经有20人回复
- 【求助】高斯频率计算死掉的原因
已经有7人回复
- 【求助】高斯优化过渡态出错
已经有14人回复
- 【求助】关于TDDFT计算激发态的问题
已经有7人回复
- 【求助】TD计算激发态时Davidson循环不过
已经有9人回复
- 【求助】单重激发态的高斯算法
已经有29人回复
- 【求助】Gaussian09 优化激发态出错
已经有8人回复
小范范1989
木虫 (著名写手)
- 应助: 123 (高中生)
- 金币: 3821.3
- 散金: 968
- 红花: 243
- 帖子: 1979
- 在线: 452.8小时
- 虫号: 2013403
- 注册: 2012-09-19
- 性别: GG
- 专业: 原子和分子物理
2楼2015-01-15 00:00:24
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
3楼2015-01-15 08:53:31
阿Q~~
至尊木虫 (文坛精英)
- 应助: 77 (初中生)
- 金币: 33343.6
- 红花: 107
- 帖子: 36186
- 在线: 1028.1小时
- 虫号: 3237212
- 注册: 2014-05-27
- 性别: MM
- 专业: 宇宙学
4楼2015-01-15 09:31:18
小范范1989
木虫 (著名写手)
- 应助: 123 (高中生)
- 金币: 3821.3
- 散金: 968
- 红花: 243
- 帖子: 1979
- 在线: 452.8小时
- 虫号: 2013403
- 注册: 2012-09-19
- 性别: GG
- 专业: 原子和分子物理
|
谢谢枪神的回复，我打算按照sob的指点试试，激发态我用cis优化，并且算频率。但是对于另外的参数，比如说绝热激发能和垂直激发能，我是不是应该按照我的b3lyp/6-31g*的形式去找？但是，激发态6-31g*出不来，3-21g出来，能否用基态是b3lyp/6-31g*，激发态优化是3-21g，来找绝热激发能？ 可是我感觉，我整个的计算的流程，基态和激发态优化的方式不一样，基泰和激发态频率的方式不一样，在进行下一步，借助清华momap去算辐射跃迁，无辐射跃迁速率，这样一系列的省略下来，不知道最后的结果和实验差不多少。一直担心这个。 谢谢神枪的指点，后面的如果运到其他我不能单独解决的问题，还希望您指点一二，thanks |
5楼2015-01-15 09:56:40
小范范1989
木虫 (著名写手)
- 应助: 123 (高中生)
- 金币: 3821.3
- 散金: 968
- 红花: 243
- 帖子: 1979
- 在线: 452.8小时
- 虫号: 2013403
- 注册: 2012-09-19
- 性别: GG
- 专业: 原子和分子物理
6楼2015-01-15 09:59:02
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
7楼2015-01-15 11:02:05
小范范1989
木虫 (著名写手)
- 应助: 123 (高中生)
- 金币: 3821.3
- 散金: 968
- 红花: 243
- 帖子: 1979
- 在线: 452.8小时
- 虫号: 2013403
- 注册: 2012-09-19
- 性别: GG
- 专业: 原子和分子物理
|
谢谢神枪指点，我还有几个问题： 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 |
8楼2015-01-15 12:05:19
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
9楼2015-01-15 13:01:27
小范范1989
木虫 (著名写手)
- 应助: 123 (高中生)
- 金币: 3821.3
- 散金: 968
- 红花: 243
- 帖子: 1979
- 在线: 452.8小时
- 虫号: 2013403
- 注册: 2012-09-19
- 性别: GG
- 专业: 原子和分子物理
10楼2015-01-15 15:02:50
10