Let's say I have a probability amplitude $\psi:\Sigma\rightarrow\mathbb{C}$ for some domain $\Sigma$ (so, $\psi$ satisfies $\int_\Sigma |\psi|^2=1$). Is there a way to use $\psi$ as initial conditions for diffusion of probability values?
In particular, is there a PDE for computing a function $\psi_t$ with $\psi_0\equiv\psi$ such that the induced probability distribution $\rho_t=|\psi_t|^2:\Sigma\rightarrow\mathbb{R}^+$ satisfies the heat equation $\Delta_\Sigma\rho=\frac{\partial\rho_t}{\partial t}$ for Laplacian $\Delta_\Sigma$ (or any other diffusive equation)?
This document seems to think that it might be an open mathematical question, at least for the "Anderson model" of diffusion (see e.g. conjecture (iii) on page 30), but I'm not sure if I am reading it properly. My hope was that it would be as easy as using the Dirac operator $D$ since it looks like the "square root" of the Laplacian in the right way, but apparently things aren't so easy!
[This website has been incredibly useful for me to get my bearings in this area of physics despite my total lack of background. I sincerely apologize for posting so many questions and really appreciate everyone's support!]