Reciprocal skin effect and its realization in a topolectrical circuit

A system is non-Hermitian when it exchanges energy with its environment and non-reciprocal when it behaves differently upon the interchange of input and response. Within the field of metamaterial research on synthetic topological matter, the skin effect describes the conspiracy of non-Hermiticity and non-reciprocity to yield extensive anomalous localization of all eigenmodes in a (quasi) one-dimensional geometry. Here, we introduce the reciprocal skin effect, which occurs in non-Hermitian but reciprocal systems in two or more dimensions: Eigenmodes with opposite longitudinal momentum exhibit opposite transverse anomalous localization. We experimentally demonstrate the reciprocal skin effect in a passive RLC circuit, suggesting convenient alternative implementations in optical, acoustic, mechanical, and related platforms. Skin mode localization brings forth potential applications in directional and polarization detectors for electromagnetic waves.

Introduction -First developed in the late 1920s by Felix Bloch and contemporaries 1 , the electronic band theory of crystals has experienced an unprecedented revival of interest and sophistication through the theoretical and experimental discovery of topological matter [2][3][4] . Whereas band theory was previously employed in the context of infinite systems, the study of topological insulators 5,6 and semimetals 7 has reshifted attention to the interplay of bulk and boundary degrees of freedom, culminating in the celebrated bulkboundary correspondence 8,9 . It allows for the prediction of universal boundary phenomena from bulk properties and relies on the fact that changing boundary conditions does not affect the bulk states at large. The predominant framework for topological crystal matter has been closed systems at zero temperature, whose band spectrum is obtained from a Hermitian Bloch Hamiltonian.
Recently, open systems that do not conserve energy or particle number have received significant attention. Within band theory, open systems are modeled by non-Hermitian Hamiltonians. The study of topological phases supported by such Hamiltonians yields qualitatively novel phenomena 10 , such as the occurrence of stable defective degeneracy points, socalled exceptional points, in two-dimensional systems 11 , and bulk Fermi arcs connecting them 12 . Another exotic feature resolved through the band theory description of open systems is the non-Hermitian skin effect 13 : all of a one-dimensional system's eigenstates can localize at one of its boundaries for open boundary conditions (OBC), which can be studied in detail since it lends itself to exact solvability despite the presence of a boundary 14 . The skin effect constitutes a breakdown of bulk-boundary correspondence, which is not found in closed systems 14,15 .
So far, the non-Hermitian skin effect has been discussed in one-dimensional systems, and was assumed to necessitate a non-reciprocal Bloch Hamiltonian usually realized through asymmetric, direction-dependent hoppings. In this work, we present a generalized skin effect in reciprocal systems, that is enabled by non-Hermiticity in dimensions higher than one. For simplicity, we constrain ourselves to two spatial dimensions. The basic idea is that, while a reciprocal model cannot localize all of its OBC eigenstates at just one boundary, it may still do so for states of a particular boundary momentum. Reciprocity then implies that the eigenstates at opposite momenta are localized at opposite boundaries.
The reciprocal skin effect dramatically expands the scenarios in which anomalous extensive skin mode localization 14 can be found in nature. As a paradigmatic experimental realization within a framework of high accessibility and tunability, we implement the reciprocal skin effect within a topolectric circuit setup [16][17][18][19][20][21] . Similar to many other classical platforms of synthetic topological matter, this roots in the insight that Berry phase phenomena 22,23 directly carry over to classical physics 24 , since they do not concern phase space, but parameter space. The realization of topological phases in photonic 25 , mechanical 26,27 and electrical systems 16,[28][29][30] explicates this connection between topological aspects of quantum and classical systems. The non-Hermitian band theory of open quantum systems can similarly be mirrored by descriptions of classical systems that involve gain and loss. Electric circuits whose circuit Laplacian 16 is modeled in analogy to a quantum Bloch Hamiltonian are particularly suited for this: loss and gain can be directly implemented by resistors and active circuit elements 18 , respectively, while Hermitian hybridization elements relate to inductive and capacitive links between circuit nodes. We design, describe, and measure 31 (1) with periodic boundary conditions. Solid and dashed lines correspond to hopping amplitudes t and −t, respectively. Bands touch in two Dirac points. Right: π-flux model with a non-Hermitian diagonal hopping (pink) included, as defined in Eq. (2). Each Dirac point splits into a pair of exceptional points (only the real part of the band eigenvalues is plotted.) b) Real and imaginary part of eigenvalues in the vicinity of an exceptional point, at which two complex eigenvalues are degenerate. c) Schematic of the reciprocal skin effect for OBC in x and PBC in y direction: near two opposite momenta, ky = π/2 and ky = 3π/2, all eigenstates are exponentially localized with localization length ξ to the left and right of the system, respectively. At ky = 0, π the modes are completely delocalized. Right: The reciprocal skin effect could serve as a direction detector for incident electromagnetic waves: dependent on the propagation direction and polarization, a voltage will build up on the left or right edge of the system. a non-Hermitian topolectrical circuit that displays, among exceptional points, the reciprocal skin effect, and is built entirely from capacitors, inductors, and resistors, without the need for non-reciprocal active elements such as operational amplifiers. Being a passive circuit network, its convenient translation into alternative platforms of synthetic topological matter promises ubiquitous realization and application in optics, mechanics, and acoustics.
Theory of the reciprocal skin effect -To underline the generality of our arguments, we will present the theory in the language of (non-Hermitian) tight-binding Hamiltonians, which may represent either a quantum system or the dynamical matrix and response function of a classical system. We start our considerations with the π-flux tight-binding model on a square lattice. It is characterized by a nearest-neighbor hopping t, where exactly one of the four sides of each plaquette has a negative hopping amplitude compared to the three others. These hoppings require a unit cell of two plaquettes of the square lattice [see Fig. 1 a)]. The Bloch Hamiltonian for a system with periodic boundary conditions (PBC) in x and y directions, yielding the momenta k x and k y , can be written as The model has two Dirac-like band touchings at momenta (k x , k y ) = (π, π/2) and (k x , k y ) = (π, 3π/2) [see Fig. 1 a) for the band structure].
We add a non-Hermitian (gain/loss) term as a diagonal hopping through one of the plaquettes in the unit cell [see Fig. 1 a)]. It assigns a complex amplitude ir (with r a real number) to the the process of a particle hopping along the diagonal toward the upper right, and the exact same amplitude ir to the reversed process. Hermiticity would require that the latter process has the complex conjugated amplitude −ir. As such, the addition violates Hermiticity for this diagonal hopping process. Despite its non-Hermiticitiy, the model still is reciprocal, as it satisfies H(k x , k y ) = H T (−k x , −k y ) under transposition. Hamiltonian (2) has two complex-valued eigenbands with remarkable properties: they touch (i.e., the two eigenvalues have equal real and imaginary part) in two pairs of points. Upon introducing a finite r, the Dirac points of the original πflux model each split into a pair of these points. Such degeneracy points in non-Hermitian systems are called exceptional points and are the generic band touchings in a space with two parameters (here, k x and k y ). The band structure in the vicinity of an exceptional point is illustrated in Fig. 1 c).
We proceed to consider Hamiltonian (2) with OBC in xdirection, which leaves k y ∈ [0, 2π] well defined as a boundary momentum, and defineH(k y ) as the Hamiltonian for the strip geometry. WhileH (k y ) =H(−k y ) guarantees reciprocity of the full system, each instance ofH(k y ) locally breaks reciprocity for a fixed k y when seen as a purely onedimensional model, except for k y = −k y mod (2π). The effective 1D model for a given k y exhibits the non-Hermitian skin effect as discussed in Refs. 13 and 32: The eigenstates of a Hermitian system form an orthonormal basis whose squared amplitudes, when summed over all states, are equal on all lattice sites. In a non-Hermitian system this need not be the case, since right (left) eigenstates of a non-Hermitian matrix do not individually form an orthonormal basis. As a result, they can all be localized at only one edge of the system, which defines the skin effect. In our model, the skin effect is realized due to the term proportional to r in Eq. (2), which renders the hopping probability for going right different from the probability for going left. This leads to an accumulation of all eigenstates towards only one edge.
The reciprocal skin effect is characterized by a k y dependence of this localization property. For k y ∈ (0, π), all OBC bulk modes localize at the right edge, while the localization switches to the left edge for k y ∈ (π, 2π). The states are delocalized at k y = 0, π, because the HamiltonianH(k y ) is locally reciprocal. We find the strongest localization at k y ∼ π/2, 3π/2, where the initial bulk π-flux model has its gapless Dirac point. Due to the global reciprocity of the Hamiltonian, each localized eigenmode at k y has a reciprocal partner at −k y mod 2π, which localizes at the opposing edge. We conclude that, depending on the boundary conditions, the Hamiltonian displays a stark difference in the spectral and eigenstate properties: reciprocal skin effect in a strip geometry versus exceptional points and extended states in a periodic geometry. This is a manifestation of the breakdown of bulk-boundary correspondence in a non-Hermitian system that contains purely reciprocal physical couplings.
Experimental topolectrical circuit realization -In principle, various non-Hermitian, i.e., lossy, classical systems could be deliberately tailored to study the physical effects outlined above. Topolectrical circuits in particular, however, offer an important advantage due to the only mild limitation imposed by local connectivity contstraints otherwise typical to many other metamaterial settings. Most importantly, a connection between the leftmost and rightmost site of a circuit can simply be toggled on/off to change between PBC and OBC. This way, the breakdown of bulk-boundary correspondence can be studied most directly.
We experimentally implemented a circuit which realizes the non-Hermitian but reciprocal model (2) through its response function within linear circuit theory, the circuit Laplacian matrix J(ω), which takes the role of the reciprocal Hamiltonian previously introduced. It connects the input currents I a (ω) at node a of the system to the voltages V b (ω) measured at b via Kirchhoff's law, represented in the frequency domain by Here we consider a fixed excitation frequency ω, the summation over the repeated index b is implied. (See Supplemental Information for the theory on circuit Laplacians.) In a Hermitian circuit, all eigenvalues of J(ω) are purely imaginary [iJ(ω) is a Hermitian matrix], while the inclusion of resistors generates complex eigenvalues. Interpreting J(ω) at a fixed frequency ω 0 as a hopping matrix, one observes that the sign change of t necessary to implement the π-flux model is achieved by connecting one of the three bonds surrounding a plaquette in a square lattice with an inductor and three with a capacitor [see Assuming translation invariance (realized to the accuracy of the circuit element specifications), we can represent the voltages and input currents in terms of Fourier modes in reciprocal space. This leads to the k-space representation of the Laplacian as well as its voltage eigenmodes and allows for the definition of a complex-valued admittance band structure. The latter can be understood as a complex mapping from wave vector k to admittance eigenvalues of the Laplacian. In a measurement, we can decompose the voltage response to an external current excitation and find the eigensystem of the Laplacian for both periodic and open boundary conditions. Figure 2 b) (left) shows the measured bands including the two Dirac-cone band touchings expected for the π-flux model.
Unlike in experimental realizations of the non-reciprocal skin effect [33][34][35] , the reciprocal skin effect can be effected without active elements. For that, we added resistive couplings across the diagonal of every other row of plaquettes, which is analogous to the non-Hermitian term proportional to r in Eq. (2). Figure 2 b) (right) shows the measured bands. The Dirac points are broadened into branch cuts, each of which spans between a pair of exceptional points. To demonstrate this, we plot the phase of the eigenvalue for the band j n (ω 0 , k x , k y ) with smaller imaginary part along a closed path in Fig. 2 c). The observed winding and phase jump by π is direct evidence that the path encircles an exceptional point in momentum space, which is topologically stable exactly through this half-integer winding number of the band eigenvalue around it.
Having confirmed that the circuit realizes the desired physics of an exceptional point band structure with PBC, we now present measurements with OBC in x-direction to demonstrate the reciprocal skin effect. Edge terminations are chosen such that the circuit grounding does not introduce undesired off-sets due to a change in the total node conductance at the edge sites, see Fig. 2 a). The measured k y -resolved spectra are shown in Fig. 2 d), in comparison to an equivalent representation of the data for PBC. The eigenvalues with small imaginary part around k y ∼ π/2, 3π/2 indeed show the expected reorganization from a spectrum with two exceptional points in the bulk towards much fewer states with OBC -a breakdown of the Hermitian bulk-boundary correspondence. The reciprocal skin effect is encoded in the coloring of the data points: red and blue dots correspond to right and left localized eigenstates. (Note that through the measurement of the full matrix J(ω 0 ) not only do we have access to its spectrum, but also to all of its eigenstates.) Remarkably, all states near k y ∼ π/2 are right-localized, while all states near k y ∼ 3π/2 are left-localized.
Discussion -We have introduced and experimentally demonstrated the concept of a reciprocal skin effect, where the breakdown of bulk-boundary correspondence occurs in the absence of any non-reciprocal coupling. Instead of having extensive mode accumulation all along one boundary, an equal number of eigenmodes localizes along opposite boundaries, with the direction of localization tied to the momentum component parallel to the boundary. Key to their realization is the gain/loss associated with couplings across different sites, which effectively behave like non-reciprocal couplings at a fixed transverse momentum. The reciprocal skin effect can in principle exist when non-Hermitian reciprocal couplings connect different internal degrees of freedom, such as transverse momenta, sublattices, or pseudospins.
We observed the reciprocal skin effect in an electric circuit with solely passive linear circuit elements. Key to its realization are serial resistors, which are natural candidates for inter-site non-Hermitian loss, instead of active non-reciprocal elements such as operation amplifiers necessary for the nonreciprocal conventional skin effect. The breakdown of bulkboundary correspondence becomes evident by comparing the PBC system with exceptional points with the markedly different OBC case, featuring oppositely localized skin modes. Our circuit, and more generally the reciprocal skin effect, facilitates novel functionalities when coupled to electromagnetic waves. For instance, it lends itself to potential applications for polarization and direction detectors for electromagnetic waves, where differently directed or polarized input signals are substantially accumulated towards opposite directions. The reciprocal skin effect can occur in any physical system in which reciprocity can be understood as a perfect cancellation of reciprocity-breaking of two oppositely non-reciprocal subsystems. As such, when separated into these subsystems, e.g., specific momentum sectors, a reciprocal system may appear effectively non-reciprocal. While this general observation pertains to most reciprocal systems, we restrict our study to non-Hermitian systems and their individual momentum subsectors, since they allow for the definition of topological invariants as key ingredients of bulk-boundary correspondence.
To illustrate the idea, we adopt the language of Hamiltonian matrix elements acting on localized degrees of freedom (nodes). On a finite lattice of at least two dimensions, a coupling λ between two nodes a, b displaced by a distance Y transverse to the boundary gives rise to hopping elements H ab = λ e ikyY and H ba = λ e −ikyY in the Hamiltonian matrix by reciprocity defined as H(k y ) = H (−k y ). If λ were real, H ab = H * ba and the Hamiltonian would not pick up any non-Hermiticity, that is necessary for the skin effect to arise. Suppose λ is complex, i.e., λ = λ 0 e iθ , where the complex phase θ corresponds to a phase difference of impedance connections in the circuit interpretation. This type of coupling gives rise to a contribution to the Hamiltonian in the (a, b) basis. For θ = mπ, m ∈ Z, it is manifestly non-Hermitian. By regarding k y as a parameter and not a variable of the theory, we induce non-reciprocity on individual slices of k y in the term proportional to σ y . In particular, the effective non-reciprocity occurs with a nonzero coupling distance Y and its complex-valued prefactor λ. Due to the oddness of sin(k y Y ) in front of σ y , states with opposite k y experience reversed effective non-reciprocity. As a result, these eigenstates at fixed opposite k y accumulate along opposite boundaries as a consequence of the skin effect. In common reciprocal system setups such as RLC circuits, this reciprocal skin effect is usually not observed because excitation signals are designed to be almost invariably reciprocal, consisting of an equal superposition of +k y and −k y mod 2π modes. Skin mode accumulation in such systems can however be physically observed, if the input signals are polarized to have unequal ±k y mod 2π weights. This is detailed in methods Sec. F, where differently directed signals, created by suitably designed phase differences between sublattices, are demonstrated to induce voltage responses with a non-reciprocal directional preference.

