Quantum dynamics of dissipative Chern insulator

<!-- metadata commented in wiki content


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>'''Quantum dynamics of dissipative Chern insulator'''</big></div>

<div id="OLE_LINK26" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Jilian Zhonga<sup> *</sup>, Xiaoyue Li<sup>a</sup></div>

''<sup>a</sup> Department of Physics, Jiangsu University, Zhenjiang 212013, People’s Republic of China''

<span id='OLE_LINK28'></span>''<sup>*</sup>Corresponding address: [mailto:zhongjilian0505@163.com zhongjilian0505@163.com]''
-->
==Abstract==

For open quantum systems, a short-time evolution is usually well described by the effective non-Hermitian Hamiltonians, while long-time dynamics requires the Lindblad master equation, in which the Liouvillian superoperators characterize the time evolution. In this paper, we constructed an open system by adding suitable gain and loss operators to the Chern insulator to investigate the time evolution of quantum states at long times by numerical simulations. Finally, we also propose a topolectrical circuits to realize the dissipative system for experimental observation. It is found that the opening and closing of the Liouvillian gap leads to different damping behaviours of the system and that the presence of non-Hermitian skin effects leads to a phenomenon of chiral damping with sharp wavefronts. Our study deepens the understanding of quantum dynamics of dissipative system.

'''Keywords''': Open quantum system, chiral damping, topolectrical circuits

==1. Introduction==

With the laboratory advances in modulating dissipation and quantum coherence,the theory of open and nonequilibrium systems has received renewed attention [1,2]. Non-Hermitian Hamiltonians have been used to describe a large number of non-conservative systems, such as classical waves with gain and loss [3-8], solids with finite quasiparticles lifetimes [9-11], and open quantum systems [12-14]. The unique features of non-Hermitian systems have been recognized in a variety of physical settings, in particular the non-Hermitian skin effect(NHSE) [15,16], where the eigenstates of the system are exponentially localized on the boundary. In recent years, the impact of NHSE has been extensively studied [17-29].

NHSE was also found in open quantum systems [30]. For open quantum systems, the non-Hermitian effective Hamiltonian describes the time evolution of the wavefunction under post-selection conditions, while the time evolution of the density matrix (without post-selection) is driven by the Liouvillian superoperator in the master equation [2,31-33]. It has been found that the Liouvillian superoperator can also exhibit non-Hermitian skin effects and that such effects can significantly affect the dynamical behaviour of the system at long times [30,34-43]. In a large class of open quantum systems, the quantum state in the long time limit converges to the steady state by algebraic damping under periodic boundary conditions and exponential damping under open boundary conditions [30].

In recent years, it has been discovered that topolectrical circuits can be used as platform to simulate the lattice systems, thus enabling the study of topological states in topolectrical circuits and gradually developing the field of topological circuitry [44-46]. Some of the early experiments and theories were extensively studied in Hermitian systems [45,47]. Since the phenomena of non-Hermitian systems are more rich than that of Hermitian systems, increasing attentions are contributed into the non-Hermitian physics, and some interesting phenomena have also been realized by topolectrical circuits [48-52].

Previous studies on open quantum dynamics and topolectrical circuits have mainly focused on one-dimensional non-Hermitian models, and relatively few studies on higher-dimensional non-Hermitian models. In this paper, we consider a two-dimensional open quantum system based on Chern insulators. Following the method developed in Song et al. [30], we study the dynamics of this system in terms of the damping matrix derived from the Liouvillian superoperator, and give a model of topolectrical circuit realization of the damping matrix based on Kirchhoff’s theory. It is found that due to the NHSE of the damping matrix, the long-time dynamics of the system under open boundary conditions is significantly different from that under periodic boundary conditions.

Our paper is organized as follows: in section 2, we briefly review the general framework on how to convert Liouvillian operators with linear jumps to non-Hermitian damping matrix. In sections 3 and 4, we compute and numerically simulate the long-time evolution of the model. In section 5, we give the circuit model of the non-Hermitian damping matrix . Finally, we conclude in section 6.

==2. General formalism of damping matrix==

An open quantum system undergoing Markovian damping satisfies the Lindblad master equation

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
|   <math>\frac{d\rho }{dt}=-i[H,\rho ]+\sum \left(2L_{\mu }\rho L_{\mu }^{\dagger }-\right. </math><math>\left. \lbrace L_{\mu }^{\dagger }L_{\mu },\rho \rbrace \right),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}

where  <math>\rho </math> is the density matrix of the system,  <math>H</math> is the Hamiltonian that represents unitary evolution of the system, and  <math>L_{\mu }</math> are Lindblad dissipation operators describing the quantum jumps induced by the coupling to the environment. The above equation can be abbreviated as  <math>\frac{d\rho }{dt}=L\rho </math>, where  <math>L</math> is called the Liouvillian superoperator. By regarding the density matrix  <math>\rho </math> as a vector that consists of matrix elements  <math>{\rho }_{i,j}</math>,  <math>L</math> is represented as a matrix whose elements are given by [53]

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>L_{ij,kl}=\sum_{\mu }2L_{\mu ;i,k}L_{\mu ;l,j}^{\dagger }-i{\left(H-i\sum_{\mu }L_{\mu }^{\dagger }L_{\mu }\right)}_{i,k}{\delta }_{l,j}+i{\left(H+i\sum_{\mu }L_{\mu }^{\dagger }L_{\mu }\right)}_{l,j}{\delta }_{ik.}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}


