Calculation of elementary diagrams using a Metropolis-like simulation method

Abstract

A method for the evaluation of higher order elementary diagrams is proposed. It utilizes symbolic algebra programming and automatic code generation. The method is used to calculate elementary diagrams for hard spheres up to the fifth. The obtained values are fitted to suitable functions of interparticle distance. The resulting truncated density expansion is compared with recent values of the bridge function calculated directly from the simulated pair distribution function. The method can be straightforwardly extended to hard-body systems and their mixtures.