I'm interested in concentration of the following random matrix sum in spectral norm

$\frac{1}{m}\sum_{k=1}^m b_k^2\mathbf{a}_k\mathbf{a}_k^*$

Here $\mathbf{a}_k\in\mathbb{R}^n$ are i.i.d. standard normal Gaussian vectors distributed as $\mathcal{N}(0,\mathbf{I}_n)$ and $b_k$ are standard normal random variables distributed as $\mathcal{N}(0,1)$ and independent of the $\mathbf{a}_k$'s. I'm interested in showing

I'm interested in showing

$\|\frac{1}{m}\sum_{k=1}^m b_k^2\mathbf{a}_k\mathbf{a}_k^*-\mathbf{I}\|\le \delta$

holds with high probability as long as $m\ge c(\delta) n$ with $c(\delta)$ a constant depending only on $\delta$. Here for a matrix $\mathbf{X}$ the spectral norm is given by $\|\mathbf{X}\|$. I can get a looser result as long as $m\ge c \frac{n\log n}{\delta^2}$. However, I would like to get rid of the extra log factor. I should also say that for the lower tail I can also get rid of the log factor i.e. prove

$\frac{1}{m}\sum_{k=1}^m b_k^2\mathbf{a}_k\mathbf{a}_k^*\succeq(1-\delta)\mathbf{I}$

holds with high probability as long as $m\ge c \frac{n}{\delta^2}$. So really need only the upper tail. Any proof that getting rid of the extra log for the upper tail is not possible would also be welcome.