Given fixed real symmetric $D\in\mathbb{R}^{n\times n}$ with $n$ distinct eigenvalues, let $U$ be a random orthogonal matrix selected uniformly from the space of $n\times n$ orthogonal matrices, and denote the first $m$ columns of $U$ as the rectangular matrix $Q\in\mathbb{R}^{n\times m}$.

We are interested in the quantity
$$\mu=\lambda_\min(D+QQ^T) - \lambda_\min(D).$$
Of course, $0 \le \mu \le 1$. But given that $Q$ is randomly selected, the decoherence between the bases of $D$ and $Q$ causes the distribution of $\mu$ to lie strictly in between the two extremes. For sufficiently large $n$, choosing $m\in\Omega(n)$ columns *should* imply that:

 1. The event $\mu=0$ occurs with probability zero.
 2. The expectation $\mathbb{E}\{\mu\}$ is bounded from below by an absolute constant.
 3. The distribution of $\mu$ concentrates about $\mathbb{E}\{\mu\}$.

Indeed, all three statements are readily confirmed using numerical simulations. But how might we go about proving these statements?

-----
Remark 1. Since $D$ is fixed, we can, WLOG, assume that $D$ is positive definite and diagonal.

Remark 2. Many of the existing work on random orthogonal projections use standard Gaussians to approximate a few orthonormal columns. But with as many columns as $m\in\Omega(n)$, the approach is no longer valid. See Jiang, Tiefeng. "How many entries of a typical orthogonal matrix can be approximated by independent normals?" *The Annals of Probability* 34.4 (2006): 1497-1529.