计算大体系的ESD(FLUO)的算力简化问题

JefferSon

老师好，
我目前在做一个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:

wzkchem5

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

JefferSon

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得到的图谱是一个在整体范围内单调减的曲线。
谢谢老师，祝工作顺利！

wzkchem5

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关键词去掉一些低频振动模式，具体用法看说明书

JefferSon

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

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