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

Light absorption in conical silicon particles

Open Access Open Access

Abstract

The problem of the absorption of light by a nanoscale dielectric cone is discussed. A simplified solution based on the analytical Mie theory of scattering and absorption by cylindrical objects is proposed and supported by the experimental observation of sharply localized holes in conical silicon tips after high-fluence irradiation. This study reveals that light couples with tapered objects dominantly at specific locations, where the local radius corresponds to one of the resonant radii of a cylindrical object, as predicted by Mie theory.

©2013 Optical Society of America

The absorption and scattering of light by objects with subwavelength dimensions is one of the most fundamental problems in electrodynamics. The strong sensitivity of the optical properties of such objects to their shape and size has also taken this problem to the forefront of many technological advances and challenges (Refs. 1,2 and references therein). More particularly, owing to their combined focusing and guiding properties, particles with a conical shape have led to many successful applications such as nanoscale optical tweezers or quasi-non-dispersive antireflection coatings [3,4], but have surprisingly received little theoretical attention comparatively to their spherical and cylindrical counterparts. This is most likely due to the problematic presence of a singularity at the cone apex [5], combined with their large aspect ratio [6]. Nevertheless, the relevance of this problem has recently arisen in the fields of metrology and microscopy [7], and, more particularly, in the context of the laser-assisted Atom Probe Tomography (APT) technique [8,9]. In this paper, we first uncover the interaction between light and infinite axisymmetric particles by extending Mie theory of light absorption by cylindrical particles [10,11], assuming that an axisymmetric particle is a mere stack of infinitesimal cylinders with a varying radius. Focusing our study on the case of a conical particle, we then deduce theoretically and show experimentally that the absorption of light occurs preferentially at specific locations on the particle, where the local radius corresponds to one of the resonance absorption/scattering radii of an illuminated cylinder at the considered wavelength. The impact of the finiteness of the particle, i.e. of the light diffracted at the apex and fed into the particle, is also discussed.

As we shall show, the light absorption by conical nanoparticles can be understood based on the analytical Mie theory of absorption by infinitely long cylindrical particles [10,11]. Since the absorbed power is proportional to the (2n)th power [n = number of photon(s) involved in the transition] of the electric field Eint induced inside the particle by the incident electromagnetic wave [12], this absorption problem can be summarized by the expression of Eint(r,ϑ), where rand ϑ are respectively the radial and azimuthal coordinates in the local axis of the cylinder cross-section. For a transverse magnetic (TM) plane wave incident normally to the axis z of symmetry of a circular cylinder of radius R, Eint(r,ϑ) is oriented along z and its amplitude reads [10,11]

Ez,int(r,ϑ)=E0[d0(Rλ0)J0(2πn˜λ0r)+k=12ikdk(Rλ0)Jk(2πn˜λ0r)cos(kϑ)],
where E0 and λ0 are respectively the amplitude and free-space wavelength of the incident electric field, n˜is the complex refractive index of the particle and Jk is thekth-order Bessel function of first kind [13]. The dk(R/λ0) coefficients are determined by the boundary conditions along the surface r=R [10,11]. A factor exp(iωt) expressing the time dependence of the wave is assumed, ω being the angular frequency of the incident wave. To illustrate Eq. (1), Fig. 1 plots the behavior of |Ez,int| in a cylindrical silicon particle as a function of the cylinder radius when illuminated with green (λ0=515 nm) and infrared (IR, λ0=1030 nm) light. The complex refractive index follows the experimental data of [14]. Figure 1 shows that, as a result of the resonant behavior of the dk(R/λ0) coefficients, |Ez,int| also presents resonant values – or maxima - in cylinders with radii of specific values. The position of these maxima as well as their spacings depend upon the incident wavelength. This resonance behavior can be understood as the result of the constructive interference between the transmitted and (multiply) internally reflected waves inside the cylinder, as also occurs in thin films of specific thicknesses [15]. Interestingly, however, the resonant thicknesses of a slab differ somewhat from the resonant diameters of a cylinder, the latter proving to be roughly 0.257×λ0 smaller than the former [16].

 figure: Fig. 1

