“…Several treatments, including Loudon's [1], regularize the potential by introducing a small cutoff a via V(x) = γ /(|x| + a), and taking the limit a → 0 at the end. This regularization procedure has led to conflicting results in which concerns the degeneracy or not of the spectrum, the existence of even wave functions and whether the associated spectrum is continuous or discrete, and the stability of the model (unboundedness of the ground state energy from below, in the a → 0 limit) [1,6,9,10]. This is, perhaps, unsurprising, since regularization procedures need to be employed with great care and often need additional input, such as symmetry, to yield meaningful results, as is well-known in quantum field theory [11] or in the study of chemical indices [12,13]-see [14,15] for rigorous treatments of the regularized 1D Coulomb potential, and [16] for a general treatment of regularized Sturm-Liouville operators.…”