计算化学公社

标题: 使用ORCA软件在外加电场后计算吸附能体系不收敛的求助 [打印本页]

作者
Author: Kaibin 时间: 2024-10-19 16:45
标题: 使用ORCA软件在外加电场后计算吸附能体系不收敛的求助
本帖最后由 Kaibin 于 2024-10-19 16:51 编辑

我的计算情况是，在一个由四个环组成的聚合物上计算一个小分子的吸附。
在没有外加的电场的情况下，可以正常计算出结果，在今天我给它施加了一个外加电场，代码如下：

! B97-3c opt freq tightSCF noautostart miniprint nopop
%maxcore 1000
%pal nprocs 4 end
%scf
efield 0.0,0.0,-0.7
end

但是出现了如下报错：
[file orca_leanscf/orca_leanscf.cpp, line 250, Process 2]: Error (ORCA_LEANSCF): unfortunately, the SCF has not converged. There may be a way out but we have to stop here
我想问的是，由于电场设置导致结构验证拉伸变形，无法收敛，我应该作何改进？

并且我目前还存在一个不知道如何设置电场的问题
我需要设置一个垂直于碳环的电场，但我不确定我的电场是否垂直于碳环。我从体系xyz中分辨不出碳环的平面，不知道有没有什么软件可以方便观看。

作者
Author: zjxitcc 时间: 2024-10-20 09:29
初始结构合理吗？是否在无电场下优化过？
“设置一个垂直于碳环的电场”这看上去没啥意义吧，体系在结构优化过程中会变化，初始的电场方向对优化第1步结构的碳环就不适用了。

作者
Author: wzkchem5 时间: 2024-10-20 13:34
仔细阅读https://www.faccts.de/docs/orca/ ... nite-electric-field后面NOTE的第一条
上来就加0.7 au的电场，十有八九是不知道0.7 au的电场是多么恐怖的大小

作者
Author: sobereva 时间: 2024-10-20 13:43
看看通常研究的电场与化学体系相互作用的电场的数量级是多少：
一篇文章深入揭示外电场对18碳环的超强调控作用
http://sobereva.com/570（http://bbs.keinsci.com/thread-19364-1-1.html）
理论设计由18碳环与锂原子构成的电场可控的光学开关
http://sobereva.com/630（http://bbs.keinsci.com/thread-27232-1-1.html）

并且搞清楚电场的原子单位，换算成常用单位有多大
1 a.u.电场=1 Hartree/(e*a0)=51.421 V/Angstrom=514.21 V/nm=5.1421E11 V/m=5.1421 E6 kV/cm

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