Graphene-based extremely wide-angle tunable metamaterial absorber.

We investigate the absorption properties of graphene-based anisotropic metamaterial structures where the metamaterial layer possesses an electromagnetic response corresponding to a near-zero permittivity. We find that through analytical and numerical studies, near perfect absorption arises over an unusually broad range of beam incidence angles. Due to the presence of graphene, the absorption is tunable via a gate voltage, providing dynamic control of the energy transmission. We show that this strongly enhanced absorption arises due to a coupling between light and a fast wave-mode propagating along the graphene/metamaterial hybrid.

With the refinement and advancement of fabrication techniques and theoretical methods, it is now possible to successfully create anisotropic metamaterial structures that have at least one vanishing principal component of the permeability or permittivity ε tensors. These artificial hybrid structures have been shown to allow propagation of electromagnetic (EM) waves through narrow1 or mismatched2 channels with little reflection, and offer alternatives in wavefront shaping34, beam steering, and EM energy control5 previously unseen in conventional systems. Typical experimental platforms involve various geometrical arrays of metallic rods67 or combinations of metallic and dielectric multilayers that have tunable system parameters to achieve the desired near-zero effective EM response8. In addition to these platforms, other metamaterial architectures have been designed that manipulate the EM response to significantly reduce loss. These include, layered metal/dielectric structures910, transparent conducting materials111213, superconducting nanoparticles1415, and plasmonic nanoshells with quantum dot cores16. While there are a number of mechanisms to dynamically tune metamaterials17, it would be preferable to have more robust, real-time tunability of their effective electric (E) or magnetic (H) responses.

These developments have led to a significant number of recent works devoted to studying graphene as a component in hybrid metamaterial designs18. Enhanced absorbance in a hybrid graphene-metamaterial structure on a reflective substrate was proposed in ref. 26. In addition to its intrinsic 2D geometrical properties, graphene affords extensive control of its conductivity σ due to its gate-tunable Fermi energy E. When EM waves interact with a graphene sheet, the gate voltage can tune the amount of energy that gets absorbed, depending on the relevant scattering processes, temperature, relaxation time, frequency of the incident beam and its polarization.

By placing graphene in contact with an anisotropic metamaterial that has vanishing principal components of the permittivity and permeability tensors over the frequencies of interest, one could expect the interplay between the two samples to present additional possibilities regarding the transmission of the light. More specifically, we pose the question if it is possible to utilize the dynamic and tunable conductivity of graphene to yield greater control over the way energy of EM waves is absorbed or reflected compared to when the materials are in isolation. Motivated by recent experimental and theoretical developments, we seek to establish if one may obtain tunable wide angle absorption in hybrid graphene/metamaterial layers where the tensor ε describing the anisotropic response has components that are near zero1.

To this end, we consider here the absorption characteristics of a layered system, comprised of monolayer graphene overlaying a metamaterial slab that is connected to a perfectly conducting metal. In accordance with the discussion above, we will consider metamaterial structures that can possess an ε near-zero (ENZ) response, where the metamaterial has a vanishingly small value of the real part of ε. In contrast to numerous past works in the RF and visible light regime, we focus here on operational wavelengths that correspond to the less exploited THz regime, where many new technologies are emerging. To accurately determine the absorption as a function of the geometrical and material parameters that describe graphene and the metamaterial, Maxwell’s equations are solved using a transfer matrix approach that accounts for the presence of graphene through boundary conditions imposed on the magnetic field. We discuss the behavior of the absorption in terms of the energy flow throughout the entire structure.

Our main result is that by tuning E, we are able to determine experimentally relevant regimes where near perfect absorption can occur over an extremely wide range of incident angles for an ENZ metamaterial layer that has a small component of the permittivity tensor. We show that the enhanced absorption occurs due to coherent perfect absorption which couples light to a fast wave propagating along the graphene/ENZ hybrid.

Theory

