Portfolio optimization techniques, such as those defined under Modern Portfolio Theory (MPT), are mildly predicated on the assumption of joint normality. Even though there will be a set of portfolio weights which minimizes variance regardless of the underlying distributions, correlation
is only a complete measure of association if the joint multivariate distribution is normal; i.e., covariance is only an exhaustive measure of co-movement if the joint distributions are themselves normal. We can see this is true because the joint distribution of X and Y is defined by joint normality:
${\frac {1}{2\pi \sigma _{X}\sigma _{Y}{\sqrt {1-\rho ^{2}}}}}\iint _{X\,Y}\exp \left[-{\frac {1}{2(1-\rho ^{2})}}\left({\frac {X^{2}}{\sigma _{X}^{2}}}+{\frac {Y^{2}}{\sigma _{Y}^{2}}}-{\frac {2\rho XY}{\sigma _{X}\sigma _{Y}}}\right)\right]\,\mathrm {d} X\,\mathrm {d} Y$
Which through a proof can be show to produce:
$\sigma _{X+Y}={\sqrt {\sigma _{X}^{2}+\sigma _{Y}^{2}+2\rho \sigma _{X}\sigma _{Y}}},$
If now, we define $\omega_i \sigma^2_i=\sigma_X$, and $\omega_j \sigma^2_j=\sigma_Y$, then we get back the equation which is used as the basis of mean variance optimization of a two asset portfolio:
$\mathbb{E}[\sigma _{p}^{2}]=\omega_{i}^{2}\sigma _{i}^{2}+\omega_{j}^{2}\sigma _{j}^{2}+2\omega_{i}\omega_{j}\sigma _{i}\sigma _{j}\rho _{ij}$
So while the portfolio covariance matrix can always be computed, to the extent that underlying assets have returns which are not normal the optimization is likely to result in spuriously optimal weights.