## Quantized electrical, thermal, and spin transports of non-Hermitian clean and dirty two-dimensional topological insulators and superconductors

Sanjib Kumar Das 1 and Bitan Roy 1

1 Department of Physics, Lehigh University, Bethlehem, Pennsylvania, 18015, USA

(Dated: December 5, 2024)

From lattice-regularized models, devoid of any non-Hermitian (NH) skin effects, here we compute the electrical ( σ xy ), thermal ( κ xy ), and spin ( σ sp xy ) Hall, and the electrical ( G xx ) and thermal ( G th xx ) longitudinal conductivities for appropriate NH planar topological insulators and superconductors related to all five non-trivial Altland-Zirbauer symmetry classes in their Hermitian limits. These models feature real eigenvalues over an extended NH parameter regime, only where the associated topological invariants remain quantized. In this regime, the NH quantum anomalous and spin Hall insulators show quantized σ xy and G xx , respectively, the NH p + ip ( p ± ip ) pairing shows halfquantized κ xy ( G th xx ), while the NH d + id pairing shows quantized κ xy and σ sp xy in the clean and weak disorder (due to random pointlike charge impurities) regimes. We compute these quantities in experimentally realizable suitable six-terminal setups using the Kwant software package. But, in the strong disorder regime, all these topological responses vanish and with the increasing nonHermiticity in the system this generic phenomenon occurs at weaker disorder.

Introduction . Topological crystals displaying the hallmark bulk-boundary correspondence in terms of robust gapless modes living on their boundaries [1-21] typically encounter the notorious non-Hermitian (NH) skin effect in open systems interacting with a bath or the environment [22-45]. It corresponds to an accumulation of all the right and left eigenvectors at opposite ends of the system with open boundary conditions. Even though topological invariants can be defined in NH systems their signatures on the the boundary modes get masked by the NH skin effect. Theoretically, this obstacle can be bypassed through the bi-orthogonal bulk-boundary correspondence, devoid of any NH skin effect [32]. Nevertheless, any direct manifestation of the topological invariants on experimentally measurable quantities even in an NH toy topological model remains illusive [46], leaving aside their actual measurements in open quantum systems. In this quest, two-dimensional NH topological models stand as cornerstones, since in the Hermitian or closed systems their topological invariants can be computed and measured in six-terminal Hall-bar setups [47-56].

Here, we consider skin effect-free NH generalization of square lattice models for all five non-trivial Hermitian Altland-Zirnbauer symmetry classes [57] that over an extended NH parameter regime feature real eigenvalue spectrum and non-vanishing topological invariants [58]. Using the Kwant software package [59-62], we numerically compute their electrical ( σ xy ), thermal ( κ xy ), and spin ( σ sp xy ) Hall, and electrical ( G xx ) and thermal ( G th xx ) longitudinal conductivities in six-terminal setups (Fig. 1), to arrive at the following conclusions.

The NH quantum anomalous Hall insulator (QAHI) and quantum spin Hall insulator (QSHI) belonging to class A and class AII in the Hermitian limit, respectively, show quantized σ xy and G xx in units of e /h 2 (Fig. 2). Here e ( h ) is the electrical charge (Planck's constant). By contrast, the NH p + ip ( p ± ip ) paired state, belonging to class D (class DIII) in Hermitian systems, supports halfquantized κ xy ( G th xx ) in units of κ 0 = π k 2 2 B T/ (3 h ) as the temperature T → 0, where k B is the Boltzmann constant (Fig. 3). Finally, the NH d + id paired state, falling in class C in the Hermitian limit, accommodates quantized κ xy (in units of κ 0 ) and σ sp xy (in units of σ sp 0 = ℏ / (8 π )). See Fig. 4. These conclusions hold in the entire topological regime in a clean and weakly disordered (due to pointlike random charge impurities) systems. In the strong disorder regime, all these responses vanish. However, with increasing t 2 or non-Hermiticity in the system, the onset of vanishing topological responses occur at weaker disorder (Fig. 5). This phenomenon takes place via quantum phase transitions at finite disorder, but the associated critical disorder strength ( W c ) can only be pinned from numerical simulations in sufficiently large systems, which is beyond the scope of our numerical resources.

Model . For all the gapped topological phases belonging to any one of the five non-trivial Altland-Zirnbauer symmetry classes in two dimensions, the model Hamiltonian on a square lattice take the following universal form

$$\underset { t \text{ over} } { \text{termi-} } \quad H = \alpha \sum _ { k } \Psi _ { k } ^ { \dagger } \left ( \sum _ { j = 1 } ^ { 3 } d _ { j } ( k ) \Gamma _ { j } \right ) \Psi _ { k } \equiv \alpha \sum _ { k } \Psi _ { k } ^ { \dagger } \mathcal { H } ( k ) \Psi _ { k }. \ ( 1 )$$

The internal structure of the spinor Ψ k with momentum k = ( k , k x y ) and representation of the Hermitian Γ matrices depend on the symmetry class. They always satisfy the anticommuting Clifford algebra { Γ j , Γ k } = 2 δ jk . The third component of the d k ( )-vector takes the form d 3 ( k ) = m 0 -t 0 [cos( k a x ) + cos( k a y )], and unless otherwise stated d 1 ( k ) = t 1 sin( k a x ) and d 2 ( k ) = t 1 sin( k a y ), where a is the lattice spacing. Within the topological regime ( | m /t 0 0 | &lt; 2), for topological insulators (superconductors) d 3 ( k ) features a band inversion (Fermi surface) near the Γ = (0 0) , and M = (1 1) , π/a points of the Brillouin zone (BZ) when 0 &lt; m /t 0 0 &lt; 2 and

FIG. 1. Six-terminal setup with a rectangular scattering region (black sites) of length L and width D connected to six semi-infinite leads (red sites), for the calculation of electrical ( σ xy ), thermal ( κ xy ), and spin ( σ sp xy ) Hall, and electrical ( G xx ) and thermal ( G th xx ) longitudinal conductivities. For the computation of σ xy and G xx , an electrical current ( I el ) flows due to a voltage bias (∆ x ≡ ∆ ) between the horizontal leads, V and we extract voltages x i ≡ V i on the vertical leads with i = 2 3 5 6. , , , Due to a temperature gradient (∆ x ≡ ∆ ), a T thermal current ( I th ) flows between the horizontal leads and temperature develops at the vertical leads with x i ≡ T i , when we compute κ xy and G th xx . While calculating σ sp xy , a spin current ( I sp ) flows between the horizontal leads, subject to a magnetic field bias with ∆ x ≡ ∆ H , yielding magnetization at the vertical leads with x i ≡ m i . The scattering region is maintained at a fixed voltage V (for σ xy and G xx ) or temperature T (for κ xy and G th xx ) or magnetic field H (for σ sp xy ). For all calculations (Figs. 2-5), we set L = 180 and D = 90.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00763v2-figure-bd6f2e2a-0dae-4741-820c-2a3cf4706d95.jpg)

