Systematic search for 𝒫𝒯-symmetric potentials with real energy spectra

Changes of coordinates represent one of the most effective ways of deriving solvable potentials from ordinary differential equations for separate special functions. Here we relax the standard Hermiticity requirement and find an innovative construction which leads to unusual, complex potentials. Their energy spectrum is shown to stay real after a weakening of the Hermiticity of the Schrödinger equation to its mere invariance under the combined (parity) and (time-reversal) symmetry. This ultimately results in richer bound-state spectra. Some of our new exactly solvable potentials generalize the current textbook models. Details are given for constructions based on the hypergeometric and confluent hypergeometric special functions.