A quantum graph approach to metamaterial design

Since the turn of the century, metamaterials have gained a large amount of attention due to their potential for possessing highly nontrivial and exotic properties—such as cloaking or perfect lensing. There has been a great push to create reliable mathematical models that accurately describe the required material composition. Here, we consider a quantum graph approach to metamaterial design. An infinite square periodic quantum graph, constructed from vertices and edges, acts as a paradigm for a 2D metamaterial. Wave transport occurs along the edges with vertices acting as scatterers modelling sub-wavelength resonant elements. These resonant elements are constructed with the help of finite quantum graphs attached to each vertex of the lattice with customisable properties controlled by a unitary scattering matrix. The metamaterial properties are understood and engineered by manipulating the band diagram of the periodic structure. The engineered properties are then demonstrated in terms of the reflection and transmission behaviour of Gaussian beam solutions at an interface between two different metamaterials. We extend this treatment to N layered metamaterials using the Transfer Matrix Method. We demonstrate both positive and negative refraction and beam steering. Our proposed quantum graph modelling technique is very flexible and can be easily adjusted making it an ideal design tool for creating metamaterials with exotic band diagram properties or testing promising multi-layer set ups and wave steering effects.


A quantum graph approach to metamaterial design
Tristan Lawrie 1 , Gregor Tanner 1* & Dimitrios Chronopoulos 2 Since the turn of the century, metamaterials have gained a large amount of attention due to their potential for possessing highly nontrivial and exotic properties-such as cloaking or perfect lensing.There has been a great push to create reliable mathematical models that accurately describe the required material composition.Here, we consider a quantum graph approach to metamaterial design.An infinite square periodic quantum graph, constructed from vertices and edges, acts as a paradigm for a 2D metamaterial.Wave transport occurs along the edges with vertices acting as scatterers modelling sub-wavelength resonant elements.These resonant elements are constructed with the help of finite quantum graphs attached to each vertex of the lattice with customisable properties controlled by a unitary scattering matrix.The metamaterial properties are understood and engineered by manipulating the band diagram of the periodic structure.The engineered properties are then demonstrated in terms of the reflection and transmission behaviour of Gaussian beam solutions at an interface between two different metamaterials.We extend this treatment to N layered metamaterials using the Transfer Matrix Method.We demonstrate both positive and negative refraction and beam steering.Our proposed quantum graph modelling technique is very flexible and can be easily adjusted making it an ideal design tool for creating metamaterials with exotic band diagram properties or testing promising multi-layer set ups and wave steering effects.
Metamaterials and quantum graphs.Since Veselago published his pioneering paper in 1968-The Electrodynamics of Substances with Simultaneously Negative Values of Permeability and Permittivity 1 -it has been understood that the manipulation of electromagnetic material properties can give rise to exotic wave effects.In this paper, Veselago showed that when a material has both permittivity and permeability less than zero, its refractive index would be negative.By balancing the wave vector and the Poynting vector, it was shown that the electromagnetic waves within would have anti-parallel phase and group velocities.As a result, Snell's law and the Doppler and Vavilov-Cherenkov effects are reversed.These results were considered purely theoretical, since no such material was known to exist.
In 2000, Pendry demonstrated a remarkable application of negative refractivity.Showing, given an ideal lossless slab of material with a refractive index of n = −1 in a medium with equal and opposite refractive index, a lens can be theoretically constructed that could focus light perfectly 2 .Crucially such a material could overcome the diffraction limit of a traditional lens.Pendry presents a practical way one could engineer such a material.Proposing the periodic arrangement of unit cells made from C-shaped metal elements or "split-ring-resonators" with wires, giving rise to an effective negative refractive index within some frequency domain.This was experimentally demonstrated by Smith 3 .Since then, these man-made materials or "metamaterials" have been investigated extensively.There have been countless proposals for different resonator designs and arrangements, giving rise to a large number of different wave effects.
Metamaterials function due to the interplay between the wavelength and the scale of the unit cell.For wavelengths of the order or less than that of the unit cell, the waves undergo Bragg scattering interacting directly with each resonant element.However, in the long wavelength regime, the material appears continuous with properties owing to the underlying structure 4,5 .By varying these structural or resonant elements, one can achieve the desired wave effects.The required constituents of a metamaterial continue to be a matter of debate and various modelling techniques are used such as, transmission line models 6 , boundary element models 7 and Finite Element Analysis 8 to gain insight into the design changes required to achieve the desired effects.These simulation techniques can be quite time consuming when setting up the models, let alone considering design modifications, so there is a need for simple and computationally cheap models for proof of principle studies.The proposed quantum graph approach can provide this simplified tool for designing metamaterial with a large range of different properties.
Since their introduction in 1997 9 , quantum graphs have become a most valuable model for studying quantum and wave effects in the context of quantum chaos and beyond.The theory of quantum graphs describe a network constructed from a set of vertices that act as point scatterers connected by a set of one-dimensional edges (bonds), wherein waves are free to travel.The wave transport through the network depends on the chosen properties of the point scatterers characterised by unitary scattering matrices.By changing the elements of the scattering matrix, one can customise the wave transport through the network.The versatility and ease of construction of quantum graphs, as well as the finite dimensionality and the exactness of 'semiclassical' expressions, make them ideal toy models to test ideas and research hypothesis'-see 10 and 11 for an overview.Besides quantum chaos, applications encompass modelling the vibrations of coupled plates 12 , formulating quantum random walks 13,14 as well as quantum search algorithms 15 .One advantage of a quantum graph formulation is, that eigenvalue conditions can be written in terms of a secular equation involving the determinant of a matrix of finite dimension.Similarly, the scattering matrix, describing the properties of open quantum graphs, can be given in terms of a closed form expression involving finite dimensional matrices.
Naturally the concise language of quantum graph theory lends itself well to describing metamaterials.The lattice nodes act as scatterers representing sub-wavelength resonant elements and the edges model phase modulation between each resonant element.For this work, the metamaterial is modeled as an infinite square periodic quantum graph embedded in a 2D real space.A generalisation to 3D is straightforward by adding additional edges and vertices in the extra dimension.With this simple quantum graph construction, a multitude of nontrivial wave effects can be modelled by adjusting the parameters of the graph metric and the node scattering matrix.By manipulating the resulting band diagrams, we will demonstrate some of these effects, namely positive and negative refraction as well as beam steering.
The paper is structured as follows: in section "Waves in finite open quantum graphs-modelling resonant element" we will describe the set up of open quantum graphs which will serve as the resonant elements in the metamaterials.We will then construct period quantum graphs in section "A quantum graph model for a metamaterial" and derive the equations for constructing plane wave solutions and the dispersion curves in section "Wave propagation in metamaterials-plane wave solutions".In section "Wave propagation in metamaterials-Gaussian beams", we consider wave energy transport based on a Gaussian beam representation of the wave solutions and focus then on the available parameter space used for band engineering in section "Band engineering".In the last two sections, section "Wave refraction between metamaterials" and section "Waves in N layered metamaterials", we focus on reflection/transmission at "interfaces" between different metamaterials both in the two and N-layer case.

