Optimization of broadband optical response of multilayer nanospheres

doi:10.1364/OE.20.018494

« Show journal navigation

Optimization of broadband optical response of multilayer nanospheres

Wenjun Qiu, Brendan G. DeLacy, Steven G. Johnson, John D. Joannopoulos, and Marin Soljačić »View Author Affiliations

Optics Express, Vol. 20, Issue 16, pp. 18494-18504 (2012)
http://dx.doi.org/10.1364/OE.20.018494

View Full Text Article

Acrobat PDF (1062 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (26)
Cited By
Figures (4)
Metrics

Abstract

We propose an optimization-based theoretical approach to tailor the optical response of silver/silica multilayer nanospheres over the visible spectrum. We show that the structure that provides the largest cross-section per volume/mass, averaged over a wide frequency range, is the silver coated silica sphere. We also show how properly chosen mixture of several species of different nanospheres can have an even larger minimal cross-section per volume/mass over the entire visible spectrum.

1. Introduction

Nanoparticles with strong optical response, characterized by scattering, absorption and total cross-sections, have wide applications in biomedical imaging, photothermal therapy, and optical obscurance [1

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]

–3

C. Noguez, “Surface plasmons on metal nanoparticles: the influence of shape and physical environment,” J. Phys. Chem. C 111, 3806–3819 (2007). [CrossRef]

]. Different applications require different optical response properties. For instance, real-time biomedical imaging is based on large scattering cross-sections, while photothermal therapy requires nanoparticles with large absorption cross-sections and small scattering cross-sections. For obscurance applications [4

S. Alyones, C. Bruce, and A. Buin, “Numerical methods for solving the problem of electromagnetic scattering by a thin finite conducting wire,” IEEE Trans. Antennas Propag. 55, 1856–1861 (2007). [CrossRef]

, 5

C. W. Bruce and S. Alyones, “Extinction efficiencies for metallic fibers in the infrared,” Appl. Opt. 48, 5095–5098 (2009). [CrossRef] [PubMed]

], the ideal nanoparticle should typically have large total cross-sections over the whole visible spectrum while keeping the volume or mass of nanoparticles as small as possible. The diversity and complexity of these requirements necessitate an engineering approach of nanoparticle design.

Previous studies on the optical response of nanoparticles are mainly based on parametric approach [6

P. K. Jain, K. S. Lee, I. H. El-Sayed, and M. A. El-Sayed, “Calculated absorption and scattering properties of gold nanoparticles of different size, shape, and composition: Applications in biological imaging and biomedicine,” J. Phys. Chem. B 110, 7238–7248 (2006). [CrossRef] [PubMed]

–8

C. E. Román-Velázquez and C. Noguez, “Designing the plasmonic response of shell nanoparticles: Spectral representation,” J. Chem. Phys. 134, 044116 (2011). [CrossRef] [PubMed]

], which works well for simple structures. However, when the structure becomes complicated and the number of design parameters increases, optimization becomes the preferable approach because it can efficiently explore the whole parameter space. Furthermore, because the optimization objective function can be an arbitrary transformation of the frequency-dependent cross-sections, this approach is very powerful in tailoring the broadband optical response of nanoparticles. For example, some applications may require that the optical resonance has both a strong peak value and a wide bandwidth. This can be achieved by maximizing the average cross-section over the bandwidth of interest. For obscurance applications, we want the total cross-section to be consistently large over the whole visible spectrum. This is equivalent to maximizing the minimal cross-section over this spectrum. In this paper, we will use an optimization tool to tailor the optical response of multilayer nanospheres over wide frequency range of interest.

Before we start, we need to select the material system. Nanoparticles composed of metal and dielectric materials support surface plasmons on the metal/dielectric interfaces and can strongly interact with light in the visible range [9

S. Oldenburg, R. Averitt, S. Westcott, and N. Halas, “Nanoengineering of optical resonances,” Chem. Phys. Lett. 288, 243 – 247 (1998). [CrossRef]

–13

R. Bardhan, S. Mukherjee, N. A. Mirin, S. D. Levit, P. Nordlander, and N. J. Halas, “Nanosphere-in-a-nanoshell: A simple nanomatryushka,” J. Phys. Chem. C 114, 7378–7383 (2010). [CrossRef]

]. At resonance, the cross-sections of these nanoparticles are much larger than their physical cross-sections, which makes them super-scatterers and super-absorbers [14

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

–17

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

]. Furthermore, the plasmon resonance frequency can be tuned by varying the physical structure of the nanoparticles. In order to be able to tailor the optical response of nanoparticles over a wide frequency range, we choose the metal/dielectric material system. For concreteness, we will focus on the silver/silica material system.

In this article, we first formulate the Transfer Matrix Method to calculate the optical response of multilayer nanospheres. After briefly reviewing the optical response of bilayer silver/silica nanospheres, we proceed to optimize the average cross-section of various silver/silica multilayer nanospheres. Our optimization results show that the structure with the maximal average cross-section is the bilayer silver/silica structure with silver as the shell. Finally, we investigate using a mixture of several species of bilayer nanospheres to enhance the minimal cross-sections over the entire visible range.

2. Calculation of optical response via transfer matrix method

We start by formulating the Transfer Matrix Method for multilayer nanosphere [18

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

, 19

C. Bohren and D. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Songs, 1983).

]. Consider a multilayer nanosphere shown in Fig. 1. Because of the spherical symmetry, the fields at a given incident frequency can be decomposed into two orthogonal polarizations: transverse electric (TE) and transverse magnetic (TM). For TE polarization, the electric fields can be written as E_TE = ∇ × rϕ_TE. For TM polarization, the magnetic fields can be written as H_TM = ∇ × rϕ_TM. The scalar potential ϕ_TE and ϕ_TM satisfy the scalar Helmholtz equation ∇²ϕ + k²ϕ = 0 where k² = ω²ε (r). Due to the spherical symmetry, ϕ can be decomposed into a discrete set of spherical modes:

