Issue 
EPJ Applied Metamaterials
Volume 2, 2015
Advanced Metamaterials in Microwaves, Optics and Mechanics



Article Number  11  
Number of page(s)  7  
DOI  https://doi.org/10.1051/epjam/2015007  
Published online  26 January 2016 
https://doi.org/10.1051/epjam/2015007
Research Article
Optic axisdriven new horizons for hyperbolic metamaterials
^{1}
Joule Physics Laboratory, Materials & Physics Research Centre, University of Salford, Greater Manchester
M5 4WT, UK
^{2}
Blackett Laboratory, Imperial College, Prince Consort Road, London
SW7 2AZ, UK
^{*} email: a.d.boardman@salford.ac.uk
Received:
1
September
2015
Accepted:
28
October
2015
Published online: 26 January 2016
The broad assertion here is that the current hyperbolic metamaterial world is only partially served by investigations that incorporate only some limited version of anisotropy. Even modest deviations of the optic axis from the main propagation axis lead to new phase shifts, which not only compete with those created by absorption but end up dominating them. Some progress has been attempted in the literature by introducing the terms “asymmetric hyperbolic media”, but it appears that this kind of asymmetry only involves an optic axis at an angle to the interface of a uniaxial crystal. From a device point of view, many new prospects should appear and the outcomes of the investigations presented here yield a new general theory. It is emphasised that the orientation of the optic axis is a significant determinant in the resulting optical properties. Whereas for conventional anisotropic waveguides homogeneous propagating waves occur over a limited range of angular dispositions of the optic axis it is shown that for a hyperbolic guide a critical angular setting exists, above which the guided waves are always homogeneous. This has significant implications for metawaveguide designs. The resulting structures are more tolerant to optic axis misalignment.
Key words: Metamaterial / Hyperbolic / Uniaxial / Optic axis
© A.D. Boardman et al., published by EDP Sciences, 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The interest in calcite (Iceland spar) dates back to 1669 and the work of Bartholin [1]. This means that the theory of propagation of electromagnetic waves in such uniaxial anisotropic crystals has naturally been of interest for some considerable time. Indeed, the early investigations are extremely useful in forming a background and vital points of contact for the creation of new horizons in the global surge of interest in what has become known as hyperbolic [2] nonmagnetic metamaterial. The current hyperbolic metamaterial world is only partially served by investigations that incorporate some limited version of anisotropy. In order to develop new pathways, it should first of all be noted that readily available uniaxial nonmetamaterial crystals, like calcite, have a direction around which rotational symmetry is apparent. This direction is called the optic axis and, given a rectangular coordinate system, it is usually set to lie along the zaxis, which in the model presented here is going to be the direction of wave propagation. Adopting this step leads to the permittivity tensor [3] emerging in diagonal form, possessing equal x and y components that are unequal to the component associated with the zdirection. Hyperbolic metamaterials have the same x and y tensor properties but possess a negative zcomponent. A completely general approach to the constructive manipulation of the optic axis position in uniaxial crystals is required if a proper description of any wave propagation within them is to be achieved. In this context, as will be emphasized here, there is a need to recognise clearly the dramatic role of anisotropy in uniaxial metamaterial investigations, as was pointed out some years ago [4].
In order to set up a full investigation of how the orientation of the optic axis affects wave propagation in hyperbolic metamaterials, it is efficient to focus upon two broad areas. Area one includes the interrogation of a slab of metamaterial, involving the transmission of light through an interface, or a number of interfaces. Area two invokes the paramount role of the optic axis when guided waves are being generated, and it is this case that will be studied here, where it will be shown that for a hyperbolic film, in contrast to the usual nonmetamaterial anisotropic guides, substantial changes in the character of guided waves occur, presenting the opportunity for novel device design.
In general, there are offdiagonal permittivity tensor elements to be considered when controlling the optic axis direction. These elements can often be brought in unintentionally at experimental level when trying to compare measurements with simulations. This is emphasised in the work by Starodubtsev [4], in which he writes
Thus, even modest deviations of the optic axis from the main zaxis leads to new phase shifts, which not only compete with those created by absorption but may end up dominating them. Some progress has been attempted by introducing the term “asymmetric hyperbolic media” [5, 6], but it appears that this kind of asymmetry only involves an optic axis at a certain angle to the interface of a uniaxial crystal. Here, a full investigation is made of the effect that the optic axis direction has on the existence of pure hybrid homogeneous waves in a wave guide consisting of a hyperbolic film, an asymmetric substrate with all positive permittivity tensor elements, and an air cover.“even small optical anisotropy is a much greater enemy of the perfect lens than absorption”.
2 The optic axis
The historic, classic work of Knoesen, Gaylord and Moharam [7] on nonmetamaterial anisotropic uniaxial waveguides introduces the geometrical concept of an optic axis in the context of a slab waveguide. It is an attractive anisotropic design, which is why it is adopted here in order to illustrate the new metamaterial concepts to be introduced in this investigation. It does not limit the generality of the overall conclusions and, for the moment, all three components of the wave number, related to right angle Cartesian coordinates will be considered.
The classic diagonal form of the permittivity tensor for a uniaxial, nonmagnetic crystal when the optic axis lies along the zaxis is given by(1)
Referring to Figure 1, the general direction of the optic axis is described through the angles θ and ϕ, where θ is the angle that the optic axis makes with the zaxis, and ϕ is the angle between the projection of the optic axis onto the xyplane measured from the xaxis (Figure 1).
Figure 1. For a rectangular coordinate system (x, y, z), the general position of the optic axis in a uniaxial crystal can be set at an angle θ to the zaxis. This axis then has a projection on the (x, y) plane, which is at an angle ϕ to the xaxis. 
In order to take into account a variation in the orientation of the optic axis, the transformation matrixcan be used to find the elements of the general permittivity tensor,(2)in the usual way, using . In the lossless cases being investigated, the tensor is Hermitian so that . Thus, for real elements, the matrix is symmetric, ε_{ij} = ε_{ji} and the elements are given by [8](3)
In special cases, for example when ϕ = 0° and the optic axis lies in xzplane [9], or when θ = 90° and the optic axis lies in the xyplane [10], the relevant permittivity tensor can be found using the appropriate values of θ and ϕ in equation (3).
For wave propagation, a useful characteristic is the dimensionless wave number with components , where i = x, y or z, c is the velocity of light in freespace and ω is the angular frequency. For the case when all the diagonal tensor elements are positive, two surfaces in (κ_{x}, κ_{y}, κ_{z}) space will appear for a given frequency. One is that of an ellipsoid and the other is that of a sphere. The former is classically associated with extraordinary modes where, in general, the power flow and phase velocity direction are not parallel. The sphere is associated with ordinary modes where the power flow and the phase velocity directions are the same. For a metamaterial described by a diagonal dielectric tensor with two equal positive elements, and one negative element, the ordinary wave sphere still emerges but the ellipsoid opens up into a hyperboloid. Figure 2 demonstrates this for the case when the optic axis lies along the zaxis [11].
Figure 2. Crosssections through the (κ_{x,}κ_{z})plane in wavenumber space when κ_{y} = 0 and the optic axis lies along the κ_{z}axis. The figure shows the transition from anisotropic (ellipse) to hyperbolic as ε_{3} changes in incremental steps from a positive to a negative value. The element ε_{1} = ε_{2} = 4.0. 
Following the magnificent lead, set many years ago by the rigorous and pioneering general work of Knoesen et al. [7], and some pioneering work on beam propagation by Fleck and Feit [12], the discussion presented below will focus upon a lossfree, planar, linear uniaxial waveguide (Figure 3), where the subscripts s, f and c are used to represent the substrate, film and cover respectively.
Figure 3. Planar, stepindex, waveguide with guided waves propagating along z, for which the central film (f) is a hyperbolic metamaterial. The cover (c) and substrate (s) can be purely isotropic, or anisotropic, with positive tensor elements. 
In this investigation the planar guide has a hyperbolic metamaterial central film with ε_{3f} < 0, an anisotropic substrate with ε_{3s} > 0, and an isotropic cover, assumed for the moment to be air. The adoption of losses can be introduced at a later stage of any device design but it is not a paramount issue here since hyperbolic metamaterials, unlike doublenegative metamaterials, are not resonant. The model here has an optic axis that is permitted to have any orientation, but is the same throughout for the film, the substrate and the cover [7]. Uncoupling of the TE and TM modes occurs when the angle ϕ = 0°. Otherwise, the waves are hybrid and will be highlighted in terms of pure guided homogeneous modes, pure guided inhomogeneous modes, or leaky guided modes [7]. A clarification of the homogeneous modes is that these waves have constant planar phase fronts, and are totally guided with constant amplitudes. Inhomogeneous waves, on the other hand, are characterised by planar phase fronts with nonconstant amplitudes. Leaky modes are not totally guided by a slab waveguide. The main thrust in this paper is in the effect that the orientation of the optic axis has on the existence of pure hybrid homogeneous modes. The term “anisotropic” will be used when all the elements of the permittivity tensor are positive, and, specifically, when ε_{1} = ε_{2} > 0; ε_{3} > 0 ≠ ε_{1}. The expression “hyperbolic” will be used to mean anisotropic, with ε_{1} = ε_{2} > 0 and ε_{3} < 0.
In the absence of attenuation, the values of κ_{z} (the component of the wavenumber in the propagating direction) are real, and, in the case of an anisotropic film, are restricted by the permittivity of the isotropic cover, and the permittivity tensors of the anisotropic substrate and the hyperbolic film. In the anisotropic film case, the maximum value that κ_{z} is can take is √ε_{1f}. However, under certain orientations of the optic axis, to be determined later, this upper boundary no longer applies when the film is hyperbolic.
If the optic axis lies in the xzplane, then the range of values of θ over which the important homogeneous waves propagate, is restricted to a specific range for both anisotropic and hyperbolic guides. For an anisotropic film, if the optic axis is also rotated so that angle ϕ has a positive value, then the range of θ over which homogeneous waves propagate decreases. In contrast, when using a hyperbolic metamaterialbased guide, increasing ϕ from 0° results in an increase in the range of θ over which pure hybrid homogeneous modes propagate until, at a certain critical value of ϕ, these modes will propagate for all angles of θ between 0° and 90°. At angles of ϕ greater than this critical angle, modes are possible with but these modes are hybrid inhomogeneous waves. This will be discussed in more detail later.
3 Wave propagation and design windows
First of all, consider a bulk nonmagnetic uniaxial crystal, with a permittivity tensor ε_{ij}. This material will support plane waves of the form A e ^{i(k⋅rωt)}, where the complex wave vector amplitude is A, and wave vector and position vector are, respectively, k = (k_{x}, k_{y}, k_{z}) and r = (x, y, z). Maxwell’s equations then yield the following equation for the electric field components (E_{x}, E_{y}, E_{z}).(4)or, in terms of the dimensionless wave vector κ_{i} defined above(5)
In the planar waveguide (Figure 2), as is usual, it is assumed that there is no variation of the fields in the ydirection so that setting κ_{y} = 0 in equation (5) yields the biquadratic(6)which, with judicious use of equations (3), reduces towhich factorises to give(7)
These equations show how the elliptic case, with only positive diagonal tensor elements, turns into the hyperbolic case as ε_{3} becomes negative (see for example Figure 2). At this stage, it is convenient to introduce the subscripts o and e to label the ordinary and extraordinary mode contributions respectively. This step results in the solutions(8)
For a waveguide, these equations give the κ_{x}/κ_{z} relationship for each layer. Since the interest here is in waves that do not attenuate in zdirection, κ_{z} must be real. Thus it is clear from (8) that for(9)the solutions for κ_{x} are real, and for(10)the solution for the ordinary mode will be imaginary, and the solution for the extraordinary mode will be complex. In the guide, these results will determine whether the modes are homogeneous, inhomogeneous or leaky modes. For propagation along the guide, fields that decay away from the surfaces into the semiinfinite media of the substrate and cover are sought. If either or both of these layers are isotropic, then, for guided waves, the first of (10) must be satisfied, and exponential decay will occur. Using the second of (10), for an anisotropic substrate, the decay will not be exponential due to the complex values of the solution for the extraordinary modes given by (8). In this case, solutions are sought which determine that the electric and magnetic fields decay to zero as x → ±∞.
For the specific cases discussed below it is assumed that that the cladding is air and in the anisotropic case ε_{1f} > ε_{1s} > ε_{3f} > ε_{3s} > ε_{1c}. In the hyperbolic film case ε_{1f} > ε_{1s} > ε_{3f} > ε_{3s} > ε_{1c} and ε_{3f} < 0. Thus, for the constituents of the slab waveguide shown in Figure 3, the necessary conditions for pure hybrid homogeneous modes are:
for the film:(11) for the anisotropic substrate:(12) and for the isotropic cover:(13)
Using (11) and (12) and the definitions in (2), pure homogeneous hybrid waves with real κ_{z} can be captured, in terms of θ and ϕ, by the following inequalities.
For the extraordinary mode:(14) and for the ordinary mode:(15)
For a substrate with parameters ε_{1s} = 3.3856, ε_{3s} = 2.25, and a thin film guide with parameters ε_{1f} = 4, ε_{3f} = 2.7889, the effect of changing the value of ϕ is very small as can be seen in Figure 4. Equations (11) and (12) give the blue curves when ϕ = 0°, and the red curves when ϕ = 90°. Applying equations (14) and (15) produces the areas in blue from which it can be seen that the change in ϕ has little effect on the range of θ over which pure homogeneous modes occur.
Figure 4. The change in the regions defining pure hybrid homogeneous waves for the anisotropic guide with allpositive diagonal tensor elements as the angles of orientation are changed. Values used: ε_{1s} = 3.3856, ε_{3s} = 2.25, ε_{1f} = 4 and ε_{3f} = 2.7889. The blue areas show how changing the value of ϕ from ϕ = 0° (the blue curves) to ϕ = 90° (the red curves) has little effect on the range of θ over which pure hybrid homogeneous modes are possible. 
Now the inequalities given in (14) and (15), permitting the existence of various types of waves, are quite general. For a hyperbolic film, where ε_{3f} < 0, the inequality on the right hand side of (14) becomes(16)
The profound effect of this change can be demonstrated, in detail, by looking analytically at the changes in the regions where the hybrid homogeneous, inhomogeneous and leaky modes can exist. In order to illustrate the specific effect of a negative value for ε_{3f}, the value ε_{3f} = −2.7889 will be selected, and the other permittivity values will remain the same. The effect of this change in sign is shown in Figure 5, and contrast greatly from the results shown in Figure 4.
Figure 5. The figure demonstrates the increasing range in θ over which homogeneous waves occur in a hyperbolic film with ε_{3f} < 0, as ϕ increases towards a certain critical value. The dark blue areas show the regions over which pure hybrid homogeneous guided wave are possible. 
From Figure 5, it is apparent that the pure hybrid homogeneous mode region depends on the point where the curves given by(16a)for the different values of ϕ cut the line . These points are given when θ = θ_{ϕ}, where, for ε_{3f} > 0,(17)which, as is apparent in Figure 4, decreases, but not very significantly, as ϕ changes from 0° to 90°. For the hyperbolic case, for which ε_{3f} < 0, (17) can be written(18)
There is a dramatic difference in the two cases. When ϕ = 0°, the ranges of θ over which pure hybrid homogeneous waves propagate extend to(19)and(20)
Comparing (19) and (20) it is clear that when ε_{3f} > 0° and ϕ = 0°, so that the optic axis is in the xzplane, the range for the anisotropic guide is wider than that for the hyperbolic guide. In the former case, as ϕ increases from zero, it is easy to see from (18) that for ε_{3f} > 0, as , the range of θ over which pure hybrid homogeneous modes can exist decreases, but only by a relatively small amount (Figures 4 and 6). In contrast to this case, when ε_{3f} < 0 and ϕ increases from zero, the range increases rapidly until, as ϕ approaches a certain critical value ϕ_{crit}, given by(21)pure hybrid homogeneous modes can exist for all values of the angle θ from 0° to 90° (Fgures 5 and 6). Note that as ϕ approaches the value ϕ_{crit}, the change in the value of θ_{max} becomes very rapid so that very small changes in the orientation of the optic axis has a profound effect, which could have important consequences for the design and manufacture of hyperbolic guides.
Figure 6. The graphs show the range of θ over which homogeneous waves can propagate for ε_{1s} = 3.3856, ε_{3s} = 2.25, ε_{1f} = 4, ε_{3f} = 2.7889. The range extends from θ = 0° up to each of the curves in the figure. As ϕ increases, there is little change when ε_{3f} > 0. When ε_{3f }< 0, there is a rapid increase in the range until, at a particular value ϕ = ϕ_{crit}, homogeneous waves can propagate at all angles of θ up to 90°. With the parameters given in the text ϕ_{crit} = 39.9° (shown by the dotted line). The upper right hand curve shows the decreasing range of θ when φ increases beyond ϕ = ϕ_{crit}. 
It is apparent from Figure 6 that for the case where ε_{3f} > 0, the critical value of ϕ has no meaning. It should be noted here that the critical value also depends crucially on the value of ε_{3f} as shown in Figure 7.
Figure 7. The dependence of the critical angle (ϕ_{c} in the figure) on the value of ε_{3f} < 0, above which the wave number κ_{z} is no longer limited by the value of ε_{3f} when this is negative. When ε_{3f} > 0 the value of must be less than ε_{3f} for propagation of hybrid homogeneous modes. If ε_{3f} < 0, there will be an angle ϕ_{crit} above which there will be propagation of inhomogeneous modes. The arrow in the diagram shows the critical angle when ε_{3f} = −2.7889. 
A second important major difference between the usual anisotropic case and the hyperbolic case is that in the latter modes can propagate at values of . However, in this region the condition for the ordinary wave in equation (2), namely, is no longer satisfied, and the waves are inhomogeneous. These waves and the associated leaky waves will be examined in detail in a future paper.
As the angle ϕ increases beyond the value of ϕ_{c}, the range of θ over which pure hybrid homogeneous waves propagate decreases again and is given by:(22)
Figure 8 demonstrates this decreasing range of range of θ over which these modes propagate as ϕ increases from a value just above ϕ_{c} to the angle ϕ = 90°.
Figure 8. Above critical angle ϕ_{c}, the normalised wave number κ_{z} is no longer limited by the value of ε_{1f}, but, in the regions above ε_{1f} and below the solid curves, the modes are inhomogeneous. The blue areas are the regions where hybrid homogeneous modes occur. The figure shows the decreasing range of θ over which these modes are possible as ϕ increases beyond the value of ϕ_{c}. The dashed line is the result for the anisotropic substrate where as ϕ changes the change in the graph is undetectable on the scale of the figure. 
The above analysis shows, for a particular data set, how varying the position of the optic axis leads to interesting conclusions about the regions of existence of different forms of waves. In particular it shows how varying the orientation of the optic axis creates, or closes, windows of opportunity to propagate the kind of waves that are certain to be part of upcoming devices. In order for a metadevice to operate in a homogeneous wave window it is important first of all to conduct an investigation into how the orientation angle ϕ can be made to influence the wave outcomes. The analysis also introduces the existence of inhomogeneous waves for values of the normalised wavenumber that lie above the maximum possibilities for the anisotropic guide with a permittivity tensor with all positive elements.
4 Conclusions
As opposed to what is revealed in the limited literature on optic axis orientations in metamaterials, the general theory given here shows the advantages of orienting the optic axis in carefully defined directions.
Through examining the dispersion relations for hyperbolic waveguides, this paper emphasises that the orientation of the optic axes is a significant determinant in the resulting optical properties. In fact, whereas for conventional anisotropic waveguides, homogeneous propagating waves occur over a limited range of the polar angle θ, a range that varies only a little with the azimuthal angle ϕ, we have shown that for a hyperbolic guide, a critical angle ϕ_{c} exists, at which value the guided waves are homogeneous for all values of θ between 0° and 90°. This has significant implications for metawaveguide designs in that it implies that, in the hyperbolic case, generation of propagating modes can be more easily achieved, and the resulting structures are more tolerant to optic axis misalignment.
The high sensitivity of the mode propagation constant to azimuthal orientation of the optic axis could also be transformative for metaguided wave technology. Small external perturbations, such as the application of mechanical pressure, or temperature, can result in changes to the effective orientation of the guide/substrate axes. These in turn, will result in modes being cut off, or allowed to propagate, opening up potential for sensor applications. Indeed, electrooptic control over the orientation of the optic axis in a hyperbolic guide could lead the way to a communications technology based on active mode control.
Also, we note the implication of our work for transformation designs that exploit the orientation of the optic axis of anisotropic media, as described in [13]. In that work, pixels of anisotropic media, where the orientation of the optic axis alternated between pixels, were shown to produce an effective transmission medium capable of a variety of transformation devices, including a wave expander and a source shifter. Our work suggests that a greater range of applications is possible through the use of pixelated hyperbolic media.
References
 E. Bartholin, Experimenta crystalli islandici disdiaclastici quibus mira & insolita refractio detegitur, Copenhagen (“Hafniæ”), Denmark, 1669. (In the text)
 L. Ferrari, C. Wu, D. Lepage, X. Zhang, Z. Liu, Progress in Quantum Electronics 40 (2015) 1. [CrossRef] (In the text)
 D.C. Kay, in Schaum’s outline of theory and problems of tensor calculus, McGrawHill, 1988. (In the text)
 E. Starodubtsev, Proc. SPIE 7353 (2009) Metamaterials IV 7353OA. (In the text)
 S.M. Hashemi, I.S. Nefedov, M. Soleimani, Photonics Letters of Poland 5 (2013) 72. (In the text)
 I.S. Nefedov, C.A. Valagiannopoulos, S.M. Hashemi, E.I. Nefedov, Scientific Reports 3 (2013) 2662. [CrossRef] (In the text)
 A. Knoesen, T.K. Gaylord, M.G. Moharam, Journal of Lightwave Technology 6 (1988) 1083. [CrossRef] (In the text)
 D.P. GiaRusso, J.H. Harris, Journal of the Optical Society A 63 (1973) 138. [CrossRef] (In the text)
 D. Marcuse, IEEE Journal of Quantum Electronics QE14 (1978) 736. [CrossRef] (In the text)
 D. Marcuse, I.P. Kaminow, IEEE Journal of Quantum Electronics QE15 (1979) 736. (In the text)
 K.E. Ballantine, P.R. Eastham, J.F. Donegan, Physical Review A 90 (2014) 013803. [CrossRef] (In the text)
 J.A. Fleck, M.D. Feit, Journal of the Optical Society A 73 (1983) 920. [CrossRef] (In the text)
 Z. Liang, X. Jiang, F. Miao, S. Guenneau, J. Li, New Journal of Physics 14 (2012) 103042. [CrossRef] (In the text)