Waves in finite open quantum graphs-modelling resonant element
Metamaterials are typically constructed from a periodic arrangement of sub-wavelength resonant elements.We consider here in particular two-dimensional metamaterials, a generalisation to three dimensions is straightforward.Using quantum graph theory, we can model a single resonant element as an open quantum graph, which is then periodically arranged to construct the metamaterial.
Consider a connected and open quantum graph Ŵ(V , E, L) as a model for a resonator.The graph is constructed from a set of bidirectional edges E = {e 0 , e 1 , . . ., e |E|−1 } with a corresponding metric L = {l e 0 , l e 1 , . . ., l e |E|−1 } , connected by a set of vertices V = {V 0 , V 1 , . . ., V |V |−1 } with corresponding boundary conditions.In addition, a set of four leads L = {L l , L r , L d , L u } (semi-infinite edges) are imposed in the left(l), right(r), down(d) and up(u) directions serving as incoming and outgoing channels in the two dimensional lattice, see the examples in Fig. 1.Both the edges in E and leads in L are endowed with a one-dimensional wave equation, (1) www.nature.com/scientificreports/Here ψ(z) = (ψ L (z), ψ E (z)) T is the vector of all lead and edge solutions and k is the wave number.The solu- tion on a given edge ψ e j (z) exists for z ∈ [0, l e j ] , while the solution on a given lead ψ L j (z) exists for z ∈ [0, ∞) .(For ease of notation, z represents the spacial coordinates on any edge or lead).The eigenfunction solution ψ(z) is expressed as a superposition of counter propagating plane waves, Here a ± = (a ± L , a ± E ) T represents the vector of all wave amplitudes heading in (−) or out ( + ) of a vertex between the leads and edges.The vertices in V are generally considered to be point scatterers, and the mapping from all incoming to all outgoing wave amplitudes at the vertices is described by a unitary matrix Ŝ , that is, We assume Neumann boundary conditions at individual vertices enforcing wave continuity and flux conservation at each vertex 10 although more general vertex scattering matrices can be considered 16 ).Enforcing Neumann boundary conditions, we obtain for the pqth element of the scattering matrix associated with vertex V j Here, δ pq is the Kronecker delta, and v V j represents the valency (number of edges and leads attached to a given vertex) at vertex V j .
To treat the connected graph as a model for a resonator, we wish to understand how waves incoming from any of the leads are redistributed and modulated into outgoing waves leaving through the leads.To do this a four dimensional graph scattering matrix ŜŴ is constructed, reducing the entire connected graph in Ŵ(V , E, L) to a single scatterer illustrated in Fig. 2.
To construct ŜŴ from the underlying connected quantum graph, one decomposes the scattering matrix into dis- crete events.First, one accounts for the prompt reflections from the leads back to the leads in a matrix ŜLL .Second, the waves on the leads are coupled to the edges in E with a matrix ŜEL .Third, the edge dynamics are described by an infinite series of scattered paths expressed as a Neumann series, ∞ n=0 ( P(k) ŜEE Here, P(k) maps the outgoing wave amplitudes on the edges a + E to the incoming wave amplitudes a − E , taking account of the phase modulation owing to the metric L , and ŜEE represents the vertex scattering between edges.Fourth and finally, the waves on the edges are coupled back onto the leads using a matrix ŜLE .Explicitly, we can write see Kottos and Smilansky 17 for details.The matrix terms in ŜŴ can be deduced from the open graph scattering matrix Ŝ in (3) using the block-matrix representation, In the next step, we will use this graph based resonator to construct a periodic metamaterial. (2) (3) (5)

A quantum graph model for a metamaterial
By periodically arranging these graph based resonant elements we will now model metamaterial behaviour.Consider placing each resonant element in a mesh, spaced by edges of length l.To keep track of the wave location in the graph, indices nm are introduced, where n representing steps in the horizontal direction and m in the vertical, as illustrated in Fig. 3.We use z again for the coordinate on the edges with z = 0 at a given resonant element and z = l at the neighbouring element for all edge (where we leave out the indices nm for convenience).
As before, all graph edges are endowed with the one dimensional wave equation (1).The vector of solutions is restricted to the edges (l, r, d, u) in each unit cell of the material at location nm, To solve for the wave solution across the metamaterial, we begin by evaluating the scattering at the central vertex n = m = 0 .Let a ± = (a ± l , a ± r , a ± d , a ± u ) T be the vector of wave amplitudes heading in (−) or out (+) of the central vertex, see Fig. 3.We have The waves once scattered, undergo phase modulation as they move along each edge allowing one to express the outgoing waves at the central vertex a + in terms of the incoming wave amplitudes c Due to the periodicity of the structure we can use Bloch's theorem to describe the wave solution in any unit cell nm in terms of the solution at the central vertex nm = 00 18 , that is, Here, κ x and κ y are the Bloch wave numbers in the horizontal and vertical direction, respectively.The incom- ing wave amplitudes at a given vertex and its neighbours are thus related by where By combining Eqs. ( 9), (10) and (12), we find the condition (8) In satisfying this condition, the incoming amplitudes at the central vertex a − can be found by solving for the eigenvector in Eq. ( 14).We then obtain the wave solutions in each unit cell nm throughout the metamaterial in the form where a + is obtained via (9).

Wave solutions for the infinite lattice
In the following, we will show how to construct plane wave and then beam-like solutions for the infinite lattice.The starting point for the construction are the solutions of Eq. ( 15) giving rise to dispersion curves of the form k(κ x , κ y ) or κ x (k, κ y ) .We encounter real or imaginary solutions of κ x (k, κ y ) corresponding to propagating or evanescent waves, respectively, the latter occurring in band gaps.The full wave solution can then be constructed from the corresponding eigenvectors in (14) extended throughout the lattice using the Bloch condition (11), as formulated in (16).
Wave propagation in metamaterials-plane wave solutions.We start by constructing plane wave solutions and will discuss a particular simple case here, where the dispersion curves can be given analytically.That is, we consider a lattice without resonators where each edge is directly connected through the vertex.Here, the vertex acts as a point scatterer, where the scattering matrix ŜŴ (k) in ( 6) is given by the Neumann scattering matrix (4) with v V = 4 , that is, A generalisation to arbitrary ŜŴ (k) is straightforward, examples will be shown in later sections.By solving Eq. ( 15), the dispersive properties of the lattice are given by the relation The resulting surface, or band, is shown in Fig. 4; it is periodic in κ-space over the Brillouin Zone (BZ), where BZ ∈ [− π l , π l ] 2 .The shape of the band is a function of the square periodic graph topology and the symmetry of the scattering matrix.To generate a solution consider taking a single value of k, resulting in a discrete ring of κ solutions on an iso-frequency contour.By picking a point on the contour, one constructs the full wave solution from the corresponding eigenfunction in Eq. ( 16) as plotted in Fig. 5.
The phase of the resulting eigenfunction is given by the chosen Bloch wave vector κ = (κ x , κ y ) T , while the energy flow is given by the Poynting vector J = (J x , J y ) T normal to the iso-frequency contour.To construct the components of J explicitly, one evaluates the flux on the graph edges.The 1D flux J of the wave function ψ nm on edge p is (15  www.nature.com/scientificreports/with ψ denoting the complex conjugate of ψ .The horizontal and vertical components of the Poynting vector can then be evaluated in terms of the waves on edges right(r)(or left(l)) and up(u)(or down(d)), respectively, that is A point on notation.In the example shown in Fig. 5, for every value of κ y , there are two corresponding values of κ x .Naturally, the choice of κ x informs the direction of energy flow.To delineate between waves traveling in opposite horizontal directions, the following notation is used: eigenfunction solutions with Poynting vector heading to the right are given index → and the corresponding eigenvector in Eq. ( 14) is expressed as a − as before where a ± := a ± (k, κ → x , κ y ) ; eigenfunction solutions with Poynting vector heading to the left are given index ← with the corresponding eigenvector in (14) relabeled as b − where b ± := a ± (k, κ ← x , κ y ) .Explicitly, In choosing the value of κ x ( κ → x or κ ← x ), one implicitly chooses the wave direction.For this reason the κ x dependence is dropped from the eigenfunction.
Wave propagation in metamaterials-Gaussian beams.Plane wave solutions travel in the direction of the wave vector κ , which does not necessarily agree with the direction of energy flow given by the Poynting vector J in (20).This can be observed when considering Gaussian beam solutions constructed via Fourier series with the eigenfunctions, shown above, as a basis.The solution of the Gaussian beam with focal point at n = n ′ , expressed in components nm (z) = (� nm,l (z), � nm,r (z), � nm,d (z), � nm,u (z)) T , is given as, ( 19) J(ψ nm,p (z)) := ℜ ψnm,p (z) where the integral is performed over the domain � = �(k) of the iso-frequency contour.Solutions outside are evanescent and so contribute nothing in the far field.The expansion coefficients α n ′ describing the beam profile at the focal point in terms of the eigenfunctions, given by the inverse transform, Here, ψ→ n ′ m,u is the complex conjugate of the right moving eigenfunction expressed for all m at horizontal location n = n ′ along the upward edges u and � n ′ m,u (z) is the beam profile for all m with focal point n = n ′ expressed on edge u.We chose the beam profile to be Gaussian, such that where σ is the width of the beam and the phase e iκ ′ y l determines the tilt angle of the beam with respect to the horizontal axis of the lattice.Since the space on the graph is discretised, so too are the spacial integrals.Figure 6 illustrates how varying the parameters σ and κ ′ y effect the shape and direction of the resulting Gaussian beam for a given iso-frequency contour �(k).
e iκ ′ y l , Figure 6.Iso-frequency contour for a metamaterial with resonant elements ŜŴ point as defined in (17) at k = 1/l together with the expansion coefficients α n ′ =0 in Eq. ( 23) and the resulting real components of the Gaussian beam profile nm with focal point set to n ′ = 0 .The beam profiles, Eq. ( 24), shown here are characterised by the parameters κ ′ y = 0 , 0 and 1/l and beam widths σ = 2.2l , 6.6l and 6.6l for (a), (b) and (c), respectively.

Band engineering
In a next step, we will use the quantum graph formulation for designing resonant elements leading to materials with "exotic" wave properties.As a simple example, we consider using the cross resonator illustrated in Fig. 1a as a resonant element for the metamaterial illustrated in Fig. 3.By setting the metric of the cross resonator to L = {l x /2, l x /2, l y /2, l y /2} as shown in Fig. 7a, one can modulate the phase in the horizontal and vertical direc- tion across the resonant element thus breaking the scattering symmetry between the four leads.By enforcing Neumann boundary conditions at the vertices as defined by Eq. ( 4), the scattering matrix between the four leads L has solution, found by Eq. (6).By solving Eq. ( 15) the dispersion relation is now given as By comparing Fig. 7b with 7c, we see that varying the phase modulation across each resonant element only in one dimension breaks the symmetry of the scattering and gives rise to a saddle shaped band, which can be exploited to yield negative refraction as shown in the next section.

Wave refraction between metamaterials
To exemplify the engineered refractive properties, consider two metamaterials represented by semi infinite square periodic quantum graphs, one with resonant elements given by ŜŴ 1 and the other with resonant elements given by ŜŴ 2 .The two materials are connected along the y direction and material 1 exists for the set of unit cells n 1 = (−∞, . . ., −2, −1} and material 2 exists for the set of unit cells n 2 = {0, 1, . . ., ∞) , as illustrated in Fig. 8.
The full wave solution nm = (� nm,l , � nm,r , � nm,d , � nm,u ) T across the two materials can be constructed by a linear superposition of counter propagating eigenfunction solutions ψ ⇋ j,nm in material j = 1 and 2, as expressed in Eq. ( 21), that is, where H(n j ) is the discretised Heavyside step function, that is, for j = 1, 2 and the coefficients A j and B j are associated with left and right moving waves, respectively.Thus, solutions with coefficients A 1 and B 2 represent waves incident on the interface, while solutions with coefficients B 1 and A 2 represent waves scattered from the interface.To determine the coefficients, we must satisfy the boundary conditions at the material interface.Naturally the wave solutions are given for a single value of k between the two materials enforcing k 1 = k 2 := k .As the system stays periodic in the y direction, the Bloch phase tangential to the interface also remain constant across the interface leading to the boundary condition κ 1,y = κ 2,y := κ y .This condition is illustrated in Fig. 9 for a chosen value of κ y = κ ′ y by a horizontal dashed grey line connecting the ) cos(κ x l)e ik(l x +l) e 2ik(l y +l) − 1 + cos(κ y l)e ik(l y +l) e 2ik(l x +l) − 1 − e 2ik(l x +l y +2l) + 1 = 0. Also plotted in the κ x , κ y plane are the iso-frequency contours of the second band for various values of k.Plot (b) is for l x = 1 , l y = 1 .Plot (c) is for l x = 0 , l y = 1.
iso-frequency contours at k = 1/l in material 1 and 2. The remaining unknowns are then the wave vectors normal to the interface κ ⇋ j,x (k, κ ′ y ) which can be obtained from the dispersion curve determined by solving Eq. ( 15) for each material, see the vertical dashed grey lines illustrated below.
To determine the scattering coefficients A j and B j , we rescale the eigenfunction solutions of each material such that the magnitude of the horizontal flux of each component is equal, that is, assuming κ ⇋ j,x ∈ R .Equally, the same condition can be enforced on the right edges, r.With the scaling choice (29), flux conservation across the interface reduces to and wave scattering at the interface can then be described in terms of a unitary scattering process.The corresponding interface scattering matrix Ŝ1,2 performing the mapping, (29)  can then be constructed by decomposing the interface scattering into two events.Event I describes a wave incident from material 1 onto material 2 with amplitude A 1 = 1 and B 2 = 0 , producing a reflected and transmitted wave with respective amplitudes r I and t I .Event II describes a wave incident from material 2 onto material 1 with amplitude B 2 = 1 and A 1 = 0 , producing a reflected and transmitted wave with respective amplitudes r II and t II , see Fig. 8.The interface scattering matrix takes on the form To evaluate the matrix elements, consider first event I.For simplicity we choose to evaluate the waves at location n = 0 , at coordinate z = 0 on edge l.Since the phase κ y is the same in both materials, it is sufficient to evaluate the solutions at m = 0 .Here, � (∀n≤−1)0,l (0)| n=0 and � (∀n≥0)0,l (0)| n=0 both exist on the same edge in the boundary region, so to stay consistent, there must be an equivalence between the incoming(−) and outgoing(+ ) wave amplitudes of these solutions, that is, This can be solved to give Exactly the same procedure can be done for event II where which yields the equivalence condition with solutions The results of this unitary interface scattering is shown in Fig. 10 for a Gaussian beam incident from material 1 with amplitude A 1 = 1 and no incident beam from material 2, that is, B 2 = 0.
The Gaussian Beam is constructed based on the method described in section "Wave propagation in metamaterials-Gaussian beams", where the basis is the full wave field nm (x, κ y ; k) as defined in Eq. ( 27), The expansion coefficients α −60 are as in Eq. ( 23) for a beam focal point n ′ = −60 and tilt given by κ ′ y = 1/l.As was shown in section "Band engineering", breaking the symmetry of the scattering matrix gives rise to a saddle shaped band resulting in the property of negative refraction at appropriate κ y values, see Fig. 7c.When choosing such a material on the right hand side, one obtains the reflection/transmission behaviour shown in Fig. 11.

Waves in N layered metamaterials
Consider now a system of N layered metamaterials as shown in Fig. 12, each with their own properties as defined previously.
Each material spans a domain defined by the set of unit cells n j .Here the wave function across all N materials is expressed as a linear superposition of counter propagating eigenfunction solutions of each material, that is, To evaluate the full wave function across all N materials, one must determine the coefficients A j and B j that satisfy the boundary conditions at all material interfaces.This is done by the Transfer Matrix Method 19 .In the previous section, the wave amplitudes between two materials were determined by constructing a scattering matrix mapping incoming to outgoing wave amplitudes between the metamaterials 1 and 2 as defined in Eq. (32).This matrix can be rearrange to give a transfer matrix that maps wave amplitudes from material 1 to material 2, Ŝ1,2 .
This procedure can be generalised to any arbitrary materials j and j + 1 giving Waves propagating across a given material j accumulate a Bloch phase which can be expressed in terms of the matrix where W j = ( max (n j ) − min (n j ))l is the width of material j.Having now formulated both scattering and propagation, one can express the wave amplitudes in material N in terms of the wave amplitudes in material 1.
We can now rearrange this transfer operator Ŝ1,N such that the entire system of layered materials acts as a single point scatterer by introducing the scattering matrix Ŝ1,N defined as Now by setting the incident wave amplitudes A 1 and B N , we determine the amplitudes of the scattered field by direct substitution into Eq.( 46).Knowing the wave amplitudes in material 1, A 1 and B 1 , it is trivial to determine all other amplitudes using Pj−1 Ŝj,j+1 .The results of this procedure is plotted in Fig. 13 for a three layered material.

Conclusion
We propose here a graph based technique for designing metamaterials allowing for a fast and flexible way to test resonant element proposals embedded periodically within a 2D square lattice.By modelling resonant elements in terms of open, graph-based scattering systems, we retain a connection with the underlying geometrical structure of the element which will inspire the engineering of physical metamaterials.This will make it possible to search for exotic band-diagrams, dispersion curves and wave effects linking these to an underlying array of resonant scattering systems.In this paper, we introduce the principles of graph based metamaterial construction including the set-up and we show how to obtain plane wave and beam-like solutions in the infinite medium.We then demonstrate the handling of boundary conditions at interfaces as well as the generalisation to N layered media.In each case, the computation can be reduced to low-dimensional matrix problems with the help of Bloch's theorem.Further results modelling various special wave effects will be presented in forthcoming publications.

Figure 1 .
Figure 1.Three examples of open quantum graphs are shown.(a) Represents a cross resonator, (b) a ring resonator and (c) some zig zag structured resonator.The edges between vertices are depicted in black and open leads in grey.

Figure 2 .
Figure 2. (a) Some arbitrary compact quantum graph Ŵ(V , E, L) . (b) The graph having been reduced to a single vertex.Noted are the lead wave amplitudes a ± L .

Figure 3 .
Figure 3.A 2D square periodic arrangement of quantum graph based resonant elements Ŵ(V , E, L) connected by edges of length l as a model for a 2D metamaterial.Location n = m = 0 shown in the subplot with labeled wave amplitudes, a ± and c − , traveling on the edges in the neighbourhood of the central vertex.

Figure 4 .
Figure 4.The band diagram for a lattice without resonant elements.Shown are iso-frequency contours for various values of k in the κ x , κ y plane.

Figure 5 .
Figure 5. Iso-frequency contour of Fig. 4 with possible wave vectors κ = (κ →x , κ ′ y ) T in (a) and κ = (κ ← x , κ ′ y ) T in (b), both shown in green, with corresponding Poynting vector J = (J x , J y ) T and J = (−J x , J y ) T normal to the contour shown in blue.(See next section for the choice of notation).The resulting real components of the eigenfunction solutions ψ ⇋ nm (z; k, κ x , κ y ) for the values of the wave vectors used in (a) and (b) are also shown.

Figure 7 .
Figure 7. (a) Shows the cross resonator plugged into the metamaterials unit cell with scattering matrix defined as ŜŴ cross (k; l x , l y ) in (25).Both (b) and (c) show the resulting band for different values of l x and l y .Also plotted in the κ x , κ y plane are the iso-frequency contours of the second band for various values of k.Plot (b) is for l x = 1 , l y = 1 .Plot (c) is for l x = 0 , l y = 1.

Figure 8 .
Figure 8. Boundary region between metamaterials 1 and 2, understood as all right(r) edges for n = −1 and all left (l) edges for n = 0 .Here, wave scattering from the boundary is divided into event I and II, where r p and t p represent reflection and transmission amplitudes for event p = I or II.

Figure 9 .
Figure 9.The iso-frequency contours for k = 1/l for a scattering matrix with ŜŴ 1 := ŜŴ point and ŜŴ 2 := ŜŴ cross (k; 0, 1) representing the two metamaterials.The horizontal dashed grey line represents a single value of κ y = κ ′ y .Intersections with contour lines give the corresponding values of κ ⇋ j,x for j = 1 and 2.

Figure 11 .
Figure 11.Top: Iso-frequency contours for k = 1/l for a scattering matrix with ŜŴ 1 := ŜŴ point and ŜŴ 2 := ŜŴ cross (k; 0, 4.45) .Bottom: The resulting real component of the Gaussian beam nm constructed from the full wave field nm across the two materials with incident wave amplitudes A 1 = 1 and B 2 = 0 .The focal point of the Gaussian Beam is set to n ′ = −60 with tilt given by κ ′ y = 1/l.

Figure 12 .
Figure 12.A system of N layered metamaterials with wave amplitudes, A j and B j , noted at each interface.

Figure 13 .
Figure13.Top: The iso-frequency contours of three metamaterials with properties defined by scattering matrices ŜŴ point , ŜŴ cross (k; 0, 4.45) and ŜŴ point for k = 1/l .Bottom: The real component of a Gaussian beam nm incident from material 1 constructed from the full wave field nm for an incident wave from metamaterial 1, A 1 = 1 , and no incident wave from material 3, B 3 = 0.