A transmutation operator
boldT transmuting harmonic functions into solutions of the radial Schrödinger equation
boldSu:=()△d−qfalse(rfalse)u=0 is studied. The potential
q is assumed to be continuously differentiable, and the Schrödinger equation is considered in a star‐shaped domain
normalΩ⊂double-struckRd (with
d⩾2). Several new properties of the transmutation operator are established including the operator relation
r2△d−q(r)T=Tr2△d,
valid on
C2 functions and its boundedness on the Bergman space. A Fourier‐Jacobi series expansion of the integral transmutation kernel is derived, and with its aid, an infinite system of solutions of the radial Schrödinger equation is obtained, which is shown to be complete with respect to the uniform norm. Explicit construction of the system is derived. In the case of
normalΩ being an open ball centered in the origin, the system of solutions represents an orthogonal basis of the corresponding Bergman space.