-2 &lt; m /t 0 0 &lt; 0, respectively, with α = 1 ( α = 1 2 due / to the Nambu doubling). Throughout, we set t 1 = t 0 = 1.

Note that d 3 ( k ) and Γ 3 do not break any crystal symmetry, and H ⊥ ( k ) = d 1 ( k )Γ 1 + d 2 ( k )Γ 2 is also invariant under all the crystal symmetries. So, H ( k ) transforms under the trivial A 1 g representation. Then the anti-Hermitian operator H AH ( k ) = Γ 3 H ⊥ ( k ) /t 1 does not break any crystallographic symmetry either. Finally, we introduce the desired NH topological operator

$$\mathcal { H } _ { \text{NH} } ( k ) = \mathcal { H } ( k ) + t _ { 2 } \mathcal { H } _ { \text{AH} } ( k ) \quad \quad ( 2 ) \quad \text{we}$$

that is symmetric under all crystal symmetries, such as reflections, four-fold rotations, and inversion. Hence, it is devoid of any NH skin effect for all parameter values [58]. Here, t 2 determines the strength of non-Hermiticity in the system. Eigenvalues of H NH ( k ) are ± E NH ( k ), where

$$E _ { \text{NH} } ( k ) = \left [ ( t _ { 1 } ^ { 2 } - t _ { 2 } ^ { 2 } ) \{ d _ { 1 } ^ { 2 } ( k ) + d _ { 2 } ^ { 2 } ( k ) \} + d _ { 3 } ^ { 2 } ( k ) \right ] ^ { 1 / 2 }. \ ( 3 ) \quad \text{(I_{eI} ) } _ { \text{the } v }$$

They are purely real when t 2 &lt; t 1 . Then H NH ( k ) fosters non-trivial topological invariants, when | m /t 0 0 | &lt; 2. Our construction from Eq. (2) is applicable to any lattice systems, such as the triangular one, where the parent Hamiltonian also assumes the universal form of Eq. (1) [14]. The resulting NH operator respects the corresponding crystal symmetries, such as three-fold rotations.

FIG. 2. Electrical Hall conductivity ( σ xy ) of an NH QAHI for (a) a fixed strength of non-Hermiticity ( t 2 ) as a function of m 0 and (b) two fixed values of m 0 with varying t 2 , showing σ xy = C in units of e /h 2 , where C is the first Chern number. Electrical longitudinal conductivity ( G xx ) of an NH QSHI for (c) a fixed t 2 as a function of m 0 and (d) two fixed values of m 0 with varying t 2 , showing G xx = C sp (in units of e /h 2 ), where C sp is the spin Chern number. For t 2 &gt; 1, σ xy = G xx = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00763v2-figure-c138f609-9eb8-472a-9a78-105457c83e19.jpg)

