Natural-orbital locking reveals hidden steady-state skin order in Gaussian open fermion chains
Pith reviewed 2026-05-08 04:16 UTC · model grok-4.3
The pith
In Gaussian open fermion chains, the dominant natural orbital locks to the Euclidean-normalized slow right eigenmode of the relaxation matrix.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The steady-state correlator admits a biorthogonal decomposition in terms of the left and right eigenmodes of the relaxation matrix X and the source matrix Y. This formula separates three ingredients: slow rapidity denominators, source loading by left eigenmodes, and real-space geometry from right eigenmodes. For a local pump, the pump position is read by the left modes, whereas the selected profile is drawn by the right modes. In a single-slow-mode regime, the dominant natural orbital locks to the Euclidean-normalized slow right mode.
What carries the argument
Biorthogonal decomposition of the steady-state correlator in left and right eigenmodes of relaxation matrix X and source matrix Y, which produces natural-orbital locking to the slow right mode.
If this is right
- The density profile is an incoherent sum over occupied natural orbitals and therefore reveals skin order less selectively than the dominant orbital.
- In a nonreciprocal Hatano-Nelson chain the locking holds and directly images the skin-localized right mode.
- In a nonreciprocal SSH chain the locked orbital crosses over from a topological edge state to a slow bulk-skin state as parameters vary.
- Natural-orbital locking supplies a mode-resolved diagnostic that distinguishes nonreciprocal localization from other steady-state features.
Where Pith is reading between the lines
- The same locking could serve as an experimental probe of hidden skin effects in engineered open systems where only steady-state correlations are accessible.
- If analogous decompositions exist in non-Gaussian or interacting open systems, natural-orbital locking might diagnose nonreciprocal order more broadly.
- The separation of left-mode source loading from right-mode geometry suggests a general way to read pump location separately from localization profile in driven-dissipative chains.
Load-bearing premise
A single-slow-mode regime exists in which the dominant natural orbital can be cleanly isolated from all other contributions.
What would settle it
Compute or measure the correlation matrix for a nonreciprocal chain in the single-slow-mode regime and check whether its leading eigenvector coincides with the right eigenmode of X after Euclidean normalization.
Figures
read the original abstract
Nonreciprocal relaxation matrices can have skin-localized right eigenmodes, but their imprint on a mixed steady state is not fixed by the density profile alone. We develop an exact steady-state theory for number-conserving Gaussian fermion chains and show that the dominant natural orbital of the correlation matrix provides a mode-resolved diagnostic of hidden skin order. The steady-state correlator admits a biorthogonal decomposition in terms of the left and right eigenmodes of the relaxation matrix $X$ and the source matrix $Y$. This formula separates three ingredients: slow rapidity denominators, source loading by left eigenmodes, and real-space geometry from right eigenmodes. For a local pump, the pump position is read by the left modes, whereas the selected profile is drawn by the right modes. In a single-slow-mode regime, the dominant natural orbital locks to the Euclidean-normalized slow right mode. The density can follow the same boundary trend, but it is a less selective incoherent sum over occupied natural orbitals. We verify this selection law in a nonreciprocal Hatano--Nelson chain and show that, in a nonreciprocal SSH chain, the selected natural orbital crosses over from a topological edge candidate to a slow bulk-skin candidate. These results identify natural-orbital locking as a steady-state diagnostic of nonreciprocal localization in Gaussian open fermion chains.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops an exact steady-state theory for number-conserving Gaussian open fermion chains with nonreciprocal relaxation. It derives a biorthogonal decomposition of the steady-state correlator C in the left and right eigenbases of the relaxation matrix X and source matrix Y, separating rapidity denominators, left-mode source loading, and right-mode spatial profiles. In the explicitly defined single-slow-mode regime, the dominant natural orbital locks to the Euclidean-normalized slow right eigenmode, providing a mode-resolved diagnostic of hidden skin order that is not captured by the density profile alone. The selection law is verified numerically in a nonreciprocal Hatano–Nelson chain and shown to exhibit a crossover in a nonreciprocal SSH chain.
Significance. If the central claims hold, the work supplies a useful diagnostic for nonreciprocal localization effects in the steady states of open quantum systems that goes beyond incoherent density profiles. The exact biorthogonal decomposition and the clean numerical tests in two concrete models constitute clear strengths; the result could inform studies of driven-dissipative fermionic systems and non-Hermitian skin physics.
major comments (1)
- The single-slow-mode regime is load-bearing for the locking claim. The manuscript should supply quantitative bounds (e.g., the ratio of the dominant rapidity to the next one) and an error estimate showing when the dominant natural orbital deviates from the normalized slow right mode; without this, the regime remains an assumption whose range of validity is not fully delimited.
minor comments (1)
- The abstract and introduction would benefit from a brief statement of the precise conditions (e.g., local pump, number conservation) under which the biorthogonal decomposition holds exactly.
Simulated Author's Rebuttal
We thank the referee for their positive evaluation and constructive feedback on our manuscript. We address the single major comment below and will revise the manuscript to incorporate quantitative bounds and error estimates as suggested.
read point-by-point responses
-
Referee: The single-slow-mode regime is load-bearing for the locking claim. The manuscript should supply quantitative bounds (e.g., the ratio of the dominant rapidity to the next one) and an error estimate showing when the dominant natural orbital deviates from the normalized slow right mode; without this, the regime remains an assumption whose range of validity is not fully delimited.
Authors: We agree that the single-slow-mode regime requires a more precise quantitative characterization to fully support the locking claim. In the revised manuscript we will add an explicit definition of the regime together with bounds on the rapidity ratio (e.g., requiring |λ₁/λ₂| ≪ 1, where λ₁ is the smallest rapidity) and a perturbative error estimate for the deviation between the dominant natural orbital and the Euclidean-normalized slow right eigenmode. The estimate follows directly from truncating the biorthogonal expansion of the steady-state correlator after the leading term and bounding the remainder by the ratio of the next rapidity and the source projections onto the left modes. We will also include additional numerical scans in both the Hatano–Nelson and SSH models that map the crossover as a function of the rapidity ratio, thereby delimiting the regime’s practical range of validity. revision: yes
Circularity Check
No significant circularity in the derivation chain
full rationale
The paper derives an exact biorthogonal decomposition of the steady-state correlator directly from the linear structure of the relaxation matrix X and source matrix Y (presumably via the Lyapunov equation for Gaussian open systems). The natural-orbital locking then follows as an algebraic consequence under the explicitly stated single-slow-mode regime, with no fitted parameters, self-referential definitions, or load-bearing self-citations. Numerical checks in the Hatano-Nelson and SSH models provide independent validation of the selection law rather than circular confirmation. The derivation remains self-contained from the model equations without reducing to its inputs by construction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The system consists of number-conserving Gaussian fermion chains
- domain assumption The steady-state correlator admits a biorthogonal decomposition in terms of left and right eigenmodes of the relaxation matrix X and source matrix Y
Reference graph
Works this paper leans on
-
[1]
Natural-orbital locking reveals hidden steady-state skin order in Gaussian open fermion chains
The density is reconstructed as nss j = X α να|φα(j)|2.(3) Thus the density is not a mode-resolved object. It is an incoherent sum over all occupied natural orbitals. When one occupation dominates, the density may look similar to the leading orbital. When subleading orbitals carry appreciable weight, the density becomes smoother. In either case, the domin...
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[2]
N. Hatano and D. R. Nelson, Localization transitions in non-Hermitian quantum mechanics, Phys. Rev. Lett.77, 570 (1996)
work page 1996
-
[3]
N. Hatano and D. R. Nelson, Vortex pinning and non- Hermitian quantum mechanics, Phys. Rev. B56, 8651 (1997)
work page 1997
- [4]
-
[5]
Z. Gong, Y. Ashida, K. Kawabata, K. Takasan, S. Hi- gashikawa, and M. Ueda, Topological phases of non- Hermitian systems, Phys. Rev. X8, 031079 (2018)
work page 2018
-
[6]
E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Ex- ceptional topology of non-Hermitian systems, Rev. Mod. Phys.93, 015005 (2021)
work page 2021
-
[7]
T. Prosen, Third quantization: a general method to solve master equations for quadratic open Fermi systems, New J. Phys.10, 043026 (2008)
work page 2008
-
[8]
Prosen, Spectral theorem for the Lindblad equation for quadratic open fermionic systems, J
T. Prosen, Spectral theorem for the Lindblad equation for quadratic open fermionic systems, J. Stat. Mech. (2010) P07020
work page 2010
-
[9]
F. Song, S. Yao, and Z. Wang, Non-Hermitian skin effect and chiral damping in open quantum systems, Phys. Rev. Lett.123, 170401 (2019)
work page 2019
-
[10]
T. Haga, M. Nakagawa, R. Hamazaki, and M. Ueda, Li- ouvillian skin effect: slowing down of relaxation processes without gap closing, Phys. Rev. Lett.127, 070402 (2021)
work page 2021
-
[11]
F. Yang, Q.-D. Jiang, and E. J. Bergholtz, Liouvillian skin effect in an exactly solvable model, Phys. Rev. Re- search4, 023160 (2022)
work page 2022
-
[12]
Natural-orbital locking reveals hidden steady-state skin order in Gaussian open fermion chains
A. McDonald, R. Hanai, and A. A. Clerk, Nonequilib- rium stationary states of quantum non-Hermitian lattice models, Phys. Rev. B105, 064302 (2022). 6 Supplemental Material for “Natural-orbital locking reveals hidden steady-state skin order in Gaussian open fermion chains” Y. T. Wang and X. Z. Zhang College of Physics and Materials Science, Tianjin Normal ...
work page 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.