arxiv: 2605.00117 · v1 · submitted 2026-04-30 · 🪐 quant-ph · math-ph· math.MP

Recognition: unknown

Topological Charge of Causality at a PT-Symmetric Exceptional Point

Kejun Liu

Authors on Pith no claims yet

Pith reviewed 2026-05-09 20:18 UTC · model grok-4.3

classification 🪐 quant-ph math-phmath.MP

keywords PT symmetryexceptional pointcausalityKramers-Kronig relationstopological chargereflection coefficientBlaschke winding numberopen quantum systems

0 comments

The pith

In a PT-symmetric open dimer, causality carries a topological charge once the exceptional point is crossed, leaving a Lorentzian residual in Kramers-Kronig relations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper argues that causality is not a simple binary property of whether a response function is analytic in the upper half-plane. In a PT-symmetric open dimer, crossing the exceptional point by tuning the gain-loss parameter causes one pole of the reflection coefficient to migrate into the upper half-plane. This migration jumps the Blaschke winding number from zero to one and adds a fixed Lorentzian residual to any attempt at standard Kramers-Kronig reconstruction. The size of the violation is largest exactly at the threshold and falls off as the distance to the critical value raised to a power near -1.08. A reader should care because the effect is sharp, protected by the structure of the exceptional point, and directly testable in a one-port reflection measurement.

Core claim

We show that in a PT-symmetric open dimer it instead carries a topological charge. As the gain-loss parameter crosses the exceptional point, a single pole of the reflection coefficient migrates into the upper half-plane, the Blaschke winding number jumps from 0 to 1, and standard Kramers-Kronig reconstruction acquires a Lorentzian residual fixed by the pole residue. The violation magnitude scales as Delta_KK ~ |gamma - gamma_c|^nu with nu ~ -1.08 in the single-port geometry. The transition is sharp, protected by the codimension-one structure of the exceptional point, and directly measurable in a one-port reflection experiment.

What carries the argument

Migration of a single pole of the reflection coefficient across the real axis at the exceptional point, which changes the Blaschke winding number from 0 to 1 and adds a residue-determined Lorentzian term to the dispersion relation.

If this is right

The breakdown of standard Kramers-Kronig relations is sharp and occurs exactly when the exceptional point is crossed.
Reconstruction of the response must include an extra Lorentzian term whose amplitude is set by the residue of the migrated pole.
The violation reaches its maximum strength at the exceptional point and weakens as the system enters deeper into the broken-symmetry phase.
The residual itself can be extracted in a THz time-domain spectroscopy experiment on a single-port device.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same pole-migration mechanism could assign topological charge to causality in other non-Hermitian scattering setups that host exceptional points.
Detection of the residual offers a concrete experimental signature for locating the PT-symmetry breaking threshold without needing full eigenvalue spectroscopy.
The observed scaling exponent near -1.08 may link to other critical exponents that appear when exceptional points control open-system dynamics.
Dispersion engineering in PT-symmetric devices could deliberately exploit the corrected relation to control group delay or absorption lineshapes.

Load-bearing premise

The reflection coefficient of the chosen PT-symmetric open dimer has exactly one pole that migrates across the real axis precisely when the exceptional point is reached.

What would settle it

Perform a one-port reflection measurement on a PT-symmetric dimer while sweeping the gain-loss parameter through the exceptional point and check whether a Lorentzian residual appears in the Kramers-Kronig reconstruction exactly at the crossing, with its magnitude scaling as the distance to threshold to the power near -1.08.

Figures

Figures reproduced from arXiv: 2605.00117 by Kejun Liu.

