KL divergence between two Gaussian distributions denoted by $\mathcal{N}(\mathbf \mu_1, \mathbf \Sigma_1)$ and $\mathcal{N}(\mathbf \mu_2, \mathbf \Sigma_2)$ is available in a closed form as:
$$\mathbb{KL}= \frac{1}{2}\left[\log\frac{|\Sigma_2|}{|\Sigma_1|} - d + \text{tr} \{ \Sigma_2^{-1}\Sigma_1 \} + (\mu_2 - \mu_1)^T \Sigma_2^{-1}(\mu_2 - \mu_1)\right]$$
from: KL divergence between two multivariate Gaussians
My question is assuming $\mathbf \mu_2=\mathbf \mu_1 +\mathbf \epsilon_\mu$ and $\mathbf \Sigma_2=\mathbf \Sigma_1 +\epsilon_\Sigma$ (each element of mean and covariance are perturbed by a small amount), and $\mathbf \Sigma_1$ and $\mathbf \Sigma_2$ have Toeplitz form, can the above formulla be simplified? If closed-from solution can not be obtained, can someone upper bound KL for this setup?