ϕ_{l m} = R_{l} (r) P_{l}^{| m |} (cos θ) exp (i m φ)

with l = 0, 1, 2,... and m = −l,...,l. Since ε(r) is a constant ε_i inside the i^th shell, R_l(r) is a linear combination of the first and second kind spherical Bessel functions within the individual shells:

{R_{l} (r) |}_{i} = A_{i} j_{l} (k_{i} r) + B_{i} y_{l} (k_{i} r)

(1)

The coefficients (A_i, B_i) of adjacent shells are linked by the transfer matrix of the interface:

[\begin{matrix} A_{i + 1} \\ B_{i + 1} \end{matrix}] = M_{i + 1, i} [\begin{matrix} A_{i} \\ B_{i} \end{matrix}]

(2)

The matrix element is determined by the boundary condition satisfied by R_l(r), which comes from the continuity of the tangent components of E and H across the boundary. For TE polarization, rR_l(r) and (rR_l(r))′ are continuous across the boundary. By writing the continuity conditions in matrix form, we get:

\begin{array}{l} M_{i + 1, i} & = & {[\begin{matrix} j_{l} (k_{i + 1} R_{i}) & y_{l} (k_{i + 1} R_{i}) \\ j_{l}^{'} (k_{i + 1} R_{i}) k_{i + 1} R_{i} + j_{l} (k_{i + 1} R_{i}) & y_{l}^{'} (k_{i + 1} R_{i}) k_{i + 1} R_{i} + y_{l} (k_{i + 1} R_{i}) \end{matrix}]}^{- 1} \times \\ [\begin{matrix} j_{l} (k_{i} R_{i}) & y_{l} (k_{i} R_{i}) \\ j_{l}^{'} (k_{i} R_{i}) k_{i} R_{i} + j_{l} (k_{i} R_{i}) & y_{l}^{'} (k_{i} R_{i}) k_{i} R_{i} + y_{l} (k_{i} R_{i}) \end{matrix}] \end{array}

(3)

For TM polarization, rR_l(r) and (ε⁻¹rR_l(r))′ are continuous across the boundary. By writing the continuity conditions in matrix form, we get:

\begin{array}{l} M_{i + 1, i} & = & {[\begin{matrix} j_{l} (k_{i + 1} R_{i}) & y_{l} (k_{i + 1} R_{i}) \\ j_{l}^{'} (k_{i + 1} R_{i}) k_{i + 1} R_{i} + j_{l} (k_{i + 1} R_{i}) & y_{l}^{'} (k_{i + 1} R_{i}) k_{i + 1} R_{i} + y_{l} (k_{i + 1} R_{i}) \end{matrix}]}^{- 1} \times \\ [\begin{matrix} j_{l} (k_{i} R_{i}) & y_{l} (k_{i} R_{i}) \\ \frac{ε_{i + 1}}{ε_{i}} (j_{l}^{'} (k_{i} R_{i}) k_{i} R_{i} + j_{l} (k_{i} R_{i})) & \frac{ε_{i + 1}}{ε_{i}} (y_{l}^{'} (k_{i} R_{i}) k_{i} R_{i} + y_{l} (k_{i} R_{i})) \end{matrix}] \end{array}

(4)

The transfer matrix of the whole system can be calculated by cascading the transfer matrices of individual interfaces.

[\begin{matrix} A_{n + 1} \\ B_{n + 1} \end{matrix}] = M_{n + 1, n} M_{n, n - 1} \dots M_{3, 2} M_{2, 1} [\begin{matrix} A_{1} \\ B_{1} \end{matrix}] = M [\begin{matrix} A_{1} \\ B_{1} \end{matrix}]

(5)

Fig. 1 Schematic of an n layer nanosphere embedded in infinite dielectric medium. The outer radius and dielectric function of individual layers are (R_i, ε_i), i = 1, 2,..., n. The dielectric function of the medium is ε_m. The solid lines represent an incident plane wave which contains incoming and outgoing waves. The dashed line represents the scattered wave which only contains outgoing wave.

Since the second kind of Bessel function is singular at the origin, we can set A₁ = 1 and B₁ = 0. So the coefficients of Bessel functions in the surrounding medium are directly given by the transfer matrix element, A_n₊₁ = M₁₁ and B_n₊₁ = M₂₁. Within the surrounding medium, it is convenient to write the radical function as a linear combination of the spherical Hankel functions:

{R_{l} (r) |}_{n + 1} = C_{n + 1} h_{l}^{1} (k_{n + 1} r) + D_{n + 1} h_{l}^{2} (k_{n + 1} r)

(6)

Taking the convention that the fields vary in time as e^−iωt,

h_{l}^{1} (k_{n + 1} r)

and

h_{l}^{2} (k_{n + 1} r)

correspond to outgoing and incoming waves respectively. The reflection coefficient of the whole system is given by

r_{l} = \frac{C_{n + 1}}{D_{n + 1}} = \frac{M_{11} - i M_{21}}{M_{11} + i M_{21}}

(7)

The reflection coefficient as a function of frequency determines the optical response of the nanoparticle under all possible illumination conditions. Specifically, when the nanoparticle is illuminated by a linearly polarized plane wave, the incident field can be decomposed into both TE and TM channels (l, m) with l = 1, 2,... and m = −1, 1. For each channel, the incident field contains both incoming and outgoing waves carrying the same power [14

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

–17

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

]

P_{l, m = \pm 1} = \frac{λ^{2}}{16 π} (2 l + 1) I_{0}

(8)

where I₀ is the incident intensity, and λ is the wavelength in the surrounding medium. The scattered field is a purely outgoing wave characterized by the scattering coefficient S_l which is related to the reflection coefficient through S_l = (r_l − 1)/2. The scattered and absorbed power in this channel is given by:

P_{l, m = \pm 1}^{sca} = \frac{λ^{2}}{16 π} (2 l + 1) I_{0} {| 1 - r_{l} |}^{2}

