calculation of DBOC

I'm using MRCC with Cfour 2.1 (I also downloaded all the available patches for MRCC). I noticed that the Diagonal Born-Oppenheimer Correction (aka DBOC or BODC or adiabatic correction) given by MRCC differ slightly from the ones computed by Cfour. 
For example, consider these input for a streched geometry of the HF molecule.
The Cfour input (using Cfour for everything) is:
It only takes a few seconds to run. It gives
RHF =       -99.849794298790385
CCSD=       -100.09088605228732
The total diagonal Born-Oppenheimer correction (DBOC) is:        627.971742 cm-1

The input for MRCC is the same but it contains the line  . I obtain:
Total CCSD energy [au]:            -100.090886055013 (same as above)
The total diagonal Born-Oppenheimer correction (DBOC) is:        631.448228 cm-1

So the value given by MRCC is larger by about 3.5 cm-1 than the one by Cfour. Not a huge difference, but it's there. Decresing the bond length this difference MRCC-Cfour also decreases, and is 0.33 cm-1 at r=1.2 a0.

If MRCC and Cfour use slightly different ways to compute the DBOC, which one should be more trustworthy? For very accurate work these small differences may be not be negligible.

BTW, the same output as above but with AUG-PVTZ for the basis set and CCSDT method doesn't run successully on my machine ("executable xvdint finished with status 17152")

Any suggestion would be appreciated! 
Lorenzo
 

Dear Lorenzo,
You should try to tighten the convergence thresholds in both programs (cctol in mrcc, CC_CONV, LINEQ_CONV, SCF_CONV in Cffour). Then the correction should come closer.

Thanks a lot! I tightened the Cfour thresholds CC_CONV, LINEQ_CONV, SCF_CONV by a factor ten (from 1e-7 to 1e-, but the results changes insignificantly (from 627.971742 cm-1 to 627.971710 cm-1).
How can I change the "cctol" parameter for MRCC when using Cfour? I tried but failed...

Set INPUT_MRCC=OFF in the ZMAT file. Then the mrcc input file is not written automatically, but you should write it manually. Create a MINP file including

iface=Cfour
calc=CCSD
dens=4
dboc=on
cctol=12

and run Cfour.

Thanks a lot for your answer! It kind-of worked, but not completely (full output attached)  . 

MRCC starts and does some part of the calculation, but eventually in the "CCSD analytic second derivative calculation" section, it stops with the message
The file fort.62 is present, but it seems MRCC is trying to get something which is not there...

 

Sorry, it doesn't work in this way. The working solution is as follows.
Run a calculation with INPUT_MRCC=ON.
Edit file fort.56, change the 16th integer in the first line, e.g., to 12 (this is cctol).
Edit ZMAT and set INPUT_MRCC=OFF
Run cfour.