Let $\varphi$ be a complex-valued function which is analytic on $\{z \in \mathbb C : |z| \leq 2\}$, let $\gamma$ be the unit circle in the complex plane, and define
$$ F_n(z) = \int_\gamma \frac{1}{s-z} \int_{-2}^{2} e^{-nt^2} \frac{\varphi(t)}{t-s}\,dt\,ds, $$
where $|z| < 1/2$.
Main Question: Is it true that $F_n(z) \to 0$ uniformly for $|z| < 1/2$ as $n \to \infty$?
We should note that the double integral exists. Indeed, by properties of Cauchy-type integrals (see Gakhov, Boundary Value Problems or here), the function
$$ g(s) = \int_{-2}^{2} e^{-nt^2} \frac{\varphi(t)}{t-s}\,dt $$
is analytic on $\mathbb C \setminus [-2,2]$ and has continuous extensions from the upper and lower half-planes to the interval $(-2,2)$ which satisfy
$$ \lim_{\epsilon \to 0^+} g(x \pm i\epsilon) = \pm i\pi e^{-nx^2} \varphi(x) + \operatorname{P.V.} \int_{-2}^{2} e^{-nt^2} \frac{\varphi(t)}{t-x}\,dt $$
for $-2 < x < 2$. Consequently, $g(s)$ is continuous on $\gamma$ except for two jump discontinuities at $s = \pm 1$.
Idea for an approach
I have an idea for an approach which I have so far been unable to make rigorous. At the end are a couple of problems that I see with it that I would greatly appreciate some feedback on.
First I'd like to split the inner integral up into
$$ \int_{-2}^{2} = \int_{|t| < 1} + \int_{1 < |t| < 2}, $$
and so write
$$ F_n(z) = \int_\gamma \frac{1}{s-z} \int_{|t| < 1} e^{-nt^2} \frac{\varphi(t)}{t-s}\,dt\,ds + \int_\gamma \frac{1}{s-z} \int_{1 < |t| < 2} e^{-nt^2} \frac{\varphi(t)}{t-s}\,dt\,ds. \tag{1} $$
Now switch the order of integration in both integrals. The first becomes
$$ \int_\gamma \frac{1}{s-z} \int_{|t| < 1} e^{-nt^2} \frac{\varphi(t)}{t-s}\,dt\,ds = \int_{|t|<1} e^{-nt^2} \varphi(t) \int_\gamma \frac{ds}{(s-z)(t-s)}\,dt, $$
and the inner integral here is
$$ \int_\gamma \frac{ds}{(s-z)(t-s)} = 2\pi i \left(\frac{1}{t-z} - \frac{1}{t-z}\right) = 0. $$
The second integral in $(1)$ becomes
$$ \begin{align} \int_\gamma \frac{1}{s-z} \int_{1 < |t| < 2} e^{-nt^2} \frac{\varphi(t)}{t-s}\,dt\,ds &= \int_{1 < |t| < 2} e^{-nt^2} \varphi(t) \int_\gamma \frac{ds}{(s-z)(t-s)}\,dt \\ &= 2\pi i \int_{1 < |t| < 2} e^{-nt^2} \frac{\varphi(t)}{t-z}\,dt, \end{align} $$
so we conclude that
$$ F_n(z) = 2\pi i \int_{1 < |t| < 2} e^{-nt^2} \frac{\varphi(t)}{t-z}\,dt. $$
Then
$$ \begin{align} e^{n} |F_n(z)| &= 2\pi \left| \int_{1 < |t| < 2} e^{-n(t^2-1)} \frac{\varphi(t)}{t-z}\,dt \right| \\ &\leq 2\pi \cdot \operatorname{Length}(\{1 < |t| < 2\}) \cdot \sup_{1 < |t| < 2} \left( e^{-n(t^2-1)} \frac{|\varphi(t)|}{|t-z|} \right) \\ &\leq 2\pi \cdot 2 \cdot \sup_{1 < |t| < 2} \left(1 \cdot \frac{|\varphi(t)|}{1/2} \right) \\ &\leq C \end{align} $$
for some constant $C > 0$, so
$$ |F_n(z)| \leq Ce^{-n} \to 0 $$
uniformly for $|z| < 1/2$ as $n \to \infty$.
There are a couple of issues that I see with this:
Can the interchange of order of integration be justified in both cases?
Are the subsequent evaluations of the inner integrals using the residue theorem valid?