Cutting and gluing with running couplings in $\mathcal{N}=2$ QCD

We consider the order parameter $u=\left<{\rm Tr}\phi^2\right>$ as function of the running coupling constant $\tau \in \mathbb{H}$ of asymptotically free $\mathcal{N}=2$ QCD with gauge group $SU(2)$ and $N_f\leq 3$ massive hypermultiplets. If the domain for $\tau$ is restricted to an appropriate fundamental domain $\mathcal{F}_{N_f}$, the function $u$ is one-to-one. We demonstrate that these domains consist of six or less images of an ${\rm SL}(2,\mathbb{Z})$ keyhole fundamental domain, with appropriate identifications of the boundaries. For special choices of the masses, $u$ does not give rise to branch points and cuts, such that $u$ is a modular function for a congruence subgroup $\Gamma$ of ${\rm SL}(2,\mathbb{Z})$ and the fundamental domain is $\Gamma\backslash\mathbb{H}$. For generic masses, however, branch points and cuts are present, and subsets of $\mathcal{F}_{N_f}$ are being cut and glued upon varying the mass. We study this mechanism for various phenomena, such as decoupling of hypermultiplets, merging of local singularities, as well as merging of non-local singularities which give rise to superconformal Argyres-Douglas theories.