简单量子化学问题答疑专帖（旧版，已关帖）

ABetaCarw · 发表于 Post on 2021-1-16 18:36:52

wzkchem5 发表于 2021-1-16 18:13
这是一个很复杂的问题，光照后可能先到一个高的激发态，再弛豫到低激发态，然后在低激发态上旋转。所以不 ...

请问有什么类似课题的研究案例吗？我目前看到的多数文章都是在瞎写甚至有些人连TD关键词默认的是S1都不知道(对普通有机分子)

算一算 · 发表于 Post on 2021-1-16 19:40:47

wzkchem5 发表于 2021-1-16 18:11
试试M06HF？HF交换成分100%，但是有相关项，原则上比HF更准

所有的泛函都不行，都是小数。算了咱们不讨论DFT了，跑题了而且浪费时间，还是回到之前我问的问题上吧，就是CIS和TD-HF在计算force的时候哪个会更准确？你之前说激发能高的时候TD-HF比CIS更准，但是多高算高？

wzkchem5 · 发表于 Post on 2021-1-16 20:44:00

算一算发表于 2021-1-16 19:40
所有的泛函都不行，都是小数。算了咱们不讨论DFT了，跑题了而且浪费时间，还是回到之前我问的问题上吧， ...

我估计3eV算高吧，这个没有一定的标准，因为说实话两者都很差
但是我的观点是，不要为了把一个性质做对，牺牲了其他的性质。固然精确的结果是每个片段的电荷是整数，但这不意味着把电荷做成整数肯定会让结果更好。你必须要证明，忽略correlation的误差对你研究的具体问题的影响比电荷非整数的影响小，才有道理用CIS或TDHF。
其实如果只算单点，CIS(D)、CC2等可能更好，但是一般程序不支持这些方法的解析梯度。再有，如果你算的激发态只涉及少数几个轨道，这些轨道做CAS能做得动的话，可以考虑CASPT2。

wzkchem5 · 发表于 Post on 2021-1-16 20:45:11

ABetaCarw 发表于 2021-1-16 18:36
请问有什么类似课题的研究案例吗？我目前看到的多数文章都是在瞎写甚至有些人连TD关键词默认的是S1都不 ...

可以搜一下rhodopsin的模型体系的研究，涉及双键的顺反异构，这个体系研究应该挺多的

strawberry · 发表于 Post on 2021-1-16 21:30:46

本帖最后由 strawberry 于 2021-1-16 21:34 编辑

C:\Users\Administrator.USER-20140818HS\Desktop

snljty · 发表于 Post on 2021-1-16 22:00:51

不要说“Link502”报错。看http://bbs.keinsci.com/thread-4829-1-1.html，直接把输出文件最后几行贴上就行，比图片看着方便还省空间。
就是自洽场不收敛。看卢老师博文《解决SCF不收敛问题的方法》
http://sobereva.com/61
用CAM-B3LYP还用这么大弥散，居然是为了算个单点能是要闹哪样，明显不合适。还有为什么要用HF做初猜？

strawberry · 发表于 Post on 2021-1-16 22:14:20

谢谢老师的指导

sobereva · 发表于 Post on 2021-1-17 03:54:47

常见问题把报错直接在公社论坛首页google框一搜就找到答案，这是老生常谈的问题

解决SCF不收敛问题的方法
http://sobereva.com/61

sobereva · 发表于 Post on 2021-1-17 03:57:42

算一算发表于 2021-1-16 15:56
那个泛函并不能解决我的问题，它给的也是小数，大约 1.6、1.4这样的。所以还是得用CIS或是TD-HF来算。现 ...

“用CIS或是TD-HF来算”，到时候审稿人看了准一脸黑线，不结合更高级别的比如CC2什么的加以检验就想糊弄过去审稿人准没戏。反正要有这么一篇稿子落我手里我肯定comment这点。但凡用一个非主流，甚至是众所周知精度较烂的方法做计算，不给出充足的证据论证合理性是不可能发表在稍微像样期刊上的。更何况HF完全不考虑动态相关，就凭这一点都会被稍微懂行的审稿人严厉抨击。

在我的认识范畴里不存在只适合用HF算的情况。不要一味地拿是否是整数当做判据，就算真怀疑DFT，也起码拿更高级的波函数方法如耦合簇也算一下电荷分布做对照。

就算体系较大，是CC2、EOM-CCSD优化做不动的情况，最最起码也该拿CIS/TDHF优化后的极小点也用CC2/EOM-CCSD算一下看受力是否已经不大。

lastzealot · 发表于 Post on 2021-1-17 22:20:03

Sob老师，请教您个问题一个Oh结构的分子
在adf里固定对称性优化得到键长为2.04，在molpro里最高到Dnh对称性，于是固定对称性为D2h、在casscf/caspt2(6,10)水平下优化，结果得到的键长不一致(1.98-2.00范围)，有细微的差异。请问在高精度水平下如何保持六个键长一致？
谢谢老师

snljty · 发表于 Post on 2021-1-18 19:16:44

请问各位，在重复卢老师博文里面Gaussian的GVPT2非谐振计算的时候，最终报错Value not found by IFind。但是输出文件里面有频率信息，请问这是什么原因导致的，是否因为无法判断振动模式的不可约表示？对结果会有影响么？谢谢！
输入文件如下
# Freq=Anharmonic wB97XD/TZVP Geom=Connectivity Integral=UltraFineGrid Guess=INDO SCF=NoVarAcc

