Path integral Monte Carlo simulations of warm dense sodium