B. Non-Hermitian circuit theory
In a passive circuit network of capacitors (C), inductors (L), and resistors (R), currents and voltages are linearly related through a discretized version of the Laplacian operator as a second order spatial derivative. However, energy in such a circuit system is not necessarily conserved, as it can dissipate in an irreversible heating process occuring in resistors. In contrast to purely capacitive or inductive couplings, which oppose the change of electric currents or voltages, a resistor opposes the flow of an electric current and is, as such, intrinsically non-Hermitian. To be able to analyze the present circuit setup, which involves all of the passive components R, L, C, we rely on non-Hermitian linear circuit theory.
Define V a and I a to be the voltage and external input current on node a of a circuit network. Using Kirchhoff's and Ohm's laws, we obtain a coupled system of differential equations for the circuit,İ where Γ ab , Σ ab and Λ ab are the reduced Laplacian matrices of capacitances, conductances and inverse inductances, and the summation over repeated indices is implied. The diagonal components a = b of the Laplacians are defined by including the circuit elements X a0 between node a and the ground. A Fourier transformation of Eq. (4) from the time to frequency domain results in where we defined J ab (ω) as the (grounded) circuit Laplacian 16 . Note that ω is treated as a parameter of the system which is fixed by the external AC driving frequency. A natural observable in a circuit is the impedance response Z a0 (ω), which is the ratio of the voltage at node a measured with respect to ground due to an input current I j = I 0 δ j,a that enters through a and exits through ground. Mathematically, Z a0 (ω) simply involves the inversion of Eq. (6) where J ab (ω) = n j n (ω)ψ n,a φ * n,b defines the spectral representation of the Laplacian with its right and left eigenvectors, ψ n and φ n . The frequency dependence of the Laplacian eigenvectors remains implicit. As the inverse of the Laplacian, the Green's function G ab (ω) = n j −1 n (ω)ψ n,a φ * n,b contains the voltage response to an external current excitation and fundamentally determines both the excitation pattern of individual eigenmodes and the circuit's impedance profile with frequency. The grounding impedances are given by the diagonal elements of G.
The off-diagonal components of the Green's function are accessible using an arrangement of measurements similar to that of the impedance response Z a0 (ω). We feed a current I a at node a and measure the voltage response V (a) b at all the other nodes. By repeating this for all input nodes, the Green's function can be reconstructed as In the circuit formalism, the Greens function as a direct observable contains full information on admittance eigenvalues and eigenmodes of the Laplacian, which can be extracted through numerical diagonalization. The procedure to measure the Green's function can be simplified in periodic models using spatial Fourier transform 31 (see Sec. D).
C. Theoretical circuit analysis

