Abstract
Based on both analytical dipole model analyses and numerical simulations, we propose a concept of coherent perfect nanoabsorbers (CPNAs) for divergent beams. This concept makes use of the properties of a slab with negative refraction and small losses. The proposed CPNA device would allow focusing radiation in nanoscale regions, and hence could be applied in optical nanodevices for such diverse purposes as reading the results of quantum computation which is based on single photon qubits.
©2012 Optical Society of America
1. Introduction
Recently, a fascinating concept of coherent perfect absorbers (CPAs) has been suggested [1, 2]. In some sense this concept is an inverse to the concept of lasing. This conception was demonstrated experimentally by the example of a symmetric Fabry-Perot cavity with small losses which was symmetrically irradiated by 2 coherent plane waves (see Fig. 1 ). In [1, 2], it was shown that under certain conditions, all the energy of incoming waves is fully absorbed by the dielectric (silicon) slab. As a result, the radiation from the system disappears for some parameters. This situation indeed is reciprocal to the case of lasing, where for high enough pumping the coherent radiation appears. More complicated examples of CPA with plane waves were considered in [3–5].
Despite generic nature of this concept, its application to more complicated geometries is a nontrivial task. In particular, since now the control of radiation of atoms and molecules with nanoparticles and metamaterials becomes very important, it is desirable to have an effective nanoabsorber of electromagnetic fields from a single atom or molecule. First of all, one may have a complicated spatial structure of light field. For example, the emission pattern of an atom would correspond to spherically divergent beams rather than plane waves. Moreover, many applications nowadays are related with making use of metamaterials for the control of light at the nanoscale [6, 7].
In the present work, we propose a concept of coherent perfect nanoabsorbers (CPNAs) using a slab made of double negative (DNG) metamaterial, that is, the metamaterial with negative refractive index [6]. Nowadays, slabs with negative refraction are widely used for many applications such as perfect lensing [8] and cloaking of small objects [9]. Our concept is based on recently discovered focusing properties of a system of sources and sinks near a slab with negative refraction [10, 11], where it was shown that if two sources and one sink were put at the ray intersection points, all the energy from the sources would go to the point sink. In this paper we propose to use a real nanoparticle instead of the hypothetic sink used in [10, 11]. The operation scheme of our CPNA is shown schematically in Fig. 2 . The device shown in Fig. 2 would allow absorbing all the energy from the 2 sources emitted in the direction towards the negative refraction slab by a nanoparticle placed inside it.
It should be emphasized that within the geometry of Fig. 2, only half of the radiated energy is going towards the slab and then absorbed. However, it is a trivial task to slightly modify this system by adding perfect magnetic mirrors just behind the x-polarized dipoles. This would result in doubling of the radiated energy going to the slab and then absorbed. For z-oriented dipoles, one should use perfect electric mirrors to absorb all the radiated energy. For simplicity, we will consider here only the systems without mirrors.
2. Analytical formulation
To describe the idea of our proposal, let us start from a 2D DNG slab without losses with ε2 = −1, µ2 = −1 placed in vacuum with ε1 = 1, µ1 = 1. Let us put into this system 3 2D dipoles which are described by line currents
where p0 is the dipole moment of the wires per unit length, and everywhere the time dependence is assumed. It can be easily verified that the magnetic fields from these currents will have the following formwhere k0 = ω/c and is the Hankel function. The electric fields can then be found with the help of Amper’s lawFrom Eq. (2) one can see that there are 3 outgoing 2D waves for the assumed time dependence. However, taking into account the fact that in the slab we have negative refraction, we immediately see that 2 of these waves (right and left) correspond to the energy sources while the third one, which is placed into the negative refraction slab, corresponds to an energy sink. The distribution of the energy flows for solution Eq. (2) is shown in Fig. 3 , and we can see from this figure that all the energy from the 2 sources radiated in the direction of the slab is absorbed by the sink dipole.It is easy to show this directly by calculating the work of the field over the sink dipole, that is the absorbed power,
Because there is no reflection, the expressions for the fields are the same in all regions [see Eq. (2)]. It is easy to see that the absorbed energy Eq. (4) is equal to the total energy radiated by one source or half of the energy radiated by 2 sourcesThus, the 3 dipole system Eq. (2) is a perfect realization of our CPNA if one approximates the nanoparticle by a dipole.However, real systems cannot be without losses and it is important to demonstrate the idea of coherent perfect nanoabsorber in more realistic situation with losses. To do this, let us consider a DNG slab with losses, that is ε2 = -1 + iδ, µ2 = -1 + iδ, and again place into this system 3 dipoles. However, now we do not know what amplitude of the dipole moment of thenanoabsorber would result in the CPNA, because now we should take into account its self-consistent interaction with the slab and other sources. Therefore, now we will consider a system of currents
where the amplitude of the dipole moment, pabs, inside the slab is an unknown quantity to be determined through the condition of perfect absorption. To prove the existence of CPNA for this case, let us again calculate the energy absorbed by the dipole inside the slabNow, however, the fields are more complicated. In particular, the field near this absorbing dipole is composed also by reflections from the walls and from the dipoles placed outside the slab, that iswhere are, respectively, the fields from first and second sources at the center of the absorbing nanoparticle, while are, respectively, the fields from in unbound space with ε2 and μ2 and corresponding reflected part at the center of the absorbing nanoparticle. If one presents electric fields through Green functions of the dipole source p0 placed at r0, that is through relationone obtainsinstead of Eq. (7). In Eq. (10), we have taken into account the fact that the fields on the z-axis have only nonzero x-component. Equation (10) is a bilinear form of the and it will definitely have a maximum which would provide the perfect absorption. To find this maximum, we should find all the Green functions entering Eq. (10).First of all, let us find the self-action part of the field,. It can be easily found because this part appears as a solution in uniform space with parameters of the slab, that is, in this case
where . Thus, at the absorber position, we haveTo calculate Eq. (12), one should use an integral presentation for Hankel function [12]
In Eq. (13) and other places, we use notion (Im kz > 0) for longitudinal wavevectors and q is its transversal part. As a result, the expression for free Green function at the origin will have the form
where integration here is in finite limits, from -1/R to 1/R, and R is the radius of the nanoparticle. The necessity to limit the integration domain is related with the fact that the absorbed or radiated power of the point source is infinite even in media with arbitrary low losses.To find other Green functions in Eq. (10), let us start from the case when the DNG slab with arbitrary dielectric constant is illuminated with line source #1 with current density
where p0 is the dipole moment of wire per unit length. For such a source, the magnetic field in free space with ε = 1, µ = 1 and the wavenumber k = ω/c has only one component which can be presented in the form of Eq. (2). Now using integral presentation for Hankel function Eq. (13), the #1 source field in unbounded space can be presented in the formUsing this presentation, one can write general expressions for the fields in all regions of space (see Fig. 2)
The coefficients A(q), B(q), C(q), D(q) can be found from the continuity of tangential components of magnetic and electric () fields. Now finding from Eq. (17) the x-component of the electric field and setting z = l and x = 0, we will obtain . Analogously, one can find and . In the fully symmetric case l = d/2 these expressions can be presented in the form
Now all the Green functions for any parameters can be calculated numerically since we know their analytical expressions. In Fig. 4 , we display the absorbed energy as a function of the relative amplitude of the dipole moment of the nanoparticle (the solid curve).
We can see from Fig. 4 that indeed for ξ» −0.6, there is a maximum of absorption which corresponds to CPNA. Here it is important to note that the energy absorbed by the nanoparticle is 5 times greater than the power of the single source in free space. Of course there is no violation of the energy conservation law here because, due to the Purcell effect [13,14], the radiated power can be substantially enhanced by the environment.
To calculate the effective emitted energy, one should again calculate the work done by the source over the field, that is
Now, the field near the dipole of the first source is formed also by reflections from the walls and from the other dipoleswhere are fields from first and second sources at the first source position r1, while is the field from the nanoparticle dipole at the first source position. If one again present the electric fields through Green functions [Eq. (9)], one obtainsinstead of Eq. (19). Green functions in Eq. (21) can be calculated in full analogy with the calculation of the Green functions of Eq. (18). The dashed line in Fig. 4 shows the dependence of normalized radiated energy [Eq. (21)] on the relative amplitude of the nanoparticle dipole moment. One can see that indeed the emitted energy grows linearly as the dipole moment of the absorber increases.Another important feature of Fig. 4 is that for relative dipole amplitude and for , the absorption become negative such that it changes to radiation, that is, some kind of lasing occurs if we provide such dipole amplitudes. However, as we will see below, it is not an easy task because these conditions need negative absorption in nanoparticles, that is gain media of the nanoparticle.
Since we have approximated the nanoparticle by a point dipole, the results of our theoretical calculations presented in this section could be only qualitative. To describe the CPNA in a quantitative manner, one should use full scale computer simulations. Nevertheless, the analytic dipole model analyses reported here give us a clear understanding of the physics underlying CPNA.
3. Numerical simulations
To illustrate the principle of operation of such type devices, we have also carried out numerical simulations of the system shown in Fig. 2, using the finite elements method (FEM) within the COMSOL Multiphysics® software. We use a nanocylinder with parameters as the absorbing element and 2 dipole nanowires with current densities
as sources. That it, we are considering a 2D geometry with p0 as the dipole moment of the wire per unit length and only nonzero Hy component of the magnetic field.Typical distribution of the energy flows is displayed in Fig. 5 . Figure 5 shows that for the chosen parameters of the cylinder, indeed, all the energy from the sources incident on the slab is totally absorbed by the cylinder. From Fig. 5, one can see also that the energy flows in this case are very close to the results obtained for the 2 point sources and one point sink in the slab without losses (see Fig. 3).
In Fig. 6 , the distribution of the squared electric fields in CPNA is shown. One can see from Fig. 6 that indeed the electric fields have maxima at the source positions and also at theabsorbing nanoparticle position. This fact strongly supports the idea of symmetry between sources and sinks which follows from the idealized solution Eq. (2). It can be also seen from Fig. 6 that the electric fields at the interface have substantially smaller amplitudes which means that the surface plasmon excitation is indeed suppressed despite very small losses in the slab metamaterial.
To find the optimal conditions for operation of CPNA, we have displayed the energy absorbed by the nanoparticle as a function of both the real and imaginary parts of its permittivity as a 2D contour plot in Fig. 7 .
From this figure one can see that at εc = 1.3 + 0.4i there is a rather broad maximum of the absorbed power and that this maximum is about 10 times more than the power radiated by the single source in free space. To understand how this happens, we have plotted in Fig. 8 both the emitted and absorbed energy as a function of the imaginary part of the permittivity of the absorbing nanoparticle.
From Fig. 8(a), one can see that the absorbed energy is even greater than the energy radiated by only one source in the region . Of course, there is no violation of energy conservation here. The reason is that due to the presence of the absorber, the radiation pattern from the sources becomes asymmetric and more than ½ of its energy is going to the slab. Note that one can increase this contribution to 1 by introducing corresponding mirrors. Another interesting feature of Fig. 8 is the dependence (growth) of energy radiated by the source on the absorber parameters. Again there is no magic. It is the well known Purcell effect [13,14] and the nanoparticle increases the radiation power of the source. Note that the feedback effects of the fields on the source current amplitude are not taken into account here. Figure 8(b) shows the absorbed and radiated energies as a function of the dipole moment of the nanoparticle. It is gratifying to see that the numerical results in Fig. 8(b) are in good agreement with that of the analytical model analyses described in the preceding section (Fig. 4). In contrast to Fig. 8(a), the dipole moment of the nanoparticle in principle can take any value. Interestingly, for and also for , the absorption becomes negative, and this corresponds to the negative values of the imaginary part of the cylinder permittivity, that is, to a gain media, as can be seen by a comparison of Fig. 8(b) with Fig. 8(a). This, excitingly, indicates that our CPNA can possibly be used to provide nanolasing as well. In short, our numerical simulations show indeed that the proposed system can be used as a CPNA for the 2 divergent sources and also as a nanolaser.
4. Discussion and conclusions
Both analytical dipole model analyses and numerical simulations have shown that an absorbing nanoparticle properly placed inside a DNG metamaterial slab can absorb up to 100% incoming energy from 2 identical divergent sources placed symmetrically outside of a slab. This allows us to put forward the concept of CPNA for diverging beams. Furthermore, in contrary to CPA with usual dielectric cavity [1–5], the CPNA device proposed is robust to small perturbations in the properties of the absorbing nanoparticle and the metamaterial as well as the position of the sources.
The generalization of this approach to 3D systems is straightforward because it has been shown that in a 3D system of the sources and sinks, analogous concentration and absorption of energy are also possible [10,11]. Besides, if we put our system in a properly tuned cavity, we can provide 100% absorption of energy emitted by 2 sources.
Importantly, our approach can also be applied to create single side perfect nanoabsorbers. Indeed due to the symmetry of our solution, one can simply cut the system along the symmetry plane and put a perfect electric conductor (PEC) or perfect magnetic conductor (PMC) mirror here (see Fig. 9 ). Furthermore, our concept can be generalized to more complicated geometries such as wedge and sphere.
Let us stress again that our system allows focusing radiation in nanoscale regions and can be applied in optical nanodevices of different purposes. For example, one may use double side or single side CPNA to arrange readout of the results of quantum computation which are based on single photon qubits.
Acknowledgments
V. Klimov would like to express his gratitude to the Russian Foundation for Basic Research (grants ## 11-02-91065, 11-02-92002, 11-02-01272, 12-02-90014) and the Presidium of the Russian Academy of Sciences for financial supports. S. Sun and G. Y. Guo would like to thank the National Science Council and National Center for Theoretical Sciences of Taiwan as well as Center for Quantum Science and Engineering, National Taiwan University (CQSE-10R1004021) for financial supports. The authors thank Konstantin Simovski and Olivier Martin for fruitful discussions and valuable comments on this work.
References and links
1. Y. D. Chong, L. Ge, H. Cao, and A. D. Stone, “Coherent perfect absorbers: time-reversed lasers,” Phys. Rev. Lett. 105(5), 053901 (2010). [CrossRef] [PubMed]
2. W. Wan, Y. D. Chong, L. Ge, H. Noh, A. D. Stone, and H. Cao, “Time-reversed lasing and interferometric control of absorption,” Science 331(6019), 889–892 (2011). [CrossRef] [PubMed]
3. S. Dutta-Gupta, O. J. F. Martin, S. D. Gupta, and G. S. Agarwal, “Controllable coherent perfect absorption in a composite film,” Opt. Express 20(2), 1330–1336 (2012). [CrossRef] [PubMed]
4. S. Longhi, “PT-symmetric laser absorber,” Phys. Rev. A 82(3), 031801 (2010). [CrossRef]
5. S. Longhi, “Coherent perfect absorption in a homogeneously broadened two-level medium,” Phys. Rev. A 83(5), 055804 (2011). [CrossRef]
6. V. Shalaev and W. Cai, Optical Metamaterials: Fundamentals and Applications (Springer, 2010).
7. Z. Jacob, J.-Y. Kim, G. V. Naik, A. Boltasseva, E. E. Narimanov, and V. M. Shalaev, “Engineering photonic density of states using metamaterials,” Appl. Phys. B 100(1), 215–218 (2010). [CrossRef]
8. J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 85(18), 3966–3969 (2000). [CrossRef] [PubMed]
9. J.-W. Dong, H.-H. Zheng, Y. Lai, H.-Z. Wang, and C. T. Chan, “Metamaterial slab as a lens, a cloak or an intermediate,” Phys. Rev. B 83(11), 115124 (2011). [CrossRef]
10. V. Klimov, J. Baudon, and M. Ducloy, “Comparative focusing of Maxwell and Dirac fields by negative refraction half-space,” Europhys. Lett. 94(2), 20006 (2011). [CrossRef]
11. V. Klimov, “Novel approach to a perfect lens,” JETP Lett. 89(5), 229–232 (2009). [CrossRef]
12. W. C. Chew, Waves and Fields in Inhomogenous Media (IEEE Press Series on Electromagnetic Wave Theory, 1995).
13. E. M. Purcell, “Spontaneous emission probabilities at radio frequencies,” Phys. Rev. 69, 681–681 (1946).
14. V. V. Klimov, M. Ducloy, and V. S. Letokhov, “Spontaneous emission of an atom in the presence of nanobodies,” Quantum Electron. 31(7), 569–586 (2001). [CrossRef]