**Figure 2.** Figure 2: FIG. 2. Phase diagram in the ( [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. KK residual at ( [PITH_FULL_IMAGE:figures/full_fig_p002_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Log-log scaling of [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗

read the original abstract

Causality in linear response is conventionally treated as a binary property: a response function is either analytic in the upper half-plane or it is not. We show that in a PT-symmetric open dimer it instead carries a topological charge. As the gain-loss parameter crosses the exceptional point, a single pole of the reflection coefficient migrates into the upper half-plane, the Blaschke winding number jumps from 0 to 1, and standard Kramers-Kronig (KK) reconstruction acquires a Lorentzian residual fixed by the pole residue. The transition is sharp, protected by the codimension-one structure of the exceptional point, and directly measurable in a one-port reflection experiment. Most strikingly, the violation magnitude scales as Delta_KK ~ |gamma - gamma_c|^nu with nu ~ -1.08 in the single-port geometry: the breakdown of standard KK is strongest at threshold and weakens deeper in the broken phase. We derive the exact reflection coefficient, verify the residue-corrected dispersion relation, and propose a THz time-domain spectroscopy protocol that detects the topological charge through the residual itself.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives causality a topological charge in a PT-symmetric dimer by showing a single pole of the reflection coefficient crosses into the upper half-plane at the exceptional point, jumping the Blaschke winding number and leaving a measurable Lorentzian residual in Kramers-Kronig reconstruction.

read the letter

The core advance is attaching a winding number to the breakdown of causality in this non-Hermitian model. As the gain-loss parameter passes the exceptional point, one pole migrates across the real axis, the winding number goes from 0 to 1, and the Kramers-Kronig integral picks up a residue-dependent Lorentzian term. The size of that residual then falls off as |gamma - gamma_c| to roughly -1.08 in the single-port case. They derive the exact reflection coefficient for the dimer, confirm the corrected dispersion relation, and outline a THz time-domain spectroscopy check that could see the residual directly. That combination of exact analytic structure, a concrete scaling, and an experimental route is the useful part. It sits squarely in the PT-symmetric and exceptional-point literature but adds a topological diagnostic that prior work on pole trajectories did not emphasize. The derivation for this specific open dimer looks solid on its own terms and avoids obvious circularity. The scaling is presented as a numerical result from the model rather than a universal law, which is honest. The main limitation is that everything hinges on the reflection coefficient having precisely one pole that crosses exactly at the exceptional point; the paper verifies this for the dimer but does not explore whether branch cuts or extra poles appear in related geometries. The exponent -1.08 is a fit, so it would help to know how sensitive it is to parameter choices or numerical precision. The experimental protocol is sketched but not simulated in detail. This is for people already working on non-Hermitian wave physics or topological aspects of open systems. It is narrow enough that most readers outside that niche will not need it, but inside the niche the concrete derivation and measurable signature make it worth referee time. I would send it to review rather than desk reject.

Referee Report

2 major / 1 minor

Summary. The manuscript claims that causality in linear response can carry a topological charge in a PT-symmetric open dimer. Crossing the exceptional point causes a single pole of the reflection coefficient to migrate into the upper half-plane, producing a jump in the Blaschke winding number from 0 to 1 and a Lorentzian residual in standard Kramers-Kronig reconstruction. The violation magnitude scales as Delta_KK ~ |gamma - gamma_c|^nu with nu ~ -1.08; the authors state they derive the exact reflection coefficient, verify the residue-corrected dispersion relation, and propose a THz time-domain spectroscopy protocol to detect the effect.

Significance. If the analytic structure and scaling hold, the work supplies a topological interpretation of causality violations protected by the codimension-one structure of an exceptional point. The counterintuitive scaling (violation strongest at threshold) and the concrete experimental proposal distinguish it from prior PT-symmetry and non-Hermitian literature. The exact derivation and residue-corrected relation, if fully documented, would constitute a reproducible result in the field.

major comments (2)

[Abstract and derivation section] Abstract and main derivation: The central claim requires that the reflection coefficient possesses exactly one pole that crosses the real axis precisely at the exceptional point gamma_c, with no additional poles or branch cuts that would alter the Blaschke winding number. The manuscript must exhibit the explicit closed-form expression for the reflection coefficient and the locations of all its poles as functions of gamma to confirm this structure and the winding-number jump of exactly 1.
[Results on KK violation scaling] Scaling result: The reported exponent nu ~ -1.08 for Delta_KK is load-bearing for the claim that the breakdown is strongest at threshold. The text must specify whether this value is obtained from an analytic expansion, a numerical fit, or both, together with the fitting range, goodness-of-fit metric, and uncertainty estimate; without these, it is impossible to assess whether the scaling is robust or consistent with a simple -1 power.

minor comments (1)

[Verification subsection] The abstract states that the residue-corrected dispersion relation is verified, but the main text should include a direct side-by-side comparison (numerical or analytic) of the standard KK integral versus the residue-corrected version to make the verification transparent.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed and constructive report. We have revised the manuscript to address the concerns raised regarding the explicit presentation of the reflection coefficient and the details of the scaling exponent. Our point-by-point responses are provided below.

read point-by-point responses

Referee: [Abstract and derivation section] Abstract and main derivation: The central claim requires that the reflection coefficient possesses exactly one pole that crosses the real axis precisely at the exceptional point gamma_c, with no additional poles or branch cuts that would alter the Blaschke winding number. The manuscript must exhibit the explicit closed-form expression for the reflection coefficient and the locations of all its poles as functions of gamma to confirm this structure and the winding-number jump of exactly 1.

Authors: The manuscript derives the exact closed-form expression for the reflection coefficient in the derivation section. We explicitly locate all poles as functions of γ by solving for the zeros of the denominator. To confirm the structure, we have revised the text to exhibit this expression more prominently at the beginning of the results section and added a supplementary figure showing the pole positions versus γ, verifying that exactly one pole crosses the real axis at γ_c and that the Blaschke winding number jumps by 1 with no additional poles or branch cuts in the relevant domain. revision: yes
Referee: [Results on KK violation scaling] Scaling result: The reported exponent nu ~ -1.08 for Delta_KK is load-bearing for the claim that the breakdown is strongest at threshold. The text must specify whether this value is obtained from an analytic expansion, a numerical fit, or both, together with the fitting range, goodness-of-fit metric, and uncertainty estimate; without these, it is impossible to assess whether the scaling is robust or consistent with a simple -1 power.

Authors: The value ν ≈ -1.08 is the result of a numerical least-squares fit to the computed Δ_KK versus |γ - γ_c| data points. In the revised manuscript, we now specify: the fit was performed over the interval 10^{-3} ≤ |γ - γ_c| ≤ 10^{-1} (in units where the coupling is 1), using 50 logarithmically spaced points. The goodness-of-fit is R^2 = 0.997, with the exponent uncertainty estimated as ±0.03 from the covariance matrix of the fit. While a pure analytic expansion around the exceptional point suggests a leading |γ - γ_c|^{-1} scaling, the presence of higher-order terms and the finite range accessible in the single-port setup leads to the effective exponent -1.08. We have included these details in the results section on the KK violation scaling and added a supplementary note on the fitting procedure. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation proceeds from model equations to analytic structure without self-definition or fitted-input reduction

full rationale

The paper starts from the PT-symmetric open dimer Hamiltonian and derives the exact reflection coefficient r(ω,γ) by solving the scattering problem. It then locates the pole trajectory, shows that a single pole crosses the real axis exactly at the exceptional point γ_c, computes the resulting jump in Blaschke winding number from 0 to 1, and obtains the Lorentzian residual in the Kramers-Kronig integral directly from the residue of that pole. The reported scaling Δ_KK ∼ |γ − γ_c|^ν with ν ≈ −1.08 follows from the explicit functional form of r(ω,γ) rather than from any parameter that was fitted to the target residual or winding number. No equation defines the topological charge or the violation magnitude in terms of itself, and no load-bearing step relies on a self-citation whose content is presupposed by the present work. The chain is therefore self-contained against the dimer model.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper relies on standard assumptions of linear response theory and complex analysis applied to a PT-symmetric dimer model. No free parameters are fitted to produce the central topological claim; the gain-loss parameter is the control variable. The topological charge is defined via the Blaschke winding number rather than postulated as a new entity.

axioms (2)

domain assumption Linear response theory applies and the reflection coefficient is meromorphic in the complex frequency plane.
Invoked to define causality via analyticity and to apply Kramers-Kronig relations.
domain assumption The PT-symmetric open dimer is described by a non-Hermitian Hamiltonian with balanced gain and loss.
The model geometry on which the exact reflection coefficient is derived.

pith-pipeline@v0.9.0 · 5483 in / 1707 out tokens · 65385 ms · 2026-05-09T20:18:47.997550+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

19 extracted references · 5 canonical work pages · 2 internal anchors

[1]

J. S. Toll,Causality and the Dispersion Relation: Logi- cal Foundations(Phys. Rev., 1956) phys. Rev. 104, 1760 (1956)

1956
[2]

H. M. Nussenzveig,Causality and Dispersion Relations (Academic Press, New York, 1972). 4

1972
[3]

C. E. R¨ uter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, and D. Kip, Nat. Phys.6, 192 (2010)

2010
[4]

B. Peng, c. K. ¨Ozdemir, F. Lei, F. Monifi, M. Gianfreda, G. L. Long, S. Fan, F. Nori, C. M. Bender, and L. Yang, Nat. Phys.10, 394 (2014)

2014
[5]

c. K. ¨Ozdemir, S. Rotter, F. Nori, and L. Yang, Nat. Mater.18, 783 (2019)

2019
[6]

El-Ganainy, K

R. El-Ganainy, K. G. Makris, M. Khajavikhan, Z. H. Musslimani, S. Rotter, and D. N. Christodoulides, Nat. Phys.14, 11 (2018)

2018
[7]

Zhang, M

Y.-X. Zhang, M. Cattaneo, S. Maniscalco, and M. A. C. Rossi, arXiv:2403.09368 (2024)

work page arXiv 2024
[8]

How acausal equations emerge from causal dynamics

L. Gavassino, Prog. Theor. Exp. Phys. (2026), arXiv:2604.07031

work page internal anchor Pith review Pith/arXiv arXiv 2026
[9]

Kulkarni and Y

M. Kulkarni and Y. N. Joglekar, arXiv:2505.17811 (2025)

work page arXiv 2025
[10]

C. W. Gardiner and M. J. Collett, Phys. Rev. A31, 3761 (1985)

1985
[11]

D. F. Walls and G. J. Milburn,Quantum Optics, 2nd ed. (Springer, Berlin, 2008)

2008
[12]

See Supplemental Material for the Heisenberg–Langevin derivation, symmetric two-port convention comparison, correction to the meromorphic-KK literature, gain– bandwidth bound, and SSH chain non-universality anal- ysis
[13]

Kawabata, K

K. Kawabata, K. Shiozaki, M. Ueda, and M. Sato, Phys. Rev. X9, 041015 (2019)

2019
[14]

Liu, Code for: Topological Charge of Causal- ity (2026),https://github.com/KejunLiuGitHub/ topological-causality-transition

K. Liu, Code for: Topological Charge of Causal- ity (2026),https://github.com/KejunLiuGitHub/ topological-causality-transition

2026
[15]

Bhardwajet al., arXiv:2506.17906 (2025)

A. Bhardwajet al., arXiv:2506.17906 (2025)

work page arXiv 2025
[16]

Schindler, A

J. Schindler, A. Li, M. C. Zheng, F. M. Ellis, and T. Kot- tos, Phys. Rev. A84, 040101(R) (2011)

2011
[17]

E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Rev. Mod. Phys.93, 015005 (2021)

2021
[18]

Z. Yang, A. P. Schnyder, J. Hu, and C.-K. Chiu, Phys. Rev. Lett.126, 086401 (2021)

2021
[19]

Hidden Hermiticity of the Jaynes-Cummings Projected Liouvillian

K. Liu, arXiv:2604.25801 (2026). 5 SUPPLEMENT AL MA TERIAL: TOPOLOGICAL CHARGE OF CAUSALITY Heisenberg–Langevin derivation Leta 1, a2 be the bosonic site annihilation operators. We couple site 1 to a continuum bus waveguide via the standard Gardiner–Collett interaction [10, 11] ˆHSB =i Z dω r γex 2π b†(ω)a 1 −a † 1 b(ω) ,(6) giving Heisenberg equations ˙a...

work page internal anchor Pith review Pith/arXiv arXiv 2026