If $B_t$ is a Brownian motion then using Hermite polynomials one can find that
$$1, B_t, B_t^2-t, B_t^3 - 3tB_t,...$$
are martingales.
If $X_t$ is a diffusion
$dX_t = \mu(X_t,t)dt + \sigma(X_t,t)dt$
Is it possible to create sequence of polynomials that would be martingales, but now with $X_t$ instead of Brownian motion?
I will be very grateful for any resources, articles that researched this question