These representations enable one to treat the Lindblad equation as a linear equation. In other words, the dynamics of the system can be understood in terms of the eigenvalue problem of the Liouvillian matrix: <math>L{\rho }^{\left(i\right)}={\lambda }_i{\rho }^{\left(i\right)}.</math> The Hamiltonian and dissipators can be expressed in terms of 2n Majorana fermions [54]

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>H=\sum_{i,j=1}^{2n}{\gamma }_iH_{{}_{ij}}^M{\gamma }_j,L_{\mu }=\sum_{i=1}^{2n}l_{{}_{\mu ,i}}^M{\gamma }_i</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

where  <math>{\gamma }_i</math> are Majorana fermions satisfying  <math>\lbrace {\gamma }_i,{\gamma }_j\rbrace =2{\delta }_{ij}</math>. The matrix  <math>H^M</math> is chosen to be an antisymmetric matrix, <math>{\left(H^M\right)}^T=-H^M</math>. Defining  <math>M_{ij}={\sum }_{\mu }l_{\mu ,i}^\ast l_{\mu ,j}^{}</math>,  <math>M_{{}_{ij}}^M={\sum }_{\mu }{\left(l_{\mu ,i}^M\right)}^\ast l_{\mu ,j}^M</math>, we have  <math>M_{{}_{}}^M=\frac{1}{4}M\otimes (1+{\sigma }_y)</math>. Under the third quantization [54,55], the Liouvillian superoperator is expressed as a quadratic form of the 2n complex fermions (4n Majorana fermions)

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>L\mbox{=}\frac{2}{i}\left(\begin{array}{cc}
c^{\dagger } & c
\end{array}\right)\left(\begin{array}{cc}
-Z^T & Y\\
0 & Z
\end{array}\right)\left(\begin{array}{c}
c\\
c^{\dagger }
\end{array}\right),</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}

where <math>Z=H^M+iRe{\left(M^M\right)}^T,Y=2Im{\left(M^M\right)}^T</math>, and  <math>c=(c_1,c_2,...,c_{2n})</math> are third quantized complex fermions. Through the above expression, we can obtain the Liouvillian eigenspectrum [54,55]

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\lambda ={\sum }_iE_iv_i</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}

with  <math>v_i\in \lbrace 0,1\rbrace </math>, where  <math>\left\{E_i\right\}</math> is the eigenspectrum of  <math>4iZ</math>. Here  <math>\lambda </math> contains valuable information of the full density-matrix dynamics, and it can be easily obtained from the damping matrix <math>X</math> with   <math>X_{ij}=ih_{ji}-{\sum }_{\mu }l_{\mu ,j}^\ast l_{\mu ,i}^{}</math> [36]. Rewriting <math>M</math> as  <math>M=M_r+iM_i</math>, where  <math>M_r,M_i</math> are real matrices, we have  <math>M^M=\frac{1}{4}(M_r+iM_i)\otimes (1+{\sigma }_y)</math>.   <math>Z</math> can be further written as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{align}
Z=&\frac{1}{4}(h_r\otimes {\sigma }_y+ih_i\otimes 1)+i\frac{1}{4}(M_r^T\otimes 1-iM_i^T\otimes {\sigma }_y)\\
=&\frac{1}{4}(h_r+M_i^T)\otimes {\sigma }_y+i\frac{1}{4}(h_i+M_r^T)\otimes 1
\end{align}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}


Therefore,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
det(4iZ-\lambda E)=det\left(\begin{array}{cc}
-(h_i+M_r^T)-\lambda  & h_r+M_i^T\\
-(h_r+M_i^T) & -(h_i+M_r^T)-\lambda 
\end{array}\right)\\
=det(X-\lambda E)det(X^\ast -\lambda E)
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}


The eigenvalue of  <math>4iZ</math> are the union of the eigenvalues of  <math>X</math> and  <math>X^\ast </math>, which gives the Liouvillian eigenspectrum.

Then we outline the general form of the Lindblad damping matrix in open quantum systems [30]. We consider tight-binding models whose Hamiltonian can generally be written as  <math>H=\sum {}_{ij}h_{ij}c_i^{\dagger }c_j</math>, where  <math>c_i^{\dagger },c_i</math> are the creation and annihilation operators on lattice site  <math>i</math>, and  <math>h_{ij}=h_{ij}^\ast </math> is the hopping amplitude between the lattice points of the system (<math>i\not =j</math>) or onsite potential (<math>i=j</math>). It is convenient to define the single-particle correlation function  <math>{\Delta }_{ij}(t)=Tr[c_i^{\dagger }c_j\rho (t)]</math> to observe the time evolution of the density matrix. Each cell is coupled to the environment through the gain jump operator  <math>L_{\mu }^g=\sum {}_iD_{\mu i}^gc_i^{\dagger }</math> and loss jump operator <math>L_{\mu }^l=\sum {}_iD_{\mu j}^lc_i</math>. Substituting the Lindblad quantum master equation into the time evolution of the single-particle correlation function, we can obtained

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\frac{d\Delta (\mbox{t})}{dt}=X\Delta (\mbox{t})+\Delta (\mbox{t})X^{\dagger }+2M_g,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}

