计算化学公社

标题: ORCA混合基组输入求助 [打印本页]

作者
Author: a617788816 时间: 2020-3-16 20:08
标题: ORCA混合基组输入求助
本帖最后由 a617788816 于 2020-3-16 22:51 编辑

ORCA手册混合基组研究了一下，不是看的很懂，论坛里关于ORCA混合基组的帖子也不多，格式是否规范很难看懂，应在molculs的Template文件怎么写，好像很麻烦.. 目标是想对22-32原子部分使用ma-弥散算单电，其他一切不变。
1.各位老师看看这样是不是就是22-32号原子带了弥散，其他没带。并且自动生成了辅助基组并加速计算。一直报错Number of NGauss expected in AddNewGTOs，实在看不懂格式，也没找到原因

2.ECP在手册说可能有必要加，我不清楚什么时候才应当要加？意义是什么？

3.另外，我发现使用Multiwfn生成ORCA输入文件时，使用PWPB95泛函，基组使用的是BJ阻尼，但输出文件仅仅写D3. 这是一个意思的嘛？还是说需要自己补充为D3BJ、在论坛sob老师文章中找到除尼苏达系列泛函,其他都应当用BJ阻尼，可是我用PM6预优化时候，只能用D3，而没办法用D3BJ，Gaussian报错没有参数，这是什么问题吗？

!PWPB95 D3 def2-TZVPP def2/J def2-TZVPP/C RIJCOSX grid4 gridx4 tightSCF miniprint nopop noautostart
%basis
NewAuxJGTO 22-32
      “AutoAux”
DelECP 22-32
end

%maxcore 1500

%pal nprocs 6 end

* xyz 0 1

[GEOMETRY]
      1 . ........
      ..
      22 .. .. ..       NEWGTO "ma-def2-TZVPP" end
      23 .........       NEWGTO "ma-def2-TZVPP" end
      ..
      32  ... .... ..... NEWGTO "ma-def2-TZVPP" end
*

作者
Author: sobereva 时间: 2020-3-17 12:13
DFT-D色散校正的使用
http://sobereva.com/210
明确说了ORCA里D3就是D3BJ

“除尼苏达系列泛函”，我限定的是泛函，PM6是另一码事

作者
Author: a617788816 时间: 2020-3-18 20:29
本帖最后由 a617788816 于 2020-3-19 18:51 编辑

sobereva 发表于 2020-3-17 12:13
DFT-D色散校正的使用
http://sobereva.com/210
明确说了ORCA里D3就是D3BJ

1.关于辅助基组，对于我的体系中部分使用弥散，1-21号不使用弥散（def2-TZVPP）。22-32号使用弥散（ma-def2-TZVPP），是否应当用autoaux生成对应的辅助基组？？？还是我使用同一辅助基组也可以吗？

2. 1-21号原子我的A部分（阳离子不带弥散算单电），22-32原子是我的B部分（阴离子带上弥散算单电），是不是说带上弥散之后，就没有相应的辅助基组，所以我在单独B部分单点时也应当使用自动生成辅助基组吗？目的是使用ORCA为了计算AB结合能（A不需弥散，B要弥散）

3.是不是因为使用了NewGTO改变了基础集，所有必须要加上DelESP 删除相应的ESP？

作者
Author: 柒月小鱼 时间: 2020-3-18 21:07
!B3LYP Def2-SVP
%bais NewECP Au "SDD" end
end

作者
Author: a617788816 时间: 2020-3-19 20:28
各位老师，def2-TZVPP，算单电能，加-ma前后，和辅助基组有没有关系啊，使用同一辅助基组有没有问题的啊？def2/J def2-TZVPP/C RIJCOSX

作者
Author: sobereva 时间: 2020-3-19 21:13

a617788816 发表于 2020-3-19 20:28
各位老师，def2-TZVPP，算单电能，加-ma前后，和辅助基组有没有关系啊，使用同一辅助基组有没有问题的啊？d ...

加了ma-显然不能用原本的def2/J辅助基组。此时应当用autoaux自动产生辅助基组。Multiwfn产生ORCA输入文件时给def2系列加弥散时也是这么处理的

作者
Author: a617788816 时间: 2020-3-20 14:26
本帖最后由 a617788816 于 2020-3-22 06:08 编辑

sobereva 发表于 2020-3-19 21:13
加了ma-显然不能用原本的def2/J辅助基组。此时应当用autoaux自动产生辅助基组。Multiwfn产生ORCA输入文件 ...

谢谢sob老师，由于我是ORCA混合基组输入算单电，所以A部分不加弥散，B部分加弥散-ma，显然B部分是应当使用autoaux的。

1.当我算AB符合体系时，A部分是不是不需要用autoaux，即两套基本基组，两套辅助基组呢？还是说两套基本基组应用同一套自动生成的辅助基组？

2. 如果是使用同一自动生成辅助基组，我单独计算计算A部分但电能时是不是也应当用autoaux呢？

麻烦各位老师解答一下，不胜感激

作者
Author: a617788816 时间: 2020-3-22 06:14
麻烦各位老师帮我看看以上问题，两三天没结果了，谢谢各位老师

作者
Author: zjxitcc 时间: 2020-3-22 08:56
本帖最后由 zjxitcc 于 2020-3-22 11:22 编辑