(9)

P_{l, m = \pm 1}^{abs} = \frac{λ^{2}}{16 π} (2 l + 1) I_{0} (1 - {| r_{l} |}^{2})

(10)

By summing over the contribution from all channels of TE and TM polarization, we get the scattering and absorption cross-sections:

σ_{sca} = \sum_{σ} \sum_{l = 1}^{\infty} \frac{λ^{2}}{8 π} (2 l + 1) {| 1 - r_{σ, l} |}^{2}

(11)

σ_{abs} = \sum_{σ} \sum_{l = 1}^{\infty} \frac{λ^{2}}{8 π} (2 l + 1) (1 - {| r_{σ, l} |}^{2})

(12)

where σ is TE or TM. The total cross-section is the sum of the scattering and absorption cross-sections. σ_tot = σ_sca +σ_abs.

3. Optical response of silver/silica bilayer nanospheres

In this section, we analyze the optical response of silver/silica bilayer nanospheres as the building elements of multilayer structures. For silver, the complex dielectric function as a function of frequency is generated by linearly interpolating the experimental data [20

E. D. Palik, Handbook of Optical Constants of Solids (Academic Press, New York, 1985).

]. For simplicity, the size dependence of Ag’s dielectric function is not taken into account. For silica, the dielectric function is taken as a constant ε = 2.1. There are two configurations of silver/silica bilayer nanoparticle depending on the core material.

First, we consider silica coated silver spheres. Colloidal suspensions of such nanoparticles were obtained from Nanocomposix (Nanocomposix Inc., San Diego, CA). Figure 2 shows the measured and calculated total cross-sections of the fabricated nanoparticles. The calculation agrees quite well with the measurement. The total cross-section peaks around 455nm. This peak comes from the l = 1 surface plasmon mode at the silver/silica boundary. The peak wavelength only varies slightly when the inner and outer radius change. For instance, consider a silica coated silver sphere suspended in air. Fixing the outer radius at 50nm, the peak wavelength varies from 410nm to 415nm when the inner radius varies from 5nm to 45nm. Fixing the aspect ratio R₁/R₂ at 0.8, the peak wavelength varies from 390nm to 480nm when the outer radius varies from 25nm to 75nm.

Fig. 2 The total cross-section of silica coated silver spheres suspended in ethanol. The cross-section is normalized by volume (the left axis) and mass (the right axis). The insert is a TEM image of the fabricated nanoparticles. The radius of the silver core has a distribution with mean 26.3nm and standard deviation 9.3nm. The thickness of the silica shell is around 25.3nm. The red line is the measured total cross-section. The black bar represents the standard deviations from eight transmission measurement on eight samples. The blue line is the Transfer Matrix calculation of the total cross-section with the radius of the silver core sampled from the measured distribution and the thickness of the silica shell fixed at 25.3nm. The dielectric function of ethanol is taken as ε_m = 1.85.

In contrast, the surface plasmon resonance of the reverse configuration has great tunability over the visible range [9

S. Oldenburg, R. Averitt, S. Westcott, and N. Halas, “Nanoengineering of optical resonances,” Chem. Phys. Lett. 288, 243 – 247 (1998). [CrossRef]

–13

R. Bardhan, S. Mukherjee, N. A. Mirin, S. D. Levit, P. Nordlander, and N. J. Halas, “Nanosphere-in-a-nanoshell: A simple nanomatryushka,” J. Phys. Chem. C 114, 7378–7383 (2010). [CrossRef]

]. Consider a silver coated silica sphere suspended in air. Fixing the outer radius at 50nm, the peak wavelength varies from 405nm to 720nm when the inner radius varies from 5nm to 45nm. Besides the peak wavelength, the relative strength of scattering and absorption cross-sections in the total cross-section also vary. For [R₁, R₂] = [5nm, 50nm], the absorption cross-section accounts for 25% of the total cross-section at resonance. For [R₁, R₂] = [45nm, 50nm], this percentage rises to 60%. The tunability of the resonance wavelength and the tunability of the total cross-section composition makes silver coated silica sphere a good candidate for achieving broadband optical response.

4. Optimization of average cross-sections over wide frequency range

The optical response of silver/silica bilayer nanosphere indicates that increasing the number of metal/dielectric interfaces can provide additional tunability in the optical response. Using this insight, we aim to design silver/silica multilayer nanosphere with a large average cross-section over wide frequency range. The figure of merit (FOM) is the scattering, absorption, and total cross-section averaged over the target frequency range, normalized by volume or mass.

F O M = \frac{1}{ω_{\max} - ω_{\min}} \int_{ω_{\min}}^{ω_{\max}} σ_{normalized} d ω

(13)

For concreteness, we take 400–600nm and 600–800nm as the target frequency range of interest. The structure under consideration is a multilayer nanosphere with alternating silver and silica layers (six layers in total). The design parameters are the thicknesses of individual layers. The lower bounds of the thicknesses are set to be zero. The upper bound of the allowed thickness is set to be a large value (1μm). Therefore, this general structure includes structures with fewer layers (one through five layers) as boundary points.

We performed the optimization using numerical optimization package NLOPT [21

S. G. Johnson, “The nlopt nonlinear optimization package, http://ab-initio.mit.edu/nlopt,”.

]. Since this problem is nonconvex, there are many local optima. To find the global optima in the design parameter space, we used Multi-Level-Single-Linkage (MLSL) algorithm. This algorithm performs a sequence of local optimization from random points by a clustering heuristic that helps it to avoid repeated searches of the same optima [22

S. Kucherenko and Y. Sytsko, “Application of deterministic low-discrepancy sequences in global optimization,” Comput. Optim. Appl. 30, 297–318 (2005). [CrossRef]

]. The local optimization algorithm used here is BOBYQA [23

M. J. D. Powell, “The bobyqa algorithm for bound constrained optimization without derivatives,” Tech. rep., Department of Applied Mathematics and Theoretical Physics, Cambridge England (2009).

]. This algorithm performs derivative-free bound-constrained optimization using an iteratively constructed quadratic approximation of the objective function.

