Time evolution of a wave packet from the time-independent Schroedinger equation [duplicate]

Question

Starting off with the time-independent Schroedinger equation (TISE)

$\quad \quad -\frac{\hbar^2}{2m} \nabla^2 \psi + V(r, \theta) \psi = E\, \psi,$

I would like to study the time evolution of an initial wave packet $\psi(0)$ dropped in the potential $V$. In particular, I want to evaluate $\langle \theta \rangle$ as a function of the time $t$.

The TISE for the potential $V$ does not have any exact solution in terms of any known functions. The wave packet $\psi (0)$ is expressed as a linear combination of a certain class of functions (infinite in number), which form a basis, but are NOT eigenstates of the Hamiltonian.

I noticed a similar question which uses the time-ordering operator. In that context, my questions were as follows:

• Given a Hamiltonian in differential form (with a particularly nasty $\nabla^2$ in polar coordinates), how does one convert it to its matrix form using Mathematica, considering that the basis used to express $\psi(0)$ has infinitely many functions?

• Is there some way to implement $\exp \big(\partial/ \partial \theta \big)$ directly as a function operating on an argument in Mathematica? As suggested by my attempts, Mathematica does not seem to recognise the dangling derivative in the exponential.

Having an answer to Part 1 of my question would enable me to generalise this answer to the problem at hand. However, I would certainly welcome any efficient approach to the direct computation of $\langle \theta \rangle (t)$ that you may suggest.

---

Edit

Just to prevent confusion about the second part of the question, by $\exp \big(\partial/ \partial \theta \big)$, I mean

$\quad \quad \exp \big(\partial/ \partial \theta \big) = 1 + \frac{\partial}{\partial \theta} + \frac{1}{2!}\frac{\partial^2}{\partial \theta^2} + \ldots$

So $\exp \big(\partial/ \partial x \big) (x^2) = 1 + 2x + x^2 \ne \exp(2x)$.

Answer 1

Try this exponential derivative operator:

expD[f_, x_] := 
 Module[{x0},
  Sum[SeriesCoefficient[f, {x, x0, i}], {i, 0, \[Infinity]}] /. {x0 -> x}
 ]

Examples:

expD[x^2, x]
(* (1 + x)^2 *)

expD[Sin[x], x]
(* Sin[1 + x] *)

expD[Exp[x], x]
(* Exp[1 + x] *)