Qsymm: Algorithmic symmetry finding and symmetric Hamiltonian generation

Symmetry is a guiding principle in physics that allows us to generalize conclusions between many physical systems. In the ongoing search for new topological phases of matter, symmetry plays a crucial role by protecting topological phases. We address two converse questions relevant to the symmetry classification of systems: is it possible to generate all possible single-body Hamiltonians compatible with a given symmetry group? Is it possible to find all the symmetries of a given family of Hamiltonians? We present numerically stable, deterministic polynomial time algorithms to solve both of these problems. Our treatment extends to all continuous or discrete symmetries of non-interacting lattice or continuum Hamiltonians. We implement the algorithms in the Qsymm Python package, and demonstrate their usefulness through applications in active research areas of condensed matter physics, including Majorana wires and Kekule graphene.