计算化学公社

标题: Multiwfn已经可以产生ORCA中最常用计算级别的输入文件 [打印本页]

作者
Author: sobereva 时间: 2019-5-5 15:43
标题: Multiwfn已经可以产生ORCA中最常用计算级别的输入文件
注：此帖已经无意义了，因为笔者之后又写了详细、深入得多的说明文章《详谈Multiwfn产生ORCA量子化学程序的输入文件的功能》（http://sobereva.com/490）

最新的Multiwfn可以在http://sobereva.com/multiwfn免费下载。

很早以前，Multiwfn的主功能100的子功能2就可以产生各种主流量子化学程序的输入文件，但是对应的级别基本上都是B3LYP/6-31G*。鉴于ORCA用的人很多，但是特别常用的级别其实也就那么些，然而ORCA的关键词没Gaussian容易写，时间长了容易忘，而且牵扯到利用RI，好多初学者不知道怎么用最佳的关键词，于是在今日更新的Multiwfn 3.6(dev)中，当使用主功能100的子功能2的选项12导出ORCA输入文件时，可以选择常见类型的任务，目前的界面是这样（以后肯定还会有改变，务必以实际情况为准）

-1 Toggle employing SMD solvation model, current: No
0 Select task, current: Single point
1 B97-3c
2 RI-BLYP-D3(BJ)/def2-TZVP
3 RI-B3LYP-D3(BJ)/def2-TZVP(-f)
4 RI-B3LYP-D3(BJ)/def2-TZVP
5 RI-wB97M-V/def2-TZVP
6 RI-PWPB95-D3(BJ)/def2-TZVPP
7 RI-PWPB95-D3(BJ)/def2-QZVPP
8 DLPNO-CCSD(T)/cc-pVTZ with tightPNO
9 CCSD(T)/cc-pVTZ
10 CCSD(T)-F12/cc-pVDZ with RI
11 CCSD(T)/CBS (cc-pVTZ->cc-pVQZ extrapolation)
20 TDA-DFT RI-PBE0/def2-SV(P) with riints_disk (much faster than 21)
21 TDDFT RI-PBE0/def2-SV(P)
22 TDDFT RI-B2GP-PLYP/def2-TZVP
23 EOM-CCSD/cc-pVTZ

选择一个级别，导出的ORCA输入文件就有相应任务最恰当的关键词。选择具体级别之前可以先用选项-1和0分别设置溶剂模型和任务类型。上面看到的这些，从1到11都是算基态的，可以算是计算精度逐渐增加。20到23，都是算激发态的，计算精度也逐渐增加。

由于导出的文件里关键词都是固定的，往往不会恰好满足需要，所以之后大家还需要修改里面的内存使用量、核数设定、激发态算的数目、基组等。反正有了这样的常用计算级别的“模板”关键词，自己改起来是很容易的事。

使用这个功能时，可以用任意Multiwfn支持的含有结构信息的文件作为Multiwfn的输入文件，比如fch/molden/wfn/wfx/mol/mo2/pdb/xyz/gjf等。

自动生成的关键词兼容ORCA 4.1.2，对以前的版本不一定兼容，对ORCA 3肯定不兼容。等ORCA以后出新版本了，如果关键词有需要调整的，在之后最新的Multiwfn里也会进行相应修改。

作者
Author: zsu007 时间: 2019-5-5 18:34
很好的办法，希望用orca的人high起来！

作者
Author: Graphite 时间: 2019-5-5 22:18
接地气又实用的功能！
实验室有一些做合成的同学，经常需要算一些不难但是批量大的东西，包括量化计算、简单的可视化分析之类的。
给脚本又没有运行环境，教写关键词三句两句也解释不清楚。
这下子Multiwfn一条龙就完事儿了。

作者
Author: ABetaCarw 时间: 2019-5-6 14:18
这是奔着把我等菜鸟惯坏了的方向去打磨Multiwfn啊。真“量化界的瑞士军刀”

作者
Author: captain 时间: 2019-5-6 15:43
本帖最后由 captain 于 2019-5-6 15:48 编辑

请问大神这两行
PWPB95 D3 def2-TZVPP def2/J def2-TZVPP/C RIJCOSX grid4 gridx4 tightSCF 对应 6 RI-PWPB95-D3(BJ)/def2-TZVPP
PWPB95 D3 def2-QZVPP def2/J def2-QZVPP/C RIJCOSX grid4 gridx4 tightSCF 对应 7 RI-PWPB95-D3(BJ)/def2-QZVPP
前面应该加上 RI- 吧？

作者
Author: sobereva 时间: 2019-5-7 08:16

captain 发表于 2019-5-6 15:43
请问大神这两行
PWPB95 D3 def2-TZVPP def2/J def2-TZVPP/C RIJCOSX grid4 gridx4 tightSCF 对应 6 RI-P ...

指定了RIJCOSX的时候默认就会用RI。凡是当前情况下默认的关键词，我都基本舍了（遵循我的计算化学剃刀原理）

作者
Author: captain 时间: 2019-5-7 09:20

sobereva 发表于 2019-5-7 08:16
指定了RIJCOSX的时候默认就会用RI。凡是当前情况下默认的关键词，我都基本舍了（遵循我的计算化学剃刀原 ...

明白了大神

作者
Author: hzamber 时间: 2019-5-21 20:49
DLPNO-CCSD(T)/cc-pVTZ with tightPNO不需要对参考态RHF使用RIJK吗？

作者
Author: sobereva 时间: 2019-5-22 03:08

hzamber 发表于 2019-5-21 20:49
DLPNO-CCSD(T)/cc-pVTZ with tightPNO不需要对参考态RHF使用RIJK吗？

这属于高精度级别，用RI对精度的降低不可完全忽略（手册里的DLPNO-CCSD(T)例子也都没用RI大抵也是考虑到这点）。而且这个级别一般也都不用于很大体系，不开RI做SCF算得动

作者
Author: blueyangliu 时间: 2019-6-16 15:38

sobereva 发表于 2019-5-7 08:16
指定了RIJCOSX的时候默认就会用RI。凡是当前情况下默认的关键词，我都基本舍了（遵循我的计算化学剃刀原 ...

为什么两个辅助基组都加上呢？def2/J def2-TZVPP/C， RIJCOSX只需要库仑部分辅助基组吧，

作者
Author: sobereva 时间: 2019-6-16 18:56

blueyangliu 发表于 2019-6-16 15:38
为什么两个辅助基组都加上呢？def2/J def2-TZVPP/C， RIJCOSX只需要库仑部分辅助基组吧，

双杂化的MP2部分需要用/C

作者
Author: blueyangliu 时间: 2019-6-16 23:02

sobereva 发表于 2019-6-16 18:56
双杂化的MP2部分需要用/C

谢谢社长。第20选项，输出为! PBE0 def2-SV(P) def2/J def2-SVP/C RIJCOSX。PBE0也不是双杂化泛函啊，怎么也用了def2-SVP/C

作者
Author: sobereva 时间: 2019-6-17 01:12

blueyangliu 发表于 2019-6-16 23:02
谢谢社长。第20选项，输出为! PBE0 def2-SV(P) def2/J def2-SVP/C RIJCOSX。PBE0也不是双杂化泛函啊，怎 ...

riints_disk要求指定/C基组

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