The Two Orbital, Interacting Hatano-Nelson Model
Pith reviewed 2026-05-10 10:29 UTC · model grok-4.3
The pith
The interacting two-orbital Hatano-Nelson model admits regions of parameter space where the full spectrum is real, set by interaction strength, hopping asymmetry, and interchain coupling.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the two-chain Hubbard-Hatano-Nelson model the spectrum is purely real inside certain regions of the three-dimensional parameter space spanned by interaction strength, non-Hermiticity, and interchain hopping. Winding-number analysis demonstrates that the spectral properties are sensitive to boundary conditions and that doublon states under periodic boundaries correspond to skin modes under open boundaries. At low filling the non-Hermitian Hamiltonian reproduces the essential features of the full Lindbladian time evolution.
What carries the argument
The two-orbital Hatano-Nelson Hamiltonian with asymmetric intra-chain hopping, symmetric inter-chain hopping, and Hubbard interaction, whose complex eigenvalues determine whether the spectrum collapses to the real axis.
If this is right
- Increasing Hubbard repulsion can drive the spectrum from complex to entirely real at fixed non-Hermiticity.
- Varying the interchain hopping strength moves the boundaries of the real-spectrum region in a systematic way.
- Winding numbers under periodic boundaries directly encode the skin-mode behavior seen under open boundaries.
- The non-Hermitian Hamiltonian supplies a qualitatively accurate shortcut for low-filling dissipative dynamics.
Where Pith is reading between the lines
- Interaction-tuned real spectra may appear generically in multi-orbital non-Hermitian lattices, offering a route to stabilize open many-body systems.
- The doublon-skin-mode correspondence suggests that bound pairs can be used as diagnostics for non-Hermitian topology in higher-dimensional or multi-band settings.
- If the low-filling agreement persists at moderate densities, effective non-Hermitian models could become practical simulators for a wider class of driven-dissipative Hubbard systems.
Load-bearing premise
Finite-lattice spectral features and the non-Hermitian effective description both survive the thermodynamic limit and remain faithful to the underlying Lindbladian dynamics at low filling.
What would settle it
A numerical diagonalization or Lindbladian simulation that finds persistent imaginary eigenvalues inside a parameter region the phase diagram labels as real-only, or that shows clear qualitative mismatch between non-Hermitian and Lindbladian dynamics even at very low particle density.
Figures
read the original abstract
The single orbital, one-dimensional, Hatano-Nelson Hamiltonian provides deep insight into the physics of non-Hermiticity, resulting from asymmetric left/right hopping, and its connections to localization. In the absence of disorder, its single particle eigenvalues $E_{\alpha}$ lie on an ellipse in the complex plane whose extent in the imaginary direction is controlled by the degree of asymmetry. When randomness is introduced, two sets of real eigenvalues emerge at the extremes of the largest and smallest real part of $E_{\alpha}$. These real eigenvalues are associated with localized eigenvectors. For spinless fermions, increasing near-neighbor interactions first cause a transition to a charge density wave phase, and ultimately, on finite lattices, a collapse of all eigenvalues to the real axis. In this paper, we explore the presence of real eigenvalues in the interacting, two-particle sector for the spinful case (Hubbard model) in a two-chain (two-band) geometry with a Hermitian interchain hopping. Our key results are to obtain the ``phase" diagrams for the existence of a purely real spectrum, as a function of the interaction strength, degree of non-Hermiticity, and interchain hopping. We study the sensitivity to boundary conditions of the spectral properties of our two-chain model with winding number analysis and explore the relationship between PBC doublon states and OBC skin modes. To address the question of stability in such non-equilibrium systems, we solve the dynamics at low filling according to Lindbladian evolution and find that the non-Hermitian description is able to qualitatively describe such systems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the two-orbital interacting Hatano-Nelson model (spinful Hubbard on two chains with asymmetric hopping and Hermitian interchain coupling). It reports phase diagrams in the two-particle sector for the region of fully real many-body spectra as functions of interaction strength U, non-Hermiticity asymmetry, and interchain hopping amplitude. Additional results include winding-number analysis of periodic- versus open-boundary spectral sensitivity, a proposed relation between PBC doublon states and OBC skin modes, and a qualitative comparison of non-Hermitian time evolution to Lindbladian dynamics at low filling.
Significance. If the reported real-spectrum regions survive extrapolation to the thermodynamic limit and the Lindbladian comparison can be made quantitative, the work would add to the understanding of interaction-tuned non-Hermitian skin effects and boundary-condition sensitivity in multi-orbital models. The direct diagonalization approach and winding-number diagnostics are standard tools, but the absence of scaling data limits the strength of the central claims.
major comments (3)
- [phase diagrams / results] The phase diagrams for purely real spectra (abstract and results section) are constructed on finite lattices with no reported system sizes L, no finite-size scaling, and no extrapolation of the critical surfaces in (U, asymmetry, interchain hopping). In non-Hermitian models the skin effect routinely produces strong L-dependence of eigenvalue trajectories; without scaling the claimed phase boundaries cannot be taken as thermodynamic-limit statements.
- [boundary-condition sensitivity] The winding-number analysis relating PBC doublons to OBC skin modes is presented without quantitative measures (e.g., winding number values, scaling of the imaginary-part gap with L, or explicit comparison of spectra for multiple L). This relation is load-bearing for the interpretation of the phase diagrams under open boundaries.
- [dynamics / Lindbladian evolution] The Lindbladian dynamics comparison is stated to be only qualitative (abstract and dynamics section) with no system sizes, no error estimates, no specific parameter values shown, and no quantitative metric (overlap, distance to real axis, etc.). The claim that the non-Hermitian effective Hamiltonian “qualitatively describes” the open-system evolution therefore rests on unverified numerical steps.
minor comments (2)
- [model definition] Notation for the two-orbital Hamiltonian (hopping terms, interaction U, asymmetry parameter) should be defined explicitly in the main text with a single equation reference rather than scattered definitions.
- [figures] Figure captions for the phase diagrams should include the lattice length(s) used and the number of disorder realizations or boundary-condition samples, if any.
Simulated Author's Rebuttal
We thank the referee for their careful reading of our manuscript and for the constructive comments, which help clarify the presentation of our results on the two-orbital interacting Hatano-Nelson model. We address each major comment point by point below, indicating the revisions planned for the next version.
read point-by-point responses
-
Referee: [phase diagrams / results] The phase diagrams for purely real spectra (abstract and results section) are constructed on finite lattices with no reported system sizes L, no finite-size scaling, and no extrapolation of the critical surfaces in (U, asymmetry, interchain hopping). In non-Hermitian models the skin effect routinely produces strong L-dependence of eigenvalue trajectories; without scaling the claimed phase boundaries cannot be taken as thermodynamic-limit statements.
Authors: We agree that the manuscript should explicitly report the lattice sizes and address finite-size effects, given the known sensitivity to the skin effect. The two-particle exact-diagonalization calculations were performed on finite chains with L ranging from 6 to 14 (depending on the parameter point to keep the Hilbert space manageable), but these values were not stated in the text or figure captions. In the revised manuscript we will add this information, include a new subsection discussing the L-dependence of the real-spectrum regions, and show representative data for at least three values of L. While a complete extrapolation of the critical surfaces to the thermodynamic limit (L→∞ at fixed particle number) would require substantial additional resources, the available data indicate that the boundaries separating real and complex spectra remain qualitatively stable for the largest accessible L; we will therefore qualify the phase diagrams as applying to finite but representative system sizes rather than claiming thermodynamic-limit validity. revision: partial
-
Referee: [boundary-condition sensitivity] The winding-number analysis relating PBC doublons to OBC skin modes is presented without quantitative measures (e.g., winding number values, scaling of the imaginary-part gap with L, or explicit comparison of spectra for multiple L). This relation is load-bearing for the interpretation of the phase diagrams under open boundaries.
Authors: We accept that the winding-number discussion would be strengthened by quantitative diagnostics. In the revised version we will report the explicit winding-number values obtained for the periodic-boundary doublon bands, plot the scaling of the imaginary-part gap versus L for several representative parameter sets, and include side-by-side spectral comparisons (PBC versus OBC) for at least two additional system sizes. These additions will make the link between PBC doublons and OBC skin modes more transparent and directly support the interpretation of the open-boundary phase diagrams. revision: yes
-
Referee: [dynamics / Lindbladian evolution] The Lindbladian dynamics comparison is stated to be only qualitative (abstract and dynamics section) with no system sizes, no error estimates, no specific parameter values shown, and no quantitative metric (overlap, distance to real axis, etc.). The claim that the non-Hermitian effective Hamiltonian “qualitatively describes” the open-system evolution therefore rests on unverified numerical steps.
Authors: The Lindbladian comparison was presented as a qualitative illustration of possible connections to open-system dynamics. We will revise the dynamics section to specify the system sizes and parameter values employed, report numerical error estimates from the time-integration routine, and introduce a simple quantitative metric (e.g., the time-averaged maximum imaginary part of the instantaneous eigenvalues or the L2 distance of the density matrix to the real-axis projection). These changes will allow readers to assess the degree of agreement more objectively while preserving the qualitative character of the original claim. revision: yes
Circularity Check
No circularity; results obtained via direct numerical diagonalization and time evolution on the defined Hamiltonian.
full rationale
The paper computes phase diagrams for real spectra in the two-particle sector by exact diagonalization of the interacting two-orbital Hatano-Nelson Hamiltonian on finite lattices, supplemented by winding-number analysis of boundary-condition sensitivity and Lindbladian dynamics at low filling. No parameters are fitted to data subsets and then repurposed as predictions of closely related quantities. No self-citations serve as load-bearing justifications for the central claims, and no ansatz or uniqueness theorem is smuggled in from prior author work. The reported spectra, phase boundaries, and qualitative dynamical agreement follow directly from the model definition and standard numerical methods without reduction to tautology or self-referential construction.
Axiom & Free-Parameter Ledger
free parameters (3)
- interaction strength
- non-Hermiticity asymmetry
- interchain hopping amplitude
axioms (2)
- domain assumption The non-Hermitian Hamiltonian provides a valid effective description of the open system
- domain assumption Finite-lattice exact diagonalization captures the relevant spectral features
Reference graph
Works this paper leans on
-
[1]
in the two-particle sector. In addition, Appendix B reviews the effect of boundary conditions for the single- chain Hatano-Nelson model and for the ladder model we investigate. B. Purely Real Spectra forN ↑ =N ↓ = 1 An equally simple analytical condition for the onset of this purely real regime in theU̸= 0 case does not exist; we therefore resort to numer...
-
[2]
J. Feinberg and A. Zee, Non-Hermitian localization and delocalization, Physical Review E59, 6433 (1999)
work page 1999
-
[3]
H. Shen, B. Zhen, and L. Fu, Topological band theory for non-Hermitian Hamiltonians, Physical Review Let- ters120, 146402 (2018)
work page 2018
-
[4]
A. Ghatak and T. Das, New topological invariants in non- Hermitian systems, Journal of Physics: Condensed Mat- ter31, 263001 (2019)
work page 2019
-
[5]
K. Kawabata, K. Shiozaki, M. Ueda, and M. Sato, Sym- metry and topology in non-Hermitian physics, Physical Review X9, 041015 (2019)
work page 2019
- [6]
-
[7]
E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Ex- ceptional topology of non-Hermitian systems, Reviews of Modern Physics93, 015005 (2021)
work page 2021
-
[8]
K. Ding, C. Fang, and G. Ma, Non-Hermitian topol- ogy and exceptional-point geometries, Nature Reviews Physics4, 745 (2022)
work page 2022
- [9]
-
[10]
C. M. Bender, Making sense of non-Hermitian Hamilto- nians, Reports on Progress in Physics70, 947 (2007)
work page 2007
-
[11]
L. Feng, R. El-Ganainy, and L. Ge, Non-Hermitian pho- tonics based on parity–time symmetry, Nature Photonics 11, 752 (2017)
work page 2017
-
[12]
R. El-Ganainy, K. G. Makris, M. Khajavikhan, Z. H. Musslimani, S. Rotter, and D. N. Christodoulides, Non- Hermitian physics and PT symmetry, Nature Physics14, 11 (2018). FIG. A10. Largest imaginary part of all eigenvalues for the 10×4 bulk system at theN ↑ =N ↓ = 1 sector, for fixed δ= 0.5, V= 0.5 and varyingU, V 0
work page 2018
-
[13]
I. Rotter, A non-Hermitian Hamilton operator and the physics of open quantum systems, Journal of Physics A: Mathematical and Theoretical42, 153001 (2009)
work page 2009
-
[14]
F. Song, S. Yao, and Z. Wang, Non-Hermitian skin effect and chiral damping in open quantum systems, Physical 16 Review Letters123, 170401 (2019)
work page 2019
-
[15]
N. Hatano and D. R. Nelson, Localization transitions in non-Hermitian quantum mechanics, Physical Review Letters77, 570 (1996)
work page 1996
-
[16]
N. Hatano and D. R. Nelson, Vortex pinning and non- Hermitian quantum mechanics, Physical Review B56, 8651 (1997)
work page 1997
-
[17]
N. Hatano and D. R. Nelson, Non-Hermitian delocal- ization and eigenfunctions, Physical Review B58, 8384 (1998)
work page 1998
-
[18]
In this paper, we will use a somewhat different convention for the hopping, namelyt±δin the two directions. See Eqs. 1-3
-
[19]
J. Feinberg and A. Zee, Spectral curves of non-Hermitian Hamiltonians, Nuclear Physics B552, 599 (1999)
work page 1999
-
[20]
P. L. Doussal, unpublished
- [21]
- [22]
-
[23]
W. N. Faugno and T. Ozawa, Interaction-induced non- Hermitian topological phases from a dynamical gauge field, Phys. Rev. Lett.129, 180401 (2022)
work page 2022
-
[24]
B. D´ ora and C. u. u. u. u. P. m. c. Moca, Full counting statistics in the many-body Hatano-Nelson model, Phys. Rev. B106, 235125 (2022)
work page 2022
-
[25]
K. Kawabata, K. Shiozaki, and S. Ryu, Many-body topology of non-Hermitian systems, Phys. Rev. B105, 165137 (2022)
work page 2022
-
[26]
F. Alsallom, L. Herviou, O. V. Yazyev, and M. Brzezi´ nska, Fate of the non-Hermitian skin ef- fect in many-body fermionic systems, Phys. Rev. Res.4, 033122 (2022)
work page 2022
-
[27]
S. Longhi, Spectral structure and doublon dissociation in the two-particle non-Hermitian Hubbard model, Annalen der Physik535, 2300291 (2023)
work page 2023
-
[28]
T. Orito and K.-I. Imura, Entanglement dynamics in the many-body Hatano-Nelson model, Phys. Rev. B108, 214308 (2023)
work page 2023
- [29]
-
[30]
A. Zee, A non-Hermitian particle in a disordered world, Physica A: Statistical Mechanics and its Applications 254, 300 (1998)
work page 1998
-
[31]
C. H. Lee, L. Li, and J. Gong, Hybrid higher-order skin- topological modes in nonreciprocal systems, Phys. Rev. Lett.123, 016805 (2019)
work page 2019
-
[32]
L. Li, C. H. Lee, S. Mu, and J. Gong, Critical non- Hermitian skin effect, Nature communications11, 5491 (2020)
work page 2020
-
[33]
T. Hayata and A. Yamamoto, Non-Hermitian Hubbard model without the sign problem, Phys. Rev. B104, 125102 (2021)
work page 2021
-
[34]
X.-J. Yu, Z. Pan, L. Xu, and Z.-X. Li, Non-hermitian strongly interacting Dirac fermions, Phys. Rev. Lett. 132, 116503 (2024)
work page 2024
-
[35]
K. Suthar, Y.-C. Wang, Y.-P. Huang, H. H. Jen, and J.-S. You, Non-Hermitian many-body localization with open boundaries, Phys. Rev. B106, 064208 (2022)
work page 2022
- [36]
- [37]
-
[38]
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
-
[39]
T. E. Lee, Anomalous edge state in a non-Hermitian lat- tice, Phys. Rev. Lett.116, 133903 (2016)
work page 2016
-
[40]
Y. Xiong, Why does bulk boundary correspondence fail in some non-Hermitian topological models, Journal of Physics Communications2, 035043 (2018)
work page 2018
- [41]
- [42]
-
[43]
Peierls, Zur Theorie des Diamagnetismus von Leitungselektronen, Zeitschrift fur Physik80, 763 (1933)
R. Peierls, Zur Theorie des Diamagnetismus von Leitungselektronen, Zeitschrift fur Physik80, 763 (1933)
work page 1933
-
[44]
J. M. Luttinger, The effect of a magnetic field on elec- trons in a periodic potential, Phys. Rev.84, 814 (1951)
work page 1951
-
[45]
S. Mu, C. H. Lee, L. Li, and J. Gong, Emergent Fermi surface in a many-body non-Hermitian fermionic chain, Phys. Rev. B102, 081115 (2020)
work page 2020
-
[46]
E. Lee, H. Lee, and B.-J. Yang, Many-body approach to non-hermitian physics in fermionic systems, Phys. Rev. B101, 121109 (2020)
work page 2020
-
[47]
For OBC, the no-jump effective Hamiltonian generated by Eq. (15) acquires a nonuniform imaginary onsite term, since edge sites participate in fewer jump operators than bulk sites. This can be remedied by introducing addi- tional boundary jump operators, equivalently viewed as couplings to fictitious vacuum sites outside the chain. With this completion, th...
-
[48]
J. Dalibard, Y. Castin, and K. Mølmer, Wave-function approach to dissipative processes in quantum optics, Phys. Rev. Lett.68, 580 (1992)
work page 1992
- [49]
-
[50]
R. Dum, P. Zoller, and H. Ritsch, Monte carlo simulation of the atomic master equation for spontaneous emission, Phys. Rev. A45, 4879 (1992)
work page 1992
-
[51]
R. Dum, A. S. Parkins, P. Zoller, and C. W. Gardiner, Monte carlo simulation of master equations in quan- tum optics for vacuum, thermal, and squeezed reservoirs, Phys. Rev. A46, 4382 (1992)
work page 1992
-
[52]
A. J. Daley, Quantum trajectories and open many-body quantum systems, Advances in Physics63, 77 (2014)
work page 2014
-
[53]
C. Wang, T.-C. Yi, J. Li, and R. Mondaini, Non- Hermitian Haldane-Hubbard model: Effective descrip- tion of one- and two-body dissipation, Phys. Rev. B108, 085134 (2023)
work page 2023
- [54]
-
[55]
K. Wang, A. Dutt, K. Y. Yang, C. C. Wojcik, J. Vuˇ ckovi´ c, and S. Fan, Generating arbitrary topologi- cal windings of a non-Hermitian band, Science371, 1240 (2021)
work page 2021
-
[56]
L. Yuan, A. Dutt, M. Qin, S. Fan, and X. Chen, Creat- ing locally interacting hamiltonians in the synthetic fre- quency dimension for photons, Photonics Research8, B8 (2020)
work page 2020
- [57]
-
[58]
G. Afek, D. Carney, and D. C. Moore, Coherent scat- tering of low mass dark matter from optically trapped sensors, Physical Review Letters128, 101301 (2022)
work page 2022
-
[59]
A. J. Brady, X. Chen, Y. Xia, J. Manley, M. Dey Chowd- hury, K. Xiao, Z. Liu, R. Harnik, D. J. Wilson, Z. Zhang, et al., Entanglement-enhanced optomechanical sensor ar- ray with application to dark matter searches, Communi- cations Physics6, 237 (2023)
work page 2023
-
[60]
M. M. Burns, J.-M. Fournier, and J. A. Golovchenko, Op- tical binding, Physical Review Letters63, 1233 (1989)
work page 1989
-
[61]
S. K. Mohanty, J. T. Andrews, and P. K. Gupta, Optical binding between dielectric particles, Optics Express12, 2746 (2004)
work page 2004
- [62]
-
[63]
J. Vijayan, J. Piotrowski, C. Gonzalez-Ballestero, K. We- ber, O. Romero-Isart, and L. Novotny, Cavity-mediated long-range interactions in levitated optomechanics, Na- ture Physics20, 1 (2024)
work page 2024
-
[64]
H. Rudolph, U. Deli´ c, K. Hornberger, and B. A. Stick- ler, Quantum theory of non-Hermitian optical binding between nanoparticles, arXiv preprint arXiv:2306.11893 (2023)
-
[65]
V. Liˇ ska, T. Zem´ ankov´ a, V. Svak, P. J´ akl, J. Jeˇ zek, M. Br´ aneck` y, S. H. Simpson, P. Zem´ anek, and O. Br- zobohat` y, Cold damping of levitated optically coupled nanoparticles, Optica10, 1203 (2023)
work page 2023
- [66]
-
[67]
M. Reisenbauer, H. Rudolph, L. Egyed, K. Hornberger, A. V. Zasedatelev, M. Abuzarli, B. A. Stickler, and U. Deli´ c, Non-Hermitian dynamics and nonreciproc- ity of optically coupled nanoparticles, arXiv preprint arXiv:2310.02610 (2023)
-
[68]
Y. Xia, A. R. Agrawal, C. M. Pluchar, A. J. Brady, Z. Liu, Q. Zhuang, D. J. Wilson, and Z. Zhang, Entanglement-enhanced optomechanical sensing, Nature Photonics17, 470 (2023)
work page 2023
-
[69]
The Two Orbital, Interacting Hatano-Nelson Model
Data and figures for “The Two Orbital, Interacting Hatano-Nelson Model” (2026), Zenodo repository con- taining all data and figures
work page 2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.