A simplified transfer function approach to beam propagation in anisotropic metamaterials
Introduction
There has been enormous interest in the development of periodic or almost periodic electromagnetic (EM) and optical structures and heterogeneous structures that can exhibit interesting properties, such as prescribed band gaps and tuning [1] sub-wavelength imaging [2], efficient harmonic generation [2], invisibility cloaking [2], [3], etc. The artificial structures that can achieve these properties are referred to as metamaterials. These metamaterials interact with EM radiation to achieve exotic material properties, such as negative permittivity, negative permeability, and negative refractive index, leading to effects such as negative refraction, etc., many of which has been experimentally demonstrated [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15]. While most of the theory is based on plane waves, the behavior of beams is more important from a practical standpoint. As is well-known, the beam propagation method (BPM) is employed to study the properties of various guided wave optics structures, such as optical waveguides with bends and tapers, gratings and spectrum analyzers, etc. [9], [10]. The propagation of beams in metamaterials can be modeled using the fast Fourier transform (FFT-BPM), or finite element and finite difference methods [9], [10], [11], [12], [13], [14], [15].
While analyzing beam propagation in anisotropic media, theanisotropy can be expressed as matrices for the electric permittivity, and magnetic permeability, [2], [13]. For any FFT based BPM, the transfer function of propagation is required to relate the angular plane wave spectrum at one plane to the spectrum at another plane during propagation [9], [10], [11], [12], [13], [14], [15]. In anisotropic media, the transfer function matrix is required.
Hyperbolic metamaterials (HMMs) are a form of anisotropic material where the dielectric tensor elements have opposite signs [16], [17]. There are two types of HMMs which can be distinguished by the signs of the principal elements of the diagonal permittivity matrix: or [15]. The dispersion occurs in these metamaterials when the principal components of diagonalized matrices change the sign [4], [15]. It has been argued that negative refraction can be achieved through the hyperbolic dispersion of these materials [16], [17], [18], [19], [20]. HMMs have many potential applications in EM and optics. HMMs are used in several optical applications, such as waveguiding, imaging, sensing, quantum, and thermal engineering [16], [17]. In this work, beam propagation through HMMs is analyzed.
Incidentally, anisotropic metamaterials can be built as a multilayer structure comprising alternating layers of dielectric and metal, and modeled as an anisotropic bulk medium (BM) with an effective dielectric permittivity tensor based on effective medium theory if the thickness of the metal and dielectric are much smaller than the wavelength [21], [22], [23], [24], [25]. Reflection and transmission from slabs of anisotropic metamaterials comprised of metallo-dielectric structures have been analyzed as a function of the angle of incidence, number of layers, and layer thicknesses, and it has been shown that effective medium theory can give results that are sufficiently close to, and computationally much faster than, rigorous computations using the transfer matrix method and finite element methods such as COMSOL [25], [26], [27], [28], [29], [30].
As stated above, FFT-BPM is a powerful tool for the computation of EM wave propagation in anisotropic media [26]. The common analysis technique is to decompose the beam into an infinite number of plane waves traveling in different directions, or plane wave spectra [31], [32], [33], [34], [35], using the Fourier transform. Propagation of the plane wave spectra is achieved by multiplying the initial angular spectrum with the transfer function for propagation, which can be done analytically or computationally [33]. The transfer function matrix for propagation in an anisotropic medium can be derived by starting from the Berreman matrix [36], [37].
Also, in related earlier work, it has been shown by Banerjee et al. [38] that propagation of radially polarized Gaussian beams can be readily analyzed rather simply by knowing the q-transformation of the underlying scalar Gaussian beam. The exact profiles of the longitudinal and transverse components of initially radially polarized lowest-order Laguerre–Gaussian beams have been derived using the transfer function approach which leads to the laws of q-transformation [39].
The objective of this work is to extend existing plane wave theory to beam propagation in anisotropic metamaterials using the notion of transfer functions in order to help understand effects such as self-focusing, negative refraction, etc. of beams in hyperbolic metamaterials from an analytical and physical standpoint. The transfer function matrix for beam propagation can be readily found from the Berreman matrix. This is summarized in Section 2, leading to the derivation of simplified transfer functions, which are subsequently used to analyze propagation of one transverse-dimensional (1TD) TE and TM polarized beams in semi-infinite HMMs. The paraxial version of the transfer function is used to analytically derive the propagation of Gaussian beams using the q-parameter approach. In Section 3, the transmission through a slab of a HMM is analyzed by extending the plane wave approach to an angular spectrum approach, by simply relating the angle of propagation to the corresponding spatial frequency to derive the transmission transfer function. The transfer function approach can be readily extended to determine the propagation of the longitudinal or z-component of the optical field. Section 4 shows illustrative examples. For instance, it is shown that in such semi-infinite HMMs, 1TD TM polarized Gaussian beams can undergo linear self-focusing, while 1TD TE polarized beams spread due to diffraction, and the results from paraxial theory match well with rigorous non-paraxial calculations. Simulations using the transfer functions include examples of realistic HMMs with complex permittivities derived from exact (complex) values of the permittivities of the metal and the dielectric; and the results are compared with those obtained from setting the imaginary part of the effective permittivity to zero. Also, it is shown that a slab of a HMM can cause negative refraction only of a 1TD TM polarized beam, evidenced from the spatial shift of the beam upon exit from the material.
Liu et al. [11] have a full vectorial description of the focusing of EM beams in media with hyperboloid and ellipsoid dispersion surfaces, and show that a paraxial beam undergoes “focusing in both x and y axes for the slab with an ellipsoid dispersion relation”, and undergoes “focusing or dispersion in different axes for the slab exhibiting a dispersion surface of a single-sheeted hyperboloid”. A nonparaxial correction to the paraxial solution is then derived. In our case, we start from the unidirectional nonparaxial solution and make the paraxial approximation which leads to the q-transformation of a Gaussian beam in HMMs. This formulation enables us to analytically predict the location of the self-focusing distance, which is also verified by numerical simulations. Comparison of non-paraxial and paraxial simulations based on the respective transfer functions can clearly show when the paraxial approximation is valid. Kivijärvi et al. [12] have a detailed theory of vectorial beam propagation in anisotropic metamaterials showing the intensity profiles of a focused beam through metamaterials composed of, for instance, a period array of silver disk and dimer metamolecules. Although based on angular spectrum decomposition, the rigorous computations are performed using COMSOL; furthermore, arbitrary radii of curvatures and oblique incidence are not examined. While complete in treatment, the physics of self-focusing and spatial shift due to propagation in HMMs for different polarizations are not clearly apparent. Zhou et al. have examined extraordinary reflection and refraction in natural hyperbolic materials for wave propagation [40]. As discussed above, our work extends explicitly to beam propagation, and the anticipated beam shift is verified using the spatial frequency dependence of the transmission function, which can be regarded as the transmission transfer function for the hyperbolic system. Incidentally, beam shift in HMMs has also been shown by Argyropoulos et al. [2], but using finite element based simulations from the outset. 2-D Gaussian beams have also been studied in the past for uniaxial crystal and uniaxial inhomogeneous media for transmission properties [41], [42]. Reflection and transmission properties have been derived using the Berreman matrix [42], [43] for a homogeneous slab, but these studies have not discussed Gaussian beam propagation in a hyperbolic medium. We feel the uniqueness of our work lies in explaining the details about beam propagation and negative refraction using BMM along with the effective medium theory along with a transfer function approach based on angular spectral decomposition of beams.
Section snippets
Mathematical derivation of transfer function matrix
To begin, a uniaxial material with only diagonal nonzero elements for , magnetic isotropy, and two transverse dimensions , with being the nominal direction of propagation is assumed (see Fig. 1). With the substitution , Maxwell’s equations can be recast into the 4 × 4 matrix equation involving the transverse components of the electric and magnetic fields, with (Eq. (2) ) given in Box I,
and where the tildes signify the spectra of the
Focusing in HMMs
In this subsection, typical simulation results for Gaussian beam propagation in a semi-infinite slab of the HMM are presented for TM and TE polarizations. The initial width and radius of curvature of the Gaussian beam in the metamaterial can be chosen independently. It is clear from the paraxial transfer functions discussed in Section 2.2 that focusing is expected for TM polarization and not for TE polarization. Different beam waists are chosen along with different (diverging) radii of
Conclusion
The transfer function for propagation TM and TE polarized beams and their corresponding z-polarizations have been derived starting from the Berreman matrix for a uniaxial hyperbolic metamaterial, along with the transmission transfer function. The transfer function concept has been used to analyze self-focusing and negative refraction in hyperbolic metamaterials. The reported work should be useful in providing an analytical tool to cross-check rigorous simulations using standard software
References (46)
- et al.
Berreman approach to optical propagation through anisotropic metamaterials: application to metallo-dielectric stacks
Opt. Commun.
(2018) From photonic band gaps to refractive index engineering
Opt. Mat.
(1999)- et al.
Negative refraction, gain and nonlinear effects in hyperbolic metamaterials
Opt. Express
(2013) - et al.
Review of left-handed metamaterials: evolution from theoretical and numerical studies to potential applications
J. Electromag. Waves
(2006) - et al.
Metamaterial electromagnetic cloak at microwave frequencies
Science
(2006) - et al.
Optical and acoustic metamaterials: superlens, negative refractive index and invisibility cloak
J. Opt.
(2017) - et al.
Experimental verification of a negative index of refraction
Science
(2001) - et al.
Experimental demonstration of electromagnetic tunneling through an epsilon-near-zero metamaterial at microwave frequencies
Phys. Rev. Lett.
(2001) The electrodynamics of substances with simultaneously negative values of and
Sov. Phys. Usp.
(1968)- et al.
Beam propagation method in anisotropic media
Appl. Opt.
(1982)
Beam propagation in uniaxial anisotropic media
J. Opt. Soc. Amer.
Focusing of vectorial fields by a slab of indefinite media
J. Opt. A: Pure Appl. Opt.
Interaction of metamaterials with optical beams
New J. Phys.
Optical negative refraction in bulk metamaterials of nanowires
Science
Midinfrared semiconductor optical metamaterials
J. Appl. Phys.
Sub-wavelength imaging at optical frequencies using a transmission device formed by a periodic layered metal–dielectric structure operating in the canalization regime
Phys. Rev. B
Hyperbolic metamaterials
Nat. Photonics
Hyperbolic metamaterials: new physics behind a classical problem
Opt. Express
Composite medium with simultaneously negative permeability and permittivity
Phys. Rev. Lett.
Photonic crystals with silver nanowires as a near-infrared superlens
Appl. Phys. Lett.
Negative refraction in semiconductor metamaterial
Nature Mater.
Far-field optical hyperlens magnifying sub-diffraction-limited objects
Science
Directed subwavelength imaging using a layered metal–dielectric system
Phys. Rev. B
Cited by (5)
Optical propagation through metamaterial structures with multilayered metallo-dielectrics: Hyperbolic dispersion and transmission filters
2021, Thin Film Nanophotonics: Conclusions from the Third International Workshop on Thin Films for Electronics, Electro-Optics, Energy and Sensors (TFE3S)An anisotropic transfer matrix approach to profiled optical field propagation through hyperbolic metamaterials
2022, Proceedings of SPIE - The International Society for Optical EngineeringVortex beam manipulation through a tunable plasma-ferrite metamaterial
2021, Scientific ReportsPhotonic Crystal for Beam Tuning Application
2020, 2020 IEEE Research and Applications of Photonics in Defense Conference, RAPID 2020 - Proceedings