Cite this article as: Boardman AD, Egan P & McCall M: Optic axisdriven new horizons for hyperbolic metamaterials. EPJ Appl. Metamat. 2015, 2, 11.
All Figures
Figure 1. For a rectangular coordinate system (x, y, z), the general position of the optic axis in a uniaxial crystal can be set at an angle θ to the zaxis. This axis then has a projection on the (x, y) plane, which is at an angle ϕ to the xaxis. 

In the text 
Figure 2. Crosssections through the (κ_{x,}κ_{z})plane in wavenumber space when κ_{y} = 0 and the optic axis lies along the κ_{z}axis. The figure shows the transition from anisotropic (ellipse) to hyperbolic as ε_{3} changes in incremental steps from a positive to a negative value. The element ε_{1} = ε_{2} = 4.0. 

In the text 
Figure 3. Planar, stepindex, waveguide with guided waves propagating along z, for which the central film (f) is a hyperbolic metamaterial. The cover (c) and substrate (s) can be purely isotropic, or anisotropic, with positive tensor elements. 

In the text 
Figure 4. The change in the regions defining pure hybrid homogeneous waves for the anisotropic guide with allpositive diagonal tensor elements as the angles of orientation are changed. Values used: ε_{1s} = 3.3856, ε_{3s} = 2.25, ε_{1f} = 4 and ε_{3f} = 2.7889. The blue areas show how changing the value of ϕ from ϕ = 0° (the blue curves) to ϕ = 90° (the red curves) has little effect on the range of θ over which pure hybrid homogeneous modes are possible. 

