Let $\mathcal H$ be a separable Hilbert space with orthonormal basis $\{e_i\}_{i\in \mathbb{N} }$, and consider the linear operators $E_{ij}\equiv e_i e_j^*$, defined as $$E_{ij}v=e_i \langle e_j,v\rangle.$$ In finite-dimensions, it's clear that these are $\dim(\mathcal H)^2$ operators that span the space of all linear operators on $\mathcal H$, call it $\operatorname{Lin}(\mathcal H)$.
Can anything similar be said in infinite dimensions? That is, is $\{E_{ij}\}_{i,j\in\mathbb{N} }$ a basis for $\operatorname{Lin}(\mathcal H)$, with respect to one of its "standard" topologies (norm, strong, weak, etc)? If not, could it maybe be a basis for the subspace of bounded/compact/Hilbert-Schmidt/something-else operators?
I'm aware that $\operatorname{Lin}(\mathcal H)$ isn't a Hilbert space, and the same is true for the set of bounded operators $B(\mathcal H)$, but I think my question should still be meaningful with the notion of convergence provided in a Banach space, which these spaces should be.
I've also found that $B(\mathcal H)$ is not separable. This tells me that most likely my question has negative answer with respect to bounded operators (and therefore also linear ones).
On the other hand Hilbert-Schmidt operators do form a Hilbert space, though I'm not sure whether it's separable. So it seems possible (or even likely) that the operators $\{E_{ij}\}$ are dense in the space of Hilbert space operators (and if with respect to the natural topology in the Hilbert-Schmidt space, they'd in fact be an orthonormal basis for it). In fact, this post discusses how finite-rank operators are dense in the Hilbert-Schmidts. Which is close enough, but also different (as far as I can tell) because their construction uses a basis comprised of all possible finite-rank operators, while here I'm focusing on a basis built from a specific basis for $\mathcal H$.