We study the unique solution $m$ of the Dyson equation \[ -m(z)^{-1} = z - a + S[m(z)] \] on a von Neumann algebra $\mathcal{A}$ with the constraint $\mathrm{Im}\,m\geq 0$. Here, $z$ lies in the complex upper half-plane, $a$ is a self-adjoint element of $\mathcal{A}$ and $S$ is a positivity-preserving linear operator on $\mathcal{A}$. We show that $m$ is the Stieltjes transform of a compactly supported $\mathcal{A}$-valued measure on $\mathbb{R}$. Under suitable assumptions, we establish that this measure has a uniformly $1/3$-Hölder continuous density with respect to the Lebesgue measure, which is supported on finitely many intervals, called bands. In fact, the density is analytic inside the bands with a square-root growth at the edges and internal cubic root cusps whenever the gap between two bands vanishes. The shape of these singularities is universal and no other singularity may occur. We give a precise asymptotic description of $m$ near the singular points. These asymptotics generalize the analysis at the regular edges given in the companion paper on the Tracy-Widom universality for the edge eigenvalue statistics for correlated random matrices [arXiv:1804.07744] and they play a key role in the proof of the Pearcey universality at the cusp for Wigner-type matrices [arXiv:1809.03971,arXiv:1811.04055]. We also extend the finite dimensional band mass formula from [arXiv:1804.07744] to the von Neumann algebra setting by showing that the spectral mass of the bands is topologically rigid under deformations and we conclude that these masses are quantized in some important cases.