 Inverse Gaussian sum closure

In the parametrization with density proportional to $y^{-3/2}\exp [-(\lambda /y+\varphi y)/2]$, the Inverse Gaussian distribution has cumulant-generating function $K(t)=\sqrt{\lambda \varphi }-\sqrt{\lambda (\varphi -2t)}$. For $n$ independent copies, the average has this same family with parameters $(n\varphi ,n\lambda )$, while the sum has parameters $(\varphi ,n^2\lambda )$. These assertions follow by multiplying moment-generating functions and scaling the density. The more usual mean parameter is $\mu =\sqrt{\lambda /\varphi }$, so the sum has mean $n\mu$ and shape $n^2\lambda$.