NH QAHI . For Γ j = τ j and Ψ ⊤ k = ( c + , c -)( k ), where c τ ( k ) is the fermionic annihilation operator on orbital with parity eigenvalue τ = ± and momentum k , we find the NH incarnation of the Qi-Wu-Zhang model for QAHI [63]. The Pauli matrices { τ j } act on the orbital index. Then the NH operator takes the form H NH ( k ) = ∑ 3 j =1 d NH j ( k ) τ j , with d NH 1 ( k ) = t 1 sin( k a x ) -it 2 sin( k a y ), d NH 2 ( k ) = t 1 sin( k a y ) + it 2 sin( k a x ), and d NH 3 ( k ) = d 3 ( k ). The first Chern number of H NH ( k ) reads [58, 64]

$$\underset { \text{we} } { \text{the} } \quad C = \int _ { \text{BZ} } \frac { d ^ { 2 } k } { 4 \pi } \left [ \partial _ { k _ { z } } \hat { d } ^ { \text{NH} } ( k ) \times \partial _ { k _ { \nu } } \hat { d } ^ { \text{NH} } ( k ) \right ] \cdot \hat { d } ^ { \text{NH} } ( k ), \ ( 4 )$$

where ˆ d NH ( k ) = d NH ( k ) / √ [ d NH ( k )] 2 . The momentum integral is performed over the first BZ. For any | t 2 | &lt; 1, we find C = -1 (+1) for 0 &lt; m /t 0 0 &lt; 2 ( -2 &lt; m /t 0 0 &lt; 0). Next, we discuss the ramification of non-trivial C on the quantized σ xy , computed in a six-terminal setup (Fig. 1), detailed in the Supplemental Material (SM) [65], with all the numerical codes available on Zenodo [66].

A voltage gradient is applied between two horizontal leads, yielding a longitudinal electrical current ( I el ) between them. The transverse leads, serving as the voltage probes, carry no electrical current. The current-voltage relation is then given by I el = G V el , with V ⊤ = ( -∆ V/ , V , V 2 2 3 , ∆ V/ , V , V 2 5 6 ) and I ⊤ el = ( I el , 0 0 , , -I el , 0 0). , Upon computing the conductance matrix G el using Kwant [59], containing only the transmission blocks of the scattering matrix, we extract different voltages from the current-voltage relation. Subsequently, we compute the transverse electrical resistance

FIG. 3. Thermal Hall conductivity ( κ xy ) of an NH p + ip paired state for (a) a fixed t 2 as a function of m 0 and (b) two fixed values of m 0 with varying t 2 , showing half-quantization κ xy = C/ 2 (in units of κ 0 ). Thermal longitudinal conductivity ( G th xx ) of an NH p ± ip superconductor for (c) a fixed t 2 with varying m 0 and (d) two fixed values of m 0 with varying t 2 , showing G th xx = C sp / 2 (in units of κ 0 ). For t 2 &gt; 1, κ xy = G th xx = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00763v2-figure-96966c0e-f6e3-46fd-9347-0a12c61043f8.jpg)

R el xy = ( V 2 + V 3 -V 5 -V 6 ) / (2 I el ), and find σ xy = 1 /R el xy = C (in units of e /h 2 ). See Figs. 2(a) and 2(b), and SM [65].

NH QSHI . An NH generalization of the BernevigHughes-Zhang model for QSHI is realized with Γ = (Γ 01 , Γ 32 , Γ 03 ), where Γ µν = σ τ µ ν and the Pauli matrices { σ µ } operate on the spin index [5]. The four-component spinor reads as Ψ ⊤ k = ( c ↑ + , c ↑ -, c ↓ + , c ↓ -)( k ), where c σ τ ( k ) is the fermionic annihilation operator with parity τ = ± , spin projection σ = ↑ ↓ , in the z direction, and momentum k . Following the same steps, highlighted for the NH QAHI model, we compute the Chern number for individual spin components, given by C ↑ and C ↓ . Within the topological regime, we find that the total Chern number of the model is C total = C ↑ + C ↓ = 0. Nevertheless, we can define an invariant, the spin Chern number C sp = | C ↑ -C ↓ | which is non-trivial in the entire topological regime ( | m /t 0 0 | &lt; 2 and | t 2 | &lt; 1), given by C sp = 2.

The non-trivial C sp leaves its footprint on the quantized G xx for which we employ the same six-terminal arrangement, previously discussed for the NH QAHI. In this case, σ xy = 0 as C total = 0 therein. But, from the same current-voltage relationship, we now compute the longitudinal electrical resistance R xx = ( V 3 -V 2 ) /I el = ( V 5 -V 6 ) /I el , in turn yielding G xx = R -1 xx = C sp (in units of e /h 2 ), as shown in Figs. 2(c) and 2(d), due to the associated counter-propagating helical edge modes.

NH p + ip pairing . For an NH p + ip paired state among spinless or spin-polarized fermions Ψ ⊤ k = ( c k , c † -k ), where c k ( c † k ) is the fermionic annihilation (creation) opera-

