PW6B95计算单点报错后加! NoTRAH会不会影响单点能？

Oliviaw · 发表于 Post on 2024-3-21 10:06:16

本帖最后由 Oliviaw 于 2024-3-22 20:20 编辑

我对一个含有过渡金属的体系算单点能，方法是PW6B95-D3/def2-tzvp，出现报错如下：
[file orca_tools/Tool-Numint/qcfxc.cpp, line 14990, Process 4]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel. Please try using LIBXC instead!

论坛里有人有相同的报错，他用的PWPB95，解决办法是加! NoTRAH。

1）我用PW6B95 加! NoTRAH能否解决报错问题，会不会影响单点能的最终结果?
2）算结合能的时候涉及三个结构，如果我对其中一个结构加了 ! NoTRAH，其他两个结构也必须一定要加 ! NoTRAH吗？
3）我还需要对比不同的结构体系，是不是其中一个结构体系单点能用了! NoTRAH，其他结构体系也要统一用 ! NoTRAH算单点？

sobereva · 发表于 Post on 2024-3-30 20:41:22

Oliviaw 发表于 2024-3-30 14:26
论坛进不去是自从我最近设置了自动登录后开始的。之后登录即使不设置自动登录，也无法进入。现在只能用流 ...

原理上跟自动登录不会有直接关系

Oliviaw · 发表于 Post on 2024-3-30 14:26:18

sobereva 发表于 2024-3-28 00:19
论坛访问问题是你的ISP的事，给ISP打电话投诉就完了

论坛进不去是自从我最近设置了自动登录后开始的。之后登录即使不设置自动登录，也无法进入。现在只能用流量，才能进论坛

sobereva · 发表于 Post on 2024-3-28 00:19:44

Oliviaw 发表于 2024-3-27 15:04
社长，请问如果对不含过渡金属的体系用pw6b95算单点，对含过渡金属的体系用wB97M-V，可以吗？
基组都用d ...

毫无理由用pw6b95，wB97M-V算什么体系都比这个强

sobereva · 发表于 Post on 2024-3-28 00:19:08

Oliviaw 发表于 2024-3-27 13:52
第一步交给gaussian算就遇到问题，我其他体系的结构算单点能都用cpcm模型，但是如果给gaussian算就不能用 ...

论坛访问问题是你的ISP的事，给ISP打电话投诉就完了

Oliviaw · 发表于 Post on 2024-3-27 20:42:53

wzkchem5 发表于 2024-3-27 17:27
在高斯里不建议用CPCM的原因是因为PCM比CPCM准，不是因为高斯的CPCM做得不好。
如果你做高斯计算的目的 ...

我现在知道为什么要在gaussian用CPCM了，谢谢

wzkchem5 · 发表于 Post on 2024-3-27 17:27:13

Oliviaw 发表于 2024-3-27 06:52
第一步交给gaussian算就遇到问题，我其他体系的结构算单点能都用cpcm模型，但是如果给gaussian算就不能用 ...

在高斯里不建议用CPCM的原因是因为PCM比CPCM准，不是因为高斯的CPCM做得不好。
如果你做高斯计算的目的是为了给ORCA的CPCM计算提供初猜，那用CPCM就比用PCM好了，因为高斯的CPCM显然比高斯的PCM更接近ORCA的CPCM。此时你的目标不是为了让高斯算准，而是让高斯收敛的波函数尽量接近ORCA收敛的波函数，那么高斯和ORCA计算的理论级别就应该尽量一致。
网上的说法是假设你直接报道高斯计算的结果，没有考虑到用高斯给其他程序准备初猜的情况。所以一定不能教条化理解，要充分理解网上这个说法背后的原因，之后自然就知道这个说法的适用范围

Oliviaw · 发表于 Post on 2024-3-27 15:06:35

zjxitcc 发表于 2024-3-27 13:56
我不对在Gaussian中用/不用cpcm 这个话题做评论。如果你一定要用cpcm，可以先用别的隐式溶剂模型（或者 ...

您说的有道理，我先试下别的隐式模型。

Oliviaw · 发表于 Post on 2024-3-27 15:04:14

sobereva 发表于 2024-3-22 20:45
PW6B95并不适合计算含有过渡金属类的体系， wB97M-V、MN15等都明显更合适
含过渡金属的体系本来就普遍更 ...

社长，请问如果对不含过渡金属的体系用pw6b95算单点，对含过渡金属的体系用wB97M-V，可以吗？
基组都用def2-TZVP.
最后得到的两个方法算的结合能能做对比吗？

zjxitcc · 发表于 Post on 2024-3-27 13:56:08

