Path Integral Monte Carlo calculation of hydrogen Hugoniot

ORAL

Abstract

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.

*Prepared by LLNL under Contract DE-AC52-07NA27344.

–

Authors

  • Saad Khairallah

    • Lawrence Livermore National Laboratory

  • John Shumway

    • Arizona State University

  • Erik Draeger

    • Lawrence Livermore National Laboratory