详谈Multiwfn产生ORCA量子化学程序的输入文件的功能

sobereva · 发表于 Post on 2021-5-24 19:33:00

为了更灵活，今日更新的Multiwfn新增加了一个选项，在文中增加了相应的介绍

此界面里有个选项-100，选了之后用户要提供一个ORCA模板文件，其中包含实际计算要用的所有设置，只不过坐标部分写成[geometry]，Multiwfn会将这部分替换为当前体系的坐标来产生ORCA输入文件。通过这个功能可以完全自定义地创建ORCA输入文件，还可以通过自写简单的脚本实现诸如把一大批xyz文件转化为带有特殊关键词的ORCA输入文件。一个模板文件的例子如下，用来产生DSD-BLYP双杂化泛函级别的自然轨道（算完后产生的以.mp2nat为后缀的文件可自行改成.gbw后缀，然后就可以用orca_2mkl转成.molden输入文件，之后可以载入Multiwfn做波函数分析）。
! DSD-BLYP def2-TZVP tightscf miniprint
%maxcore 6000
%pal nprocs 10 end
%mp2
density relaxed
NatOrbs true
end
* xyz 0 1
[geometry]
*

喵星大佬 · 发表于 Post on 2021-5-29 10:47:27

激发态计算可以加入一个CIS(D)的选项，对于很多吸收700nm以上的染料精度大大好于TDDFT和双杂化，比如

! RHF def2-TZVP Def2-TZVP/C noautostart miniprint nopop
%scf
Tole 1e-10
TolRMSP 1e-10
TolMaxP 1e-8
end
%cis
dcorr 2
Nroots 5
MaxDim 20
end

sobereva · 发表于 Post on 2021-5-30 06:14:46

喵星大佬发表于 2021-5-29 10:47
激发态计算可以加入一个CIS(D)的选项，对于很多吸收700nm以上的染料精度大大好于TDDFT和双杂化，比如

双杂化试过DSD-PBEP86-D3(BJ)和wB2GP-PLYP没有？

喵星大佬 · 发表于 Post on 2021-5-30 15:03:21

本帖最后由喵星大佬于 2021-5-30 15:09 编辑

sobereva 发表于 2021-5-30 06:14
双杂化试过DSD-PBEP86-D3(BJ)和wB2GP-PLYP没有？

试了，基本对于750nm以上吸收的染料，激发能都高估了接近100nm了，和TDDFT普通泛函没有太他区别，CIS(D)基本在30nm以内，接近100个原子的体系，def2-SVP def2-SVP/C很容易算，def2-TZVP def2-TZVP/C就很容易挂掉

wzkchem5 · 发表于 Post on 2021-5-30 17:17:59

喵星大佬发表于 2021-5-30 08:03
试了，基本对于750nm以上吸收的染料，激发能都高估了接近100nm了，和TDDFT普通泛函没有太他区别，CIS(D) ...

会不会不是双杂化泛函的问题，而是你没有考虑vibronic effect的关系？

喵星大佬 · 发表于 Post on 2021-5-30 17:58:54

wzkchem5 发表于 2021-5-30 17:17
会不会不是双杂化泛函的问题，而是你没有考虑vibronic effect的关系？

刚性比较强的大共轭染料分子，激发态和基态结构几乎没啥区别(实验光谱的位移也很小)，这个问题应该不是主要问题吧

sobereva · 发表于 Post on 2021-5-31 05:05:03

喵星大佬发表于 2021-5-30 17:58
刚性比较强的大共轭染料分子，激发态和基态结构几乎没啥区别(实验光谱的位移也很小)，这个问题应该不是主 ...

建议用DLPNO-STEOM-CCSD验证一下到底是理论方法的问题还是error cancellation

我不觉得CIS(D)表现更优秀是原理上说得通而且是普遍的现象

sobereva · 发表于 Post on 2021-7-2 08:51:38

今日更新了Multiwfn程序以使得产生ORCA输入文件的功能完全兼容ORCA 5.0，并同时更新了本文。新版界面里专门加入了从ORCA 5.0开始可以用的wB97X-2-D3(BJ)和r2SCAN-3c的选项，并在本文中予以了说明

Y30181009 · 发表于 Post on 2021-7-10 18:56:32

sob老师，请问wB2GP-PLYP在ORCA中计算TD-DFT该如何加DFT-D3，直接写D3最后报错显示无色散校正参数：
Warning: (DFTVDWD): Non-parameterized functional used for dispersion correction !

Warning: (D3BJ): Non-parameterized functional used for dispersion correction !

Please define all necessary parameters (s6, s8, a1, and a2) in the %method block!
[file orca_tools/Tool-Misc/qcdftvdw.cpp, line 36127]: ERROR in DFT dispersion correction!

sobereva · 发表于 Post on 2021-7-11 12:51:04

Y30181009 发表于 2021-7-10 18:56
sob老师，请问wB2GP-PLYP在ORCA中计算TD-DFT该如何加DFT-D3，直接写D3最后报错显示无色散校正参数：
War ...

没必要加D3，色散校正不影响激发能

wuzhiyi · 发表于 Post on 2021-7-23 03:49:21

卢老师好，感谢卢老师每次都非常快速的更新软件。
我还没有测试过，但看说明好像在ORCA>5里，溶剂模型默认是使用gaussian charge scheme了，这貌似导致了SMD模型失去了解析梯度，这样的话，将SMD用在opt+freq的任务就不是很合适了。
虽然可以手动在这种情况下替换成CPCM(Water)，但感觉在非单点能的情况下，当选择 -1 Toggle employing implicit solvation model, current: No的时候，默认使用CPCM是不是对新手更加友好。

sobereva · 发表于 Post on 2021-7-23 04:10:22

wuzhiyi 发表于 2021-7-23 03:49
卢老师好，感谢卢老师每次都非常快速的更新软件。
我还没有测试过，但看说明好像在ORCA>5里，溶剂模型默认 ...

ORCA 5.0 SMD有解析梯度，只不过没有解析Hessian
在-1 Toggle employing implicit solvation model里输入c之后就可以输入溶剂名用CPCM，我觉得不算麻烦
另外极有可能下一个版本就有SMD的解析Hessian了，所以这块就不动它了

wzkchem5 · 发表于 Post on 2021-7-23 04:45:13

sobereva 发表于 2021-7-22 21:10
ORCA 5.0 SMD有解析梯度，只不过没有解析Hessian
在-1 Toggle employing implicit solvation model里输 ...

下一个版本（5.0.1）不会，但是5.1.0有一定概率会有。
SMD没有解析Hessian的原因是非静电项的Hessian没做，而不是gaussian charge scheme的原因。现在我们内部已经基本做完了，可能还需要一些调试

喵星大佬 · 发表于 Post on 2021-7-23 07:48:03

wzkchem5 发表于 2021-7-23 04:45
下一个版本（5.0.1）不会，但是5.1.0有一定概率会有。
SMD没有解析Hessian的原因是非静电项的Hessian没 ...

不过说实话基本不会在SMD下优化+振动分析

qlyfe · 发表于 Post on 2021-8-31 15:09:54

老师，用最新版Multiwfn产生的sTD-DFT任务中除了%maxcore一项外，在end前面还多了一个maxcore，请问这里的maxcore作用是什么呢？

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

[ORCA] 详谈Multiwfn产生ORCA量子化学程序的输入文件的功能