Figure 3 summarizes the optimization results. In all cases, the optimal structure returned by the optimization engine is always a silver coated silica sphere. Although multilayer structures can offer greater tunability of the optical response, bilayer structures already maximize the average cross-section over wide frequency range. From Table 1, we can see that nanospheres with an outer radius around 70nm have the largest normalized average scattering cross-section. The wavelength of the scattering peak can be further tuned by varying the aspect ratio. For absorption cross-sections, our optimization engine found many local optima with approximately equal FOM’s. These local optima have the same aspect ratio, and the thickness of the silver layer varies from zero to several nanometers. This can be explained by the quasi-static approximation. The absorption cross-section of a nanoparticle can be written as σ_abs(ω) = ωIm[α(ω)], where α(ω) is the polarizability [24

N. T. Fofang, T.-H. Park, O. Neumann, N. A. Mirin, P. Nordlander, and N. J. Halas, “Plexcitonic nanoparticles: Plasmonexciton coupling in nanoshell j-aggregate complexes,” Nano Lett. 8, 3481–3487 (2008). [CrossRef] [PubMed]

–26

V. S. Lebedev, A. S. Medvedev, D. N. Vasil’ev, D. A. Chubich, and A. G. Vitukhnovsky, “Optical properties of noble-metal nanoparticles coated with a dye j-aggregate monolayer,” Quantum Electron. 40, 246–248 (2010). [CrossRef]

]. When the nanoparticle diameter is much smaller than the wavelength, the quasi-static approximation is valid, under which the polarizability is proportional to the volume of the nanoparticle with the proportionality coefficient dependent on frequency and aspect ratio [24

]. Therefore, the normalized absorption cross-section averaged over a frequency range is only determined by the aspect ratio and independent of the nanoparticle diameter as long as the nanoparticle diameter is much smaller than the wavelength. In Table 2, the thickness of the silver layer is set to be ≥ 2nm. Nanospheres with thinner silver layers have approximately the same normalized average absorption cross-section but are difficult to fabricate. The structures that give the largest average scattering and absorption cross-sections are quite different. Since the FOM of the absorption cross-section is about twice as large as the FOM of the scattering cross-section, the structure that gives the largest average total cross-section is essentially identical to the structure providing the largest average absorption cross-section. The structures given in Table 1 are super-scatters and super-absorbers. For instance, the structure with [R, T] = [60.40nm, 8.68nm] has an average scattering cross-section of 8.65m²/g over 600–800nm, while its physical cross-section is only 2.07m²/g. The structure with [R, T] = [18.09nm, 2.00nm] provides an average total cross-section of 17.52m²/g over 600–800nm, which means that only 1g of such nanoparticles, when fully dispersed, can obscure an area as large as 17.52m².

Fig. 3 Optimization of average cross-sections over wide frequency range. The structure under consideration is a silver/silica multilayer nanosphere. The optimal structure found by the optimization engine is always silver coated silica sphere. For all subfigures, blue (red) lines show the optimized average cross-sections over the blue (red) shaded frequency range. (a)(b)(c) correspond to scattering, absorption and total cross-sections per volume respectively. (d)(e)(f) correspond to scattering, absorption and total cross-sections per mass respectively. The radius of the silica cores and the thickness of silver shells exhibiting the cross-sections shown above are given in Table 1.

Table 1 Optimization of average cross-sections

Normalized by volume
cross-section	range (nm)	silica (nm)	silver (nm)	FOM (1/nm)
σ_sca	400–600 600–800	31.25 60.32	26.65 9.65	0.0486 0.0464
σ_abs	400–600 600–800	6.07 14.80	2.00 2.00	0.0767 0.0817
σ_tot	400–600 600–800	6.09 16.12	2.00 2.00	0.0773 0.0846
Normalized by mass
cross-section	range (nm)	silica (nm)	silver (nm)	FOM (m²/g)
σ_sca	400–600 600–800	38.53 60.40	17.88 8.68	5.64 8.65
σ_abs	400–600 600–800	6.54 17.89	2.00 2.00	10.87 16.71
σ_tot	400–600 600–800	8.71 18.09	2.00 2.00	10.93 17.52

Table 2 Optimization of minimal cross-sections

Normalized by volume
cross-section	# of species	silica (nm)	silver (nm)	weights (%)	FOM (1/nm)
σ_sca	1	68.09	55.44	100	0.0192
	2	52.54	44.21	82.8	0.0229
	2	73.41	7.63	17.2	0.0229
	3	52.85	47.27	74.4	0.0235
		75.61	8.31	20.8
		35.55	23.21	4.8
σ_abs	1	60.14	10.00	100	0.0032
	2	37.03	4.00	57.5	0.0075
	2	44.02	15.99	42.5	0.0075
	3	22.95	10.26	37.4	0.0114
		26.08	4.39	34.3
		21.96	2.00	28.3
σ_tot	1	80.15	41.87	100	0.0201
	2	59.96	34.12	82.8	0.0270
	2	61.96	6.29	17.2	0.0270
	3	50.21	32.98	75.2	0.0285
		52.97	4.90	14.4
		55.76	7.32	10.4
	10	...	...	...	0.0310
Normalized by mass
cross-section	# of species	silica (nm)	silver (nm)	weights (%)	FOM (m²/g)
σ_sca	1	67.95	52.38	100	2.10
	2	52.77	38.70	87.4	2.80
	2	75.72	8.17	12.6	2.80
	3	76.71	21.93	46.1	2.88
		21.96	34.54	43.2
		70.45	7.31	10.7
σ_abs	1	42.80	6.55	100	0.58
	2	20.56	2.15	51.3	1.31
	2	37.16	13.61	48.7	1.31
	3	27.43	14.22	36.9	1.59
		49.06	11.02	34.9
		36.66	3.66	28.2
