number of configurations $\Phi_p$, and $C_p$ a variational coefficient. All configurations are symmetry-adapted (this is expressed in the last equation with the operator $\hat{O}(\hat{L}^2)$). The operator $\hat{\mathcal{A}}$ is the antisymmetrization operator and $\chi$ is the spin eigenfunction: $$\chi =\left[ (\alpha \beta -\beta \alpha )\alpha \right].$$ The Hartree products are multiplied by up to one interelectronic coordinate $r_{ij}$ $$\psi _p=r_{ij}^\nu \prod_{k=1}^n\phi _k(r_k,\theta _k,\varphi _k),$$ where the choice $\nu = 0$ corresponds to the CI wave function, while for $\nu = 1$ we have the Hy-CI wave function. To calculate the bound (doublet) states in the C$^{3+}$ and F$^{6+}$ ions we have optimized the exponents of the wave function expansions for the Li and Be$^{+}$ atoms of Ref.[@Bel1] for the nuclear charges $Z = 6$ and $Z = 9$. The Hy-CI calculations have been conducted using quadruple precision arithmetic and we have employed $\approx$ 1000 configurations, while in the CI calculations double precision has been used and we have employed between 1000 and 1500 configurations. The construction of the configurations is described in [@Bel1]. The configurations of the newly calculated states in this work are similar to the ones of [@Bel1] employing higher $l$ quantum number for the outer electron. The CI 