I noticed that
$$t^2 \frac{\partial^3}{\partial t^3}\Delta_t(s)+s^2 \frac{\partial}{\partial s} \Delta_t(s)=0 $$
is satisfied by
$$\Delta_t(s)= - \sqrt{\frac{t}{s}}Y_1{(4\pi\sqrt{ts})}+ \sqrt{\frac{t}{s}}K_1(4\pi \sqrt{ts})$$
for $Y_1$ and $K_1$ Bessel functions.
I noticed this looked similar to the "error term" in the Dirichlet divisor problem
$$\Delta(t)= -\sum_{s \in \Bbb N} d(s) \sqrt{\frac{t}{s}}Y_1{(4\pi\sqrt{ts})}-\sum_{s \in \Bbb N} d(s)\sqrt{\frac{t}{s}}K_1(4\pi \sqrt{ts})$$
I know from number theory that the order of magnitude of $\Delta(t)$ as $t\to \infty$ is close to $O(t^{1/4+\epsilon}).$ Making this exact is the notorious Dirichlet divsor problem.
I wonder if there are techniques from the theory of pde's that I can use to show this asymptotic holds or at least get close.
I observed that
$$t^2\frac{\partial^3}{\partial t^3}\sum_{n=1}^\infty \Psi_n(t,s)=s^2\frac{\partial}{\partial s}\sum_{n=1}^\infty \Psi_n(t,s)+\sum_{n=2}^\infty b(n)\frac{\partial}{\partial s} \Psi_n(t,s)$$
has a solution
$$\Psi_n(t,s)=2 \sqrt{\frac{tn}{s}}K_1(2\sqrt{t ns})$$
where the coefficients are
$$b(n)=\big\lbrace 2,12,24,40,60,84, \cdot\cdot\cdot\big\rbrace.$$
which I noticed correspond to twice the pronic/oblong numbers (if you throw out the first term $2$).
I thought about writing out the coefficients of the sequence $r(n)$ (to be determined) for the PDE corresponding to $Y_1$ as we did for $b(n)$ corresponding to $K_1$.
Then through linearity of solutions we can obtain a PDE involving both $r(n)$ and $b(n)$:
$$t^2\frac{\partial^3}{\partial t^3}\sum_{n=1}^\infty \chi_n(t,s)+s^2\frac{\partial}{\partial s}\sum_{n=1}^\infty \chi_n(t,s)=\sum_{n=2}^\infty r(n)\frac{\partial}{\partial s} \chi_n(t,s)+\sum_{n=2}^\infty b(n)\frac{\partial}{\partial s} \chi_n(t,s)$$
where $$\chi_n(t,s)=-2 \sqrt{\frac{tn}{s}}Y_1{(2\sqrt{tns})}+2 \sqrt{\frac{tn}{s}}K_1(2 \sqrt{tns})$$
Is $r(n)$ a known arithmetic sequence?