计算化学公社

标题: Mg-Cl surface scan [打印本页]

作者
Author: fhh2626 时间: 2015-11-13 19:05
标题: Mg-Cl surface scan
我最近做了一个小测试，以Mg2+和Cl-的距离为反应坐标做了一个surface scan，结果发现不同方法所得的结果相差巨大，似乎无法用误差来解释，而且曲线形状也不一样，不知道这是什么原因呢？使用开闭壳层似乎没什么影响。忘大家赐教，谢谢！
(, 下载次数 Times of downloads: 57)
输入文件：
! ccsd(t) aug-cc-pvtz tightscf nopop
%pal nprocs 8
end
%paras R = 2, 8, 61
end
% maxcore 1000

* int 1 1
Mg 0 0 0 0 0 0
Cl 1 0 0 {R} 0 0
*

作者
Author: sobereva 时间: 2015-11-13 21:33
可以用高斯重复一下看看，如果也是这样，对比较远的点检验下波函数稳定性，没准收敛到了不同电子态

作者
Author: fhh2626 时间: 2015-11-13 21:46

sobereva 发表于 2015-11-13 21:33
可以用高斯重复一下看看，如果也是这样，对比较远的点检验下波函数稳定性，没准收敛到了不同电子态

用高斯重复过，结果差不多，奇怪的是对NaCl的计算结果却没问题

作者
Author: fhh2626 时间: 2015-11-13 22:47

sobereva 发表于 2015-11-13 21:33
可以用高斯重复一下看看，如果也是这样，对比较远的点检验下波函数稳定性，没准收敛到了不同电子态

感谢指教~我做了波函数稳定性测试，并且使用了开壳层和guess=mix，但是结果还是
The wavefunction has an internal instability.

使用stable=opt程序不能正常结束，请问怎么办呢？

作者
Author: sobereva 时间: 2015-11-14 00:48

fhh2626 发表于 2015-11-13 22:47
感谢指教~我做了波函数稳定性测试，并且使用了开壳层和guess=mix，但是结果还是
The wavefunction has a ...

不清楚stable=opt怎么不能结束，但应该是初猜波函数轨道占据方式问题，初猜不是基态排布，而不在于非限制性开壳层的使用上。用guess=alter调换调换初猜时对不同不可约表示轨道的占据方式，或许能收敛到稳定波函数通过stable检测。也注意看看HF计算时不同不可约表示占据轨道的数目和DFT时是否不同。

作者
Author: zyj19831206 时间: 2015-11-14 17:28
这个是IRC曲线吧？？

作者
Author: fhh2626 时间: 2015-11-14 18:11

sobereva 发表于 2015-11-14 00:48
不清楚stable=opt怎么不能结束，但应该是初猜波函数轨道占据方式问题，初猜不是基态排布，而不在于非限制 ...

多谢指教，我现在解决了波函数稳定性的问题了，不过加上stable=opt以后似乎就不能做scan了

看来只能写个脚本在每个位置都算单点了？

作者
Author: fhh2626 时间: 2015-11-14 18:11

zyj19831206 发表于 2015-11-14 17:28
这个是IRC曲线吧？？

不是啊。。。

作者
Author: sobereva 时间: 2015-11-14 18:43

fhh2626 发表于 2015-11-14 18:11
多谢指教，我现在解决了波函数稳定性的问题了，不过加上stable=opt以后似乎就不能做scan了

看来只能写 ...

这个scan可能会涉及势能面交叉。如果你就是想得到从始至终的基态势能面，每个点都stable=opt是最稳妥的做法。

作者
Author: fhh2626 时间: 2015-11-15 22:20

sobereva 发表于 2015-11-14 18:43
这个scan可能会涉及势能面交叉。如果你就是想得到从始至终的基态势能面，每个点都stable=opt是最稳妥的做 ...

感谢~！

不过请问一下，同样的问题在orca中有办法解决吗

作者
Author: sobereva 时间: 2015-11-16 02:59

fhh2626 发表于 2015-11-15 22:20
感谢~！

不过请问一下，同样的问题在orca中有办法解决吗

貌似没有。或者用高斯找到不同距离下基态的轨道占据方式，在orca里也用同样的

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