Elastic properties of dense solid phases of hard cyclic pentamers and heptamers in two dimensions

Systems of model planar, nonconvex, hard-body "molecules" of fivefold and sevenfold symmetry axes are studied by constant pressure Monte Carlo simulations with variable shape of the periodic box. The molecules, referred to as pentamers (heptamers), are composed of five (seven) identical hard disks "atoms" with centers forming regular pentagons (heptagons) of sides equal to the disk diameter. The elastic compliances of defect-free solid phases are computed by analysis of strain fluctuations and the reference (equilibrium) state is determined within the same run in which the elastic properties are computed. Results obtained by using pseudorandom number generators based on the idea proposed by Holian and co-workers [Holian et al., Phys. Rev. E 50, 1607 (1994)] are in good agreement with the results generated by DRAND48. It is shown that singular behavior of the elastic constants near close packing is in agreement with the free volume approximation; the coefficients of the leading singularities are estimated. The simulations prove that the highest density structures of heptamers (in which the molecules cannot rotate) are auxetic, i.e., show negative Poisson ratios.