[guided]We need a genuine derivative formula, not just a formal calculation. Let $(\Omega,\mathcal F,\mathbb P)$ be the probability space on which $F:(\Omega,\mathcal F)\to(\mathbb R,\mathcal B(\mathbb R))$ is defined. Define the map $M_F:I\to(0,\infty)$ by assigning to each $\gamma\in I$ the value
\begin{align*}
M_F(\gamma)=\mathbb E[e^{\gamma F}].
\end{align*}
The codomain is $(0,\infty)$ because $e^{\gamma F}>0$ almost surely and the expectation is finite by hypothesis.
Fix $\gamma_0\in I$. Since $I$ is open, there is $\varepsilon>0$ such that the closed interval $[\gamma_0-\varepsilon,\gamma_0+\varepsilon]$ is contained in $I$. The purpose of choosing this smaller interval is to create integrable exponential functions on both sides of $\gamma_0$, which will dominate the derivative of $e^{\gamma F}$.
For $\gamma\in[\gamma_0-\varepsilon/2,\gamma_0+\varepsilon/2]$, the function $x\mapsto |x|e^{\gamma x}$ is controlled by the two endpoint exponentials with an explicit constant. Define $C_{\gamma_0,\varepsilon}:=2/(e\varepsilon)$. If $x\ge0$, then $\gamma\le\gamma_0+\varepsilon/2$, so
\begin{align*}
x e^{\gamma x}\le x e^{(\gamma_0+\varepsilon/2)x}=x e^{-\varepsilon x/2}e^{(\gamma_0+\varepsilon)x}\le C_{\gamma_0,\varepsilon}e^{(\gamma_0+\varepsilon)x}.
\end{align*}
If $x<0$ and $y=-x>0$, then $\gamma\ge\gamma_0-\varepsilon/2$, so
\begin{align*}
|x|e^{\gamma x}=y e^{-\gamma y}\le y e^{-(\gamma_0-\varepsilon/2)y}=y e^{-\varepsilon y/2}e^{-(\gamma_0-\varepsilon)y}\le C_{\gamma_0,\varepsilon}e^{(\gamma_0-\varepsilon)x}.
\end{align*}
Combining the two cases gives, for every $x\in\mathbb R$,
\begin{align*}
|x|e^{\gamma x}\le C_{\gamma_0,\varepsilon}\left(e^{(\gamma_0+\varepsilon)x}+e^{(\gamma_0-\varepsilon)x}\right).
\end{align*}
Substituting $x=F(\omega)$ gives the pointwise bound
\begin{align*}
|F(\omega)|e^{\gamma F(\omega)}\le C_{\gamma_0,\varepsilon}\left(e^{(\gamma_0+\varepsilon)F(\omega)}+e^{(\gamma_0-\varepsilon)F(\omega)}\right)
\end{align*}
for every $\omega\in\Omega$ outside a null set on which $F$ may be represented arbitrarily. The right-hand side is integrable with respect to $\mathbb P$ because $\gamma_0+\varepsilon$ and $\gamma_0-\varepsilon$ belong to $I$.
Now differentiate the map $\gamma\mapsto e^{\gamma F(\omega)}$ pointwise. Its derivative is $F(\omega)e^{\gamma F(\omega)}$, and the domination just obtained is integrable locally in $\gamma$. Therefore the [Dominated Convergence Theorem](/theorems/4) permits differentiation to pass through the expectation, giving
\begin{align*}
M_F'(\gamma)=\mathbb E[F e^{\gamma F}]
\end{align*}
for every $\gamma\in I$. Since $M_F(\gamma)>0$, the logarithm is differentiable at $M_F(\gamma)$, and the chain rule gives
\begin{align*}
\Lambda_F'(\gamma)=\frac{M_F'(\gamma)}{M_F(\gamma)}=\frac{\mathbb E[F e^{\gamma F}]}{\mathbb E[e^{\gamma F}]}.
\end{align*}[/guided]