计算化学公社

标题: 铁配合物片段初猜设置求助 [打印本页]

作者
Author: ZH123456 时间: 2022-9-17 22:56
标题: 铁配合物片段初猜设置求助
本帖最后由 ZH123456 于 2022-9-19 09:20 编辑

我采用片段初猜的方法设置了铁配合物断键来产生自由基片段，设置两部分如下 0 2 a 和 2 3 a （第一个片段带一个电子）输入文件如下。计算结果文件查看Mulliken charges and spin densities condensed to fragments:    1    0.272564    -0.552047
   2    1.727436    3.552047
发现第一部分并不是携带一个电子，请问这个问题该怎么解决啊？我是设置断键产生自由基，那么我的初设结构键长应该是按照断键的要求设置吗？第一步得到稳定的波函数再进行过渡态计算吗？请老师指教，谢谢老师！
输入文件：

%chk=R-Fe3+.chk
%mem=4GB
%nprocshared=40
#p ub3lyp/6-31+g(d,p) nosymm guess(fragment=2) scf=qc

E-R-Fe3+

2 4 0 2 2 3
H(Fragment=1) 0             0             0
H(Fragment=1)    0             0             0
C(Fragment=2)    - 0             0             0
O(Fragment=2)    0             0             0
O(Fragment=2) 0             0             0
H(Fragment=1)    0             0             0
H(Fragment=1) 0             0             0
H(Fragment=1)    0             0             0
H(Fragment=1)    0             0             0
O(Fragment=2)    0             0             0
Fe(Fragment=2) 0             0             0
H(Fragment=1)    0             0             0

--Link1--
%chk=E-R-Fe3+.chk
%mem=4GB
%nprocshared=40
#p ub3lyp chkbasis guess=read geom=allcheck scrf=(cpcm,solvent=water) nosymm stable=opt
scf(xqc,maxcycle=128)

作者
Author: wzkchem5 时间: 2022-9-18 03:14
参见http://sobereva.com/82最后一节，计算出来的自旋布居本来就不应该精确等于设定的单电子数

作者
Author: sobereva 时间: 2022-9-18 04:54
片段初猜只能设置初猜波函数的状态，没法决定收敛完了什么样

没事甭乱加弥散函数，看
谈谈弥散函数和“月份”基组
http://sobereva.com/119
而且带弥散函数的时候Mulliken分析结果有误导性

没事甭写scf=qc，正常情况下能收敛的情况时候这又浪费时间又增加收敛到不稳定波函数的概率

没事甭用cpcm，仔细看
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）

不要做出给几十核却只给几G内存这种事，仔细看
Gaussian的安装方法及运行时的相关问题
http://sobereva.com/439（http://bbs.keinsci.com/thread-10814-1-1.html）

上传完整输入文件或者给出结构图别人才能给出更确切的回复

发帖时标题必须具体，这次给你改了，别就写个“片段初猜”，以后注意。仔细阅读http://bbs.keinsci.com/thread-9348-1-1.html

作者
Author: ZH123456 时间: 2022-9-18 20:40

sobereva 发表于 2022-9-18 04:54
片段初猜只能设置初猜波函数的状态，没法决定收敛完了什么样

没事甭乱加弥散函数，看

明白了，谢谢老师的指导，我会认真阅读的

作者
Author: ZH123456 时间: 2022-9-18 20:41

wzkchem5 发表于 2022-9-18 03:14
参见http://sobereva.com/82最后一节，计算出来的自旋布居本来就不应该精确等于设定的单电子数

明白了，谢谢老师的指导

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