Abstract
We present a detailed theoretical description of a broadband omnidirectional light concentrator and absorber with cylinder geometry. The proposed optical “trap” captures nearly all the incident light within its geometric cross-section, leading to a broad range of possible applications – from solar energy harvesting to thermal light emitters and optoelectronic components. We have demonstrated that an approximate lamellar black-hole with a moderate number of homogeneous layers, while giving the desired ray-optical performance, can provide absorption efficiencies comparable to those of ideal devices with a smooth gradient in index.
©2010 Optical Society of America
1. Introduction
Recent advances in the development of optical metamaterials, artificial structures with the desired local electromagnetic response determined by artificial subwavelength patterning of the media, have opened a wide range of new concepts and applications, from super-resolution imaging to optoelectronics and photovoltaics [1–8]. In particular, the recently proposed [8] “optical black hole”, an omnidirectional light concentrator and absorber, could significantly enhance the performance of solar energy harvesting systems and lead to novel nonlinear-optical devices. In this work, we present a detailed theoretical description of this device, and develop appropriate numerical methods, which – in addition to the accurate analysis of the optical black hole, - could be applied to a broad range of metamaterial-based optical systems with centrally symmetric gradient index. We note that in the context of this paper, the term “black hole” simply refers to the highly efficient “capture” of the light incident on the device, and does not imply any profound analogy to General Relativity.
The remaining part of the paper is organized as follows. In Section 2, we introduce the initial assumptions, state the geometry of the cylindrical optical system with an axially symmetric graded index, and show basic examples of ray-tracing inside such devices, including the fundamental optical black-hole [8]. Section 3 deals with the ray optical description of the system followed by its full-wave description in Section 4. The analytical and numerical apparatus for approximating the ideal fundamental black-hole with concentric lamellar systems of uniform isotropic layers is then discussed in Section 5. Section 5 also covers derivations of the scattering and absorption cross-sections, which are then compared with their semiclassical approximations in Section 6. Finally, we summarize our work in Section 7.
2. The system
In its original proposal [8], the optical “black hole” (see Fig. 1(a) ) consists of the “payload” core with the radius , supporting the desired functionality (such as e.g. a solar cell for concentrated photovoltaics applications), surrounded by the concentrating shell of radius with an axially symmetric dielectric function, , changing in the radial direction as
whereand is a constant. In order to match the real parts of relative permittivity at the interface the radius of the core is defined by the ratio of the permittivity of the surrounding medium () to that of the core (),Note that we assume that the metamaterial forming the shell, is essentially non-magnetic, so that the refractive index is .3. Ray-optical description
In the semi-classical limit, the radial variation of the refractive index of the shell (where ) corresponds to the effective potential
which for the dielectric permittivity defined by (2) with traps all the incident light trajectories – see Figs. 1(b)-1(e). For these ray trajectories, a straightforward integration of the Hamiltonian equations with the effective potential defined by (4) yields in the polar coordinateswhere m is the (conserved) angular momentum, (with being ), and the constants and are set by the initial conditions of the given trajectory. In particular, for the permittivity defined in the shell (2), Eq. (5) after integration giveswhere is yet another constant defined by .Case corresponds to the gravitational potential (or electrostatic potential) and emulates the motion of a non-accelerated particle in central gravitational field; so that the parabolic trajectories (Fig. 1c) are expected. Case gives the potential of fatal attraction [9]. In fact for all (6) corresponds to the ray trajectories “falling” into the core of the “black hole” – see Fig. 1(b). Figures 1(b)-1(e) show the traces of light rays for p equal to −1, 1, 2, and 3 respectively; the trivial case of the beam propagating in free-space () is omitted.
This concept can be realized in both the spherical and cylindrical geometries, with the latter suitable for optoelectronic devices where the light propagation is generally realized in a waveguide geometry. In the present paper we will focus on this cylindrical realization of the optical black hole, leaving its spherical counterpart to a forthcoming work.
4. The wave optical description
In the present section of the paper, we develop the wave-optical description of the “black hole” light concentrator. We define a given monochromatic electromagnetic vector field by the vectors of electric field and magnetic field , and by its angular frequency ω. We start with the monochromatic Maxwell’s equations , arriving first at [10]:
Getting , from , we then rewrite (7)-(8) as here , , and , where is a smooth spatially-dependent dielectric function, the free-space parameters , c, , and respectively denote the wavenumber, the velocity of light, the permittivity, and the permeability.Here we assume only transverse electric (TE, ) and transverse magnetic (TM, ) waves within a nonmagnetic medium (), so we reduce the above vector equations and use in-plane scalars e and h, which are further on collectively denoted as f.
4.1 Scalar wave equations and separation of variables
Hence, both TE and TM wave Eqs. (9) and (10) reduce to a common scalar wave equation
For any TM-wave () the operator L is given by and does not vanish for spatially dispersive ε, while it vanishes for any TE-wave ().The cylindrical symmetry of the device allows a straightforward separation of radial and angular variables in the cylindrical coordinates, (r, ϕ, z), linked to the Cartesian coordinates (x, y, z) through , . Exploiting the cylindrical symmetry of the dielectric function, , and decomposing the solution of (11) into the cylindrical waves, , we arrive at equation for scalar modes
where the prime denotes the radial derivative, and or for TE or TM modes respectively.4.2 TE and TM solutions for
The solvable cases that we are using here include the trivial instance of for both TE and TM modes. In this case, (12) reduces to the Bessel equation for , and the general solution is a linear combination of two linearly independent Bessel functions [11].
In most of the cases it is more convenient to use the following combination
since it represents the unique decomposition (usually found in scattering problems) into the incident and scattered field if the domain is infinite. On the other hand if the solution is written in the domain, which contains the origin, then the additional boundary condition should be provided to avoid singularity ().For simulations we use the normalized variant of (13), which improves the condition number of the linear system matrix (Section 5.1):
However in some applications, including the calculation of cross sections, a better option is an equivalent combination of the Hankel functions which gives expressions that are easier to handle for mathematical analyses (Section 5.3). Back and forth translations between formulas (13) and (15) are performed by and4.3 TE and TM solutions for the ‘quadratic decay’ of permittivity,
For the cylindrical realization of the black hole, which in its simplest form corresponds to an infinite cylinder with a radial distribution of the dielectric permittivity given byor , all optical modes of the system can also be separated into the TE and TM categories. For the TM wave, while for TE we always have .
For a TM mode, (12) can be therefore reduced as
and for a TE mode asBoth (16) and (17) are the Euler differential equations, which are solved by substitution . Then, since and we arrive at the characteristic equation for with roots , where the general solution readsand the characteristic roots in (18) areor equivalently, , where for TM, or for TE, and in both cases.4.4 Boundary conditions and the mode matching
In previous sections the general analytical solutions were written in the domains, where the dielectric function is defined by either the Bessel equation or the Euler equation. At each interface between such domains, the electromagnetic field should satisfy the standard boundary conditions for the continuity of the tangential components of the magnetic and electric fields. Furthermore, the physical requirement for the behavior at the infinity and at the origin, introduce additional constrains to the solution, as described earlier in Section 4.2
In particular for a nonmagnetic, cylindrically symmetric medium without any free charges the condition of continuity of electromagnetic fields at a given interface reads
TE:
TM:where . Since we assume that both the inner and outer sub-domains are of constant permittivity, with the general solution given by Eq. (13), which should have a physical solution at , then andIf the function f describes the scattered field (i.e. does not include the incident wave), it should obey the Sommerfeld radiation condition, . The asymptotic behavior of the Hankel functions for forces us to choose only one term in (13); hence, for the scattered field we have
4.5 Incident light: the Gaussian beams and plane wave
While the angular momentum formalism introduced above, represents the “natural” choice for a system with rotational symmetry, to complete the description of the scattering problem, we also need to develop the angular momentum expansion of the incident field. Here, we will consider the most common examples of the illumination field: the plane wave and the Gaussian beam.
Starting from the classical generating function for the Bessel functions [12]:
for the incident plane wave we obtain ()While plane wave illumination represents perhaps the most common setup of a scattering problem, it does not convey the case when the target is illuminated by a focused beam. The latter is generally adequately described by the Gaussian beam [13], which yields an accurate description of the beam focused by linear optical elements. However, one needs to keep in mind that the Gaussian beam in its standard form is not an exact solution of Maxwell’s equations in the free space, but instead corresponds to the paraxial approximation [13]. To avoid the approximate nature of this approach and resulting inaccuracy, we consider that the field, while following the Gaussian profile at the waist cross section, is in fact the exact solution of the Helmholtz wave equation with the given direction of propagation [14]:
The resulting solution of (26) in each of half-spaces , can be written as a composition of y-plane waveswhere the amplitudes satisfy the condition . Here, the field is separated into the contributions of the propagating (the first integral in (27)) and the evanescent waves (the second term in (27)). The contribution of evanescent waves, can be estimated asAs a result of (28), for a Gaussian beam with waist propagating in media with neglecting the evanescent waves integral in (27) gives a reasonable approximation, since the relative error of the field magnitude is less than .Adding the arbitrary shift of the focus location by straightforward translation, , the final expression of a focused beam in Cartesian coordinates is
or (29) can be equivalently written in the cylindrical coordinates with and givingwhere .Using the generating function (24) for and we obtain from (30) the final serieswhere .5. Layered Systems
Practical realizations of concentric cylinder optical and optoelectronic devices often involve a design that includes a number of (homogeneous) layers. In particular, external layers may be incorporated as a protective cover, while an internal layer between the absorber and the shell of the optical concentrator provides additional mechanical support or works as a channel for a liquid absorber. In this Section, we develop the theoretical description of light scattering and absorption by such layered systems.
5.1 Cascading the cylindrical layers
The boundary value problem of the scattered cylindrical modes inside a set of concentric cylindrical layers leads to a system of linear equations. Utilizing a direct linear algebra method for obtaining the scattered cylindrical modes inside layers with losses or axially symmetric distribution of the dielectric function could be problematic due to typically large condition numbers of the system matrices. The algorithm that we use includes a normalization step, making it possible to improve poor conditioning. However, if the algorithm were used for getting analytical formulae the normalization would not be necessary.
We consider a concentric cylindrical device having layers, , , with the “outer” () layer filled with a host media, as shown in Fig. 2 .
The general solution of the wave equation in i th layer can be represented as a sum of angular harmonics, where the radial part consists of two partial solutions and with yet undefined coefficients. These coefficients can be obtained from appropriate boundary conditions. All cylindrical modes are linearly independent, thus each of the boundary conditions can be separated into the individual equations for each m th mode, similar to that shown in (14) presented in a more general form,
where for the simplicity of notation (and from now on in this section) we suppress the index m. To satisfy the condition (22) in the core layer we put and .Furthermore, the total field in the “outer” layer can be separated into the “incident” and “scattered” parts. The Sommerfeld radiation condition is imposed on the scattered field (cf. Eq. (23)), while for the incident wave , with the values of the coefficients taken from the expansion of the incident wave (25), (31).
To obtain the rest of equations we use the standard boundary condition (20), (21) for TE polarization at each interface ,
whereFor the TM polarization the only change is that . Inverting the matrix at the left-hand-side of (33) and evaluating its determinant as we obtain
whereand is the Wronskian [12], which is taken for the partial solutions of the layer. In particular, for homogeneous layers we exploit the known value of the Wronskians for the Bessel functions [15],For a layer with the inverse-square radial dependence of the dielectric we have , which yields . This approach is also known as the ‘cascading layers’ method [16], since the solution in each layer (34) is then a composition of transfer matrices applied to the first layer ) [13].The resulting system obtained after stating all 2l boundary conditions has the form:
The system (37) can be solved in two sweeps, similar to the standard Thomas algorithm [17]. In this approach, the auxiliary normalized coefficients , and , are calculated by forward sweep using the recurrence (38)
followed by the backward sweep of the reflection and transmission coefficients (39), while using the initial condition at the external interface:where is known and .The above algorithm includes the normalization of the partial solution in each layer by its value on the inner interface (the field in the core layer is an exception and is normalized by its value on the outer interface , the latter is achieved by substitution in Eq. (32)). The normalization significantly improves the accuracy of the mathematical operations, performed numerically with a limitation of floating-point representations and arithmetic. Finally, the procedure yields the m th mode of the field by using (32); the solution to the entire problem in the i th layer is given as a sum of the angular momentum modes, .
5.2 Scattering and absorption efficiencies
For any optical system, scattering and absorption efficiencies have the physical meaning of the relative size of the equivalent “black body” target which results in the same amount of scattered and absorbed power [18].
First, we recapitulate that the total field (TE) or (TM) outside the device at the external boundary of the shell, is defined as ( is omitted)
Then, we assume that the m th mode of the total field is given by,orso that the total field (40) is decomposed into a plane wave illumination sourceand the scattered fieldwhere coefficients are determined from the boundary value problem. In (41)-(44) and further on we use abbreviations , for , .Equations (40)-(43) are used for the time-averaged Poynting vector , and its normalized radial component , where W is the total power transmitted by the incident wave through –wide strip per its unit length:
and since the absorption efficiency per unit length is taken as a normalized flux of power flowing into the device, , we arrive atAfter applying (41) to (46), we note that includes a vanishing sum of a complex number with its conjugate, , and another term, , so that
which is obtained by substituting with (cf. Eq. (36)).Similarly the scattering efficiency per unit length is calculated by taking the flux of normalized power of the scattered field (44) moving away from the device, , where likewise (45), ; and we obtain
Note that if the solution of the same scattering problem is given for an arbitrary source of illumination as , then due to the linearity of the problem the coefficients relate to those from (41) as . In particular this gives immediate expressions for the absorption and scattering efficiencies, if the total field is given as for which ; then substitution and the use of gives5.3 The ideal electromagnetic black hole
We consider the ideal “fundamental” (p = 2) electromagnetic black hole, see (1), (2), (3). It is a three layered system with a gradient-index shell and absorbing core, for which the formalism of Section 5.1 can be applied to calculate the exact full-wave solution. Fig. 3 shows an example simulation of an ideal black hole with , , and. The device is illuminated with a Gaussian beam (free-space wavelength and full width ), which is focused at fixed , and different (where , and , and 0).
Furthermore the approach is also used to determine the final form of the absorption efficiency. For this we use Eqs. (47), (48) and derive the scattering coefficients. We write the total field outside the device as, where , . Then the exact formula for the scattering coefficients is given as
The coefficients are readily available from (35), then writing the explicit expressions for the matrix elements to (51) we obtain:
whereThe coefficients are the characteristic roots of the Euler equation (cf. Eq. (19)), and factor A is obtained from the boundary conditions,For reshaping we use the identity . Similarly in the expression for X, we use the identity so that coefficients can be written for both polarizations as
5.4 The lamellar electromagnetic black hole
The formalism of a layered system for example can be used to study a non-ideal lamellar “black hole” optical concentrator and absorber, which approximates an ideal device (, , and) with a different number of homogeneous layers, l. The device is illuminated by a plane wave with free-space wavelength ; for l equal to 3, 5, 9, and 17. Scattering and absorption efficiencies are given in (49)-(50), the latter, , yields 72%, 84%, 90%, and 94% respectively. Fig. 4(e) depicts the field map of the ideal black hole with smooth gradient, where 99% absorption efficiency is achieved.
6. Semiclassical description
The primary objective of the present section is the derivation of the semiclassical expansion () of the scattering and absorption cross sections of the “fundamental” (p = 2) electromagnetic black hole, introduced in Section 2.
Intuitively, from the ray dynamics of Section 2, we expect the absorption cross section to be close to the “geometrical” limit, so that . The ray-optical treatment of Section 5 however gives no quantitative measure of how closely this limit is actually approached in a practical device. On the other hand, the exact description of Section 5, while giving a precise quantitative account, provides no compact lucid description of the system behavior. In contrast to this behavior, the semiclassical approach, mathematically corresponding to taking the formal limit of the exact solution, combines the clear physical picture of the actual dynamics in terms of the underlying ray trajectories, with the quantitative description of the interference and diffraction in the system.
The semiclassical approximation for the reflection coefficients and the absorption cross-section can be obtained by taking the limit in (52) and (54). Mathematically, this corresponds to the famous Debye expansion [12] of the Bessel functions
where the “phase” .Physically, the angular momenta summation in absorption efficiency (47) can be separated into two parts – the sum over and the sum over . As the dimensionless angular momentum m has physical meaning of the product of the (fixed) wavenumber and the impact parameter – see Fig. 2 – the sum over corresponds to the rays that actually hit the target and as a result experience strong scattering. On the other hand, the sum for takes into account the rays that miss the object. In the geometrical-optics limit, those rays experience no scattering and thus do not contribute to the cross-section. Beyond ray optics, the waves with , do experience some amount of scattering normally associated with diffraction, however in the limit the resulting contribution is exponentially small. Therefore, in the semiclassical limit the sum over large angular momenta () can be neglected, and we obtain
Furthermore, for the relative change in , when one goes from m to , is reasonably small – leading to the relative error – and thus the summation in Eq. (56) can be approximated by the integral
Note that in this approximation the values of the parameter that are different for the TE and TM polarization (cf. Eq. (53)), should be considered equal, i.e. , then, applying Debye approximation to (52) and using identity we obtain,
Furthermore, for the term can be neglected, in addition, from (54) we haveand since the fraction is very small for the TM polarization and is exactly equal to zero for TE polarization, using the Debye expansion for Bessel function in the limit and low losses in (59) we obtain,Using a similar approximation for in (60) we write,where phase factor is already defined in (55), and .Then, dropping the small term and taking (due to small losses) in (61), we are back to (58) and obtain,
with the magnitude of (62) beingwhere we used, .For the absorption efficiency (57) after substitution (63) we then obtain,
and after variable substitution (64) can be written as a standard integralwhere and are the modified Bessel and Struve functions of the first order, andEquations (65) and (66) finalize the derivation of the semiclassical expansion () of the scattering and absorption cross sections of the “fundamental” () electromagnetic black hole, which is the specific aim of this Section. As expected, for both polarizations in the absence of losses we find zero absorption cross-section. However, even for moderate losses the absorption efficiency is close to one. Thus, as predicted by the semiclassical theory, the device does indeed absorb nearly all the incident light from every direction.
To further illustrate our results we first qualitatively portray the decay of performance as we move away from the semiclassical limit. Fig. 5 shows field patterns of the black hole device ( , , r s = 20 µm, and ) illuminated by a TE-polarized plane wave for a free-space wavelength changing from 1.5 to 6.0 µm, where increasing scattering is observed with increase of wavelength. Then, we perform a quantitative comparison of the absorption efficiency, , versus the ratio .
The efficiencies, shown in Fig. 6a , are separately calculated for the TE- and TM-polarized plane wave using the exact method of Section 5, and the semiclassical result of Section 6, which is valid for both polarizations. Figures 6a also indicates good quality of the semiclassical approximation of the absorption cross-section even far beyond the semiclassical limit. Figure 6b depicts a more detailed plot for changing from 0.075 to 0.15.
7. Summary
In summary, we have developed a theoretical description of wave propagation in cylindrically-symmetric gradient-index systems, and applied this approach to light “trapping” in the recently proposed electromagnetic black hole concentrator and absorber.
In particular we have demonstrated that an approximate lamellar black-hole with a relatively limited number of homogeneous layers, while giving the desired ray-optical performance, can provide absorption efficiencies comparable with those of ideal devices with smooth gradients in index. It is forecasted that the use of non-uniform layers could further improve the performance of the lamellar device with a given number of homogeneous layers. The numerical validation of the mathematical apparatus for designing the black-hole devices opens the path for the development of test beds that could further advance the application of this light-trapping technique into a viable commercial option.
Acknowledgments
This work was partially supported by ARO MURI56154-PH-MUR, 50342-PH-MUR, and also through IIP SB RAS No 113, and RFFS No 09-01-00352-а.
References and links
1. V. M. Shalaev, W. S. Cai, U. K. Chettiar, H. K. Yuan, A. K. Sarychev, V. P. Drachev, and A. V. Kildishev, “Negative index of refraction in optical metamaterials,” Opt. Lett. 30(24), 3356–3358 (2005). [CrossRef]
2. J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 85(18), 3966–3969 (2000). [CrossRef] [PubMed]
3. Z. Jacob, L. V. Alekseyev, and E. Narimanov, “Optical Hyperlens: Far-field imaging beyond the diffraction limit,” Opt. Express 14(18), 8247–8256 (2006). [CrossRef] [PubMed]
4. A. Salandrino and N. Engheta, “Far-Field Subdiffraction Optical Microscopy Using Metamaterial Crystals: Theory and Simulations,” Phys. Rev. B 74(7), 075103 (2006). [CrossRef]
5. N. M. Litchinitser and V. M. Shalaev, “Metamaterials: transforming theory into reality,” J. Opt. Soc. Am. B 26(12), B161–B169 (2009). [CrossRef]
6. N. I. Landy, S. Sajuyigbe, J. J. Mock, D. R. Smith, and W. J. Padilla, “Perfect metamaterial absorber,” Phys. Rev. Lett. 100(20), 207402 (2008). [CrossRef] [PubMed]
7. T. V. Teperik, F. J. García de Abajo, A. G. Borisov, M. Abdelsalam, P. N. Bartlett, Y. Sugawara, and J. J. Baumberg, “Omnidirectional absorption in nanostructured metal surfaces,” Nat. Photonics 2(5), 299–301 (2008). [CrossRef]
8. E. E. Narimanov and A. V. Kildishev, “Optical black hole: broadband omnidirectional light absorber,” Appl. Phys. Lett. 95(4), 041106 (2009). [CrossRef]
9. L. D. Landau, and E. M. Lifshitz, Mechanics, Pergamon Press, Oxford, (1976).
10. R. K. Luneburg, Mathematical Theory of Optics, University of California Press, Berkeley, 1964, p. 12.
11. E. J. Post, Formal Structure of Electromagnetics: General Covariance and Electromagnetics, 1962, p. 152.
12. S. Gradshteyn, and I. M. Ryzhik, Tables of integrals, series and products, Academic Press, New York, CD-ROM Edition, 1994, Eq. (8).511.4.
13. M. Born, and E. Wolf, Principles of Optics, Electromagnetic Theory of Propagation, Interference and Diffraction of Light, 7th ed. (Cambridge University Press, New York, 1999).
14. G. P. Agrawal and D. N. Pattanayak, “Gaussian beam propagation beyond the paraxial approximation,” J. Opt. Soc. Am. 69(4), 575 (1979). [CrossRef]
15. M. Abramowitz, and I. A. Stegun, eds., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, New York, Dover, (1972).
16. V. Kildishev, U. K. Chettiar, Z. Jacob, V. M. Shalaev, and E. E. Narimanov, “Materializing a binary hyperlens design,” Appl. Phys. Lett. 94(7), 071102 (2009). [CrossRef]
17. S. D. Conte, and C. W. Boor, “Elementary Numerical Analysis: An Algorithmic Approach, 3rd edition,” 1980.
18. H. Hulst, Light scattering by small particles, Dover, 1981, p. 309.