计算化学公社

标题: orca 和 gaussian 算出来的极化率不同 [打印本页]

作者
Author: ulosggs 时间: 2017-12-5 12:28
标题: orca 和 gaussian 算出来的极化率不同
计算乙醚的极化率，版本分别是 gaussian 16 和 orca 4.0.1.2MP2能量相同: gauss: -233.1856357 orca: -233.1856341
极化率结果不同: gauss: 59.98 orca: 56.56

gauss输入文件

%nprocshared=16
%mem=8GB
# MP2 aug-cc-pVTZ polar

Title Card Required

0 1
C 4.620  30.690 7.940
C 5.890  30.080 8.600
O 5.810  29.910  10.010
C 6.980  29.320  10.680
C 6.690  29.120  12.180
H 4.810  30.870 6.880
H 4.330  31.610 8.460
H 3.760  30.010 8.030
H 6.130  29.120 8.130
H 6.660  30.830 8.390
H 7.120  28.350  10.190
H 7.890  29.910  10.550
H 6.470  30.080  12.660
H 5.880  28.400  12.300
H 7.580  28.620  12.590

orca输入文件

! MP2 aug-cc-pVTZ

%maxcore 800

%pal
nprocs 16
end

%elprop
Polar 1
end

*xyz 0 1
C 4.620  30.690 7.940
C 5.890  30.080 8.600
O 5.810  29.910  10.010
C 6.980  29.320  10.680
C 6.690  29.120  12.180
H 4.810  30.870 6.880
H 4.330  31.610 8.460
H 3.760  30.010 8.030
H 6.130  29.120 8.130
H 6.660  30.830 8.390
H 7.120  28.350  10.190
H 7.890  29.910  10.550
H 6.470  30.080  12.660
H 5.880  28.400  12.300
H 7.580  28.620  12.590
*

作者
Author: 冰释之川 时间: 2017-12-5 12:48
本帖最后由冰释之川于 2017-12-5 12:49 编辑

注意：ORCA里的SCF默认收敛限很宽松。在高斯里SCF=conver=8 为默认；x相同条件对应于ORCA应该是：

%scf
TolE 1e-10
TolRMSP 1e-10
TolMaxP 1e-8
end

复制代码

作者
Author: beefly 时间: 2017-12-5 12:53
看你想计算什么级别的偶极矩。根据你的输入，gaussian偶极矩是hf级别的，orca需要仔细读输出信息

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