24小时热门版块排行榜

>论坛更新日志 (4022)
>虫友互识 (644)
>文献求助 (491)
>导师招生 (369)
>考博 (207)
>休闲灌水 (202)
>论文投稿 (188)
>硕博家园 (162)
>基金申请 (119)
>招聘信息布告栏 (101)
>博后之家 (83)
>教师之家 (72)
>论文道贺祈福 (61)
>公派出国 (53)
>考研 (48)
>找工作 (43)

返回列表

【悬赏金币】回答本帖问题，作者1466036890将赠送您 110 个金币

1466036890

铜虫 (正式写手)

应助: 0 (幼儿园)
金币: 1591.7
散金: 83
帖子: 353
在线: 34.1小时
虫号: 3830547
注册: 2015-04-24
性别: MM
专业: 物理化学

[求助] 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}

回复此楼

» 猜你喜欢

吉林大学材料物理本科生求问调剂信息已经有23人回复
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！已经有1人回复
物理化学论文润色/翻译怎么收费? 已经有263人回复
光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有0人回复
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes 已经有2人回复
帮导师招2026CSC博士（巴塞罗那自治大学UAB-CSC博士项目）已经有21人回复
有没有LIMS客户的资源，求推荐已经有0人回复
福州大学梁宇航副教授招收2026年申核制博士研究生/硕士研究生（理论计算方向）已经有20人回复
福州大学梁宇航副教授招收2026年申核制博士研究生/硕士研究生（理论计算方向）已经有18人回复

» 本主题相关商家推荐: (我也要在这里推广)

心若在，梦就在！！！

1楼 2023-06-09 23:28:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

1234zou

金虫 (小有名气)

应助: 7 (幼儿园)
金币: 1941.5
红花: 5
帖子: 92
在线: 31.6小时
虫号: 2914121
注册: 2014-01-04
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★
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