Path Integral Monte Carlo calculation of hydrogen Hugoniot

We use restricted path integral Monte Carlo to calculate the compression of liquid Hydrogen/Deuterium between 50 and 600 GPa for a temperature ranging from 10,000 to 1 Million Kelvin. Different theories and experiments place the compression between 4 and 6. We attempt to tackle this inconsistency with a new fermion fixed node algorithm and careful finite size effect study. Our preliminary results place the compression at 4.3.