In the Derivation of CAVI, it rewrites the ELBO using iterated expectation. Then, it becomes the negative KL divengence between $q_j(z_j)$ and $q^*_j(z_j)$. [Variational Inference: A Review for Statisticians][1] [1]: https://arxiv.org/abs/1601.00670