Abstract
In this paper, we conducted analysis of p-V-T relation of cryogenic hydrogen using classical Molecular Dynamics (MD) and path integral Centroid MD (CMD) method to understand an effect of quantum nature of hydrogen molecules. We performed NVE constant MD simulation across a wide density-temperature region to obtain an Equation Of State (EOS). Simulation results were compared with experimental data. As a result, it was confirmed that classical MD cannot reproduce the experimental data at the high density region. On the other hand, CMD well reproduces the thermodynamic properties of liquid hydrogen. Moreover, it was clarified that taking the quantum effect into account makes repulsion force larger and the potential well smaller. Because of this mechanism, the intermolecular interaction of hydrogen diminishes and the virial pressure increases.