In the text 
Figure 5. The figure demonstrates the increasing range in θ over which homogeneous waves occur in a hyperbolic film with ε_{3f} < 0, as ϕ increases towards a certain critical value. The dark blue areas show the regions over which pure hybrid homogeneous guided wave are possible. 

In the text 
Figure 6. The graphs show the range of θ over which homogeneous waves can propagate for ε_{1s} = 3.3856, ε_{3s} = 2.25, ε_{1f} = 4, ε_{3f} = 2.7889. The range extends from θ = 0° up to each of the curves in the figure. As ϕ increases, there is little change when ε_{3f} > 0. When ε_{3f }< 0, there is a rapid increase in the range until, at a particular value ϕ = ϕ_{crit}, homogeneous waves can propagate at all angles of θ up to 90°. With the parameters given in the text ϕ_{crit} = 39.9° (shown by the dotted line). The upper right hand curve shows the decreasing range of θ when φ increases beyond ϕ = ϕ_{crit}. 

In the text 
Figure 7. The dependence of the critical angle (ϕ_{c} in the figure) on the value of ε_{3f} < 0, above which the wave number κ_{z} is no longer limited by the value of ε_{3f} when this is negative. When ε_{3f} > 0 the value of must be less than ε_{3f} for propagation of hybrid homogeneous modes. If ε_{3f} < 0, there will be an angle ϕ_{crit} above which there will be propagation of inhomogeneous modes. The arrow in the diagram shows the critical angle when ε_{3f} = −2.7889. 

In the text 
Figure 8. Above critical angle ϕ_{c}, the normalised wave number κ_{z} is no longer limited by the value of ε_{1f}, but, in the regions above ε_{1f} and below the solid curves, the modes are inhomogeneous. The blue areas are the regions where hybrid homogeneous modes occur. The figure shows the decreasing range of θ over which these modes are possible as ϕ increases beyond the value of ϕ_{c}. The dashed line is the result for the anisotropic substrate where as ϕ changes the change in the graph is undetectable on the scale of the figure. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.