This is not a complete answer but a summary of progress I have made. I suspect the primary author was aware of results related to those I am about to present and decided to include only a special case, perhaps in the expectation that resolving that case would help advance the full conjecture. I now state a conjecture and provide further details and intuition.
Let the power-sum symmetric polynomial of degree $k$ in $n$ variables $x_1,\ldots,x_n$ be denoted by $P_k$, defined by
$$
P_k(x_1,\ldots,x_n)=\sum_{i=1}^n x_i^k.
$$
Below I state a conjecture which I believe may hold under additional assumptions. I remain open to the possibility that the conjecture, as stated, may be false; however, I expect that suitable extra hypotheses will make it valid.
Conjecture 1. Let $(A_i)_{i=1}^{m}$ be a collection of $n \times n$ positive definite Hermitian matrices such that
$\sum_{i=1}^{m} A_i = I_n$. For each matrix $A_i$, let $\lambda(A_i) = (\lambda_1(A_i), \ldots, \lambda_n(A_i))$ denote the corresponding eigenvalues of the matrix, ordered decreasingly.
For $ \alpha_i,\beta_i \in \mathbb{R} $,let
\begin{align*}
f(A_1,\ldots,A_m) = \sum_{i=1}^{m} \alpha_i P_i(A_1,\ldots,A_m) \\
g(A_1,\ldots,A_m) = \sum_{i=1}^{m} \beta_i P_i(A_1,\ldots,A_m)
\end{align*}
be symmetric $m$ variable matrix functions satisfying the scalar inequality, meaning that for any $x_1,\ldots,x_m \in \mathbb{R}^{+}$ such $\sum_{i=1}^{m} x_i = 1$ we have
$$
0 \le g(x_1,\ldots,x_m) \le f(x_1,\ldots,x_m).
$$
Then
$$
\lambda\left( g(A_1, \ldots, A_m) \right) \prec_{w} \lambda\left( f(A_1, \ldots, A_m) \right).
$$
By $\prec_{w}$ we mean weakly majorization.
If the conjecture above is valid, then a suitable modification of the following result yields a stronger statement.
Theorem 1. Let $n\ge 1$, let $d\ge 0$ be even, and let $x_1,\dots,x_n$ be real numbers. If $(x_1,\dots,x_n)$ majorizes $(y_1,\dots,y_n)$, then
$$
h_d(x_1,\dots,x_n)\ge h_d(y_1,\dots,y_n),
$$
with equality if and only if $(y_1,\dots,y_n)$ is a permutation of $(x_1,\dots,x_n)$.
Moreover, if one restricts the variables to be positive and replaces majorization by weak majorization, the same argument leads to stronger matrix inequalities (in particular, one obtains corresponding inequalities for determinants and other spectral quantities).
To support these claims I have tested the conjecture numerically for the original question and have not found counterexamples so far. Another example I examined is the following
$$
f(A,B,C)=\frac{116}{81}I_n + 64\big(A^5+B^5+C^5\big),\qquad
g(A,B,C)=20\big(A^3+B^3+C^3\big),
$$
under the usual assumptions. For these choices, $f-g$ need not be positive semidefinite in general; an explicit $2\times2$ counterexample is given by
$$
A=\begin{bmatrix}0.363 & 0\\[4pt]0 & 0.767\end{bmatrix},\quad
B=\begin{bmatrix}0.306 & -0.13534244755789848\\[4pt]-0.13534244755789848 & 0.075\end{bmatrix},\quad
C=\begin{bmatrix}0.331 & 0.13534244755789848\\[4pt]0.13534244755789848 & 0.158\end{bmatrix},
$$
for which $f-g$ is not necessarily positive semidefinite.
One natural tool to attempt here is the following variational characterization of partial sums of eigenvalues.
Theorem 2 (Ky-Fan).
Let $A$ be an $n\times n$ Hermitian matrix with eigenvalues $\lambda_1(A)\ge\lambda_2(A)\ge\cdots\ge\lambda_n(A)$. Then, for each $1\le k\le n$,
$$
\max_{Y\in\mathbb{C}^{n\times k},\;Y^*Y=I_k}\operatorname{Tr}(Y^* A Y)=\sum_{i=1}^k \lambda_i(A),
$$
where the maximum is taken over all $n\times k$ matrices $Y$ whose columns are orthonormal.
The advantage of this theorem is that, for a given choice of $Y$, the compressed matrices $(Y^*A_iY)_{i=1}^m$ again sum to $I_k$, so one can reduce dimension and attempt an inductive or variational approach. I was, however, unable to make substantial progress along this line. One important point to note is that for any minimizer or maximizer $Y$, its columns span an invariant subspace of $A$ associated with the first $k$ eigenvalues, check On Generalizing Trace Minimization.
For the final part, I will show that the following theorem is correct. It represents a very special case of the material discussed above.
Theorem 3. By the stated assuptions one have
$$
\text{Tr}(5(A^2+B^2+C^2)) \le \text{Tr}(I_n+6(A^3+B^3+C^3))
$$
Proof . Since $ \text{Tr}(A)+\text{Tr}(B)+\text{Tr}(C)=\text{Tr}(I_n) $, a matrix exists such that it's trace is strictly less than $\text{Tr}( I_n/3 )$. Otherwise
$$
\text{Tr}(I_n) = \text{Tr}(\frac{I_n}{3})+\text{Tr}(\frac{I_n}{3})+\text{Tr}(\frac{I_n}{3}) < \text{Tr}(A)+\text{Tr}(B)+\text{Tr}(C) .
$$
Which is not possible. Hence, without loss of generality we can assume that $\text{Tr}(C) \le \text{Tr}(\frac{I_n}{3})$. Additionally consider the Eigendecomposition of $C$ as follow
$$
C=U^{*}C^{'}U,U^{*}U=I_n.
$$
Where $C^{'}$ is a diagoanl matrix. Set
$$
X_n = I_n+6(A^3+B^3+C^3), Y_n = 5(A^2+B^2+C^2) .
$$
It's straightforward to verify that
$$
\text{Tr}(X_n)=\text{Tr}(U^{*}X_nU),\text{Tr}(Y_n)=\text{Tr}(U^{*}Y_nU) .
$$
We additionally establish
$$
U^{*}AU+U^{*}BU+C{'} = U^{*}(A+B+C)U = U^{*}U = I_n .
$$
As a result, it's enough to prove that
\begin{align*}
0 \le \text{Tr}(U^{*}X_nU) - \text{Tr}(U^{*}Y_nU) .
\end{align*}
Let us define
\begin{align*}
A^{'} &= \text{diag}(U^{*}AU) \\
B^{'} &= \text{diag}(U^{*}BU) \\
K &= U^{*}AU-A^{'} .
\end{align*}
Clearly
$$
(A^{'}+K) + (B^{'}-K) + C^{'} = I_n .
$$
From the preceding results, we obtain
\begin{align*}
& \text{Tr}(U^{*}X_nU)-\text{Tr}(U^{*}Y_nU) \\
& = \text{Tr}(I_n+6((U^{*}AU)^3+(U^{*}BU)^3+C'^3)-5((U^{*}AU)^2+(U^{*}BU)^2+C'^2)) \\
& = \text{Tr}(I_n+6((A^{'}+K)^3+(B^{'}-K)^3+C'^3) -5((A^{'}+K)^2+(B^{'}-K)^2+C'^2)) \\
& = \text{Tr}(I_n+6(A'^3+B'^3+C'^3)-5(A'^2+B'^2+C'^2))+6\text{Tr}(K^3-K^3)+6\text{Tr}(3(A'^2-B'^2)K) \\
& + 6\text{Tr}(3(A'+B')K^2)-10\text{Tr}(K^2)-10\text{Tr}((A'-B')K) .
\end{align*}
Since $K$ is a Hermitian matrix and it's diagonal elements are zero and matrices $A^{'}$ and $B^{'}$ are diagonal, we can deduce
$$
6\text{Tr}(3(A'^2-B'^2)K)=10\text{Tr}((A'-B')K)=6\text{Tr}(K^3-K^3)=0 .
$$
Consequently
\begin{align*}
& \text{Tr}(U^{*}X_nU)-\text{Tr}(U^{*}Y_nU) \\
& =\text{Tr}(I_n+6(A'^3+B'^3+C'^3)-5(A'^2+B'^2+C'^2))+\text{Tr}((18(A'+B')-10I_n)K^2).
\end{align*}
The next step is to show that
$$
0 \le \text{Tr}((18(A'+B')-10I_n)K^2) .
$$
To prove this, please note that $A^{'}+B^{'} = I_n- C^{'}$ and $\text{Tr}(C^{'}) \le \text{Tr}(I_n/3)$. Resulting in
\begin{align*}
\text{Tr}((18(A'+B')-10I_n)K^2) & \ge \\
\text{Tr}((18(I_n-C')-10I_n)K^2) & \ge\\
\text{Tr}((18\times \frac{2}{3}I_n-10I_n)K^2) & \ge 0 .
\end{align*}
It's left to demonstrate that
$$
\text{Tr}(I_n+6(A'^3+B'^3+C'^3)-5(A'^2+B'^2+C'^2)) \ge 0 .
$$
Which is equivalent to the scalar case, this completes the proof.
$$\Delta \stackrel{def}{=} 6A\left(A-\frac{5}{12} I_n\right)^2 + 6B\left(B-\frac{5}{12} I_n\right)^2 + 6C\left(C-\frac{5}{12} I_n\right)^2$$
is positive semi-definite. This leads to
$$\det\left(6(A^3+B^3+C^3)+\color{red}{\frac{25}{24}}I_{n}\right)\ge 5^n\det(A^2+B^2+C^2)$$
Since $A + B + C = I_n$, the $3$ squares in $\Delta$ cannot vanish at the same time. It sounds pausible to me $\Delta$ is actually "bounded below" by $\frac{1}{24}I_n$. If this is indeed the case, then the original inequality is true.
– achille hui Aug 21 '14 at 19:19