Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Phase diagram for passive electromagnetic scatterers

Open Access Open Access

Abstract

With the conservation of power, a phase diagram defined by amplitude square and phase of scattering coefficients for each spherical harmonic channel is introduced as a universal map for any passive electromagnetic scatterers. Physically allowable solutions for scattering coefficients in this diagram clearly show power competitions among scattering and absorption. It also illustrates a variety of exotic scattering or absorption phenomena, from resonant scattering, invisible cloaking, to coherent perfect absorber. With electrically small core-shell scatterers as an example, we demonstrate a systematic method to design field-controllable structures based on the allowed trajectories in this diagram. The proposed phase diagram and inverse design can provide tools to design functional electromagnetic devices.

© 2016 Optical Society of America

1. Introduction

With state-of-the-art nano-optical technologies, there are renewed interests on having an efficient way to manipulate and design nanostructures with unusual electromagnetic properties [1, 2]. In particular, resonant scattering [3–6], coherent perfect absorption [7–10], invisible cloaking [11–16], subwavelength superscattering [17–19], and minimum-scattering superabsorbers [20, 21], are revealed in a single, isotropic, and multi-layered scatterer. These meta-structures could be expected to serve as functional nano-devices with promising applications in light harvesting [22–24], heat generation by metal nanoparticles [25, 26], optical nanocircuits and nonlinear optical processes [27, 28].

To have exotic electromagnetic properties at subwavelength scale, a variety of specific conditions are asked to be satisfied. Undoubtedly, a better understanding in the scattering coefficients could provide an access to design nanostructures. In general, we need to consider information about scattering limitation, power assignment, scattered radiation pattern, and robustness on the corresponding extrinsic field response of real scatterers. For working frequency of interests, for example, many metals are associated with a strong dispersion in the visible spectra, that introduces real loss effects and suppresses desired functions [29, 30]. As possible mismatching in physical parameters may occur, it is natural to seek optimized invisible cloaks or performance boundary in a cloaked sensor with the consideration of intrinsic loss in reality [30].

In this paper, we study the general relation between amplitude and phase in the scattering coefficients for any passive electromagnetic scatterers. A phase diagram is introduced by imposing the power conservation on absorption cross section for each partial wave channel, which acts as a universal map to design passive scatterers. Not only all physically allowed regions can be defined to satisfy the intrinsic power conservation, but also all exotic electromagnetic properties in the literature can be illustrated in this phase diagram. Moreover, we take electrically small core-shell scatterers as an example to illustrate a systematic way in designing the composition of subwavelength-structures with required scattering and absorption properties.

2. Phase diagram for a passive scatterer

We consider a linearly polarized plane wave with time evolution e−iωt at the angular frequency ω, which is illuminating on a single spherical object. The object could be made of multiple layers of uniform and isotropic media with complex permittivity and permeability, denoted as ε = ε′+ ″ and μ = μ′+ ″, respectively. Here, ε″ and μ″ are both assumed to be positive real numbers for a passive medium. Without loss of generality, the surrounding environment is taken as non-absorptive, non-magnetic, and free of external sources or currents, i.e., ε0 = μ0 = 1. In the following, we express the electric field E and magnetic field H in the environment by two auxiliary vector potentials, i.e., the transverse magnetic (TM) and transverse electric (TE) modes, which are respectively generated by two scalar spherical wave equations. Each scalar functions can be built by an infinite series with unknown coefficients determined through boundary conditions. By following the conventional notations, let the scattering coefficients be CnTM and CnTE for the transverse magnetic (TM) and transverse electric (TE) modes in each spherical harmonic channel labeled by the index n, respectively [31–34]. The corresponding absorption and scattering cross-sections, σabs and σscat, defined as the total power absorbed and scattered by a single scatterer with respect to the unit intensity of incident plane wave, can be expressed as

