Simulating Z 2 topological insulators via a one-dimensional cavity optomechanical cells array

We propose a novel scheme to simulate Z2 topological insulators via one-dimensional (1D) cavity optomechanical cells array. The direct mapping between 1D cavity optomechanical cells array and 2D quantum spin Hall (QSH) system can be achieved by using diagonalization and dimensional reduction methods. We show that the topological features of the present model can be captured using a 1D generalized Harper equation with an additional SU(2) guage structure. Interestingly, spin pumping of effective photon-phonon bosons can be naturally derived after scanning the additional periodic parameter, which means that we can realize the transition between different QSH edge states. © 2018 Optical Society of America OCIS codes: (270.0270) Quantum optics; (270.5585) Quantum information and processing; (120.4880) Optomechanics. References and links 1. M. Z. Hasan and C. L. Kane, “Colloquium: Topological insulators,” Rev. Mod. Phys. 82, 3045 (2010). 2. X. L. Qi and S. C. Zhang, “Topological insulators and superconductors,” Rev. Mod. Phys. 83, 1057 (2011). 3. P. Delplace, D. Ullmo, and G. Montambaux, “Zak phase and the existence of edge states in graphene,” Phys. Rev. B 84, 195452 (2011). 4. L. Ge, L. Wang, M. Xiao, W. Wen, C. T. Chan, and D. Han, “Topological edge modes in multilayer graphene systems,” Opt. Express 23, 21585-21595 (2015). 5. Y. J. Lin, R. L. Compton, K. Jiménez-García, J. V. Porto, and I. B. Spielman, “Synthetic magnetic fields for ultracold neutral atoms,” Nature (London) 462, 628-632 (2009). 6. M. S. Rudner, N. H. Lindner, E. Berg, and M. Levin, “Anomalous Edge States and the Bulk-Edge Correspondence for Periodically Driven Two-Dimensional Systems,” Phys. Rev. X 3, 031005 (2013). 7. K. Jimenéz-García, L. J. LeBlanc, R. A.Williams, M. C. Beeler, A. R. Perry, and I. B. Spielman, “Peierls Substitution in an Engineered Lattice Potential,” Phys. Rev. Lett. 108, 225303 (2012). 8. F. Mei, J. You, D. Zhang, X. C. Yang, R. Fazi, S. L. Zhu, and L. C. Kwek, “Topological insulator and particle pumping in a one-dimensional shaken optical lattice,” Phys. Rev. A 90, 063638 (2014). 9. A. He, W. Luo, Y. Wang, and C. Gong, “Wave functions for fractional Chern insulators in disk geometry,” New J. Phys 17, 125005 (2015). 10. H. Z. Shen, M. Qin, X. Q. Shao, and X. X. Yi, “General response formula and application to topological insulator in quantum open system,” Phys. Rev. E 92, 052122 (2015). 11. Z. C. Shi, H. Z. Shen, W. Wang, and X. X. Yi, “Response of two-band systems to a single-mode quantized field,” Phys. Rev. E 93, 032120 (2016). 12. H. Z. Shen, W. Wang, and X. X. Yi, “Hall conductance and topological invariant for open systems,” Sci. Rep. 4, 6455 (2014). 13. X. Li, E. Zhao, and W. V. Liu, “Topological states in a ladder-like optical lattice containing ultracold atoms in higher orbital bands,” Nat. Commun. 4, 1523 (2013). 14. S. Ganeshan, K. Sun, and S. Das Sarma, “Topological Zero-Energy Modes in Gapless Commensurate Aubry-AndréHarper Models,” Phys. Rev. Lett. 110, 180403 (2013). 15. J. Koch, A. A. Houck, K. L. Hur, and S. M. Girvin, “Time-reversal-symmetry breaking in circuit-QED-based photon lattices,” Phys. Rev. A 82, 043811 (2010). 16. W. L. Yang, Z. Q. Yin, Z. X. Chen, S. P. Kou, M. Feng, and C. H. Oh, “Quantum simulation of an artificial Abelian gauge field using nitrogen-vacancy-center ensembles coupled to superconducting resonators,” Phys. Rev. A 86, 012307 (2012). 17. B. Peropadre, D. Zueco, F. Wulschner, F. Deppe, A. Marx, “Tunable coupling engineering between superconducting resonators: From sidebands to effective gauge fields,” Phys. Rev. B 87, 134504 (2013). ar X iv :1 81 1. 06 16 7v 1 [ qu an tph ] 1 5 N ov 2 01 8 18. F. Mei, J. You, W. Nie, R. Fazio, S. Zhu, and L. C. Kwek, “Simulation and detection of photonic Chern insulators in a one-dimensional circuit-QED lattice,” Phys. Rev. A 92, 041805(R) (2015). 19. N. Goldman, I. Satija, P. Nikolic, A. Bermudez, M. A. Martin-Delgado, M. Lewenstein, and I. B. Spielman, “Realistic Time-Reversal Invariant Topological Insulators with Neutral Atoms,” Phys. Rev. Lett. 105, 255302 (2010). 20. F. Mei, S. L. Zhu, Z. M. Zhang, C. H. Oh, and N. Goldman, “Simulating Z2 topological insulators with cold atoms in a one-dimensional optical lattice,” Phys. Rev. A 85, 013638 (2012). 21. G. Y. Wang, D. Y. Wang, W. X. Cui, H. F. Wang, A. D. Zhu, and S. Zhang, “Direct conversion of a three-atomW state to a Greenberger-Horne-Zeilinger state in spatially separated cavities,” J. Phys. B: At. Mol. Phys. 49, 065501 (2016). 22. H. F. Wang, A. D. Zhu and S. Zhang, “One-step implementation of a multiqubit phase gate with one control qubit and multiple target qubits in coupled cavities,” Opt. Lett. 39, 1489-1492 (2014). 23. H. F. Wang, A. D. Zhu and S. Zhang, “Physical optimization of quantum error correction circuits with spatially separated quantum dot spins,” Opt. Express 21, 12484-12494 (2013). 24. H. F. Wang, A. D. Zhu, S. Zhang, and K. H. Yeon, “Optically controlled phase gate and teleportation of a controlled-NOT gate for spin qubits in quantum dot-microcavity coupled system,” Phys. Rev. A 87, 062337 (2013). 25. H. F. Wang, A. D. Zhu, S. Zhang, and K. H. Yeon “Deterministic CNOT gate and entanglement swapping for photonic qubits using a quantum-dot spin in a double-sided optical microcavity,” Phys. Lett. A. 377, 2870-2876 (2013). 26. H. F. Wang, A. D. Zhu, S. Zhang, X. X. Yi, and K. H. Yeon, “Local conversion of four Einstein-Podolsky-Rosen photon pairs into four-photon polarization-entangled decoherence-free states with non-photon-number-resolving detectors,” Opt. Express 19, 25433-25440 (2011). 27. W. X. Cui, S. Hu, H. F. Wang, A. D. Zhu, and S. Zhang, “Deterministic conversion of a four-photon GHZ state to a W state via homodyne measurement,” Opt. Express 24, 15319-15327 (2016). 28. H. F. Wang, A. D. Zhu, S. Zhang, and K. H. Yeon, “Quantum computation and entangled state generation via long-range off-resonant Raman coupling,” Quantum Inf. Process. 12, 2207 (2013). 29. H. F. Wang, A. D. Zhu, S. Zhang, and K. H. Yeon, “Simple implementation of discrete quantum Fourier transform via cavity quantum electrodynamics,” New J. Phys. 13, 013021 (2011). 30. T. J. Kippenberg, and K. J. Vahala, “Cavity optomechanics: back-action at the mesoscale,” Science 321, 1172-1176 (2008). 31. F. Marquardt, and S. M. Girvin, “Trend: Optomechanics,” Physics 2, 40 (2009). 32. I. Favero, and K. Karrai, “Optomechanics of deformable optical cavities,” Nat. Photonics 3, 201-205 (2009). 33. M. Aspelmeyer, S. Gröblacher, K. Hammerer, and N.Kiesel, “Quantum Optomechanics-throwing a glance,” J. Opt. Soc. Am. B 27, A189-A197 (2010). 34. I. Wilson-Rae, N. Nooshi, W. Zwerger, and T. J. Kippenberg, “Theory of Ground State Cooling of a Mechanical Oscillator Using Dynamical Backaction,” Phys. Rev. Lett. 99, 093901 (2007). 35. Y. C. Liu, Y. F. Xiao, X. S. Luan, Q. H. Gong, and C. W. Wong, “Coupled cavities for motional ground-state cooling and strong optomechanical coupling,” Phys. Rev. A 91, 033818 (2015). 36. X. Chen, Y. C. Liu P. Peng, Y. Zhi, and Y. F. Xiao, “Cooling of macroscopic mechanical resonators in hybrid atom-optomechanical systems,” Phys. Rev. A 92, 033841 (2015). 37. Y. Guo, K. Li, W. Nie, and Y. Li, “Electromagnetically-induced-transparency-like ground-state cooling in a double-cavity optomechanical system,” Phys. Rev. A 90, 053841 (2014). 38. F. Marquardt, J. P. Chen, A. A. Clerk, and S. M. Girvin, “Quantum Theory of Cavity-Assisted Sideband Cooling of Mechanical Motion,” Phys. Rev. Lett. 99, 093902 (2007). 39. Z. Q. Yin, T. Li, and M. Feng, “Three-dimensional cooling and detection of a nanosphere with a single cavity,” Phys. Rev. A 83, 013816 (2011). 40. J. B. Clark, F. Lecocq, R. W. Simmonds, J. Aumentado, and J. D. Teufel, “Sideband Cooling Beyond the Quantum Limit with Squeezed Light,” arXiv: 1605.08795 (2016). 41. M. Underwood, D. Mason, D. Lee, H. Xu, L. Jiang, A. B. Shkarin, K. Børkje, S. M. Girvin, and J. G. E. Harris, “Measurement of the motional sidebands of a nanogram-scale oscillator in the quantum regime,” Phys. Rev. A, 92, 061801 (2015). 42. J. M. Dobrindt, I. Wilson-Rae, and T. J. Kippenberg, “Parametric Normal-Mode Splitting in Cavity Optomechanics,” Phys. Rev. Lett. 101, 263602 (2008). 43. S. Grölacher, K. Hammerer, M. R. Vanner, and M. Aspelmeyer, “Observation of strong coupling between a micromechanical resonator and an optical cavity field,” Nature (London) 460, 724-727 (2009). 44. H. Xu, D. Mason, L. Jiang and J. G. E. Harris, “Topological energy transfer in an optomechanical system with exceptional points,” Nature (London), 537, 80-83 (2016). 45. L. Tian, “Cavity cooling of a mechanical resonator in the presence of a two-level-system defect,” Phys. Rev. B 84, 035417 (2011). 46. Y. D. Wang, and A. A. Clerk, “Using Interference for High Fidelity Quantum State Transfer in Optomechanics,” Phys. Rev. Lett. 108, 153603 (2012). 47. H. Xu, U. Kemiktarak, J. Fan, S. Ragole, J. Lawall and J. M. Taylor, “Observation of optomechanical buckling transitions,” Nat. Commun., 8, 14481 (2017). 48. Z. Q. Yin, W. L. Yang, L. Sun, and L. M. Duan, “Quantum network of superconducting qubits through an optomechanical interface,” Phys. Rev. A 91, 012333 (2015). 49. A. A.Clerk, F. Marquardt, and K. Jacobs, “Back-action evasion and squeezing of a mechanical resonator using a cavity detector,” New J. Phys. 10, 095010 (2008). 50. Q. Wu, J. Q. Zhang, J. H. Wu, M. Feng, and Z. M. Zhang, “Tunable multi-channel inverse optomechanically induced transparency and its applications,” Opt. Express 23, 18534-18547 (2015). 51. W. L. Li, Y. F. Jiang, C. Li, and H. S. Song, “Parity-time-symmetry enhanced optomechanically-inducedtransparency,” Sci. Rep. 6, 31095 (2016). 52. C. H. Bai, D. Y. Wang, H. F. Wang, A. D. Zhu, and S. Zhang, “Robust entanglement between a movable mirror and atomic ensemble and entanglement transfer in coupled optomechanical system,” Sci. Rep. 6, 33404 (2016). 53. D. Vitali, S. Gigan, A. Ferreira, H. R. 


