Semiclassical expansions up to $\hbar^4$-order in relativistic nuclear   physics
  We present the first calculation of the $\hbar^4$-Wigner--Kirkwood corrections to a relativistic system of fermions in the presence of external scalar and vector potentials. The method we propose allows to compute efficiently semiclassical corrections to one body operators such as mean energies and local densities. It also preserves gauge invariance and produces explicitly convergent results despite some apparent divergencies at the classical turning points. As a byproduct we obtain the $\hbar^4$ corrections stemming from the polarization of the Dirac sea. We apply our results to the relativistic $\sigma$-$\omega$ Lagrangian in the Hartree valence approximation. We compare the semiclassical expansion with the exact result and with the Strutinsky average whenever it can be obtained. We find that the $\hbar^4 $ corrections are much smaller than typical shell effects. Our results provide convincing arguments to neglect higher than second order $\hbar$ effects in the Wigner--Kirkwood scheme to relativistic nuclear physics in the mean field approximation. 
J. Caro, E. Ruiz Arriola and L.L. Salcedo