Let $\mathcal{A}$ be a (unital) $C^\star$ algebra, and $a,b$ be two positive elements of $\mathcal{A}$. We know that $ab$ is not positive in general (in fact, it is positive iff $a$ and $b$ commute), but it is always true that $\sigma(ab)\geq 0$. I was wondering if it would be possible to approximate $ab$ by a positive element in $\mathcal{A}$, that is, if for all $\epsilon>0$ there exists some $c>0$ in $\mathcal{A}$ such that $\lVert ab - c \rVert < \epsilon$ .
My obvious first attempt was with $c=\frac{ab+ba}{2}$ . But the noncommutativity of $a$ and $b$ blocks this from working. I was wondering if functional calculus techniques might work, but I am not sure. I know that the set of hermitian elements is not dense in $\mathcal{A}$, so the approximation might actually not work as well, but I don't have a counterexample at hand. I would be ineterested in both finite and infinite dimensional settings.
(Of course, the problem is trivial if $\mathcal{A}$ is commutative.)
If you had your approximation, you would have $ab$ as a limit of positive elements. But a limit of positive elements is positive.
This can be seen most easily if you work on a faithful representation $A\subset B(H)$: if $c_n\to ab$ with $c_n\geq0$, then $$ \langle ab\xi,\xi\rangle=\lim\langle c_n\xi,\xi\rangle\geq0,\qquad\qquad \xi\in H. $$
If you want to stay abstract, an element $a\in A$ is positive if and only if $\varphi(a)\geq0$ for all states $\varphi$ (proof below). And for any state you have $$ \varphi(ab)=\lim\varphi(c_n)\geq0, $$ so $ab\geq0$.
Proof that positivity can be tested on states.
Let $a\in A$ such that $\varphi(a)\geq0$ for all states $\varphi$. As $\varphi(\operatorname{Im} a)=\operatorname{Im} \varphi(a)=0$, we get $\operatorname{Im} a=0$ (by functional calculus, as any selfadjoint attains its norm at a state) and so $a=a^*$. Knowing that $a$ is selfadjoint we can do the following. If it were not positive, it has $\lambda<0$ in its spectrum; using functional calculus, we can get a state (a character, actually) $\varphi$ with $\varphi(a)=\lambda<0$.