Graph Symmetry Reveals Hidden Exceptional Points in Open Quantum Systems

Physicists have introduced a symmetry-based method that finds exceptional points—singularities where a quantum system’s eigenvalues and eigenvectors collapse together—directly from the full equations governing realistic open quantum systems, rather than from simplified models engineered to display them. The approach, described in a study published in Discover Physics on 7 October 2026 by Eric Bittner, Bhavay Tyagi and Kevin Bassler, uses the symmetry of correlated dissipation to break a large quantum dynamics problem into small, exactly solvable pieces in which exceptional behavior can be located precisely.

Exceptional points are among the most striking phenomena in non-Hermitian physics, the branch of physics that describes systems which lose energy or information to their surroundings. At an exceptional point, two or more eigenvalues of the system’s governing operator merge, and at the same instant their eigenvectors coalesce. The operator becomes defective: it no longer possesses a complete set of independent eigenvectors. This coalescence produces non-analytic parameter dependence, transient amplification of signals, and abrupt qualitative changes in dynamical behavior, which is why exceptional points have become central organizing structures in modern non-Hermitian physics and have been pursued in photonics, wave systems and sensing technologies.

Most theoretical studies of exceptional points begin from deliberately constructed low-dimensional effective Hamiltonians whose parameters are tuned to produce the desired spectral degeneracies. Real open quantum systems, however, are governed not by Hamiltonians alone but by Liouvillian superoperators, which unify coherent evolution, irreversible dissipation and quantum jumps within a single dynamical generator. For a quantum system whose Hilbert space has dimension d, the Liouvillian is generally a non-Hermitian matrix of size d squared by d squared, acting on the space of density operators. Its eigenvalues encode both decay rates and coherent oscillation frequencies. Exceptional points of the Liouvillian therefore arise as intrinsic properties of the complete open-system dynamics rather than of an assumed effective Hamiltonian, yet a general framework for identifying such exceptional dynamics directly from microscopic dissipative models has been lacking.

The new work supplies that framework by exploiting a structural feature of trace-preserving quantum dynamics. Because the trace of the density matrix is conserved, the identity operator is always a left zero-eigenvector of the Liouvillian, and the operator space splits cleanly into a traceful part and a traceless part. The traceless subspace is invariant under the dynamics, and all nontrivial decay, oscillation, degeneracy and Jordan-block structure relevant to measurable observables is governed by the restriction of the Liouvillian to that sector. Crucially, the authors distinguish rank deficiency, which reflects a nontrivial kernel associated with steady states or conserved quantities, from spectral defectiveness, which arises when the geometric multiplicity of an eigenvalue falls below its algebraic multiplicity. A trace-preserving Liouvillian may possess a unique stable stationary state and still exhibit defective dynamics within nonzero spectral sectors, producing algebraic time dependence and enhanced transient sensitivity.

The key insight concerns correlated dissipation. When the noise acting on a network of quantum sites is correlated, the noise-correlation matrix can be read as the adjacency matrix of a graph, with its elements quantifying the degree of bath correlation between sites. The eigenvectors of this matrix furnish irreducible representations of the dissipation network, and in this basis the dissipator decomposes into independent symmetry sectors. The Liouvillian inherits the symmetry group of the underlying graph, so the full dynamics block-diagonalizes into a direct sum of low-dimensional invariant blocks. Exceptional points then emerge within specific symmetry-selected subspaces rather than as global singularities of the entire spectrum. Observables supported predominantly within a defective sector may show algebraic prefactors or enhanced transient amplification even while the global steady state remains unique and stable, whereas observables in diagonalizable sectors remain insensitive to the non-Hermitian criticality.

To make these ideas concrete, the authors analyzed a two-site dimer coupled to a common reservoir, treating two distinct dissipation mechanisms that share the same noise-correlation matrix. In the correlated relaxation case, where local operators remove excitations, projection onto the symmetric and antisymmetric sector yields a reduced two-by-two generator whose eigenvalues depend on the site-energy imbalance and the correlation strength. An exceptional point occurs when the spectral discriminant vanishes, at a critical detuning proportional to half the absolute value of one minus twice the correlation parameter. In the correlated dephasing case, where the environment randomly modulates local transition frequencies without exchanging populations, the reduced generator instead depends on the tunneling amplitude, and the exceptional point falls at a critical correlation equal to one minus the ratio of tunneling to dephasing rate.

