If the Hamiltonian $H(t)=-\Delta+V(t)$ is time-dependent, the approach involving evolution semigroups has to be modified and below we give a proof for the existence of solutions to the Schrödinger equation with potential $V(t)$. An obvious restriction is $V(t) \in K$ from the class of Kato perturbations from before for all times. By the Kato–Rellich theorem this guarantees that the family $\{H(t)\}_t$ of Hamiltonians has a joint domain $D(H)=D(-\Delta)$. But it remains open what time-related restriction on the potential we have to choose to get a well-defined solution, i.e. a family of evolution operators $U(t,s)$ forming an evolution system.

Definition. The unitary operators $U(t,s), 0\leq s,t \leq T$, form a evolution system if

1. $U(t,t) = \mathrm{id}$
2. $U(t,r)U(r,s) = U(t,s)$
3. $U(t,s)^{-1} = U(t,s)^* = U(s,t)$
4. $\partial_t U(t,s) = -i H(t) U(t,s)$ and $\partial_s U(t,s) = i U(t,s) H(s)$ on $D(H)$

Note that the members of the family $\{H(t)\}_t$ are in general not commuting among each other or with the evolution operator. We use a technique called the stepwise static approximation¹⁾ that divides the time interval $[0,T]$ into short subintervals and takes $H(t)$ constant over each such subinterval, then passing to smaller and smaller subintervals. The relevant potential classes, also with regard to regularity issues later, are the following:

Definition. (“Sobolev–Kato space”) $K^m = \{ V \in K \mid \partial^{\alpha}V \in K, |\alpha|\leq m \}$ with the respective norm where multi-index notation is used for the (weak) partial derivative.

Definition. (“Sobolev–Kato–Lipschitz space”) $\mathrm{Lip}([0,T], K^m)$ with norm \[ \max_{t\in[0,T]} \|V(t)\|_{K^m} + \max_{s<t}\frac{\|V(t)-V(s)\|_{K^m}}{|t-s|}. \]

Proof sketch. We show that potentials $V \in \mathrm{Lip}([0,T], K)$ give a well defined Schrödinger solutions but the proof omits a few important steps that would take up too much space here.²⁾

1. Take the time interval $[0,t]$ with $t\leq T$ divided into $k$ subintervals $[t_j,t_{j+1}]$ that all decrease in length with increasing number $k$. We define the stepwise static approximation by using Stone's theorem for $k$ different static Hamiltonians: \[ U_k(t,0) = \exp(-iH(t_{k-1})(t_k-t_{k-1})) \ldots \exp(-iH(t_{0})(t_1-t_0)) =: U_k^{(k)} \ldots U_k^{(1)}. \] We show now that for $k\rightarrow \infty$ this so-defined evolution operator converges uniformly in $t$ in the strong topology of the Hilbert space $\mathcal{H}$.

2. So we show the Cauchy property of this sequence with a similar trick as in the proof of Stone's theorem: \[ (U_k(t,0)-U_l(t,0))\psi_0 = i\int_0^t U_l(t,s) (V(s_{(l)}) - V(s_{(k)})) U_k(s,0)\psi_0 \,ds \] Here the time $s_{(k)}$ is the largest time step $t_k$ in the partitioning belonging to $U_k$ smaller than $s$. Clearly it holds $s_{(l)} - s_{(k)} \rightarrow 0$ if both indices increase. So it seems plausible that $V(s_{(l)}) - V(s_{(k)})$ in the $K$-norm, but we need also $U_k(s,0)\psi_0$ converging strongly in $D(H)$ to be able to even multiply with the potential difference and use the norm estimate (here bottom).

3. By adding a big enough constant one first makes all $H(t) \geq 1$ (i.e. $\langle \varphi,H(t)\varphi \rangle \geq 1$ for all $\varphi \in D(H)$) and thus has a well-defined inverse $H(t)^{-1}:\mathcal{H} \rightarrow D(H)$ bounded by 1.

4. We add $D(H)$-identities to the definition of $U_k$: \begin{align*} U_k(t,0) &= U_k^{(k)} \ldots U_k^{(1)} \\ &= H(t_{k-1})^{-1}H(t_{k-1}) U_k^{(k)} \ldots H(t_{0})^{-1}H(t_{0}) U_k^{(1)} \\ &= H(t_{k-1})^{-1} U_k^{(k)} H(t_{k-1}) \ldots H(t_{0})^{-1} U_k^{(1)} H(t_{0}) \end{align*}