ORCA混合基组输入求助

a617788816 · 发表于 Post on 2020-3-16 20:08:42

本帖最后由 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
*

abin · 发表于 Post on 2020-6-11 16:29:54

本帖最后由 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

手册上肯定有这些内容。

13104468322 · 发表于 Post on 2020-6-11 15:56:13

! 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)

sobereva · 发表于 Post on 2020-4-4 05:27:23

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

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

sobereva · 发表于 Post on 2020-4-4 05:25:48

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

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

a617788816 · 发表于 Post on 2020-4-3 13:37:48

本帖最后由 a617788816 于 2020-4-3 13:38 编辑

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

a617788816 · 发表于 Post on 2020-4-2 13:09:50

本帖最后由 a617788816 于 2020-4-2 13:12 编辑

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

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

snljty · 发表于 Post on 2020-4-2 12:53:56

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关键词了嘛？

a617788816 · 发表于 Post on 2020-4-2 12:44:45

本帖最后由 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的差异，正常吗？还是只要保证没虚频就正常??

sobereva · 发表于 Post on 2020-3-22 15:27:04

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

所有都用autoaux

zjxitcc · 发表于 Post on 2020-3-22 08:56:25

本帖最后由 zjxitcc 于 2020-3-22 11:22 编辑

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

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

a617788816 · 发表于 Post on 2020-3-22 06:14:46

麻烦各位老师帮我看看以上问题，两三天没结果了，谢谢各位老师

a617788816 · 发表于 Post on 2020-3-20 14:26:56

本帖最后由 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呢？

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

sobereva · 发表于 Post on 2020-3-19 21:13:49

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

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

a617788816 · 发表于 Post on 2020-3-19 20:28:39

各位老师，def2-TZVPP，算单电能，加-ma前后，和辅助基组有没有关系啊，使用同一辅助基组有没有问题的啊？def2/J def2-TZVPP/C RIJCOSX

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

[ORCA] ORCA混合基组输入求助

浏览过的版块