For $x\in\mathbb{R}^{n}$ and $t\in[0,\infty)$, consider the linear PDE initial value problem $$\dfrac{\partial u}{\partial t} = \left(a \Delta - \dfrac{b}{|x|}\right)u, \quad u(x,0) = u_0(x)\quad\text{(given)},$$ where $\Delta$ is the standard Laplacian, $\vert\cdot\vert$ denotes the Euclidean 2-norm, and the constants $a,b>0$. I want to write its solution as $$u(x,t) = \int_{\mathbb{R}^n} K(x,y;t)u_0(y)\:{\rm{d}}y.$$
Is there a known explicit (in terms of $x,y,t,n$) expression for this kernel/Green's function $K$?
It feels like something well-known but I am having difficulty finding a reference stating/deriving such $K$. The solution may be related to time-dependent Schrödinger equation with Coulomb potential but I do not know enough background to be sure.