Bound states and critical behavior of the Yukawa potential

We investigate the bound states of the Yukawa potential V (r)=−λexp(−αr)/r, using different algorithms: solving the Schr?dinger equation numerically and our Monte Carlo Hamiltonian approach. There is a critical α = α_C, above which no bound state exists. We study the relation between α_C and λ for various angular momentum quantum number l, and find in atomic units, α_C(l) = λ[A ₁ exp(−l/B ₁) + A ₂ exp(−l/B ₂)], with A ₁ = 1.020(18), B ₁ = 0.443(14), A ₂ = 0.170(17), and B ₂ = 2.490(180).