计算化学公社

标题: 计算MOF的紫外吸收光谱，和实验数据差的有点离谱 [打印本页]

作者
Author: njust-lbc 时间: 2024-6-19 17:37
标题: 计算MOF的紫外吸收光谱，和实验数据差的有点离谱
本帖最后由 njust-lbc 于 2024-6-25 20:55 编辑

请问各位老师，我对某MOF进行TDDFT计算并希望得到其紫外吸收光谱。但是得到的数据感觉和实验数据差的有点离谱，不知道问题出现在哪里，希望得到老师们的指点。
1. 我先将实验测得的cif中所有重原子固定，在PBE/DZVP-MOLOPT-SR-GTH条件下对氢原子做结构优化；
2. 对上一步得到的结构在PBE0/DZVP-MOLOPT-SR-GTH的条件下做TDDFT计算，读取上一步的波函数文件，考虑所有空轨道，只算十个激发态。
3. 这里面锌原子是+2价的，所以MOF中没有未成对的电子，体系的自选多重度设的是1
结果算出来发现第一激发能是负值，改用PBE/DZVP-MOLOPT-SR-GTH条件做TDDFT计算，第一激发能只有0.00924 eV，实验上紫外测得的光学带隙有2.4 eV，明显算出来的结果是有问题的。
是因为我选的基组有问题么，还是我做TDDFT计算的cif有问题， :(
PS：我已经上传了优化后的MOF cif文件，还有我用PBE0做TDDFT的输入输出文件，希望老师们不吝赐教！！

作者
Author: sobereva 时间: 2024-6-19 23:23
跟基组没关系
没准OT收敛到了不稳定波函数。改用对角化再试。并且检查结构，我发现扩胞后在原晶胞边缘处能明显看到位置不合理的原子

作者
Author: njust-lbc 时间: 2024-6-20 20:26
本帖最后由 njust-lbc 于 2024-6-20 20:42 编辑

sobereva 发表于 2024-6-19 23:23
跟基组没关系
没准OT收敛到了不稳定波函数。改用对角化再试。并且检查结构，我发现扩胞后在原晶胞边缘处能 ...

请问老师，我白天试了一下用对角化优化氢原子位置，结果发现SCF根本没有收敛的趋势，Convergence一直在0.几和1.几之间来回震荡，这种是什么问题呢？明明用OT是可以收敛的啊。我读取OT的波函数做初猜也不行

是因为我用的2022.1版本的原因么？

cif是测单晶如何解出来的，结构的话我也看了，我没太看出来哪里有不合理的地方，，，

这个不合理是指原子间的距离太近吗？

作者
Author: wzkchem5 时间: 2024-6-21 00:57
检查该补氢原子的地方是不是都补氢原子了。不能盲目相信cif文件补氢的结果，需要根据化学常识手动检查

作者
Author: sobereva 时间: 2024-6-21 06:39

njust-lbc 发表于 2024-6-20 20:26
请问老师，我白天试了一下用对角化优化氢原子位置，结果发现SCF根本没有收敛的趋势，Convergence一直在0. ...

对角化本来就往往比OT难收敛

跟版本没关系

扩成超胞，检查不合理接触

作者
Author: njust-lbc 时间: 2024-6-21 08:56

wzkchem5 发表于 2024-6-21 00:57
检查该补氢原子的地方是不是都补氢原子了。不能盲目相信cif文件补氢的结果，需要根据化学常识手动检查

好的，谢谢老师

作者
Author: njust-lbc 时间: 2024-6-21 08:56

sobereva 发表于 2024-6-21 06:39
对角化本来就往往比OT难收敛

跟版本没关系

好的，谢谢老师

作者
Author: njust-lbc 时间: 2024-6-25 20:54
该楼的问题已经解决，cif中原来有二甲胺阳离子我以为没什么用就给删了，加回来后结构优化和后续的TDDFT计算就正常了。。。

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