Let $X$ and $Y$ be independent random variables such that $M_X(t) = \exp(e^t - 1)$ and $M_Y(t) = \exp(3(e^t - 1))$. What is the characteristic function $\phi_{X+Y}(t)$?