σ_tot	1	92.95	39.04	100	2.46
	2	55.54	29.60	83.5	3.59
	2	67.77	7.39	16.5	3.59
	3	53.36	29.01	60.1	3.79
		79.49	10.25	31.1
		24.64	25.04	8.8
	10	...	...	...	4.62

5. Optimization of the minimal cross-sections over wide frequency range

In this section, we aim to design nanoparticles with consistently large cross-sections over a wide frequency range. This is equivalent to maximizing the minimal cross-section over the target frequency range. From the previous section, we see that silver coated silica spheres have strong surface plasmon resonances, and the peaks of their resonances are highly tunable. Therefore, we expect that a mixture of several species of silver coated silica spheres with different resonance frequencies can effectively cover a wide frequency range.

This intuition can be formalized by optimization language. We take the target wavelength range to be 400–800nm. The FOM is the minimal cross-section over this range, normalized by either volume or mass.

F O M = min_{ω} σ_{normalized}

(14)

The structure under consideration is a mixture of N species of silver coated silica spheres where N = 1, 2, 3,.... The design parameters are the size parameters of individual species and relative weights (i.e. proportions) of each species in the mixture. The weights represent the relative weights in volume (mass) when the normalization is over volume (mass).

σ_{normalized} = \sum_{i = 1}^{N} w_{i} σ_{i, normalized}

(15)

When the size parameters of individual species are fixed, the problem of finding the optimal weights turns out to be a Linear Programming (LP) program. Therefore, we employed a two-level optimization structure. In the lower level, we used a standard LP solver to find out the optimal weights given the current size parameters. The resulting FOM as a function of size parameters is further optimized in the upper level with the same nonlinear method we used in the previous section. This separation into linear and nonlinear parts of the original optimization problem reduces the dimension of the parameter space and helps MLSL algorithm to find the global optimal in less iterations.

Figure 4 summarizes the results. We can see that the optimization engine tries to build a plateau over the target range to maximize the minimal cross-section. For scattering cross-sections, it is relatively easy to build such a plateau since the scattering peaks have wide bandwidth. When there is only one species (the blue lines in Fig. 4(a)(d)), the nanoparticle has a large size so that its l = 1 resonance can cover the large wavelength region and its l = 2 resonance can cover the small wavelength region. When the number of species increases, the new species try to cover the dips in the original scattering spectra with their resonant peaks. On the other hand, it is relatively difficult to build an absorption plateau because absorption peaks have narrow bandwidth. The FOM of the scattering cross-section is about twice as large as that of the absorption cross-section. The total cross-section can be enhanced significantly (35% for volume normalization and 46% for mass normalization) when N increases from 1 to 2. The enhancement when N increases from 2 to 3 is only moderate. The benefit of adding more species gradually saturates.

Fig. 4 Optimization of minimal cross-sections over a wide frequency range. The structure under consideration is a mixture of several species of silver coated silica spheres. The target frequency range is shaded in yellow. For all subfigures, blue, red, black lines corresponds to one, two, and three species of nanospheres. The black dashed lines in (c) and (f) correspond to ten species of nanospheres. (a)(b)(c) correspond to scattering, absorption, and total cross-sections per volume respectively. (d)(e)(f) correspond to scattering, absorption and total cross-sections per mass respectively. The radius of the silica cores and the thickness of silver shells corresponding to these cross-sections are given in Table 2.

6. Concluding remarks

In this article, we used optimization tools to tailor the optical response of silver/silica multilayer spheres. We show that the structure that gives the largest average cross-section over wide frequency range is the bilayer structure with a silver shell. We also show that using several species of nanoparticles can significantly enhance the minimal cross-section over the whole visible range although this enhancement saturates when the number of species increase. Because the FOM can be an arbitrary function of the frequency dependent cross-sections, the optimization approach described here can be used to design nanoparticles with more complicated optical response. We are investigating efficient numerical methods of calculating and optimizing the optical response of non-spherical geometries as future directions.

Acknowledgments

This research was funded by the Department of the Army Basic Research Program and sponsored by the Edgewood Chemical Biological Center. Support was also provided by the U.S. Army Research Office under contract W911NF-07-D-004.

References and links

