This question arose when I am considering the disturbed Generating Moment Function (GMF) of a Bernoulli random variable. Mathematically speaking, given a Bernoulli random variable $X$ with mean $p$, we define its disturbed GMF as $$ \varphi_\varepsilon (t) \;=\; \mathbb{E} (e^{t(X - p - \varepsilon)}) \;=\; p e^{t(1-p-\varepsilon)} + (1-p) e^{-t(p+\varepsilon)}. $$ I want to prove for any $p\in [0, 1]$ (since it is the probability), for any $\varepsilon > 0$, $$ \varphi_\varepsilon (4 \varepsilon) \le 1, $$ which is equivalent to $$ f(p, \varepsilon, \theta) \;=\; e^{\theta\varepsilon (p+\varepsilon)} - p e^{\theta\varepsilon} + p - 1 \;\ge\; 0 $$ with $\theta = 4$, $p\in [0,1]$, and $\varepsilon > 0$.
When I tried to solve this, I found that it seems for any $t\ge 0$, for any $x,y\in \mathbb{R}$, $f(p, \varepsilon, t) \ge 0$, and the equality holds when $x=0$. I wonder if there is any way to solve this problem, or the former specific case. Any comment helps.