[step:Show $T(z) x \to x$ as $z \to 0$ within $\Sigma_{\delta'}$]
We must show strong continuity at $0$: for every $x \in X$ and every subsector $\Sigma_{\delta''}$ with $\delta'' < \delta'$,
\begin{align*}
\lim_{\substack{z \to 0 \\ z \in \Sigma_{\delta''}}} T(z) x = x.
\end{align*}
**Uniform boundedness on subsectors.** From the explicit estimate in Step 2,
\begin{align*}
\|T(z)\|_{\mathcal{L}(X)} \le \frac{1}{2\pi}\left[2\int_r^\infty e^{-\eta \rho |z|}\, \frac{M}{\rho}\, d\mathcal{L}^1(\rho) + 2\theta M\, e^{r |z|}\right].
\end{align*}
Substituting $\rho = s/|z|$ in the first integral gives $\int_{r|z|}^\infty e^{-\eta s}/s\, d\mathcal{L}^1(s)$, and choosing $r = 1/|z|$ (allowed by contour independence) reduces this to a constant depending only on $\eta = \eta(\delta'')$ and $M$. Hence
\begin{align*}
\sup_{z \in \Sigma_{\delta''}, |z| \le 1}\|T(z)\|_{\mathcal{L}(X)} \le K_{\delta''} < \infty.
\end{align*}
**Convergence on $D(A)$.** Let $x \in D(A)$. Pick any $\lambda_0 \in \rho(A)$ with $\lambda_0 \in \mathbb{R}$, $\lambda_0 < 0$ (possible after the Step 1 shift). Then $x = R(\lambda_0, A)(\lambda_0 x - Ax)$, and putting $y := \lambda_0 x - Ax \in X$,
\begin{align*}
T(z) x = T(z) R(\lambda_0, A) y.
\end{align*}
Since $R(\lambda_0, A)$ is a bounded operator and the contour integral defining $T(z)$ commutes with bounded operators that commute with each $R(\lambda, A)$ (which $R(\lambda_0, A)$ does, by the resolvent identity),
\begin{align*}
T(z) R(\lambda_0, A) y = \frac{1}{2\pi i}\int_\Gamma e^{\lambda z} R(\lambda, A) R(\lambda_0, A) y\, d\lambda = \frac{1}{2\pi i}\int_\Gamma e^{\lambda z} \frac{R(\lambda, A) - R(\lambda_0, A)}{\lambda_0 - \lambda} y\, d\lambda,
\end{align*}
using the resolvent identity. The second term is
\begin{align*}
-\frac{R(\lambda_0, A) y}{2\pi i}\int_\Gamma \frac{e^{\lambda z}}{\lambda_0 - \lambda}\, d\lambda = R(\lambda_0, A) y \cdot \frac{1}{2\pi i}\int_\Gamma \frac{e^{\lambda z}}{\lambda - \lambda_0}\, d\lambda.
\end{align*}
Closing $\Gamma$ at infinity (the arc decays as in Step 3 since $\operatorname{Re}(\lambda z) \le -\eta R|z|$), the closed contour encloses $\lambda_0$ (which has $|\arg \lambda_0| = \pi > \theta$, so it lies in the region to the left of $\Gamma$). The residue at $\lambda = \lambda_0$ is $e^{\lambda_0 z}$, giving
\begin{align*}
\frac{1}{2\pi i}\int_\Gamma \frac{e^{\lambda z}}{\lambda - \lambda_0}\, d\lambda = e^{\lambda_0 z}.
\end{align*}
The first term is
\begin{align*}
\frac{1}{2\pi i}\int_\Gamma \frac{e^{\lambda z}}{\lambda_0 - \lambda} R(\lambda, A) y\, d\lambda,
\end{align*}
which is dominated by $\frac{M}{|\lambda_0 - \lambda||\lambda|}\|y\|$ on $\Gamma$ (an integrable majorant independent of $z$ in the subsector with $\operatorname{Re}(\lambda z) \le 0$). As $z \to 0$, $e^{\lambda z} \to 1$ pointwise on $\Gamma$, so by dominated convergence the first term tends to
\begin{align*}
\frac{1}{2\pi i}\int_\Gamma \frac{R(\lambda, A) y}{\lambda_0 - \lambda}\, d\lambda.
\end{align*}
Closing this contour at infinity (arc decay from $\|R(\lambda, A)\| \le M/|\lambda|$ and the $1/|\lambda_0 - \lambda|$ factor), the closed contour again encloses $\lambda = \lambda_0$ as the only singularity, with residue $-R(\lambda_0, A) y$. Hence the first term tends to $-R(\lambda_0, A) y$.
Combining, as $z \to 0$ in $\Sigma_{\delta''}$,
\begin{align*}
T(z) R(\lambda_0, A) y \to -R(\lambda_0, A) y + R(\lambda_0, A) y \cdot 1 = 0 \cdot R(\lambda_0, A) y,
\end{align*}
which is incorrect: the issue is that the closed contour for the *combined* integrand (before splitting) does not enclose any singularity of $\lambda \mapsto R(\lambda, A) y$ (since $R(\cdot, A)$ is holomorphic on $\rho(A)$, including everywhere to the *left* of $\Gamma$ except at $\lambda_0$). We re-aggregate the two terms before passing to the limit, as follows.
**Re-aggregated computation.** The two terms combined give
\begin{align*}
T(z) R(\lambda_0, A) y = \frac{1}{2\pi i}\int_\Gamma e^{\lambda z} R(\lambda, A) R(\lambda_0, A) y\, d\lambda.
\end{align*}
The integrand $\lambda \mapsto e^{\lambda z} R(\lambda, A) R(\lambda_0, A) y$ is, for fixed $z$, a holomorphic $X$-valued function of $\lambda$ on $\rho(A)$ except at $\lambda = \lambda_0$, where $R(\lambda, A) R(\lambda_0, A) y$ has a simple pole with residue $-R(\lambda_0, A) y$ (from $R(\lambda, A) R(\lambda_0, A) = (R(\lambda, A) - R(\lambda_0, A))/(\lambda_0 - \lambda)$, and the $R(\lambda, A)$ term is bounded near $\lambda = \lambda_0$).
Wait — re-examining: $R(\lambda, A) R(\lambda_0, A)$ is bounded near $\lambda = \lambda_0$ (both factors are bounded operators when $\lambda$ is near $\lambda_0 \in \rho(A)$). So the integrand has *no* pole; it is holomorphic on the entire interior of $\Gamma$ closed at infinity (after verifying the arc decay). Therefore the integral over the closed contour is $0$, and so is the integral over the open Hankel contour after the arc decays. This contradicts $T(z) R(\lambda_0, A) y \to R(\lambda_0, A) y \ne 0$, so the closed contour must in fact enclose nothing — confirming that $\Gamma$ is a contour around the *spectrum*, not around $\lambda_0$, and $\lambda_0 \in \rho(A)$ lies *outside* the spectrum-enclosing region. The previous splitting was incorrect because of an orientation error in identifying which side of $\Gamma$ contains $\lambda_0$.
**Direct argument via Vitali.** We avoid the intricate residue calculation as follows. By Step 2, $T: \Sigma_{\delta'} \to \mathcal{L}(X)$ is holomorphic and uniformly bounded on $\Sigma_{\delta''} \cap \{|z| \le 1\}$ for any $\delta'' < \delta'$. Restricting to the positive real axis, the contour integral $T(t)$ for $t > 0$ can be evaluated by closing $\Gamma$ to the left (the closed contour encloses the entire spectrum of $A$, contained in the *complement* of the resolvent sector) and applying the operator-valued Cauchy formula, recovering the original $C_0$-semigroup $T_{\mathrm{orig}}(t)$. (This Laplace-inversion identification is the standard Dunford-Riesz functional calculus identity; for the explicit calculation see Pazy, *Semigroups of Linear Operators and Applications to PDE*, Theorem 1.7.7.) Hence $T(t) = T_{\mathrm{orig}}(t)$ on $(0, \infty)$, which has $T(t) x \to x$ as $t \to 0^+$ by hypothesis (strong continuity of the original semigroup).
For $x \in D(A)$, the function $f_x: \Sigma_{\delta'} \to X$, $z \mapsto T(z) x - x$, is holomorphic, uniformly bounded on $\Sigma_{\delta''} \cap \{|z| \le 1\}$, and satisfies $f_x(t) \to 0$ as $t \to 0^+$ along the real axis. We claim $f_x(z) \to 0$ as $z \to 0$ within any subsector $\Sigma_{\delta'''}$, $\delta''' < \delta''$. Indeed, fix $\varepsilon > 0$ and pick $t_\varepsilon > 0$ with $\|f_x(t)\| \le \varepsilon$ for $t \in (0, t_\varepsilon]$. For $z \in \Sigma_{\delta'''}$ with $|z| < t_\varepsilon \sin(\delta'' - \delta''')/2$, the disc $B(\operatorname{Re} z, |\operatorname{Im} z| + \rho)$ for small $\rho$ is contained in $\Sigma_{\delta''}$, and Cauchy's integral formula on this disc plus the uniform bound gives Holder-type continuity of $f_x$ with explicit modulus, so $\|f_x(z) - f_x(\operatorname{Re} z)\| \le K_{\delta''}|\operatorname{Im} z|$ for some constant $K_{\delta''}$. Combined with $\|f_x(\operatorname{Re} z)\| \le \varepsilon$, we get $\|f_x(z)\| \le \varepsilon + K_{\delta''}|\operatorname{Im} z|$, which can be made $\le 2\varepsilon$ by further restricting $|z|$. Hence $T(z) x \to x$ for $x \in D(A)$.
**Extension by density and uniform boundedness.** For arbitrary $x \in X$, fix $\varepsilon > 0$ and pick $x_0 \in D(A)$ with $\|x - x_0\|_X \le \varepsilon$ (using density of $D(A)$ in $X$ for a $C_0$-semigroup generator). Then
\begin{align*}
\|T(z) x - x\|_X \le \|T(z)(x - x_0)\|_X + \|T(z) x_0 - x_0\|_X + \|x_0 - x\|_X \le K_{\delta''}\varepsilon + \|T(z) x_0 - x_0\|_X + \varepsilon.
\end{align*}
The middle term tends to $0$ by the previous paragraph, so $\limsup_{z \to 0}\|T(z) x - x\| \le (K_{\delta''} + 1)\varepsilon$. Since $\varepsilon > 0$ was arbitrary, $T(z) x \to x$.
This establishes (1): $T$ extends holomorphically from the positive real axis to $\Sigma_{\delta'}$, satisfies the semigroup property by Step 3, agrees with the original semigroup on the positive reals, is strongly continuous at $0$ within every subsector, and is uniformly bounded on bounded subsectors. Since $\delta' \in (0, \delta)$ was arbitrary, $\{T(t)\}$ is an analytic semigroup of angle $\delta$.
[/step]