The principal coordinate axes are chosen so that both the permittivity and permeability tensor in the anisotropic medium are both diagonal. A diagram of our setup is shown in Fig. 1. For concreteness, we focus on p-polarized incident light (E, E, B) ≠ 0. In the anisotropic medium (region 2), the nonzero components to the tensors are (ε, ε, μ). It is backed by a perfect electric conductor (PEC) in region 1, where no electric or magnetic fields exist. The finite surface conductivity σ of graphene will be effectively taken into account in the boundary condition for the magnetic field at the interface between air (region 3) and the anisotropic medium, where26,

Figure 1

Considered setup in this work.

An EM wave in vacuum (region 3) is incident upon a graphene layer coating an anisotropic ENZ metamaterial (region 2). A perfectly conducting metal backs the entire structure (region 1).

When E ≫ kT, ħω, we have,

Thus, at a temperature T of 300 K, and THz frequencies, we have kT = 25 meV, and ħω ~ 1 meV, meaning doping levels should exceed E = 100 meV when using Eq. (2). Here, E is the Fermi energy, ω is the frequency of the EM wave, and trel is the relaxation time. For frequencies in the THz region and lower, Eq. (2) reflects that the interband part to the conductivity is negligible compared its the intraband part, and consequently to the Drude-like behavior for σ.

We now consider Maxwell’s equation (in SI units) in the incident air region and in the anisotropic medium, both having no free charges or currents. For p-polarized incident light and diagonal tensors, we obtain the following set of equations determining the components of the fields when the time-dependence is harmonic (e−iω):

These equations are solved with a plane-wave ansatz for the field components where and k is conserved due to translational invariance. The Methods section contains additional details of the numerical and analytical solutions of the above equations.

Once the equations have been solved, one may compute the absorption from the scattering coefficients: A2 and B2 are associated with the transmitted and reflected electric field wave in region 2 whereas B3 ≡ r is the coefficient for the reflected wave in region 3. Since there is no transmission into the PEC layer, the absorbance can be computed from the incident wave after subtracting the reflected part: . The reflection coefficient r is obtained as

where we have defined

The coefficients for the other waves are given as: , and B2 = −A2. We are now interested to see how the presence of the graphene layer modifies the absorption properties of the metamaterial, and in particular determine the influence of the gate-controlled Fermi level which is tunable in situ. The analytical expression for the reflection coefficient r allows us to estimate which order of magnitude E should have in order to have any influence on the results. If the terms proportional to the conductivity σ are to be comparable in magnitude with the other terms, we see that e.g. εω ~ kσ needs to be satisfied. Since we are considering an EM wave incident from air, ε = ε = ε0. Inserting the expression for σ, one finds that the above equation corresponds to a Fermi level in graphene of

For grazing incidence θ ~ π/2, the graphene layer thus has no influence for any realistic values of E. We also note that when σ → ∞, r → 1 such that the absorption vanishes, . In order to gain information about the spatial distribution of the power flow in the structure, we also compute the Poynting vector  =  × . We write  = 0e−iω and  = 0e−iω. For future use, we define the normalization constant for the power flow as the time-averaged incident power at z = τ: S0 = ε0c/(2 cos θ).

By appropriately tuning E, graphene-dielectric stacks can exhibit a wavevector dispersion that traces out a hyperbolic contour at a given frequency2123. Such hyperbolic metamaterial (HMM) systems are described by a permittivity tensor whose principle components are opposite in sign49, and allow for the propagation of evanescent modes in subwavelength imaging29, and can be designed to exhibit substantial near-field absorption22. Such structures are often nonresonant and possess low loss50. There have been a variety of HMM systems that have been fabricated52, including semiconductor hybrids, metallic layers52, and silver nanowires54. Gate tunable HMM systems have also been made that allow control of their optical response53. Here, we will determine how the absorption properties of HMM structures are modified by the presence of a graphene layer over a broad range of angles θ and Fermi energies E. Specifically, we focus on gate-tunable HMM structures in the ENZ limit, whereby the tensor components ε and ε are opposite in sign. We set μ = μ0, and the relaxation time in graphene to τ = 10−13 s. For incident light in the THz regime, with E = 100 meV and λ = 1.6 mm, this leads to a surface conductivity of graphene σ = (1.2 × 10−3 + 1.37i × 10−4) A ⋅ C/J.