FIG. 4. Thermal Hall conductivity ( κ xy ) of an NH d + id pairing for (a) a fixed t 2 with varying m 0 and (b) two fixed values of m 0 with varying t 2 , showing κ xy = C (in units of κ 0 ), where C = 2. Its spin Hall conductivity ( σ sp xy ) for (c) a fixed t 2 as a function of m 0 and (d) two fixed values of m 0 with varying t 2 shows σ sp xy = C (in units of σ sp 0 ). For t 2 &gt; 1, κ xy = σ sp xy = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00763v2-figure-93ea680b-823c-4e31-9988-97ec1c30cbc0.jpg)

tor with momentum k . The Pauli matrices { τ j } operate on the Nambu or particle-hole index [21]. The effective single-particle NH Bogoliubov de-Gennes operator H NH ( k ) has the topological invariant C [Eq. (4)]. It manifests via a half-quantized κ xy . Now, a thermal current ( I th ) flows between two horizontal leads, held at fixed but different temperatures (Fig. 1). Four vertical or transverse leads serve as the temperature probe, carry no thermal current. The thermal current-temperature relation is a matrix equation I th = A T th , where I ⊤ th = ( I th , 0 0 , , -I th , 0 0) , and T ⊤ = ( -∆ T/ , T , T 2 2 3 , ∆ T/ , T , T 2 5 6 ). The elements of A th are

$$\intertext { n \, the } \int _ { 0 } ^ { \infty } \frac { E ^ { 2 } } { T } \left ( - \frac { \partial f ( E, T ) } { \partial E } \right ) \left [ \delta _ { i j } \mu _ { j } - \text{Tr} ( \mathfrak { t } _ { i j } ^ { \dagger } \mathfrak { t } _ { i j } ) \right ] d E, \\ \intertext { r t h e - \frac { \dots } { \dots } \, \ e / \varepsilon \, \varepsilon \, \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsi \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdot \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e / \varepsilon \, \cdots \, \ e }$$

where f ( E,T ) = 1 (1 + exp[ / E/ k ( B T )]) is the FermiDirac distribution function, µ j denotes the number of propagating modes in the j th lead, t ij is the transmission part of the scattering matrix between the leads i and j , and the trace (Tr) is taken over the conducting channels. From A th , we extract temperatures at various leads. The transverse thermal resistance is then given by R th xy = ( T 2 + T 3 -T 5 -T 6 ) / (2 I th ) [60-62, 67], from which we compute κ xy = 1 /R th xy at T = 0 01, showing . k xy = αC in units of κ 0 (Figs. 3(a) and 3(b), and SM [65]), where α = 1 2 accounts for the Nambu doubling. Thus, model / NH p + ip paired state features a half-quantized κ xy in the entire topological regime.

NH p ± ip pairing . This paired state occurs among spin-1/2 fermions, with p + ip and p -ip pairing symmetries for opposite spin projections. The associated

FIG. 5. Disorder averaged (a) electrical Hall conductivity ⟨ σ xy ⟩ in QAHI, (b) electrical longitudinal conductivity ⟨ G xx ⟩ in QSHI, (c) thermal Hall conductivity ⟨ κ xy ⟩ in p + ip paired state, (d) thermal longitudinal conductivity ⟨ G th xx ⟩ in p ± ip superconductor, and (e) ⟨ κ xy ⟩ and (f) spin Hall conductivity ⟨ σ sp xy ⟩ in d + id paired state in Hermitian ( t 2 = 0) and NH ( t 2 = 0 1 0 5) systems. The . , . error bars correspond to the saturated (after averaging over a large number of disorder realizations) standard deviations [65].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00763v2-figure-a403485b-2fa1-43e1-9454-00ccbd731938.jpg)

four-component spinor is Ψ k = ( ψ , σ k 2 ψ ⋆ -k ) with ψ k = ( c ↑ , c ↓ )( k ), where c σ ( k ) is the fermionic annihilation operator with spin projection σ = ↑ ↓ , and momentum k . The involved Γ matrices are Γ 1 = σ τ 0 1 , Γ 2 = σ τ 3 2 , and Γ 3 = σ τ 0 3 . Similar to the situation in NH QSHI, C total = 0 but C sp = 2 for the NH p ± ip paired. We numerically confirm that now κ xy = 0. Nevertheless, from the longitudinal thermal resistance R th xx = ( T 3 -T 2 ) /I th = ( T 6 -T 5 ) /I th , we find that G th xx = ( R th xx ) -1 = αC sp κ 0 with α = 1 2, as shown in Figs. 3(c) and 3(d). /

compute from the spin current-magnetization matrix relation I sp = G M sp , where I ⊤ sp = ( I sp , 0 0 , , -I sp , 0 0) and , M ⊤ = ( -∆ H/ ,m ,m , 2 2 3 ∆ H/ ,m ,m 2 5 6 ), and G sp is the spin conductance matrix, extracted using Kwant [62]. From its solutions, we compute the spin Hall resistance R sp xy = ( m 2 + m 3 -m 5 -m / 6 ) (2 I sp ), yielding σ sp xy = ( R sp xy ) -1 = C (in units of σ sp 0 ), see Figs. 4(c) and 4(d).

