Exactly solvable $\mathcal{PT}$-symmetric models in two dimensions

Non-Hermitian, -symmetric Hamiltonians, experimentally realized in optical systems, accurately model the properties of open, bosonic systems with balanced, spatially separated gain and loss. We present a family of exactly solvable, two-dimensional, potentials for a non-relativistic particle confined in a circular geometry. We show that the -symmetry threshold can be tuned by introducing a second gain-loss potential or its Hermitian counterpart. Our results explicitly demonstrate that breaking in two dimensions has a rich phase diagram, with multiple re-entrant -symmetric phases.