1.	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]
2.	X. Huang, S. Neretina, and M. A. El-Sayed, “Gold nanorods: From synthesis and properties to biological and biomedical applications,” Adv. Mater. 21, 4880–4910 (2009). [CrossRef]
3.	C. Noguez, “Surface plasmons on metal nanoparticles: the influence of shape and physical environment,” J. Phys. Chem. C 111, 3806–3819 (2007). [CrossRef]
4.	S. Alyones, C. Bruce, and A. Buin, “Numerical methods for solving the problem of electromagnetic scattering by a thin finite conducting wire,” IEEE Trans. Antennas Propag. 55, 1856–1861 (2007). [CrossRef]
5.	C. W. Bruce and S. Alyones, “Extinction efficiencies for metallic fibers in the infrared,” Appl. Opt. 48, 5095–5098 (2009). [CrossRef] [PubMed]
6.	P. K. Jain, K. S. Lee, I. H. El-Sayed, and M. A. El-Sayed, “Calculated absorption and scattering properties of gold nanoparticles of different size, shape, and composition: Applications in biological imaging and biomedicine,” J. Phys. Chem. B 110, 7238–7248 (2006). [CrossRef] [PubMed]
7.	J. Zhu, J. Li, and J. Zhao, “Tuning the wavelength drift between resonance light absorption and scattering of plasmonic nanoparticle,” Appl. Phys. lett. 99, 101901 (2011). [CrossRef]
8.	C. E. Román-Velázquez and C. Noguez, “Designing the plasmonic response of shell nanoparticles: Spectral representation,” J. Chem. Phys. 134, 044116 (2011). [CrossRef] [PubMed]
9.	S. Oldenburg, R. Averitt, S. Westcott, and N. Halas, “Nanoengineering of optical resonances,” Chem. Phys. Lett. 288, 243 – 247 (1998). [CrossRef]
10.	E. Prodan and P. Nordlander, “Structural tunability of the plasmon resonances in metallic nanoshells,” Nano Lett. 3, 543–547 (2003). [CrossRef]
11.	E. Prodan, C. Radloff, N. J. Halas, and P. Nordlander, “A hybridization model for the plasmon response of complex nanostructures,” Science 302, 419–422 (2003). [CrossRef] [PubMed]
12.	R. Bardhan, N. K. Grady, T. Ali, and N. J. Halas, “Metallic nanoshells with semiconductor cores: Optical characteristics modified by core medium properties,” ACS Nano 4, 6169–6179 (2010). [CrossRef] [PubMed]
13.	R. Bardhan, S. Mukherjee, N. A. Mirin, S. D. Levit, P. Nordlander, and N. J. Halas, “Nanosphere-in-a-nanoshell: A simple nanomatryushka,” J. Phys. Chem. C 114, 7378–7383 (2010). [CrossRef]
14.	R. E. Hamam, A. Karalis, J. D. Joannopoulos, and M. Soljačić, “Coupled-mode theory for general free-space resonant scattering of waves,” Phys. Rev. A 75, 053801 (2007). [CrossRef]
15.	Z. Ruan and S. Fan, “Superscattering of light from subwavelength nanostructures,” Phys. Rev. Lett. 105, 013901 (2010). [CrossRef] [PubMed]
16.	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–7329 (2010). [CrossRef]
17.	Z. Ruan and S. Fan, “Design of subwavelength superscattering nanospheres,” Appl. Phys. lett. 98, 043101 (2011). [CrossRef]
18.	H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).
19.	C. Bohren and D. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Songs, 1983).
20.	E. D. Palik, Handbook of Optical Constants of Solids (Academic Press, New York, 1985).
21.	S. G. Johnson, “The nlopt nonlinear optimization package, http://ab-initio.mit.edu/nlopt,”.
22.	S. Kucherenko and Y. Sytsko, “Application of deterministic low-discrepancy sequences in global optimization,” Comput. Optim. Appl. 30, 297–318 (2005). [CrossRef]
23.	M. J. D. Powell, “The bobyqa algorithm for bound constrained optimization without derivatives,” Tech. rep., Department of Applied Mathematics and Theoretical Physics, Cambridge England (2009).
24.	N. T. Fofang, T.-H. Park, O. Neumann, N. A. Mirin, P. Nordlander, and N. J. Halas, “Plexcitonic nanoparticles: Plasmonexciton coupling in nanoshell j-aggregate complexes,” Nano Lett. 8, 3481–3487 (2008). [CrossRef] [PubMed]
25.	A. Yoshida and N. Kometani, “Effect of the interaction between molecular exciton and localized surface plasmon on the spectroscopic properties of silver nanoparticles coated with cyanine dye j-aggregates,” J. Chem. Phys. C 114, 2867–2872 (2010). [CrossRef]
26.	V. S. Lebedev, A. S. Medvedev, D. N. Vasil’ev, D. A. Chubich, and A. G. Vitukhnovsky, “Optical properties of noble-metal nanoparticles coated with a dye j-aggregate monolayer,” Quantum Electron. 40, 246–248 (2010). [CrossRef]

OCIS Codes
(290.2200) Scattering : Extinction
(290.5825) Scattering : Scattering theory

ToC Category:
Scattering

History
Original Manuscript: June 4, 2012
Revised Manuscript: July 20, 2012
Manuscript Accepted: July 23, 2012
Published: July 27, 2012

Virtual Issues
Vol. 7, Iss. 9 Virtual Journal for Biomedical Optics

Citation
Wenjun Qiu, Brendan G. DeLacy, Steven G. Johnson, John D. Joannopoulos, and Marin Soljačić, "Optimization of broadband optical response of multilayer nanospheres," Opt. Express 20, 18494-18504 (2012)
http://www.opticsinfobase.org/vjbo/abstract.cfm?URI=oe-20-16-18494

Sort: Author | Year | Journal | Reset

