Let $B$ be an $n\times n$ matrix, and define $f$ to be the function that maps positive semidefinite (PSD) $n\times n$ matrices $A$ to real numbers by
$$ f(A) = \mathrm{trace}( (B^*A^2B)^{1/3}). $$
In other words, $f$ maps $A$ to the sum of $1/3$-powers of the eigenvalues of the PSD matrix $B^*A^2B$.
Is the function $f$ concave over the PSD cone? I.e. is it true that for any two PSD matrices $X$ and $Y$, $f((X + Y)/2) \ge f(X)/2 + f(Y)/2$?
A more general question is whether the function $f(A) = \mathrm{trace}( (B^*A^2B)^{p})$ is concave over the PSD cone for $0< p < 1/2$.
Carlen and Lieb have some closely related results, but I could not find the particular combination of matrix powers in their paper or in other related work on trace inequalities.