Tunable three-way topological energy-splitter

Strategically combining four structured domains creates the first ever three-way topological energy-splitter; remarkably, this is only possible using a square, or rectangular, lattice, and not the graphene-like structures more commonly used in valleytronics. To achieve this effect, the two mirror symmetries, present within all fully-symmetric square structures, are broken; this leads to two nondistinct interfaces upon which valley-Hall states reside. These interfaces are related to each other via the time-reversal operator and it is this subtlety that allows us to ignite the third outgoing lead. The geometrical construction of our structured medium allows for the three-way splitter to be adiabatically converted into a wave steerer around sharp bends. Due to the tunability of the energies directionality by geometry, our results have far-reaching implications for applications such as beam-splitters, switches and filters across wave physics.


Results
formulation. The group theoretic and topological concepts foundational to our approach hold irrespective of any specific two-dimensional scalar wave system. We choose to illustrate them here using a structured thin elastic Kirchhoff-Love (K-L) plate 41 for which many results for point scatterers are explicitly available 42 ; the geometrical ideas themselves carry across to photonics, phononics and plasmonics. Displacement Bloch eigenstates ψ | 〉 for Bloch-wavevector κ, n denoting the eigenmodes and ω κ n, the non-dimensionalised frequency; reaction forces at the point constraints, F x ( ), introduce dependence upon the direct lattice. The most straightforward constraints, sufficient for our purposes, are point mass-loading F x ( ) with the reaction forces proportional to the displacement via an effective impedance coefficient. Here l labels each elementary cell that repeats periodically to create the infinite physical plate crystal, and each cell contains = … p P 1 constraints. In an infinite medium the displacements are Bloch eigenfunctions n n n where | 〉 κ u n is a periodic eigenstate. The displacements satisfy the following completeness and orthogonality relations: up to first-order in κ |∆ |. This expansion will be used in the subsequent section, alongside symmetry considerations, to engineer the Dirac cones. A categorisation of planar structures that yield non-symmetry induced Dirac cones was shown in 43 . C 2v cellular structure. In this section we examine the cellular structure shown in Fig. 2(a). Spatially, this structure solely has σ v reflectional symmetry; however, in Fourier space, it has C 2v symmetry, due to the presence of time-reversal symmetry. In subsection 3.1, we utilise the expansion (5) and group theoretical considerations to demonstrate how an accidental Dirac cone is engineered. The effects of σ v symmetry breaking, on the bulk bandstructure, are discussed in subsection 3.2. Subsection 3.3 demonstrates how the strategic stacking of geometrically distinct media results in valley-Hall edge states 10,15 . A ZLM connected to a valley-Hall edge state is shown in subsection 3.4 alongside a justification for why this particular C 2v model does not allow for propagation around corners. Since the valley-Hall state is a weak topological state protected solely by symmetry, care must be taken to prohibit backscattering hence knowledge of the long-scale envelope is especially useful for finite length interfaces as it can used to minimise the backscattering. An asymptotic method, more commonly known as high-frequency www.nature.com/scientificreports www.nature.com/scientificreports/ homogenisation (HFH) allows for the characterisation of this long-scale envelope (see Methods); this is applied to a C 2v ZLM in subsection 3.4.
Engineering an accidental Dirac cone. Band coupling at high-symmetry point for C 2v structure. The point group symmetry of the structure, shown in Fig. 2(a) this is also the point group symmetry at N, Table 1). The C 2v point group arises from a combination of spatial (reflectional) and time-reversal symmetries; the latter relates κ κ → − . The group theoretical arguments used throughout this subsection, are reminiscent of those found in 21 although in our calculations we have applied an actual asymptotic scheme whereby we have judiciously chosen a small parameter with a distinguished limit.
The irreducible representations (IRs) at N are one-dimensional hence there is no symmetry induced degeneracy. Despite this, we shall demonstrate in this subsection how two of the IRs can be tuned such that an accidental degeneracy (that is not symmetry repelled) forms. The four solid bands in Fig. 2(c) (bands numbered 3-6 inclusive) are associated with the eigensolutions, shown in Fig. 3, these match the basis function symmetries of the the C 2v group (Table 1); hence this indicates that bands 3-6 are symmetry induced and the sequential ordering of them (lowest to highest) is deduced numerically, via the eigensolutions, as : B A B A  { , , , }   2  1  1 2 . It is expected, from the dispersion curves ( Fig. 2(c)) that the two bands that form the accidental degeneracy, namely A 1 , B 1 , have a strong influence on each other whilst the other two symmetry induced bands, B 2 , A 2 , will have a limited effect on the local curvature or slope of the A 1 , B 1 bands 44 ; the effect, by the bands that lie outside of bands 3-6, on the A 1 , B 1 bands, is expected to be negligible; to see these points mathematically we initially separate out the eigenket expansion Eq. (4) into two sets of bands; namely, the symmetry set eigensolutions (SSE), bands 3-6, and those that lie outside the SSE, Motivated by the orthogonality condition (3) and the expansion (6), we multiply Eq. (4), by (where ∈ l SSE, β ∉ SSE) before integrating over the primitive cell to obtain the following two equations,   where we have neglected terms which couple states outside the SSE to other states outside the SSE. If we let = ∈ n l SSE and κ κ κ = + ∆ 0 then the frequency term on the left-hand side is expanded to yield, Hence, from this expansion it is easy to see that the second summation in (11), that couples states within the SSE to those outside, falls into second-order hence the effective first-order equation is, l 0 and ∈ l SSE. Notably, the higher-order corrections, that encompass the coupling between bands within the SSE to those outside, provide the band curvature details away from a locally linear point. In this instance, Eq. (13) is a 4 × 4 matrix eigenvalue problem, where the Hamiltonian, with components H lj , is Hermitian. If, for a particular κ 0 , the first-order term is zero we would have to proceed to second-order; here additional terms would come from the fourth-ordered derivative, the ω κ l expansion and band coupling between outside SSE and inside SSE bands.
Compatibility relations and band tunability along NX. Bands tend to vary continuously except possibly at accidental degeneracies where modal inversion may occur which in turn leads to a discontinuity of the intersecting surfaces. Hence, the eigenfunctions continuously transform as you progress along a continuous IBZ path of simple eigenvalues. The associated IRs, that describe the transformation properties of the eigenfunctions, themselves smoothly transition into IRs that belong to the point groups along N Γ or NX.
In physical space the cellular structure only has σ v spatial symmetry, this is equivalent to σ h symmetry in Fourier space, Fig. 4(b). Recall the definition of a point group symmetry, i.e. any symmetry operator ∈ Γ R G that satisfies, κ κ = R G mod , where G is a reciprocal lattice basis vector; this implies that κ ∈ NX solely has the mirror symmetry operator, σ h within its point group. Similarly, for a κ ∈ ΓN , only the vertical mirror symmetry operator, σ ∈ C v v 2 , satisfies the point group criterion. The symmetries of the eigenfunctions, for a κ belonging to either of the paths, NX and NΓ, are shown within the basis functions column of the σ v,h Table 2. If we solely consider the two strongly coupled bands, represented by the IRs A 1 and B 1 , then the associated eigenfunctions transitional behaviour, away from N, is described by the σ v h , character table. Due to the continuity of the bands the A B , 1 1 IRs belonging to the C 2v table will transform into the IRs, of the σ v h , table, as we move away from N; the relationships between different IRs are more commonly referred to as compatibility relations 44,45 . Initially, we consider symmetry σ h , the eigenstates at N and along NX satisfy the following, Physically, this transition is also evident from the eigensolutions; as → B B 1 the eigensolution may also satisfy oddness relative to the x-axis (see Fig. 3). Similarly, at N and along NΓ, the eigenstates transform under σ v as, This implies that the bands A B ( , ) 1 1 (at N) are compatible with A A ( , ) (along NΓ). These compatibility relations are summarised pictorially in the unfurled IBZ path (Fig. 4(c)). Importantly, note that, in deriving Eq. (13) we have only assumed that κ 0 belongs to a particular symmetry set band (surfaces 3-6) (the band at κ 0 must be continuously connected to the same band at N). Therefore, the compatibility relations allow us to choose any expansion point along the the path ΓNX where the eigenfunction basis set, Eq. (4), transforms accordingly i.e.
In order to solve the 2-band eigenvalue problem, Eq. (13), we compute the determinant of the truncated matrix, (2019) 9:18939 | https://doi.org/10.1038/s41598-019-55485-0 www.nature.com/scientificreports www.nature.com/scientificreports/ where parity considerations 44,45 allows for simplification of the Hermitian matrix; the eigensolutions are evaluated at κ 0 . Solving the eigenvalue problem yields the following result, where the ± corresponds to the A B , bands, respectively. This result implies that the A B , bands have an identical slope, albeit with opposite gradients; hence, if, at N an instance can be found where ω ω > B A 1 1 then the bands will invariably cross along the path NX. The parametric variation afforded to us, and encompassed in the variable Λ A B , , merely increases or decreases the slope thereby increasing or decreasing the distance between N and the Dirac point. Note that the Dirac cone occurs along the spatial symmetry path, σ h , of the structure due to the opposite parities of the A B , bands; band repulsion occurs along the NΓ path 22 thereby resulting in a partial band gap along NΓ. If ω ω > B A 1 1 , then the partial gap along NΓ isolates the Dirac cone along a portion of the IBZ path, ΓNX. The distance between the Dirac and high-symmetry point is highly relevant for the transmission properties of the topological guide 26 stated that the transmission is better for short wave envelopes, as opposed to long wave  Fig. 2; σ v symmetry present within physical space cell. In this instance, B 1 curve lies above A 1 curve at N, hence there is band crossing along NX. www.nature.com/scientificreports www.nature.com/scientificreports/ envelopes, hence, for transmission post the nodal region, it is desirable to increase the distance between the Dirac cone and N. The latter is true due to the connection between the bulk and projected bandstructures 46 ; the bulk BZ is reduced to a one-dimensional BZ because the only relevant wavevector component for a straight guide is the one parallel to the ZLM. All wavevectors are projected onto the ΓM line in Fourier space, hence if the distance between N and the Dirac cone is increased then the Fourier separation between oppositely propagating modes, along the topological guide, would be increased. A mechanism to do this would be by altering the system parameters; Fig. 4(d,e) and Eq. (16) demonstrate that the slopes of the A and B bands can be increased or decreased by the system parameters thereby altering the position of the band intersection.
Breaking σ v symmetry. From the previous subsection we know that when the σ v symmetry is preserved an accidental Dirac degeneracy can be created; the bands coalescing along NX in Fig. 5(b) are parametrically engineered to do so. An important nuance is that the Dirac points are solely located along the two high-symmetry lines (HSLs), parallel to σ h , and not along the perpendicular HSLs (see Fig. 4(b)); this is critical when it comes to energy-splitting. The σ v symmetry is lost in Fourier space when the internal set of inclusions is rotated and this breaks open the Dirac point to create a band-gap, Fig. 5(c). The locally quadratic curves, in the vicinity of the former Dirac cones, are commonly referred to as "valleys" and they carry nonzero valley Chern numbers ( Fig. 5(a)) which in turn leads to the generation of valley-Hall edge modes 47,48 . The valley Chern number (C v ) is formally defined as follows,  15 . Both, C v and ϕ are dependent upon the perturbation strength and hence are not quantized quantities. Despite this, a topological integer, C sgn( ) v , exists that yields a bulk-boundary correspondence valid for specific interfaces 49 . These interfaces must ensure that regions of opposite C sgn( ) v are not projected onto the same point in Fourier space; hence for the C 2v case, from Fig. 5(a), it is evident that an interface associated with the ΓM direction will lead to well-defined topological integers. In the next subsection, we shall show how, the locations of nonzero C v , dictates how the geometrically distinct media are stacked.
C 2v adjoining ribbons. Attaching two topological media, with opposite C sgn( ) v yields broadband chiral edge states 10 . This is achieved by placing one gapped medium, above its σ v reflected twin; in essence, the stacking in Fourier space results in regions of opposite C sgn( ) v overlaying each other, this local disparity ensures the presence of valley-Hall edge modes. The two distinct orderings of the media create two distinct interfaces, as seen in Fig. 6(a) one of which supports only the even modes and the other only the odd modes. This evenness and oddness of the edge modes is inherited from the even and odd bulk modes, Fig. 3. The gapless curves are a symptom of the topologically nontrivial nature of the edge states; this is akin to the valley-Hall modes seen for the zigzag interface within hexagonal structures.
The simplicity of this construction, the apriori knowledge of how to tessellate the two media to produce these broadband edge states, and the added robustness 26 are the main benefits of these topological valley-Hall modes. The additional functionality of having a three-way topological splitter (Fig. 1) comes with a caveat: The Fourier separation between the valleys controls the intervalley scattering and the smaller separation in the square lattice, Figure 5. (a) IBZ (shaded region) within the Brillouin zone (BZ); circles indicate Dirac cone locations preperturbation, whilst ± denotes the signum of the post-perturbation C v . Dirac cones solely along single set of parallel HSLs, not both. Dispersion curves for the C v 2 case. Panel (b) pre-perturbation curves, when structure in Fig. 2(a) possesses σ v symmetry. Rotation of inclusion set in Fig. 2(a) removes σ v symmetry and yields the postperturbation curves shown in panel (c). Fig. 6(a), vis-a-vis that for graphene-like structures 26 leads to increased scattering. This can be mitigated as the Fourier separation can be artificially increased by parametrically increasing the distance between the Dirac cone and N in Fig. 5(b) thereby acting to increase the robustness of the edge states against shorter-range defects.
C 2v ZLM and absence of post-bend propagation. The property of the C 2v case that prohibits propagation around the bend is the absence of well-defined valleys with nonzero C v , along the vertical HSLs of the BZ, see Fig. 5(a). Hence, there is no arrangement that can be placed to the right of either stacking in Fig. 6(a) to obtain a ZLM perpendicular to the blue-orange interface, Fig. 6(b). The ZLM, Fig. 6(c), has a long-scale periodic envelope that can be captured using an effective medium theory 50 (see Methods). Knowledge of the long-scale envelope is especially useful for these finite length interfaces as it can used to minimise the backscattering as one has, in effect, a Fabry-Pérot resonator.
To summarise, for this C 2v case, there are ZLMs along straight interfaces, however the energy cannot navigate around a π/2 bend because there is no post-bend mode to couple with. C 4v cellular structure. We now extend the concepts illustrated in Sec. 3 to a cellular structure that possesses C 4v point group symmetry at Γ (see Fig. 7(a)). Due to this structure possessing two perpendicular mirror symmetries, as opposed to one, there exists regions of nonzero C v along both edges of the square BZ, subsection 4.1 (Fig. 7(b)); this will be shown to yield propagation around a corner (subsection 4.2) aswell as a three-way splitting of energy (subsection 4.3). An extensive comparison between the square structures, discussed within this article, and the earlier valleytronics models based upon graphene-like structures [16][17][18][19][20][26][27][28][29][30][31][32][33][34][35][36][37][38][39] will be pictorially shown at the end of subsection 4.3.
Breaking σ v,h symmetries. The C 4v case, Fig. 7(a), is reminiscent of the C 2v case but now with the addition of σ h symmetry in physical space. This reflectional symmetry yields additional Dirac cones along, a parallel set of HSLs, perpendicular to those connected with the σ v symmetry (Fig. 7(b)). This is evident, for the unperturbed C 4v case, in its dispersion curves, Fig. 7(c); note that we have plotted around the C 2v IBZ to clearly illustrate the correspondence between the two sets of dispersion curves, Figs. 4(b) and 7(c). The additional Dirac cone, for the C 4v case, along XM is due to the additional σ h reflectional symmetry in physical space. Rotating the inclusion set ( Fig. 7(a)) results in the breaking of both σ h,v symmetries thereby opening up a band-gap (Fig. 7(d)). Importantly, an identical band-gap is present whether we're plotting along the C v 2 or C v 4 IBZ's. www.nature.com/scientificreports www.nature.com/scientificreports/ In the subsequent section, we demonstrate how the additional reflectional symmetry enables mode coupling from the pre-bend to post-bend ZLM thereby allowing for energy navigation around a corner.
Propagation around a bend. A crucial property that allows for wave steering for the C v 4 case is the presence of Dirac cones along both edges of the BZ. Another important property is that, like the C v 2 case, both, even and odd edge modes exist, however they are now present along the same interface as opposed to different interfaces. The orthogonality of these opposite-parity modes ensures that they do not couple along the same edge. The presence of both parity modes along the same interface (for the C v 4 case) arises from the relationship between the orange over blue stacking and its reverse ( Fig. 8(a,b)). Specifically, it is clearly evident from Fig. 8(a,b) that a right propagating mode for one stacking is a left propagating mode on the other and vice versa. This special property is also what allows for the three-way splitting of energy (see subsection 4.3).
We now move onto deriving an edge mode for the C v 4 case. Due to there being only a single unique interface, we choose to use a Fourier-Hermite spectral method 51 , that purely finds the decaying solution along a single interface, as opposed to simultaneously along both; the latter occurs when the PWE method is used in Figure 7. Dispersion curves C v 4 case -(a) Cellular structure shown; maroon mass value of 1, blue mass value of 2, lattice constant of 2, centroid to vertex mass distance of 0.70. Pre-perturbation structure has σ v and σ h symmetries, both of these symmetries are broken in the post-perturbation structure. (b) Shows IBZ (shaded region) within the BZ; circles indicate Dirac cone locations pre-perturbation, whilst ± denotes C sgn( ) v , postperturbation. Unlike the C v 2 case (Fig. 4(b)), Dirac cones now present along both sets of parallel HSLs. (c) Preperturbation dispersion curves. We have opted to plot around the IBZ of the C v 2 case (Fig. 2(b)) in order to explicitly show the Dirac cone that arises from the added σ h symmetry. (d) Post-perturbation dispersion curves. If we were to plot along the C v 4 IBZ an identical band-gap, in location and width, would be present.

