Optimal calculation of the pair correlation function for an orthorhombic system

We present a new computational method to calculate arbitrary pair correlation functions of an orthorombic system in the most efficient way. The algorithm is demonstrated by the calculation of the radial distribution function of shock compressed liquid hydrogen.