While working on the journal version of our paper, we encountered the following problem, which seems to be fairly simple but we could not find the answer:
Suppose $A_i \in \mathbb{R}^{p\times n}$ are i.i.d. Gaussian matrices ($p<n$), and $A_i = V_i \Sigma_i U_i^T$ is the reduced SVD of the matrix. Then taking these $U_i$, we can form $\sum_{i=1}^m U_i U_i^T$. How large $m$ needs to be for this to be full rank?
I think the $U_i$ formed this way follow a Haar measure (meaning its uniformly distributed on Stiefold manifold).
To get a full rank matrix, the size of m should be at least $p$.