本帖最后由 huangyoyo 于 2019-3-25 09:19 编辑 sobereva 发表于 2019-3-22 19:22 那么例子中优化所有的激发态的目的是什么? |
huangyoyo 发表于 2019-3-22 16:25 ORCA都不支持TDDFT二阶解析导数,按手册里的方式去算你的体系,猴年马月去 |
Novice 发表于 2019-3-21 11:03 如果要模拟振动分辨的精细光谱,不是要输入每个激发态的Hess文件么? |
huangyoyo 发表于 2019-3-21 09:37 这个任务我也做了,你这么写估计能算一个月,这是ORCA模拟的振动分辨的精细光谱,Sob老师说的很明白了。 |
huangyoyo 发表于 2019-3-21 09:37 我最近正好也在看ORCA的ESD模块,下面是我的看法: 在社长写的“振动分辨的电子光谱的计算”(http://sobereva.com/223)提到了FC近似、HT近似,考虑了HT效应,就是考虑Vibronic coupling,是需要电子的波函数对核坐标的导数信息,所以这里需要提供一个Hessian。 因为我也是刚开始接触ORCA的ESD模块,以上是个人看法,如有错误,还请sob老师指出。 |
本帖最后由 huangyoyo 于 2019-3-21 10:01 编辑 sobereva 发表于 2019-3-20 16:51 再比如这个输入文件,仅仅考虑S0和S1的振动耦合?为什么要优化所有的激发态? |
QQ截图20190321092511.png (84.88 KB, 下载次数 Times of downloads: 49)
huangyoyo 发表于 2019-3-20 15:12 前面说了,不考虑振动耦合就根本不需要hess |
本帖最后由 huangyoyo 于 2019-3-20 15:14 编辑 sobereva 发表于 2019-3-20 14:46 这是那一节的所有信息,hess 文件用于计算吸收速率?如果问我仅仅需要吸收谱,只需要用tddft就可以了。 |
QQ截图20190320151051.png (179.65 KB, 下载次数 Times of downloads: 55)
QQ截图20190320151156.png (221.53 KB, 下载次数 Times of downloads: 53)
huangyoyo 发表于 2019-3-20 14:31 你没有正确理解这个例子的意图 |
sobereva 发表于 2019-3-20 11:28 我是看他手册里那个ESD模块的输入文件中要求的。 |
QQ截图20190320143037.png (78.16 KB, 下载次数 Times of downloads: 45)
huangyoyo 发表于 2019-3-20 09:09 “ORCA中计算吸收谱需要激发态的hess文件”这是错误认识。只有计算振动分辨的电子光谱才需要激发态的hess 只对某个结构计算吸收谱而且不考虑振动耦合的话,直接做TDDFT,然后用Multiwfn绘图就完了,见 使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱 http://sobereva.com/462(http://bbs.keinsci.com/thread-12114-1-1.html) 使用Multiwfn绘制红外、拉曼、UV-Vis、ECD、VCD和ROA光谱图 http://sobereva.com/224 可以对动力学过程的所有帧(或有间隔地提取)计算光谱,然后求平均。这在Multiwfn里相当于绘制构象权重平均光谱,只不过权重都相同。 |
sobereva 发表于 2019-3-19 18:29 恩,谢谢! 另外我想用NPT系综做分子动力学再将结果在ORCA中计算吸收谱。 首先我的结果算频率分析得到的基态hess文件是有虚频的(因为T和p都较高),其次,ORCA中计算吸收谱需要激发态的hess文件,而您之前提过我的结果是不能进行激发态几何优化的(因为激发态几何优化是在真空下进行的,温度和压力的影响不能反映出来),综合这两个因素,我想问NPT之后的结果能在ORCA中计算吸收谱么?若不能,可否推荐一个可行的方法? |
| 一般情况都不建议用RIJK,限制多,体系稍微大点效率就显著低于RIJCOSX |
sobereva 发表于 2019-3-19 17:15 换成RICOSX就能算了 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-19 08:43 , Processed in 0.178328 second(s), 25 queries , Gzip On.