A new class of orthogonal polynomials associated with Coulomb wave functions is introduced. These polynomials play a role analogous to that the Lommel polynomials have in the theory of Bessel functions. The orthogonality measure for this new class is described in detail. In addition, the orthogonality measure problem is discussed on a more general level. Apart from this, various identities derived for the new orthogonal polynomials may be viewed as generalizations of certain formulas known from the theory of Bessel functions. A key role in these derivations is played by a Jacobi (tridiagonal) matrix J(L) whose eigenvalues coincide with the reciprocal values of the zeros of the regular Coulomb wave function F-L(eta, rho). The spectral zeta function corresponding to the regular Coulomb wave function or, more precisely, to the respective tridiagonal matrix is studied as well. (C) 2014 Elsevier Inc. All rights reserved.