Introduction
Topological insulators have become a very active area of research and raised widespread attention since the discovery of quantum Hall (QH) effects [1,2] in 1980.Normally, the topological insulators exhibit as insulator in the bulk and possess gapless edge states in surface or boundary.These surface states or boundary states are very robust against the disorder and perturbation, as long as the energy gaps of the system remain open, the topologically protected edge states cannot vanish.Besides, the topological classification and the essentially topological features of the system can be captured by a topological invariant.These novel properties stimulate that many theoretical models and experimental schemes have been proposed to investigate topological insulators, including graphene ribbons [3,4], cold atoms trapped in optical lattices [5][6][7][8][9], open system [10][11][12], p-orbit optical ladder systems [13], off-diagonal bichromatic optical lattices [14], and circuit-QED lattices [15][16][17][18].It has been verified that these systems can be used to simulate and investigate topological insulators exhibiting QH edge states via the synthetic gauge field, dimensional reduction, etc. methods.The QSH effect, moreover, which is found in materials displaying strong spin-orbit coupling, has opened another new area to investigate Z 2 topological insulators, a new family of topological states.The intrinsic spin-orbit coupling, as the origin of the QSH effect, determines that the QSH system contains two spin- 1  2 fermions where the two-component spin fermions can be described as QH states at equal but opposite "magnetic fields".Interestingly, this kind of topological insulator state possesses two pairs of edge states with opposite spin components at the boundary of the system.Recently, a number of models and proposals have been put forward to mimic Z 2 topological insulators.In [19], a specific 2D tight-binding model has been introduced to simulate a Z 2 topological insulator by engineering a synthetic gauge field to subject cold atoms.Afterward, Mei et al. [20] proposed a scheme to simulate Z 2 topological insulators using dimension reduction method with cold atoms trapped in a 1D optical lattice.
In our knowledge, although both cavity optomechanical system and topological insulators have rapidly developed, the simulation of topological features for the topological insulators based on the optomechanical system is rarely investigated yet.Motivated by this, we propose a conceptually simple and theoretically feasible method to simulate Z 2 topological insulators with a 1D cavity optomechanical cells array.We demonstrate that the original cavity optomechanical array model can be decoupled to two independent bosonic chains, which can be used to simulate Z 2 topological insulators by introducing a periodic modulation appropriately.The topological features of our model can be described by a generalized 1D Harper equation, and we find that the system exhibits trivial and nontrivial topological phases when the system parameters are chosen to different values.Particularly, the analogous spin pumping of effective photon-phonon bosons can be spontaneously derived with the periodic parameter varying at the range of (0, 2π).
The paper is organized as follows.In Section II, we derive the effective Hamiltonian of 1D cavity optomechanical arrays and realize the mapping between our model and 2D QSH system.In Section III, two examples are presented to simulate trivial and nontrivial topological insulators as the system parameters take different forms.We also discuss the Z 2 topological index which is used to distinguish different QSH phases.And we discuss the experimental feasibility of our scheme before conclusion.A conclusion is given in Section IV.

