Free probability theory was developed precisely to deal with noncommuting operators, represented by matrices $X$, $Y$. If $XY\neq YX$, the eigenvalues $\sigma$ of the sum $X+Y$ are not given by the sum of the eigenvalues $\lambda,\mu$ of $X$ and $Y$, so you cannot use the usual construction of the probability distribution $P_\sigma$ of $\sigma$ as the convolution of $P_\mu$ and $P_\lambda$. Free probability theory tells you how to construct, under certain conditions, $P_\sigma$ in terms of $P_\mu$ and $P_\lambda$. If this is the case, $X$ and $Y$ are called "free independent random variables".
Most of the physics applications I know of involve the quantization of classically chaotic systems. Suppose you have a billiard geometry of a certain shape, and you vary the shape a bit, parameterized by some variable $s$. If the billiard has classically chaotic dynamics, then the Hamiltonians $H(s_1)$ and $H(s_2)$ of that system evaluated at two sufficiently different values of $s$ are free independent.