简单量子化学问题答疑专帖

sobereva · 发表于 Post on 2023-5-15 09:13:56

奶啵啵发表于 2023-5-15 08:49
#p pbe1pbe/genecp em=gd3bj int=fine freq scrf 计算的为吸附过程，opt和freq最开始一起计算，OPT不收敛, ...

内存给得太少
仔细看
Gaussian的安装方法及运行时的相关问题
http://sobereva.com/439（http://bbs.keinsci.com/thread-10814-1-1.html）

sobereva · 发表于 Post on 2023-5-15 09:14:34

cokie 发表于 2023-5-15 08:21
请问各位老师，目前除了帅老师的MOMAP软件外，还有没有其他途径可以较准确的计算Kisc和Krisc呢？谢谢！

ORCA可以计算kisc，FCclasses可以计算kisc，我在论坛里都回复过。都不要钱

化小白 · 发表于 Post on 2023-5-16 12:25:19

老师，我想找一个五元环成环的过渡态，环内含过渡金属，想要的是两个碳成键。我优化所用的是b3lyp/6-31G(d)-SDD,PCM溶剂模型，加了色散矫正，要不一直震荡，要不就是成反应物，要不就是成产物。我就把震荡的最后结构做了振动分析，也确实是想要的振动方向，做IRC分析，两边都走到了产物。然后我就用M06/def2-SVP优化震荡的结构，直接就变成了产物；又用了b3lyp/6-31G(d,p),震荡；又用了b3lyp/6-311G(d,p)震荡；又用了b3lyp/6-31+G(d)震荡。那个不收敛的博文也都试了很多种方法，各种坐标下优化，放宽收敛限，maxstep=3,5都试了。想问一下老师，还有没有其他办法让它收敛？

sxh3231796 · 发表于 Post on 2023-5-16 16:04:15

sobereva 发表于 2023-5-13 23:01
Gaussian直接就会给出各个态中的组态函数的系数，求平方乘100%就是贡献量

谢谢sob老师，看到文献里写的计算的组态系数,分成使用natural orbital和localized orbitals，请问这种结果是要在高斯关键词中加些特定关键词吗？还是使用后续其他软件分析结果得到的呀？还有后面这个2200的电子占据轨道的结果是直接在高斯计算结果中读取的吗？

wzkchem5 · 发表于 Post on 2023-5-16 17:41:35

化小白发表于 2023-5-16 05:25
老师，我想找一个五元环成环的过渡态，环内含过渡金属，想要的是两个碳成键。我优化所用的是b3lyp/6-31G(d) ...

如果IRC两边都走到产物，说明当前结构离产物太近，把结构往反应物方向调整一下，再优化

化小白 · 发表于 Post on 2023-5-16 19:14:09

wzkchem5 发表于 2023-5-16 17:41
如果IRC两边都走到产物，说明当前结构离产物太近，把结构往反应物方向调整一下，再优化

老师，有一个过渡态结构像反应物的做了IRC之后也是一边走向了产物，另一边就走了一个点就结束了，我感觉那个过渡态也不对

sobereva · 发表于 Post on 2023-5-16 20:56:01

sxh3231796 发表于 2023-5-16 16:04
谢谢sob老师，看到文献里写的计算的组态系数,分成使用natural orbital和localized orbitals，请问这种结 ...

Gaussian输出文件里有各个组态函数的定义，以及各自的系数

文献里这个应该估计不是标准的CASSCF，怀疑给的是基于不同轨道做CASCI得到的系数，没考虑轨道优化

sxh3231796 · 发表于 Post on 2023-5-16 21:17:38

sobereva 发表于 2023-5-16 20:56
Gaussian输出文件里有各个组态函数的定义，以及各自的系数

文献里这个应该估计不是标准的CASSCF，怀疑 ...