These two conditions reveal a striking inversion. Under dephasing, correlations suppress the effective dissipation in the protected antisymmetric channel and stabilize coherent oscillations, enlarging the region of parameter space in which the dynamics remains PT-preserving; correlations act protectively. Under relaxation, correlations enhance the imbalance between collective decay channels and promote non-Hermitian mixing, expanding the PT-broken region instead. With uncorrelated dissipation, any finite energy imbalance immediately generates a complex-conjugate eigenvalue pair in the symmetry-resolved block, so finite bath correlations are required to stabilize a preserved regime at all. Identical noise correlations therefore generate inequivalent exceptional-point landscapes depending on the physical dissipation mechanism, a result the authors identify as the central finding of their minimal-model analysis. The terminology of PT-preserved and PT-broken regimes follows their earlier work on correlated-noise synchronization, although they neither require nor assume a formally defined parity-time symmetry of the complete Liouvillian.

Beyond analytic reduction, the study introduces a numerical diagnostic called the exceptional-point strength, denoted by the symbol E, defined as the inverse of the smallest singular value of the matrix of right eigenvectors of the Liouvillian. At an exceptional point the eigenvectors become linearly dependent and this matrix loses rank, so the measure diverges; for diagonalizable generators it remains finite. The definition is invariant under similarity transformations and independent of eigenvalue ordering, making it a basis-independent probe of proximity to defective dynamics that requires no analytic reduction of the Liouvillian. Perturbation theory predicts that near a generic second-order exceptional point the eigenvalue splitting and eigenvector coalescence scale as the square root of the distance to the singularity, implying that the exceptional-point strength grows as the inverse square root of that distance.

Applied to tight-binding models with correlated dephasing and relaxation, the diagnostic performed as predicted. In one-dimensional parameter scans of the full Liouvillian, the exceptional-point strength exhibited a sharp divergence precisely at the analytically predicted exceptional-point locations, and on logarithmic axes its growth near the peak followed the expected power law, confirming the square-root coalescence of eigenvectors at second-order exceptional points. Away from the seam the measure remained finite and slowly varying, indicating a well-conditioned eigenbasis. Two-dimensional scans reproduced the full exceptional seams, and the method also revealed universal scaling of the strength near the singularities. The same rank-and-nullity test applied after symmetry reduction allows defectiveness to be localized to specific invariant sectors of any Liouvillian superoperator once its symmetry-resolved blocks are identified.

The framework also clarifies how long-lived collective oscillations arise in dissipative networks. When the correlation strength approaches its physical bound, set by complete positivity of the Lindblad generator, the dissipation eigenvalue of a protected sector vanishes and the corresponding Liouvillian mode becomes strictly marginal. Time-domain propagation then shows the dynamics converging onto a persistent oscillatory trajectory in operator space rather than relaxing to a static fixed point—a limit cycle of the density matrix stabilized by symmetry-induced suppression of dissipation rather than by fine tuning. For connected k-regular graphs, the fully symmetric sector is the only irreducible representation that can become truly protected under anti-correlated noise, while all other adjacency eigenmodes retain strictly positive decay rates even at the positivity boundary. The authors note that the approach inverts the conventional logic of the field: instead of assuming a low-dimensional non-Hermitian reduction and tuning its parameters, the full Lindbladian is treated as the fundamental object and exceptional behavior is detected through numerical conditioning diagnostics. Because the diagnostics rely only on repeated application of the Liouvillian rather than dense diagonalization, they are compatible with tensor-network representations, offering a practical route to hunting exceptional points in many-body open systems such as excitonic aggregates, spin networks, superconducting and photonic lattices, and cavity-QED arrays where correlated noise and network structure are essential.

Subject of Research: Symmetry-resolved identification of exceptional points in open quantum systems governed by Lindblad dynamics

Article Title: Graph symmetry organizes exceptional dynamics in open quantum systems

Article References: Bittner, E. R., Tyagi, B., & Bassler, K. E. (2026). Graph symmetry organizes exceptional dynamics in open quantum systems. Discover Physics, 2(1), Article 14. https://doi.org/10.1007/s44418-026-00018-8

Image Credits: AI Generated

DOI: 10.1007/s44418-026-00018-8

Keywords: exceptional points, open quantum systems, Liouvillian, Lindblad dynamics, graph symmetry, correlated dissipation, non-Hermitian physics, eigenvector coalescence, dephasing, relaxation, limit cycles, tensor networks