Broadband converging plasmon resonance at a conical nanotip

We propose an analytical theory which predicts that Converging Plasmon Resonance (CPR) at conical nanotips exhibits a red-shifted and continuous band of resonant frequencies and suggests potential application of conical nanotips in various fields, such as plasmonic solar cells, photothermal therapy, tip-enhanced Raman and other spectroscopies. The CPR modes exhibit superior confinement and ten times broader scattering bandwidth over the entire solar spectrum than smooth nano-structures. The theory also explicitly connects the optimal angles and resonant optical frequencies to the material permittivities, with a specific optimum half angle that depends only on the real permittivity for high-permittivity and low-loss materials. © 2013 Optical Society of America OCIS codes: (240.6680) Surface plasmons; (310.6628) Subwavelength structures, nanostructures; (350.6050) Solar energy; (350.5340) Photothermal effects. References and links 1. V. Ferry, J. Munday, and H. Atwater, “Design considerations for plasmonic photovoltaics,” Adv. Mater. 22, 4794–4808 (2010). 2. X. Lu, M. Rycenga, S. Skrabalak, B. Wiley, and Y. Xia, “Chemical synthesis of novel plasmonic nanoparticles,” Annu. Rev. Phys. Chem. 60, 167–192 (2009). 3. C. Lee, S. Bae, S. Mobasser, and H. Manohara, “A novel silicon nanotips antireflection surface for the micro sun sensor,” Nano Lett. 5, 2438–2442 (2005). 4. A. R. Parker and H. E. Townley, “Biomimetics of photonic nanostructures,” Nat. Nanotechnol. 2, 347–353 (2007). 5. A. V. Goncharenko, H. C. Chang, and J. K. Wang, “Electric near-field enhancing properties of a finite-size metal conical nano-tip,” Ultramicroscopy 107, 151–157 (2007). 6. X. Huang, I. H. El-Sayed, W. Qian, and M. A. El-Sayed, “Cancer cell imaging and photothermal therapy in the near-infrared region by using gold nanorods,” J. Am. Chem. Soc 128, 2115–2120 (2006). 7. R. Rodriguez-Oliveros and J. A. Sanchez-Gil, “Gold nanostars as thermoplasmonic nanoparticles for optical heating,” Opt. Express 20, 621–626 (2012). 8. A. Mohammadi, F. Kaminski, V. Sandoghdar, and M. Agio, “Fluorescence enhancement with the optical (bi-) conical antenna,” J. Phys. Chem. C 114, 7372–7377 (2010). 9. H. Shen, N. Guillot, J. Rouxel, M. de la Chapelle, and T. Toury, “Optimized plasmonic nanostructures for improved sensing activities,” Opt. Express 20, 21278–21290 (2012). 10. A. F. Stevenson, “Solution of electromagnetic scattering problems as power series in the ratio (dimension of scatter)/wavelength,” Appl. Phys. Lett. 24, 1134–1141 (1953). 11. M. I. Stockman, “Nanofocusing of optical energy in tapered plasmonic waveguides,” Phys. Rev. Lett. 93, 137404 (2004). 12. A. Goncharenko, J. K. Wang, and Y. C. Chang, “Electric near-field enhancement of a sharp semi-infinite conical probe: Material and cone angle dependence,” Phys.Rev.B 74, 235442 (2006). 13. J. M. Pitarke, V.M. Silkin, E. V. Chulkov, and P. M. Echenique, “Theory of surface plasmons and surface-plasmon polaritons,” Rep. Prog. Phys 70, 1–87 (2007). !"#$%&'()(*"+,--(./


Introduction
Surface plasmon resonance on metallic nanostructures has been extensively studied recently in their potential application in solar cell, photothermal therapy, etc. Conical nanotips have been shown to achieve better performance in the above mentioned fields, compared to smooth nanostructures like nanospheres.In solar cell application, metallic nanospheres are commonly used to enhance light absorption [1].However, the narrow resonance spectrum of metallic nanospheres, with a limited bandwidth of about 100 nm [2], restricts the frequency range of light that can be absorbed from the entire solar spectrum (400 nm to 1200 nm).Conical nanotips have realized broadband absorption in experimental and numerical efforts [3][4][5], thus improving solar cell efficiency.Although broad bandwidth adsorption at the gold tip does not contribute to current generation directly, the plasmonics are scattered from the tip into the photovoltaic material.Metallic nanostructures, due to the intense and localized heat generated at plasmon resonant frequency, are suitable for cancer photothermal therapy.And it is desirable to shift resonant frequency from visible range to near infrared (NIR) range to reduce light absorbed by intrinsic chromophores in native tissue [6].Simulation and experimental results have shown that nanotips can excite red-shifted plasmon resonance without compromising the intensity compared to gold nanospheres [7][8][9].The same red-shifted broad bandwidth plasmonic excitation also suggests that multi-spectral Raman and other optical near field spectroscopy techniques can be enhanced by a conical tip.
However, a closed-form analytical theory for broadband resonance at conical nanotips that is able to quantify the bandwidth and capture the underlying mechanisms is still lacking.Existing theories failed to capture the intrinsic features of CPR due to their invalid assumptions when dealing with singular geometries.Quasi-static (long wavelength perturbation) theory [10] is frequently used to determine the plasmonic frequency of nanostructures in the literature, particularly for nanospheres, but this theory does not apply for sharp geometries with significant wavelength compression.Slender body WKB theory [11] that models a singular cone with a cylinder with a slowly varying diameter does not describe the crucial EM azimuthal field interaction around the tip of the cone and fails to capture its broadband resonance spectrum.Although a closed-form dispersion relation was given in [12], further analysis is lacking to explicitly relate broadband resonance and red-shifted resonance wavelength to material dielectric properties.
This Letter reports a rigorous analysis of the CPR resonance spectrum on conical nanotips.Our theory unveils the singular behavior of CPR on conical nanotips in the near field, captures its unique broadband resonance property and explicitly relates the resonance frequency and optimum angle to the material dielectric property.With our theory, an optimized cone is shown to have resonance frequency red-shifted and scattering spectra that span the entire solar spectrum, ten times that of a nanosphere.

Derivation of dispersion relationship
We solve the full Maxwell equations for localized EM field at the tip of an infinite long cone in spherical coordinate, where the origin of the coordinate resides at cone tip, as shown in the inset of Fig. 1(a).The apex angle of the cone is 2α and the cone is made of a material with a complex permittivity ε m and a permeability µ m surrounded by a dielectric medium with a real permittivity ε o and a permeability µ o , as depicted in the inset of Fig. 1(a).For the TM (transverse magnetic) mode with a magnetic field in the angular direction, the non-zero components of the EM fields are H φ , E r , E θ , and the full set of Maxwell equations are simplified to a single scalar equation ∇ 2 H φ − H φ /(r 2 sin 2 θ )+β 2 j H φ = 0, where β 2 j = ω 2 µ j ε j ( j = m, o), ω is the frequency of the incoming light, which is the same for the CPR resonance.The radially confined harmonics within the cone-like tips (0 (βr) is the Bessel function of the first kind and P 1 υ (cos(θ )) is the associate Legendre polynomial of the first kind.The corresponding harmonics for the dielectric Note that two types of Legendre polynomials are chosen for the two media to satisfy the symmetry conditions.The θ -dependence of the EM field suggests that the EM field does not decay away from the interface, in contrast to the logarithm decay in WKB theory [11].Applying the continuity of the field and displacement on the interface ) and omitting higher order less singular plasmonic modes, a transcendental CPR dispersion relationship is obtained for the dominant eigenvalue υ, This dispersion relationship, derived using different approaches, is identical to Eq. ( 6) in [12].
But new underlying mechanism are to be unveiled by conformal map techniques and expansion near singular point, which will be extensively discussed in this letter.
As expected, the symmetry of Eq. ( 1) at α = π/2 yields the classical planar plasmonic resonance condition ε m /ε o + 1 = 0 and the frequency dependence of the complex permittivity ε m defines the planar plasmon resonant frequency ω s with this planar resonance condition [13].For TE (transverse electric) mode, the permittivity ratio in Eq. ( 1) is replaced by the permeability ratio (This ratio is close to unity for most materials in the optical range and since Re[g(α, υ)] grows rapidly as α decreases below π/2, TE resonance can never be realized [14]).In near field, the radial component of the electric field E r scales as r υ−1 (r is the radial distance from the apex of the cone), while the azimuthal magnetic field H φ is less singular since it scales as r υ (υ = υ r + iυ i ).The restriction that the total EM energy be finite over the entire domain requires that − 1 2 < υ r < 1 for cones [15].The position-dependent wave-number κ = υ i /r increases dra- matically as CPR waves propagate toward the tip, which is found to be negative (υ i < 0) for the converging waves with speed retardation.

Conformal mapping of dispersion relationship
Equation (1) defines a conformal map −g(α, υ) from the complex υ space to the complex permittivity ε m /ε o (ε m = ε r + iε i ) space for cones (The dielectric permittivity ε o of the outside medium will be set to unity for simplicity).The CPR sustainable region (grey area in Fig. 1(a)) is confined within the bold curve (the image of υ r = 1) and the horizontal axis.The point 2 )) on the real line of the complex permittivity plane, which is the Grey regions shows CPR sustainable region for a cone with half angle α = π/6 and the contour lines correspond to constant υ r .The branch point dg(α, υ)/dυ = 0 of the conformal map is indicated as a solid dot and the limit point (υ i → ∞)i si n d i c a t e da sa no p e nc i r c l eo nt h en e g a t i v eε r axis.The bold circular contour is the image of υ r = 1.The hollow circles are empirical permittivities for gold ranging from infrared to ultraviolet wavelength [16] and the dashed curve is the Drude permittivity model for gold [13].The band of excited frequencies for the Drude model fall at the "resonance" band between the branch point and the limit point on the real axis.(b) Intensification exponents for the Drude gold cone with half angle 5, 10, 15 and 20 degrees are plotted as solid curves.The dotted line marks the resonance peak position of a single nanosphere for comparison [2].The broad bandwidth extends beyond the visible range for the Drude cone and the resonance frequency is red-shifted compared to gold nanospheres.image of υ →− 1 2 + 0i, corresponding to the most singular CPR mode with vanishing wave number κ, is a branch point of the conformal map (dg(α, υ)/dυ = 0, marked as a solid circle in Fig. 1(a)).Another limit point (ε r = −1, ε i = 0,marked as an open circle in Fig. 1(a)) of the conformal map corresponds to the common large υ r limit for the pure dielectric (the image of the most singular mode with infinite imaginary part υ →− 1 2 + i∞).Maximum EM field intensification occurs at the real permittivity axis from the branch point to the limit point (ε ro < ε r < −1).The gold permittivity data of classical Drude free electron theory [13] is plotted as a dashed line in Fig. 1(a), whose imaginary permittivity decreases monotonically with frequency.At high frequency where ε i ∼ 0 and when maximum possible intensification occurs, CPR of Drude conical nanotips hence exhibits a continuous band of resonant frequencies, whose corresponding real permittivies vary from ε ro to -1.This continuous broad band of resonant frequency is unique to conical nanotips due to their continuous change of length scale, while finite length scale of smooth nanoparticles only exhibits discrete and sharp resonance frequency.The intensification exponents 1 − υ r computed from Eq. ( 1) for a Drude gold cone are shown as solid curves in Fig. 1(b) for different angles and wavelengths λ of incident light.The entire solar spectrum is excited at small angles.In contrast, the extinction spectrum for a single nanosphere [2] has a bandwidth that is 1/10 of the broadband CPR spectrum for conical nanotips.Empirical gold permittivity data which includes conduction loss due to electron interband transition [16] are also inserted as open circles in Fig. 1(a).The imaginary (loss) permittivity of gold exhibits a minimum with respect to frequency which suggests a global resonant frequency for these metals instead of the continuous band of resonant frequen-wavelength (nm) half angle (degree) wavelength (nm) (a) wavelength (nm) Fig. 2. Constant intensification exponent contours lines (1 − υ r = 1.428, 1.4, 1.3, 1.1) for a gold cone for empirical permittivities from [16], with different angles and incident light wavelengths, where the dashed curve connects the local resonant conditions from large angles until the global optimal angle (circle).The inset (a) plots the intensification exponent for a gold cone with different cone angles (α = 5 o , 10 o , 15 o , 20 o , 30 o ) and shows the CPR spectrum broadens and exhibits a red shift away from the planar plasmon resonant wavelength λ s = c/ω s .The inset (b) plots the ratio of the imaginary to real permittivity for gold as a function of wavelength [16], where maximum intensification occurs at λ = 892nm.
cies for Drude-type metals.Nevertheless, the scattering feature of the cone still suggests a large bandwidth of excited/scattered light around the global resonant frequency.

