Boundary-induced dynamics in 1D topological systems and memory effects of edge modes

Boundary-induced dynamics in 1D topological systems and memory effects of edge modes

Yan He College of Physical Science and Technology, Sichuan University, Chengdu, Sichuan 610064, China    Chih-Chun Chien School of Natural Sciences, University of California, Merced, CA 95343, USA.
July 20, 2019

Dynamics induced by a change of boundary conditions reveals rate-dependent signatures associated with topological properties in one-dimensional Kitaev chain and SSH model. While the perturbation from a change of the boundary propagates into the bulk, the density of topological edge modes in the case of transforming to open boundary condition reaches steady states. The steady-state density depends on the transformation rate of the boundary and serves as an illustration of quantum memory effects in topological systems. Moreover, while a link is physically broken as the boundary condition changes, some correlation functions can remain finite across the broken link and keep a record of the initial condition. By testing those phenomena in the non-topological regimes of the two models, none of the interesting signatures of memory effects can be observed. Our results thus contrast the importance of topological properties in boundary-induced dynamics.

I Introduction

The discovery of topological insulators and other materials exhibiting topological properties in their band structures has opened a growing research field (see Refs. Hasan and Kane, 2010; Qi and Zhang, 2011; Bernevig and Hughes, 2013; Shen, 2012; Asboth et al., 2016 for a review). Due to the bulk-edge correspondence, a nontrivial topological invariant in the bulk indicates localized edge modes that only appear in the presence of boundaries. Interesting transport properties can then arise from those edge modes. In addition to electronic systems, topological systems have been realized in ultracold atoms in optical potentials Miyake et al. (2013); Aidelsburger et al. (2013); Jotzu et al. (2014) and photonic and phononic systems Lu et al. (2014); Yang et al. (2015).

Recently there has been growing interest in studying nonequilibrium behavior of topological systems. For example, dynamics of topological superconductors in one and two dimensions undergoing a global quench of interactions reveals decaying topological quantities Sacramento (2014), and a Chern insulator experiencing a global interaction quench is expected to exhibit edge currents Caio et al. (2015). There have also been studies on other topological systems after global quenches of parameters Wang et al. (2015); Zeng et al. (2015); Sacramento (2016) or piecewise local quenches Grushin et al. (2015). Here we explore dynamics of topological models after boundary conditions are changed from periodic to open and vice versa. In the absence of topological properties, different boundary conditions should not cause observable effects in a large system. This is not the case for topological systems possessing edge modes that only emerge if there are boundaries. The emergence and disappearance of edge modes when boundary condition changes is the main theme of this work.

Two paradigmatic one-dimensional (1D) topological models will be implemented to demonstrate interesting dynamics induced by a change of boundary conditions. The Kitaev model Kitaev (2001) exhibits Majorana-fermion states and describes a topological superconductor. The Su-Schrieffer-Heeger (SSH) model Su et al. (1979) was originally proposed as a model for electronic transport in polyacetylene, and its topological properties are summarized in Refs. Shen, 2012; Asboth et al., 2016. The Zak phase of the SSH model has been measured using cold-atoms Atala et al. (2013). According to the classification of topological insulators Ryu et al. (2010), the Kitaev model belongs to the symmetry class D and the SSH model belongs to the class AIII. Therefore, the Kitaev model has particle-hole symmetry with -wave pairing while the SSH model has sublattice (or chiral) symmetry.

In both models we found the densities of edge modes reach steady-state values after the boundary condition transforms from periodic to open. Steady-state behavior in non-topological systems lacking interaction and dissipation has been reported Lai and Chien (2016), and here it is observable in the edge-mode dynamics of topological models. Another interesting feature of the edge mode dynamics studied here is the emergence of memory effects, where the rate of boundary transformation is recorded in the steady-state density of edge modes. Here we emphasize the quantum nature of the memory effects and will call them quantum memory effects. In different models the dependence of the steady-state density on the rate of boundary transformation can be different. Memory effects have been proposed in noninteracting quantum systems possessing interesting properties such as a tunable bound state Cornean et al. (2013), a geometric flat-band Lai and Chien (2016), or rate-dependent hysteresis Metcalf et al. (2016a). Here we show that memory effects from pure quantum dynamics should also be observable in edge-mode dynamics of topological models.

