一个Fortran写的电荷转移积分程序

Warm_Cloud · 发表于 Post on 2019-8-23 22:01:26

本帖最后由 Warm_Cloud 于 2020-7-19 13:31 编辑

attach://20792.x电荷转移积分的计算是很多小伙伴的需求，这两天基于高斯的输出文件做了一个，细节可以参考http://bbs.keinsci.com/thread-11369-1-1.html
使用方法很简单，需要提供dimer和monomers的高斯输出文件和fchk文件，以此按照提示输入就可以。输入文件要加 nosymm iop(3/33=1)。
顺便说一下，以后更强大和好用的电荷转移功能会加入到我自己写的量化软件中。

CT.x (8.96 MB, 下载次数 Times of downloads: 197)
代码地址：https://github.com/zhangyingfeng1993/Charge-Transport-Integral
有问题的小伙伴可以留言讨论。

kaiden · 发表于 Post on 2022-4-26 10:43:31

zjxitcc 发表于 2022-4-26 10:37
不收敛总是相对好解决的《解决SCF不收敛问题的方法》http://sobereva.com/61

最怕的是保留线性相关基 ...

谢谢zjxitcc！我有尝试用sob这个博文的建议解决不收敛，但是暂时还没有解决成功。另外想问下用J_from_g03那个程序算的电荷转移积分会不会受基函数线性依赖的影响呢？（我测试了这个程序，输出文件显示有基函数线性依赖，但是用J_from_g03计算过程没有报错，不知道结果是否可信。）

zjxitcc · 发表于 Post on 2022-4-26 10:37:40

kaiden 发表于 2022-4-26 10:35
我有尝试过加IOp(3/32=2)让Gaussian保留线性相关的基函数，但是算单点总是不收敛。

不收敛总是相对好解决的《解决SCF不收敛问题的方法》http://sobereva.com/61

最怕的是保留线性相关基函数后算出的电子能量很离谱

zjxitcc · 发表于 Post on 2022-4-26 10:35:57

本帖最后由 zjxitcc 于 2022-4-26 10:39 编辑

Warm_Cloud 发表于 2022-4-26 09:35
加上IOp(3/32=2)，这样Gaussian就不会处理线性依赖的基函数了。

这样容易造成电子能量异常，后续分析结果可信度降低。如果只有几个线性相关基函数，估计问题不大。

kaiden · 发表于 Post on 2022-4-26 10:35:06

Warm_Cloud 发表于 2022-4-26 09:35
加上IOp(3/32=2)，这样Gaussian就不会处理线性依赖的基函数了。

我有尝试过加IOp(3/32=2)让Gaussian保留线性相关的基函数，但是算单点总是不收敛。

snljty · 发表于 Post on 2022-4-26 10:23:24

Warm_Cloud 发表于 2022-4-26 09:35
加上IOp(3/32=2)，这样Gaussian就不会处理线性依赖的基函数了。

大师，请问然后会不会出现线性依赖导致系数矩阵不可逆，LAPACK报错的的问题？

Warm_Cloud · 发表于 Post on 2022-4-26 09:35:59

kaiden 发表于 2022-4-25 16:17
感谢楼主回复！因为测试失败的文件太大被我删掉了。后来查看fchk文件发现是因为Number of basis function ...

加上IOp(3/32=2)，这样Gaussian就不会处理线性依赖的基函数了。

kaiden · 发表于 Post on 2022-4-25 16:17:56

Warm_Cloud 发表于 2022-4-24 20:28
文件上传一下我试试。

感谢楼主回复！因为测试失败的文件太大被我删掉了。后来查看fchk文件发现是因为Number of basis functions和Number of independent functions不一致（应该是算dimer单点出现基函数线性依赖）。算单点用的m062x/ma-SVP。后来用B3LYP-D3BJ/def2svp重新算的单点就可以了。不知道后者缺乏带弥散函数的组合算出来的结果能不能接受呢？谢谢楼主！

Warm_Cloud · 发表于 Post on 2022-4-24 20:28:19

kaiden 发表于 2022-4-24 12:15
楼主您好，最近我用您写的程序测试电荷转移积分遇到一点问题，请教下面报错会是什么问题呢？谢谢！（在m062 ...

文件上传一下我试试。

kaiden · 发表于 Post on 2022-4-24 12:15:05

楼主您好，最近我用您写的程序测试电荷转移积分遇到一点问题，请教下面报错会是什么问题呢？谢谢！（在m062x/ma-SVP下计算的dimer，monomer单点。dimer.out载入时出错，但是把monomer.out当作dimer载入又没问题。）

nput for Dimer:
  Please input the Gaussian log file :<e.g. out.log>
pPDM_Dimer1_SP.out

  Spin Multiplicity:          1
  Number of basis functions:  2560
  Readind Orbit...
forrtl: severe (59): list-directed I/O syntax error, unit 1101, file /root/Documents/Study/Charge-Transport-Integral-master/pPDM_Dimer1_SP.fchk
Image             PC             Routine          Line       Source
CT.x             000000000094B016  Unknown             Unknown  Unknown
CT.x             0000000000967727  Unknown             Unknown  Unknown
CT.x             0000000000965985  Unknown             Unknown  Unknown
CT.x             0000000000404E8B  Unknown             Unknown  Unknown
CT.x             0000000000401AE0  Unknown             Unknown  Unknown
CT.x             0000000000400D29  Unknown             Unknown  Unknown
CT.x             0000000000400CDE  Unknown             Unknown  Unknown
CT.x             00000000009DE5E9  Unknown             Unknown  Unknown
CT.x             0000000000400BBA  Unknown             Unknown  Unknown

tulingzhi · 发表于 Post on 2021-12-3 14:32:00

请问楼主，您计算转移积分所采用的理论方法是基于这篇 JACS,128(30):9882-9886吗？是属于能级劈裂法还是在位能修正法啊？我个人的理解是作者是在能级劈裂法的基础上进行修正？

Warm_Cloud · 发表于 Post on 2020-7-19 13:32:27

更新了一下输出。

Puying · 发表于 Post on 2019-11-30 21:50:47

太牛批了！

小白的暖暖 · 发表于 Post on 2019-10-28 10:02:51

Warm_Cloud 发表于 2019-10-27 13:42
Je

好的，谢谢您

Warm_Cloud · 发表于 Post on 2019-10-27 13:42:49

小白的暖暖发表于 2019-10-26 13:17
请问，这四项哪一个是我想要的电荷转移积分啊
J(meV) S
Je(meV) Je(kcal/mol)

Je

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[辅助/分析程序] 一个Fortran写的电荷转移积分程序

评分 Rate

评分 Rate