orca 和 gaussian 算出来的极化率不同

ulosggs

计算乙醚的极化率，版本分别是 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
*

冰释之川

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

1. %scf
   
2.     TolE 1e-10
   
3.     TolRMSP 1e-10
   
4.     TolMaxP 1e-8
   
5. end

复制代码

beefly

看你想计算什么级别的偶极矩。根据你的输入，gaussian偶极矩是hf级别的，orca需要仔细读输出信息