简单量子化学问题答疑专帖

sobereva

cxh19990601 发表于 2023-3-24 00:35
老师好，请问一下ADC(2)和RI加速的RI-ADC(2)除了速度不同外，精度等会有所差异吗？

用RI总会引入多多少少的误差，没有差异是不可能的，只不过是大小的问题。二者的差异可以控制在一般研究完全能接受的程度以内。

wzkchem5

cxh19990601 发表于 2023-3-23 17:35
老师好，请问一下ADC(2)和RI加速的RI-ADC(2)除了速度不同外，精度等会有所差异吗？

RI的精度稍低，但引入的误差远低于方法本身的误差，一般低2个数量级左右。
这个结论不仅适用于ADC(2)，也适用于几乎所有其他方法

wzkchem5

498746012 发表于 2023-3-23 16:17
各位老师好，请教一下苯环上连有吸电子基团如卤素或者三氟甲基，它们的诱导效应可以通过共轭体系传递吗？其 ...

另一种方法是查一下这些基团的Hammett取代基常数

for_you

sobereva 发表于 2023-3-24 00:35
基组定义了做量子化学计算时，各种原子都用多少以及什么特征的基函数。所以有了一个给定的化学体系，再告 ...

感谢老师解答

糖糖糖豆9988

请教各位老师，集群计算产生的chk文件，利用高斯转化成fch文件后，multiwfn在读取时会出现图片内容，之后闪退；这是因为集群产生的chk文件在读取方面有问题吗？谢谢老师

xiayu1985

wzkchem5 发表于 2023-3-13 16:57
试过把结构优化收敛限设严了吗？另外也可以试试加大格点

您好，还是之前咨询过得，中间体结构无法收敛至一样的问题。我已经使用如下的设置进行了计算：opt=(calcfc,tight,gdiis) int=superfine freq scrf(SMD,solvent=anisole) pbe1pbe/GenEcp Temperature=323.15。但是还是无法完全一样，达不到能量值小数点后五位一样的标准。具体的差异见第一张图，苯环等部位都有略微差异，无法完全重叠。如果不加溶剂化效应，如第二张图所示，只有叔丁基(红色部分)的位置无法重叠，其他部分都是可以完美重叠的。不知道还有什么办法吗？

糖糖糖豆9988

糖糖糖豆9988 发表于 2023-3-24 09:41
请教各位老师，集群计算产生的chk文件，利用高斯转化成fch文件后，multiwfn在读取时会出现图片内容，之后闪 ...

问题解决啦，在linux下用命令行转换才可以，打扰大家

wzkchem5

xiayu1985 发表于 2023-3-24 07:15
您好，还是之前咨询过得，中间体结构无法收敛至一样的问题。我已经使用如下的设置进行了计算：opt=(calcf ...

那可能本来就是两个相距极近的局部极小值点。另外注意一下有没有虚频

桃子量子

计算出只有一个虚频的过渡态，做IRC，设置maxpoints＝50，stepsize＝5，lqa
跑出来的IRC只有三个点，还正确的结束了，
想请教一下这是什么原因
输入文件如下
%nprocshared=40
%mem=4GB
%chk=C4F9SO3--TS-1-IRC.chk
#p irc=(calcfc,maxpoints=50,stepsize=5,lqa) m062x/6-311+g(2d,p) scrf=(smd,solvent=water) nosymm


C4F9SO3--TS-1-IRC


-2 2
C                 -3.33878900   -0.03226000    0.28230800
C                 -2.08569600    0.78355800   -0.09658800
C                 -0.75080600    0.01003700    0.02854500
C                  0.52189100    0.82155100    0.03881800
S                  2.85190100   -0.24258600    0.04356000
O                  3.77063600   -0.00660500    1.22683100
O                  2.45925400   -1.70329000    0.03592700
O                  3.61382200    0.05145500   -1.23201400
F                 -3.31797200   -1.22113000   -0.30757900
F                 -3.40496500   -0.20456300    1.59471100
F                 -0.70624600   -0.84775500   -1.01483000
F                 -2.23511000    1.19411900   -1.36496700
F                 -0.80822500   -0.70751600    1.17347900
F                  0.51763500    1.66447400   -1.02977100
F                 -4.42841800    0.61697100   -0.11176500
F                 -2.06023500    1.86249900    0.70288100
F                  0.55779000    1.59241100    1.15615600
IRC图如下

wzkchem5

仔细阅读http://sobereva.com/400

桃子量子

好的，谢谢您

sobereva

糖糖糖豆9988 发表于 2023-3-24 09:41
请教各位老师，集群计算产生的chk文件，利用高斯转化成fch文件后，multiwfn在读取时会出现图片内容，之后闪 ...

服务器上chk没转换完整或者转换失败，或者下载的时候没下载全。如果你确认不存在这些问题，把gjf文件和molden文件压缩后传网盘上给出下载链接。

糖糖糖豆9988

请教各位老师，用orca计算三重态的O2，如下的设置正确吗？
%pal nprocs 8
end
! UHF DLPNO-CCSD(T) tightPNO RIJK cc-pVTZ cc-pVTZ/JK cc-pVTZ/C tightSCF
%maxcore 8000
* xyz 0 3
O                  0.00000000    0.00000000    0.59430300
O                  0.00000000    0.00000000   -0.59430300
*

糖糖糖豆9988

sobereva 发表于 2023-3-25 03:54
服务器上chk没转换完整或者转换失败，或者下载的时候没下载全。如果你确认不存在这些问题，把gjf文件和mo ...

谢谢社长，解决啦~

wzkchem5

糖糖糖豆9988 发表于 2023-3-25 04:12
请教各位老师，用orca计算三重态的O2，如下的设置正确吗？
%pal nprocs 8
end

如果因为体系里其他物种太大，DLPNO-CCSD(T)/cc-pVTZ以上的级别算不动的话，可以