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

kaiden

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

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

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

snljty

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

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

kaiden

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

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

zjxitcc

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

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

zjxitcc

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

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

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

kaiden

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

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

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

ljh123

大佬您好，我想请问下，如果单电子转移反应前后自旋多重度发生变化了的话，该怎么计算电子转移积分呢？

Warm_Cloud

ljh123 发表于 2024-11-27 19:46
大佬您好，我想请问下，如果单电子转移反应前后自旋多重度发生变化了的话，该怎么计算电子转移积分呢？

Amesp可以直接算不同自旋多重度的电子转移积分。

ljh123

Warm_Cloud 发表于 2024-11-30 08:52
Amesp可以直接算不同自旋多重度的电子转移积分。

大佬我有两个问题想咨询下您，（1）我想算的是自旋多重度5和自旋多重度7的电子转移积分。我之前了解到的似乎都是根据TDDFT计算，并且只能计算基态和三重态。所以我不太清楚我这种情况是否能算。（2）我的体系是周期性体系，总共有400多个原子，1个Fe，其他都是C,H,O,N。如果可以计算电子转移积分的画，不清楚这种体系能不能算得动

cokie

ljh123 发表于 2024-11-30 14:09
大佬我有两个问题想咨询下您，（1）我想算的是自旋多重度5和自旋多重度7的电子转移积分。我之前了解到的 ...

Worm_Cloud 大佬开发的Amesp算转移积分很好用。你的第一个问题，按Amesp手册，对两个不同多重度片段定义好自旋多重度即可。

(2) “周期性体系” 和 “总共400个原子” 的描述似乎有点矛盾。
电荷转移积分的程序可以算【你感兴趣的片段】之间的转移积分。
比如周期性体系，你可以把感兴趣的部分拆分出两个能发生电荷转移的重复单元片段，其他不感兴趣的部分可以删掉并用H把边缘做饱和即可。

ljh123

cokie 发表于 2024-11-30 21:23
Worm_Cloud 大佬开发的Amesp算转移积分很好用。你的第一个问题，按Amesp手册，对两个不同多重度片段定义 ...

关于周期性和400个原子的表述，因为周期性是用pbc盒子表现的，实际盒子里就是400个原子，相当于一个单元。可以请问下用H把边缘做饱和是怎么操作吗，我的体系是水溶液和溶质小分子。最后就是想请问下耗时问题，因为我可能要算几百个构象下的转移积分，希望能用便宜一些的方法计算。
此外，我想问一下如果把这种自旋多重度会发生变化的过程的电子转移积分，简化为同一个自旋多重度的电子转移积分，这样是否可以（这样简化就可以直接在CP2K里算，速度快）
十分感谢大佬

cokie

ljh123 发表于 2024-11-30 21:50
关于周期性和400个原子的表述，因为周期性是用pbc盒子表现的，实际盒子里就是400个原子，相当于一个单元 ...

我有看到你23L的表述，事实上发生电子转移时(前)，自旋多重度是一样的对吧？自旋多重度发生变化是指两个产物，而这两个产物之间并不发生电子转移。

如果是的话，我个人认为算电荷转移积分应该使用反应前的自旋多重度和几何构型（或者说更准确一点是发生反应那一时刻的自旋多重度和几何构型）

另外至于体系应该怎么做具体处理，我不是很明白，因为这个表述让我还是有点迷糊

ljh123

cokie 发表于 2024-11-30 21:58
我有看到你23L的表述，事实上发生电子转移时(前)，自旋多重度是一样的对吧？自旋多重度发生变化是指两个 ...

反应是：A+B→(A+)+(B-)。其中反应物态（A+B）总自旋多重度为5，产物态总自旋多重度（(A+)+(B-)）为7。我想这应该是您说的意思？
关于体系，我不清楚大佬有没有用过CP2K做周期性体系的AIMD，就比如一个溶液体系，我们不可能把所有的水分子和溶质都建模出来，于是就可以像材料的晶胞一样，只建模一个重复性单元，并用pbc周期性边界条件表达出整体的效果。而这一个重复性单元中，一共有400个原子（而实际溶液中，肯定不止这么点原子，而是400*N倍个原子）

cokie

ljh123 发表于 2024-11-30 22:20
反应是：A+B→(A+)+(B-)。其中反应物态（A+B）总自旋多重度为5，产物态总自旋多重度（(A+)+(B-)）为7。我 ...

我一看周期性体系这个表述，以为发生反应的物质是两个周期性的晶体或者聚合物。。其实你只是建了个盒子跑AIMD搜构象对吗？实际电子转移反应还是发生在溶质之间对吗？溶质是小分子对吗？

如果反应物是A和B的话，我个人认为你要算转移积分，应该是A和B之间的转移积分，而不是A+和B-之间，或者（A+B）和（(A+) + (B-)）之间的转移积分。