计算化学公社

标题: 求助片段组合波函数、自旋多重度的设置 [打印本页]

作者
Author: ZH123456 时间: 2022-9-14 17:55
标题: 求助片段组合波函数、自旋多重度的设置
我想采用片段组合波函数计算金属配合物（二价体系，自旋多重度为6）过渡态产生自由基片段，输入文件第一部分设置0，2 a，第二部分设置2,5,a,输出文件总是第一部分0 2 a 第二部分总是0 1 a，请问怎么解决这个问题呢？片段组合波函数可以计算多重度的带电体系吗？

作者
Author: zjxitcc 时间: 2022-9-15 00:33
本帖最后由 zjxitcc 于 2022-9-15 00:37 编辑

我猜你是用GaussView划分片段，然后保存gjf再打开发现不是你想要的片段划分方式。GaussView划分片段是一种辅助手段，但不是必须手段，如果你觉得你的划分是合理的，直接使用文本编辑器（如Sublime Text3, UltraEdit）打开gjf文件编辑片段电荷和自旋多重度即可。例如这个帖子里11L我的回复，[Fe(H2O)6]^(3+)片段组合波函数输入文件我就是手动添加片段信息的。

作者
Author: sobereva 时间: 2022-9-15 05:14
总要把问题交代清楚，最起码让别人知道拿什么程序算的。说不清楚就上传输入文件

a指什么说清楚

“片段组合波函数可以计算多重度的带电体系吗？”明显是病句

把此文看了
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

作者
Author: ZH123456 时间: 2022-9-15 09:50

zjxitcc 发表于 2022-9-15 00:33
我猜你是用GaussView划分片段，然后保存gjf再打开发现不是你想要的片段划分方式。GaussView划分片段是一种 ...

谢谢您，我的体系是2+，自旋多重度为6，我要设计一部分生成自由基，用gaussview设置好片段（0 2 a 和 2 3 a），提交Linux系统计算完成后查看log文件，用GaussView打开查看了atom groups 发现不是我当时设计的电荷和自旋多重度，变成了0 2 a 和 0 1 a，怎么计算这个体系肯定不会变成0价，不知道应该加什么关键词解决这个问题，这方面的介绍也很少

作者
Author: ZH123456 时间: 2022-9-15 10:15
本帖最后由 ZH123456 于 2022-9-19 09:18 编辑

sobereva 发表于 2022-9-15 05:14
总要把问题交代清楚，最起码让别人知道拿什么程序算的。说不清楚就上传输入文件

a指什么说清楚

拿GaussView设置的atom group 然后保存gjf文件提交在linux系统计算的，a指Alpha电子，我看这个文章写的是中性体系自旋多重度1的，我要算的是2+体系自旋多重度为4，分层设置然后计算过渡态产生自由基片段，所以我不知道可以算不，这方面的介绍很少，下面是我输入文件，log文件查看是2 4 0 2 0 1，感谢您的解答，新手这方面有很多不懂得
#p opt=(calcfc,ts,noeigentest,gdiis,maxcycle=300) freq ub3lyp/6-31+g(d,p)
scrf=(cpcm,solvent=water) nosymm scf=qc int=ultrafine

TS-1-3

2 4 0 2 2 3
H(Fragment=1)    000
H(Fragment=1)    000
S(Fragment=2) 000
O(Fragment=2) 000
O(Fragment=2)    000
F(Fragment=1) 000
F(Fragment=1) 000
F(Fragment=1) 000
F(Fragment=1) 000
O(Fragment=2) 000
Fe(Fragment=2)  000
F(Fragment=1) 000

作者
Author: zjxitcc 时间: 2022-9-15 12:35

ZH123456 发表于 2022-9-15 10:15
拿GaussView设置的atom group 然后保存gjf文件提交在linux系统计算的，a指Alpha电子，我看这个文章写的是 ...

这里面问题很多，有比“log文件查看是2 4 0 2 0 1”这种次要问题重要的多得多的问题：

（1）关键词那行没写guess(fragment=片段数目)，底下无论怎么写都没用，都不是片段组合波函数，再次强烈建议反复、仔细阅读《谈谈片段组合波函数与自旋极化单重态》http://sobereva.com/82

（2）一般来讲不建议使用scf=qc，除非你很清楚自己在干什么。如果你是想防止scf不收敛，scf=qc是优先级十分靠后的一种手段，绝不应该作为首选。

（3）过渡金属用B3LYP不是明智的行为，强烈建议阅读《简谈量子化学计算中DFT泛函的选择》http://sobereva.com/272

（4）复杂体系应先正确计算单点、检验波函数稳定性，然后再读取轨道进行结构优化。强烈建议阅读这个帖子11L的回复http://bbs.keinsci.com/thread-31621-1-1.html

作者
Author: ZH123456 时间: 2022-9-16 21:53

zjxitcc 发表于 2022-9-15 12:35
这里面问题很多，有比“log文件查看是2 4 0 2 0 1”这种次要问题重要的多得多的问题：

（1）关键词那 ...

谢谢老师的指导

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