Puiseux series about exceptional singularities dictated by symmetry-allowed Hessenberg forms of perturbation matrices
Phys. Rev. A 114, 042212 – Published 16 October, 2026
DOI: https://doi.org/10.1103/3fxh-7y36
Abstract
We develop a systematic framework for determining the nature of exceptional points of order () in non-Hermitian (NH) systems, represented by complex square matrices. By expressing symmetry-preserving perturbations () in the Jordan-normal basis of the defective matrix at an , we show that the upper- Hessenberg structure of the perturbation matrix directly dictates the leading-order eigenvalue- and eigenvector-splitting to be , when expanded in a Puiseux series. Applying this to three-band NH models invariant under parity (P), charge-conjugation (C), or parity-time-reversal (PT), we find that in P- and C-symmetric systems are restricted to at most branch points, while PT-symmetric systems generically support with the strongest possible singularities (viz., ). We illustrate these results with concrete three-dimensional models in which exceptional curves and surfaces emerge. We further show that fine-tuned perturbations can suppress the leading-order branch point to a less-singular splitting, which have implications for designing direction-dependent EP-based sensors. The Appendices extend the analysis to four-band P- and C-symmetric models, establishing the existence of with and singularities, respectively.