I am trying to reason to myself why is it that VAEs can approximate arbitrary probability distributions even though 𝑞𝜙(𝑧|𝑥) and 𝑝𝜃(𝑥|𝑧) are Gaussian.
I understand that the parameters are typically learned using neural networks 𝜙=𝑓𝑤(𝑥) and 𝜃=𝑔𝑤(𝑧) , where 𝑓 and 𝑔 denote arbitrary neural networks. Therefore, some ideas I had is perhaps the entire VAE is able to learn differentiable, change of variables mappings, which allows the VAE to ensure that, although 𝜙 and 𝜃 both characterize Gaussian distributions, they do so using some transformation of the original variables 𝑥 and/or 𝑧 and so effectively correspond to some arbitrary distribution rather than a Gaussian in 𝑥 or 𝑧 .
However, and I may just be unnecessarily picky here, but I can't quite convince myself why the Gaussian PDF itself, even with the change of variables mapping, will not act as a "bottleneck" in learning some rather complex non-Gaussian distributions. I was wondering if there is some universal approximation theorem in the case of VAEs, as there in feed-forward neural networks / RNNs etc.