where  <math>X=ih^T-(M_l^T+M_g)</math> is the damping matrix with  <math>{\left(M_g\right)}_i{}_j=\sum {}_{\mu }D_{\mu i}^{g\ast }D_{\mu j}^g</math> and  <math>{\left(M_l\right)}_i{}_j=\sum {}_{\mu }D_{\mu i}^{l\ast }D_{\mu j}^l</math>. The steady state correlation  <math>{\Delta }_s=\Delta (\infty )</math>, to which the long-time evolution of any initial state converges, is determined by  <math>d{\Delta }_s/dt=0</math> or  <math>X{\Delta }_s+{\Delta }_sX^{\dagger }+2M_g=0</math>.  Focusing on the deviation towards the steady state  <math>\tilde{\Delta }(t)=\Delta (t)-{\Delta }_s</math>, whose time evolution is  <math>d\tilde{\Delta }(t)/dt=X\tilde{\Delta }(t)+\tilde{\Delta }(t)X^{\dagger }</math>, we can integrate it with Eq. (1) to obtain

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\tilde{\Delta }(\mbox{t})=e^{Xt}\tilde{\Delta }(\mbox{0})e^{X\dagger t}.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}


Therefore, the dynamical behaviour of the system can be characterized by the damping matrix.

==3. Model==

In this paper, we consider the Chern insulator model with the Hamiltonian in momentum space as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>h(k)=l_xsink_x{\sigma }_x+l_ysink_y{\sigma }_y+{\epsilon }_k{\sigma }_z,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
|}

where  <math>{\epsilon }_k=m+t_xcosk_x+t_ycosk_y</math>. Let each unit cell contain a single loss and gain dissipator,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{align}
L_x^l= & \frac{\sqrt{2\gamma }}{2}\left(e^{-\displaystyle\frac{\pi }{4}i}c_{xA}+e^{\displaystyle\frac{\pi }{4}i}c_{xB}\right)\\
L_x^g= & \frac{\sqrt{2\gamma }}{2}\left(e^{\displaystyle\frac{\pi }{4}i}c_{xA}^{\dagger }+e^{-\displaystyle\frac{\pi }{4}i}c_{xB}\right),
\end{align}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (11)
|}

where  <math>x</math> denotes the lattice site,  <math>A,B</math> refer to the sublattice. The Fourier transformation of  <math>X</math> is  <math>X(k)=ih^T(-k)-M_l^T(-k)-M_g(k)</math>. The gain and loss dissipators are intra-cell, so these  <math>M(k)</math> matrices are independent of  <math>k</math>,  <math>M_l(k)=\frac{\sqrt{2}}{2}\lambda +\frac{1}{2}{\sigma }_x-</math><math>\frac{1}{2}{\sigma }_y,M_g(k)=\frac{\sqrt{2}}{2}\lambda +</math><math>\frac{1}{2}{\sigma }_x+\frac{1}{2}{\sigma }_y</math>. Then, the damping matrix in momentum space is

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>X(k)=i[l_xsink_x{\sigma }_x+l_ysink_y{\sigma }_y+{\epsilon }_k{\sigma }_z+</math><math>i[\lambda {\sigma }_x+\lambda {\sigma }_y]-\sqrt{2}\lambda .</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
|}


It can be written in the form of left and right eigenvectors,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>X=\sum_n{\lambda }_n\vert u_{Rn}\rangle \langle u_{Ln}\vert ,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)
|}

where  <math>X^{\dagger }\vert u_{Ln}\rangle ={\lambda }_n^\ast \vert u_{Ln}\rangle ,X\vert u_{Rn}\rangle =</math><math>{\lambda }_n\vert u_{Rn}\rangle</math>. It is worth noting that our  <math>M_l</math> and   <math>M_g</math> satisfy  <math>M_l^T+M_g=2M_g</math>, guaranteeing that  <math>{\Delta }_S=\frac{1}{2}I_{2L\times 2L}</math> is a steady state solution, where  <math>L=N_x\times N_y</math>, <math display="inline"> L </math> is the system size, and  <math>N_x,N_y</math> are the size in <math display="inline"> x,y </math> direction, respectively. We assume that the initial state of the system is the completely filled state, i.e., <math>\Delta (0)</math> is an identity matrix. Therefore, Eq.(9) can be re-expressed as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{aligned}
\tilde{\Delta}(\mathrm{t})& =\frac{1}{2}\displaystyle\sum_{n,n^{'},l}\exp[(\lambda_{n}+\lambda_{n^{'}}^{*})t]\widetilde{u_{R}}(i,n)\widetilde{u_{L}}(l,n)\widetilde{u_{L}^{*}}(l,n^{'})\widetilde{u_{R}^{*}}(j,n^{'})  \\
&=\displaystyle\frac{1}{2}\displaystyle\sum_{n,n^{'}}\displaystyle\frac{\displaystyle\sum_{l}\exp[(\lambda_{n}+\lambda_{n^{'}}^{*})t]u_{R}(i,n)u_{L}(l,n)u_{L}^{*}(l,n^{'})u_{R}^{*}(j,n^{'})}{\displaystyle\sum_{k}u_{R}(k,n)u_{L}(k,n)\displaystyle\sum_{m}u_{L}^{*}(m,n^{'})u_{R}^{*}(m,n^{'})}
\end{aligned}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (14)
|}


