24小时热门版块排行榜    

查看: 1051  |  回复: 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}
回复此楼
心若在,梦就在!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1234zou

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1466036890 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 入职高校3年发表10+SCI,尽人事听天命 +34 kaoyan250 2024-05-27 52/2600 2024-06-03 13:49 by kaoyan250
[论文投稿] 投稿的时候编辑发来邮件说作者机构信息没填写怎么办 5+3 飞飞0421 2024-06-02 7/350 2024-06-03 13:39 by topedit
[基金申请] 为什么美国要封锁华为而不封锁丘C同呢? +13 zju2000 2024-06-02 16/800 2024-06-03 11:54 by chemhua
[基金申请] B口人才项目 +11 WOWO159357 2024-05-29 23/1150 2024-06-03 11:18 by XUAAAAA
[基金申请] 化学B02口青基 代表作都是什么水平的?向大佬求助 +14 arthas_007 2024-06-01 19/950 2024-06-03 11:14 by prisonrabbit
[基金申请] 离职后国自然项目不能变更单位,在新单位还能申请新的国自然项目吗 +3 宜兰 2024-05-31 11/550 2024-06-03 10:26 by lilyxjw
[硕博家园] 材料博士毕业 +3 112233ssg 2024-06-02 3/150 2024-06-03 09:59 by Lwyy1106
[电化学] 电化学传感器 +3 哈哈嘿哟 2024-06-01 5/250 2024-06-03 08:46 by mecart2021
[基金申请] 化学口B0109(高分子合成),拿青年基金一般需要怎样的文章水平? +14 salmon95 2024-05-30 24/1200 2024-06-03 08:31 by LNP@mRNA
[硕博家园] 论大家对6070后普通教授导师的看法 +8 SNaiL1995 2024-05-28 14/700 2024-06-03 05:07 by SNaiL1995
[论文投稿] 求助大神,Fe和Al离子对MOF都有淬灭,当两种离子共存时,怎么区分两种离子? 10+3 maoxiao 2024-06-02 3/150 2024-06-03 01:43 by nono2009
[硕博家园] 各位同学能否分享一下实验室的学生劳务发放标准呀? +19 ma3252788 2024-05-30 19/950 2024-06-02 18:15 by Lcy6666
[基金申请] 九部门发文:不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件 +17 sjtu2012 2024-05-28 21/1050 2024-06-02 13:43 by 欢乐颂叶蓁
[考博] 求25博导,金属增材制造方向 +3 22机械 2024-06-01 3/150 2024-06-02 11:17 by Napoleonsky
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天,完结那天晚上几乎失眠 +11 瞬息宇宙 2024-05-27 13/650 2024-06-01 22:15 by otani
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[论文投稿] 纠结选哪一个期刊,电化学领域 50+8 Freya163 2024-05-28 10/500 2024-05-31 15:09 by wzykobe
[硕博家园] 求问endnote上Wiley引用格式在哪下载 +3 Elon. 2024-05-30 3/150 2024-05-31 10:13 by crazy peng
[材料综合] 真空封石英管 北京 +4 dessha 2024-05-29 5/250 2024-05-30 16:40 by mpdfwxgui
[论文投稿] 审稿专家比较坚定的让补充实验,但实在没法补充实验,修回还有希望吗? (EPI+1) 3+3 qweasd12345 2024-05-29 6/300 2024-05-30 08:11 by qweasd12345
信息提示
请填处理意见