Global resonance frequency
As field focusing at small cone angles can increase both the field intensity and also the conductive loss, an optimal angle is expected.Conversely, the ratio −ε i /ε r of the material determines the energy loss-storage ratio of the system [17], thus should determine the global resonant frequency discussed above.Contour lines of constant EM field intensification exponent 1 − υ r computed from Eq. ( 1) for a gold cone with the permittivity model of [16] are shown in Fig. 2 for different angles and wavelengths λ of incident light.A global maximum in the intensification exponent exists at λ = 892nm and α = 11 o .As expected, this global resonant frequency resides at the minimum energy loss-storage ratio of gold, as shown in Fig. 2(b).The dash line connects local resonant frequency for different angles which ends at the global optimal angle.The resonance spectrum for a gold cone at different cone angles is shown in Fig. 2(a).

Branch point
Asymptotic behavior of g (α, υ) in the limit of low conduction loss (ε i /ε r → 0) is most relevant for both Drude-type and other metals at the IR and visible frequency range.The five curves in the complex υ plane shown in the inset of Fig. 3 are the loci of the dominant eigenvalues as a function of the half-angle α for a cone with ε r = −182 and for five different imaginary permittivities ε i = 0.01, 0.5, 1, 3 and 10.Optimal angles for each ε i , marked as solid circles on the root loci, follows a horizontal line, which indicates that the optimal angles are only a function of ε r for metal with low conduction loss.This invariance to conduction loss is due to a curious asymptotic behavior of the dominant eigenvalues in the no-loss dielectric limit.Their smooth loci in Fig. 3 approach the real line on left half plane (υ r < 0, υ i = 0) and a vertical line (υ r = − 1 2 and υ i < 0) at the pure dielectric limit of ε i → 0, with the intersection of the two orthogonal lines occurring at α = α * .This intersection of the two lines at υ = −1/2 and α = α * is the image of the branch point (ε r , ε i )=(ε ro , 0).Given the wavelength of incident light and the material property, the branch point angle α = α * can be obtained from dg/dα(α * , −1/2)= 0. A leading order asymptotic estimate of the branch-point angle α * for pure dielectric with |ε r | 1, leads to the following explicit relation However, α * is not the optimum half angle for the zero-loss dielectric.The above asymptotic limit, shown as a bold line in Fig. 3, is actually the upper bound of a continuous range of optimum angles for metal with zero loss-all angles smaller than α * can sustain plasmon waves with the maximum amplitude -1/2 for the eigenvalue υ r .Also, the resonance bandwidth (related to υ r ) can be estimated, given this range of half cone angles.The change of optimal angle from a discrete value for finite ε i to a continuous band for ε i = 0 is reflected by the contortion of the smooth loci to a kink singularity connecting two straight lines in Fig. 3. Specific and discrete optimum angles for metal with low but finite loss (ε i > 0) are depicted as solid circles in the inset and they can be better estimated by another critical point υ o = −1/2 + υ io i at α = α o with a constant asymptote for the imaginary exponent υ i < 0.