Model and mapping
The system of 1D cavity optomechanical array we study here as shown in Fig. 1.In this array, each unit cell is composed of a mechanical resonator with angular frequency ω n b and a cavity with angular frequency ω n a .The mechanical resonator couples to the cavity which is driven by an external optical field with angular frequency ω n p and rabi frequency Ω n via radiation pressure, forming a standard optomechanical subsystem.The total Hamiltonian of the system can be Fig. 1.Schematic of the 1D cavity optomechanical cells array.The nth unit cell contains a cavity mode with frequency ω n a labelled by a n and a oscillator with frequency ω n b labelled by b n , each cavity mode is driven by a classical pumping field with frequency ω n p .The coupling strengths between the adjacent cavities can be modulated to an appropriate value J n .
expressed as where a † n (b † n ) is the creation operator of the optical cavity (mechanical resonator) for the nth unit cell while a n (b n ) is its corresponding annihilation operator, and g n is the single-phonon optomechanical coupling strength resulted from radiation pressure.The first two terms represent the free energy of the nth unit cell, the third term describes the interaction between the cavity and mechanical resonator in nth optomechanical unit cell, the forth term stands for the interaction between the cavity field and the external classical laser field, and the last term expresses the hopping between the adjacent optomechanical unit cells with the hopping strength J n .
In the rotating frame with respect to the driving frequency ω n p , the above Hamiltonian becomes where −∆ n = ω n a − ω n p is the detuning between the cavity field frequency and pumping laser pulse.Under the condition of strong laser driving, we rewrite the operators as the sum of mean values and small quantum fluctuations to linearize Eq. ( 2), which means Substituting the above formulas into Eq.( 2), and ignoring the first and third order terms, and further dropping the notation "δ" for all the fluctuation operators for the sake of simplicity, the standard linearization Hamiltonian can be obtained as where G n = g n α n is the effective optomechanical coupling strength.For Hamiltonian in Eq. ( 3), the 1D cavity optomechanical array can be easily decoupled to two irrelevant bosonic chains by proceeding the diagonalization using the usual Bogoliubov transformation [67].
To derive the type of the beam splitter Hamiltonian, we consider the case of red-detuned regime, namely −∆ n ≈ ω n b , in which the Hamiltonian in Eq. ( 3) becomes To obtain the diagonal form of the above Hamiltonian, we make the bose mode transformation as follows For the case of strongly off-resonant regime G n J n , substituting Eq. ( 5) into Eq.( 4) and together with the rotating wave approximation, the effective diagonalization Hamiltonian can be written as where ω n A = −∆ n − G n and ω n B = −∆ n + G n are the eigenenergy of the effective photon-phonon bosonic modes A n and B n , respectively.Equation (6) shows that the original full Hamiltonian of the 1D cavity optomechanical array can be equivalent to a Hamiltonian consisting of two decoupled bosonic chains, as shown in Fig 2 .The first two terms and last four terms of Eq. ( 6) represent the on-site eigenenergy and adjacent hopping interaction of the two bosonic chains respectively.Obviously, the form of Eq. ( 6) is consistent with the tight-binding Hamiltonian that is investigated in topological system.It is natural to think that our present 1D optomechanical model can be used to simulate a system that exhibits nontrivial topological phase.
It is worth highlighting that, different from a Fermi system with spin-1 2 electron, our model is a standard Bose system, which contains two kinds of bosonic photon and phonon simultaneously.However, the quantized phenomenon in topological insulators originates from its topological properties, and is independent on Fermi statistical distribution of the electron with semi-odd spin.Besides, from the points of energy spectrum and Bloch theorem, electron, photon, and phonon present certain similarity.So, in the field of bosons with integral spin, the quantized topological states of light, sound [68], mechanical motion and so on should also be implemented similarly.If we regard the two decoupled bosonic chains as a two-component atomic gas mentioned in [20], our system with two kinds of effective photon-phonon bosons can be naturally used to simulate the Z 2 topological insulator with two-compont fermions, in which one of two effective photon-phonon bosonic chains A n can be equivalent to spin-up atomic gas while another can be equivalent to spin-down atomic gas.
In order to simulate the 2D tight-binding model [19] via the 1D cavity optomechanical cells array, it is significantly necessary to introduce another periodic parameter to map the second dimension of 2D model in present 1D cavity optomechanical cells array.Thus choosing the system parameters ω n A , ω n B , and J n 2 as ω n A = λ cos(2π βn + φ), ω n B = λ cos(2π βn − φ), and J n 2 = t.Here, parameter λ is the strength of periodic modulation term introduced into the system, β Fig. 2. Two decoupled bosonic chains.The initial optomechanical array system can be equivalent to two effective irrelevant bosonic chains via diagonalization under the condition of red-detuned regime.
is a real parameter, φ is an additional phase which can be changed continuously from 0 to 2π, and t represents the hopping strength.Noticing that this can be achieved by modulating system detuning with respect to driving field ∆ n = −λ cos(2π βn) cos(φ) and effective coupling strength G n = λ sin(2π βn) sin(φ).It has been verified that the manipulation of optomechanical coupling parameter can be relized by acting periodicly modulated driving light or microwave field on the optomechanical system.We can achieve the periodic modulation of system by changing the frequency of driving field in the way of periodic control [61].Substituting the above periodic parameters into Eq.( 6), the Hamiltonian becomes Obviously, if we associate the operators A † n (A n ) with the creation (annihilation) operators of the spin-up component fermions (represented by ↑) and B † n (B n ) with the creation (annihilation) operators of the spin-down component (represented by ↓) fermions, regard the parameter β as the magnetic flux, and treat the additional phase φ as the quasimomentum k y , our 1D cavity optomechanical cells array can be directly mapped into the 2D setup as in [19] to realize Z 2 topological insulators.In this configuration, the system can be described by a generalized 1D Harper equation with an additional SU(2) gauge structure as where Based on this equation, the topological phases of our 1D cavity optomechanical array can be captured via the energy spectrum and the distribution of the edge states.

