Let $a_1, \cdots, a_n\in\mathbb{R}^k$ be independent random vectors sampled from $N(0,\Sigma)$, where $\Sigma = \operatorname{diag}(\lambda_1, \cdots, \lambda_k)$ and $\lambda_1 \ge \cdots \ge \lambda_k$. We aim to establish a high probability bound on the largest eigenvalue or, equivalently, the spectral norm of the matrix $\sum_{i=1}^n \|a_i\|^2 a_ia_i^T$.
Can we show there exist positive constants $c_1$ and $c_2$ such that $\mathbb{P}(\|\sum_{i=1}^n \|a_i\|^2 a_ia_i^T\| \le c_1n\lambda_1) \ge c_2$?
For example, we can show the above result for the matrix $\sum_{i=1}^n a_ia_i^T$ using the bound on the norm of matrices with Gaussian entries.
Thanks!
Update: Can we get a bound of order $O(n\lambda_1)$ without $k$ if $k = O(n)$?