According to the dissipative property,  <math>Re\lbrace {\lambda }_n\rbrace \leq 0</math> always holds. The Liouvillian gap  <math>\Lambda =min[2Re(-{\lambda }_n)]</math> plays a decisive role in long-time dynamics. The opening gap (<math>\Lambda \not =0</math>) implies an exponential rate of convergence to the steady state, while the closing gap (<math>\Lambda =0</math>) implies algebraic convergence [34].

==4. chiraldamping==

For simplicity, the parameters of our model are taken as  <math>l_x=l_y=1,t_x=t_y=-1.</math> We first study the dynamical behaviour under the periodic boundary conditions. Diagonalizing  <math>X(k)</math> , we obtain the energy spectrum as shown in Fig. 1. It is found that the Liouvillian gap vanishes at  <math>m=1.5</math> , while the gap opens at  <math>m=2.5</math> . So we expect the damping rate to be algebraic and exponential in each case, respectively.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Zhong_847600978-image84.png|600px]] </div>

FIG. 1. The eigenvalues of the damping matrix X. Blue: periodic boundary; Red: open boundary . The Liouvillian gap under periodic boundary condition vanishes for (a) and (b) while it is nonzero for (c) and (d). Under open boundary condition, the Liouvillian gap is nonzero in all four cases. This significant difference between open and periodic boundary comes from the NHSE of X. (a) <math>\lambda =0.1,m=1.5</math> ;(b) <math>\lambda =0.5,m=1.5</math> ; (c)  <math>\lambda =0.1,m=2.5</math> ; (d) <math>\lambda =0.5,</math>  <math>m=2.5</math> .

To verify this, we define the site-averaged fermion number deviation from the steady state  <math>R(t)=\sqrt{\frac{1}{N_XN_Y}{\sum_xR_x(t)}^2}=\sqrt{\frac{1}{N_XN_Y}{\sum_x\left(\frac{n_x(t)-n_x(t-{\delta }_t)}{{\delta }_t}\right)}^2},</math> where  <math>R_x(t)=n_x(t)-n_x(\infty ),</math>  <math>n_x(t)={\Delta }_{xA,xA}(t)+{\Delta }_{xB,xB}(t)</math> . The numerical results are shown in Fig. 2. As anticipated, it is observed that the damping of R(t) is algebraic for cases black and red lines with m=1.5 while exponential for blue and green lines with m=2.5 under the periodic boundary condition.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Zhong_847600978-image93.png|348px]] </div>

FIG. 2. The damping of site-averaged fermion number towards the steady state under periodic boundary condition with size  <math>L=30\times 30</math> . m = 1.5 (black and red) exhibits a slow algebraic damping while m = 2.5 (blue and green) is an exponential damping. The  initial state is the completely filled state.

Next we turn to the open boundary conditions. Since the damping matrix X has  NHSE, its energy spectrum is no longer that of the periodic boundary conditions. At this point all the energy spectrums have a non-zero energy gap (red part of Fig. 1) ,therefore, we expect an exponential long-time damping of  <math>\tilde{\Delta }(t)</math> . The numerical simulation in Fig. 3 confirms this exponential behaviour with R(t) having a period of algebraic damping before entering into the exponential damping. The time of the algebraic damping increases with the size L[Fig.3(a)]. To better understand this feature, we plot the damping in several unit cells in the same x dimension (ix = 1), as shown in Fig. 3(b). It can be seen that the left end ( iy = 1) enters the exponential damping immediately, and the other sites enter the exponential damping in turn according to their different distances to the left end.Due to a process of algebraic damping that occurs before entering the exponential stage,,there is a "damping wavefront" from left ( <math>ix=1,iy=1</math> ) to right ( <math>ix=1,iy=N_y</math> ). This phenomenon is known as "chiral damping".

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Zhong_847600978-image98.png|600px]] </div>

FIG. 3. (a) The site-averaged particle number damping under periodic boundary conditions (solid line) and open boundary conditions (dashed line) for several sizes L. The long-time damping of R(t) follows a power law under periodic boundary condition, while the damping follows an exponential law after an initial power law stage under open boundary condition. (b) The particle number damping on several sites. The system size is 30×30, and the left end ( <math>ix=1,iy=1</math> ) enters the exponential phase from the beginning, followed by the other sites in turn. For (a) and (b), the initial state is completely filled state. <math>m=1.5,\lambda =0.1</math> .