谢谢sob老师，https://onlinelibrary.wiley.com/doi/10.1002/anie.201309767之前看过的就是这篇文献，这种系数用高斯可以算吗？还是需要用其他软件，类似于mokit之类的呀？

sobereva · 发表于 Post on 2023-5-16 21:25:08

sxh3231796 发表于 2023-5-16 21:17
谢谢sob老师，https://onlinelibrary.wiley.com/doi/10.1002/anie.201309767之前看过的就是这篇文献，这 ...

如果你就是要做CASCI，用iop(5/7=-1)，此时不做轨道优化

自然轨道和定域化轨道都可以用Multiwfn产生，参考下文，之后用主功能100的子功能2导出成fch，再用unfchk转成chk，就可以在CASSCF计算时用guess=read当初猜轨道
Multiwfn的轨道定域化功能的使用以及与NBO、AdNDP分析的对比
http://sobereva.com/380（http://bbs.keinsci.com/thread-6053-1-1.html）
在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例
http://sobereva.com/403（http://bbs.keinsci.com/thread-9001-1-1.html）

用作CASSCF初猜的自然轨道用Gaussian自己的功能产生也可以，看下文相关部分
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379（http://bbs.keinsci.com/thread-6020-1-1.html）

sxh3231796 · 发表于 Post on 2023-5-16 22:00:27

sobereva 发表于 2023-5-16 21:25
如果你就是要做CASCI，用iop(5/7=-1)，此时不做轨道优化

自然轨道和定域化轨道都可以用Multiwfn产生， ...

多谢sob老师，之前在做完casscf之后就用Multiwfn做了轨道定域化分析，查看轨道都没有问题，但是没有文献里显示的相关组态系数，是做完casci在高斯结果中读取吗？还是Multiwfn也有相应的功能啊？

sobereva · 发表于 Post on 2023-5-16 22:27:03

sxh3231796 发表于 2023-5-16 22:00
多谢sob老师，之前在做完casscf之后就用Multiwfn做了轨道定域化分析，查看轨道都没有问题，但是没有文献 ...

Multiwfn基于CASSCF任务产生的波函数文件做轨道定域化没有丝毫意义。轨道定域化只能用于单行列式波函数。
前面说了，当前情况下定域化轨道是作为CASSCF的初猜轨道用的

sxh3231796 · 发表于 Post on 2023-5-16 22:39:05

sobereva 发表于 2023-5-16 22:27
Multiwfn基于CASSCF任务产生的波函数文件做轨道定域化没有丝毫意义。轨道定域化只能用于单行列式波函数。 ...

哦哦哦，老师的意思是说，做常规的结构优化得到轨道波函数，用Multiwfn做轨道定域化分析生成新的波函数，以此作为初猜再进行计算casscf吗？还是做casci的计算呢？

sxh3231796 · 发表于 Post on 2023-5-16 23:08:31

wzkchem5 发表于 2023-5-4 20:14
一般来说，研究过渡金属配合物，应该把过渡金属的价层d轨道全部选进来。因为即便有些d轨道不参与成键，这 ...

还想再跟老师请教一下，怎么判断那个轨道时金属的价层d轨道，从轨道形状来看吗？还是高斯计算之后会有明确的定义呢？

wzkchem5 · 发表于 Post on 2023-5-16 23:45:20

sxh3231796 发表于 2023-5-16 16:08
还想再跟老师请教一下，怎么判断那个轨道时金属的价层d轨道，从轨道形状来看吗？还是高斯计算之后会有明 ...

用轨道形状或轨道成分（http://sobereva.com/131）

sobereva · 发表于 Post on 2023-5-17 10:48:36

sxh3231796 发表于 2023-5-16 22:39
哦哦哦，老师的意思是说，做常规的结构优化得到轨道波函数，用Multiwfn做轨道定域化分析生成新的波函数， ...

是
对Gaussian而言是CASSCF任务，但通过IOp让Gaussian实际做的是CASCI

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

[综合交流] 简单量子化学问题答疑专帖