References

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]
X. Huang, S. Neretina, and M. A. El-Sayed, “Gold nanorods: From synthesis and properties to biological and biomedical applications,” Adv. Mater.21, 4880–4910 (2009). [CrossRef]
C. Noguez, “Surface plasmons on metal nanoparticles: the influence of shape and physical environment,” J. Phys. Chem. C111, 3806–3819 (2007). [CrossRef]
S. Alyones, C. Bruce, and A. Buin, “Numerical methods for solving the problem of electromagnetic scattering by a thin finite conducting wire,” IEEE Trans. Antennas Propag.55, 1856–1861 (2007). [CrossRef]
C. W. Bruce and S. Alyones, “Extinction efficiencies for metallic fibers in the infrared,” Appl. Opt.48, 5095–5098 (2009). [CrossRef] [PubMed]
P. K. Jain, K. S. Lee, I. H. El-Sayed, and M. A. El-Sayed, “Calculated absorption and scattering properties of gold nanoparticles of different size, shape, and composition: Applications in biological imaging and biomedicine,” J. Phys. Chem. B110, 7238–7248 (2006). [CrossRef] [PubMed]
J. Zhu, J. Li, and J. Zhao, “Tuning the wavelength drift between resonance light absorption and scattering of plasmonic nanoparticle,” Appl. Phys. lett.99, 101901 (2011). [CrossRef]
C. E. Román-Velázquez and C. Noguez, “Designing the plasmonic response of shell nanoparticles: Spectral representation,” J. Chem. Phys.134, 044116 (2011). [CrossRef] [PubMed]
S. Oldenburg, R. Averitt, S. Westcott, and N. Halas, “Nanoengineering of optical resonances,” Chem. Phys. Lett.288, 243 – 247 (1998). [CrossRef]
E. Prodan and P. Nordlander, “Structural tunability of the plasmon resonances in metallic nanoshells,” Nano Lett.3, 543–547 (2003). [CrossRef]
E. Prodan, C. Radloff, N. J. Halas, and P. Nordlander, “A hybridization model for the plasmon response of complex nanostructures,” Science302, 419–422 (2003). [CrossRef] [PubMed]
R. Bardhan, N. K. Grady, T. Ali, and N. J. Halas, “Metallic nanoshells with semiconductor cores: Optical characteristics modified by core medium properties,” ACS Nano4, 6169–6179 (2010). [CrossRef] [PubMed]
R. Bardhan, S. Mukherjee, N. A. Mirin, S. D. Levit, P. Nordlander, and N. J. Halas, “Nanosphere-in-a-nanoshell: A simple nanomatryushka,” J. Phys. Chem. C114, 7378–7383 (2010). [CrossRef]
R. E. Hamam, A. Karalis, J. D. Joannopoulos, and M. Soljačić, “Coupled-mode theory for general free-space resonant scattering of waves,” Phys. Rev. A75, 053801 (2007). [CrossRef]
Z. Ruan and S. Fan, “Superscattering of light from subwavelength nanostructures,” Phys. Rev. Lett.105, 013901 (2010). [CrossRef] [PubMed]
Z. Ruan and S. Fan, “Temporal coupled-mode theory for fano resonance in light scattering by a single obstacle,” J. Phys. Chem. C114, 7324–7329 (2010). [CrossRef]
Z. Ruan and S. Fan, “Design of subwavelength superscattering nanospheres,” Appl. Phys. lett.98, 043101 (2011). [CrossRef]
H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).
C. Bohren and D. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Songs, 1983).
E. D. Palik, Handbook of Optical Constants of Solids (Academic Press, New York, 1985).
S. G. Johnson, “The nlopt nonlinear optimization package, http://ab-initio.mit.edu/nlopt ,”.
S. Kucherenko and Y. Sytsko, “Application of deterministic low-discrepancy sequences in global optimization,” Comput. Optim. Appl.30, 297–318 (2005). [CrossRef]
M. J. D. Powell, “The bobyqa algorithm for bound constrained optimization without derivatives,” Tech. rep., Department of Applied Mathematics and Theoretical Physics, Cambridge England (2009).
N. T. Fofang, T.-H. Park, O. Neumann, N. A. Mirin, P. Nordlander, and N. J. Halas, “Plexcitonic nanoparticles: Plasmonexciton coupling in nanoshell j-aggregate complexes,” Nano Lett.8, 3481–3487 (2008). [CrossRef] [PubMed]
A. Yoshida and N. Kometani, “Effect of the interaction between molecular exciton and localized surface plasmon on the spectroscopic properties of silver nanoparticles coated with cyanine dye j-aggregates,” J. Chem. Phys. C114, 2867–2872 (2010). [CrossRef]
V. S. Lebedev, A. S. Medvedev, D. N. Vasil’ev, D. A. Chubich, and A. G. Vitukhnovsky, “Optical properties of noble-metal nanoparticles coated with a dye j-aggregate monolayer,” Quantum Electron.40, 246–248 (2010). [CrossRef]

Ali, T.
Alyones, S.
Averitt, R.
Bardhan, R.
Bohren, C.
Bruce, C.
Bruce, C. W.
Buin, A.
Chubich, D. A.
El-Sayed, I. H.
El-Sayed, M. A.
Fan, S.
Fofang, N. T.
Grady, N. K.
Grange, R.
Halas, N.
Halas, N. J.
Hamam, R. E.
Hsieh, C.-L.
Huang, X.
Huffman, D.
Jain, P. K.
Joannopoulos, J. D.
Karalis, A.
Kometani, N.
Kucherenko, S.
Lebedev, V. S.
Lee, K. S.
Levit, S. D.
Li, J.
Medvedev, A. S.
Mirin, N. A.
Mukherjee, S.
Neretina, S.
Neumann, O.
Noguez, C.
Nordlander, P.
Oldenburg, S.
Palik, E. D.
Park, T.-H.
Powell, M. J. D.
Prodan, E.
Psaltis, D.
Pu, Y.
Radloff, C.
Román-Velázquez, C. E.
Ruan, Z.
Soljacic, M.
Sytsko, Y.
van de Hulst, H. C.
Vasil’ev, D. N.
Vitukhnovsky, A. G.
Westcott, S.
Yoshida, A.
Zhao, J.
Zhu, J.

ACS Nano

R. Bardhan, N. K. Grady, T. Ali, and N. J. Halas, “Metallic nanoshells with semiconductor cores: Optical characteristics modified by core medium properties,” ACS Nano4, 6169–6179 (2010). [CrossRef] [PubMed]

Adv. Mater.

X. Huang, S. Neretina, and M. A. El-Sayed, “Gold nanorods: From synthesis and properties to biological and biomedical applications,” Adv. Mater.21, 4880–4910 (2009). [CrossRef]

Appl. Opt.

C. W. Bruce and S. Alyones, “Extinction efficiencies for metallic fibers in the infrared,” Appl. Opt.48, 5095–5098 (2009). [CrossRef] [PubMed]

Appl. Phys. lett.

J. Zhu, J. Li, and J. Zhao, “Tuning the wavelength drift between resonance light absorption and scattering of plasmonic nanoparticle,” Appl. Phys. lett.99, 101901 (2011). [CrossRef]
Z. Ruan and S. Fan, “Design of subwavelength superscattering nanospheres,” Appl. Phys. lett.98, 043101 (2011). [CrossRef]

Chem. Phys. Lett.

S. Oldenburg, R. Averitt, S. Westcott, and N. Halas, “Nanoengineering of optical resonances,” Chem. Phys. Lett.288, 243 – 247 (1998). [CrossRef]

Comput. Optim. Appl.