σabsn=1σnabs(TE)+σnabs(TM)=n=1(2n+1)λ22π(Re{CnTM}+|CnTM|2+Re{CnTE}+|CnTE|2),
σscat=n=1n=(2n+1)λ22π(|CnTM|2+|CnTE|2),
where λ is the wavelength of incident wave in vacuum. For a given radius of particle, denoted as a, the value of size parameter 2πa/λ determines how many terms in these two convergent series to be dominant [32]. Here, we define the partial absorption cross section for each spherical harmonic channel, labeled by n, as σnabs(TE,TM)(2n+1)λ22π(Re{Cn(TE,TM)}+|Cn(TE,TM)|2), for TE or TM mode. We further express the scattering coefficient as Cn(TE,TM)=|Cn(TE,TM)|exp{iθn(TE,TM)}, where the magnitude |Cn(TE,TM)| is a positive real value and θn(TE,TM) is the corresponding phase. Due to the conservation of power, these partial absorption cross-sections would be equal or larger than zero in each spherical harmonic channel, σnabs(TE,TM)0.

By decomposing into partial waves, in terms of each spherical channel, we can have a universal phase diagram for any passive electromagnetic scatterers, as shown in Fig. 1. Interestingly, even though we do not write down any exact formulas for the scattering coefficient, the range to support physical values for the amplitude square only exists within 0|Cn(TE,TM)|21; while the phase is bound within π/2θn(TE,TM)3π/2, as illustrated in Fig. 1. It is worth to remark that the range to support allowable solutions of scattering coefficients results from the intrinsic power conservation regardless of any specific scattering events. We depict the allowed solutions in colors for passive cases, which correspond to σnabs(TE,TM)0 and depict the forbidden ones in white color for gain cases, which correspond to σnabs(TE,TM)<0. Along this lossless contour, there exist a family of solutions with the same value on σnabs(TE,TM)=0, but with different scattering coefficients in amplitude and phase. It is known that an ideally localized surface plasmon in the subwavelength structure relies on lossless resonance condition [3, 5], which corresponds to the point |Cn(TE,TM)|=1 and θn(TE,TM)=π in our phase diagram. As for invisible cloaks [11–16], one can look for the solutions near bottom of the phase diagram, i.e., |Cn(TE,TM)|=0, for dominant channels.

 figure: Fig. 1

Fig. 1 A phase diagram for each spherical harmonic channel, labeled by n, is generated by imposing the power conservation on the partial absorption cross section, for TE or TM mode separately. Marked numbers shown in the contour lines correspond to the values of normalized absorption cross section in the individual channel: 2π(2n+1)λ2σnabs(TE,TM). Colored regions are physically allowed solutions; while uncolored regions represent forbidden solutions. It is noted that the amplitude square is bounded within the range [0, 1]; while the allowed phase is within [π/2, 3π/2]. The Green cross-marker, localed at (θn(TE,TM)=π,|Cn(TE,TM)|2=0.25), indicates the maximum value, 0.25, in the normalized absorption cross-section.

Download Full Size | PDF

Once the composited material in a scatterer has intrinsic loss, the scattering coefficients move to reside inside the colored region. For each channel, the maximum value in the normalized absorption cross section is 2πσnabs/(2n+1)λ2=1/4, i.e., the Green cross-marker shown in Fig. 1, corresponding to coherent perfect absorbers [7–10], but which is also associated with the same amount of electromagnetic scattering power. The phase and amplitude of scattering coefficients to achieve a maximum absorption power is π and 1/2, respectively. Moreover, along the contour for a constant absorption power, there exist a maximum and a minimum values in the scattering amplitude, both at the phase θn(TE,TM)=π. It implies that one may design a scatterer possessing the same absorption power, but with different scattering signals. As for seeking an optimized “cloaking a sensor” [30], i.e., to have a better absorption efficiency with a constant absorption power, in the phase diagram one can follow the constant absorption contour to find the corresponding solution located at θn(TE,TM)=π (conjugate-matched condition) with a minimum scattering amplitude.

3. From the trajectory in the phase diagram to design passive scatterers

Through above examples, the phase diagram provides a universal map to display all possible solutions for any passive scatterers. In principle, without knowing the composition in a scatterer, one can also have the same scattering coefficients. In this way, one may design specific scatterers with the required scattering and absorption properties by the choice of allowed trajectories in the phase diagram.

In the following, we introduce a systematic way to do the inverse design for the scatterer by specifying the required scattering or absorption properties. For a scatterer with N layers made of isotropic and homogeneous media, the corresponding scattering coefficient can be expressed in a compact form:

