 Weighted Neumann heat kernel with constant drift

On $(0,L)$, the differential operator $\mathcal{L}={\partial }_x^2+\beta {\partial }_x=e^{-\beta x}{\partial }_x(e^{\beta x}{\partial }_x)$ with homogeneous Neumann boundary conditions is self-adjoint in the weighted inner product with weight $e^{\beta x}$. Its eigenfunctions and nonnegative decay rates areTheir squared norms are $N_0=(e^{\beta L}-1)/\beta$ and $N_k=L[1+{\beta }^2/(4q_k^2)]/2$. The Sturm-Liouville eigenfunction expansion givesThis kernel acts against the measure $e^{\beta y}dy$, not Lebesgue measure alone. It is symmetric in $x,y$; the transition density against Lebesgue measure is $e^{\beta y}{K}_{\beta }(x,y,t)$.