Not sure this is the best place to ask for verification, but I can't seem to find a derivation anywhere else. I want to calculate $\mathbb{E}[e^{\sigma(W_t + W_s)}]$, where $W_t$ and $W_s$ are two values of a standard Brownian motion and $\sigma$ is some positive constant. WLOG assume $s < t$ and let $\mathcal{F}$ be the filtration generated by the Brownian motion. Then using standard results I get
$$ \begin{align*} \mathbb{E}[e^{\sigma(W_t + W_s)}] & = \mathbb{E}[\mathbb{E}[e^{\sigma(W_t - W_s)}e^{2\sigma W_s}| \mathcal{F}_s]] \\ & = \mathbb{E}[\mathbb{E}[e^{\sigma(W_t - W_s)}]e^{2\sigma W_s}] \\ & = \mathbb{E}[e^{2\sigma W_s}]e^{\frac{1}{2}\sigma^2(t-s)} \\ & = e^{2\sigma^2 s + \frac{1}{2}\sigma^2(t-s)} \\ & = e^{\frac{1}{2}\sigma^2(3s + t)}. \end{align*} $$
Thanks!