关于CCSD(T)的一个小问题

圣人惠

请问各位，我想用CCSD（T）算单点，可是这个有赝势的不知道该如何处理，看了sob老师的文章也还是不太清楚，想请各位指教，谢谢大家

Warm_Cloud

%mem=2000MB
%chk=augHIO3NH3LW2.chk
%nproc=4
#p CCSD(T)/genecp
   
Title Card Required

0 1
53           0       -0.146672    0.013541   -0.203279
  8           0        0.159851    1.459190    0.829371
  8           0        0.094495   -1.467559    0.801469
  8           0       -2.112714    0.010299   -0.169953
  1           0       -2.398631   -0.272404    0.708647
  7           0        2.407672   -0.044012   -0.198535
  1           0        2.564434   -0.900194    0.323162
  1           0        2.566150    0.744979    0.420318
  1           0        3.054884    0.002568   -0.975708

O,H,N,0
aug-cc-pVDZ
****
I 0                           
S 9 1.00                     
0.24497900E+04 0.41900000E-03     
0.35980800E+03 0.22400000E-02     
0.14405800E+02 0.39722300E+00     
0.90763200E+01 -0.93224900E+00   
0.20881000E+01 0.93713800E+00     
0.10349800E+01 0.39208600E+00     
0.31628400E+00 0.12485000E-01
0.12171900E+00 -0.13290000E-02   
0.42000000E-01 0.00000000E+00     
S 9 1.00                          
0.24497900E+04 0.17500000E-03     
0.35980800E+03 0.10570000E-02     
0.14405800E+02 0.16900000E+00     
0.90763200E+01 -0.42179300E+00
0.20881000E+01 0.63886400E+00     
0.10349800E+01 0.32011500E+00
0.31628400E+00 -0.81442800E+00   
0.12171900E+00 -0.48979800E+00
0.42000000E-01 0.00000000E+00     
S 1 1.00                     
0.31628400E+00 1.0               
S 1 1.00                          
0.12171900E+00 1.0               
S 1 1.00                          
0.42000000E-01 1.0               
P 7 1.00                     
0.19530100E+02 0.58934000E-01     
0.11088200E+02 -0.23093000E+00   
0.27156300E+01 0.66480100E+00     
0.12043000E+01 0.45067300E+00     
0.33994500E+00 0.28980000E-01     
0.11088100E+00 -0.28890000E-02
0.33800000E-01 0.00000000E+00     
P 7 1.00                     
0.19530100E+02 -0.18836000E-01   
0.11088200E+02 0.80006000E-01
0.27156300E+01 -0.30665200E+00   
0.12043000E+01 -0.14759400E+00   
0.33994500E+00 0.60750600E+00     
0.11088100E+00 0.54704900E+00     
0.33800000E-01 0.00000000E+00     
P 1 1.00                     
0.11088100E+00 1.0               
P 1 1.00                     
0.33800000E-01 1.0            
D 6 1.00                           
0.45476500E+02 0.42660000E-02      
0.13192800E+02 -0.13625000E-01      
0.42274100E+01 0.30975600E+00      
0.19428000E+01 0.50977200E+00      
0.83977100E+00 0.29746100E+00      
0.30000000E+00 0.40164000E-01      
D 1 1.00                           
0.30000000E+00 1.0                  
D 1 1.00                           
0.11910000E+00 1.0                              
****

I 0                  
ECP28MDF 4 28         
G-Komponente         
1                     
2 1.000000 0.000000         
S-G                  
3                     
2 40.033376 49.989649        
2 17.300576 281.006556      
2 8.851720 61.416739         
P-G                  
4                     
2 15.720141 67.416239        
2 15.208222 134.807696      
2 8.294186 14.566548         
2 7.753949 28.968422         
D-G                  
4                     
2 13.817751 35.538756        
2 13.587805 53.339759        
2 6.947630 9.716466         
2 6.960099 14.977500         
F-G                  
4                     
2 18.522950 -20.176618      
2 18.251035 -26.088077      
2 7.557901 -0.220434         
2 7.597404 -0.221646      

圣人惠

Warm_Cloud 发表于 2017-5-24 21:34
%mem=2000MB
%chk=augHIO3NH3LW2.chk
%nproc=4

谢谢，非常感谢

圣人惠

Warm_Cloud 发表于 2017-5-24 21:34
%mem=2000MB
%chk=augHIO3NH3LW2.chk
%nproc=4

这个是gaussian的吧，我想请教你会不会写orca的输入呢

Warm_Cloud

molpro对avdz-pp(aug-cc-pvdz-pp)支持的好。

圣人惠

Warm_Cloud 发表于 2017-5-24 22:17
molpro对avdz-pp(aug-cc-pvdz-pp)支持的好。

老师认准了orca了，我回头和老师建议建议

liyuanhe211

看ORCA 3.0.3 手册 6.3.5 、6.3.6 节（或4.0.0手册的对应段落），很清楚。

圣人惠

liyuanhe211 发表于 2017-5-25 00:13
看ORCA 3.0.3 手册 6.3.5 、6.3.6 节（或4.0.0手册的对应段落），很清楚。

谢谢，非常感谢