Topological phases
Considering that we have used the condition G n J n at the moment of decoupling initial Hamiltonian, corresponding to that λ sin(2π βn) sin(φ) 2t.Therefore, the effective strength of the periodic modulation must satisfy λ 2 √ 2t, which is valid as revealed in [69].In the following we will investigate the spectral properties of Eq. (8).

Trivial topological phase
Firstly, we set the parameters as β = 1 2 t, λ = 15t, t = 1, and φ ∈ (0, 2π).In this case, each spin component, which is simulated by an effective photon-phonon bosonic chain (A n or B n ), can be mapped onto a 2D Hofstader model with π flux in each plaquette in Eq. ( 8), which implies that the system is invariant under the time-reversal transformation for the magnetic flux term only well-defined modulo 2π.Therefore, the system exhibits no QH effect and the QH edge states don't emerge.To demonstrate this point, we plot the energy spectrum of the system under the open boundary condition, as shown in Fig. 3.
Figure 3(a) clearly shows that the energy spectrum exists no band gap but contains two Dirac points with a liner dispersion connecting the valence band and conduction band.These results are exactly consistent with the viewpoints mentioned above.To further demonstrate this point, the distributions of wave functions around the Dirac point (0.5π, 0) corresponding to different effective photon-phonon bosonic "spin" modes (A n and B n ) are also plotted varying with the number of sites in Figs.3(b) and 3(c).The results indicate that no matter at the Dirac point or at other place in Brillouin zone, the distributions of wave functions are extended, which shows that there exists no local QH edge modes.
Furthermore, we find that the distributions of the states are completely same in Figs.3(b) and 3(c), which means that the motor directions of the two "spin" components A n and B n are from left to right simultaneously (or from right to left) with φ varying from 0 to 2π.The reason is that the matrix in Eq. ( 9) becomes M(n, φ) = Diag[∓λ cos(φ), ∓λ cos(−φ)] when β takes 1  2 t.The periodical modulation term of "spin-up" component A n becomes identical with the "spin-down" component B n due to the parity of cosine function, which implies that two irrelevant chains degenerate one chain and the system is a topologically trivial Z 2 insulator.

