All calculations were performed using DIRAC10 code interfaced with our in-house coupled-cluster code . We have solved the DHF equation and generated the necessary integrals using the DIRAC code, and the extracted integrals were then used for further molecular calculations using our in-house code. The finite-size nuclear model described by the Gaussian charge distribution and the default values of the nuclear parameters of DIRAC were considered in this work unless otherwise mentioned. We have used the following uncontracted basis sets: double-zeta (DZ): dyall.cv2z for Sr and cc-pCVDZ for F, triple-zeta (TZ) basis: dyall.cv3z for Sr and cc-pCVTZ and aug-cc-pCVTZ for F. Small basis functions are generated from the large basis using the restricted kinetic balance condition. We have projected out the negative energy spectrum by employing the "no virtual pair approximation" . The bond length of SrF used for calculation is 2.075 Å . The spin and the magnetic moment of the 87 Sr nucleus are considered to be 9/2 and -1.09283 nuclear magneton, respectively. We have correlated all the electrons in the molecular calculations but excluded the virtual spinors with energy higher than a certain cutoff value and the details are given in Table .