Cn(TE,TM)=11+iVn(TE,TM)Un(TE,TM)11+i[αn(TE,TM)+iβn(TE,TM)],
where Un(TE,TM) and Vn(TE,TM) are determinants of a 2N × 2N matrix constituted by spherical harmonic functions [5,11,15,24]. Here, we also rewrite this scattering coefficient by introducing two auxiliary real numbers: αn(TE,TM) and βn(TE,TM) in Eq. (3). By substituting Eq. (3) into the scattering coefficients, one can have different trajectories in the phase diagram, as illustrated in Fig. 2. Supported trajectories for the parameter sets with a constant value of αn(TE,TM) or βn(TE,TM) are plotted in Red dotted-dashed- and Blue dotted-dashed-curves, respectively. The ranges to support allowable scattering solutions for passive scatterers on these numbers would be αn(TE,TM)=[,] and βn(TE,TM)=[,0].

 figure: Fig. 2

Fig. 2 Supported trajectories in the phase diagram are shown for different sets of the parameters: αn(TE,TM) and βn(TE,TM) defined in Eq. (3). Here, trajectories with a constant βn(TE,TM) are shown in Blue dotted-dashed-curves; while trajectories with a constant αn(TE,TM) are shown in Red dotted-dashed-curves. Two contours for a constant absorption power are also depicted in the Black color.

Download Full Size | PDF

In particular, for the contour with a constant absorption power in the phase diagram, as shown in the Black color in Fig. 2, one can use the parametric representation of the curve to describe this supported trajectory. The corresponding parameter sets of αn(TE,TM) and βn(TE,TM) can be found as

αn(TE,TM)(t)=14[qn(TE,TM)]21qn(TE,TM)sin(t),
βn(TE,TM)(t)=[112qn(TE,TM)]+14[qn(TE,TM)]21qn(TE,TM)cos(t).

Here, the independent variable, t, is used for the parametric equation, which is also bounded within the range t = [0,2π]. For a given normalized absorption power, qn(TE,TM)2π(2n+1)λ2σnabs(TE,TM) in each channel, from Eqs. (4)(5) one can see that the trajectory of αn(TE,TM) and βn(TE,TM) in the phase diagram is following the request for a constant absorption power.

As an example, we consider a passive scatterer in the configuration of a core-shell sphere, as illustrated in the insect of Fig. 3, which is composed by two concentric layers of isotropic and homogeneous materials. The geometrical parameters and material properties for this core-shell scatterer are the radius of core, ac, the radius of whole particle, a, and εs(μs)/εc(μc) for the permittivity (permeability) in the shell/core regions, respectively. We limit our system to nonmagnetic case, so μs = μc = 1. If the electrically small approximation is satisfied for such a core-shell scatterer, it is known that the main contribution dominantly comes from the electric dipole-wave scattering, i.e., n = 1 and TM mode. We choose the constant absorption power with q1(TM)=0.2 in Eqs. (4)(5). For such a two-layered scatterer, the corresponding scattering coefficients are conducted from a 4 × 4 matrix by tracking TE and TM modes. By applying the continuity of electric and magnetic fields established at the two boundaries of shell-environment and core-shell, one can approximately express the term V1TM/U1TM as

V1TMU1TM=3λ32(2πa)32γ3(1εs)(εcεs)(2+εs)(εc+2εs)γ3(εsεc)(2εs+1)+(1εs)(εc+2εs),
where γ is defined as the ratio between the core radius to the whole particle radius, γ ≡ ac/a. If one replaces ε by μ, then we can obtain the other term V1TE/U1TE. However for non-magnetic media it is automatically zero for μ0 = μs = μc = 1. By taking γ = 1 or εs = εc, above result can be reduced to the electric dipole equation for a solid sphere.

 figure: Fig. 3

Fig. 3 Absorption and scattering cross sections correspond to the contour shown in Fig. 2, which are depicted in terms of the parametric variable t defined in Eqs. (4)(5). Here, a constant absorption power is requested by setting q1TM=0.2; while there is a degree of freedom in the scattering power. The insect illustrates the core-shell scatterer used as an example to design a passive electromagnetic devices with the constant absorption power.

