Darboux Transformation of Diffusion Processes
Alexey Kuznetsov, Minjian Yuan
TL;DR
This work develops a diffusion-theoretic incarnation of the Darboux transformation by combining Doob's $h$-transform, Siegmund duality, and a subsequent $h$-transform to produce a transformed diffusion $\tilde Y$ whose transition density is explicitly linked to the original $Y$ via a concise formula that depends on a positive seed $h$ solving ${\rm L}_Y h=\lambda h$. The main theoretical result provides a precise identity relating $p_t^{\tilde Y}$ to $p_t^Y$ through $h$ and $\,\lambda$, under suitable boundary and finiteness assumptions, thereby extending Siegmund duality to a broader class of boundary conditions. The paper then applies the construction to five Brownian-motion–based examples, deriving closed-form transition densities and, in four cases, spectral representations, and revealing how the Darboux transform shifts or augments the spectrum (e.g., adding a new eigenfunction $1/h$ or shifting eigenvalues by $-2\lambda$). These examples connect the probabilistic transform to Sturm–Liouville theory, Krein dual strings, and propagators for Schrödinger-type problems, highlighting potential applications in mathematical finance and physics. Overall, the results significantly broaden the set of diffusions with tractable transition densities by translating classical operator techniques into stochastic-analytic tools.
Abstract
Darboux transformation of a second-order linear differential operator is a well-known technique with many applications in mathematics and physics. We study Darboux transformation from the point of view of Markov semigroups of diffusion processes. We construct the Darboux transform of a diffusion process through a combination of Doob's $h$-transform and a version of Siegmund duality. Our main result is a simple formula that connects transition probability densities of the two processes. We provide several examples of Darboux transformed diffusion processes related to Brownian motion and Ornstein-Uhlenbeck process. For these examples, we compute explicitly the transition probability density and derive its spectral representation.
