Inhibition of tunnelling and edge state control in polariton topological insulators

We address the inhibition of tunnelling in polariton condensates confined in a potential landscape created by a honeycomb array of microcavity pillars in the presence of spin-orbit coupling and Zeeman splitting in the external magnetic field. The coupling rate between the microcavity pillars can be strongly impacted even by weak out-of-phase temporal modulations of the depths of the corresponding potential wells. When such a modulation is implemented in truncated honeycomb arrays that realize a polariton topological insulator, which supports unidirectional edge states in the presence of spin-orbit coupling and Zeeman splitting, it allows controlling the velocity of the states. The origin of the phenomenon is the dynamical modulation with a proper frequency, which notably changes the dispersion of the system and the group velocity of edge states. We show that such a control is possible for modulation frequencies close to resonances for inhibition of tunnelling in a two-well configuration. Edge states considerably slow down, and even stop completely, when the modulation frequency approaches a resonant value, while above such frequency splitting of the edge states into wavepackets moving with different velocities occurs.

The uncommon dispersion of the edge states, whose unidirectional motion becomes apparent when wavepackets with localized envelope are used for edge state excitation, is determined by the physical factors leading to opening of the topological gap and, most importantly, by the internal structure (symmetry) of the insulator. Thus, photonic Floquet [20][21][22] and polariton 25,29 insulators are most frequently constructed using honeycomb arrays of waveguides or microcavity pillars, respectively. The symmetry of the underlying array determines the momentum intervals where edge states can appear for a particular truncation type. Elucidation of new mechanisms that allow controlling the dispersion of the topological edge states, hence the transport properties of the system, is an important problem from both fundamental and practical points of view.
Here we expose one such mechanism that keeps the internal structure of the topological insulator unchanged and that is based on shallow periodic temporal modulations of the depths of the potential wells in the honeycomb array. It relies on the phenomenon of inhibition of tunnelling that was introduced in quantum mechanics and that leads to a considerable suppression of tunnelling for certain modulation frequencies in time-modulated potentials 41,42 . Nowadays, the control of tunnelling by periodic forces is routinely used in non-topological systems 43,44 . In photonics, this phenomenon has been observed in waveguide arrays modulated in the direction of light propagation 45,46 (see also the recent review 47 ). Honeycomb structures frequently used for construction of topological insulators are especially suitable for observation of inhibition of tunnelling, since this effect requires out-of-phase modulation of the depths of all neighbouring potential wells 48 , a property that is easily realizable in honeycomb array.
We consider a polariton topological insulator built as a honeycomb array of microcavity pillars that can be readily fabricated experimentally and that has been already used for demonstration of nontopological polaritonic edge states 49,50 . Existence of topological edge states in such a system relies on the combination of polarization-dependent tunnelling between neighbouring pillars, analogous to spin-orbit coupling effect 51,52 , and Zeeman splitting in a magnetic field through the excitonic component of the polariton condensate. We show that the group velocity of the edge states can be controlled by a weak periodic out-of-phase temporal modulation of the neighbouring microcavity pillars, which can be realized, e.g., using electro-optic 53 , acoustic 54 , free-carrier based 55 , and other 56,57 techniques. When the modulation frequency approaches one of the resonant values, the group velocity of the edge states can be decreased practically to zero, thus leading to a decreasing displacement along the edge for states with a localized envelope. The phenomenon takes places with a suppressed emission into the bulk due to the topological protection, which ensures motion of the wavepacket only along the edge with a velocity controlled by the modulation.
We adopt a circular polarization basis, where the relation between the wavefunctions of the spinpositive and spin-negative components and wavefunctions of conventional TE (subscript y ) and with the width d and depth p , the separation between neighbouring potential wells in the resulting array  is given by a . Each potential well in  oscillates in time out-of-phase with its three nearest neighbours, which is a crucial ingredient for tunnelling inhibition. Here we assume that the array is truncated along the x -axis and that it is periodic along the y -axis, so that two zigzag edges appear for - profile, as depicted in Fig 1(a). Results that we obtained for the case of bearded edges are qualitative similar, thus they are not shown here. The y period of the structure is In Eq. (1) all spatial coordinates are normalized to the characteristic length L , all energy parameters (including potential depth and Zeeman splitting) are normalized to the characteristic energy is the effective polariton mass, while time is normalized to was used, which gives 0 0.7 meV   . Topological edge states appear due to simultaneous action of spin-orbit coupling and Zeeman splitting. These two effects, acting together, break the time-reversal invariance in the Eq. (1) that, according to work [3], should lead to opening of the topological gap between former Dirac points in the spectrum of the honeycomb array 25,29 . To illustrate the effect for the truncated array from Fig. 1(a), we first switch off the temporal modulation by setting 0 m = and search for eigenmodes of the structure in the form of Bloch waves ( , , ) + , which are periodic ( , ) ( , ) u x y u x y T   = + along the y -axis and localized ( , ) 0 u x y   ¥ = along the x -axis. Such Bloch modes are characterized by the dependence of their energy  on Bloch momentum k , which should be periodic with a period K 2 /T p = according to the Floquet theorem. Such dependence is found from the linear eigenvalue problem: and presented in Fig. 1(b), which shows a zoom of the dispersion relation between Dirac points. We studied a structure containing 44 microcavity pillars per unit cell (one y -period of the array): four such periods are shown in Fig. 1(a). The presence of the topological gap with unidirectional edge states in it (red and blue lines between remnants of Dirac points at K/3 k = and 2 K/3 k = ) is readily visible. Since our system is spinor and transforms at 0 b = into two independent equations, there are always two similar groups of bands shifted by 2W and characterized by the dominance of different components, y + or y -, in spinor wavefunction 29 . We deliberately selected a sufficiently large Zeeman splitting