NH d + id pairing . In the same Nambu-doubled spinor basis, a spin-singlet NH d + id paired state is realized with d 1 ( k ) = t 1 [cos( k a x ) -cos( k a y )], d 2 ( k ) = t 1 sin( k a x ) sin( k a y ), and Γ j = σ τ 0 j [21, 68]. Within the topological regime ( | m /t 0 0 | &lt; 2 and | t 2 | &lt; 1), the Chern number of the corresponding NH operator for each spin projection is C = 2, and thus its total Chern number is 2 C = 4. This NH paired state accommodates κ xy = C (in units of κ 0 ), as shown in Fig. 4(a) and 4(b). The spin-charge separation allows this NH paired state to harbor quantized spin Hall conductivity ( σ sp xy ), for which a spin current ( I sp ) passes through the scattering region due to a difference in the magnetic field bias between the two horizontal leads. Then magnetization ( m j ) develops in the vertical leads, which we

Disorder . Finally, we unfold the effects of disorder on all the topological responses, discussed so far in clean NH systems. We consider only pointlike charge impurities, the dominant source of elastic scattering in any real material. In the above NH systems, we introduce the terms V ( r ) τ 0 , V ( r ) σ τ 0 0 , V ( r ) τ 3 , V ( r ) σ τ 0 3 , and V ( r ) σ τ 0 3 , respectively, where V ( r ) is uniformly and randomly distributed within the range [ -W/ ,W/ 2 2] at each site of the scattering region and W denotes the strength of disorder. The results are shown in Fig. 5. All the topological transport quantities (disorder averaged) retain robust (half-)quantized values in the weak disorder regime, while they all vanish in the strong disorder regime. Furthermore, with the introduction of non-Hermiticity in the system as t 1 → √ t 2 1 -t 2 2 [Eq. (3)], fermions become weakly dispersive and the latter event occurs at weaker disorder, which can be seen by comparing the results for t 2 = 0

and 0 5. . However, even L = 2 D = 180 systems (Fig. 1) turn out to be insufficient to compare the results for different finite t 2 values. Numerical simulations in larger systems become too time expensive for us to pursue at this stage. Nevertheless, most of the results from Fig. 5 strongly suggest that with increasing t 2 , the disappearance of the topological responses occurs at weaker disorder. It can be seen by comparing the results for t 2 = 0 1 . and 0 5. Then the transport quantities deviate from their . (half-)quantized values at weaker disorder for t 2 = 0 5 . than for t 2 = 0 1, except for the . d + id paired state, for which the curves are very close to each other. The prominent finite size effects in the d + id paired state due to longer range hopping encoded in d 2 ( k ) is confirmed by numerically computing disorder averaged σ sp xy in a smaller L = 2 D = 120 system [65].

Discussions . Here we identify a family of NH operators, devoid of NH skin effects, for planar topological insulators and superconductors that manifests the topological invariant through (half-)quantized electrical, thermal, and spin transport quantities, closely mimicking the ones previously found in Hermitian systems [61, 62]. For any | t 2 | &gt; 1, when the bulk topological invariants ( C and C sp ) vanish, there is no quantization of any transport quantity. In the SM, we analytically show that the topological bound state exists only when t 2 &lt; 1 [65]. However, in order to observe vanishing transport responses we need to set t 2 to be slightly bigger than unity to bypass finite size effects [Figs. 2-4]. In model NH operators, displaying the NH skin effect, besides the topological edge modes all the left and right eigenvectors reside near the opposite edges of the scattering region by definition, ruining any (half-)quantized transport responses therein. Even though the requisite six-terminal Hall bar arrangement is well-developed by now, controlled synthesis of NH quantum crystals remains far from reality. Nevertheless, optical lattices of neutral atoms constitute a promising testbed for our theoretical predictions, where a plethora of topological band engineering protocols has been proposed [69] and realized [70]. Simplicity of our construction, in which the NH operators result from nearest-neighbor hopping modulations, makes them achievable on such highly tunable platforms that also harbor superfluids (charge-neutral superconductors). In such systems, the Hall conductivity can be obtained from the 'heating effect' [71], for example, which can also be employed to measure quantized longitudinal transports.

Acknowledgments. S.K.D. was supported by the Startup Grant of B.R. from Lehigh University. B.R. was supported by NSF CAREER Grant No. DMR-2238679. Portions of this research were conducted on Lehigh University's Research Computing infrastructure partially supported by NSF Award No. 2019035. We thank Vladimir Juriˇi´ for critical reading of the manuscript. c c

