| 查看: 1049 | 回复: 1 | |||
| 【悬赏金币】回答本帖问题，作者1466036890将赠送您 110 个金币 | |||
1466036890铜虫 (正式写手)
|
[求助]
molpro计算中SCF能量错误地收敛到较高的态已有1人参与
|
||
|
molpro计算某个过渡态的SCF能量收敛到了其他态，并非基态能量。与orca和gaussian等软件计算的SCF能量对不上，本人试了给个初猜轨道和减小能量收敛限，仍就收敛不到基态能量，这个分子没有对称性。虫友们有什么建议吗？输入如下： ***,ts5 memory,5000,M print,orbitals,basis symmetry,nosym geomtyp=xyz geometry={ 17 ts5 C -0.58190338 -0.29113736 -0.26498053 H -0.65442938 -0.15726736 -1.35048253 C -1.88987438 -0.51838836 0.39308347 H -1.87495838 -0.93184336 1.39649247 C -3.18321738 -0.10897436 -0.20109953 H -3.75259338 -1.00267736 -0.49832153 H -3.80709738 0.42903064 0.51910347 H -3.04383438 0.51122464 -1.08762653 O -0.53130838 0.95524264 0.39509547 O 0.73677062 1.89780064 -0.29047753 H 0.55621762 2.67405064 0.25874547 C 0.51258362 -1.29805736 0.08935247 H 0.62157062 -1.39013436 1.17482747 H 0.23000762 -2.27088536 -0.33078453 O 1.77006262 -1.00920736 -0.49590353 O 2.49056562 -0.09430436 0.36571047 H 2.07483362 0.76433664 0.12919547 } basis=aug-cc-pVDZ {hf,energy=1E-12;wf,65,1,1} |
回复此楼» 猜你喜欢
纳米银粉漂浮在溶剂上部，混合不均匀
已经有3人回复
-
一站式科研采购平台来啦！！电池相关专业快过来~
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有282人回复
-
全电池问题
已经有6人回复
-
避坑- 北京清析技术研究院 - 一个典型的坑人检测机构
已经有16人回复
-
催化材料信息学博士后招聘 总包40-60w/年 松山湖材料实验室刘利峰研究员、张磊教授
已经有5人回复
-
从容量瓶转移溶液
已经有2人回复
-
想写一篇关于气液相变---临界点理论的文章，含气液界面张力。邀请感兴趣的同仁加入！
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
1234zou
金虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 1938.5
- 红花: 4
- 帖子: 92
- 在线: 31.5小时
- 虫号: 2914121
- 注册: 2014-01-04
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★
paramecium86: 金币+3, 应助指数+1, 感谢回复。 2023-08-14 12:30:14
paramecium86: 金币+3, 应助指数+1, 感谢回复。 2023-08-14 12:30:14
|
极其简单。 （1）使用fch2com小程序从Gaussian传轨道给Molpro，使后者1圈收敛。注意，高斯gjf文件里最好写上nosymm int=nobasistransform，即不适用点群对称性，不要让程序旋转分子。 （2）使用mkl2com小程序从ORCA传轨道给Molpro，使后者几乎1圈收敛。注意，ORCA里要写TightSCF或VeryTightSCF，否则收敛过松。 更严谨地，你还可以在Gaussian或ORCA里检验好波函数稳定性，然后再传轨道给Molpro。就连Gaussian->ORCA之间也有fch2mkl小程序。其实你只要算1个程序，其他程序都可以立即收敛，节约时间。 fch2com和mkl2com都是MOKIT里的小程序。MOKIT免费、开源，可联网在线安装，亦可下载Linux预编译版解压写好环境变量就用。 程序仓库主页：https://gitlab.com/jxzou/mokit 中文简介：https://gitlab.com/jxzou/mokit/-/blob/master/README_zh.md fch2com文档：https://jeanwsr.gitlab.io/mokit- ... 5.html#4520-fch2com |
2楼2023-08-12 16:41:53
10