Figure 8.
Interface comparison between C v 4 case (a,b) and graphene-like structure (c,d) -representative hexagonal structure taken from 26 . Evidently, the two hexagonal zigzag interfaces that host ZLMs are distinct whilst, the two square interfaces, are identical under TRS. Even and odd-parity edge modes exist along the same interface for the C v 4 case and different interfaces for the graphene-like structures and C v 2 cases (Fig. 6(a)). This latter point is what allows for coupling between the pre-bend and post-bend modes, Fig. 10, for the C v 4 case but not the C v 2 . Crucially, this property is also what yields three-way splitting for the C v 4 case, Fig. 1 www.nature.com/scientificreports www.nature.com/scientificreports/ conjunction with two-dimensional periodic Bloch conditions. Hence, from Fig. 9, we clearly see that, for a variant of the C v 4 case, the orange over blue ( Fig. 9(a)) or blue over orange ( Fig. 9(b)) stacking yields an even-parity decaying mode. More specifically, the orange over blue stacking gives solutions to the right of Γ whilst the blue over orange yields solutions to the left of Γ; this implies that the two stackings host the same mode and are TRS pairs of each other. Note that a parametric variant of the C v 4 case was used to ensure faster convergence of the Fourier-Hermite spectral method. The local curvature, and thereby the characterisation of the envelope, is obtained for modes in the vicinity of Γ (see asymptotics in the bottom panels of Fig. 9). Similar to the earlier C v 2 structure the edge states that arise are topologically nontrivial and gapless 22 .
Transmission around a bend. The perturbed C v 4 system has valleys of nonzero C v along all HSLs of the BZ (Fig. 7(b)). This allows for the strategic arrangement of four structured media such that valleys of opposite C sgn( ) v overlay each other along, both, horizontal and vertical interfaces, Fig. 10(a). This strategic arrangement necessitates the existence of broadband ZLMs along both of these interfaces simultaneously; therefore, unlike the C v 2 case, energy is navigable around bends.
The four-cell arrangement shown in Fig. 10(a) encompasses the design of the nodal region (and by extending it outwards, the entire region) for the π/2 wave steerer and three-way energy-splitter. If the bottom-right inclusion set is rotated clockwise then a wave incident along the leftmost interface will follow the red arrows around the π/2 bend. The indistinguishable, pre-and post-bend interfaces, ensure that, as the energy traverses the turning point, an even-parity mode will couple into itself. An example of, topological wave steering around a bend, is shown in Fig. 10(b-d). Notably, the π/2 wave steerers observed within hexagonal structures require coupling between a zigzag mode with an armchair mode. The latter termination hosts topologically trivial edge states due to the overlaying regions of identical C sgn( ) v resulting in gapped states. Contrastingly, the structure shown in Fig. 10(a) allows for topologically nontrivial π/2 wave steering.
Similar to the C v 2 ZLM, the short-scale oscillations are discernible from the long-scale modulation. The importance of this long-scale modulation is numerically elucidated in the subsequent section.
Relevance of envelope to transmission around a bend. The characterisation of the energy-carrying envelope is important, as the tuning of it can lead to higher transmission along finite length interfaces. This principle is elucidated by examining the wave-steering example, Fig. 10(e,f); using finite element integration the intensity of the wave-field in each arm of such a steerer is calculated. The ratio of these intensities is the measure of the transmission of the wave steerer ( Fig. 10(f)); this quantity can be seen to oscillate rapidly across the band-gap. This is similar to the behaviour of conventional Fabry-Pérot resonators, where for maximal transmission an integer number of wavelengths must be completely contained in each lead. Thus the length of the interfaces is of importance for optimising the transmission. This effect is clearly seen by the contrast in transmission between Fig. 10(e)  and (b-d). Despite the paradigm utilising the valley-Hall topological phase, the robustness and bandwidth of the effect can be further increased, by parametic variation, introducing a TRS-breaking active component, nonlinearity and/or resonators within the nodal region.
Topological 3-way splitter. We now move onto the construction of the three-way energy-splitter; rotating the bottom-right inclusion set anti-clockwise, in Fig. 10(a); results in four partitions of geometrically distinct media. A wave incoming, from the leftmost interface, will now follow, both, the red and green arrows thereby splitting www.nature.com/scientificreports www.nature.com/scientificreports/ the energy three-ways. The resulting scattering solution, for a monopolar source, is shown in Fig. 1; the topological nature of the modes is demonstrated by the chiral fluxes. The three-way splitter can be tuned to a wave steerer by rotating the cellular structures in lower-right quadrant.
Importantly, the right-sided interface (blue over orange) is the reverse of the left-sided interface, hence a right propagating mode on the right-sided interface is identical to a left propagating mode on the left-sided interface. Mathematically, the latter implies that κ κ = − L R where κ L R , are the wavevectors along the left and right-sided interfaces respectively. Due to the chiral relationship between the interfaces the κ L and κ R modes conveniently have matching chirality. This allows for the seamless coupling between the left-sided and right-sided modes however there remains a phase difference to account for. In order for the left-sided mode to couple with the right-sided mode a phase must be acquired as the incident wave passes through the nodal point; this is similar to how an incident wave acquires a phase when it passes through a gratings coupler 52 . This phase difference between the modes implies that κ κ + = 0 L R thereby resulting in a zero imaginary component along the rightward lead, Fig. 11(b). Despite this phase difference, the chiral relationship between the two interfaces, Fig. 8(a,b), ensures that there is conservation of chirality (a necessary condition for topological mode coupling) throughout this four-region structured domain. The relationship between the interfaces is crucial in allowing the third lead to be triggered. This provides further evidence for why only two-way splitting has thus far been obtained for TRS-breaking topological systems [53][54][55][56] .

