Dots and boxes algorithm for Peierls substitution: application to multidomain topological insulators.

Magnetic fields can be introduced into discrete models of quantum systems by the Peierls substitution. For tight-binding Hamiltonians, the substitution results in a set of (Peierls) phases that are usually calculated from the magnetic vector potential. As the potential is not unique, a convenient gauge can be chosen to fit the geometry and simplify calculations.