Nontrivial topological phase
Since the case of β = 1 2 t cannot generate the topologically nontrivial phase, we consider β = 1 3 t and keep λ = 15t, t = 1.Interestingly, we find that the energy spectrum splits into three bulk subbands and appears two band gaps, as shown in Fig. 4(a).One can see that there exists two degenerate states with opposite "spins" that are localized at two ends of the system at each gap, as shown in Fig. 4(b) (corresponding to the first energy gap).In more detail, for the "spin-up" component A n , there are two states which localized at the sites n = 1 and n = N respectively corresponding to φ = 0.5π and φ = 1.5π, as depicted in Fig. 4(b).Similarly, for the case of "spin-down", we find that a state is localized at n = N when φ = 0.5π and a state is localized at n = 1 when φ = 1.5π, which is just contrary to the "spin-up" case.These results indicate that each end of the system exists two degenerate edge states with opposite spins.All these features indicate that our system exhibits QSH effect in analogy to the 2D system which possesses two pairs of helical edge states.
Besides, Fig. 4(b) also depicts that the "spin-up" localized edge state which is at site 1 will transit to site N (the "spin-down" emerges the transition from N to 1) by adiabatically changing the parameter φ from 0 to 2π.We stress that these states are still localized at each edge during the whole process except for the case of φ = π in which the edge states integrate into the bulk state.In the analogous the QSH effect, we find that the above interesting process is the so-called spin pumping.More specifically, the topologically protected edge states with the opposite "spin" components can be pumped through the cavity optomechanical arrays by scanning the pumping laser phase adiabatically.This process is an effective way to realize the transition between different topologically protected edge states, which is significant to quantum information processing.

Fig. 3 .
Fig. 3. Energy spectrum and density distribution.(a) Energy varys with φ.The insets show two Dirac points with a liner dispersion.(b) Distribution of states corresponding to spin up component represented by symbol "•".The up panel corresponding to φ = 0, the middle panel corresponding to φ = 0.5π, and the bottom corresponding to φ = π.(b) Distribution of states corresponding to spin down component represented by symbol "+".The up panel corresponding to φ = 0, the middle panel corresponding to φ = 0.5π, and the bottom corresponding to φ = π.The system parameters are chosen as β = 1 2 t, λ = 15t, t = 1, and the lattice size N = 80.

Fig. 4 .
Fig. 4. Energy spectrum and density distributions.(a) Energy varys with φ.The insets exhibit fine structures corresponding two energy gaps.(b) Density distributions of spin up component represented by symbol "•" and spin down component represented by symbol "+" vary with φ.The up panel corresponding to φ = 0.5π, the middle panel corresponding to φ = π, and the bottom corresponding to φ = 1.5π.The system parameters are chosen as β = 1 3 , λ = 15t, t = 1, and the lattice size N = 80.