用RI加辅助基的目的是为了保证精度的前提下，加速计算。因此这个问题你可以自己回答的：用一个稍微小一点的类似体系，或稍微小一点的类似基组，加/不加RI（加RI的话用何种辅助基），算个单点看一下，如果电子能量只有小于0.001 a.u.的变化，那就是完全可以的。你只需相信自己的数据，这个结论不依赖于他人给的建议。
两天足够测试出这个结论。

另外，我觉得grid4 gridx4还是low了点。。。既然都TightSCF了，建议grid6 gridx6

作者
Author: sobereva 时间: 2020-3-22 15:27

a617788816 发表于 2020-3-20 14:26
谢谢sob老师，由于我是ORCA混合基组输入算单电，所以A部分不加弥散，B部分加弥散-ma，显然B部分是应当使 ...

所有都用autoaux

作者
Author: a617788816 时间: 2020-4-2 12:44
本帖最后由 a617788816 于 2020-4-3 13:37 编辑

The energy has been calculated at the MP2 level but the densities used in the property calculations are still SCF densities MP2 response densities have not been calculated use %mp2 Density relaxed end or %mp2 Density unrelaxed end

1. 请问我用PWPB95 QZVPP算单电能，没报错但出来这个警告，是什么问题? 目的是计算复合物结合能
2.另外对比了了一下，B3LYP-D3/def2-TZVP同级别同结构下，分别优化，和频率计算，发现电子能量有约20KJ/mol的差异，正常吗？还是只要保证没虚频就正常??

作者
Author: snljty 时间: 2020-4-2 12:53

a617788816 发表于 2020-4-2 12:44
The energy has been calculated at the MP2 level but the densities used in the property calculations ...

因为你没算MP2密度，这个不影响单点能结果，只有你想获得更好的密度和获得比如后HF(多组态)的自然轨道的时候才需要用。计算这个要额外花一定时间。
没看明白你说的差了20 kJ/mol是谁和谁比，不过你要是算复合物，用B3LYP的话D3色散校正几乎必不可少。你加D3关键词了嘛？

作者
Author: a617788816 时间: 2020-4-2 13:09
本帖最后由 a617788816 于 2020-4-2 13:12 编辑

snljty 发表于 2020-4-2 12:53
因为你没算MP2密度，这个不影响单点能结果，只有你想获得更好的密度和获得比如后HF(多组态)的自然轨道的 ...

使用了D3，已补充。就是同一个复合物，同方法同级别，先算了优化，使用优化后的结构再计算频率，无虚频，但发现电子能量有20kJ/mol差异，担心算错，不知道是否属于合理范围
虽然这个能量似乎不重要，因为后面还进行了高精度的单电能计算。上面补充了一下第三个问题，麻烦您也帮忙看一下

作者
Author: a617788816 时间: 2020-4-3 13:37
本帖最后由 a617788816 于 2020-4-3 13:38 编辑

请问应用于Multiwfn波函数分析，输入文件，应该用低精度（B3LYP-D3/ma-def2-TZVP）下的优化的波函数文件，还是应该是用高精度下(PWPB95/ma-def2-QZVPP)的单电能的波函数文件

作者
Author: sobereva 时间: 2020-4-4 05:25

a617788816 发表于 2020-4-3 13:37
请问应用于Multiwfn波函数分析，输入文件，应该用低精度（B3LYP-D3/ma-def2-TZVP）下的优化的波函数文件， ...

用前者
后者的话实际上给出的波函数只是其中的杂化泛函部分，真需要在双杂化下做波函数分析还得产生自然轨道，麻烦。
而且TZ级别基组做波函数分析就已经足够了，比用QZ时间省很多。

作者
Author: sobereva 时间: 2020-4-4 05:27

a617788816 发表于 2020-4-2 13:09
使用了D3，已补充。就是同一个复合物，同方法同级别，先算了优化，使用优化后的结构再计算频率，无虚频， ...

差20 kJ/mol，绝对是哪里的设置或者结构和之前不同，或者操作犯了什么错误

作者
Author: 13104468322 时间: 2020-6-11 15:56
! PBE0 DKH2 DKH-def2-TZVP SARC/J RIJCOSX cpcm(CH2Cl2) tightSCF miniprint grid4 gridx4
%basis
NewGTO Cu "SARC-DKH-TZVP" end
end
%pal nprocs 24 end
%maxcore 2500
%tddft nroots=10 TDA false dosoc true printlevel 3 end
************************************************

我对Cu的基组进行了设置，输出文件为什还会出现如下报错，望各位大佬指教

**********************************************

There are no main basis functions on atom number 0 (Cu)

作者
Author: abin 时间: 2020-6-11 16:29
本帖最后由 abin 于 2020-6-11 16:34 编辑

13104468322 发表于 2020-6-11 15:56
! PBE0 DKH2 DKH-def2-TZVP SARC/J RIJCOSX cpcm(CH2Cl2) tightSCF miniprint grid4 gridx4
%basis
NewGT ...

你有读ORCA手册嘛？

SARC-DKH-TZVP, Segmented all-electron relativistically contracted basis sets for use with the DKH2 and ZORA Hamiltonians. Available for elements beyond Xe.

在手册中，
同样提到，
For use in DKH or ZORA calculations we provide adapted versions of the def2 basis setsfor the elements H–Kr.
.....
DKH-def2-TZVP and ZORA-def2-TZVP

手册上肯定有这些内容。

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