Download Full Size | PDF

Now, for our core-shell system with the geometric size fixed, we provide a systematic way to find out the corresponding material properties with a constant absorption power, as specified by the contour in the phase diagram shown in Fig. 2. To give a clear illustration, first, one may fix the material property in the shell or in the core region. If we assume that the composition for the shell region is given, i.e., εs is fixed, then, based on Eqs. (4)(6), the corresponding solution for the permittivity in the core region is found to satisfy:

εc=εs3(2εs42γ3+2γ3εs)2(α1TM+iβ1TM)(2πa/λ)3(22εs+2εsγ3+γ3)3(εs+2+2γ3εs2γ3)+2(α1TM+iβ1TM)(2πa/λ)3(1εs2εsγ3γ3).

Solutions obtained from the analytical formula in Eq. (7) are shown in Figs. 4(a) and 4(b) for the real and imaginary parts of the permittivity in the core region, respectively. In terms of the parametric variable, t, we can have a wide rang in selecting core materials, and all of them have the same absorption power. Moreover, based on these found parameters, the corresponding absorption and scattering cross sections would satisfy our request for a constant absorption power, as shown in Fig. 3. From the comparison between Fig. 3 and Figs. 4(a)–4(b), we find that when t = π/2 the scattering power reaches a maximum value; while the εc for the required material has a minimum value, due to the reason that dissipative loss is proportional to the local electric field. In this scenario, with the help of a strong electric field, it becomes possible to maintain the same absorption power simultaneously.

 figure: Fig. 4

Fig. 4 The permittivities to support a constant absorption power are shown as a function of the parametric variable t. For a given material in the shell region, εs = 3.12, found solutions for the real and imaginary parts of the permittivity in the core region are shown in (a) and (b), respectively. For a given material in the core region, εc = 5, two families of found solutions shown in Eq. (8) are denoted as εs+ and εs for the shell region, with the corresponding real and imaginary parts of the permittivity shown in (c, e) and (d, f), respectively. Results obtained from analytical formulas are depicted in solid-curves; while exact solutions from scattering theory are depicted in dashed-curves. In all cases, the core-shell geometries are fixed with a = 1/24λ and γ = 0.9.

Download Full Size | PDF

On other hand, if the material property in the core region is specified, based on Eqs. (4)(6) the corresponding solutions for the permittivity in the shell region εs to support a constant absorption power are governed by

εs±=g±g24fh2f,
with the shorthanded notations:
f=2(1γ3)[32(α1TM+iβ1TM)(2πa/λ)3],
g=2(α1TM+iβ1TM)(2πa/λ)3[γ3(12εc)+2εc]+3(2γ3+2γ3εc+εc+4),
h=εc(1γ3)[6+2(α1TM+iβ1TM)(2πa/λ)3].

Derivations to have solutions in Eqs. (7) and (8) are shown in Appendix.

From Eq. (8), there exist two families to support the materials in the shell region, denoted as εs±. We show the real and imaginary parts of the permittivity in the shell region for these two families in Figs. 4(c)–4(d) and 4(e)–4(f), as well as the exact solutions in dashed-curves, respectively. We reveal a good agreement between our analytical solutions and the numerical ones obtained by the exact scattering theory. Again, based on these results, the corresponding absorption and scattering cross sections would satisfy our request for a constant absorption power, as shown in Fig. 3. It is surprised to find that there exist a variety of choices for the material properties even a specific scattering or absorption property is given in the beginning. Without the introduction of this phase diagram and the inverse design, it is not only difficult to find out the required material properties, but also complicated to recognize the power competition and limitation among these cross sections for each channel.

Before concluding, we remark that when the electrically small approximation is not valid, multiple channels for TE and TM modes may be excited as expected. In this scenario, one can also apply our phase diagram, but not for a single channel only. By embedding multi-layered coatings to excite multiple channels, the intrinsic single channel limitation can be broken to generate superscattering or superabsorber phenomena [17–21, 29]. Scattering coefficients from several dominant channels are just a natural extension by considering all of them onto the phase diagrams simultaneously. In addition, although the proposed phase diagram is based on the well-known scattering formulas for symmetrically spherical scatterers, the concept of our phase diagram can be applied to non-spherical scatterers as well. Although in our example of inverse design we use well-known electric dipole formulas for core-shell systems, our approach to find out the corresponding materials for a given request on the scattering and absorption properties is always non-trivial. Through this universal phase diagram, one can have a systematic way to design functional passive electromagnetic scatterers.

