Hybrid Interference Induced Flat Band Localization in Bipartite Optomechanical Lattices

The flat band localization, as an important phenomenon in solid state physics, is fundamentally interesting in the exploration of exotic ground property of many-body system. Here we demonstrate the appearance of a flat band in a general bipartite optomechanical lattice, which could have one or two dimensional framework. Physically, it is induced by the hybrid interference between the photon and phonon modes in optomechanical lattice, which is quite different from the destructive interference resulted from the special geometry structure in the normal lattice (e.g., Lieb lattice). Moreover, this novel flat band is controllable and features a special local density of states (LDOS) pattern, which makes it is detectable in experiments. This work offers an alternative approach to control the flat band localization with optomechanical interaction, which may substantially advance the fields of cavity optomechanics and solid state physics.


I. INTRODUCTION
The electron localization in a crystal is an important phenomenon in solid state physics, which relates to many fundamental problems, e.g., the metal-insulator transition.Normally, the localization phenomenon is due to the presence of disorder, i.e., the celebrated Anderson localization [1].But in some special crystal lattices, electrons can be localized without any disorder and form a completely flat band in the whole Brillouin zone, the reason of which is the destructive wave interference resulted from the lattice geometry [2][3][4].This is just the flat band localization.Lieb, Kagome, Diamond, stub, and sawtooth lattices are some examples of the flat band lattices, and some general methods are proposed to design more lattice structures with flat bands [5].That the flat band electrons are of special interest is because that, due to the quenched kinetic energy, tiny interaction can induce some exotic correlated ground states, such as ferromagnetism [2,[6][7][8], superconductivity [9], and Wigner crystals [10][11][12].However, though the flat band electrons have been intensively studied in last three decades, it has not been experimentally confirmed in natural materials due to the complexity of real materials.
Most recently, a essential progress about the flat band lattice has been achieved in artificial quantum lattice systems, where flat band lattices have been realized in experiment, and flat band localization is observed as well [13][14][15].For example, it is reported that the Lieb lattice has been realized in various quantum systems, such as photonic crystals [14][15][16][17], cold atoms [18], artificial electron lattice on metal surface [19,20].Interestingly, the flat band localization are clearly demonstrated in both the photonic crystals [14,15,21] and the artificial electron lattice on metal surface [22][23][24], while the flat band superconductivity is realized in the cold atom system [18].
Here, we investigate the flat band localization in the bipartite optomechanical lattice.Our main finding is that, instead of the destructive interference resulted from the lattice geometry, the hybrid photon-phonon-interference in optomechanical lattice can also induce the flat band localization.An immediate consequence is that, in the optomechanical lattice, the structure of the flat band lattice, as well as the corresponding band structures, can be different from that in other lattices.This actually reflects the quantum interference characteristic of the optomechanical lattice.
Concretely, we propose the model of hybrid bipartite optomechanical lattice with one or two dimensional framework, whose unit cell consists of one optical sublattice and one optomechanical sublattice.A novel flat band is found, which is resulted from the hybrid interference between phonon and photon modes, and does not exist in other artificial quantum lattice systems.More interestingly, this new flat band corresponds to a special photon-phonon-localization pattern, i.e., photons are only localized on the optical sublattice, while phonons are only localized on the optomechanical sublattice.In the previous lattices which exist flat-band localization, excitations (photons or phonons) are only localized on one sublattice in comparison.This property can be used to identify this new flat band state in experiment.Furthermore, this new flat band together its special localization pattern is controllable by the driving laser applied into the bipartite optomechanical lattices, which is an obvious advantage compared with the previous flat band localization induced by the lattice geometry.Note that, for the case of two dimensional optomechanical lattice, we choose the optomechanical honeycomb and Lieb lattice as the examples.In the optomechanical Lieb lattice, we find three flat bands.We demonstrate that the three flat bands have two different origins.The middle one is the new flat band and it is because of the hybrid interference, and the others are due to the lattice geometry.The two kinds of flat band also have distinct photon-phononlocalization patterns.

II. MODEL AND METHODS
We consider a bipartite optomechanical lattice, whose unit cell consists of an optical sublattice (i.e., site A) and an optomechanical sublattice (i.e., site B).As shown in Fig. 1, this bipartite optomechanical lattice can have one or two dimensional framework.Applying a strong laser with frequency ω L on the bipartite optomechanical lattice, the system Hamiltonian in a frame rotating with ω L reads [48][49][50][51] where a (a † ) and b (b † ) are the annihilation (creation) operators of the optical and mechanical modes, respectively.The joint index (i, α) and (i, β) contain the subindexes α and β.The Hamiltonian can represent a one or two dimensional lattice by properly defining i, α and β.For the one dimension (1D) lattice shown in Fig. 1(a), α = A, B, and β = B, denoting A-sites and B-sites, respectively.And i indexes the unit cell.In words, for the case each site has a localized optical mode, which is evanescently coupled to the optical modes at adjacent sites.Its unit cell can be separated into two sublattice, the optical sublattice, site A and the optomechanical sublattice, site B.
In the optomechanical sublattice, a localized mechanical mode couples to the optical mode in the same site.
Extending to the case of two dimension (2D), two kinds of bipartite lattices are considered here, i.e., the hybrid honeycomb lattice and Lieb lattice, shown in Fig. 1(b,c).Now i denotes the unit cell, α = A, B, β = B correspond to the optomechanical honeycomb lattice.And α = A, B, C and β = B, C correspond to the optomechanical Lieb lattice, denoting A-sites, B-sites and Csites.Similar as the 1D case, the unit cell includes an optical sublattice and an optomechanical sublattice.Then the nearest neighbor optical hopping with strength J is considered and it is denoted by i, α, j, α ′ .The frequency detuning ∆ ≈ ω c − ω L with the cavity frequency ω c , and g is the linearized optomechanical interaction strength, which is much smaller than the mechanical frequency Ω.
In principal, the proposed hybrid bipartite optomechanical lattice is general and could be implemented in cavity (or circuit) QED system in the optical (or microwave) frequency range [52].As shown in Fig. 1(a), the 1D optomechanical lattice could be realized in the optomechanical crystals [36][37][38][39][40].The defects are generated by a appropriate local modification of the pattern of holes, which localizes the optical and mechanical modes on the crystals.The accessible system parameters for our model could be λ = 1, 550 nm, Ω/2π = 3.75 GHz, κ/2π ≈ 900 MHz, γ/2π ≈ 250 KHz, and g ≈ κ.Here λ, g are the optical wavelength and linearized optomechanical coupling strength under the condition of strongly optical driving, respectively, and κ, γ are the optical and mechanical decay rate.
Here we will examine the local density of states (LDOS) of lattice sites for both photon ρ O (ω; jα) and phonon ρ M (ω; jβ).In experiments, the LDOS of the photon at each site can be directly measured via a auxiliary probe laser.The photon and phonon LDOSes are formally defined as where G R OO (ω; jα, jα) and G R MM (ω; jβ, jβ) are the retarded Green's function of photons and phonons in real space, respectively.And the definitions are To calculate the Green's functions, we start from the Heisenberg-Langevin equation of motion in momentum space, where ψ is the vector of the photonic and phononic annihilation operators of lattices, and ξ is the vector of the noise operators of baths.Their specific formula depends on the lattice considered.For example, ψ = (a kA , a kB , b kB ) T and ξ = (a in kA , a in kB , b in kB ) T for the 1D case we consider.In k space, the retarded Green's function Then the retarded Green's function can be obtained via Fourier transformation: The diagonal components of G R (ω; k) matrix give the photonic and phononic retarded Green's functions.With Fourier transformation, the photonic and phononic LDOSes can be expressed as: where N is the number of unit cells of lattices.

A. hybrid one-dimension optomechanical lattice
Hybrid interference between the optical and mechanical modes exists in the proposed bipartite optomechanical lattice, which ultimately induces a new flat band together with the photon and phonon localization.
In the case of 1D array, shown in Fig. 1(a), transforming to the momentum space, the Hamiltonian becomes The eigenvalue E 0 (k) = 0 corresponds to the appearance of flat band, which is clearly exhibited in Fig. 2(a).Nor-mally, the lattice geometry in a pure photon or phonon 1D lattice will not induce the destructive interference, and hence no flat band appears in the normal 1D lattice.Moreover, Fig. 2 also shows a gap between the middle band and the up (or down) band, which is holden even when the middle flat band disappears under the condition ∆ = Ω.This gap is induced by the optomechanical interaction and its width is decided by the interaction strength.
It should be noticed that, formally, one may naively think this model is a stub lattice [53][54][55] due to the similar Hamiltonian.However, the hybrid 1D optomechanical lattice has three fundamental distinctions comparing with the stub lattice.First, the flat bands in the two systems have different physical origins.In the stub lattice, the flat band results from the lattice geometry induced destructive interference.But in the 1D bipartite optomechanical lattice discussed above, the flat band is induced by the hybrid interference between two different species of mode, i.e., the optical and the mechanical modes.Second, in other artificial quantum lattice systems, e.g., the photonic crystal, once the stub lattice is constructed, the energy dispersion is fixed.However, as we mentioned above, the energy dispersion of the optomechanical lattice can be tuned by adjusting the laser detuning ∆.So, as illustrated in Fig. 2, the flat band here can be changed into a dispersive band, and vice versa.Finally, featuring the optomechanical system, the bipartite optomechanical lattice is a one-dimensional lattice.It is contrary to the stub lattice, which is a quasi-one-dimensional lattice.
Different from the previous flat band in the lattice with special geometry structure (e.g., 2D Lieb lattice), this flat band is induced by the photon-phonon hybrid interference between transitions a A ↔ a B and b B ↔ a B , which does not exist in the pure photon (or phonon) 1D lattice.This ultimately leads to the result that it has distinguishable photon-phonon-localization property charactered by the special LDOS pattern (see Figs. 3 and 4).
Specifically, the photonic LDOS of sites A, B, and the phononic LDOS of site B are where , and κ (γ) is the dissipation of the optical (mechanical) mode.Here the subscripts O and M denote the photon and phonon, respectively.It is shown from Figs. 3 and 4 that, when the system energy is at the flat band [corresponding to ω = Ω in Figs.3(a,c,e)] under the condition ∆ = Ω, photons are only localized in the optical sublattice (i.e., A-sites), while phonons are only localized in the optomechanical sublattice (i.e., B-sites).This special LDOS pat- tern is detectable experimentally by probing the photon and phonon excitations in the lattices, and it offers a simple method to prove the emergence of this new flat band in our model.When the resonant condition ∆ = Ω is violated, the hybrid photon-phonon-interference is destroyed, leading that the flat band localization disappear [see Figs.In principle, the presented flat band localization is general and it also could be realized in a two-dimension bipartite optomechanical lattice.Here we choose the 2D honeycomb and Lieb lattices as the examples.Now the lattice periodicity leads to f (k) = 1 + e ik•a1 + e ik•a2 [with the basis vector a 1 = (1, 0), 2 )] for the honeycomb lattice, and f j (k) = 1 + e ik•aj , j = 1, 2 [with a 1 = (1, 0), a 2 = (0, 1)] for the Lieb lattice.(a,c).The parameters are same as that in Fig. 5 except for γ = 10 −3 J.
In Figs.5-8, we plot the energy structure and the LDOS pattern of the hybrid honeycomb and Lieb optomechanical lattices by numerically solving system Hamiltonian.Firstly, the flat bands are exhibited under the resonant condition ∆ = Ω, as shown in Figs.5(a,c) and the middle flat band in Figs.7(a,c).Note that the lattice geometry in a normal honeycomb lattice will not induce the destructive interference, and hence no flat band appears in the normal photon (or phonon) honeycomb lattice.Even for the normal Lieb lattice, there only is one flat band induced by the destructive interference resulted from its lattice geometry, and it will be holden when the geometry structure is not changed.
Secondly, similar as the case of 1D lattice, the flat band and the middle flat band, respectively, appearing in the optomechanical honeycomb and Lieb lattices are induced by the hybrid photon-phonon-interference.Because they This also can not be applied into the flat band localizations in the normal Lieb lattice induced by its special geometry structure, as shown in Figs.7(b,d) and Figs.8(b,d,f).
In addition, Fig. 9 plots the LDOS pattern in real space.It can been seen that excitations are dispersive, meaning that the localization does not exist for the case of ω = Ω ± J, as shown in Fig. 9(a,d).And Fig. 9(b,c) show two flat band localizations resulted from distinct origins, which corresponding to ω = Ω, and ω = Ω + g.Photons are localized at optical sublattice (i.e., A-sites), and phonons are localized at optomechanical sublattice (i.e., B,C-sites) at ω = Ω; i.e., the photon-phonon flat band localization in hybrid Lieb lattice.And Fig. 9(c) shows the intrinsic flat-band localization due to its geometry, in which photons and phonons are localized at optomechanical sublattices (i.e., B, C-sites).

III. CONCLUSION
We have investigated the flat band localization in the bipartite optomechanical lattice both in the cases of 1D and 2D, including a pure optical sublattice and an optomechanical sublattice.We shown that a new flat band together with a special photon-phonon-localization property is exhibited under the optimal photon-phononresonant condition i.e., ∆ = Ω.This leads to the results that the present flat band localization is detectable experimentally, and can be easily controlled by tuning the frequency of driving laser applied into the optomechanical sublattice.This study might inspire further explorations regarding the connection bewteen cavity optomechanics and the many-body physics.

FIG. 1 :
FIG. 1: (Color online) The sketch of hybrid bipartite optomechanical lattices with (a) one dimension framework, and two dimension (b) honeycomb, (c) Lieb frameworks.They both are bipartite lattices, which can be separated into two sublattices, denoted A and B, C. The sublattice A only has the optical mode (pure optical site), while sublattice B or C has coupled optical and mechanical modes (optomechanical site).The lower illustration in (a) denotes the implementation of one dimensional hybrid lattice with the optomechanical crystals.

FIG. 4 :
FIG. 4: (Color online) (a,b) correspond to the case of ω = Ω, ω = Ω − J, respectively, under ∆ = Ω.And (c) describes the LDOS pattern at ω = Ω, ∆ = Ω.The cycles denote the photonic LDOS pattern and the rectangles denote phononic LDOS pattern.(a) shows the special photon-phonon flat band localization, while the localization does not exist in (b) at ω = Ω − J under ∆ = Ω.(c) shows excitations are dispersive, meaning the disappearance of the special localization.Other parameters are same as that in Fig. 3 except for κ = 0.1J.
2(b), 3(b,d,f), and 4(c)].This demonstrates that the flat band localization in our model is controllable via adjusting the frequency of driving laser.Otherwise, it can be seen that the localization would not emerge when the system energy is not at ω = Ω, even if the condition ∆ = Ω is satisfied, as shown in Fig. 4(b).

B.
Flat band localization in two-dimension optomechanical lattice.

FIG. 6 :
FIG. 6: (Color online) The LDOS of 2D optomechanical honeycomb lattice for different decay rate κ.The local photon or phonon DOS of sites A and B under the condition (a,c,e) ∆ = Ω and (b,d,f) ∆ = Ω − J.Here the point ω = Ω corresponds to the flat band in Figs.5(a,c).The parameters are same as that in Fig.5except for γ = 10 −3 J.