The phenomenon of chiral damping can be observed more intuitively as shown in Fig. 4(a) where the colour shades indicate the value of R(t). Under the periodic boundary condition, the time evolution follows a slow power law while under the open boundary condition, a wavefront moving to the upper right is observed. This can be intuitively linked to the phenomenon that all eigenstates of X are localized in the upper right corner, which arises from the non-Hermitian skin effect of the damping matrix X. If the matrix X does not have NHSE under the open boundary condition, the fermion number of the system should have a similar behaviour of damping under different boundary conditions. Therefore, the non-Hermitian skin effect plays an important role in open quantum systems and significantly affects the dynamical behaviour of open quantum systems.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Zhong_847600978-image101-c.png|456px]] </div>

FIG. 4. The evolution of R(t) at each lattice site under open boundary conditions (a) and periodic boundary conditions (b).

==5.Experiment realized==

Next we give the scheme of topolectrical circuits to simulate the damping matrix. Based on the similarity between the Kirchhoff equation and the Schrödinger equation, it is possible to simulate the Hamiltonian of the system using different circuit components, and the different parameters in the Hamiltonian can be adjusted independently by various components. The circuit Laplacian corresponding to the Hamiltonian can be written as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>J=D-C+W,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (15)
|}


where W and D are diagonal matrices containing the total conductance from each node to the ground and to the rest of the circuit, respectively. C is the adjacency matrix of conductances [44].

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Zhong_847600978-image103.png|600px]] </div>

FIG. 5. The structure of topolectrical circuit to realize damping matrix under periodic boundary conditions.  (a) Connection relations between the nodes. The blue solid line box containing two “sublattice” nodes A (red) and B (blue) simulates a unit cell of X. The black(grey) solid line indicates the coupling between nodes in the x(y)-direction. (b)The circuit element structure is detailed for the green dashed framed rectangle in (a).    (c) Internal circuit diagram of the INIC element, consisting of an operational amplifier and impedances  <math>Z_{\pm }</math> . The impedance  <math>Z</math> is the target element, and different conductance in different directions of  <math>V_{i,r}</math> can be achieved by connecting the INIC in series.  <math>Z_{\pm }</math> satisfies  <math>Z_+=Z_{-}</math> . (d) Grounding module of the nodes. The resistances <math>R_{A,}R_B</math> and capacitance C are used to simulate the onsite potential, and inductance L allows the Laplacian eigenvalue spectrum to be shifted uniformly as desired.

Considering the periodic boundary conditions first, the topolectrical circuit for realizing the damping matrix X is illustrated in Fig.5.  Fig. 5(a) depicts the schematic diagram of the overall circuit structure, which gives the connection relationship between the nodes. Fig. 5(b) shows the detailed circuit component of the unit which is the green dashed box in Fig. 5(a). The blue box in Fig. 5(a) represents a unit cell in the system, and the two nodes inside it correspond to sublattices A (red) and B (blue). The circuit connections in the x and y directions are distinguished by black and gray. From Fig. 5(b) we can obtain the matrices C and D in Eq. (9), so that

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
D-C=-i\omega \left(\begin{array}{cc}
-i\frac{2}{\omega R_2}cosk_x-i\frac{2}{\omega R_4}cosk_y-C_1-i\frac{1}{\omega R_0} & i2C_2sink_x+\frac{2}{\omega R_6}sink_y-i\frac{1}{\omega R_1}+C_1\\
i2C_2sink_x+\frac{2}{\omega R_5}sink_y+i\frac{1}{\omega R_1}+C_1 & -i\frac{2}{\omega R_3}cosk_x-i\frac{2}{\omega R_7}cosk_y-C_1+i\frac{1}{\omega R_0}
\end{array}\right),\\
\mbox{    }
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (16)
|}


with <math>C_1=-\lambda ,C_2=l_x/2,\frac{1}{R_0}=\frac{1}{R_1}-</math><math>2(\frac{1}{R_3}+\frac{1}{R_4}),R_1=-1/\left(\omega \lambda \right),R_2=</math><math>-R_3=-2/t_x,</math>  <math>R_4=-R_7=-2/t_y</math> ， <math>R_5=-R_6=-2/l_y.</math>

Comparing it with the damping matrix, we need to add grounding elements to match the onsite potential. The grounding elements of nodes A and B are shown in  Fig. 5(d), where the resistors  <math>R_{A,}R_B</math> and capacitors  <math>C</math> simulate the lattice potential , and  <math>R_{A,}R_B</math> satisfies  <math>R_A=-R_B</math> . So the diagonal matrix W is

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>W=\left(\begin{array}{cc}
i\omega C+\frac{1}{R_A}+\frac{1}{i\omega L} & \\
 & i\omega C+\frac{1}{R_B}+\frac{1}{i\omega L}
\end{array}\right).</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (17)
|}