4. Conclusion

In summary, we introduce a phase diagram as a compact tool to link the scattering and absorption powers for each spherical harmonic channel. Intrinsically, the power conservation for any passive scatterers gives the physically allowable solutions in the scattering coefficient. Not only the known exotic scattering and absorption phenomena can be illustrated in this diagram, but supported trajectories are also demonstrated to design extrinsic-field-controllable scatterers. With the core-shell scatterers at the subwavelength scale as an example, we reveal a systematic way to find out a variety of solutions in the composited materials to possess the same absorption power. In general, one can easily extend this methodology to go beyond the small particle size limitation, by considering interferences from several channels in the map. With the analogy among wave phenomena, the concept of this phase diagram and our inverse design method can be ready applied to acoustic systems as well as quantum scattering system.

Appendix

In this Appendix, we provide the formula to derive Eqs. (7) and (8) in detail. When the electrically small condition (quasi-static limit) is satisfied in a non-magnetic core-shell structure, the dominant scattered wave can be approximated by Eq. (6) [5, 11, 15, 24]. Then, by defining Z1TMα1TM+iβ1TM as a shorthanded notation, Eq. (6) becomes

V1TMU1TM=Z1TM=32(λ2πa)32γ3(1εs)(εcεs)(2+εs)(εc+2εs)γ3(εsεc)(2εs+1)+(1εs)(εc+2εs).

With the help of Eq. (12), one has

εc{2(2πaλ)3Z1TM[2εsγ3γ3+1εs]3[2γ32γ3εs2εs]}=3[4εs2εs22γ3εs+2γ3εs2]2(2πaλ)3Z1TM[γ3(2εs2+εs)+(2εs2εs2)],
which gives us the result in Eq. (7). Similar process can be applied to derive Eq. (8) by expanding every terms and collecting coefficients of εs, i.e.,
εs2(2γ32)[2(2πaλ)3Z1TM3]+[2(2πaλ)3Z1TMεc(1γ3)6εc(γ31)]+εs[2(2πaλ)3Z1TM(γ32εcγ3+2εc)+3(2γ3+2γ3εc+4+εc)]=0.

Finally, with the shorthanded notations introduced in Eqs. (9)(11), we have

fεs2+gεs+h=0.

As a result, one can easily solve εs in Eq. (15), which gives the solutions shown in Eq. (8).

Acknowledgments

This work is supported in part by the Ministry of Science and Technologies, Taiwan, under the contract No. 101-2628-M-00 7-003-MY4.

References and links

1. J. Shi, F. Monticone, S. Elias, Y. Wu, D. Ratchford, X. Li, and A. Alú, “Modular assembly of optical nanocircuits,” Nat. Commun. 5, 3896 (2014). [CrossRef]   [PubMed]  

2. N. Engheta, “Circuits with light at nanoscales: optical nanocircuits inspired by metamaterials,” Science 317, 1698–1702 (2007). [CrossRef]   [PubMed]  

3. M. I. Tribelsky and B. S. Lukyanchuk, “Anomalous light scattering by small particles,” Phys. Rev. Lett. 97, 263902 (2006). [CrossRef]  

4. F. Monticone, C. Argyropoulos, and A. Alú, “Multi-layered plasmonic covers for comblike scattering response and optical tagging,” Phys. Rev. Lett. 110, 113901 (2013). [CrossRef]  

5. A. Alú and N. Engheta, “Polarizabilities and effective parameters for collections of spherical nanoparticles formed by pairs of concentric double-negative, single-negative, and/or double-positive metamaterial layers,” J. Appl. Phys. 97, 094310 (2005). [CrossRef]  

