A comprehensive examination was conducted on the ring-opening polymerization (ROP) mechanism of ε-caprolactone (ε-CL) employing bis-lithium N-heterocyclic carbene complexes (bisLiNHCs) in a theoretical framework. The modeled reactions were carried out utilizing density functional theory at B3LYP level. An activated-monomer mechanism was observed, characterized by occurrence of two distinct transition states. Initiating with activation of ε-CL via two lithium atoms from bisLiNHCs catalysts, the ensuing step involved nucleophilic addition of carbonyl group stemming from the lithium benzoxide (LiBnO) initiator. The ring opening ε-CL was accomplished via classical cleavage of acyl-oxygen bond. The determination of reaction barrier heights for ε-caprolactone with bisLiNHCs catalysts involved the examination of potential energy profiles. The computed energy profiles revealed exothermic characteristics for each ROP reaction, with the identification of a rate-determining step occurring at the second transition state, which entailed cleavage of acyl-oxygen bond from monomer. Furthermore, through comparisons catalytic capabilities of bisLiNHCs derivatives and energy barrier heights, it was observed that the length and flexibility of linker chain as well as steric effects from substituent groups, significantly influenced the geometric configuration. These specific variations induce alterations in both distance and angular orientation, which are pivotal aspects influencing the chelation dynamics between lithium and oxygen atoms.