| 查看: 3059 | 回复: 2 | ||
| 【悬赏金币】回答本帖问题，作者1466036890将赠送您 110 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
1466036890铜虫 (正式写手)
|
[求助]
molpro计算中SCF能量错误地收敛到较高的态已有2人参与
|
|
|
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} |
回复此楼» 猜你喜欢
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
物理化学论文润色/翻译怎么收费?
已经有104人回复
-
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
-
书籍求助：汽车市场营销理论与实务（电子版）——章小平
已经有0人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有26人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有45人回复
-
东方理-中科大联合博士生招聘
已经有0人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
yongleli
木虫 (正式写手)
- 应助: 93 (初中生)
- 金币: 3769.1
- 散金: 369
- 红花: 24
- 帖子: 764
- 在线: 101.3小时
- 虫号: 303595
- 注册: 2006-12-03
- 专业: 原子和分子物理
3楼2024-07-10 13:28:28
1234zou
金虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 1941.5
- 红花: 5
- 帖子: 92
- 在线: 31.6小时
- 虫号: 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