计算化学公社

标题: 计算大体系的ESD(FLUO)的算力简化问题 [打印本页]

作者
Author: JefferSon 时间: 2023-5-17 16:56
标题: 计算大体系的ESD(FLUO)的算力简化问题
老师好，
我目前在做一个200个原子体系的荧光激发振动分辨谱，我在计算的过程中发现在最后一步计算correlation function时其点的数量是100个原子的平方倍，所需时间太大了。我想咨询一下有没有减少这个点的计算方法，如用VG代替AH（计算点数会下降嘛，手册说大体系精度可能不佳），还有就是如何降低图谱分辨率。谢谢老师的答疑解惑，祝老师工作顺利！

      ***Everything is set, now computing the correlation function***

Homogeneous linewidth is: 75.00 cm-1
Inhomogeneous linewidth is: 200.00 cm-1
Number of points: 2147483648
Maximum time: 1061.77 fs
Spectral resolution: 5.00 cm-1
Temperature used: 298.15 K
Z value: 1.601387e-386
Cutoff for the correlation function: 1.00e-06
Adiabatic energy difference: 18272.94 cm-1
0-0 energy difference: 17658.33 cm-1
Reference transition dipole (x,y,z): (0.09621 0.00000),
                                    (-0.07548 0.00000),
                                    (-0.01533 0.00000)
Calculating correlation function:

作者
Author: wzkchem5 时间: 2023-5-17 17:05
VG代替AH只节省correlation function之前的计算步骤的时间
correlation function的点数可以用npoints关键词修改，具体做法在手册里搜。谱图分辨率用maxtime控制，根据傅里叶变换的基本理论易知，maxtime越短，分辨率越低

作者
Author: JefferSon 时间: 2023-5-21 07:27

wzkchem5 发表于 2023-5-17 17:05
VG代替AH只节省correlation function之前的计算步骤的时间
correlation function的点数可以用npoints关键 ...

老师您好，那如果我已经计算得到了ES和GS的hess文件，我下次计算可以直接进行最后一步的“Everything is set, now computing the correlation function”吗，因为这样好像能减少很多时间。
还有就是可能是因为我点设太少了，我看他的0-0峰应该是在550nm左右，但是从spectrum得到的图谱是一个在整体范围内单调减的曲线。
谢谢老师，祝工作顺利！

作者
Author: wzkchem5 时间: 2023-5-21 16:01

JefferSon 发表于 2023-5-21 00:27
老师您好，那如果我已经计算得到了ES和GS的hess文件，我下次计算可以直接进行最后一步的“Everything is ...

当不考虑H-T效应的时候：可以，虽然ESD前面的TDDFT计算可能不能不算，但是那个耗时不大。当考虑H-T效应的时候，我不确定能不能只算correlation function，原理上可以但我不确定程序是否支持，建议自己看一下说明书，说明书没写就是不能，如果那样的话建议在orca论坛的feature request版面提一个request。
假设波长范围正确的话，那么曲线整体单调递减不是因为点数太少，而是因为低频振动模式的非谐性效应太大导致correlation function发散。应该用tcutfreq关键词去掉一些低频振动模式，具体用法看说明书

作者
Author: JefferSon 时间: 2023-5-21 19:44

wzkchem5 发表于 2023-5-21 16:01
当不考虑H-T效应的时候：可以，虽然ESD前面的TDDFT计算可能不能不算，但是那个耗时不大。当考虑H-T效应的 ...

好的谢谢老师，我这就去了解一下，祝工作顺利！

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