Non-Hermitian Rayleigh-Schrödinger-like Perturbation Theory at Exceptional Point

We develop a Rayleigh--Schr\"{o}dinger-like perturbation theory for non-Hermitian quantum systems at an exceptional point of order $N$. Working in the Jordan basis of the unperturbed Hamiltonian and employing a Puiseux expansion of the perturbed eigenvalues and eigenstates, we derive explicit recursion relations for the expansion coefficients. The corrections to the unperturbed eigenvalue in the Puiseux expansion govern the splitting near the exceptional point; the first two are obtained iteratively in two equivalent forms. One is given in terms of the perturbation Hamiltonian in the Jordan basis, and the other in terms of the generator that drives the eigenvalue evolution with respect to the perturbation. The latter constitutes the exceptional-point counterpart of a geometric perturbation method recently developed for the non-exceptional-point regime. Both representations are verified explicitly for the $N = 2$ and $N = 3$ cases.