The commutative Gelfand-Naimark theorem tells us that every unital commutative C* algebra is isometrically isomorphic to the space of continuous functions on its maximal ideal space. The non- commutative GNS construction, on the other hand, tells us that every C* algebra is isometrically embedded into a * subalgebra of $\mathfrak{B}(\mathbb{H})$ for some Hilbert space $\mathbb{H}$.
My question is the following: how does the commutative case fit into the non commutative one?
I could see that given a commutative C* algebra $C(X)$, where $X$ is a locally compact Hausdorff space, one can find a non zero positive linear functional $F$ on $C(X)$ such that $F(e)=1$, where $e$ is the identity function on $C(X)$. Then using Riesz Representation theorem, one can find a corresponding Borel $\sigma$ algebra on $X$ and a unique probability measure $\mu$ such that $F(f)= \int f\,d\mu \quad \forall f\in C(X)$. From here it is easy to deduce that $C(X)$ is isometrically embedded into a * subalgebra of $\mathfrak{B}(\mathcal{L^2}(X,\mu))$ via multiplication operators.
My question is: Is there a canonical way of finding such a probability measure? Clearly if we choose a different positive linear functional, we get a different probability measure. So is there some unique/ natural way of doing this? Also, is this where the Hahn-Hellinger theorem comes in?
I would be grateful for some clarity on this, and also references if possible.
The commutative case fits into the noncommutative case quite simply; in both cases there is not a canonical choice of linear functional which leads to an isometric representation of the C$^*$-algebra.
Let's consider the GNS representation more carefully.
It is not the case that there is always one linear functional $f$ whose GNS representation gives you an isometric representation of your C$^*$-algebra $A$ in some $B(H)$. If the linear functional $f$ is faithful, i.e., $f(x^*x)=x$ implies $x=0$, then the corresponding representation $\pi_f$ is isometric. If your $C^*$-algebra is separable, then there will always be such a faithful positive linear functional $f$, but there will be many different choices for $f$. Of course, the C$^*$-algebras $\pi_f(A)$ will all be isometrically isomorphic for all these different choices of faithful states $f$.
In general you will need to consider a direct sum of many different GNS representations using many different states in order to obtain the isometric inclusion of $A$ into some $B(H)$. For example, one can take the "universal representation" formed by taking the direct sum of the GNS representations for all positive linear functionals in $A^*$.
Similarly, one may need more than one measure to obtain an isometric representation of your commutative C$^*$-algebra. Instead of fixing one $F \in C(X)^*$, $F \geq 0$, one could take a universal representation formed by taking the direct sum over all such $F$. So there is not a canonical way to pick one such $F$ or corresponding measure $\mu$.
Also, if $X$ is a locally compact space that is not compact, then one should instead consider $C_0(X)$, the continuous functions on $X$ that vanish at infinity. In this case, the function $e$ in your question (given by $e(x)=1$ for all $x$) will not be an element of your C$^*$-algebra.
Yes, this lack of a canonical representation of self-adjoint operators by multiplication operators makes it more difficult to determine when two such operators are unitarily equivalent and it is for this reason that we need results like the Hahn-Hellinger theorem.