24小时热门版块排行榜    

查看: 3176  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 274求调剂 +3 时间点 2026-03-13 3/150 2026-03-14 22:27 by 钓鱼佬as
[考研] 一志愿中国科学技术大学,086000,生物与医药 +3 12341234789 2026-03-08 3/150 2026-03-14 04:08 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分,求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
[基金申请] 面上项目还需要AI说明吗? +3 liyundong 2026-03-08 3/150 2026-03-09 22:30 by kingkocxr
信息提示
请填处理意见