Me gustaría agregar algunas cosas para complementar las respuestas de Arturo y Fabian.
Si tomas el producto externo de un vector unitario, \hat{e} (\lvert e \rangle, notación de Dirac), con su dual, \hat{e}^\ast (\langle e \rvert), obtienes una matriz que proyecta vectores en el espacio definido por el vector unitario, es decir,
$$\begin{aligned}
\mathbf{P}_e &= \hat{e} \otimes \hat{e}^\ast \ &= \lvert e \rangle \langle e \rvert \ &= \left( \begin{array}{ccc} \lVert e_{1} \rVert^2 & e_{1} e_2^\ast & \ldots \\ e_2 e_1^\ast & \lVert e_{2} \rVert^2 & \ldots \\ \vdots & \vdots & \ddots \end{array} \right)
\end{aligned} $$
donde
\mathbf{P}_e \vec{v} = (\vec{v} \cdot \hat{e}) \hat{e}.
En otras palabras, \mathbf{P}_e es un operador de proyección. Utilizando esto, puedes reescribir una matriz, \mathbf{A}, en términos de sus valores propios, \lambda_i, y eigenvectores, \lvert e_i \rangle,
\mathbf{A} = \sum_i \lambda_i \lvert e_i \rangle \langle e_i \rvert
lo cual se conoce como la descomposición espectral. A partir de esto, es evidente que cualquier eigenvector, \hat{e}_i, con un valor propio de cero no contribuye a la matriz, y para cualquier componente de vector en uno de esos espacios, \mathbf{P}_{e_i}\vec{v} = \vec{v}_i, \mathbf{A} \vec{v}_i = 0.
Esto implica que la dimensionalidad del espacio en el que \mathbf{A} opera es más pequeña que la dimensión del espacio utilizado para describirlo. En otras palabras, \mathbf{A} no tiene rango completo, y no es invertible.