Critical point
The optimum angle at a particular frequency can be obtained by solving Eq. ( 1) and optimal condition dυ r /dα = 0 simultaneously at the limit of ε i /ε r → 0 and ε r /ε o →−∞.Restricting ourselves to regions close to the critical point, an explicit Taylor expansion about the critical point can be carried out to yield υ r + 1/2 = −iε i /g υ (α o , υ o ).The optimal condition then re- duces to g o α − g o υ g o υα /g o υυ = 0. Combination of this criterion and Eq.(1) (ε i = 0) gives rise to the optimal angle α o and the corresponding υ io for any given ε r .Good agreement between the optimal angles predicted by Eq. (1) (solid curve in Fig. 3), critical local resonance points (dotted curve in Fig. 3) and numerical data from [18] is observed in a large range of optical excitation frequency that yields a range of ε r between -200 and -2.Independent of ε i , both the branch point angles and critical point angles approach zero at the large ε r limit.
It can be shown that, for finite loss and at the limit of ε r /ε o →−∞, iε r /g υ (α o , υ o ) approaches a constant and hence υ r + 1/2 ∼−0.693(ε i /ε r ).This estimate of the dominant exponent at optimal angles for each frequency confirms that the global resonant frequency occurs at the minimum of the loss-storage ratio in Fig. 2(b).It also allows a convenient estimate of the local and global resonant frequencies from the material properties (ε i /ε r )(ω).As seen in Fig. 3, both the predicted exponent and optimum angle at a particular real permittivity are in good agreement with numerical values from Eq. (1).