Derivation of the circuit Laplacian
We realize the reciprocal skin effect in an electrical circuit whose unit cell is shown in Fig. 2 a), with the conceptually important elements of each plaquette shown in Fig. 2 b). The two-point Laplacian of a capacitor is given by Under omission of the prefactor iω, the Laplacian matrix is symmetric and Hermitian. The Laplacian for an inductor takes a similar form that differs in the frequency-dependent admittance prefactor and, importantly, in the overall sign. At the resonance frequency ω 0 = 1/ √ LC of an LC resonator, the inductor effectively acts a negative capacitor, where −1/(ω 2 0 L) = −C in front of the Laplacian matrix. Therefore, one circuit plaquette of the π-flux model consists of three capacitors (representing t) and an inductor (representing −t at resonance), all of which are Hermitian and reciprocal. At nodes of sublattice A, the diagonal contributions in the Laplacian add to 4C, whereas for ω = ω 0 the capacitive and inductive contributions cancel each other in the diagonal term at sublattice B. To avoid the sublattice asymmetry contributing to a σ z term, we ground nodes of type A with an inductor L g ≈ L/4 such that the diagonal contribution to the Laplacian vanishes for both sublattices at resonance frequency. In total, the circuit Laplacian in momentum space is given by resembling the π-flux tight-binding model for ω → ω 0 , where the second term vanishes.
Additional to J π (k x , k y ), we introduce a resistor connecting nodes A and B of adjacent unit cells in y-direction [see Fig. 2 a), b)]. Its admittance representation reads The total Laplacian is given by J = J π + J r . With iω factored out, the added Laplacian J r is non-Hermitian and breaks time-reversal symmetry. For an arbitrary circuit network described by J(k y , k y ), reciprocity is defined as J (k x , k y ) = J(−k x , −k y ). Thus, a resistor is a reciprocal circuit element. However, if one considers a fixed k y slice of the model, reciprocity is broken as e iky = e −iky for arbitrary k y = 0, π. Consider OBC in x-direction, such that the y-direction preserves translational invariance and k y remains to be welldefined. We find the non-Hermitian skin effect with an extensive number of boundary localized modes, for a range of momenta k y . The emergence of the skin effect on the level of the tight-binding model is analyzed in the following.