Fig. 1 Resonance behavior of the electric field in a cylindrical silicon particle of radiusR, illuminated with a TM wave of amplitude E0 and wavelengths λ0=515 nm (full line) and λ0=1030 nm (dotted line), incident normally to the cylinder axis. The shown values are on the illuminated side (r=R,ϑ=π).

Download Full Size | PDF

We propose to generalize Eq. (1) to the case of an infinite axisymmetric particle with an arbitrary surface R(z) illuminated normally to its axis (referred to as normal incidence in the rest of this paper). This is physically acceptable provided the resulting z-dependence of the dk[R(z)/λ0] coefficients is such that |2Ez,int/z2|<<|2πn˜/λ0|2|Ez,int|. When this condition is verified, the wave equation determining Ez,int is indeed hardly perturbed by the z-dependence of the dk coefficients. In other words, if the particle surface is smooth enough along its symmetry axis, |2Ez,int/z2| is small and the particle can be considered as locally cylindrical, i.e. the coupling between cross-sections with different radii can be neglected. An extra benefit of this condition for the normally illuminated cones considered in this paper is that each local cylinder is also illuminated normally to its axis. Consequently, the contributions to the internal field oriented along the r and ϑ directions can be expected to be small relative toEz,intand thus be ignored. Considering a conical particle with cone angle α and a surface defined by

R(z)=ztan(α/2),
the above condition reduces to
|Ez,int|>>|λ0tan(α/2)2πn˜|2|2Ez,intR2|.
Figure 2 shows the value of the right-hand side of Eq. (3) in the case of a conical silicon particle with cone angle α=10° illuminated respectively by green and IR TM waves of amplitudesE0 incident normally to the cone axis. Comparison with Fig. 1 shows that Eq. (3) is verified everywhere for green illumination. For IR light, however, Eq. (3) is violated in very narrow regions of the cone close to some of the resonant cross-sections (R(z)175,240,300,...nm). These peaked perturbations generate localized source terms in the wave equation, which should invalidate the use of Eq. (1) for the conical particle. These perturbations, nonetheless, are drastically reduced when the cone angle decreases and, as we show experimentally in the sequel, do not seem to critically impact the applicability of Eq. (1) to conical particles.

 figure: Fig. 2

Fig. 2 Variation in the right-hand side of Eq. (3) as a function of the local radius R(z)in a conical silicon particle with cone angle α=10° illuminated respectively by green (full line) and IR (dotted line) TM waves of amplitudeE0, incident normally to the cone axis. The shown values are on the illuminated side (r=R,ϑ=π).

Download Full Size | PDF

To experimentally validate our theory summarized by Eqs. (1) and (2), conical silicon tips have been irradiated with high-fluence (~100 mJ/cm2) TM IR light in a Laser-Assisted Wide-Angle Tomographic Atom Probe (LAWATAP) tool of Cameca [17]. The Scanning Electron Microscope (SEM) image of Fig. 3 shows that, after irradiation, the tips present sharp holes, indicating strongly localized absorption. The calculated value of |Ez,int[r=R(z),ϑ=π,z]|4 obtained with Eqs. (1) and (2) for silicon irradiated with IR light is reproduced on the left of Fig. 3, making it obvious that these holes are precisely located at the theoretically predicted maxima of |Ez,int|4. Similar holes are also observed in the case of green irradiation but are packed much more tightly together (Fig. 1), making their discussion more difficult (not shown). The sharp localization of the holes at the resonance radii predicted by Mie theory for the cylindrical particle proves that the light is coupled into the cone dominantly at these locations, hence confirming the validity of the present theory. The relatively large cone angle (α~10°) of the tip in Fig. 3, also suggests that the theory remains valid, at least in first order, even when perturbations of the wave equation induced by a non-negligible |2Ez,int/z2| are present.

 figure: Fig. 3

