This submodularity property is known to not hold, as the OP has found out already. However, I'd like to mention the following observation:

Let $f$ be defined on $(0,\infty)$ such that $-f'$ is operator monotone (i.e., for $A\le B \implies f'(A) \ge f'(B)$), then
\begin{equation*}
\text{trace} f(A+B+C)+\text{trace} f(A) \le \text{trace}f(A+B)+\text{trace}(A+C).
\end{equation*}

**Example:** the above inequality holds for $f(t)=t^p$ for $p\in (0,1)$.