Entropy and chirality in sphinx tilings

As a toy model of chiral interactions in crowded spaces, we consider sphinx tilings in finite regions of the triangular lattice. The sphinx tiles, hexiamonds composed of six equilateral triangles in the shape of a stylized sphinx, come in left and right enantiomorphs. Regions scaled up from the unit sphinx by an integer factor ("Sphinx frames") require tiles of both chiral forms to produce tilings, including crystalline, quasicrystalline, and fully disordered tilings. For frames up to order 13, we describe methods that permit exact enumeration and computation of partition functions using "accelerated backtracking," "seam," and "dangler" algorithms. For larger frames, we introduce a Monte Carlo (MC) method to sample typical tilings. Key to the latter is the identification of fundamental shapes (polyads) that admit multiple tilings and which allow a rejection-free MC simulation.