[step:Construct resolvents on the positive and negative half-lines]
For $\lambda>0$, define the bounded [linear map](/page/Linear%20Map) $R_\lambda:H\to H$ by
\begin{align*}
R_\lambda\psi:=\int_0^\infty e^{-\lambda t}U(t)\psi\,d\mathcal{L}^1(t).
\end{align*}
The integral is Bochner integrable because $\|e^{-\lambda t}U(t)\psi\|_H=e^{-\lambda t}\|\psi\|_H$ and $e^{-\lambda t}\in L^1((0,\infty),\mathcal{L}^1)$. For $h\neq 0$, translation by $h$ in the Bochner integral gives
\begin{align*}
\frac{U(h)R_\lambda\psi-R_\lambda\psi}{h}=\frac{e^{\lambda h}-1}{h}\int_h^\infty e^{-\lambda r}U(r)\psi\,d\mathcal{L}^1(r)-\frac{1}{h}\int_0^h e^{-\lambda r}U(r)\psi\,d\mathcal{L}^1(r).
\end{align*}
Here the substitution is $r=t+h$, the measure is unchanged by translation, and the interval $(0,\infty)$ becomes $(h,\infty)$. The formula is interpreted with oriented Bochner integrals: if $h<0$, then $\int_0^h=-\int_h^0$. We separate the two limiting terms. First,
\begin{align*}
\left\|\int_h^\infty e^{-\lambda r}U(r)\psi\,d\mathcal{L}^1(r)-R_\lambda\psi\right\|_H\leq \int_0^{|h|} e^{\lambda |h|} \|\psi\|_H\,d\mathcal{L}^1(r)=|h|e^{\lambda |h|}\|\psi\|_H,
\end{align*}
so the translated tail converges in $H$ to $R_\lambda\psi$. Since $(e^{\lambda h}-1)/h\to\lambda$, the first term converges to $\lambda R_\lambda\psi$. Second, the averaged interval term satisfies
\begin{align*}
\frac{1}{h}\int_0^h e^{-\lambda r}U(r)\psi\,d\mathcal{L}^1(r)\to \psi
\end{align*}
by strong continuity of $r\mapsto e^{-\lambda r}U(r)\psi$ at $0$. Hence $R_\lambda\psi\in D(K)$ and
\begin{align*}
K R_\lambda\psi=\lambda R_\lambda\psi-\psi.
\end{align*}
Therefore
\begin{align*}
(\lambda I-K)R_\lambda\psi=\psi.
\end{align*}
So $\operatorname{Range}(\lambda I-K)=H$ for every $\lambda>0$.
For $\lambda<0$, define the bounded linear map $R_\lambda:H\to H$ by
\begin{align*}
R_\lambda\psi:=-\int_{-\infty}^0 e^{-\lambda t}U(t)\psi\,d\mathcal{L}^1(t).
\end{align*}
Since $\lambda<0$, the scalar function $t\mapsto e^{-\lambda t}$ belongs to $L^1((-\infty,0),\mathcal{L}^1)$. The same translation computation, now with the interval $(-\infty,0)$ translated to $(-\infty,h)$ and with $\int_0^h$ understood as an oriented Bochner integral for $h<0$, yields
\begin{align*}
\frac{U(h)R_\lambda\psi-R_\lambda\psi}{h}=\frac{e^{\lambda h}-1}{h}R_\lambda\psi-\frac{e^{\lambda h}}{h}\int_0^h e^{-\lambda r}U(r)\psi\,d\mathcal{L}^1(r).
\end{align*}
The first term converges to $\lambda R_\lambda\psi$ by the same dominated-convergence argument on the translated half-line, and the second term converges to $\psi$ because the oriented averages over $[0,h]$ converge to the value at $0$ by strong continuity of $U$. Thus $R_\lambda\psi\in D(K)$ and
\begin{align*}
(\lambda I-K)R_\lambda\psi=\psi.
\end{align*}
So $\operatorname{Range}(\lambda I-K)=H$ for every $\lambda<0$. In particular,
\begin{align*}
\operatorname{Range}(I-K)=H,\qquad \operatorname{Range}(I+K)=H.
\end{align*}
[/step]