From Eq. (10) we get the conductance matrix of the circuit of Fig. 5(a) at  <math>\omega </math> frequency

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
J(\omega )=-i\omega \left(\begin{array}{cc}
-i\frac{2}{\omega R_2}cosk_x-i\frac{2}{\omega R_4}cosk_y-(C+C_1)-i(\frac{1}{\omega R_0}+\frac{1}{\omega R_A}) & i2C_2sink_x+\frac{2}{\omega R_6}sink_y-i\frac{1}{\omega R_1}+C_1\\
i2C_2sink_x+\frac{2}{\omega R_5}sink_y+i\frac{1}{\omega R_1}+C_1 & -i\frac{2}{\omega R_3}cosk_x-i\frac{2}{\omega R_7}cosk_y-(C+C_1)+i(\frac{1}{\omega R_0}+\frac{1}{\omega R_A})
\end{array}\right)+\frac{1}{i\omega L}\epsilon \\
=-i\omega J_P+\frac{1}{i\omega L}\epsilon ,
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (18)
|}


where  <math>C=(1-\sqrt{2})\lambda </math> ， <math>\frac{1}{R_A}=-\frac{1}{R_B}=-\omega m-\frac{1}{R_0}.</math> Comparing this Laplacian matrix with the damping matrix, the mapping relationship can be established by  <math>J_P\Leftrightarrow X</math> .

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Zhong_847600978-image123-c.png|600px]] </div>

FIG. 6. Negative impedance module [22]. (a) A single-port circuit to ground. The input impedance is  <math>Z_g=-Z</math> ; (b) Free-port circuit. Its input impedance at both ends is  <math>Z_{ij}=Z_{ji}=-Z</math> .  The markings on the ideal amplifier indicate the output voltage versus the input voltage.

Notice that the circuit requires a negative component,which is implemented as shown in Fig.6. Fig.6(a) and (b) show the equivalent negative impedance modules for a single port to ground and a free two-terminal port, respectively. They achieve the equivalent negative impedance through an amplifier. According to Kirchhoff's law, the input impedance of the single-port circuit to ground [Fig. 6(a)] can be obtained as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>Z_g=\frac{V_g}{I_g}=\frac{V_g}{(V_g-2V_g)/Z}=-Z.</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (19)
|}


The input impedance at both ends of the free port circuit [Fig. 6(b)] are

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\begin{array}{c}
Z_{ij}=\frac{V_i-V_j}{I_i}=\frac{V_i-V_j}{[V_i-2(V_i-V_j)]/Z}=-Z,\\
Z_{ji}=\frac{V_j-V_i}{I_j}=\frac{V_j-V_i}{[V_j-2(V_j-V_i)]/Z}=-Z.
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (20)
|}


That is  <math>Z_{ij}=Z_{ji}=-Z</math> .

Under the open boundary condition, the hopping amplitude of the cells located at the boundary weakens, leading to fewer branches connected to the boundary nodes in the circuit model, as shown in Fig.7(a).  Fig.7(a) gives the connection relationship between the nodes of the circuit under the open boundary condition, and the circuit nodes can be classified into body nodes (in the black dashed box), edge nodes (in yellow) and corner nodes (in green). Changes in the branch circuit of the nodes at the boundary will cause variations of the matrices C and D. The matrix C corresponds to the hopping amplitude between the lattice points, which is allowed to change. Whereas the change of D is not desired due to the same onsite potential under different boundary condition.

Therefore, we need to design specific grounding elements to eliminate the effects of variations in D. Owing to the asymmetry of the coupling strengths under periodic boundary condition, the types of the edge and corner nodes are different for each of the four orientations, so there are a total of 16 different grounding modules, as shown in Fig. 7(b). The additional grounding elements keep the diagonal matrix D+W unchanged, i.e., the onsite potential is unchanged, which achieves the mapping of the circuit Laplacian in Fig. 7 to the damping matrix X under the open boundary condition.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Zhong_847600978-image129-c.png|600px]] </div>

FIG. 7. Schematic diagram of the circuit of the damping matrix X under open boundary conditions. (a) Schematic diagram of the connection relations among the nodes. The black, yellow and green dashed boxes correspond to the body, edge and corner nodes, respectively. The circuit connections of the body node are the same as those of the periodic boundary, while the edge and corner nodes require additional grounding elements to regulate the onsite potential. (b) Grounding modules for edge and corner nodes . The grounding elements for the edge and corner nodes are different for each of the four orientations, where the negative impedance elements can be realized by Fig.6(a). Note that in addition to these grounding elements, all nodes need to be connected to the elements in Fig. 5(d).

==6. Conclusion==

In summary, we study the dynamical properties of a two-dimensional open system. The open quantum system is constructed by introducing appropriate gain and loss to the Chern insulator, and then using the damping matrix derived from the Liouvillian superoperator explore its long-time evolution. It is found that  the site-averaged fermion number deviation from the steady state under periodic boundary conditions shows a slow algebraic damping when the energy gap closes and an exponential damping when the energy gap opens. Under open boundary conditions, due to the non-Hermitian skin effect of the damping matrix, the system exhibits the phenomenon of chiral damping that the fermion number at each site undergoes a period of algebraic damping before entering an exponential damping, and the transition time that is proportional to the distance from that site to the boundary. Finally, we map the damping matrix in terms of the circuit Laplacian to give a model diagram of the topolectrical circuit implementation of the system.