6. R. E. Hamam, A. Karalis, J. D. Joannopoulos, and M. Soljacic, “Coupled-mode theory for general free-space resonant scattering of waves,” Phys. Rev. A 75, 053801 (2007). [CrossRef]  

7. H. Noh, Y. Chong, A.D. Stone, and H. Cao, “Perfect coupling of light to surface plasmons by coherent absorption,” Phys. Rev. Lett. 108, 186805 (2012). [CrossRef]   [PubMed]  

8. H. Noh, S.M. Popoff, and H. Cao, “Broadband subwavelength focusing of light using a passive sink,” Opt. Express 21, 17435–17446 (2013). [CrossRef]   [PubMed]  

9. M. I. Tribelsky, “Anomalous light absorption by small particles,” Europhys. Lett. 94, 14004 (2011). [CrossRef]  

10. V. Grigoriev, N. Bonod, J. Wenger, and B. Stout, “Optimizing nanoparticle designs for ideal absorption of light,” ACS Photonics 2, 263–270 (2015). [CrossRef]  

11. A. Alú and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E 72, 016623 (2005). [CrossRef]  

12. A. Alú and N. Engheta, “Plasmonic materials in transparency and cloaking problems: mechanism, robustness, and physical insights,” Opt. Express 15, 3318–3332 (2007). [CrossRef]   [PubMed]  

13. S. Muhlig, M. Farhat, C. Rockstuhl, and F. Lederer, “Cloaking dielectric spherical objects by a shell of metallic nanoparticles,” Phys. Rev. B 83, 195116 (2011). [CrossRef]  

14. S. Muhlig, A. Cunningham, J. Dintinger, M. Farhat, S. B. Hasan, T. Scharf, T. Burgi, F. Lederer, and C. Rockstuhl, “A self-assembled three-dimensional cloak in the visible,” Sci. Rep. 3, 2328 (2013). [CrossRef]   [PubMed]  

15. M. Farhat, S. Mhlig, C. Rockstuhl, and F. Lederer, “Scattering cancellation of the magnetic dipole field from macroscopic spheres,” Opt. Express 20, 13896–13906 (2012). [CrossRef]   [PubMed]  

16. A. Alú and N. Engheta, “Cloaking a sensor,” Phys. Rev. Lett. 102, 233901 (2009). [CrossRef]   [PubMed]  

17. Z. Ruan and S. Fan, “Superscattering of light from subwavelength nanostructures,” Phys. Rev. Lett. 105, 013901 (2010) [CrossRef]   [PubMed]  

18. Z. Ruan and S. Fan, “Design of subwavelength superscattering nanospheres,” Appl. Phys. Lett. 98, 043101 (2011). [CrossRef]  

19. A. Mirzaei, A. E. Miroshnichenko, I. V. Shadrivov, and Y. S. Kivshar, “Superscattering of light optimized by a genetic algorithm,” Appl. Phys. Lett. 105, 011109 (2014). [CrossRef]  

20. N. M. Estakhri and A. Alú, “Minimum-scattering superabsorbers,” Phys. Rev. B. 89, 121416 (2014). [CrossRef]  

21. Z. Ruan and S. Fan, “Temporal coupled-mode theory for Fano resonance in light scattering by a single obstacle,” J. Phys. Chem. C 114, 7324 (2010). [CrossRef]  

22. K. R. Catchpole and A. Polman, “Plasmonic solar cells,” Opt. Express 16, 21793–21800 (2008). [CrossRef]   [PubMed]  

23. X. Zhang, Y. L. Chen, R. S. Liu, and D. P. Tsai, “Plasmonic photocatalysis,” Rep. Prog. Phys. 76, 046401 (2013). [CrossRef]   [PubMed]  

24. J.-Y. Lee, M.-C. Tsai, P.-C. Chen, T.-T. Chen, K.-L. Chan, C.-Y. Lee, and R.-K. Lee, “Thickness effects on light absorption and scattering for nano-particles in shape of hollow-spheres,” J. Phys. Chem. C 119, 25754–25760 (2015). [CrossRef]  