Fig. 3 (left) Theoretical value of |Ez,int[r=R(z),ϑ=π,z]|4 along the illuminated side of an IR-irradiated conical silicon tip as a function of the local radius R(z). (right) Silicon tip irradiated with a high-fluence (~100 mJ/cm2) IR pulsed laser. Resonance absorption occurs in the cross-sections with R(z)~113, 178 and 241 nm. The tip is tilted 45° with respect to the primary beam of the SEM.

Download Full Size | PDF

Although the agreement between theory and experiment in Fig. 3 is clear, three important remarks need to be made for its complete and unambiguous understanding. Firstly, it has to be stressed that the holes are sharply located inside each resonant cross-section of the cone. Their apparent elliptical shape in Fig. 3 is merely due to the 45° tilt angle of the primary beam of the SEM with respect to the tip axis. This tilt is also responsible for the small (downward) apparent misalignment of the holes with respect to the corresponding resonance peaks. Actually, only the leftmost point of each hole is reliable for measurement of the local radius R(z) of the holes (locations of the 3 measured values of 113, 178 and 241 nm). Secondly, the use of the fourth power of |Ez,int| is justified by the dominant two-photon absorption process of high-fluence IR radiation by silicon [18]. For the sake of completeness, note also that the inclusion of this nonlinear effect in the value of the complex refractive index in Eq. (1) proves to slightly reduce the heights of the resonance |Ez,int| peaks while keeping them at the same position (not shown). Thirdly and finally, the reason for the one-sided shape of the observed holes is still under discussion, with two competing possible explanations. This shape may indeed be due to a transient nonlinear absorption phenomenon common to pulsed-laser ablation of Si, wherein free-carrier absorption becomes very intense and dominant in the course of the laser pulse, hence strongly decreasing the light penetration depth [19]. It is, however, also suspected that the hole shape may result from an optical interference between Ez,int and the perturbation field induced by the local |2Ez,int/z2| highlighted in Fig. 2. This question is still under investigation.

The impact of the finiteness of a particle such as shown in Fig. 3 remains to be determined. The agreement between the developed theory for infinitely long particles and the experimental observation of holes summarized in Fig. 3 seems to infer that the apex does not have much impact upon its absorption characteristics in the case of incidence normal to the cone axis. The same conclusion is reached in recent theoretical studies of the electromagnetic field inside an illuminated truncated cylinder with planar endings [20]. Let us, however, investigate the amplitude Ez,apex of the light diffracted at a rounded apex, considered here as an additional independent contribution by invoking the linearity of Maxwell’s wave equations. To understand why Ez,apex is likely to be negligible, two determining factors must be investigated, namely the coupling efficiency of the apex and the guiding ability of the cone. First, to determine how efficiently light can be coupled into the apex, let us consider the hemispherical apex as an isolated silicon sphere. It is then easy to show, based on Mie theory for spherical particles [10], that only a few percent of the incident IR light is coupled into the apex when its radius Rapex<<145 nm, i.e. far from the first resonance peak of the spherical particle. Note, furthermore, that similar conclusions would be reached using Bethe-Bouwkamp diffraction theory by assimilating the apex as a hole in a conductive screen [21]. In other words, it can be understood that the approximate 35-nm apex radius of the silicon nanoparticle shown in Fig. 3 is too small to significantly couple to IR light. Second, let us study how much electromagnetic power a conical particle, acting as a dielectric fiber [22], can guide away from the apex. In such a waveguiding structure, indeed, only bound modes propagate with a significant electric field inside the guide and therefore have a chance of being partially absorbed. The relevant parameter to be tested is the normalized frequency V=2πRapex/λ0((n˜))21 [22]. Noting that the condition for single-mode operation V<2.4 is reached for Rapex<120 nm when the guide is excited by IR light, it is clear that only the fundamental HE11 mode can be excited at the apex and be guided through the nanoparticle of Fig. 3 [22,23]. This obviously further reduces the power actually fed into the particle via the apex, in particular far from the incidence parallel to the axis, measured to be optimal for coupling to HE11 [24]. Similar conclusions can be reached for green light although the 35-nm apex radius is closer to the first resonance peak of the spherical silicon particle (Rapex=59 nm) and to the single-mode condition (Rapex<50 nm). It can thus be expected that more electromagnetic power can be fed into the cone through its apex in case of green illumination. In summary, the light diffracted at the apex can be ignored because the coupling efficiency and guiding ability of a normally illuminated cone with small apex radius are low. It is important to note that the situation would be different for light incident (quasi-) parallelly to the axis, where the coupling through the ending/apex is indeed significant [20,24]. This coupling is e.g. used in recent successful applications where arrays of cylindrical and conical Mie resonators serve as antireflection coating [25,26].

