“…The single-configurational methods include Hartree−Fock (HF), DFT with local density approximation (LDA) VWN functional V (donated as VWN5), with generalized gradient approximation (GGA) BLYP, , and with hybrid functional B3LYP , (in which VWN functional III and V 38 are used for local correlation, denoted, respectively, as B3LYP3 and B3LYP5), coupled-cluster with single and double excitations (CCSD), and CCSD with perturbative triples corrections, that is, CCSD(T) . The single-reference configuration of (I) 3 Δ g corresponds to nickel (3d δ ) 3 (3d σ ), while for (I) 1 A 1 , the lowest five a 1 , two b 1 and b 2 , and one a 2 molecular orbitals are fully occupied. The multi-configurational method used is internal-contracted multireference average quadratic single and double coupled-cluster (ic-MR-AQCC) with the active space of (12,8).…”