a Coulomb scalar, vector and tensor interaction. They were successful in work- ing out the exact solutions as inferred from the 1D-radial Schr¨ odinger Coulomb problem as E 2 $_{−}$m 2 = − b $^{2}$( E + m ) 2 4 N 2 ; N = n r + k$_{−}$ a + 1 , (1) for the spin symmetric case (see their Eq.(12a)), and E 2 $_{−}$m 2 = − b $^{2}$( E$_{−}$ m ) 2 4 ˜ N 2 ; ˜ N = n r + k$_{−}$ a, (2) for the pseudo-spin symmetric case (see their Eq.(17a)). A proiri, these two equations are not quadratic in E (cf. e.g., Mustafa [1,2]) to imply two branches of eigenvalues. obviously, equation (1) yields E$_{−}$ m = − b $^{2}$( E + m ) 4 N 2 , (3) and equation (2) implies E + m = − b $^{2}$( E$_{−}$ m ) 4 ˜ N 2 . (4) Therefore, the energy eigenvalues should read E = ( 4 N 2 $_{−}$b 2 ) m 4 N 2 + b 2 , (5) for the spin symmetric case, and E = − ( 4 ˜ N 2 $_{−}$b 2 ) m 4 ˜ N 2 + b 2 , (6) for the pseudo-spin case. 2