TABLE I: The point-coupling constants and pairing strengths of PC-PK1 set. The corresponding QCD-scaled coupling constants c$_{ln}$ are given in the last column as well. Coupling Constant Value Dimension c$_{ln}$ α$_{S}$ − 3 . 96291 × 10 − 4 MeV − 2 -1.695 β$_{S}$ 8 . 6653 × 10 − 11 MeV − 5 1.628 γ$_{S}$ − 3 . 80724 × 10 − 17 MeV − 8 -3.535 δ$_{S}$ − 1 . 09108 × 10 − 10 MeV − 4 -0.277 α$_{V}$ 2 . 6904 × 10 − 4 MeV − 2 1.151 γ$_{V}$ − 3 . 64219 × 10 − 18 MeV − 8 -0.338 δ$_{V}$ − 4 . 32619 × 10 − 10 MeV − 4 -1.097 α$_{TV}$ 2 . 95018 × 10 − 5 MeV − 2 0.505 δ$_{TV}$ − 4 . 11112 × 10 − 10 MeV − 4 -4.171 V$_{n}$ − 349 . 5 MeV fm 3 V$_{p}$ − 330 MeV fm 3 TABLE II: The calculated binding energies (in MeV) for selected spherical nuclei by PC- PK1 in comparison with the data [23] and those by DD-PC1 [15], PC-F1 [11], PC-LA [10], and NL3* [28]. The bold-faced quantities denote that the experimental values of the corresponding nuclei are used in the parametrization fitting. The root mean square (rms) deviation Δ = N ∑ i √ ( E exp i − E cal i ) $^{2}$/N and the root of relative square (rrs) deviation δ = N ∑ i √ ( E exp i − E cal i ) $^{2}$/ ( E exp i ) 2 N are respectively listed in the last two rows. Nuclei Exp. PC-PK1 DD-PC1 PC-F1 PC-LA NL3* $^{16}$O 127.619 127.280 128.527 127.691 127.407 128.112 $^{18}$O 139.806 140.223 141.145 140.028 140.356 140.504 $^{20}$O 151.370 151.962 152.790 151.606 152.228 151.955 $^{22}$O 162.026 162.285 163.141 162.054 162.665 161.990 $^{18}$Ne 132.143 132.088 132.923 132.216 132.317 132.494 22