Let $\lambda \geq 1$ and $\varepsilon \leq \frac{1}{2}\sqrt{\lambda}$.
Why is it the case that $\exists C$ a universal constant (i.e. that does not depend on $\lambda, \varepsilon$) such that $\forall i \in \mathbb{N}, \forall t \in [\lambda - \varepsilon , \lambda + \varepsilon ]$, $$e^{-t} t^i \leq \frac{1}{2} C \left( e^{-(\lambda - \varepsilon)}(\lambda - \varepsilon)^i + e^{-(\lambda + \varepsilon)}(\lambda + \varepsilon)^i \right) $$
The careful reader will have noticed that this is a property about the pmf of the Poisson distribution. For information, this statement appears in "An Automatic Inequality Prover and Instance Optimal Identity testing -- G. Valiant & P. Valiant, Lemma 7" , but it isn't as obvious to me as it is to them.
Here is were we are at so far:
Write $f_i(t) = e^{-t}t^i$ for simplicity of notation, we then want to show that there exists $C$ universal such that $f_i(t) \leq \frac{1}{2} C \left( f_i(\lambda - \varepsilon) + f_i(\lambda + \varepsilon) \right) $.
Notice that since $f_i'(t) = f_{(i-1)}(t)(i-t)$, $f_i$ is increasing on $[0, i]$ and decreasing on $[i, \infty]$.
- For $i < \lambda - \varepsilon$ it is the case that $f_i(t) \leq f_i(\lambda - \varepsilon)$ and $C = 2$ is sufficient.
- For $i > \lambda + \varepsilon$ it is the case that $f_i(t) \leq f_i(\lambda + \varepsilon)$ and $C = 2$ is sufficient as well.
- For $i \in [\lambda - \varepsilon, \lambda + \varepsilon]$, $\max_{t \in [\lambda - \varepsilon, \lambda + \varepsilon]} f_i(t) = f_i(i) = e^{-i}i^i = e^{i(\ln(i) - 1)} := M_i$. Now I can bound $e^{-i}i^i \leq e^{-(\lambda - \varepsilon)}(\lambda+\varepsilon)^i$ from the definition of the range of $i$, but $e^{-(\lambda - \varepsilon)}(\lambda+\varepsilon)^i = e^{3\varepsilon}e^{-(\lambda + \varepsilon)}(\lambda+\varepsilon)^i$ so this does not control it by a constant, rather a function of $\varepsilon$.
Could this come from some convexity argument ? Consider the function $h: x \mapsto e^{-x}x^x$, and compute its second derivative $h''(x) = e^{-x} x^{x-1} [x \ln^2(x) + 1]$ which is positive on the range of interest, yielding convexity of $h$.
I'll change the notation up a bit, otherwise I'll get confused.
For $n=1,2,\dots,$ let $f_n(x) = e^{-x}x^n, x\ge 0.$ Verify that $f_n$ is increasing on $[0,n],$ decreasing thereafter, with a global maximum value of $e^{-n}n^n.$
As you indicated, we only need to worry about $n\in[\lambda-\epsilon, \lambda+\epsilon].$ Suppose first $n\in[\lambda-\epsilon, \lambda].$ I'll show
$$\tag 1\frac{e^{-n}n^n}{e^{-(\lambda-\epsilon)}(\lambda-\epsilon)^n}$$
is bounded independent of $n\in [\lambda-\epsilon,\lambda],$ $\lambda\ge 1$ and $\epsilon \le \sqrt \lambda /2.$
Now $n= \lambda-\epsilon +\delta$ for some $\delta \in [0,\epsilon].$ Substitute that into $(1)$ to get
$$\frac{e^{-(\lambda-\epsilon +\delta)}(\lambda-\epsilon +\delta)^{\lambda-\epsilon +\delta}}{ {e^{-(\lambda-\epsilon)}(\lambda-\epsilon)^{\lambda-\epsilon +\delta}}}$$
$$\tag 2= e^{-\delta} \cdot \left (1+\frac{\delta}{\lambda -\epsilon}\right)^{\lambda -\epsilon}\cdot\left (1+\frac{\delta}{\lambda -\epsilon}\right)^{\delta}.$$
We need a few lemmas:
Lemma 1: If $0\le u\le v,$ then
$$ \left(1+\frac{u}{v} \right )^v \le e^u.$$
Lemma 2: If $u\ge 1,$ then
$$\left(1+\frac{u/2}{u^2-u/2}\right)^{u/2} \le e^{1/2}.$$
I'll leave the proofs of these to you for now. Now Lemma 1 implies the second factor in $(2)$ is $\le e^{\delta}.$ As for the third factor, it only gets larger if we take $\delta = \epsilon=\sqrt \lambda/2.$ Lemma 2 then implies, with $u=\sqrt \lambda,$ that the third factor is no more than $e^{1/2}.$
So we have shown the desired boundedness in $(1),$ assuming $n\in[\lambda-\epsilon, \lambda].$ Surely the same result, using similar ideas, holds for $n\in[\lambda,\lambda+\epsilon].$ I haven't checked the latter super carefully, so it might be a good exercise for you. Assuming everything works, we have shown the desired universal $C$ exists. (And probably there is a simple estimate for $C$ that can be found.)