Moreover, certain correlations can remain finite across a link that is physically broken after the boundary condition changes if the systems are in the topological regimes. This is because of the initial intra-cell correlations that survives the change of boundary condition. In contrast, the correlations will be shown to simply decay to zero across a broken link if the system is topologically trivial. Therefore, topological models are capable of retaining correlations in dynamics. Since experimental realizations of the Kitaev and SSH models are possible, we will briefly summarize their experimental implications.

The paper is organized as follows. Sec. II briefly summarizes the Kitaev and SSH models and the formalism for investigating their dynamics, along with suitable initial conditions. Sec. III presents the dynamics of the Kitaev and SSH models after a change of boundary conditions. Evidence of quantum memory effects will be found in the edge modes, and topological effects will be contrasted by the results from the topologically trivial counterparts. Sec. IV concludes our study. Details of the two models and their analyses are summarized in the Appendix.

Ii Theoretical background

ii.1 Kitaev model

The Hamiltonian of the 1D Kitaev model is given by


Here we assume the system is arranged so that the -th site coincides with the first site. In our study we choose to be a uniform constant. The lattice constant is taken as the unit of length. We assume that


Here and are constants while and are functions of time allowing one to change the boundary conditions. In this work we take the linear form and for and and for to model a transformation from periodic boundary condition (PBC) to open boundary condition (OBC) with a characteristic time . Similar considerations also allow the system to transform from OBC to PBC. We note that Ref. Fagotti, has studied global reactions to local changes in spin chains by switching a few links, but the system remains periodic.

It is known Kitaev (2001); Bernevig and Hughes (2013) that the Kitaev chain is topologically nontrivial when . In this case, an open chain with an even number of sites has two zero-energy modes located at the two ends of the chain. On the other hand, when , the Kitaev chain is topologically trivial and there is no zero-energy edge modes in an even-numbered open chain. Details of the topological property of the Kitaev model are summarized in Appendix A.

ii.1.1 Time evolution and initial condition

The quantum dynamics can be obtained from the equation of motion in the Heisenberg picture


We define the following correlation functions to characterize the dynamics


Here is the initial quantum state, and in the following we will skip the label if there is no ambiguity. It is more convenient to transform to the Majorana fermion representations. Assuming that is real for any , we introduce the Majorana fermion operators and the Hamiltonian becomes


The evolution of the equal time correlations can be obtained accordingly. The density on site and current are


Here the current is from site to site . From the details shown in Appendix A, all those physical quantities can be expressed in terms of .

The equations of motion need to be supplemented with proper initial conditions. In Appendix A we present analytic and numerical methods for obtaining suitable initial quantum states with open or periodic boundary condition. The initial state is half-filled and in the ground state, and then the system evolves accordingly.

ii.1.2 Topological edge mode

Since the Kitaev model with is topologically nontrivial, an edge mode emerges as an eigenstate with OBC. The edge mode wavefunction satisfies the eigen-equation, Eq. (23) in the Appendix, with and . Thus,


Here are normalization factors and . The edge-mode creation and annihilation operators can be constructed by using Eq. (A) with . Note that in terms of Majorana fermions there are two edge modes located at the two ends of the chain. However, the two Majorana fermions will combine into one ordinary-fermion edge mode. The edge mode density corresponds to the occupation of the edge mode and can be obtain from Eq. (A) with .

ii.2 SSH model

The Hamiltonian of the SSH model Su et al. (1979) is given by


Here label the two sublattices. When , , and for all , the system is topologically non-trivial if because the Zak phase Atala et al. (2013), or equivalently the winding number Asboth et al. (2016), is finite. Moreover, the Zak phase changes signs from to . By including additional hopping or onsite terms, Ref. Li et al., 2014 shows rich topological properties of generalized SSH models.

