计算化学公社

标题: 如何在orca中使用M052X关键词在隐式溶剂模型计算溶解自由能 [打印本页]

作者
Author: 牧生 时间: 2025-7-17 15:13
标题: 如何在orca中使用M052X关键词在隐式溶剂模型计算溶解自由能
本帖最后由牧生于 2025-7-17 15:16 编辑

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

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

(, 下载次数 Times of downloads: 26)

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

! 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，就可以正常算完。

作者
Author: zjxitcc 时间: 2025-7-17 15:41
本帖最后由 zjxitcc 于 2025-7-17 15:43 编辑

For RIJCOSX, try

! RKS 6-31G(d) VeryTightSCF RIJCOSX def2/J defgrid3
%method
method dft
exchange hyb_mgga_x_m05_2x
correlation mgga_c_m05_2x
end

复制代码

For NoRI, try

! RKS 6-31G(d) VeryTightSCF defgrid3 noRI
%method
method dft
exchange hyb_mgga_x_m05_2x
correlation mgga_c_m05_2x
end

复制代码

作者
Author: 牧生 时间: 2025-7-18 15:41
本帖最后由牧生于 2025-7-18 15:43 编辑

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

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

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

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-7-18 15:55

牧生发表于 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就行。

作者
Author: 牧生 时间: 2025-7-18 19:16

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

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

作者
Author: adamm 时间: 2025-9-1 10:11

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

感谢！看说明还需仔细啊

作者
Author: 加加加 时间: 2025-12-29 18:13

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

大佬，我对比了一下高斯# m052x/6-31g(d) scrf=SMD 和您提供的关键词加上SMD(water)，结果两者数值不一致，请问如何调节LibXC参数能使这两个软件的值一致啊

作者
Author: zjxitcc 时间: 2025-12-29 18:19
本帖最后由 zjxitcc 于 2025-12-29 18:21 编辑

加加加发表于 2025-12-29 18:13
大佬，我对比了一下高斯# m052x/6-31g(d) scrf=SMD 和您提供的关键词加上SMD(water)，结果两者数值不一 ...

无法一致。不同量化软件的DFT格点实现细节不同，隐式溶剂模型的实现细节也不同，无法通过关键词调整。只有在气相下、关闭RI近似、不使用DFT，使用HF等波函数方法，这些严格限制下才能算出一致的结果。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)