I'm trying to find the gradient of the joint MGF of the multivariate normal distribution, $M_{\pmb{X}}(\pmb{t}) = \exp\left[\pmb{t}^\mathsf{T} \pmb{\mu} + \frac{1}{2}\pmb{t}^\mathsf{T}\Sigma \pmb{t}\right]$, and am having trouble finding clues online.
I understand how to derive the MGF and that the gradient should be
$\nabla M_{\pmb{X}}(\pmb{t}) = \begin{bmatrix}\frac{\delta M_{\pmb{X}}(\pmb{t})}{\delta t_1}\\ \vdots \\ \frac{\delta M_{\pmb{X}}(\pmb{t})}{\delta t_n} \end{bmatrix}$
but can't figure out how to express the partial derivatives in terms of $\mu$, $\Sigma$, and $t$. Differentiating $\exp\left[\pmb{t}^\mathsf{T} \pmb{\mu} + \frac{1}{2}\pmb{t}^\mathsf{T}\Sigma \pmb{t}\right]$ with respect to each $t$ seems to lead to some really messy algebra.
Since $\pmb{t}^\mathsf{T} \pmb{\mu}$, and $\pmb{t}^\mathsf{T}\Sigma \pmb{t}$ are scalar is it possible to write $M_{\pmb{X}}(\pmb{t}) = \exp\left[\pmb{t}^\mathsf{T} \pmb{\mu}\right] \cdot \exp\left[\frac{1}{2}\pmb{t}^\mathsf{T}\Sigma \pmb{t}\right]$ and use the product rule for two scalar valued $n$-variable functions? I.e.
$\nabla M_{\pmb{X}}(\pmb{t}) = \nabla [f(\pmb{t}) \cdot g(\pmb{t})] = f(\pmb{t}) \cdot \nabla g(\pmb{t}) + g(\pmb{t}) \cdot \nabla f(\pmb{t})$
where $f(\pmb{t}) = \exp\left[\pmb{t}^\mathsf{T} \pmb{\mu}\right]$ and $g(\pmb{t}) = \exp\left[\frac{1}{2}\pmb{t}^\mathsf{T}\Sigma \pmb{t}\right]$.
However this also leads to some mess with $\nabla g$ and the covariance matrix $\Sigma$. Though I can keep trying down this road if this is the best approach.
Any help/hints would be greatly appreciated, thanks!