We read "everywhere" that, because of Lorentz invariance, the energy-momentum tensor of any field in the vacuum state should reduce to a constant multiplying the metric tensor (I'm using the metric signature $\eta = (1, -1, -1, -1)$ and units such that $c = 1$): $$\tag{1} T_{ab}^{(0)} = \rho_{\text{vac}} \, \eta_{ab}. $$ I agree with that. I'm considering the case of the electromagnetic field: $F_{ab} = \partial_a A_b - \partial_b A_a$, of energy-momentum (up to a constant factor) $$\tag{2} T_{ab} = F_a^{\; c} \, F_{cb} + \frac{1}{4} \, \eta_{ab} \, F_{cd} \, F^{cd}. $$ These components are such that their trace is 0, by construction: $\mathrm{Tr} (T_{ab}) \equiv T_a^{\; a} \equiv 0$. For an homogeneous and isotropic state, this tensor reduces to the one representing a perfect fluid: $$\tag{3} T_{ab} = (\rho + p) \, u_a \, u_b - \eta_{ab} \, p. $$ So the condition $\mathrm{Tr} (T_{ab}) = \rho - 3 p = 0$ imposes $p = \frac{1}{3}\, \rho$, as it should for random electromagnetic radiation. Whatever the energy density $\rho \ne 0$ (and $\rho < \infty$), this energy-momentum cannot reduces to (1). The classical solution to the vacuum state is to have $\rho_{\text{vac}} = 0$. My issue appears when I consider the "quantum mechanics vacuum" of the radiation field, according to which $\rho_{\text{vac}} \ne 0$ (the usual "naive" calculations actually give a divergent integral: $\rho_{\text{vac}} \rightarrow \infty$). In that case, how can we reconcile the 0 trace condition of (2)-(3) (which can't reduce to (1)) and the Lorentz invariance of the vacuum (i.e energy-momentum tensor (1)), which implies the vacuum pressure $p_{\text{vac}} = -\, \rho_{\text{vac}}$ instead of $p_{\text{vac}} = \frac{1}{3} \, \rho_{\text{vac}}$?
EDIT: Maybe my query isn't clear enough, so I'll try to be more precise in what I want to know. The vacuum is Lorentz invariant, so we need to get the tensor (1) with some quantum calculations. But then, (2) has a 0 trace, by construction. So it appears to be mathematically impossible to get (1) from (2), except in the special and trivial (non-quantum) case $\rho = 0$. The usual quantum calculation gives (this is the "naive" calculation) $$\tag{4} \rho_{\text{vac}} \propto \int_0^{\infty} \omega^3 \, d\omega = \infty. $$ Authors frequently use an "hard cutoff" at $\omega_{\text{max}}$ to get a finite integral, but this explicitly destroys the Lorentz invariance of the vacuum, so that $p = \frac{1}{3} \, \rho$, instead of the proper relation $p = -\, \rho$. So it is clear that the hard cutoff isn't right. If really $\rho_{\text{vac}}$ is finite and not 0, then the energy-momentum (2) can't be the right expression to be used in the quantum calculation.