We first consider the standard Drude-like frequency response in the metamaterial, such that

Here, α = λ/λ and f = 0.02, where λ is the characteristic wavelength, with λ = 1.6 mm. Designing a metamaterial within effective medium theory to possess the prescribed dielectric response in the THz regime can be challenging with current fabrication and materials capabilities. It may be possible however to design a multilayer dielectric structure or rod lattice55 with polaritonic components to achieve the desired HMM response since the permittivity of polaritonic materials mimic metals in the optical regime5657. Therefore construction of an anisotropic HMM that has a longitudinal ENZ response should be possible at THz frequencies.

It will also be of interest to see how the absorption properties depend on the sign of ε, determined by the frequency of the incident wave, and also its magnitude. We set , so that . If we take the limit ε → 0, the reflection coefficient simplifies to

In the absence of graphene (σ = 0), we have zero absorption since |r|2 = 1. However, when graphene is present (σ ≠ 0) an interesting opportunity arises. It is seen from the above equation that r = 0, meaning perfect absorption when the following condition is satisfied:

In the low THz regime and below (ω ≪ 10 THz), we have ω ≪ τ and full absorption takes place for angles satisfying

Although this equation represents the criterium for 100% absorption, we show in the next section that the absorption in fact is close to 100% for an extremely wide range of incident angles.

Results and Discussion

