Time-independent scattering methods are widely used to analyze transport in periodic $\mathcal{PT}$-symmetric systems. However, their predictions become unphysical when the system supports time-growing bound states (TGBSs), which manifest as $S$-matrix poles in the first quadrant of the complex wave-number plane. Here, we analytically delineate the region of physical relevance for a $\mathcal{PT}$-symmetric chain of $N$ unit cells with gain/loss strength $γ$. We derive the TGBS onset threshold $γ_c = 2\sin[π/(4N)]$, which scales as $π/(2N)$ for large $N$ and vanishes in the thermodynamic limit. Enlarging the structure thus enriches stationary scattering phenomenology but inevitably triggers TGBSs at weaker gain/loss. Time-dependent wave-packet simulations confirm this analytical boundary quantitatively. Applying this criterion, we show that many previously reported predictions of gain-loss-induced localization, reflectionless transport, and coherent perfect absorbers and lasers in large periodic structures fall outside the physically relevant regime. $S$-matrix pole analysis is therefore an indispensable prerequisite for interpreting time-independent scattering predictions in periodic non-Hermitian systems.
Physical relevance of time-independent scattering predictions in periodic PT-symmetric chains
Time-independent scattering methods are widely used to analyze transport in periodic PT-symmetric systems. However, their predictions become unphysical when the system supports time-growing bound states (TGBSs), which manifest as S-matrix poles in the first quadrant of the complex wave-number plane. Here, we analytically delineate the region of physical relevance for a PT-symmetric chain of N unit cells with gain/loss strength γ. We derive the TGBS onset threshold γc=2sin[π/(4N)], which scales as π/(2N) for large N and vanishes in the thermodynamic limit. Enlarging the structure thus enriches stationary scattering phenomenology but inevitably triggers TGBSs at weaker gain/loss. Time-dependent wave-packet simulations confirm this analytical boundary quantitatively. Applying this criterion, we show that many previously reported predictions of gain-loss-induced localization, reflectionless transport, and coherent perfect absorbers and lasers in large periodic structures fall outside the physically relevant regime. S-matrix pole analysis is therefore an indispensable prerequisite for interpreting time-independent scattering predictions in periodic non-Hermitian systems.
Thresholdless dynamical instability in transverse PT-symmetric scattering systems: The hidden role of bound states in the continuum
The stationary scattering properties of transverse parity-time (PT) symmetric systems have been extensively studied, yet their dynamical stability, a prerequisite for any stationary description, remains largely unexplored. Here we uncover a thresholdless dynamical instability in such systems, driven by symmetry-protected bound states in the continuum (BICs). Without gain and loss, the up-down mirror symmetry of the transverse geometry generically protects a BIC, which manifests as the coalescence of an S-matrix pole and zero on the real axis of the complex wave-number plane. A Hermitian symmetry-breaking perturbation shifts the poles into the lower half-plane, converting the BIC into a resonance with a Fermi-golden-rule decay width. An anti-Hermitian PT-symmetric perturbation instead reverses the sign of the second-order energy shift, driving the poles into the upper half-plane and producing a time-growing bound state with ImE=γ2ΓV/2+O(γ3), where γ is the gain-loss strength and ΓV is the golden-rule coupling of the BIC to the continuum. The instability therefore sets in at arbitrarily small γ whenever ΓV>0. We confirm this mechanism in a three-site side-coupled model that is unstable despite possessing a unitary scattering matrix, and in a four-site rhombic model where transverse and longitudinal gain-loss placements yield vanishing and finite thresholds, respectively. These results establish a microscopic stability criterion and a design principle for stable PT-symmetric scattering devices.
Time-independent scattering methods are widely employed to analyze transport in non-Hermitian systems. Their application, however, rests on a critical yet often overlooked assumption: that an incident wave is a pure superposition of scattering states. In practice, any physically realistic, spatially localized wave packet will generally have a nonzero overlap with the system's bound states, thereby violating this premise. While this violation is inconsequential in Hermitian systems, it can invalidate the conventional scattering picture in their non-Hermitian counterparts. The underlying cause is the emergence of time-growing bound states, which manifest as poles of the scattering matrix (S matrix) in the first quadrant of the complex wave-number plane. Any initial overlap with these states becomes exponentially amplified, eventually dominating the long-time dynamics. Consequently, the actual evolution of a wave packet diverges dramatically from the conventional scattering picture, rendering the transmission and reflection coefficients derived from time-independent scattering methods unphysical. Using tight-binding models with non-Hermiticity introduced via imaginary on-site potentials or asymmetric hopping, we demonstrate that parameter regimes supporting such growing states are common. We therefore conclude that an analysis of S-matrix poles is an indispensable step to confirm the physical relevance of time-independent scattering calculations in non-Hermitian systems.
Parity-time (PT) symmetric systems are classical, gain-loss systems whose dynamics are governed by non-Hermitian Hamiltonians with exceptional-point (EP) degeneracies. The eigenvalues of a PT-symmetric Hamiltonian change from real to complex conjugates at a critical value of gain-loss strength that is called the PT breaking threshold. Here, we obtain the PT threshold for a one-dimensional, finite Kitaev chain---a prototype for a p-wave superconductor---in the presence of a single pair of gain and loss potentials as a function of the superconducting order parameter, on-site potential, and the distance between the gain and loss sites. In addition to a robust, nonlocal threshold, we find a rich phase diagram for the threshold that can be qualitatively understood in terms of the band structure of the Hermitian Kitaev model. In particular, for an even chain with zero on-site potential, we find a re-entrant PT-symmetric phase bounded by second-order EP contours. Our numerical results are supplemented by analytical calculations for small system sizes.
In non-Hermitian systems, spectra can be maximally boundary-sensitive, yet bulk physics need not be. Here we experimentally show that spectral moments provide boundary-robust bulk observables in finite non-Hermitian lattices, even when the spectra undergo dramatic geometry-dependent reshaping due to the skin effect. Using a unified acoustic platform with full spectral reconstruction and time-domain access, we probe one-, two- and three-dimensional lattices and demonstrate that spectral moments remain nearly invariant across distinct boundary geometries while the corresponding complex spectra differ strongly. To connect the thermodynamic theorem to realistic finite systems, we develop a loop-counting theory that identifies the physical origin of finite-size deviations in terms of missing boundary loops, quantitatively captures the corrections, and predicts a scaling law, which we verify experimentally. Beyond acoustic spectroscopy, we reveal a counterintuitive dynamical consequence of moment invariance: a dispersive-to-proliferative bulk transition governed by bulk moment structure rather than spectral boundary sensitivity. As a result, local bulk dynamics can remain stable (dispersive) even in a PT-broken spectral regime, challenging the conventional expectation that PT breaking necessarily implies feedback-induced dynamical instability (proliferation) through exponentially amplifying spectral components. These results establish spectral moments as practical bulk descriptors for finite non-Hermitian matter and open a route to extracting and controlling intrinsic bulk behavior in realistic wave-based non-Hermitian devices.
Analysis of fundamentally open systems which exhibit both spatial reflection (parity) and time-reversal symmetry, called PT-symmetric [1] systems, has played an important role in understanding the dynamics of a broad class of non-Hermitian systems. Such systems support distinct regimes of stability and instability respectively, and the transition between these plays an important role in the system dynamics. We investigate a variety of non-Hermitian systems which exhibit a PT-symmetry breaking transition in both classical and quantum contexts using an array of theoretical and numerical approaches which are directly amenable to physical experiments. Organizationally, we divide these studies into those done on classical and quantum systems. Specifically for classical systems we analyze several interesting phenomena in the context of optical waveguide arrays. These provide an important testing ground for much of the discrete, PT-symmetric theory. In our case, we analyze PT-symmetric extensions of several interesting lattice models. Firstly, we explore single gain/loss pair, PT-symmetric extensions of the Su-Schreiffer-Heeger (SSH) model [2] and the Aubry-Andre-Harper (AAH) model [3,4]. In both cases, we find a veiled symmetry which is responsible for the PT-symmetry. Secondly we expand upon the Rudner-Levitov (RL) model [5,6], a non-Hermitian extension of the SSH model with multiple gain/loss pairs, and we explore the dependence of topological phenomena on changes in the initial state as well as their robustness in the presence of gain-induced nonlinearity. For quantum systems, we explore the rich PT-symmetric phase diagram resulting from time-periodic, Floquet PT-symmetry [7,8]. This is experimentally implementable in ultracold atomic systems [9], and features a reentrant PT-symmetric phase. Secondly, we examine quantum multi-photon statistics for a simple non-Hermitian model which is transformed and mapped to a system with unitary evolution via a unitary dilation, and is fully realizable on a small-scale quantum simulator [10].
A PT-Symmetric XXZ Chain with an External Constraint Field: Exceptional Points, Finite-Size Scaling, and Falsifiable Predictions This work introduces a non-Hermitian quantum-spin model that combines three established areas of modern theoretical physics: The XXZ Heisenberg spin-½ chain; PT-symmetric gain/loss dynamics; An externally coupled constraint field capable of modifying the location of exceptional points. The resulting Hamiltonian provides a mathematically explicit framework for studying how externally driven constraints interact with open quantum systems exhibiting parity-time symmetry. Rather than proposing a speculative unified theory, the objective of this work is to formulate a concrete, experimentally testable model that produces quantitative predictions capable of confirmation or falsification through laboratory measurement. Scientific Background Parity-time (PT) symmetric systems occupy an important position within contemporary non-Hermitian physics. Although such systems are not Hermitian, they can nevertheless exhibit purely real energy spectra within specific parameter regimes. As gain-loss asymmetry increases, the system eventually reaches a critical transition known as an Exceptional Point (EP). At an Exceptional Point: eigenvalues coalesce; eigenvectors coalesce; sensitivity to perturbations increases dramatically; system behavior changes qualitatively. Exceptional points have become a major area of research in: photonics; optics; condensed matter physics; open quantum systems; topological materials; waveguide lattices; sensing technologies. The present work extends this framework by introducing an explicit constraint-field coupling and examining its impact on exceptional-point behavior. Model Overview The total Hamiltonian consists of three components: Hermitian XXZ Backbone The first contribution is the XXZ spin-chain Hamiltonian. This term represents: coherent spin transport; exchange interactions; anisotropic coupling between neighboring spins. In the absence of gain/loss and external constraints, the model reduces to the conventional XXZ chain studied throughout quantum many-body physics. PT-Symmetric Gain/Loss Sector The second contribution introduces alternating gain and loss across the lattice. This term produces: non-Hermitian dynamics; PT-symmetric behavior; exceptional-point transitions. The gain/loss amplitude controls the transition between: PT-unbroken phases; PT-broken phases. External Constraint Field The third contribution is a spatially structured driving field. This field is modeled as a superposition of harmonically varying modes possessing: amplitude; frequency; phase; spatial localization. The field couples directly to the spin degrees of freedom and modifies the effective dynamics of the system. Exact Two-Site Solution To establish a mathematically transparent foundation, the paper solves the two-site case exactly. The resulting spectrum yields a closed-form exceptional-point condition: γEP=J\gamma_{EP}=J where: JJ is the exchange coupling; γ\gamma is the gain/loss strength. This analytical result separates three physically distinct regimes: PT-Unbroken Phase γ<J\gamma < J Properties: real eigenvalues; coherent transport; stable dynamics. Exceptional Point γ=J\gamma = J Properties: eigenvalue coalescence; eigenvector coalescence; maximal sensitivity. PT-Broken Phase γ>J\gamma > J Properties: complex eigenvalues; dissipative dominance; broken PT symmetry. Central Prediction The principal result of the paper is the prediction that the exceptional point becomes tunable when the external constraint field is activated. Specifically: γEP(λ)=J1+λ⟨G⟩J\gamma_{EP}(\lambda) = J\sqrt{ 1+ \frac{\lambda\langle G\rangle}{J} } where: λ\lambda is the coupling strength; ⟨G⟩\langle G\rangle is the mean constraint-field contribution. For weak coupling: ΔγEP≈λ⟨G⟩2\Delta\gamma_{EP} \approx \frac{\lambda\langle G\rangle}{2} The prediction is therefore simple: increasing constraint coupling should systematically shift the exceptional point. Falsifiability A central design goal of this work is scientific falsifiability. The model can be rejected by experiment. The theory is considered falsified if: no exceptional point is observed; the exceptional point remains fixed when the constraint field is varied; the observed shift fails to scale with the predicted dependence; an alternative model explains the data more successfully. These conditions are intentionally defined prior to experimentation. Finite-Size Scaling Program While the exact solution is obtained for the two-site system, the broader objective is to determine whether the phenomenon survives in larger many-body systems. The proposed finite-size scaling program investigates: N=4,6,8,10,12N = 4,6,8,10,12 site chains through exact diagonalization. The goal is to determine whether the exceptional-point shift converges toward a well-defined thermodynamic limit. A successful scaling analysis would provide evidence that the mechanism is not merely a finite-size artifact. Feedback Dynamics The paper further proposes a feedback mechanism in which the constraint field is partially sourced by the state of the spin chain itself. This introduces: nonlinear behavior; effective gain/loss renormalization; hysteresis effects; potentially emergent collective behavior. A predicted observable is a measurable hysteresis loop in the exceptional-point response. Experimental Realization A practical realization is proposed using PT-symmetric photonic waveguide arrays. In such systems: waveguides represent lattice sites; optical coupling represents exchange interactions; alternating amplification and attenuation implement gain/loss; an externally driven resonator grid provides the constraint field. The experimentally measurable quantity is the eigenvalue splitting: ΔE=4J2−γ2\Delta E = 4\sqrt{J^2-\gamma^2} which collapses at the exceptional point. By measuring how the collapse location changes with constraint strength, the core prediction can be directly tested. Research Roadmap The project is organized into successive stages: Stage 1 Exact analytical treatment of the two-site model. Stage 2 Finite-size numerical diagonalization. Stage 3 Photonic-lattice simulation. Stage 4 Laboratory implementation. Stage 5 Independent experimental replication. Contribution This work should be understood as a contribution to: non-Hermitian quantum physics; PT-symmetric systems; exceptional-point dynamics; open quantum systems. The paper introduces: an explicitly defined Hamiltonian; a solvable toy model; a finite-size scaling program; a laboratory realization pathway; a quantitative prediction subject to experimental verification. Its scientific value rests not on the novelty of the concept alone, but on the existence of a clear pathway by which the prediction may be tested, reproduced, supported, revised, or rejected. Keywords: PT Symmetry, Exceptional Points, XXZ Spin Chain, Non-Hermitian Physics, Open Quantum Systems, Quantum Many-Body Systems, Finite-Size Scaling, Photonic Lattices, Gain-Loss Dynamics, Condensed Matter Physics, Mathematical Physics, Quantum Information, Experimental Quantum Optics, Falsifiable Prediction, Project Prometheus.