Repulsion integrals involving Slater-type functions and Yukawa potential

In range-separated functionals using basis sets of exponential functions, the Yukawa potential appears as the natural choice for the attenuated electron–electron repulsion. The development of these methods requires one- and two-center integrals involving two exponential functions and the Yukawa potential. In this work, general expressions are derived for these integrals, and stable algorithms for their efficient calculation, covering a wide range of screening parameters, are developed and tested.