24小时热门版块排行榜    

查看: 3392  |  回复: 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}
回复此楼
心若在,梦就在!!!
已阅   回复此楼   关注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的回帖

yongleli

木虫 (正式写手)

【答案】应助回帖

建议检查一下多参考性大小。有可能这个分子需要多参考才能算准。我之前算过一个体系,要算7个量子态,找能量最低的那个电子态作为初始猜测,然后做MRCI才能算对。
3楼2024-07-10 13:28:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1466036890 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗?20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见