We will consider the last inter-cell hopping coefficient to be time dependent with


Here we consider to be a linear transformation with a characteristic time . The dynamic equations can be derived from . The exact time evolution of the equal-time correlation function can be obtained and monitored in computer simulations. The initial state is half-filled and placed in the ground state with the corresponding boundary condition. Details of the initial condition is summarized in Appendix B.

ii.2.1 Topological edge modes

For the SSH model with an even number of sites and open boundary condition Shen (2012); Asboth et al. (2016), there can be two edge modes if or no edge mode if . Moreover, the edge modes survive in the presence of alternating onsite energies, but their energies are no longer pinned at zero. In the topological regime with alternating onsite energies and , the wavefunctions of the two edge modes expressed in terms of the annihilation operators are given by


They correspond to the eigen-energies and , respectively. Here are normalization factors.

Iii Result and discussion

iii.1 Boundary-induced dynamics of Kitaev model

For the 1D Kitaev model in the topological regime, we choose and , and use as the time unit. We present the results for sites, and the results are not sensitive to the number of sites if it is reasonably large and even. In Figs. 1 and 2 we show the results for the time evolution from closed boundary condition to open boundary condition. The transformation of the link between site- and site- is assumed to be linear in time with a tunable characteristic time .

Figure 1: Time evolution of the density ((a) for a 3D view and (b) for a top view) and the correlation Re (c) of the 1D Kitaev model as the system transforms from periodic to open boundary condition. The propagation speed is estimated as . Here , , , and .
Figure 2: The left panel shows the time evolution of the edge mode density for selected ramping times , respectively. Here , , and . The right panel shows the final edge mode density as a function of the ramping time . The inset of the right panel plots vs. and the red curve is a linear fit.

As shown in Fig. 1 (a) and (b), the perturbation due to the change of boundary condition propagates with a light-cone structure into the bulk. We found that the propagation speed is close to the maximal group velocity of the Kitaev model. The excitation energy follows the spectrum with . Following Refs. Ashcroft and Mermin, 1976; Metcalf et al., 2016b, the group velocity can be found from . For the parameters , , and lattice constant , the maximal group velocity is , which is very close to the slope of the light cone observed in Figure 1 (b). As increases, decreases. We have verified that the observed light-cone propagation speed of the density and correlation profiles also decreases with increasing , and the value agrees reasonably with the maximal group velocity.

After a transformation from periodic to open boundary condition, two end-points appear and as a consequence, an edge mode should arise at the boundaries. Figure 2 shows the rising of the edge mode density. There are two important features. First, the density of the edge mode reaches a steady state exhibiting a plateau after the transformation is completed. The steady state allows us to unambiguously identify memory effects in the growth of the edge mode, which is the second feature. As the ramping time gets longer, the edge mode has higher steady-state density. Such a dependence shows that the steady state of the edge mode is sensitive to the rate of the boundary change, and this is a manifestation of memory effects.

We emphasize that the isolated and noninteracting system considered here is not expected to equilibrate. The edge mode reaches a steady state because it is an eigenstate of the final Hamiltonian with open boundary condition, so its population remains after the boundary transformation is completed. In the right panel of Figure 2, we plot the steady-state value of the edge mode density as function of the ramping time . Within the range we explored, is exponentially approaching as increases. In the inset, we plot against , and it basically follows a straight line and confirms the exponential dependence.

Here we give pictorial explanations of the steady state and memory effect. A full analytical analysis is hindered by a lack of the full expressions of the eigenstates with open boundary condition. The steady state of the edge mode is because, after the boundary transformation is completed, the edge mode is an exact eigenstate of the final Hamiltonian with open boundary condition. The unitary evolution afterwards thus keeps the occupation of the edge mode intact, so a steady state emerges. We can consider two extreme cases: For a very short ramping time, one would expect, from the sudden approximation, that the the initial state is almost intact and the densities of the final eigenstates are found to be projections of the initial state. Since the edge mode was absent in the initial configuration, the projection leads to a small overlap and this implies a small steady-state density. On the other hand, for a very large ramping time , we are practically in the adiabatic limit where . Then the initially occupied eigenstate at half filling will smoothly evolve into the edge mode and the final density approaches . For arbitrary finite , the resulting will be somewhere in between those two extremes.