==References==

[1] E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Exceptional topology of non-Hermitian systems,Rev. Mod. Phys. 2021, 93, 015005.

[2] Y. Ashida, Z. Gong, and M. Ueda, Non-Hermitian physics, Adv. Phys. 2020, 69, 249-435, 3.

[3] X. Zhu, H. Ramezani, C. Shi, J. Zhu, and X. Zhang, PT- symmetric acoustics, Phys. Rev. X 2014, 4, 031042.

[4] B. I. Popa and S. A. Cummer, Non-reciprocal and highly nonlinear active acoustic metamaterials, Nat. Commun. 2014, 5, 3398.

[5] A. Regensburger, C. Bersch, M.-A. Miri, G. Onishchukov, D. N. Christodoulides, and U. Peschel, Parity{time synthetic photonic lattices, Nature (London) 2012, 488, 167-171.

[6] L. Feng, Z. J. Wong, R.-M. Ma, Y. Wang, and X. Zhang, Single-mode laser by parity-time symmetry breaking, Science 2014, 346, 972-975.

[7] H. Zhou, C. Peng, Y. Yoon, C. W. Hsu, K. A. Nelson, L. Fu, J. D. Joannopoulos, M. Soljacic, and B. Zhen, Observation of bulk Fermi arc and polarization half charge from paired exceptional points, Science 2018, 359, 1009-1012.

[8] A. Cerjan, S. Huang, M. Wang, K. P. Chen, Y. Chong, and M. C. Rechtsman, Experimental realization of a Weyl exceptional ring, Nat. Photonics 2019,13, 623-628.

[9] Micha‛l Papaj, Hiroki Isobe, and Liang Fu, “Nodal arc of disordered dirac fermions and non-hermitian band theory,” Phys. Rev. B 2019, 99, 201107.

[10] Huitao Shen and Liang Fu, “Quantum oscillation from in-gap states and a non-hermitian landau level problem,” Phys. Rev. Lett. 2018, 121, 026403.

[11] Y. Cao, Y. Li, and X. Yang, “Non-hermitian bulkboundary correspondence in a periodically driven system,” Physical Review B, 2021, 103, 075126

[12] Ingrid Rotter, “A non-hermitian hamilton operator and the physics of open quantum systems,” Journal of Physics A: Mathematical and Theoretical 2009, 42, 153001.

[13] L. Xiao, T. Deng, K. Wang, G. Zhu, Z. Wang, W. Yi, and P. Xue, Observation of non-Hermitian bulk-boundary correspondence in quantum dynamics, Nat. Phys. 2020, 16, 761-768.

[14] L. Xiao, T. Deng, K. Wang, Zhong Wang, Wei Yi, and Peng Xue, Observation of non-Bloch parity-time symmetry and exceptional points, Phys. Rev. Lett.  2021, 126, 230402.

[15] Shunyu Yao and Zhong Wang, “Edge states and topological invariants of non-hermitian systems,”Phys. Rev. Lett.  2018, 121, 086803.

[16] S. Yao, F. Song, and Z. Wang, “Non-hermitian chern bands,” Physical review letters, 2018, 121, 136802.

[17] Ching Hua Lee and Ronny Thomale, “Anatomy of skin modes and topology in non-hermitian systems,” Phys. Rev. B 2019, 99, 201103.

[18] Huaiqiang Wang, Jiawei Ruan, and Haijun Zhang, “Non-hermitian nodal-line semimetals with an anomalous bulk-boundary correspondence,”Phys. Rev. B 2019, 99, 075130.

[19] Ching Hua Lee, Linhu Li, and Jiangbin Gong, “Hybrid higher-order skin-topological modes in nonreciprocal systems,” Phys. Rev. Lett. 2019, 123, 016805.

[20] Flore K. Kunst and Vatsal Dwivedi, “Non-hermitian systems and topology: A transfer-matrix perspective,” Phys. Rev. B 2019, 99, 245116.

[21] Borgnia D S , Kruchkov A J , Slager R J .Non-Hermitian Boundary Modes, PHYS. REV. LETT. 2019, 124, 056802.

[22] B. X. Wang and C. Y. Zhao, “Topological phonon polaritons in one-dimensional nonhermitian silicon carbide nanoparticle chains,” Phys. Rev. B 2018, 98, 165435.

[23] Y. Fu, J. Hu and S. Wan, "Non-Hermitian second-order skin and topological modes." Physical Review B 2021, 103(4), 045420.

[24] Li, Yang, et al. "Topological energy braiding of non-Bloch bands." Physical Review B 2022,10, 195425.

[25] Motohiko Ezawa, “Braiding of majorana-like corner states in electric circuits and its non-hermitian generalization,” Phys. Rev. B 2019, 100, 045407.

[26] Yang Xiaosen, Yang Cao, and Yunjia Zhai. "Non-Hermitian Weyl semimetals: Non-Hermitian skin effect and non-Bloch bulk–boundary correspondence." Chinese Physics B 2022, 31, 010308.

[27] Zi-Yong Ge, Yu-Ran Zhang, Tao Liu, Si-Wen Li, Heng Fan, and Franco Nori, “Topological band theory for non-hermitian systems from the dirac equation,” Phys. Rev. B 2019, 100, 054105.

[28] X. Ji and X. Yang, "Generalized bulk-boundary correspondence in periodically driven non-Hermitian systems." Journal of Physics: Condensed Matter 2024, 36, 243001.

[29] Li, Yang, et al. "Universal characteristics of one-dimensional non-Hermitian superconductors." Journal of Physics: Condensed Matter 2022, 35, 055401.

[30] Song, Fei, Shunyu Yao, and Zhong Wang. "Non-Hermitian skin effect and chiral damping in open quantum systems." Physical review letters 2019,123, 170401.

[31] J. Dalibard, Y. Castin, and K. Mølmer, Wave-Function Approach to Dissipative Processes in Quantum Optics, Phys. Rev. Lett. 1992, 68, 580.

[32] H. J. Carmichael, Quantum Trajectory Theory for Cascaded Open Systems, Phys. Rev. Lett. 1993, 70, 2273.

[33] H. Weimer, A. Kshetrimayum, and R. Or´us, Simulation methods for open quantum many-body systems, Rev. Mod. Phys. 2021, 93, 015008.

[34] Cai, Zi, and Thomas Barthel. "Algebraic versus exponential decoherence in dissipative many-particle systems." Physical review letters  2013, 111, 150403.

[35] Zhou, Ziheng, and Zhenhua Yu. "Non-Hermitian skin effect in quadratic Lindbladian systems: An adjoint fermion approach." Physical Review A 2022, 106, 032216.

[36] Li, Tianyu, Yong-Sheng Zhang, and Wei Yi. "Engineering dissipative quasicrystals." Physical Review B 202, 105, 125111.

[37] He, Peng, et al. "Damping transition in an open generalized Aubry-André-Harper model." Physical Review A  2022,105, 023311.

[38] Liu, Chun-Hui, et al. "Helical damping and dynamical critical skin effect in open quantum systems." Physical Review Research 2020, 2, 043167.

[39] Xi-Wang Luo and Chuanwei Zhang, “Higher-order topological corner states induced by gain and loss,”Phys. Rev. Lett. 2019, 123, 073601.

[40] Manzano, Daniel. "A short introduction to the Lindblad master equation." Aip Advances  2020,10, 025106.

[41] Haga, Taiki, et al. "Liouvillian skin effect: Slowing down of relaxation processes without gap closing." Physical Review Letters  2021, 127, 070402.

[42] McDonald, Alexander, and Aashish A. Clerk. "Exponentially-enhanced quantum sensing with non-Hermitian lattice dynamics." Nature communications  2020, 11, 5382.

[43] McDonald, Alexander, Ryo Hanai, and Aashish A. Clerk. "Nonequilibrium stationary states of quantum non-Hermitian lattice models." Physical Review B  2022, 105, 064302.

[44] Lee, Ching Hua, et al. "Topolectrical circuits." Communications Physics 2018, 1, 39.

[45] Imhof, Stefan, et al. "Topolectrical-circuit realization of topological corner modes." Nature Physics  2018, 14, 925-929.

[46] Albert, Victor V., Leonid I. Glazman, and Liang Jiang. "Topological properties of linear circuit lattices." Physical review letters 2015, 114, 173902.

[47] Yang, Huanhuan, et al. "Observation of symmetry-protected zero modes in topolectrical circuits." Physical Review Research 2020, 2, 022028.

[48] Helbig, Tobias, et al. "Generalized bulk–boundary correspondence in non-Hermitian topolectrical circuits." Nature Physics  2020, 16, 747-750.

[49] Hofmann, Tobias, et al. "Chiral voltage propagation and calibration in a topolectrical Chern circuit." Physical review letters  2019, 122, 247702.

[50] Ezawa, Motohiko. "Electric circuits for non-Hermitian Chern insulators." Physical Review B 2019, 100, 081401.

[51] Ezawa, Motohiko. "Non-Hermitian boundary and interface states in nonreciprocal higher-order topological metals and electrical circuits." Physical Review B 2019, 99, 121411.

[52] Schindler, Joseph, et al. "PT-symmetric electronics." Journal of Physics A: Mathematical and Theoretical 2012, 45, 444029.

[53] N. Okuma and M. Sato, “Quantum anomaly, non-hermitian skin effects, and entanglement entropy in open systems,” Physical Review B, 2021, 103, 085428.

[54] T. Prosen, “Third quantization: a general method to solve master equations for quadratic open fermi systems,” New Journal of Physics, 2008, 10, 043026.

[55] S. Lieu, M. McGinley, and N. R. Cooper, “Tenfold way for quadratic lindbladians,” Physical Review Letters, 2020, 124, 040401.