25. L. R. Hirsch, R. J. Stafford, J. A. Bankson, S. R. Sershen, B. Rivera, R. E. Price, J. D. Hazle, N. J. Halas, and J. L. West, “Nanoshell-mediated near-infrared thermal therapy of tumors under magnetic resonance guidance,” Proc. Natl. Acad. Sci. U. S. A. 100, 13549–13554 (2003). [CrossRef]   [PubMed]  

26. M. I. Tribelsky, A. E. Miroshnichenko, Y. S. Kivshar, B. S. Lukyanchuk, and A. R. Khokhlov, “Laser pulse heating of spherical metal particles,” Phys. Rev. X 1, 021024 (2011).

27. Y. Pu, R. Grange, C. L. Hsieh, and D. Psaltis, “Nonlinear optical properties of core-shell nanocavities for enhanced second-harmonic generation,” Phys. Rev. Lett. 104, 207402 (2010). [CrossRef]   [PubMed]  

28. S.-W. Chu, T.-Y. Su, R. Oketani, Y.-T. Huang, H.-Y. Wu, Y. Yonemaru, M. Yamanaka, H. Lee, G.-Y. Zhuo, M.-Y. Lee, S. Kawata, and K. Fujita, “Measurement of a saturated emission of optical radiation from gold nanoparticles: Application to an ultrahigh resolution microscope,” Phys. Rev. Lett. 112, 017402 (2014). [CrossRef]   [PubMed]  

29. A. Mirzaei, I. V. Shadrivov, A. E. Miroshnichenko, and Y. S. Kivshar, “Cloaking and enhanced scattering of core-shell plasmonic nanowires,” Opt. Express 21, 10454–10459 (2013). [CrossRef]   [PubMed]  

30. R. Fleury, J. Soric, and A. Alú, “Physical bounds on absorption and scattering for cloaked sensors,” Phys. Rev. B. 89, 045122 (2014). [CrossRef]  

31. M. Born and E. Wolf, Principle of Optics, 7th ed. (Cambridge University, 1999). [CrossRef]  

32. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).

33. S. A. Maier, Plasmonics: Fundamentals and Applications (Springer, 2007).

34. H.C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1
Fig. 1 A phase diagram for each spherical harmonic channel, labeled by n, is generated by imposing the power conservation on the partial absorption cross section, for TE or TM mode separately. Marked numbers shown in the contour lines correspond to the values of normalized absorption cross section in the individual channel: 2 π ( 2 n + 1 ) λ 2 σ n a b s ( T E , T M ). Colored regions are physically allowed solutions; while uncolored regions represent forbidden solutions. It is noted that the amplitude square is bounded within the range [0, 1]; while the allowed phase is within [π/2, 3π/2]. The Green cross-marker, localed at ( θ n ( T E , T M ) = π , | C n ( T E , T M ) | 2 = 0.25 ), indicates the maximum value, 0.25, in the normalized absorption cross-section.
Fig. 2
Fig. 2 Supported trajectories in the phase diagram are shown for different sets of the parameters: α n ( T E , T M ) and β n ( T E , T M ) defined in Eq. (3). Here, trajectories with a constant β n ( T E , T M ) are shown in Blue dotted-dashed-curves; while trajectories with a constant α n ( T E , T M ) are shown in Red dotted-dashed-curves. Two contours for a constant absorption power are also depicted in the Black color.
Fig. 3
Fig. 3 Absorption and scattering cross sections correspond to the contour shown in Fig. 2, which are depicted in terms of the parametric variable t defined in Eqs. (4)(5). Here, a constant absorption power is requested by setting q 1 T M = 0.2; while there is a degree of freedom in the scattering power. The insect illustrates the core-shell scatterer used as an example to design a passive electromagnetic devices with the constant absorption power.
Fig. 4
Fig. 4 The permittivities to support a constant absorption power are shown as a function of the parametric variable t. For a given material in the shell region, εs = 3.12, found solutions for the real and imaginary parts of the permittivity in the core region are shown in (a) and (b), respectively. For a given material in the core region, εc = 5, two families of found solutions shown in Eq. (8) are denoted as ε s + and ε s for the shell region, with the corresponding real and imaginary parts of the permittivity shown in (c, e) and (d, f), respectively. Results obtained from analytical formulas are depicted in solid-curves; while exact solutions from scattering theory are depicted in dashed-curves. In all cases, the core-shell geometries are fixed with a = 1/24λ and γ = 0.9.

