使用orca做tddft/post-HF如何产生对应的自然轨道 |
上善若水 发表于 2022-8-8 21:47 你可以说吸收了某个频率的光,激发到某激发态后,键被削弱了 |
卢老师,请问您:由本文中方法与基态中所计算的结果相对比得出的结论,如键级的降低,可以描述为分子(基态)在某波长光(跃迁到某个态)下的响应吗? |
Novice 发表于 2020-11-6 14:44 明明是mwfn格式有效位数更多,所以说mwfn是记录波函数最完美的格式 差异更大可能是来自于产生自然轨道的细节上。这点差异完全可以忽略。 |
sobereva 发表于 2020-11-6 03:33 Sob老师,我刚才迫不及待试了下,发现分析ADCH电荷时基于gaussian直接输出的wfn文件计算出的Ir的hirsfeld电荷为0.28,而基于mwfn文件计算出的hirsfeld电荷为0.30。其它原子的差别却小一些(大概是在小数点后四位以后)。最终的ADCH电荷也有一定不同。是因为mwfn文件的数字有效位数较少,还是有什么bug? |
sobereva 发表于 2020-11-6 02:27 谢谢sob老师,我下载最新的Multiwfn软件进行分析! |
wzkchem5 发表于 2020-11-5 22:19 好的,我检查下,谢谢您! |
Novice 发表于 2020-11-3 16:48 目前最新版本已改成输出mwfn格式了 还有一些原本输出molden文件的功能我也都改成输出mwfn格式了,完全杜绝了这个问题 |
参与人数Participants 1 | eV +5 | 收起 理由Reason |
---|---|---|
| + 5 | 赞! |
陈丢丢 发表于 2020-11-5 21:17 今天我更新了官网上的Multiwfn最新版本,此功能输出.mwfn格式代替了老版本的.molden格式,完全没这个问题了。 Gaussian计算时候写Pop=spinnatural是完全多余的 |
陈丢丢 发表于 2020-11-5 21:17 检查一下molden文件的轨道占据数对不对? |
本帖最后由 陈丢丢 于 2020-11-5 21:24 编辑 sobereva 发表于 2020-11-4 14:46 sob老师好,我现在算一个金属体系,fe、K用的LANL2DZ,其他原子用的6-31G*。目前想对这个体系做自旋自然轨道分析,同样遇到了“molden文件不记录原子的实际核电荷数信息”这个问题,我修改了molden文件里面fe和K的核电子数,将修改后的molden文件重新拉入Multiwfn中,但是还是报错,是因为我修改数据错误吗?我将fe和K的核电子数改成16和9C:\Users\dell\Desktop\计算SNO报错.png。此外如果我用高斯算自旋自然轨道,产生wfn文件,用wfn文件进行SNO分析的话,是用Multiwfn直接看轨道吧?用的关键词是: ub3lyp/genecp Pop=spinnatural out=wfn geom=connectivity empiricaldispersion=gd3 ***** ***** ts5-p-t-k-top-sno.wfn |
计算SNO报错.png (198.9 KB, 下载次数 Times of downloads: 144)
sobereva 发表于 2020-11-4 14:46 明白了,多谢社长回复 ![]() |
Novice 发表于 2020-11-3 16:48 molden文件不记录原子的实际核电荷数信息(原子序数扣除被赝化的电子数),所以重新载入molden后程序就不知道实际原子的核电荷数是多少了。解决办法在Multiwfn手册2.5节的molden格式介绍部分提了,需要自己改一下molden文件。 这是个molden格式的缺陷,而非Multiwfn的bug。molden格式的这个不足在后来我提出的mwfn格式(http://bbs.keinsci.com/thread-16363-1-1.html)中避免了,我以后会让这个功能导出mwfn格式而非molden格式。 实际上,对于你当前的分析,由于不需要基函数信息,还不如直接用Gaussian产生的wfn/wfx文件当Multiwfn输入文件,明显更省事。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-12 23:28 , Processed in 0.522689 second(s), 26 queries , Gzip On.