Oliviaw 发表于 2024-3-27 13:52
第一步交给gaussian算就遇到问题，我其他体系的结构算单点能都用cpcm模型，但是如果给gaussian算就不能用 ...

我不对在Gaussian中用/不用cpcm 这个话题做评论。如果你一定要用cpcm，可以先用别的隐式溶剂模型（或者甚至是气相）把SCF算收敛了，然后在cpcm下读取收敛的轨道作为初猜应该可以很快收敛。

Oliviaw · 发表于 Post on 2024-3-27 13:52:40

zjxitcc 发表于 2024-3-27 13:23
MOKIT提供了大量传轨道的小程序，用户想用哪个就用哪个；你也可以从ORCA往Gaussian传。
Gaussian->ORCA ...

第一步交给gaussian算就遇到问题，我其他体系的结构算单点能都用cpcm模型，但是如果给gaussian算就不能用cpcm，网上都说不建议在gaussian用cpcm。

我打算先看gaussian能不能收敛，再考虑要不要用MOKIT或者multiwfn.

最近不知道怎么回事，科音论坛的网页进不去，总是提示要域名白名单，我在火绒里加了科音论坛的白名单，还是不行。

zjxitcc · 发表于 Post on 2024-3-27 13:23:22

本帖最后由 zjxitcc 于 2024-3-27 13:29 编辑

Oliviaw 发表于 2024-3-27 12:56
其实，我开始用orca是因为gaussian算类似的体系不收敛，或者收敛很慢。换orca后就很快收敛了。

MOKIT提供了大量传轨道的小程序，用户想用哪个就用哪个；你也可以从ORCA往Gaussian传。
Gaussian->ORCA可以用fch2mkl，详细示例在本帖4L中有。
ORCA->Gaussian可以用mkl2fch或molden2fch。
也可以从/向很多其他量化程序（如PySCF, Molpro, OpenMolcas, GAEMSS, PSI4, Dalton, BDF, AMESP）传轨道，看图可了解概况。

我只能说如果你哪一个具体例子有SCF收敛困难，愿意发上来的话，我可以用MOKIT的小程序演示给你看。至于你想怎么做（一个个关键词尝试去解决SCF不收敛、还是ORCA传Gaussian、还是采取其他做法），最终还是取决于你自己。

zjxitcc · 发表于 Post on 2024-3-27 13:19:01

Oliviaw 发表于 2024-3-27 12:54
您好！
“这里不需要写基组名称def2TZVP，因为在inp文件下方已经写了。”这句说inp文件下方写了，是说gb ...

fch2mkl小程序生成的inp和mkl文件里都有def2-TZVP基组数据信息，由mkl生成的gbw文件里也有。而在生成的inp文件里不用写def2-TZVP基组名称，这是因为inp文件里有def2-TZVP基组数据信息（夹在坐标中，打开inp文件一看便知）。如果你非要在!这一行中写def2-TZVP，一般是没有影响的（我没观察到过有影响的情况）。

Oliviaw · 发表于 Post on 2024-3-27 12:56:38

wzkchem5 发表于 2024-3-22 18:19
TRAH算得慢是没错的，但是”结果通常不收敛“有多大比例是TRAH本身的问题，有多大比例是只有对于难收敛的 ...

其实，我开始用orca是因为gaussian算类似的体系不收敛，或者收敛很慢。换orca后就很快收敛了。

Oliviaw · 发表于 Post on 2024-3-27 12:54:39

zjxitcc 发表于 2024-3-22 17:40
加! NoTRAH不一定能解决报错，它只是让ORCA不要用TRAH收敛方法，而是多跑一些常规SCF迭代，以期能侥幸收 ...

您好！
“这里不需要写基组名称def2TZVP，因为在inp文件下方已经写了。”这句说inp文件下方写了，是说gbw文件是用def2-TZVP，然后inp就不用写def2-TZVP了吗？我不懂这句话指的是哪个位置写了基组。
我如果在！后面自己写基组会有影响吗？

sobereva · 发表于 Post on 2024-3-22 20:45:32

Oliviaw 发表于 2024-3-22 17:14
您好社长，我加了NoTRAH和defgrid3算单点，不收敛报错如下：
********************************* ...

PW6B95并不适合计算含有过渡金属类的体系， wB97M-V、MN15等都明显更合适
含过渡金属的体系本来就普遍更难收敛
参考下文的解决思路。针对ORCA也有专门的帮助收敛的做法可以尝试，手册里有提，本论坛里也都讨论过，首页google框搜
解决SCF不收敛问题的方法
http://sobereva.com/61

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] PW6B95计算单点报错后加! NoTRAH会不会影响单点能？