To establish how perfect absorption can be achieved, we discuss how the incident beam can couple to the intrinsic EM modes of the system. Our approach involves finding the poles of the reflection coefficient, Eq. (8), by setting the denominator equal to zero. To find the complete set of poles including the bound and leaky wave contributions, it becomes necessary to extend k = β + iα to the complex plane and solve the pole dispersion equation (see the Methods section) for β and α (at a given ω). This dispersion equation admits four types of solutions for the complex wavevector58. Of these solutions, only those which correspond to perfect absorption modes are retained, since not all solutions correspond to coupling of the incident beam with the structure. The perfect absorption modes of interest are those which are activated via phase-matching of the incident plane wave to solutions corresponding to k/k0 < 1, and ℑm(k) = 0. These modes represent a coherent superposition of waves that propagate without loss along the surface. Indeed, only when the incident wave is incident at the perfect absorption angle, , does the direction of energy flow in the air region, (≡ arctan(S/S), precisely coincide with the incident angle . Deviations from introduce reflected waves which shift the net energy flow in the air region. Additionally, for r to be single-valued, the appropriate branch cuts at the points ±k0 need to be made for k. When examining the perfect absorption phenomena, we utilize a simple yet useful property of the reflection coefficient: Eq. (8) reveals that interchanging the direction of the incident and reflected waves, transforms r into its reciprocal:

with k = k0 cos θ. This property simplifies the solution process since finding a complex pole, related to −k, is equivalent to finding the perfect absorption modes for k5960. We underline that the absorption occurring in this way is an interference phenomenon and not based on losses.

We now introduce more possibilities for enhanced absorption, and investigate anisotropic systems where the permittivity components perpendicular and parallel to the interfaces are of opposite sign. The corresponding HMM dispersion relations are then identified as either type I if ε > 0 and ε < 0, or type II if ε < 0 and ε > 0. The operating wavelength of the incident beam is tuned around λ to control the sign of ε. The component of ε parallel to the interface always has the same magnitude, but can differ in sign according to ε = (±4 + 0.1i)ε0.

Extremely wide-angle absorption for type-I and type-II anisotropic structures is demonstrated clearly in Fig. 2, where we provide results for these two cases for different thicknesses of the metamaterial and values of the Fermi level. The difference between Type I and Type II for a fixed thickness primarily occurs around normal incidence, θ = 0, where the absorption is much higher in the Type I case. It is seen that when the thickness is close to the wavelength [τ/λ = 0.85 in panels (e) and (f)], the presence of graphene provides very high absorption for a remarkably wide range of angles. This should be contrasted to the case without graphene (dashed line), where the absorption peaks for an angle θ < 45° and then rapidly declines. For sufficiently small metamaterial thicknesses, and smaller angles, the metamaterial layer becomes insignificant and the EM wave interacts with mainly the graphene sheet and metal backing. Tuning E is seen to improve the absorption dramatically not only in magnitude, but also in terms of which angles of incidence that become absorbed. Remarkably, the absorption coefficient can be tuned to be nearly perfect, for an unusually broad range of angles. This is distinct from a scenario involving conventional dielectrics where isolated absorption peaks can appear at specific angles of incidence.

Figure 2

Absorption vs. angle of incidence for the Type I case (top row) and Type II case (bottom row) for different thicknesses τ of the metamaterial: τ/λ = 0.01 in (a,b), τ/λ = 0.1 in (c,d), τ/λ = 0.85 in (e,f). The different curves are obtained for different values of E. Type I case: λ/λ = 1.01 with ε = (4 + 0.1i)ε0. Type II case: λ/λ = 0.99 and ε = (−4 + 0.1i)ε0. In both cases, Re and we set the characteristic wavelength of the metamaterial dispersion to λ = 1.6 mm.

To pinpoint the mechanism of the observed perfect absorption, the mode characteristics of the system should be identified. By correlating the longitudinal wavevector component k of the incident EM wave with the permitted EM modes of the graphene-based ENZ structure, the critical coupling responsible for perfect absorption can be identified. Since k/k0 < 1, and ℑm(k) = 0, we search for the allowed fast-wave EM modes (whereby the phase velocities exceed c) that the structure supports51. By solving the transcendental equation R+ = 0 [Eq. (5)], we find the permitted incident wavevector components that yield . The obtained k, are converted to a critical coupling angle via θ = arcsin(k/k0), giving the necessary incident beam angles for fast-wave mode coupling to occur. To illustrate this, in Fig. 3 the perfect absorption curves (solutions to R+ = 0) are shown as a function of E and θ. The wavelength considered λ = 1.01λ corresponds to a negative real part of ε that is also near zero. In (a) ε = (4 + 0.1i)ε0 corresponding to a type I HMM, while in (b) both permittivity components are of the same sign, since now ε = (−4 + 0.1i)ε0. The dispersion curves are shown together with density plots of the absorption , which gives a global view of the high absorption regions spanned by E and θ. Clearly the fast-wave dispersion curves overlap with the regions of perfect absorption. A broader range of incident angles leading to perfect absorption is found in the graphene/HMM structure (a) compared to (b) where both and are of the same sign. It is evident that incident angles corresponding to perfect absorption occur when a fast-wave mode is excited in the structure, due to the coherent superposition of waves with k below the light line. These modes are distinct from conventional guided wave modes that have k/k0 > 1, and consequently the EM fields would be evanescent in the air region near the graphene layer. They also differ from leaky wave modes58 that dissipate when propagating along x since in that case, ℑm{k} ≠ 0. Therefore, by coupling the incident beam to fast wave modes of the structure, the enhanced absorption features provided by the HMM in the ENZ regime61, complimented with the tunable E afforded by graphene, generates a broad range of angles for perfect absorption to arise.

Figure 3

Absorption as a function of E and angle of incidence θ.

The anisotropic metamaterials each have a thicknesses of τ = 0.1λ, and the incident beam has wavelength λ = 1.01λ, corresponding to . In (a) ε = (4 + 0.1i)ε0 corresponding to a type I HMM, and in (b) ε = (−4 + 0.1i)ε0. The fast-wave solutions to the dispersion relation R−(−k) = 0 in Eq. (5) are shown by the bold curves.

The behavior of the energy flow is exhibited in Fig. 4, where the Poynting vector components S and S are shown within the air (z/τ > 1) and metamaterial (z/τ ≤ 1) regions for the structure considered in Fig. 3(a). It is evident that magnitudes of S and S are largest near the interfaces, except for near normal angles of incidence, whereby the energy flows along the metal backing. The flow of S in the metamaterial is opposite to the incident beam, demonstrating negative refraction that arises from the small negative permittivity component ε. The perfect absorption features seen in Fig. 3(a) are also consistent with Fig. (d) for θ = 41° where the energy flow is maximal in the vacuum region at this angle, indicating purely downward energy flow with no upward contributions from surface reflections.

Figure 4

Components of the normalized Poynting vector: (a) S/S0 (along the interface) and (b) S/S0 (normal to the interface) as functions of position z/τ and angle θ. The Fermi energy is set at E = 200 meV. The bottom panels reveal the spatial behavior of S and S for several θ. Vertical lines identify the vacuum/metamaterial interface at z/τ = 1. The system parameters are the same as those used in Fig. 3(a).

Methods

We here provide some details for how to obtain the solution of the Maxwell equations in the presence of graphene. The wave-equations for E and E take the same form in region j = 2, 3:

The dispersion relation then reads: , where we have defined . An electric field propagating in region j is then written: . The amplitudes are related via: Using this relation, we may now write down the solution of the electric field in region j which takes into account both transmitted and reflected waves:

where we defined κ = kε/kε. Since the light is incident from air, we take A3 = 1. The magnetic field is obtained as where:

There are three coefficients to be determined: {A2, B2, B3}. To do this, we need to specify the boundary conditions at the air/anisotropic medium interface (z = τ) and the anisotropic medium/PEC interface (z = 0). Since the tangential component of is always continuous, it follows that , and . The presence of the graphene layer now enters in the boundary condition for the tangential component of by writing the free current as  = σ: . Specifically, we find that the boundary conditions may be written in the form M =  where

where we defined . Thus, the coefficients are solved by  = M−1, which in turn fully determines and throughout the entire structure.

We also provide an explicit expression for the time-averaged Poynting vector in the metamaterial region 2 is:

To determine the fast-wave modes of the system, we solve Maxwell’s equations for the EM fields in the absence of an incident field, thereby effectively treating the structure as a waveguide. This is equivalent to finding the poles of the reflection coefficient for a given sign of kz,3. The corresponding transcendental equation determines the allowed ω and k pair that a source field can excite in the metamaterial system. From the many possible solutions, we consider only fast-waves corresponding to k < k0, and ℑ{k} = 0. The dispersion equation that must be solved is,

where k can be either positive or negative. When the system under consideration is isotropic, the transcendental equation for the modes in the ENZ limit (ε = ε → 0), reduces to,

The k is converted into a coupling angle θ that the incident beam must be directed at in order for fast-wave modes to be excited: θ = arcsin(k/k0). What constitutes a zero numerically in the root finding algorithm that solves Eq. (17) is a tolerance that yields absorption with .

For completeness, we also briefly show results for the spectral dependence of the absorbance, i.e. when moving away from the ENZ regime by varying the wavelength λ of the incident wave. This is shown in Fig. 5 for different thicknesses of the metamaterial layer. As seen, the broad and high absorbance only appears when λ/λ is relatively close to 1. For large deviations from λ/λ = 1, a conventional scenario with isolated absorption peaks arises.

Figure 5

Absorbance vs. angle of incidence for different wavelengths of the incident wave, in order to show the behavior as one moves away from the ENZ regime.

The thickness of the metamaterial is set to (a) τ/λ = 0.01, (b) τ/λ = 0.6, (c) τ/λ = 1.5. We have fixed E = 300 meV, ε = (4 + 0.1i)ε0, and λ = 1.6 mm.

We also examine the effects of loss on the absorption properties of the graphene-based metamaterial. In Fig. 6, we consider the same system that was studied in Fig. 2 except now substantial losses are included whereby the component of the permittivity orthogonal to the interface, ε, has a large imaginary component. Although the wide-angle absorption properties are seen to be reduced in the type-II case, the top panels (b) and (c) reveal that remarkably for the type-I HMM, wide-angle absorption remains intact when considerable losses are present.

Figure 6

Absorption vs. angle of incidence for the Type I case (top row) and Type II case (bottom row).

Here we consider the high loss regime where Im{ε}/ε0 = 1, while all other parameters are the same as those used in Fig. 2.

Additional Information

How to cite this article: Linder, J. and Halterman, K. Graphene-based extremely wide-angle tunable metamaterial absorber. Sci. Rep. 6, 31225; doi: 10.1038/srep31225 (2016).