关于sob-MECP计算问题的咨询，谢谢指点

小范范1989

利用sob-MECP寻找S1-T2的交叉点，最后runMECP，看着也正常结束了，但是执行到了99步，然后ReportFile中能量都一样。比如说：

Energy of First State:    -2219.5971489500
Energy of Second State:   -2219.5971489500

Convergence Check (Actual Value, then Threshold, then Status):
Max Gradient El.:   0.052833 (0.000700)   NO
RMS Gradient El.:   0.008798 (0.000500)   NO
Max Change of X:    0.036983 (0.004000)   NO
RMS Change of X:    0.006158 (0.002500)   NO
Difference in E:    0.000000 (0.000050)  YES

但是我看例子中的，都是能量最后才几乎接近。我这里ReportFile中的能量时刻都是一样的。然后最后99步才说收敛，我猜，是我是我哪里弄错了？输入文件还是什么？
谢谢老师指点
具体的计算输出的文件我已经上传，请查看，谢谢

sobereva

应该是脚本对数据解析得不对。不知道为什么你的A计算用TD，B计算用TDA。你看看reportfile里解析出来的能量到底对应于输出文件里的哪项，如果有误可以自行改改脚本。

小范范1989

sobereva 发表于 2017-9-18 11:30
应该是脚本对数据解析得不对。不知道为什么你的A计算用TD，B计算用TDA。你看看reportfile里解析出来的能量 ...

谢谢sob老师，因为我计算单重态用的td，三重态用的tda。
我查阅了这个reportfile中的能量，对应的是高斯计算log输出文件中文件最后HF=.....这块能量。也就是基态势能面上的能量呗，所以，计算出来reportfile中两个能量一直一样。
老师说的“如果有误可以自行改改脚本。”修改哪个脚本呢？如何修改呢？谢谢老师指点

sobereva

小范范1989 发表于 2017-9-18 15:09
谢谢sob老师，因为我计算单重态用的td，三重态用的tda。
我查阅了这个reportfile中的能量，对应的是高斯 ...

sobMECP源码没多少，仔细看一遍就知道了

小范范1989

sobereva 发表于 2017-9-18 19:34
sobMECP源码没多少，仔细看一遍就知道了

恩，好的，谢谢sob老师

beefly

可能是Gaussian版本的问题。据测试，用G09.A做DFT正常，一旦用到TDDFT，结构/能量就固定不变。G09.D，G09.E没有这个问题。

小范范1989

beefly 发表于 2017-11-18 21:11
可能是Gaussian版本的问题。据测试，用G09.A做DFT正常，一旦用到TDDFT，结构/能量就固定不变。G09.D，G09.E ...

偶，还有这一点，谢谢beefly大神指点。

kaegi

需要新建一个读取TD文件能量的脚本，比如取名为extract energyTD，仿造原来读取能量文件的格式，修改为读取TD计算的激发态能量，同时要修改sub_script，调用对应的extract energyTD文件。