24小时热门版块排行榜    

查看: 3174  |  回复: 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的回帖

yongleli

木虫 (正式写手)

【答案】应助回帖

建议检查一下多参考性大小。有可能这个分子需要多参考才能算准。我之前算过一个体系,要算7个量子态,找能量最低的那个电子态作为初始猜测,然后做MRCI才能算对。
3楼2024-07-10 13:28:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自科面上基金字体 +3 iwuli 2026-03-12 4/200 2026-03-14 22:33 by iwuli
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 310 070300化学求调剂 +5 扑风铃的猫 2026-03-08 6/300 2026-03-14 03:22 by JourneyLucky
[考研] 0703化学求调剂 +5 很老实人 2026-03-09 5/250 2026-03-14 02:57 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
信息提示
请填处理意见