Need help for a CBS extrapolation
Hi. I need help with a CBS extrapolation. Can't ask anyone for help bc nobody I know can help me. Therefore:
I did DLPNO-CCSD(T1) calculations with orca (see sull inputs below). My system is a C-C dimer which breaks into two identical C radicals and I want to benchmark the electronic reaction energy.
I used cc-pVDZ and cc-pVQZ basis sets for the extrapolation and will extrapolate SCF energy (parameter 4.42; exponential fit) and correlation energy (parameter 2.46; polynomial fit) separately. I will use the formulas and fitting parameters of
Frank Neese, Edward F. Valeev; Revisiting the Atomic Natural Orbital Approach for Basis Sets: Robust Systematic Basis Sets for Explicitly Correlated and Conventional Correlated ab initio Methods?. J. Chem. Theory Comput. 11 January 2011; 7 (1): 33–43. https://doi.org/10.1021/ct100396y.
Now to the problem. I did the calculations with ORCA (input see below) and when I checked if E(SCF)+E(CCSD)+E(T)=E(final single point) I found that the value was off by -0.03378. It seems that for the MDCI ORCA did not use the final SCF energy but a "reference energy" which is exactly -0.03378 off of the SCF energy. (SCF energy: -585.730392772957; E(0) of the MDCI: -585.696614349) My LLM of choice could not sufficiently explain, where this comes from. The whole ORCA output file does not contain this exact value (I searched for 0.03378 and -0.03378).
Now to the questions:
- What is this energy difference? Why is it here and where does it come from?
- Which correlation energy should I use for the extrapolation? I would choose E(CCSD)+E(T) (values see below).
- Which SCF energy should I use for the extrapolation? I would have used the "reference energy" E(0) which was used in the MDCI (values see below).
Examples for a radical monomer calculation:
Total SCF Energy:
----------------
TOTAL SCF ENERGY
----------------
Total Energy : -585.73039277295754 Eh -15938.53429 eV
...
Correlation energy part 1 (CCSD iterations):
--- The CCSD iterations have converged ---
E(0) ... -585.696614349
E(CORR)(strong-pairs) ... -1.982761334
E(CORR)(weak-pairs) ... -0.002944788
E(CORR)(corrected) ... -1.985706122
E(TOT) ... -587.682320471
Singles norm <S|S>**1/2 ... 0.154748668 ( 0.062079079, 0.092669589)
T1 diagnostic ... 0.018629550
<S**2>(linearized) ... 0.7533235 (ideal value: 0.7500000)
Correlation energy part 2 (triples correction):
Triples Correction (T) ... -0.086625089
alpha-alpha-alpha ... -0.002628018 ( 3.0%)
alpha-alpha-beta ... -0.041032413 ( 47.4%)
alpha-beta -beta ... -0.040372285 ( 46.6%)
beta -beta -beta ... -0.002592373 ( 3.0%)
Final correlation energy ... -2.072331210
E(CCSD) ... -587.682320471
E(CCSD(T)) ... -587.768945560
ORCA Input for a radical monomer:
!UHF DLPNO-CCSD(T1) cc-pVDZ cc-pVDZ/C TightPNO
%scf
maxiter 1000
end
%pal
nprocs 12
end
%maxcore 8500
*xyzfile 0 2 a0of.xyz
ORCA Input for a closed shell dimer:
!DLPNO-CCSD(T1) cc-pVDZ cc-pVDZ/C TightPNO
%mdci
UseFullLMP2Guess false
end
%scf
maxiter 1000
end
%pal
nprocs 12
end
%maxcore 8500
*xyzfile 0 1 A0of.xyz
Feel free to criticize my input files and point out bad practices. I like to learn and improve. And thanky for the help :)