24小时热门版块排行榜    

查看: 3173  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 332分材料工程调剂 +3 莓好时光海苔 2026-03-09 3/150 2026-03-14 02:03 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 285求调剂 +4 ytter 2026-03-12 4/200 2026-03-13 14:48 by jxchenghu
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 083000环境科学与工程调剂 +8 mingmingry 2026-03-09 9/450 2026-03-11 10:23 by 沙漠之狐994
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见