Hermitian part of the model
We omit the unit matrix contribution to the circuit Laplacian J and interchangeably consider the tight-binding model that it induces. The Hermitian part of the model is based on the π-flux model with the Hamiltonian H π = [1 + cos(k x )] σ x + sin(k x ) σ y + 2 cos(k y ) σ z . (13) It features two Dirac band touchings, which are located on the boundary of the Brillouin zone at (k x , k y ) = (π, π/2), (π, 3π/2). Applying open boundary conditions to the π-flux model in y-direction in a cylindrical geometry results in a pair of counter-propagating modes in the boundary Brillouin zone, as the two Dirac cones are projected onto the same point. There are no topological boundary modes in this configuration. Open boundaries in x-direction preserve the separation of the Dirac crossings at k y = π/2, 3π/2. If we treat k y as a parameter instead of as a variable in Eq. (13), we end up with an effectively one-dimensional hopping model, that is form-invariant to a hopping model on a 1D chain with an additional chiral symmetry breaking mass of m = 2 cos(k y ) parametrized by k y . As the amplitudes of the intracell and intercell hoppings of the chain are identical, there will be no topological edge states for open boundary conditions and in particular no flat band joining the two projected Dirac cones in the k y boundary Brillouin zone.

Full non-Hermitian model
We now introduce a non-Hermitian extension to the π-flux model by adding a diagonal coupling across one of the two plaquettes in the unit cell and arrive at the total Hamiltonian which is related to the circuit Laplacian by H = [J(ω 0 ) − ir1]/(iω 0 C) and r = 1/(ω 0 RC). It is reciprocal, H(k x , k y ) = H(−k x , −k y ), in agreement with the fact that it was constructed out of passive elements only. We can tune the strength of the non-Hermiticity by the resistance R, which is chosen such that r = 1 in our implementation. The eigenvalues of the Hamiltonian H are complex-valued and the Dirac cones each split into two exceptional points, which are located at k x = π and k y = 1/2 arccos(r 2 /2 − 1). The exceptional points are band closings of a complex-valued band structure where the corresponding Hamiltonian is defective, i.e., the two eigenvectors coalesce to one and the matrix is a non-diagonalizable Jordan block. The band closing point for k x remains unchanged upon the added non-Hermiticity. For r = 1, the four PBC exceptional points are pinned to k y = ±π/3, ±2π/3 mod 2π.
The effectively one-dimensional model obtained by treating k y as a parameter is extended by two non-Hermitian terms resulting in with r x = r cos(k y ), r y = r sin(k y ) and m = 2 cos(k y ).
If we fix k y to one particular value, which means taking r x , r y and m as constants, H ky (k x ) breaks reciprocity resulting from the term ir y σ y . The combined breaking of reciprocity and Hermiticity in this model gives rise to the onedimensional skin effect for r y = 0 [13][14][15] .
The localization length ξ is positive, if k y ∈ (0, π) and negative if k y ∈ (π, 2π) in the boundary Brillouin zone with open boundaries along x. This leads to left edge localized modes, if k y ∈ (0, π) and to right edge localization, if k y ∈ (π, 2π). The strongest localization is found at the original positions of the Dirac crossings in the boundary Brillouin zone, k y = π/2, 3π/2. At those points, ξ vanishes for r → 1 leading to infinite localization, which is accompanied by exceptional points in the OBC spectrum of the model. The localization inherits a reciprocal symmetry from the full model as ξ(−k y ) = −ξ(k y ). A state at k y has a reciprocal partner mode at −k y mod 2π, which possesses the same absolute localization length, but is localized on the opposite edge. Combining those partners preserves reciprocity in the full model. For k y = 0 or π, the localization length diverges. The corresponding bulk eigenstates are delocalized and bulk-boundary correspondence is restored. This roots in the fact that r y = 0 at k y = 0 or π and H 0,π (k x ) is reciprocal leading to a matching of the PBC and OBC spectrum. Due to the inclusion of non-Hermitian terms, topological states can emerge in the full Hamiltonian, that have not been present in the initial Hermitian π-flux model. Our discussion closely follows Ref. 13. To study the topological aspects in the system with N unit cells, we perform a non-unitary transformation H ky = S −1 H ky S with S = diag(1, a, a, a 2 , a 2 , . . . , a N −1 , a N ) and a = (1 − ir x − r y )/(1 − ir x + r y ), where open boundary conditions are implied. The transformed Hamiltonian under application of PBCs can be rewritten as a reciprocal Su-Schrieffer-Heeger (SSH) chain with complex coefficients, where t 0 = (1 − ir x ) 2 − r 2 y ∈ C. SSH-type edge modes exist on both edges of the system, if |t 0 | < 1, which translates to a topological regime for parametrized by the momentum k y . The bulk spectrum of H ky is gapped, if either t 0 = 1 or m = 0. The gap closes at r = √ 2 for k y = π/2, which marks the topological phase transition in r. If r ≥ √ 2, there are no topological states. As the spectra of H ky and H ky are identical, the total Hamiltonian with OBCs in x-direction exhibits a bulk gap for all k y in the boundary Brillouin zone for r = √ 2. The topological phase transition does in general not occur at the PBC band closings, but is instead parametrically displaced due to the non-Hermiticity. Only for the special choice of r = 1, the two conditions coincide, resulting in topological transitions at the PBC exceptional points for k y = ±π/3, ±2π/3 mod 2π. The topological boundary modes disperse with ±m = ±2 cos(k y ) and assume a linear dispersion at the original location of the Dirac crossings in the non-Hermitian model at k y = π/2, 3π/2 in the boundary Brillouin zone. Note that the bulk is gapped at those points, while the edge states cross at zero energy. The analysis shows that SSHtype modes can coexist with skin modes in the non-Hermitian model with open boundaries along x.