- [1] M. Z. Hasan and C. L. Kane, Colloquium: Topological insulators, Rev. Mod. Phys. 82 , 3045 (2010).
- [2] X.-L. Qi and S.-C. Zhang, Topological insulators and superconductors, Rev. Mod. Phys. 83 , 1057 (2011).
- [3] C.-K. Chiu, J. C. Y. Teo, A. P. Schnyder, and S. Ryu, Classification of topological quantum matter with symmetries, Rev. Mod. Phys. 88 , 035005 (2016).
- [4] C. L. Kane and E. J. Mele, Z 2 Topological Order and the Quantum Spin Hall Effect, Phys. Rev. Lett. 95 , 146802 (2005).
- [5] B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Quantum Spin Hall Effect and Topological Phase Transition in HgTe Quantum Wells, Science 314 , 1757 (2006).
- [6] L. Fu, C. L. Kane, and E. J. Mele, Topological Insulators in Three Dimensions, Phys. Rev. Lett. 98 , 106803 (2007).
- [7] L. Fu and C. L. Kane, Topological insulators with inversion symmetry, Phys. Rev. B 76 , 045302 (2007).
- [8] J. E. Moore and L. Balents, Topological invariants of time-reversal-invariant band structures, Phys. Rev. B 75 , 121306 (2007).
- [9] R. Roy, Topological phases and the quantum spin Hall effect in three dimensions, Phys. Rev. B 79 , 195322 (2009).
- [10] S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. Ludwig, Topological insulators and superconductors: tenfold way and dimensional hierarchy, New J. Phys. 12 , 065010 (2010).
- [11] A. Kitaev, Periodic table for topological insulators and superconductors, AIP Conf. Proc. 1134 , 22 (2009).
- [12] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, Classification of topological insulators and superconductors in three spatial dimensions, Phys. Rev. B 78 , 195125 (2008).
- [13] L. Fu, Topological Crystalline Insulators, Phys. Rev. Lett. 106 , 106802 (2011).
- [14] R.-J. Slager, A. Mesaros, V. Juriˇi´, and J. Zaanen, The c c space group classification of topological band-insulators, Nat. Phys. 9 , 98 (2012).
- [15] K. Shiozaki and M. Sato, Topology of crystalline insulators and superconductors, Phys. Rev. B 90 , 165114 (2014).
- [16] B. Bradlyn, L. Elcoro, J. Cano, M. G. Vergniory, Z. Wang, C. Felser, M. I. Aroyo, and B. A. Bernevig, Topological quantum chemistry, Nature (London) 547 , 298 (2017).
- [17] G. E. Volovik, The Universe in a Helium Droplet (Oxford University Press, Oxford, UK, 2009).
- [18] Y. Ando and L. Fu, Topological Crystalline Insulators and Topological Superconductors: From Concepts to Materials, Annu. Rev. Condens. Matter Phys. 6 , 361 (2015).
- [19] M. Sato and Y. Ando, Topological superconductors: a review, Rep. Prog. Phys. 80 , 076501 (2017).
- [20] T. Senthil, J. B. Marston, and M. P. A. Fisher, Spin quantum Hall effect in unconventional superconductors, Phys. Rev. B 60 , 4245 (1999).
- [21] N. Read and D. Green, Paired states of fermions in two dimensions with breaking of parity and time-reversal symmetries and the fractional quantum Hall effect, Phys. Rev. B 61 , 10267 (2000).
- [22] L. E. F. F. Torres, Perspective on topological states of non-Hermitian lattices, J. Phys. Mater. 3 , 014002 (2019).
- [23] A. Ghatak and T. Das, New topological invariants in