To conclude, we have developed a simplified light-absorption theory for conical nanoparticles based on the assumption that an axisymmetric particle is a stack of cylinders with varying radius. This theory predicts that light couples and is therefore absorbed in these objects at specific locations, where the local radius corresponds to a resonance of the Mie theory developed for cylindrical particles. We have shown theoretically that, for this approximation to be good, the particle surface should be sufficiently smooth and the wavelength sufficiently short. Further, we have shown that the finiteness of the particle, i.e. the presence of an apex, has little impact on its absorption efficiency provided the radius-to-wavelength ratio is small and the incidence is close to normal to the cone axis. An experimental confirmation of the proposed theory has been found in the observation of sharply localized holes in nanoscale conical silicon particles illuminated with high-fluence IR light. The location of these holes is in excellent agreement with the developed theory.

References and links

1. M. I. Mishchenko, G. Videen, V. A. Babenko, N. G. Khlebtsov, and T. Wriedt, “T-matrix theory of electromagnetic scattering by particles and its applications: a comprehensive reference database,” J. Quant. Spectrosc. Radiat. Transf. 88, 357–406 (2004).

2. M. I. Mishchenko, N. T. Zakharova, G. Videen, N. G. Khlebtsov, and T. Wriedt, “Comprehensive T-matrix reference database: a 2007-2009 update,” J. Quant. Spectrosc. Radiat. Transf. 111(4), 650–658 (2010). [CrossRef]  

3. L. Novotny, R. X. Bian, and X. S. Xie, “Theory of nanometric optical tweezers,” Phys. Rev. Lett. 79(4), 645–648 (1997). [CrossRef]  

4. J. Zhu, Z. Yu, G. F. Burkhard, C. M. Hsu, S. T. Connor, Y. Xu, Q. Wang, M. McGehee, S. Fan, and Y. Cui, “Optical absorption enhancement in amorphous silicon nanowire and nanocone arrays,” Nano Lett. 9(1), 279–282 (2009). [CrossRef]   [PubMed]  

5. M. Idemen, “Confluent tip singularity of the electromagnetic field at the apex of a material cone,” Wave Motion 38(3), 251–277 (2003). [CrossRef]  

6. M. I. Mishchenko, L. D. Travis, and D. W. Mackowski, “T-matrix computations of light scattering by nonspherical particles: A review,” J. Quant. Spectrosc. Radiat. Transf. 55(5), 535–575 (1996). [CrossRef]  

7. L. Novotny and B. Hecht, Principles of Nano-Optics (Cambridge, 2006).

8. A. Vella, B. Deconihout, L. Marrucci, and E. Santamato, “Femtosecond field ion emission by surface optical rectification,” Phys. Rev. Lett. 99(4), 046103 (2007). [CrossRef]   [PubMed]  

9. A. Vella, N. Sevelin-Radiguet, J. Houard, and B. Deconihout, “Polarization dependence of the optical absorption of a subwavelength tip,” Appl. Surf. Sci. 258(23), 9202–9207 (2012). [CrossRef]  

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

11. P. W. Barber and S. C. Hill, Light Scattering by Particles (World Scientific, 1990).

12. R. Menzel, Photonics, 2nd ed. (Springer, 2007).

13. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Courier Dover Publications, 1965)

14. T. Yasuda and D. E. Aspnes, “Optical-standard surfaces of single-crystal silicon for calibrating ellipsometers and reflectometers,” Appl. Opt. 33(31), 7435–7438 (1994). [CrossRef]   [PubMed]  

