计算化学公社

标题: 关于sob-MECP计算问题的咨询，谢谢指点 [打印本页]

作者
Author: 小范范1989 时间: 2017-9-18 09:01
标题: 关于sob-MECP计算问题的咨询，谢谢指点
利用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步才说收敛，我猜，是我是我哪里弄错了？输入文件还是什么？
谢谢老师指点
具体的计算输出的文件我已经上传，请查看，谢谢

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

作者
Author: 小范范1989 时间: 2017-9-18 15:09

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

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

作者
Author: sobereva 时间: 2017-9-18 19:34

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

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

作者
Author: 小范范1989 时间: 2017-9-18 20:24

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

恩，好的，谢谢sob老师

作者
Author: beefly 时间: 2017-11-18 21:11
可能是Gaussian版本的问题。据测试，用G09.A做DFT正常，一旦用到TDDFT，结构/能量就固定不变。G09.D，G09.E没有这个问题。

作者
Author: 小范范1989 时间: 2017-11-19 07:31

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

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

作者
Author: kaegi 时间: 2023-10-31 21:25
需要新建一个读取TD文件能量的脚本，比如取名为extract energyTD，仿造原来读取能量文件的格式，修改为读取TD计算的激发态能量，同时要修改sub_script，调用对应的extract energyTD文件。

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