- non-Hermitian systems, J. Phys.: Condens. Matter 31 , 263001 (2019).
- [24] E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Exceptional topology of non-Hermitian systems, Rev. Mod. Phys. 93 , 015005 (2021).
- [25] K. Esaki, M. Sato, K. Hasebe, and M. Kohmoto, Edge states and topological phases in non-Hermitian systems, Phys. Rev. B 84 , 205128 (2011).
- [26] S.-D. Liang and G.-Y. Huang, Topological invariance and global Berry phase in non-Hermitian systems, Phys. Rev. A 87 , 012118 (2013).
- [27] S. Yao and Z. Wang, Edge States and Topological Invariants of Non-Hermitian Systems, Phys. Rev. Lett. 121 , 086803 (2018).
- [28] S. Yao, F. Song, and Z. Wang, Non-Hermitian Chern Bands, Phys. Rev. Lett. 121 , 136802 (2018).
- [29] Z. Gong, Y. Ashida, K. Kawabata, K. Takasan, S. Higashikawa, and M. Ueda, Topological Phases of NonHermitian Systems, Phys. Rev. X 8 , 031079 (2018).
- [30] K. Kawabata, K. Shiozaki, and M. Ueda, Anomalous helical edge states in a non-Hermitian Chern insulator, Phys. Rev. B 98 , 165148 (2018).
- [31] H. Shen, B. Zhen, and L. Fu, Topological Band Theory for Non-Hermitian Hamiltonians, Phys. Rev. Lett. 120 , 146402 (2018).
- [32] F. K. Kunst, E. Edvardsson, J. C. Budich, and E. J. Bergholtz, Biorthogonal Bulk-Boundary Correspondence in Non-Hermitian Systems, Phys. Rev. Lett. 121 , 026808 (2018).
- [33] K. Kawabata, S. Higashikawa, Z. Gong, Y. Ashida, and M. Ueda, Topological unification of time-reversal and particle-hole symmetries in non-Hermitian physics, Nat. Commun. 10 , 297 (2019).
- [34] H. Zhou and J. Y. Lee, Periodic table for topological bands with non-Hermitian symmetries, Phys. Rev. B 99 , 235112 (2019).
- [35] K. Yokomizo and S. Murakami, Non-Bloch Band Theory of Non-Hermitian Systems, Phys. Rev. Lett. 123 , 066404 (2019).
- [36] K. Kawabata, K. Shiozaki, M. Ueda, and M. Sato, Symmetry and topology in non-hermitian physics, Phys. Rev. X 9 , 041015 (2019).
- [37] C. H. Lee and R. Thomale, Anatomy of skin modes and topology in non-Hermitian systems, Phys. Rev. B 99 , 201103 (2019).
- [38] N. Okuma, K. Kawabata, K. Shiozaki, and M. Sato, Topological Origin of Non-Hermitian Skin Effects, Phys. Rev. Lett. 124 , 086801 (2020).
- [39] K. Zhang, Z. Yang, and C. Fang, Correspondence between Winding Numbers and Skin Modes in NonHermitian Systems, Phys. Rev. Lett. 125 , 126402 (2020).
- [40] D. S. Borgnia, A. J. Kruchkov, and R.-J. Slager, NonHermitian Boundary Modes and Topology, Phys. Rev. Lett. 124 , 056802 (2020).
- [41] H. C. Wu, L. Jin, and Z. Song, Nontrivial topological phase with a zero Chern number, Phys. Rev. B 102 , 035145 (2020).
- [42] X.-Q. Sun, P. Zhu, and T. L. Hughes, Geometric Response and Disclination-Induced Skin Effects in NonHermitian Systems, Phys. Rev. Lett. 127 , 066401 (2021).
- [43] W. Zhu, W. X. Teo, L. Li, and J. Gong, Delocalization of topological edge states, Phys. Rev. B 103 , 195414 (2021).
- [44] A. Panigrahi, R. Moessner, and B. Roy, Non-Hermitian dislocation modes: Stability and melting across excep-
- tional points, Phys. Rev. B 106 , L041302 (2022).
- [45] S. Manna and B. Roy, Inner skin effects on non-Hermitian topological fractals, Commun. Phys. 6 , 10 (2023).
- [46] R. Chaturvedi, V. K¨nye, E. M. Hankiewicz, J. van den o Brink, and I. Cosma Fulga, Non-Hermitian topology of transport in Chern insulators, arXiv:2406.14303.
- [47] K. v. Klitzing, G. Dorda, and M. Pepper, New Method for High-Accuracy Determination of the Fine-Structure Constant Based on Quantized Hall Resistance, Phys. Rev. Lett. 45 , 494 (1980).
- [48] M. Ko¨ onig, S. Wiedmann, C. Br¨ne, A. Roth, H. Buhu mann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Quantum Spin Hall Insulator State in HgTe Quantum Wells, Science 318 , 766 (2007).
- [49] I. Knez, R.-R. Du, and G. Sullivan, Evidence for Helical Edge Modes in Inverted InAs / GaSb Quantum Wells, Phys. Rev. Lett. 107 , 136603 (2011).
- [50] R. Yu, W. Zhang, H.-J. Zhang, S.-C. Zhang, X. Dai, and Z. Fang, Quantized Anomalous Hall Effect in Magnetic Topological Insulators, Science 329 , 61 (2010).
- [51] C.-Z. Chang, J. Zhang, X. Feng, J. Shen, Z. Zhang, M. Guo, K. Li, Y. Ou, P. Wei, L.-L. Wang, Z.-Q. Ji, Y. Feng, S. Ji, X. Chen, J. Jia, X. Dai, Z. Fang, S.-C. Zhang, K. He, Y. Wang, L. Lu, X.-C. Ma, and Q.-K. Xue, Experimental Observation of the Quantum Anomalous Hall Effect in a Magnetic Topological Insulator, Science 340 , 167 (2013).
- [52] C.-Z. Chang, W. Zhao, D. Y. Kim, H. Zhang, B. A. Assaf, D. Heiman, S.-C. Zhang, C. Liu, M. H. W. Chan, and J. S. Moodera, High-precision realization of robust quantum anomalous Hall state in a hard ferromagnetic topological insulator, Nature Mater 14 , 473 (2015).
- [53] M. Banerjee, M. Heiblum, V. Umansky, D. E. Feldman, Y. Oreg, and A. Stern, Observation of half-integer thermal Hall conductance, Nature (London) 559 , 205 (2018).
- [54] Y. Kasahara, T. Ohnishi, Y. Mizukami, O. Tanaka, S. Ma, K. Sugii, N. Kurita, H. Tanaka, J. Nasu, Y. Motome, T. Shibauchi, and Y. Matsuda, Majorana quantization and half-integer thermal quantum Hall effect in a Kitaev spin liquid, Nature (London) 559 , 227 (2018).
- [55] S. K. Srivastav, M. R. Sahu, K. Watanabe, T. Taniguchi, S. Banerjee, and A. Das, Universal quantized thermal conductance in graphene, Sci. Adv. 5 , eaaw5798 (2019).
- [56] G. Le Breton, R. Delagrange, Y. Hong, M. Garg, K. Watanabe, T. Taniguchi, R. Ribeiro-Palau, P. Roulleau, P. Roche, and F. D. Parmentier, Heat Equilibration of Integer and Fractional Quantum Hall Edge Modes in Graphene, Phys. Rev. Lett. 129 , 116803 (2022).
- [57] A. Altland and M. R. Zirnbauer, Nonstandard symmetry classes in mesoscopic normal-superconducting hybrid structures, Phys. Rev. B 55 , 1142 (1997).
- [58] D. J. Salib, S. K. Das, and B. Roy, Model nonHermitian topological operators without skin effect, arXiv:2309.12310.
- [59] C. W. Groth, M. Wimmer, A. R. Akhmerov, and X. Waintal, Kwant: A software package for quantum transport, New J. Phys. 16 , 063065 (2014).
- [60] I. C. Fulga, Y. Oreg, A. D. Mirlin, A. Stern, and D. F. Mross, Temperature Enhancement of Thermal Hall Conductance Quantization, Phys. Rev. Lett. 125 , 236802 (2020).
- [61] S. K. Das, S. Manna, and B. Roy, Topologically distinct atomic insulators, Phys. Rev. B 108 , L041301 (2023).
- [62] S. K. Das and B. Roy, Quantized thermal and spin trans-