15. C. F. Bohren, “Scattering by a sphere and reflection by a slab: some notable similarities,” Appl. Opt. 27(2), 205–206 (1988). [CrossRef]   [PubMed]  

16. K. G. Ayappa, H. T. Davis, S. A. Barringer, and E. A. Davis, “Resonant microwave power absorption in slabs and cylinders,” AIChE J. 43(3), 615–624 (1997). [CrossRef]  

17. B. Gault, F. Vurpillot, A. Vella, M. Gilbert, A. Menand, D. Blavette, and B. Deconihout, “Design of a femtosecond laser assisted tomographic atom probe,” Rev. Sci. Instrum. 77(4), 043705 (2006). [CrossRef]  

18. A. D. Bristow, N. Rotenberg, and H. M. Van Driel, “Two-photon absorption and Kerr coefficients of silicon for 850-2200 nm,” Appl. Phys. Lett. 90(19), 191104 (2007). [CrossRef]  

19. B. R. Tull, “Femtosecond laser ablation of silicon: nanoparticles, doping and photovoltaics,” PhD thesis, Harvard University (2007)

20. A. R. Jones, “Some calculations on the internal electromagnetic fluxes of cylinders,” J. Opt. 14(6), 065702 (2012). [CrossRef]  

21. H. A. Bethe, “Theory of diffraction by small holes,” Phys. Rev. 66(7-8), 163–182 (1944). [CrossRef]  

22. A. W. Snyder and J. D. Love, Optical waveguide theory (Springer, 1983)

23. L. Tong, J. Lou, and E. Mazur, “Single-mode guiding properties of subwavelength-diameter silica and silicon wire waveguides,” Opt. Express 12(6), 1025–1035 (2004). [CrossRef]   [PubMed]  

24. A. W. Snyder, C. Pask, and D. J. Mitchell, “Light-acceptance property of an optical fiber,” J. Opt. Soc. Am. 63(1), 59–64 (1973). [CrossRef]   [PubMed]  

25. P. Spinelli, M. A. Verschuuren, and A. Polman, “Broadband omnidirectional antireflection coating based on subwavelength surface Mie resonators,” Nat Commun 3, 692 (2012). [CrossRef]   [PubMed]  

26. Y. Qiu, H.-C. Hao, J. Zhou, and M. Lu, “A close to unity and all-solar-spectrum absorption by ion-sputtering induced Si nanocone arrays,” Opt. Express 20(20), 22087–22094 (2012). [CrossRef]   [PubMed]  

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 (3)

Fig. 1
Fig. 1 Resonance behavior of the electric field in a cylindrical silicon particle of radius R , illuminated with a TM wave of amplitude E 0 and wavelengths λ 0 =515 nm (full line) and λ 0 =1030 nm (dotted line), incident normally to the cylinder axis. The shown values are on the illuminated side ( r=R,ϑ=π ).
Fig. 2
Fig. 2 Variation in the right-hand side of Eq. (3) as a function of the local radius R(z) in a conical silicon particle with cone angle α=10° illuminated respectively by green (full line) and IR (dotted line) TM waves of amplitude E 0 , incident normally to the cone axis. The shown values are on the illuminated side ( r=R,ϑ=π ).
Fig. 3
Fig. 3 (left) Theoretical value of | E z,int [r=R(z),ϑ=π,z] | 4 along the illuminated side of an IR-irradiated conical silicon tip as a function of the local radius R(z) . (right) Silicon tip irradiated with a high-fluence (~100 mJ/cm2) IR pulsed laser. Resonance absorption occurs in the cross-sections with R(z) ~113, 178 and 241 nm. The tip is tilted 45° with respect to the primary beam of the SEM.

Equations (3)

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

E z,int (r,ϑ)= E 0 [ d 0 ( R λ 0 ) J 0 ( 2π n ˜ λ 0 r) + k=1 2 i k d k ( R λ 0 ) J k ( 2π n ˜ λ 0 r)cos(kϑ) ] ,
R(z)=ztan(α/2),
| E z,int |>> | λ 0 tan(α/2) 2π n ˜ | 2 | 2 E z,int R 2 | .
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.