By using modular functions on the upper complex half-plane, we study a class of strain energies for crystalline materials whose global invariance originates from the full symmetry group of the underlying lattice. This follows Ericksen’s suggestion which aimed at extending the Landau-type theories to encompass the behavior of crystals undergoing structural phase transformation, with twinning, microstructure formation, and possibly associated plasticity effects. Here we investigate such Ericksen-Landau strain energies for the modelling of reconstructive transformations, focusing on the prototypical case of the square-hexagonal phase change in 2D crystals. We study the bifurcation and valley-floor network of these potentials, and use one in th...