ports of dirty planar topological superconductors, Phys. Rev. B 109 , 195403 (2024).

- [63] X.-L. Qi, Y.-S. Wu, and S.-C. Zhang, Topological quantization of the spin Hall effect in two-dimensional paramagnetic semiconductors, Phys. Rev. B 74 , 085308 (2006).
- [64] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Quantized Hall Conductance in a TwoDimensional Periodic Potential, Phys. Rev. Lett. 49 , 405 (1982).
- [65] See Supplemental Material at XXX-XXXX for the details of transport calculations, saturation of transport responses with the number of disorder realizations, calculations of electrical and thermal Hall conductivities, finite size effects on transport in d + id paired state, and analytical solution of the topological zero mode.
- [66] S. K. Das and B. Roy, Quantized electrical, thermal, and spin transports of non-Hermitian clean and dirty two-
- dimensional topological insulators and superconductors 10.5281/zenodo.13328890 (2024).
- [67] W. Long, H. Zhang, and Q.-f. Sun, Quantum thermal Hall effect in graphene, Phys. Rev. B 84 , 075416 (2011).
- [68] R. B. Laughlin, Magnetic induction of d x 2 -y 2 + id xy order in highT c superconductors, Phys. Rev. Lett. 80 , 5188 (1998).
- [69] D.-W. Zhang, Y.-Q. Zhu, Y. X. Zhao, H. Yan, and S.-L. Zhu, Topological quantum matter with cold atoms, Adv. Phys. 67 , 253 (2018).
- [70] N. R. Cooper, J. Dalibard, and I. B. Spielman, Topological bands for ultracold atoms, Rev. Mod. Phys. 91 , 015005 (2019).
- [71] L. Asteria, D. T. Tran, T. Ozawa, M. Tarnowski, B. S. Rem, N. Fl¨schner, a K. Sengstock, N. Goldman, and C. Weitenberg, Measuring quantized circular dichroism in ultracold topological matter, Nat. Phys. 15 , 449 (2019).
Sanjib Kumar Das Bitan Roy
Quantized electrical, thermal, and spin transports of non-Hermitian clean and dirty two-dimensional topological insulators and superconductors
From lattice-regularized models, devoid of any non-Hermitian (NH) skin
effects, here we compute the electrical ($\sigma_{xy}$), thermal
($\kappa_{xy}$), and spin ($\sigma^{sp}_{xy}$) Hall, and the electrical
($G_{xx}$) and thermal ($G^{th}_{xx}$) longitudinal conductivities for
appropriate NH planar topological insulators and superconductors related to all
five non-trivial Altland-Zirbauer symmetry classes in their Hermitian limits.
These models feature real eigenvalues over an extended NH parameter regime,
only where the associated topological invariants remain quantized. In this
regime, the NH quantum anomalous and spin Hall insulators show quantized
$\sigma_{xy}$ and $G_{xx}$, respectively, the NH $p+ip$ ($p \pm ip$) pairing
shows half-quantized $\kappa_{xy}$ ($G^{th}_{xx}$), while the NH $d+id$ pairing
shows quantized $\kappa_{xy}$ and $\sigma^{sp}_{xy}$ in the clean and weak
disorder (due to random pointlike charge impurities) regimes. We compute these
quantities in experimentally realizable suitable six-terminal setups using the
Kwant software package. But, in the strong disorder regime, all these
topological responses vanish and with the increasing non-Hermiticity in the
system this generic phenomenon occurs at weaker disorder.