D. Experimental circuit implementation
The circuit board was fitted with metallized polypropylene film capacitors (Vishay MKP1837410161G, 100nF 160V) implementing C and SMD power inductors of 33µH (Bourns SRR0604-330KL, max. 250mΩ DC resistance) for L and 8.2µH (TDK MLF2012E8R2JT000, max. 700 mΩ DC resistance) for L g . For R we use the 18 Ω resistor (Yageo RC1206FR-0718RL). The capacitors and resistors have a tolerance of 1% and were therefore not characterized. All inductances were pre-characterized with the Hameg HM8118 LCR Brige to obtain tolerances of below 1% of the nominal component values.
The measurement was conducted at the frequency of 87.25 kHz, which was identified as the resonance frequency f 0 = ω 0 /(2π). This frequency differs from the nominal value of 1/ √ LC as parasitic effects in the inductors and capacitors shift the resonance. To find the measurement frequency, a frequency sweep of the voltage at a node A and a node B due to a current excitation was done and compared to an LTSPICE simulation of the circuit including parasitic resistances. On the basis of features such as local minima and maxima in the frequency sweep the resonance frequency could be identified.
The circuit was fed by the sine wave generator of a lockin amplifier (Zurich Instruments MFLI series). We worked with an input signal with an amplitude of 900 mV. The current was measured with an Oscilloscope (Rhode & Schwarz RTM 1054) over a resistor (11 Ω). In order to measure the band structure of the circuit, we fed the current into a node via the resistor and we measured the voltage response due to the excitation at all nodes. For the voltage measurements, the same lock-in amplifier was used. We perform a Fourier transformation of the input current and voltage vector, for PBC in x and y direction and for OBC in y direction only. Using Eq. (8), we can reconstruct one column of the Green's function for all wave vectors. By repeating this procedure for all other inequivalent nodes, the complete Green's function is accessible. For PBC only the different nodes in the unit cell are inequivalent, whereas for OBC all nodes in the non-periodic direction serve as current input for a set of measurements. After the Green's function is determined, we obtain the eigenvalues and eigenvectors by diagonalizing it. For details on the measurement process see Ref. 31.

