Sums of random symmetric matrices and quadratic optimization under orthogonality constraints