  Let $(\kappa_n(a))_{n\geq 1}$ denote the sequence of free cumulants of a
random variable $a$ in a non-commutative probability space
$(\mathcal{A},\varphi)$. Based on some considerations on bipartite graphs, we
provide a formula to compute the cumulants $(\kappa_n(ab+ba))_{n\geq 1}$ in
terms of $(\kappa_n(a))_{n\geq 1}$ and $(\kappa_n(b))_{n\geq 1}$, where $a$ and
$b$ are freely independent. Our formula expresses the $n$-th free cumulant of
$ab+ba$ as a sum indexed by partitions in the set $\mathcal{Y}_{2n}$ of
non-crossing partitions of the form
  \[ \sigma=\{B_1,B_3,\dots, B_{2n-1},E_1,\dots,E_r\}, \quad \text{with }r\geq
0, \]
  such that $i\in B_{i}$ for $i=1,3,\dots,2n-1$ and $|E_j|$ even for $j\leq r$.
Therefore, by studying the sets $\mathcal{Y}_{2n}$ we obtain new results
regarding the distribution of $ab+ba$. For instance, the size
$|\mathcal{Y}_{2n}|$ is closely related to the case when $a,b$ are free Poisson
random variables of parameter 1. Our formula can also be expressed in terms of
cacti graphs. This graph theoretic approach suggests a natural generalization
that allows us to study quadratic forms in $k$ free random variables.
