Choice of computational protocol for carbon-lithium spin-spin coupling constants ¹ J_CLi

In this work, we tested various computational schemes for calculations of ¹ J_CLi constants with a high accuracy. On the example of six organolithium reagents (phenyllithium monomer and dimer, monomer s-butyllithium, monomers of 1- and 2-lithionaphthalenes, and a methyllithium tetramer), the following aspects are discussed: (i) the role of a model system geometry, (ii) influence of solvent effects, and (iii) the choice of a functional and basis set. Practical recommendations for calculation of ¹ J_CLi with an accuracy ±2 Hz are formulated.