Equations of motion governing the dynamics of the exceptional points of parameterically dependent nonhermitian Hamiltonians

TL;DR

This paper derives equations of motion for the trajectories of exceptional points in parameter-dependent non-Hermitian Hamiltonians, providing a theoretical framework and a practical numerical method to study their response to parameter changes.

Contribution

It introduces a self-contained set of equations of motion for EPs' trajectories, enabling analysis of their dynamics and response to external perturbations in complex quantum systems.

Findings

Derived equations of motion for EP trajectories.

Validated the approach with numerical tests on a toy model.

Provided a practical numerical method for locating EPs.

Abstract

We study exceptional points (EPs) of a nonhermitian Hamiltonian $\hat{H}(\lambda,\delta)$ whose parameters $\lambda \in {\mathbb C}$ and $\delta \in {\mathbb R}$. As the real control parameter $\delta$ is varied, the $k$-th EP (or $k$-th cluster of simultaneously existing EPs) of $\hat{H}(\lambda,\delta)$ moves in the complex plane of $\lambda$ along a continuous trajectory, $\lambda_k(\delta)$. We derive a self contained set of equations of motion (EOM) for the trajectory $\lambda_k(\delta)$, while interpreting $\delta$ as the propagation time. Such EOM become of interest whenever one wishes to study the response of EPs to external perturbations or continuous parametric changes of the pertinent Hamiltonian. This is e.g.~the case of EPs emanating from hermitian curve crossings/degeneracies (which turn into avoided crossings/near-degeneracies when the Hamiltonian parameters are continuously varied). The presented EOM for EPs have not only their theoretical merits, they possess also a substantial practical relevance. Namely, the just presented approach can be regarded even as an efficient numerical method, useful for generating EPs for a broad class of complex quantum systems encountered in atomic, nuclear and condensed matter physics. Performance of such a method is tested here numerically on a simple yet nontrivial toy model.