W=
to have two upper bands from one of these groups well isolated from the rest of the spectrum [ Fig. 1(b)]. The wavefunctions emerging from these bands have a dominating ycomponent. The width of the topological gap generally increases with increase of the spin-orbit coupling strength b . From now on we set The topological mode associated to the red branch (superscript r ) resides on the right edge of the insulator, with a group velocity that is negative. Therefore, the mode moves along the negative direction of the y -axis. In contrast, the mode from the blue branch (superscript b ) resides on the left edge of the insulator and at 0.55 K k = moves in the positive direction of the y -axis. One observes the appearance of charge-2 vortices in each well in the y + component, a phenomenon that is a consequence of the spin-orbit coupling specific for polaritons.
As the group velocity of the edge states is defined by / k ¶ ¶  , their dispersion is defined by The motion of the edge state becomes apparent in the presence of a broad localized envelope along the y -axis. An illustrative example of the evolution of the edge state from the right edge [red dot in Fig. 1(b)] with the Gaussian envelope in the unmodulated array is shown in Fig. 2. According to Fig. 1(b) the velocity of such state is 0.1 v »-, therefore it traverses about 60 periods of the structure and notably broadens over the time interval shown. Notice that due to topological protection, the wavepacket moves along the edge and does not penetrate into the bulk. Our main goal here is showing that the velocity of the wavepacket along the edge can be efficiently controlled, while its broadening can be practically cancelled, by a suitable weak temporal modulation of the potential (microcavity pillars). To gain insight into the physics underlying the mechanism, we first address the impact on tunnelling of the out-of-phase modulation of potential wells in a simple two-pillar structure. Note that we take into account all physical effects that are included into model (1), i.e. Zeeman splitting and spin-orbit coupling. The latter dictates that the mode of a single potential well exhibits a bell-shaped profile in the ycomponent and carries a charge-2 vortex in the substantially less localized y + component, as visible in the insets of Fig.  3(a). When the pillars are brought in contact with each other one observes Josephson oscillations, namely periodic switching between the pillars. For selected parameters of the structure (that we optimized to get sufficiently large velocity v of the edge states), the period of oscillations at 0 m = is b 13.15 T » . If an out-of-phase temporal modulation of the potential wells is imposed, the evolution of the t -periodic system is governed by quasi-energies, rather than by usual energy  as in the static system. The solutions in this case can be found in the form is the time-periodic function whose period is defined by the period of the temporal modulation. As we found numerically, in modulated two-pillar system quasi-energies  of localized modes practically (but not exactly) coincide at a set of resonance frequencies w that depend on the system parameters, such as depths of potential wells, separation between them, and spin-orbit coupling strength. At these frequencies the modes of the system do not experience dephasing 41,45,46 , and therefore the input wavepacket that may excite an arbitrary combination of such modes remains practically unchanged upon evolution, i.e., tunnelling is practically arrested. For the parameters studied in this paper, the largest modulation frequency at which quasi-energies coincide in two-pillar system is 0.75 w » . To illustrate the above-described effect using a single-pillar excitation, we calculate the timeaveraged norm in the excited (right) potential well defined as: where 0 U is the norm in the excited well in the initial moment of time 0 t = , and averaging is performed over . The dependence av ( ) U w is shown in Fig. 3(a). Maxima in this dependence are indicators of inhibited tunnelling -as one can see, there are multiple resonances, each of which corresponds to a minimum of the difference between quasi-energies. Inhibition is not complete, because in the system under study the quasi-energies do not coincide exactly (otherwise one would have av 1 U = at the maxima), but coupling is delayed drastically in resonances. The phenomenon is visible in the evolution dynamics in 0 y = cross-sections for 0.65 w = [out of resonance, blue dot in Fig. 3(a)] and for 0.75 w = [primary resonance, red dot in Fig. 3(a)], in Figs. 3(b) and 3(c), respectively. At resonance [ Fig. 3(c)] the excitation remains in the excited well at all times and the fraction of norm switching to neighbouring well is negligible. The same picture is observed in secondary resonances, but quality of the tunnelling inhibition reduces with the decrease of the resonant frequency. In Fig. 3 we use optimal modulation depth 0.2 m = , because too shallow modulations cannot inhibit coupling, while too deep modulations result in strong radiation. We observed that resonance frequencies grow almost linearly with m . Notice that inhibition of tunnelling persists even if the presence of a phase shift between the modulations of depths of the two potential wells. However, its efficiency is maximal when the modulation is out-of-phase. This is the case for both, the two-pillar system and the extended topological insulator.  Since tunnelling between potential wells determines the evolution of the excitations also in complex extended systems, one expects that a similar resonant phenomenon may occur in the topological insulator in the presence of out-of-phase modulation of the potential wells. However, importantly, in that case the modes of the system are extended edge states having nonzero transverse momentum k . To study the impact of the potential modulation on their evolution we use the same input conditions as in Fig. 2 Fig. 4(b). There is a clear correlation between the frequencies at which the displacement is minimal (practically zero around 0.75 w = ) and the resonances for tunnelling inhibition in Fig. 3(a). For frequencies between the resonant ones the edge states traverse considerable distances ( 60  periods) along the edge. The displacement is practically zero in the principal (i.e. right) resonance and in the second (next to principal) resonance. It is also relevant to show the expansion rate of the edge states during evolution, and to quantify it we introduced the formfactor, defined as: 58 The quantity 1 ccharacterizes the width of the wavepacket. In Figs. 4(c) and 4(d), the dependence 1 ( ) c w is shown during evolution and in the final moment of time 1500 t = , respectively. There is a clear correlation between the width and the displacement. For example, the minimal width is achieved close to the resonance frequencies, where displacement is also minimal.  Figure 5 shows the output ydistributions at 1500 t = for modulation frequencies around the principal resonance. Remarkably, already for frequencies 0.59 w » the broadening of the wavepacket is notably reduced in comparison with the case of the unmodulated topological insulators. The displacement along the edge monotonically decreases (the interval of variation of wavepacket position is enormous and can be further increased for longer evolution times) when the modulation frequency approaches the primary resonance ( 0.59 0.68) w = -. This constitutes the central result of our paper, i.e., the possibility to control displacements in the topological system by a weak modulation of the structure. To gain further insight, we calculated the quasi-energy bands of the modulated topological insulator for different modulation frequencies w . Examples of quasi-energy bands are presented in the bottom row of Fig. 3. In the y -periodic topological insulator, quasi-energies  are functions of the Bloch momentum k . One can see that the topological gap with edge states in it persists even for the modulation amplitude 0.2 m = used here. However, the width of the gap notably decreases when the modulation frequency approaches the resonant values. Because the width of the gap decreases, the derivatives / k ¶ ¶  and 2 2 / k ¶ ¶  for topological branches determining the velocity and effective dispersion of the corresponding topological edge states decrease too, hence the reduced displacement and dispersion close to resonant frequencies. Figure 3(g) shows that the topological gap shrinks at 0.70 w » (it should be stressed that in periodic system the resonant frequency may depart from resonant frequency of two-pillar system). This is consistent with a strong expansion of the wavepacket into the bulk of the array observed in direct simulations shown in Fig. 5. Thus, in resonances the topological properties of the system degrade and as a result in a narrow range of frequencies one observes radiation into the bulk. When modulation frequency increases above the resonant value, the wavepacket splits into two wavepackets that may move in opposite directions, while radiation into the bulk cancels ( 0.78) w = . The amplitude of one of these wavepackets gradually decreases and for high frequencies one again observes considerable shifts in the negative direction of the y -axis ( 0.84) w = .
Similar dynamics was observed for other transverse momenta k and not only around principal, but also around secondary resonances in w . Figure 6 shows transformation of the output state with increase of modulation frequency, when it approaches the second resonance. Just as in the primary resonance, the displacement of the wavepacket monotonically decreases with w , until the wavepacket experiences an abrupt destruction around the resonance at 0.304 w = .
Summarizing, we studied the inhibition of tunnelling of topological edge states by introducing out-of-phase modulations in microcavity pillars constituting a polariton topological insulator. The modulation was found to drastically affect the displacement of the states, without destroying their internal structure. Also, we showed that efficient control of the displacement is possible for modulation frequencies close to frequencies at which inhibition of tunnelling occurs in two-pillar system. Our findings are relevant to optoelectronic devices based on topologically protected edge states. The concept may be also relevant in the case of nonlinear excitations and gap solitons 59,60 .