Comparison with numerical results
Intensification exponent calculated from Eq. ( 1) can be compared to simulation data published in the literature [18].Knowing |E r | ∼ r υ r −1 , we obtain the dependence of the ratio of the electric field intensity on the exponents and the geometry E tip /E o =(r 1 /r 2 ) υ r −1 , where r 1 is the radial distance from the apex of the cone to the point where E tip is measured and r 2 is the radial distance from the apex to the point where laser with electric field E o impinges on the cone.
Branch point angle approximated by Eq. ( 2) Optimal angle predicted by Eq. ( 1) Critical point angle Optimal angle from literature optimal half angle (degree)

Critical point
Fig. 3. Optimum angles of a gold cone at different ε r corresponding to different excitation wavelength.Simulation data from the solution of the Maxwell equation [18] are compared to estimates from the branch point, the critical point and the exact solution of the CPR dispersion relationship Eq. ( 1).The inset depicts the dominant CPR intensification exponent as a function of the half angle for five imaginary permittivities ε i = 0.01, 0.5, 1, 3 and 10.0 at ε r = −182.The optimum angle yielding the most negative υ r is computed numerically for each ε i and is marked as a solid circle on each root locus.Analytical estimate, with a ε i /ε r scaling, from an expansion about the critical point is shown as open triangles.At the limit of ε i /ε r → 0, the eigenvalue of this optimum angle approaches the critical point at υ o = −1/2 + υ io i.
From simple geometry, r 1 = R 1 sin α and r 2 = R 2 sin α where R 1 =5 nm is the radius of curvature of the apex and R 2 =300 nm is the radius of the cylindrical wire at the far end.Figure 4 favorably compares the literature data log 10 (E tip /E o ) to our analytical estimate (υ r − 1) log 10 (R 1 /R 2 ), where the exponents are evaluated from Eq. ( 1).[18] (scattered points) for gold cones of different angles (R1=5 nm, R2=300 nm).(b) Comparison against literature data [18] for three different materials at excitation wavelength 488 nm with the same R1 and R2.

