计算化学公社

 找回密码 Forget password
 注册 Register

Mg-Cl surface scan

查看数: 16829 | 评论数: 10 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2015-11-13 19:05

正文摘要:

我最近做了一个小测试,以Mg2+和Cl-的距离为反应坐标做了一个surface scan,结果发现不同方法所得的结果相差巨大,似乎无法用误差来解释,而且曲线形状也不一样,不知道这是什么原因呢? 使用开闭壳层似乎没什么影响 ...

回复 Reply

sobereva 发表于 Post on 2015-11-16 02:59:40
fhh2626 发表于 2015-11-15 22:20
感谢~!

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


貌似没有。或者用高斯找到不同距离下基态的轨道占据方式,在orca里也用同样的
fhh2626 发表于 Post on 2015-11-15 22:20:31
sobereva 发表于 2015-11-14 18:43
这个scan可能会涉及势能面交叉。如果你就是想得到从始至终的基态势能面,每个点都stable=opt是最稳妥的做 ...

感谢~!

不过请问一下,同样的问题在orca中有办法解决吗
sobereva 发表于 Post on 2015-11-14 18:43:13
fhh2626 发表于 2015-11-14 18:11
多谢指教,我现在解决了波函数稳定性的问题了,不过加上stable=opt以后似乎就不能做scan了

看来只能写 ...

这个scan可能会涉及势能面交叉。如果你就是想得到从始至终的基态势能面,每个点都stable=opt是最稳妥的做法。
fhh2626 发表于 Post on 2015-11-14 18:11:25
zyj19831206 发表于 2015-11-14 17:28
这个是IRC曲线吧??

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

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

看来只能写个脚本在每个位置都算单点了?
zyj19831206 发表于 Post on 2015-11-14 17:28:39
这个是IRC曲线吧??
sobereva 发表于 Post on 2015-11-14 00:48:05
fhh2626 发表于 2015-11-13 22:47
感谢指教~我做了波函数稳定性测试,并且使用了开壳层和guess=mix,但是结果还是
The wavefunction has a ...

不清楚stable=opt怎么不能结束,但应该是初猜波函数轨道占据方式问题,初猜不是基态排布,而不在于非限制性开壳层的使用上。用guess=alter调换调换初猜时对不同不可约表示轨道的占据方式,或许能收敛到稳定波函数通过stable检测。也注意看看HF计算时不同不可约表示占据轨道的数目和DFT时是否不同。
fhh2626 发表于 Post on 2015-11-13 22:47:20
sobereva 发表于 2015-11-13 21:33
可以用高斯重复一下看看,如果也是这样,对比较远的点检验下波函数稳定性,没准收敛到了不同电子态

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

使用stable=opt程序不能正常结束,请问怎么办呢?
fhh2626 发表于 Post on 2015-11-13 21:46:37
sobereva 发表于 2015-11-13 21:33
可以用高斯重复一下看看,如果也是这样,对比较远的点检验下波函数稳定性,没准收敛到了不同电子态

用高斯重复过,结果差不多,奇怪的是对NaCl的计算结果却没问题
sobereva 发表于 Post on 2015-11-13 21:33:11
可以用高斯重复一下看看,如果也是这样,对比较远的点检验下波函数稳定性,没准收敛到了不同电子态

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 23:26 , Processed in 0.166541 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list