S. Kucherenko and Y. Sytsko, “Application of deterministic low-discrepancy sequences in global optimization,” Comput. Optim. Appl.30, 297–318 (2005). [CrossRef]

IEEE Trans. Antennas Propag.

S. Alyones, C. Bruce, and A. Buin, “Numerical methods for solving the problem of electromagnetic scattering by a thin finite conducting wire,” IEEE Trans. Antennas Propag.55, 1856–1861 (2007). [CrossRef]

J. Chem. Phys.

C. E. Román-Velázquez and C. Noguez, “Designing the plasmonic response of shell nanoparticles: Spectral representation,” J. Chem. Phys.134, 044116 (2011). [CrossRef] [PubMed]

J. Chem. Phys. C

A. Yoshida and N. Kometani, “Effect of the interaction between molecular exciton and localized surface plasmon on the spectroscopic properties of silver nanoparticles coated with cyanine dye j-aggregates,” J. Chem. Phys. C114, 2867–2872 (2010). [CrossRef]

J. Phys. Chem. B

P. K. Jain, K. S. Lee, I. H. El-Sayed, and M. A. El-Sayed, “Calculated absorption and scattering properties of gold nanoparticles of different size, shape, and composition: Applications in biological imaging and biomedicine,” J. Phys. Chem. B110, 7238–7248 (2006). [CrossRef] [PubMed]

J. Phys. Chem. C

C. Noguez, “Surface plasmons on metal nanoparticles: the influence of shape and physical environment,” J. Phys. Chem. C111, 3806–3819 (2007). [CrossRef]
R. Bardhan, S. Mukherjee, N. A. Mirin, S. D. Levit, P. Nordlander, and N. J. Halas, “Nanosphere-in-a-nanoshell: A simple nanomatryushka,” J. Phys. Chem. C114, 7378–7383 (2010). [CrossRef]
Z. Ruan and S. Fan, “Temporal coupled-mode theory for fano resonance in light scattering by a single obstacle,” J. Phys. Chem. C114, 7324–7329 (2010). [CrossRef]

Nano Lett.

N. T. Fofang, T.-H. Park, O. Neumann, N. A. Mirin, P. Nordlander, and N. J. Halas, “Plexcitonic nanoparticles: Plasmonexciton coupling in nanoshell j-aggregate complexes,” Nano Lett.8, 3481–3487 (2008). [CrossRef] [PubMed]
E. Prodan and P. Nordlander, “Structural tunability of the plasmon resonances in metallic nanoshells,” Nano Lett.3, 543–547 (2003). [CrossRef]

Phys. Rev. A

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

Phys. Rev. Lett.

Z. Ruan and S. Fan, “Superscattering of light from subwavelength nanostructures,” Phys. Rev. Lett.105, 013901 (2010). [CrossRef] [PubMed]
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]

Quantum Electron.

V. S. Lebedev, A. S. Medvedev, D. N. Vasil’ev, D. A. Chubich, and A. G. Vitukhnovsky, “Optical properties of noble-metal nanoparticles coated with a dye j-aggregate monolayer,” Quantum Electron.40, 246–248 (2010). [CrossRef]

Science

E. Prodan, C. Radloff, N. J. Halas, and P. Nordlander, “A hybridization model for the plasmon response of complex nanostructures,” Science302, 419–422 (2003). [CrossRef] [PubMed]

Other

M. J. D. Powell, “The bobyqa algorithm for bound constrained optimization without derivatives,” Tech. rep., Department of Applied Mathematics and Theoretical Physics, Cambridge England (2009).
H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).
C. Bohren and D. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Songs, 1983).
E. D. Palik, Handbook of Optical Constants of Solids (Academic Press, New York, 1985).
S. G. Johnson, “The nlopt nonlinear optimization package, http://ab-initio.mit.edu/nlopt ,”.

2011, Zhu, Appl. Phys. lett.
2011, Román-Velázquez, J. Chem. Phys.
2011, Ruan, Appl. Phys. lett.
2010, Ruan, Phys. Rev. Lett.
2010, Ruan, J. Phys. Chem. C
2010, Bardhan, ACS Nano
2010, Bardhan, J. Phys. Chem. C
2010, Pu, Phys. Rev. Lett.
2010, Yoshida, J. Chem. Phys. C
2010, Lebedev, Quantum Electron.
2009, Huang, Adv. Mater.
2009, Bruce, Appl. Opt.
2008, Fofang, Nano Lett.
2007, Noguez, J. Phys. Chem. C
2007, Alyones, IEEE Trans. Antennas Propag.
2007, Hamam, Phys. Rev. A
2006, Jain, J. Phys. Chem. B
2005, Kucherenko, Comput. Optim. Appl.
2003, Prodan, Nano Lett.
2003, Prodan, Science
1998, Oldenburg, Chem. Phys. Lett.

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper

Figures


Fig. 1	Fig. 2	Fig. 3


Fig. 4

« Previous Article | Next Article »

OSA is a member of CrossRef.

OpticsInfoBase

Virtual Journal for Biomedical Optics

| EXPLORING THE INTERFACE OF LIGHT AND BIOMEDICINE

Optimization of broadband optical response of multilayer nanospheres

Article Tools

Article Navigation

Abstract

1. Introduction

2. Calculation of optical response via transfer matrix method

3. Optical response of silver/silica bilayer nanospheres

4. Optimization of average cross-sections over wide frequency range

5. Optimization of the minimal cross-sections over wide frequency range

6. Concluding remarks

Acknowledgments

References and links

References

ACS Nano

Adv. Mater.

Appl. Opt.

Appl. Phys. lett.

Chem. Phys. Lett.

Comput. Optim. Appl.

IEEE Trans. Antennas Propag.

J. Chem. Phys.

J. Chem. Phys. C

J. Phys. Chem. B

J. Phys. Chem. C

Nano Lett.

Phys. Rev. A

Phys. Rev. Lett.

Quantum Electron.

Science

Other

Cited By

Figures

Select a Conference
Session or Paper #
Year