Conclusion
In conclusion, we have provided an analytical theory that explains the underlying mechanism of red-shifted resonance wavelength and broadband CPR spectrum at conical nanotips.Our theory predicts a continuous band of resonant frequencies for Drude-type metals and a global resonant frequency for metals with electron interband transition.Asymptotic local analysis of the CPR dispersion relation offers simple estimates of the optimal angles and resonant frequencies for materials with low loss.Intensification exponents calculated from the CPR dispersion relation are favorably compared with published simulation results.Our theory suggests that conical nanotips are ideal scattering centers for enhancing light absorption in solar cells, due to its broadband resonance and strong EM field intensity.Its red-shifted resonate bandwidth are also desirable for cancer photothermal therapy, multi-spectral Raman and other near-field spectroscopy techniques.

Fig. 1 .
Fig.1.(a) Conformal map in complex permittivities domain.Grey regions shows CPR sustainable region for a cone with half angle α = π/6 and the contour lines correspond to constant υ r .The branch point dg(α, υ)/dυ = 0 of the conformal map is indicated as a solid dot and the limit point (υ i → ∞)i si n d i c a t e da sa no p e nc i r c l eo nt h en e g a t i v eε r axis.The bold circular contour is the image of υ r = 1.The hollow circles are empirical permittivities for gold ranging from infrared to ultraviolet wavelength[16] and the dashed curve is the Drude permittivity model for gold[13].The band of excited frequencies for the Drude model fall at the "resonance" band between the branch point and the limit point on the real axis.(b) Intensification exponents for the Drude gold cone with half angle 5, 10, 15 and 20 degrees are plotted as solid curves.The dotted line marks the resonance peak position of a single nanosphere for comparison[2].The broad bandwidth extends beyond the visible range for the Drude cone and the resonance frequency is red-shifted compared to gold nanospheres.