如何在orca中使用M052X关键词在隐式溶剂模型计算溶解自由能

牧生

社长的博文中 http://sobereva.com/327 ，使用了gaussian隐式溶剂模型计算溶解自由能，其中说了“M05-2X/6-31G*这个很low的级别，除了结合SMD算溶解自由能很适合外，没有任何其它实际用处！”。

我想用orca在相同级别下计算，查了手册，关键词是这样的




比如一个简单的甲烷，隐式溶剂用水，我的输入文件是这样写的

! M052X 6-31G* def2/J RIJCOSX noautostart miniprint
%maxcore     1000
%pal nprocs   4 end
%cpcm
smd true
SMDsolvent "water"
end
* xyz   0   1
C      0.00000000    0.00000000    0.00000000
H      0.62911800    0.62911800    0.62911800
H     -0.62911800   -0.62911800    0.62911800
H     -0.62911800    0.62911800   -0.62911800
H      0.62911800   -0.62911800   -0.62911800
*


结果报错，发现orca不认M052X，同样的，也不认M05

            !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                                 INPUT ERROR
            UNRECOGNIZED OR DUPLICATED KEYWORD(S) IN SIMPLE INPUT LINE
              M052X
            !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[file orca_main/main_input_keywordline.cpp, line 13872]:

[file orca_main/main_input_keywordline.cpp, line 13872]:




但是如果改成M062X，就可以正常算完。

zjxitcc

For RIJCOSX, try

1. ! RKS 6-31G(d) VeryTightSCF RIJCOSX def2/J defgrid3
   
2. %method
   
3. method dft
   
4. exchange hyb_mgga_x_m05_2x
   
5. correlation mgga_c_m05_2x
   
6. end

复制代码

For NoRI, try

1. ! RKS 6-31G(d) VeryTightSCF defgrid3 noRI
   
2. %method
   
3. method dft
   
4. exchange hyb_mgga_x_m05_2x
   
5. correlation mgga_c_m05_2x
   
6. end

复制代码

牧生

zjxitcc 发表于 2025-7-17 15:41
For RIJCOSX, try
For NoRI, try

谢谢，你给的这些关键词是可行的。。社长的博文中也说了，M062X的结果也是可以的，于是我就用的M062X。

但是很想了解一下，为什么M062X就可以，但手册中写出的关键词M052X就不行呢？

Uus/pMeC6H4-/キ

牧生 发表于 2025-7-18 15:41
谢谢，你给的这些关键词是可行的。。社长的博文中也说了，M062X的结果也是可以的，于是我就用的M062X。
...

1楼截图是来自3.3.6.3节的吧，这节一开始就说了“部分LibXC泛函可以用!LibXC(Keyword)调用，如!LibXC(TPSS)”，这跟如3.3.4.1节说的内置的泛函当然写法不同。
M05-2X未内置于ORCA 6.1.0，所以试图些在!关键词行直接写M052X就没法识别，得写成LibXC(M052X)，等价于2楼用%method分别指定交换与相关泛函的写法；而M06-2X是内置的，所以直接在!关键词行写M062X就行。

牧生

Uus/pMeC6H4-/キ 发表于 2025-7-18 15:55
1楼截图是来自3.3.6.3节的吧，这节一开始就说了“部分LibXC泛函可以用!LibXC(Keyword)调用，如!LibXC(TPS ...

谢谢。看着M05和M06都在一起，压根就没想到它们有啥区别。