E. Representation of the Laplacian spectra in the complex plane
We present an alternative representation of the complex spectra of the circuit Laplacian, both theoretical and experimental. For this, we consider k y as a parameter and plot for each fixed k y the set of eigenvalues in the complex plane Re(j n )-Im(j n ). This generates a two-dimensional closed surface, as shown in Fig. 3 a) for PBC. As a function of k y , the eigenvalues trace out either a single circle or two circles in the Re(j n )-Im(j n ) plane. The transitions between these two topologically distinct situations are exceptional points. The planes k y = 0 and k y = π are special, because the path traced out in the complex plane as a function of k x is reciprocal on them. Therefore, it collapses into lines on which each non-end point is visited two times as k x is varied 14 .
In summary, since we have two pairs of exceptional points, the two-dimensional eigenvalue surface in k y -Re(j n )-Im(j n ) space has genus three 36 . We can plot the experimental spectra of J(ω 0 ) in a similar fashion, shown in Fig. 3 c). We observe, for PBC, a clear from slices at constant k y with a single circle and two circles in the complex plane of eigenvalues. This strongly indicates the presence of exceptional points between these two slices.
F. Response of the reciprocal skin-effect circuit to bulk perturbations A distinct feature of the reciprocal skin effect, which may lead to applications for polarization detection of electromagnetic waves, is the response to the circuit to perturbations in the bulk. From Eq. (6) we deduce that the voltage response at site a to a driving current at frequency ω is given by Since our circuit has finite resistivity, J ab (ω) has no zero eigenvalues at real frequencies, and can therefore be readily inverted. To model the experimental situation, we consider PBC in y-direction and OBC in x-direction. We then apply a minimal driving current to two bulk sites aligned along the y-direction, where we assign the current at one site a +π/2 (−π/2) phase shift with respect to the current at the other site.
We therefore expect the current to excite the circuit Laplacian eigenstates at k y = π/2 (k y = 3π/2), which are right (left) localized in x-direction due to the reciprocal skin effect. Figure 4 shows the theoretically calculated response of the circuit that is obtained by inverting the real-space version of the Hamiltonian given by Eq. (2). We find that the two driving current patterns indeed lead to nonlocal voltage responses at the far right (left) of the circuit. The reciprocal skin effect could therefore be the basis of a polarization detection device for electromagnetic radiation. Reciprocal skin effect voltage response due to a localized bulk driving current with phase shift. The sites where the current is applied are framed in black, the zoom-ins show how the relative phase is implemented. Importantly, for a phase shift of +π/2 (−π/2) we find a nonlocal voltage response at the far right (left) side of the circuit with OBC, which is absent for phase shifts 0 and π. The arrows show the direction of voltage accumulation.