计算化学公社

标题: 请老师看看我的晶体过渡金属配合物分子算吸收光谱和荧光发射光谱方案有无毛病 [打印本页]

作者
Author: Priscilla 时间: 2023-9-23 20:17
标题: 请老师看看我的晶体过渡金属配合物分子算吸收光谱和荧光发射光谱方案有无毛病
本帖最后由 Priscilla 于 2023-9-23 20:21 编辑

如题，之前关于这个问题的每个细节都已经发帖求助过了，但是还是有些疑问，现在列出本菜鸟的最后方案。
第一个是算垂直激发能的UV-Vis光谱，第二个是算荧光光谱的，体系是在晶体环境中olex2解出来的，又参考sob老师的579号博文，我私以为重原子在晶体环境下应该已经是最优不会有太大变化，所以我没有对基态进行优化，只是用MOPAC的PM7进行了H的优化然后加入了用Multiwfn算出来的chelpg背景电荷，但是我并不是很确定这点有问题所以向各位老师请教。
另外就是请老师们看看我最终的两个输入文件，给一些建议，非常感谢！

补充：之前的帖子里，一直把吸收能当成了发射能误会了PBE0的精度，真是低级错误，现在换了个def2tzvp算，目前的精度（即第一个文件的结果）经过对比和实验值只差了5nm，因此选择PBE0继续算第二个文件

作者
Author: sobereva 时间: 2023-9-24 04:35
算垂直激发能完全没必要写density，你又不对激发态的弛豫密度感兴趣，写这个只会白增加几倍的耗时。
之前我就反复说过，激发态优化没法光用背景电荷表现周围分子环境，因为背景电荷位置是固定的，没法自发地响应被优化的那个分子的结构而变化。只有QM/MM或ONIOM(QM:MM)才能做到令周围分子结构自发地响应。另外，228个原子做算TDDFT的受力本身就是巨耗时的事，你还用def2-TZVP这么昂贵的基组，根本不可能算得动

作者
Author: Priscilla 时间: 2023-9-24 10:49
本帖最后由 Priscilla 于 2023-9-24 10:50 编辑

原来如此学到了

作者
Author: Priscilla 时间: 2023-9-24 10:50
本帖最后由 Priscilla 于 2023-9-24 11:46 编辑

sobereva 发表于 2023-9-24 04:35
算垂直激发能完全没必要写density，你又不对激发态的弛豫密度感兴趣，写这个只会白增加几倍的耗时。
之前 ...

感谢老师不吝赐教，本菜新完全不知道这些问题，再向您请教一下，弛豫密度是指激发态没有优化并不是最低能量状态下的密度吗，考察这部分的波函数density是没有意义吗？您说的QM/MM在这里具体是怎么操作的，可否详述或有相关博文？您在579博文说oniom做了也无法对波函数分析那又如何解决想分析波函数的问题呢？还有"228个原子做算TDDFT的受力"是指对S1的优化吗，def2svp试过了误差较大，我应该用什么基组比较好呢？或者说，在晶体环境应该怎么去优化结构呢？

作者
Author: sobereva 时间: 2023-9-24 12:19

Priscilla 发表于 2023-9-24 10:50
感谢老师不吝赐教，本菜新完全不知道这些问题，再向您请教一下，弛豫密度是指激发态没有优化并不是最低能 ...

看下文了解什么叫弛豫和非弛豫的激发态密度
使用Multiwfn做空穴-电子分析全面考察电子激发特征
http://sobereva.com/434（http://bbs.keinsci.com/thread-10775-1-1.html）

QM/MM或ONIOM(QM:MM)几句话说不清楚，细节超级多。北京科音基础量子化学培训班（http://www.keinsci.com/workshop/KBQC_content.html）专门有一节做详细讲解

Multiwfn可以对ONIOM(QM:MM)产生的波函数分析，formchk转换的时候特别注意就完了

显然优化的每一步都得算所有原子受力。误差较大怎么判断的说清楚，说不定完全是你的误解

作者
Author: Priscilla 时间: 2023-9-24 13:30

sobereva 发表于 2023-9-24 12:19
看下文了解什么叫弛豫和非弛豫的激发态密度
使用Multiwfn做空穴-电子分析全面考察电子激发特征
http:// ...

那请问老师为什么垂直激发能不需要density呢？因为我想探究一下激发后的波函数分布和电子转移，所以我写了density
另因为课题组不允许频繁外出培训，之前报了初级量化最近又报了cp2k，想自己报名参加但是手头并不宽裕所以基础班得等来年，问题又很紧张需要处理，我看老师在一篇帖子http://bbs.keinsci.com/thread-942-1-1.html下面放了一些oniom的资源，想请问自学这个可以解决当前的晶体激发态优化问题吗
我总是把和实验值的拟合程度作为衡量误差判断的标准，这次我觉得较大是因为换了大基组（尽管还是背景电荷的错误做法）后，带隙能量上升0.3ev的原因，不过之前您也回复过这并不是很大的误差。但是由于这次相比于垂直激发能在def2tzvp下和实验值几乎没有误差，在这种误差校正下，所以觉得的def2svp误差非常大。希望老师指教

作者
Author: sobereva 时间: 2023-9-24 13:40

Priscilla 发表于 2023-9-24 13:30
那请问老师为什么垂直激发能不需要density呢？因为我想探究一下激发后的波函数分布和电子转移，所以我写 ...

写不写density又不影响你算的激发能，你就说了算激发能，事先又不说还要研究激发态波函数，我当然不让你加density浪费时间
有特殊目的在提问的时候事先说清楚

搞清楚所有可能导致误差的来源，搞清楚什么叫误差的巧合抵消，别把所有和实验对比的误差都赖到某一个因素上

欢迎光临计算化学公社 (http://bbs.keinsci.com/)