Based on the canonical correlation analysis, we derive series representations of the probability density function (PDF) and the cumulative distribution function (CDF) of the information density of arbitrary Gaussian random vectors as well as a general formula to calculate the central moments. Using the general results, we give closed-form expressions of the PDF and CDF and explicit formulas of the central moments for important special cases. Furthermore, we derive recurrence formulas and tight approximations of the general series representations, which allow efficient numerical cal-culations with an arbitrarily high accuracy as demonstrated with an implementation in PYTHON publicly available on GITLAB. Finally, we discuss the (in)validity of Gaussian approximations of the information density.