Optimized under wB97X-D/def-TZVP

0 1
C                0.61058200 3.64191600 0.00000000
C                1.87324500 3.18234400 0.00000000
C                2.80871200 2.39739500 0.00000000
C                3.48056100 1.23371800 0.00000000
C                3.69261400 0.03110600 0.00000000
C                3.45928300 -1.29217800 0.00000000
C                2.84870100 -2.34973800 0.00000000
C                1.81936900 -3.21345000 0.00000000
C                0.67184900 -3.63111300 0.00000000
C                -0.67184900 -3.63111300 0.00000000
C                -1.81936900 -3.21345000 0.00000000
C                -2.84870100 -2.34973800 0.00000000
C                -3.45928300 -1.29217800 0.00000000
C                -3.69261400 0.03110600 0.00000000
C                -3.48056100 1.23371800 0.00000000
C                -2.80871200 2.39739500 0.00000000
C                -1.87324500 3.18234400 0.00000000
C                -0.61058200 3.64191600 0.00000000

输出文件最后几行如下：
（前面还有47个）
WARNING: Irreducible representation not found for mode    48.
WARNING: Irreducible representation not found for rot/transl. num. 1.
WARNING: Irreducible representation not found for rot/transl. num. 2.
WARNING: Irreducible representation not found for rot/transl. num. 3.
WARNING: Irreducible representation not found for rot/transl. num. 4.
WARNING: Irreducible representation not found for rot/transl. num. 5.
WARNING: Irreducible representation not found for rot/transl. num. 6.

Definition of the model system: Active modes
--------------------------------------------
The 48 Active Modes are:
   1    2    3    4    5    6    7    8    9 10 11 12 13
14 15 16 17 18 19 20 21 22 23 24 25 26
27 28 29 30 31 32 33 34 35 36 37 38 39
40 41 42 43 44 45 46 47 48

Data for Anharmonic Potential Energy Surface
--------------------------------------------
WARNING: Unreliable CUBIC force constant i= 21,j= 18,k= 24
- Fjik =    -0.99573289 while Fijk is NULL
WARNING: Unreliable CUBIC force constant i= 18,j= 21,k= 24
- Fjik =    -0.99560968 while MEAN(Fijk,Fjik) is NULL
WARNING: Unreliable CUBIC force constant i=    8,j=    9,k= 13
- Fjik is NULL while MEAN(Fijk,Fjik) =    0.14951167
WARNING: Unreliable CUBIC force constant i=    7,j=    7,k= 10
- Fjik is NULL while Fijk =    0.12873608
WARNING: Unreliable CUBIC force constant i=    7,j=    3,k=    9
- Fjik is NULL while Fijk =    -0.11419075
WARNING: Unreliable CUBIC force constant i=    3,j=    7,k=    9
- Fjik =    -0.11463660 while MEAN(Fijk,Fjik) is NULL
WARNING: Unreliable CUBIC force constant i=    1,j=    1,k= 15
- Fjik is NULL while Fijk =    -0.15406589

Data for Electric Dipole
------------------------
Property available.

   ==================================================
      Vibro-Rotational Analysis Based on Symmetry
   ==================================================

Input/Output information
------------------------
Normal modes will be PRINTED in DESCENDING order (imag. freq. first)
and sorted by irreducible representation
The connection between this new numbering (A) and the one used before
(H) is reported in the present equivalency table:
Value not found by IFind.
Error termination via Lnk1e in /home/zz/g16/l717.exe at Mon Jan 18 18:40:20 2021.

C18_freq_anharm.out.zip (181.33 KB, 下载次数 Times of downloads: 1)

谢谢各位！

wzkchem5 · 发表于 Post on 2021-1-19 11:07:58

lastzealot 发表于 2021-1-17 22:20
Sob老师，请教您个问题一个Oh结构的分子
在adf里固定对称性优化得到键长为2.04，在molpro里最高到Dnh对称 ...

你确定活性空间选对了吗？可能活性空间破坏了对称性

lastzealot · 发表于 Post on 2021-1-19 21:31:00

wzkchem5 发表于 2021-1-19 11:07
你确定活性空间选对了吗？可能活性空间破坏了对称性

活性空间为3个双占据p轨道和7个f轨道，也就是10个前线轨道，6个电子，没有调整轨道顺序。三个占据轨道能量是简并的，且能量比下面的明显高。

wzkchem5 · 发表于 Post on 2021-1-20 10:25:47

lastzealot 发表于 2021-1-19 21:31
活性空间为3个双占据p轨道和7个f轨道，也就是10个前线轨道，6个电子，没有调整轨道顺序。三个占据轨道能 ...

f轨道是简并的吗？会不会某个f轨道在CASSCF过程中被转出去了

329678001 · 发表于 Post on 2021-1-20 10:58:24

sob老师您好，我用B2PLYPD3/def2TZVP算H质子的单点能，任务正常结束，但是没有mp2的能量呢，这个是怎么回事呀，是不是得换计算级别了呀？

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

[综合交流] 简单量子化学问题答疑专帖（旧版，已关帖）

高斯求单点能出错L502