Memory effects usually arise if there are competing time scales. For example, a competition of driving period and relaxation time leads to rate-dependent hysteresis, which is another manifestation of memory effects (see Ref. Metcalf et al., 2016a and references therein). Here the intrinsic time scale is the hopping time determined by . The edge mode is special in the sense that it is inside the energy gap between the two bands and brings a different energy scale. Filling or emptying the edge mode is at a different time scale compared to the continuum, so different transformation times lead to different final edge-mode densities.

Figure 3: Time evolution of (a) density and (b) Re of the 1D Kitaev model from open to periodic boundary condition for . Here , , and .

In Figs. 3 and 4, we show the time evolution from open to periodic boundary conditions. The transformation of the boundary is again assumed to be linear in time with a characteristic time . In the initial state with open boundary condition, there is an edge mode and it should decay after the system becomes periodic. Interestingly, the edge mode density exhibits oscillatory behavior as shown in Fig 4. The lack of steady-state behavior inhibits a search for memory effects in this case because the edge mode is not an eigenstate of the final Hamiltonian with periodic boundary condition. However, we found that the variance of the total current in the whole system is relatively small and this indicates the system closely follows the averaged behavior.

Figure 4: The left and right panels show the time evolution of the edge mode density and the variance of the total current (in units of of the 1D Kitaev model from open to periodic boundary conditions for , respectively. Here , , and .

To contrast the effects from topological properties, we also present the evolution of the 1D Kitaev model in the topologically trivial regime by taking and . In Figure 5 we show the results for a transformation from periodic to open boundary conditions following a linear time dependence with a tunable characteristic time . One important difference between the topological and non-topological cases is the different behavior of the correlation function across the broken link. As shown in Fig. 5 (b), in the non-topological regime the correlation vanishes completely when the link is broken. In contrast, the correlation decays but remains finite in the topological regime even after the link is physically broken, and one can observe this in Fig. 1 (c).

One may understand the different behavior of the selected correlation by resorting to the dimer picture of Majorana fermions, where pairs of Majorana fermions form. Firstly, . When we focus on the nearest-neighbor correlation between the Majorana fermions, reveals the correlation between and . Next, one may use two chosen sets of parameters to contrast the topological and topologically trivial cases. By choosing with finite in the topological regime, the Hamiltonian shown in Eq. (6) only has Majorana pairs of the type , which includes the pair . Therefore, an initially periodic system has correlations between the Majorana fermions on the two side of the link that will be broken after a boundary condition change. In contrast, in the topologically trivial regime one may choose with a finite . Then, Eq. (6) only has terms like , so is not directly correlated. For general parameters in the two regime, our numerical simulations confirm that the correlation in the steady state is finite in the topological regime and vanishes in the topologically trivial regime, which is another example of quantum memory effects of initial correlations.

Figure 5: The time evolution of (a) density and (b) Re of the Kitaev model from periodic to open boundary conditions for . Note the correlation Re decays to zero as shown in (b) (the rising part shown at the rightmost site). The system is in the topologically trivial regime with , , and .

iii.2 Dynamics of SSH model

For the SSH model, is the hopping between the A site and B site within one unit cell, and is the hopping between the B site and A site from two neighboring unit cells. The onsite energies for A and B sites are and , respectively. Here for a chain with even number of sites we take and in the topological regime and use as the time unit. We present results for systems with site, and the conclusions are insensitive to the total site number. In contrast to the edge mode of the 1D Kitaev model with open boundary condition, the edge modes of the SSH model are highly oscillatory in space in the initial state.

When , the eigenstates with open boundary condition are close to those with periodic boundary condition. Dynamic signals from a boundary transformation are barely observable on the density profile or correlation functions. In contrast, we found clear dynamic signatures of topological properties with finite values of . Due to the alternating on the two sublattices, the density and current distributions for the A and B sites are showing oscillating behavior, so here we show the result on A sites to avoid over-crowded plots. The results for B sites are very similar. In Figure 6 (a), we show the results from periodic to open boundary conditions with a linear transformation of the link connecting site and site with a characteristic time . Similar to the Kitaev model shown in Fig. 1 (b), the propagation of the perturbation into the bulk also shows a light-cone structure on the evolution of density profile, and we found that the propagation speed is roughly the maximal group velocity of the system.

Figure 6: Time evolution of the density profile on the A-sublattice for transforming (a) from periodic to open boundary condition and (b) from open to periodic boundary condition. Here , , , and . The density profile on the B-sublattice has a similar structure.

Figure 6 (b) shows the evolution of the density profile from open to periodic boundary conditions. In this case the edge modes are initially present, but they decay away in the time evolution. By inspecting the density and current profiles, it is not easy to identify features distinguishing the two types of boundary changes. However, if we focus on the dynamics of the edge modes, visible differences can be observed.

From Eq. (12), the edge mode creation operator can be expressed as a superposition of the fermion operators on different sites. The edge mode density, which reflects the occupation of the edge mode, can be expressed as a superposition of the time-evolved correlation matrix . By analyzing the density of the edge mode, we found the edge-mode density exhibits steady-state behavior as shown in Fig. 7 for both types of boundary changes. The steady-state values of the edge-mode density allow us to identify memory effects from the boundary-induced dynamics although the oscillations around the steady-state values are smaller in the case transforming to open boundary condition. The steady-state values of the case transforming from periodic to boundary conditions shown in Fig. 7 (a) depend explicitly on , while the steady-state values of the other transformation shown in Fig. 7 (b) are insensitive to . The rate-dependent edge-mode density in the former case again serves as evidence that quantum memory effects can be found in topological systems, while there is no memory effect in the latter case.

When compared with the Kitaev model transforming from periodic to open boundary conditions, one can see that the steady-state value of the edge mode density increases monotonically with in the Kitaev model, but the dependence is non-monotonic in the SSH model as one can see in the inset of Fig. 7 (a). One may view the alternating onsite energies of the SSH model as an internal bias, which can further tune the dynamics. Such a feature is absent in the Kitaev mode because the onsite energy is uniform. We have checked that the behavior of the edge-mode density of the SSH model depends on , and as the steady-state values are too small to be discerned.

Figure 7: Time evolution of the left edge mode density when the system transforms from periodic to open boundary conditions (left panel) and from open to periodic boundary conditions (right panel) for , respectively. The inset shows the steady-state density (s. s. density) as a function of for the case from periodic to open boundary conditions. Here , , and .

In Figure 8 we show the time evolution of the correlation function Im when the system transforms from periodic to open boundary conditions with two different ratios of . An open chain with even number of sites and supports two edge modes and is in the topological regime while one with does not have any edge mode and is in the non-topological regime. For the former case, we choose while for the latter case, we choose in Fig. 8. The onsite energies are for both cases. The link between site and site is broken after the transformation, and the correlation Im for the non-topological case (the case without any edge mode) decays to zero as expected, but the correlation remains finite for the topological case. We have checked the correlation remains the same for much larger system sizes and confirmed the finite correlation is not a finite-size effect.

The finite correlation in the topological regime is again due to the initial inter-cell pair, and it is another manifestation of memory effects in topological systems. One may understand this by analyzing the limit with a finite in the topological regime. In this case, pairs are localized across each -link and we can focus on the pair that will be separated in the transformation. We may treat the left and right sites of the pair as the two components of a pseudo-spinor. Then the system is similar to the problem of a quantum particle in a double-well potential and the Hamiltonian is , where are the Pauli matrices. The solution to is formally written as , where is a two-component spinor and denotes time-ordering. Then is associated with the accumulated dynamic phase when is reduced to zero. In contrast, in the other limit with a finite in the non-topological regime, the pairs are across those -links. In this case cutting a link has no effects on and it should remain zero after the transformation. Our numerical results confirm the emergence of memory effect in the selected correlation function.

Figure 8: Time evolution of Im of the SSH model from periodic to open boundary conditions for the case with two edge modes (black) and the case with no edge mode (red). The parameters are () for the topological (non-topological) case with . The topological case retains a finite correlation while the correlation in the non-topological case decays to zero. Here and .

The observation of quantum memory effects in topological models provides another example of how noninteracting quantum systems can support interesting dynamic phenomena. Quantum memory effects, however, are not unique to topological systems. When compared to previous studies utilizing a tunable bound state Cornean et al. (2013), a flat-band of dispersionless states Lai and Chien (2016), a competition between the characteristic times of driving and dissipation Metcalf et al. (2016a), the memory effects in the topological models discussed here offer an additional realization for retaining information of dynamics.

iii.3 Experimental implications

The Kitaev model has inspired experimental searches of Majorana fermions in condensed matter systems. For example, superconductor-semiconductor nanowire devices Mourik et al. (2012) and superconductor-ferromagnetic chain systems Nadj-Perge et al. (2014) are promising systems for realizing the model. Since those solid-state devices can be engineered, we envision a ring-shape device with a tunable link may be fabricated for simulating the phenomena discussed here.

The SSH model, on the other hand, has been realized using cold-atoms and optical superlattices Atala et al. (2013). However, cold-atom systems are usually tubes of atomic clouds and may not be suitable for investigating properties with periodic boundary condition. By using a rapid painting potential, ring-shape optical lattices have been realized Henderson et al. (2009). Future modifications of the painting-potential technique may allow the dynamics discussed here to be simulated. As discussed in Ref. Chien et al., 2016, direct observations of edge modes in cold-atom systems can be challenging, but techniques like the quantum-gas microscope Bakr et al. (2009) may help resolve the density distribution.

Since correlation functions such as of the Kitaev model or of the SSH model are useful in identifying memory effects, measurements of those correlations are highly desired but they can be challenging as well. One possible way is to measure the time evolution of the density correlation . From Wick decomposition Fetter and Walecka (2003), pairwise correlations such as and will be included in , and the dominant density-density contribution is subtracted. The correlations without memory effects contribute to a background while the correlations with memory effects may exhibit observable rate dependence.

Therefore, boundary-induced dynamics in topological models may still be testable in future experiments with quantum technologies. We caution that dissipation from the environment will eventually relax the system into thermal equilibrium and wash out any dynamic signature in the long-time limit. Measurement of the steady state and memory effect should be performed immediately after the boundary transformation.

Iv Conclusion

We have shown, by tuning only one link in 1D topological models like the Kitaev chain or the SSH model to change boundary conditions, interesting dynamical phenomena emerge. The edge modes exhibit steady-state behavior and allow us to identify quantum memory effects in topological systems when the boundary condition transforms from periodic to open. Correlations across the link broken in the transformation can remain finite in the topological regimes and vanish in topologically trivial regimes. The boundary-induced dynamics thus complements global quench dynamics and reveals non-equilibrium behavior due to topological properties.

Acknowledgment: Y. H. thanks the support of NSFC via Grant No. 11404228. C. C. C. thanks the Hellman Family Foundation for partially supporting this work.

Appendix A Details of Kitaev model

This topological property of the Kitaev model can be characterized by a bulk topological invariant called the Majorana number, which is defined by Kitaev in his original paper Kitaev (2001). For a one-dimensional -band non-interacting Kitaev chain, the Hamiltonian can be expressed in terms of the Majorana fermions as . Here and are annihilation operators of Majorana fermions with denoting momentum and labeling the bands (Note that grouping two Majorana fermions leads to one ordinary fermion), and the Majorana fermions are their own antiparticles. Then the Majorana number is given by . Here is the Pfaffian of the matrix . It can be shown that the Majorana number is equivalent to the number of the edge modes Kitaev (2001).

The correlations of the original fermions can be expressed in terms of the Majorana-fermion correlations . One can verify that and is pure imaginary if . Thus, we define and satisfies . For ,




Periodic boundary condition — When the boundary condition is periodic, the energy spectrum and eigenfunctions can be found analytically in momentum space by introducing the lattice Fourier transform and . Here with integer in the range . For convenience, we take as an even number. Then the Hamiltonian with for all becomes


with . It can be diagonalized by the Bogoliubov transformation


Here we assume that is real and is complex. Then,


Here is the phase of and sgn is the sign of . We also define .

Then the ground state satisfying leads to the correlation functions ( denote the sites)


Open boundary condition — For the open boundary case, there is no translational invariance. This makes it difficult to simplify the problem by transforming to momentum space. Instead, we will closely follow the method in Ref. E. Lieb, 1961, which implements a Bogoliubov transformation in real space. The Hamiltonian can be written in the form


Here is a Hermitian matrix and is an anti-symmetric matrix. For the Kitaev model with uniform parameter, and . We introduce the quasi-particle annihilation and creation operators


and define and . By requiring that diagonalizes the Hamiltonian, we find that and satisfy


Here is the eigenvalues of the above matrices, not to be confused with the superfluid quasi-particle dispersion. Therefore, we can find from . can be chosen to be real and satisfy , and the same is true for . The Hamiltonian becomes .

The Majorana fermion operators can be expressed as and . Then the correlation functions are given by and . From the correlations of Majorana fermions, one can find the correlation of original fermions as


and similar expressions for and .

We can also express the quasi-particles creation and annihilation operators in terms of the Majorana fermion operators as


The density of quasi-particles is given by

Appendix B Details of SSH model

The initial conditions of the SSH model can also be analyzed in a similar fashion.

Periodic boundary condition — Using lattice Fourier transforms and , the Hamiltonian with and becomes


Here we define , , and . The Hamiltonian can be diagonalized by introducing


Here . Then the Hamiltonian becomes


We assume that the SSH chain is half filled, so the lower band is fully filled and the upper band is empty. The equal-time correlations are given by


Open boundary condition — For the SSH model with open boundary condition, the Hamiltonian in the matrix form is


Here odd numbers label the A-sites and even numbers label the B-sites. Then we can numerically find the -th eigenvector corresponding to eigenvalue , where labels the components. Since the Hamiltonian matrix is symmetric, we can choose all the eigenvector to be real vectors.

In terms of the quasi-particle operator , the Hamiltonian becomes


Then it is straightforward to evaluate the equal-time correlation functions as



  • Hasan and Kane (2010) M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010).
  • Qi and Zhang (2011) X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83, 1057 (2011).
  • Bernevig and Hughes (2013) B. A. Bernevig and T. L. Hughes, Topological insulators and topological superconductors (Princeton University Press, 2013).
  • Shen (2012) S. Q. Shen, Topological Insulators: Dirac Equation in Condensed Matters (Springer-Verlag, 2012).
  • Asboth et al. (2016) J. K. Asboth, L. Oroszlany,  and A. Palyi, A short course on topological insulators: Band-structure topology and edge states in one and two dimensions (Springer, 2016).
  • Miyake et al. (2013) H. Miyake, G. A. Siviloglou, C. J. Kennedy, W. C. Burton,  and W. Ketterle, Phys. Rev. Lett. 111, 185302 (2013).
  • Aidelsburger et al. (2013) M. Aidelsburger, M. Atala, M. Lohse, J. T. Barreiro, B. Paredes,  and I. Bloch, Phys. Rev. Lett. 111, 185301 (2013).
  • Jotzu et al. (2014) G. Jotzu, M. Messer, R. Desbuquois, M. Lebrat, T. Uehlinger, D. Greif,  and T. Esslinger, Nature 515, 237 (2014).
  • Lu et al. (2014) L. Lu, J. D. Joannopoulos,  and M. Soljacic, Nat. Photonics 8, 821 (2014).
  • Yang et al. (2015) Z. Yang, F. Gao, X. Shi, X. Lin, Z. Gao, Y. Chong,  and B. Zhang, Phys. Rev. Lett. 114, 114301 (2015).
  • Sacramento (2014) P. D. Sacramento, Phys. Rev. E 90, 032138 (2014).
  • Caio et al. (2015) M. D. Caio, N. R. Cooper,  and M. J. Bhaseen, Phys. Rev. Lett. 115, 236403 (2015).
  • Wang et al. (2015) P. Wang, M. Schmitt,  and S. Kehrein, “Universal nonanalytic behavior of the hall conductance in a chern insulator at the topologically driven nonequilibrium phase transition,”  (2015), arXiv: 1511.05255.
  • Zeng et al. (2015) Y. Zeng, A. Hamma,  and H. Fan, “Thermalization of topological entropy after a quantum quench,”  (2015), arXiv: 1509.08613.
  • Sacramento (2016) P. D. Sacramento, “Edge mode dynamics of quenched topological wires,”  (2016), arXiv: 1601.05476.
  • Grushin et al. (2015) A. G. Grushin, S. Roy,  and M. Haque, “Response of fermions in chern bands to spatially local quenches,”  (2015), arXiv: 1508.04778.
  • Kitaev (2001) A. Y. Kitaev, Phys. -Usp. 44, 131 (2001).
  • Su et al. (1979) W. P. Su, J. R. Schrieffer,  and A. J. Heeger, Phys. Rev. Lett. 42, 1698 (1979).
  • Atala et al. (2013) M. Atala, M. Aidelsburger, J. T. Barreiro, D. Abanin, T. Kitagawa, E. Demler,  and I. Bloch, Nat. Phys. 9, 795 (2013).
  • Ryu et al. (2010) S. Ryu, A. P. Schnyder, A. Furusaki,  and A. W. W. Ludwig, New J. Phys. 12, 065010 (2010).
  • Lai and Chien (2016) C. Y. Lai and C. C. Chien, Phys. Rev. Applied. 5, 034001 (2016).
  • Cornean et al. (2013) H. D. Cornean, A. Jensen,  and G. Nenciu, Ann. Henri Poincaré 15, 1919 (2013).
  • Metcalf et al. (2016a) M. Metcalf, C. Y. Lai,  and C. C. Chien, Phys. Rev. A 93, 053617 (2016a).
  • (24) M. Fagotti, “Control of global properties in a closed many-body quantum system by means of a local switch,” ArXiv: 1508.04401.
  • Li et al. (2014) L. Li, Z. Xu,  and S. Chen, Phys. Rev. B 89, 085111 (2014).
  • Ashcroft and Mermin (1976) N. W. Ashcroft and N. D. Mermin, Solid state physics (Thomson Learning, Stamford, USA, 1976).
  • Metcalf et al. (2016b) M. Metcalf, G. W. Chern, M. Di Ventra,  and C. C. Chien, J. Phys. B: At. Mol. Opt. Phys. 49, 075301 (2016b).
  • Mourik et al. (2012) V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers,  and L. P. Kouwenhoven, Science 336, 1003 (2012).
  • Nadj-Perge et al. (2014) S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, A. H. MacDonald, B. A. Bernevig,  and A. Yazdani, Science 346, 602 (2014).
  • Henderson et al. (2009) K. Henderson, C. Ryu, C. MacCormick,  and M. G. Boshier, New J. Phys. 11, 043030 (2009).
  • Chien et al. (2016) C. C. Chien, S. Peotta,  and M. Di Ventra, Nat. Phys. 11, 998 (2016).
  • Bakr et al. (2009) W. S. Bakr, J. I. Gillen, A. Peng, S. Folling,  and M. Greiner, Nature 462, 74 (2009).
  • Fetter and Walecka (2003) A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems (Dover Publications, 2003).
  • E. Lieb (1961) D. M. E. Lieb, T. Schultz, Ann. Phys. 16, 407 (1961).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description