Methods
Dispersion curves and scattering solutions. The dispersion curves throughout the article were obtained using a combination of standard spectral methods as well as the Galerkin method. The standard scheme utilised an adaptation of the plane-wave expansion method to determine the eigenstates 50 . Specifically, a doubly periodic Fourier series expansion is employed by applying Floquet-Bloch conditions on opposite sides of the unit cell, resulting in a generalised eigenvalue problem, that is solved for the non-dimensionalised eigenstates presented throughout.
The dispersive behaviour of the edge states is obtained in a similar manner whereby we consider a stretched unit cell centred on the interface between the two media. Floquet-Bloch conditions are applied to both edges of the ribbon and the height of the ribbon is taken such that it is much greater than the decay length of the localised edge states. The relatively slower decay of the C v 4 edge mode (Fig. 9) required an alternative method; namely, a Galerkin method where the rate of the decay is built into the expansion of the wavefield through scaled orthonormal Hermite functions 51 .
The scattering solutions for point forcings are obtained through the solution of a system of linear equations by standard methods 42,68 . When a forcing is applied we utilise a Green's function approach where the total wavefield is given for P scatterers by   Figure 11. Real (a) and imaginary (b) components of the displacement field shown in Fig. 1. Notably in panel (a), the monopolar source triggers an even-parity ZLM along the left-sided interface; this mode couples into identical parity modes along all three outgoing leads. The absence of excitation along the right-hand lead, in panel (b), indicates that there is a phase difference between this lead and the other three excited leads. Panel (c) Two-way energy-splitting for hexagonal structure and three-way energy-splitting for C v 4 case -The suppression of intervalley scattering restricts graphene-like structures to two-way splitting of energy. The incoming ZLM, that has group velocity > v 0 g and wavevector +κ, is unable to couple to the post-nodal region ZLM, κ > − v 0, g , due to their differing valley indices.
www.nature.com/scientificreports www.nature.com/scientificreports/ Characterising energy-carrying envelope and relevance to robustness. The efficacy of transmitting energy around a bend, coupling modes between different leads within a network or even transmission through a straight ZLM is contingent upon the displacement of the mode at the turning, nodal or end point. Knowledge of the long-scale envelope is especially useful for these finite length interfaces as it can used to minimise the backscattering as one has, in effect, a Fabry-Pérot resonator. Examples of the characterisation of the energy-carrying envelope, using high-frequency homogenisation (HFH) 50 , for the C v 2 and C v 4 cases are shown. In addition to this, the interfacial dispersion curves for a variation of the C v 4 case are derived using a Fourier-Hermite spectral method 51 . www.nature.com/scientificreports www.nature.com/scientificreports/ To fully characterise the long-scale periodic behaviour of topological edge states along a crystal interface we utilise HFH, applying the methodology directly in reciprocal space 51 , to further bolster the plane wave expansion (PWE) method that was used to obtain the dispersion curves. This technique is a multiple scale asymptotic method, that (for non-degenerate curves with locally quadratic curvature) results in the following homogenised PDE, ij X X 0, 2 2 0 i j where f 0 is the long-scale envelope defined on the coordinate system X X ( , ) i j ; whilst the T ij coefficients fully encapsulates the short-scale behaviour (similar analysis can be carried over to any scalar and vectorial system 51,69-73 . The tensor coefficients T ij are geometrically dependant and, from the simple solution of the homogenised PDE, determine the envelope wavelength for a given frequency. These coefficients are determined entirely from integrated quantities of the wave-field in physical space. To avoid the need for regularisation (higher order corrections) we work in reciprocal space and calculate the T ij 's directly, using the PWE method. Our eigenvalue problem is recast into matrix form, where ω 0 and W 0 are the solutions obtained from the PWE method, and … + [ ] denotes the pseudoinverse. The explicit characterisation of the envelope is shown in Fig. 6(c).