Equations (15)

Equations on this page are rendered with MathJax. Learn more.

σ abs n = 1 σ n abs ( TE ) + σ n abs ( TM ) = n = 1 ( 2 n + 1 ) λ 2 2 π ( Re { C n TM } + | C n TM | 2 + Re { C n TE } + | C n TE | 2 ) ,
σ scat = n = 1 n = ( 2 n + 1 ) λ 2 2 π ( | C n TM | 2 + | C n TE | 2 ) ,
C n ( T E , T M ) = 1 1 + i V n ( T E , T M ) U n ( T E , T M ) 1 1 + i [ α n ( T E , T M ) + i β n ( T E , T M ) ] ,
α n ( T E , T M ) ( t ) = 1 4 [ q n ( T E , T M ) ] 2 1 q n ( T E , T M ) sin ( t ) ,
β n ( T E , T M ) ( t ) = [ 1 1 2 q n ( T E , T M ) ] + 1 4 [ q n ( T E , T M ) ] 2 1 q n ( T E , T M ) cos ( t ) .
V 1 T M U 1 T M = 3 λ 3 2 ( 2 π a ) 3 2 γ 3 ( 1 ε s ) ( ε c ε s ) ( 2 + ε s ) ( ε c + 2 ε s ) γ 3 ( ε s ε c ) ( 2 ε s + 1 ) + ( 1 ε s ) ( ε c + 2 ε s ) ,
ε c = ε s 3 ( 2 ε s 4 2 γ 3 + 2 γ 3 ε s ) 2 ( α 1 T M + i β 1 T M ) ( 2 π a / λ ) 3 ( 2 2 ε s + 2 ε s γ 3 + γ 3 ) 3 ( ε s + 2 + 2 γ 3 ε s 2 γ 3 ) + 2 ( α 1 T M + i β 1 T M ) ( 2 π a / λ ) 3 ( 1 ε s 2 ε s γ 3 γ 3 ) .
ε s ± = g ± g 2 4 f h 2 f ,
f = 2 ( 1 γ 3 ) [ 3 2 ( α 1 T M + i β 1 T M ) ( 2 π a / λ ) 3 ] ,
g = 2 ( α 1 T M + i β 1 T M ) ( 2 π a / λ ) 3 [ γ 3 ( 1 2 ε c ) + 2 ε c ] + 3 ( 2 γ 3 + 2 γ 3 ε c + ε c + 4 ) ,
h = ε c ( 1 γ 3 ) [ 6 + 2 ( α 1 T M + i β 1 T M ) ( 2 π a / λ ) 3 ] .
V 1 T M U 1 T M = Z 1 T M = 3 2 ( λ 2 π a ) 3 2 γ 3 ( 1 ε s ) ( ε c ε s ) ( 2 + ε s ) ( ε c + 2 ε s ) γ 3 ( ε s ε c ) ( 2 ε s + 1 ) + ( 1 ε s ) ( ε c + 2 ε s ) .
ε c { 2 ( 2 π a λ ) 3 Z 1 T M [ 2 ε s γ 3 γ 3 + 1 ε s ] 3 [ 2 γ 3 2 γ 3 ε s 2 ε s ] } = 3 [ 4 ε s 2 ε s 2 2 γ 3 ε s + 2 γ 3 ε s 2 ] 2 ( 2 π a λ ) 3 Z 1 T M [ γ 3 ( 2 ε s 2 + ε s ) + ( 2 ε s 2 ε s 2 ) ] ,
ε s 2 ( 2 γ 3 2 ) [ 2 ( 2 π a λ ) 3 Z 1 T M 3 ] + [ 2 ( 2 π a λ ) 3 Z 1 T M ε c ( 1 γ 3 ) 6 ε c ( γ 3 1 ) ] + ε s [ 2 ( 2 π a λ ) 3 Z 1 T M ( γ 3 2 ε c γ 3 + 2 ε c ) + 3 ( 2 γ 3 + 2 γ 3 ε c + 4 + ε c ) ] = 0.
f ε s 2 + g ε s + h = 0.
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.