Home | About | Journals | Submit | Contact Us | Français |

**|**Biomed Opt Express**|**v.2(6); 2011 June 1**|**PMC3114226

Formats

Article sections

- Abstract
- 1. Introduction
- 2. Absorption efficiency computation
- 3. The adaptive Monte-Carlo method
- 4. Sensitivity design of nanoshells for photothermal therapy
- 5. Conclusion
- References and links

Authors

Related links

Biomed Opt Express. 2011 June 1; 2(6): 1584–1596.

Published online 2011 May 17. doi: 10.1364/BOE.2.001584

PMCID: PMC3114226

Thomas Grosges,^{1,}^{*} Dominique Barchiesi,^{1} Sameh Kessentini,^{1} Gérard Gréhan,^{2} and Marc Lamy de la Chapelle^{3}

Received 2011 February 2; Revised 2011 April 22; Accepted 2011 May 6.

Copyright ©2011 Optical Society
of America

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-No Derivative Works 3.0 Unported License, which permits download and redistribution, provided that the original work is properly cited. This license restricts the article from being modified or used commercially.

This article has been cited by other articles in PMC.

The optimization of the coated metallic nanoparticles and nanoshells is a current challenge for biological applications, especially for cancer photothermal therapy, considering both the continuous improvement of their fabrication and the increasing requirement of efficiency. The efficiency of the coupling between illumination with such nanostructures for burning purposes depends unevenly on their geometrical parameters (radius, thickness of the shell) and material parameters (permittivities which depend on the illumination wavelength). Through a Monte-Carlo method, we propose a numerical study of such nanodevice, to evaluate tolerances (or uncertainty) on these parameters, given a threshold of efficiency, to facilitate the design of nanoparticles. The results could help to focus on the relevant parameters of the engineering process for which the absorbed energy is the most dependant. The Monte-Carlo method confirms that the best burning efficiency are obtained for hollow nanospheres and exhibit the sensitivity of the absorbed electromagnetic energy as a function of each parameter. The proposed method is general and could be applied in design and development of new embedded coated nanomaterials used in biomedicine applications.

Extensive theoretical and experimental studies of the interaction between an incident wave with metallic nanostructures, such as nanoparticles, nanodots [1], nanorings [2], nanocubes [3], nanoshells [4], have been achieved in physics and chemistry. Under illumination, these structures, due to local plasmon resonances, are known to exhibit a high enhancement of the electromagnetic field at their surfaces. This field strength depends on the size parameters of the particle and on two related quantities: the wavelength and the permittivity of the used materials (in the case of nanoshell, the core could be made of silica or just a vacuum). Due to this structural tunability of the plasmon resonances, such nanostructures can be of interest in a wide range of applications in biomedicine [5–7].

Therefore, these researches permit to develop applications in nanoscale plasmon waveguide [8] and nanosensing modalities [9–11]. In the past, the choice of the metallic materials and synthesized alloys, size and shape geometries of materials were studied experimentally [12–14] and interpreted as plasmon hybridization [15]. Therefore, the development of metallic nanodevices, where the relative role of electromagnetic and resonant enhancement response could be precisely delineated, necessitates an accurate control of the local electromagnetic field enhancement near the metallic surface [16, 17].

For cancer photothermal therapy, nanodevices like coated metallic nanoparticles and nanoshells are currently used to burn cancer cells. Indeed, most biological tissues have a relatively low light absorption coefficient in the visible (VIS) and near infrared (NIR) regions (600–1300 nm) known as the tissue optical window or therapeutic window. Over this window, organic molecules have limited absorption [18], whereas gold nanodevices, which are biocompatible nontoxic and easily conjugated to antibodies, absorb light millions of times stronger than the organic molecules. Then, almost all the absorbed light is converted to heat via series of nonradiative processes [18] to burn cells in which they are embedded. The plasmon resonance tuning helps to increase the absorbed light and therefore, the nanodevice capability to burn diseased cells.

Until recently, in the absence of optimization of such devices, the absorption efficiency (i.e. the rate of absorbed energy by the particle relatively to the incoming illumination intensity and to the surface of the particle) is only of one order of magnitude [19, 20]. In a recent study, it has been shown that it becomes possible to obtain an absorption efficiency of two order of magnitude in the most convenient wavelength domain for biomedical applications (around 900 nm) [21]. Such a design optimization permits to obtain an increase in the absorption up to two order of magnitude in comparison with the classical coated nanoparticles. However, the sensitivity of such devices to manufacturing uncertainties has never been investigated.

In this study, we focus on the sensitivity design of nanoshells used in photothermal therapy (i.e. the absorption of these particles). The main principle of the method is to select the particles parameters that keep the absorption efficiency greater than a threshold. These parameters are within a range that defines the manufacturing tolerance.

The paper is organized as follows. Section 2 is devoted to the description of the absorption efficiency computation and the requirements on absorption efficiency computation are defined. Then, in section 3, the optimisation scheme by the Monte-Carlo method is presented and explained. The numerical setting related to the considered biomedical application (photothermal therapy for deep cancer using nanoshells) and results will be given and discussed in section 4. Finally, section 5 concludes the paper.

The total electromagnetic field (near-field and far-field) properties of a metallic nanoobject depends intrinsically on the geometry and on the optical properties of the involved materials. Nanoshells are composed of a core with radius *r*_{1} and of a metallic coating or shell of thickness *e* (see Fig. 1). The core could be made of silica or just vacuum (the commonly used name of such nanoshells is hollow nanospheres with *n*_{1} = 1.0) whereas the shell is made of gold. The permittivity of the core and the permittivity of the coating are denoted *ε*_{1} and *ε*_{2}, respectively. The metallic material is characterized by a complex permittivity *ε*_{2} = *ε _{r}* +

$${W}_{\mathit{\text{abs}}}\hspace{0.17em}=\hspace{0.17em}\frac{1}{2}\hspace{0.17em}\Re \hspace{0.17em}\left\{{\int}_{\Omega}[{\mathbf{\text{E}}}_{\mathbf{\text{i}}}\hspace{0.17em}\times \text{\hspace{0.17em}}{\mathbf{\text{H}}}_{\mathbf{\text{s}}}^{*}\hspace{0.17em}+\hspace{0.17em}{\mathbf{\text{E}}}_{\mathbf{\text{s}}}\hspace{0.17em}\times \text{\hspace{0.17em}}{\mathbf{\text{H}}}_{\mathbf{\text{i}}}^{*}]\hspace{0.17em}d\Omega \hspace{0.17em}\right\}\hspace{0.17em}-\hspace{0.17em}\frac{1}{2}\hspace{0.17em}\Re \hspace{0.17em}\left\{{\int}_{\Omega}[{\mathbf{\text{E}}}_{\mathbf{\text{s}}}\hspace{0.17em}\times \text{\hspace{0.17em}}{\mathbf{\text{H}}}_{\mathbf{\text{s}}}^{*}]\hspace{0.17em}d\Omega \right\}\hspace{0.17em},$$

(1)

where **E _{i}**,

$${C}_{\mathit{\text{abs}}}\hspace{0.17em}=\hspace{0.17em}\frac{{W}_{\mathit{\text{abs}}}}{{I}_{i}}\hspace{0.17em}=\hspace{0.17em}\frac{2\pi}{{k}^{2}}\hspace{0.17em}\sum _{n=1}^{\infty}(2n\hspace{0.17em}+\hspace{0.17em}1)\hspace{0.17em}\left\{\Re \hspace{0.17em}[{a}_{n}\hspace{0.17em}+\hspace{0.17em}{b}_{n}]\hspace{0.17em}-\hspace{0.17em}\left[|{a}_{n}{|}^{2}\hspace{0.17em}+\hspace{0.17em}|{b}_{n}{|}^{2}\right]\right\}\hspace{0.17em},$$

(2)

where *I _{i}* is the incident intensity,
$k\hspace{0.17em}=\hspace{0.17em}2\pi \sqrt{{\varepsilon}_{m}}/\lambda $ being the wave vector in the surrounding medium, and

$${Q}_{\mathit{\text{abs}}}\hspace{0.17em}=\hspace{0.17em}\frac{{C}_{\mathit{\text{abs}}}}{S}\hspace{0.17em}=\hspace{0.17em}\frac{2}{{k}^{2}\hspace{0.17em}{({r}_{1}\hspace{0.17em}+\hspace{0.17em}e)}^{2}}\hspace{0.17em}\sum _{n=1}^{\infty}(2n\hspace{0.17em}+\hspace{0.17em}1)\hspace{0.17em}\left\{\Re \hspace{0.17em}[{a}_{n}\hspace{0.17em}+\hspace{0.17em}{b}_{n}]\hspace{0.17em}-\hspace{0.17em}\left[|{a}_{n}{|}^{2}\hspace{0.17em}+\hspace{0.17em}|{b}_{n}{|}^{2}\right]\right\}\hspace{0.17em}.$$

(3)

In the limit of small particles (*r*_{2}*/λ* << 1) [22], the computation of *Q _{abs}* is often achieved by neglecting the order

We show in Fig. 2 that even for small particles as nanoshells, this dipolar approximation is not valid. We summary in Table 1 the acceptable intervals of parameters for inner radius *r*_{1}, thickness *e*, illuminating wavelength *λ* and the relative real and imaginary parts of the bulk permittivity of gold *ε*_{2}(*λ*) [31]. Figures 2(a) and 2(b) show the relative error between the absorption efficiency *Q _{abs}* computed from the dipolar approximation and the series (Eq. (3), with 60 terms, ensuring a convergence of the series, better than 10

Relative error between the exact computation of *Q*_{abs} and its approximation in the small particle limit for radii, as a function of the real *ε*_{r} and imaginary part *ε*_{i} of the material permittivity for particle radii: (a) *r*_{1}*/λ* = 0.089 **...**

Therefore, for sensitivity design, the multipolar computation of *Q _{abs}* must be achieved (i.e. more than one term must be computed in the series). In the following, the computations of the

The main purpose is to determine the parameters of the nanoparticles that are critical for medical application. For this, we introduce a boundary adaptive method based on a Monte-Carlo scheme. The proposed approach differs strongly from the metaheuristic optimization algorithms we applied in the previous studies, either to the inverse problem solving in near-field optics [33], or to optimize the Surface Plasmon Resonance planar biosensor [17, 34–37]. In those papers, the goal was to compute the best set of parameters to get the most efficient biosensor. Nevertheless, no propagation of uncertainty has been determined, despite the known uncertainties on permittivities and geometrical parameters. Recently, we have proposed a method to design the nanoshells with silica core for imaging and burning applications [21]. In that study, an evolutionary method (i.e. developping similar operational principles based on the evolution of the searched parameters as objective variables to reach a target) was used to optimize the size parameters and the permittivity of the shell using existing permittivities. However, the optical index of the core was fixed to that of silica (which is not the case in this study). The goal was to obtain at least 80% of the mathematical maximum of the absorption *Q _{abs}* (resp. scattering) efficiency, for the burning (resp. imaging) purpose. This method is now used to find the best set of parameters to reach

Otherwise, in the present paper, the proposed method consists in fixing a tolerance on the optimum of the absorption efficiency, to determine a class of acceptable parameters. Then, the tolerance on the various parameters can be deduced and interpreted either in terms of sensitivity of the model or acceptable uncertainty in the process of fabrication. For this, a family of parameters of the nanoshells (size and permittivities) is randomly generated, and the algorithm retains the maximal and minimal values (i.e. the boundaries) of the parameters for which the efficiency remains greater than a given fraction *α* of the optimal efficiency. The process is repeated (by regenerating randomly another set of nanoshells given the new boundary constraints and looking for additional “good” nanoshells) until a given number of “good” nanoshells are obtained and analyzed to depict the sensitivity to each size or optical parameter. This approach is derived from a classical metaheuristic engineering approach of the study of uncertainties propagation through a numerical model of the system, to compute the engineering tolerances. Actually, the Monte-Carlo method relies on repeated random sampling to compute their results (e.g. probability of fitting some conditions) in a fixed domain (or boundaries) given usually by physical constraints. It is known to be an effective substitute to the traditional method of experimental plan, indeed to solve inverse problems [38].

In this paper, we present the method permitting to compute the intervals of optimized size parameters and core index that guarantee absorption efficiency *Q _{abs}*, given an acceptable percentage

This numerical method is summarized as following, considering a threshold *α* = 90%, delineating an acceptable value of *Q _{abs}*:

- Define an initial domain of possible parameters (size parameters ${r}_{1}^{0}$ and
*e*^{0}and permittivities ${\varepsilon}_{1}^{0}$, ${\varepsilon}_{2}^{0}(\lambda )$) by setting the boundaries according to physical constraints. The number of parameters set*npt*is fixed as the maximum of the domain of variation of the initial parameters divided by their uncertainty*u*in the process of fabrication [38]:*npt*= max { $\left[\text{max}\left({r}_{1}^{0}\right)-\hspace{0.17em}\text{min}\left({r}_{1}^{0}\right)\right]\hspace{0.17em}/u({r}_{1})$; [max(*e*^{0}) − min(*e*^{0})]*/u*(*e*); $\left[\text{max}({n}_{1}^{0})\hspace{0.17em}-\hspace{0.17em}\text{min}\left({n}_{1}^{0}\right)\right]/u({n}_{1})$; [max(*λ*^{0}) − min(*λ*^{0})]*/u*(*λ*)}. Within the range of parameters, the maximum of*Q*is computed with the above mentioned evolutionary method [21]._{abs} - Random generation of a family of
*npt*unknown parameters sets**p**(*t*), using uniform laws in the considered domain.*t*is the iteration index initially set to 1. - Computation of the optimum geometrical parameters with Nelder-Mead method: this method requires an initial values for the optimum search. These initial values are the parameters set that gives the maximum value max
(_{t}*Q*), within the family: [_{abs}*r*_{1}*, e,λ,ε*_{2}(*λ*)],*ε*_{1}being unchanged. This maximum is max(_{NM}*Q*), and varies slightly with_{abs}*t*. Actually, the basic Nelder-Mead method is a multidimensional unconstrained nonlinear method that cannot handle with boundaries in the search domain and therefore, if*n*_{1}would be also a parameter to be optimized, the output of the algorithm would give*n*_{1}= 0, which is not compatible with the physical constraint on the optical index of a dielectric medium*n*_{1}≥ 1. The optimum set generated by this method will be denoted*p*(_{max}*t*) and the corresponding absorption efficiency max(_{NM}*Q*)._{abs} - Selection and storage of the best parameters sets of the family
**p**(*t*) for which*Q*>_{abs}*α*max(_{NM}*Q*). The retained sets form a family denoted_{abs}**p**(_{α}*t*). The set*p*(_{max}*t*) is not added to this family**p**(_{α}*t*) but is used only to determine max(_{NM}*Q*)._{abs} - Updating boundaries of the domain using the minimum and the maximum of each parameter in
**p**(_{α}*t*). - Increment
*t*and loop on step 2, until the size of all the selected sets of parameters_{t}**p**(_{α}*t*) is greater than*npt*. The final number of iteration is*N*._{t}

The associated tolerance in fabrication will be deduced from the boundaries of this last family.

At this stage, two points should be clarified. First, the advantage of the adaptive method v.s. the classical Monte-Carlo method. To reach an accuracy on the boundaries determination lower than 1%, the number of random parameters should be greater than 100,000, the convergence of the classical Monte-Carlo method being in $1/\sqrt{\mathit{\text{npt}}}$. The advantage of our adaptive method lies in the adaptation of the domain of search, keeping constant the number of random parameters at each iteration, from a generation to the following. Therefore the convergence is better at each new iteration and the precision on the boundaries is increased.

Second, the numerical cost of the method has to be compared to the cost of the four loops required in a systematic study. For this, we use the numerical parameters, that are physically acceptable for the initialization of the algorithm. We consider
${r}_{1}\hspace{0.17em}\in \hspace{0.17em}[1;\hspace{0.17em}100]\hspace{0.17em}=\hspace{0.17em}[\text{min}\left({r}_{1}^{0}\right);\hspace{0.17em}\text{max}\left({r}_{1}^{0}\right)]\hspace{0.17em}\text{nm}$; *e* [1;50] nm as small particles are preferred in burning applications.

For the benchmark of the algorithm, we used the same parameters as in Ref. [21] where non heuristic method was used to compute the tolerance on the parameters. Table 2 shows the various parameters of the study and the obtained results. In that study the optical index of the silica core (resp. external medium) was *ε*_{1} = 2.4 (resp. *ε*_{3} = 1.77). For illustration, we focus on a specific result for gold nanoshell and *λ* = 850 nm. The threshold was *α* = 80%. The computation of the standard deviation of the boundaries obtained for all iterations can be held as an indicator of the confidence on the computed boundaries. Indeed, the standard deviation of the boundaries can be considered as uncertainty and is an indicator for limiting the significative digits of the results. Increasing *npt* results in a
$1/\sqrt{\mathit{\text{npt}}}$ decreasing of the standard deviations of the boundaries, the number of iterations being almost constant (*N _{t}* around 40).

The Monte-Carlo method, gives the tolerance on the radius *r*_{1} and the thickness *e*. These tolerances are deduced from the intervals of acceptable values of each parameter. The tolerances correspond to those obtained from a systematic study like experience plane [21], but the number of required evaluations of *Q _{abs}* is reduced. The computational scheme requires less than 40 iterations (with 990 evaluations of

These results are validating the Monte-Carlo approach which helps to compute sensitivity of the model to the input parameters. From Table 2, we can deduce that the critical parameter is the thickness of gold, the tolerance on *r*_{1} being around four times greater. After convergence of the Monte-Carlo scheme, the boundaries of the last hypercube can be used to define the tolerance on each parameter.

A last bench has been made before the use of this method for more general problems. Indeed, it is well known that the reiteration of metaheuristic methods has to be made to characterize stability and dispersion of the results. Thousand realizations of the proposed Monte-Carlo code, with different initializations of the random generator, have been made to check the stability of the algorithm. The obtained results show a standard deviation of the boundaries, over the thousand realizations, which is smaller than the initial tolerance on each parameter. Therefore, the method could be applied with sufficient confidence to the sensitivity design of nanoshells. Let us investigate the case of hollow nanoshells (*n*_{1} ≈ 1) and coated silica (*n*_{1} = 1.54) nanospheres in fat instead of water [39].

We focus on the design of gold coated nanoparticles (or nanoshells) for applications in the biomedical domain of wavelengths (*λ* [800;1000] nm), for cell burning purpose. The control of absorption efficiency by tuning the wavelength and the size in visible and infrared region has been proved experimentally [19, 34]. Nevertheless, the optical constant of gold nanoparticles is not known exactly. Moreover, strong variations have been found for bulk and obviously for thin layers, depending on the mode of deposition. With such variations in the values of the permittivity for gold, the question is: is the numerical optimization using the bulk permittivity for the gold shell adapted? In fact, despite the variations between values of the permittivity for gold nanoparticles and bulk values, the results of computation of cross sections are in good agreements with experiment results, especially in the case of spectroscopic studies of such nano-devices [40–42]. Hopefully the following results confirm that the geometrical parameters are more critical in the design than the optical index of gold.

The permittivity of the core *ε*_{1} [1;16] as the commonly used core is made of silica or vacuum (with vacuum core, the corresponding nanoshells are commonly known as hollow nanospheres). The uncertainty on the optical core index
${n}_{1}\hspace{0.17em}=\hspace{0.17em}\sqrt{{\varepsilon}_{1}}$ is fixed to 0.01. The permittivity of the gold coating (or shell) is obtained from the Palik’s data [31], in the [800;1000] nm range of wavelength *λ : ε*_{2}(*λ*). Therefore, the real part of the gold permittivity lays between −42 and −24. The imaginary part of this permittivity is between 1.5 and 3. A fit of this permittivity with *λ*’s step equal to 1 nm is used. On the contrary of Ref. [21], the material of the shell is fixed and our results reveal the shallow sensitivity on this parameter.

Let us investigate first the most general case, where the index of the core must also be determined, the initial set of parameters is summarized in Table 3. Therefore, the size of each family of parameters at each iteration, *npt*, can be evaluated: *npt* = max{(100 – 1)/0.1,(50 – 1)/0.1,(1000 – 800)/1,(4 – 1)/0.01,(800 – 400)/1} = 990.

In the following, we present the histograms for each parameter (see Fig. 3). These results are obtained for a value of *Q _{abs}* which falls within the interval [0.9 max

- hopefully, the sensitivity of
*Q*to the wavelength (and consequently to the permittivity of gold) is low. In the investigated range of wavelengthes ([800;1000] nm), the permittivity has monotonic variations. Therefore the nanoshell can be efficient in this whole domain of wavelengthes. Moreover, the exact knowledge of permittivity of the shell seems not to be critical. This behavior can be explained by the broadening of the multipolar resonance of the nanoshell due to the non negligible imaginary part of the permittivity._{abs} - The less the inner permittivity is, the better is the efficiency of the nanoshell (see Fig. 3(b) with
*n*_{1}≈ 1.0 and Fig. 3(e) with 12.7 ≤*Q*≤ 14.0). This confirms that the hollow nanosphere seems to be the more efficient. Schwartzberg_{abs}*et al.*have shown that the absorption band of these particles can be tuned by adjusting the thickness of the gold shell and the inner radius and thus would enable both strong scattering and absorption efficiency [34]. - The computed tolerances on radii are within these of fabrication processes (0.6 nm [34, 35]). But the thickness of the gold shell is the most critical parameter:
*e*= 2.1 ± 0.6.*r*_{1}= 22.5 ± 6.5 and the best parameters are*r*_{1}= 22.2 nm,*e*= 2.0 nm, to get max(_{NM}*Q*) = 14.0. The best parameters for_{abs}*λ*,*r*_{1}and*e*are close the center of the intervals of tolerance and*n*_{1}≈ 1.0. The standard deviation of the boundaries of intervals are respectively 0.5 for*r*_{1}, 0.05 for*e*and 10^{−3}for*n*_{1}. These values are around ten times lower than the size of each interval, ensure the validity of the above intervals of tolerance, and confirm the sensitivity of*Q*to each parameter. In particular, the high sensitivity of_{abs}*Q*to the thickness of gold is observed and can be assimilated to the same behavior observed for Surface Plasmon Resonance based planar biosensors [17]._{abs} - The histogram of
*Q*shows a maximum of 14.0 and of course a minimum of 12.9 which is related to the threshold_{abs}*α*= 0.9.

Histograms of (a) the wavelength, (b) the optical index of the core (
${n}_{1}\hspace{0.17em}=\hspace{0.17em}\sqrt{{\varepsilon}_{1}}$) (c) the radius of the core, (d) the thickness of the shell, (e) the absorption efficiency. The relative frequency is plotted in percents, and the number **...**

With such parameters, the absorption efficiency, relatively to the reference particle size (see Loo et al. [19]), corresponds to a gain of 60 to 100.

Figures 4(a)–4(d) show the convergence of the boundaries of the domain. In this case, the maximum number of iterations required to get at least *npt* parameters sets _{t}**p*** _{α}* (

A second application is carried to study one of the commonly manufactured nanoparticles: the nanoshells with silica core. For this, the optical index *n*_{1} is fixed to 1.54. This is possible as its variation with the wavelength within the range [800;1000] nm is negligible. The proposed Monte-Method results are *r*_{1} [18.2(0.6);28.2(0.6)] nm and *e* [2.1(0.06);3.3(0.06)] nm. The standard deviation of each boundary in the family _{t}**p*** _{α}* (

Moreover, max* _{NM}*(

We proposed an iterative Monte-Carlo method with boundary adaptation, to compute the sensitivity of the absorption efficiency of nanoparticles. The goal was to determine the sensitivity of the geometrical and material parameters, that could maintain a sufficiently high absorption efficiency, to get an elevation of its temperature for burning purposes. The proposed algorithm enables to spare evaluations of the target (*Q _{abs}*) with regards to a systematic study, through simple loops. The obtained results are coherent with former experimental and theoretical results and are statistically significant. In particular, the method shows that the hollow nanoshells (no material in the core) are more efficient. The extent of the domain of possible parameters shows that the shell should be synthesized with more precision to guarantee a maximal absorption. On the other hand, the efficiency of the nanoparticles also depends on the gold process of deposition and on the inner radius. The method introduced in this study is general and can be applied to a wide range of problems and especially to nanoparticles in order to obtain an optimized design of the size, the shape, the geometry and to select the most adapted for a specific application. The extension to non-spherical particles with complex geometry may give perspectives for both research and applicative purposes, especially in the optical imaging and in the biomedical domain. The extension to complex or non-spherical geometries will necessitate the definition of additional geometrical parameters (like roughness), the development of an adapted numerical method (e.g. FEM, DDA,..) to compute the target to optimize (Qabs, or absorbed enery) for the considered geometry. Nevertheless, this study opens a way in the easier manufacturing and development of a new family of sophisticate nanostructures (nanoshells, nanorings, nanodots, etc.) with optimized design.

Authors thank the Centre de Calcul Intensif ROMEO2 for computational facilities, the Région Champagne-Ardennes, the Conseil Régional de l’Aube and the *Nanoantenna* European Project (FP7 Health-F5-2009-241818) for financial supports.

1. Link S., El-Sayed M. A., “Spectral properties and relaxation dynamics of surface plasmon electronic oscillations in gold and silver nanodots and nanorods,” J. Phys. Chem. B 103, 8410–8426 (1999).10.1021/jp9917648 [Cross Ref]

2. Aizpurua J., Hanarp P., Sutherland D. S., Kall M., Bryant G. W., de Abajo F. J. G., “Optical properties of gold nanorings,” Phys. Rev. Lett. 90, 057401 (2003).10.1103/PhysRevLett.90.057401 [PubMed] [Cross Ref]

3. Sun Y. G., Xia Y. N., “Shape-controlled synthesis of gold and silver nanoparticles,” Science 298, 2176–2179 (2002).10.1126/science.1077229 [PubMed] [Cross Ref]

4. Oldenburg S. J., Averitt R. D., Westcott S. L., Halas N. J., “Nanoengineering of optical resonances,” Chem. Phys. Lett. 288, 243–247 (1998).10.1016/S0009-2614(98)00277-2 [Cross Ref]

5. Loo C., Lowery A., Halas N. J., West J. L., Drezek R., “Immunotargeted nanoshells for integrated cancer imaging and therapy,” Nano Lett. 5, 709–711 (2005).10.1021/nl050127s [PubMed] [Cross Ref]

6. Sershen S. R., Westcott S. L., Halas N. J., West J. L., “Temperature-sensitive polymer-nanoshell composites for photothermally modulated drug delivery,” J. Biomed. Mater. Res. 51, 293–298 (2000).10.1002/1097-4636(20000905)51:3<293::AID-JBM1>3.0.CO;2-T [PubMed] [Cross Ref]

7. O’Neal D. P., Hirsch L. R., Halas N. J., Payne J. D., West J. L., “Photo-thermal tumor ablation in mice using near infrared-absorbing nanoparticles,” Cancer Lett. 209, 171–176 (2004).10.1016/j.canlet.2004.02.004 [PubMed] [Cross Ref]

8. Brongersma M. L., Hartman J. W., Atwater H. A., “Electromagnetic energy transfer and switching in nanoparticle chain arrays below the diffraction limit,” Phys. Rev. B 62, 16356–16359 (2000).10.1103/PhysRevB.62.R16356 [Cross Ref]

9. Okamoto T., Yamaguchi I., Kobayashi T., “Local plasmon sensor with gold colloid monolayers deposited upon glass substrates,” Opt. Lett. 25, 372–374 (2000).10.1364/OL.25.000372 [PubMed] [Cross Ref]

10. Jackson J. B., Westcott S. L., Hirsch L. R., West J. L., Halas N. J., “Controlling the surface enhanced Raman effect via the nanoshell geometry,” Appl. Phys. Lett. 82, 257–259 (2003).10.1063/1.1534916 [Cross Ref]

11. Tam F., Halas N. J., “Plasmon response of nanoshell dopants in organic films: a simulation study,” Prog. Org. Coat. 47, 275–278 (2003).10.1016/j.porgcoat.2003.08.001 [Cross Ref]

12. Link S., Wang Z.L., El-Sayed M.A., “Alloy formation of gold-silver nanoparticles and the dependence of the plasmon absorption on their composition,” J. Phys. Chem. B 103, 3529–3533 (1999).10.1021/jp990387w [Cross Ref]

13. Kreibig U., Vollmer M., Optical Properties of Metal Clusters (Springer, 1995).

14. Kelly K. L., Eduardo C., Zhao L .L., Schatz G. C., “The optical properties of metal nanoparticles: the influence of size, shape, and dielectric environment,” J. Chem. Phys. B 107 (3), 668–677 (2003).10.1021/jp026731y [Cross Ref]

15. Prodan E., Nordlander P., “Plasmon hybridization in spherical nanoparticles,” J. Chem. Phys. 120(11), 5444–5454 (2004).10.1063/1.1647518 [PubMed] [Cross Ref]

16. Barchiesi D., Kremer E., Mai V. P., Grosges T., “A Poincaré’s approach for plasmonics: the plasmon localization,” J. Microsc. 229(3), 525–532 (2008).10.1111/j.1365-2818.2008.01938.x [PubMed] [Cross Ref]

17. Barchiesi D., Macias D., Belmar-Letellier L., van Labeke D., Lamy de la Chapelle M., Toury T., Kremer E., Moreau L., Grosges T., “Plasmonics: influence of the intermediate (or stick) layer on the efficiency of sensors,” Appl. Phys. B, Lasers Opt. 93(1), 177–181 (2008).10.1007/s00340-008-3173-5 [Cross Ref]

18. Liu C., Mi C.C., Li B.Q., “Energy Absorption of Gold Nanoshells in Hyperthermia Therapy,” IEEE Transactions on Nanobioscience 7, 206–214 (2008).10.1109/TNB.2008.2002284 [PubMed] [Cross Ref]

19. Loo C., Hirsch L., Lee M. H., Chang E., West J., Halas N., Drezek R., “Gold nanoshell bioconjugates for molecular imaging in living cells,” Opt. Lett. 30(9), 1012–1014 (2005).10.1364/OL.30.001012 [PubMed] [Cross Ref]

20. Grady N. K., Halas N. J., Nordlander P., “Influence of dielectric function properties on the optical response of plasmon resonant metallic nanoparticles,” Chem. Phys. Lett. 399, 167–171 (2004).10.1016/j.cplett.2004.09.154 [Cross Ref]

21. Grosges T., Barchesi D., Toury T., Gréhan G., “Design of nanostructures for imaging and biomedical applications by plasmonic optimization,” Opt. Lett. 33(23), 2812–2814 (2008).10.1364/OL.33.002812 [PubMed] [Cross Ref]

22. Bohren C. F., Huffman D. R., Absorption and Scattering of Light by Small Particles (John Willey & Sons, Inc., 2003).

23. Mie G., “Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen,” Ann. Phys. 25, 377–445 (1908).10.1002/andp.19083300302 [Cross Ref]

24. Gouesbet G., Maheu B., Gréhan G., “Light scattering from a sphere arbitrarily located in a Gaussian beam, using a Bromwich formulation,” J. Opt. Soc. Am. A 5, 1427–1443 (1988).10.1364/JOSAA.5.001427 [Cross Ref]

25. Gréhan C., Gouesbet G., Guilloteau F., “Comparison of the diffraction theory and the generalized lorenz-mie theory for a sphere arbitrarily located into a laser beam,” Opt. Commun. 90, 1–6 (1992).10.1016/0030-4018(92)90315-I [Cross Ref]

26. Du H., “Mie-scattering calculation,” Appl. Opt. 43, 1951–1956 (2004).10.1364/AO.43.001951 [PubMed] [Cross Ref]

27. Grosges T., Vial A., Barchiesi D., “Models of near-field spectroscopic studies: comparison between Finite-Element and Finite-Difference methods,” Opt. Express 13(21), 8483–8497 (2005), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-13-21-848310.1364/OPEX.13.008483 [PubMed] [Cross Ref]

28. Barchiesi D., Guizal B., Grosges T., “Accuracy of local field enhancement models: toward predictive models?,” Appl. Phys. B, Lasers Opt. 84(1–2), 55–60 (2006).

29. Grosges T., Borouchaki H., Barchiesi D., “Improved scheme for accurate computation of high electric near-field gradients,” Opt. Express 15(3), 1307–1321 (2007), http://www.opticsinfobase.org/abstract.cfm?URI=oe-15-3-130710.1364/OE.15.001307 [PubMed] [Cross Ref]

30. Borouchaki H., Grosges T., Barchiesi D., “Improved 3D adaptive remeshing scheme applied in high electromagnetic field gradient computation,” Finite Elem. Anal. Des. 46(1–2), 84–95 (2010).10.1016/j.finel.2009.06.026 [Cross Ref]

31. Palik E. D., Handbook of Optical Constants of Solid I (Academic Press, 1985).

32. Barchiesi D., van Labeke D., “Application of Mie scattering of evanescent waves to scanning optical microscopy theory,” J. Mod. Opt. 40(7), 1239–1254 (1993).10.1080/09500349314551331 [Cross Ref]

33. Macias D., Vial A., Barchiesi D., “Application of evolution strategies for the solution of an inverse problem in near-field optics,” J. Opt. Soc. Am. A 21, 1465–1471 (2004).10.1364/JOSAA.21.001465 [PubMed] [Cross Ref]

34. Schwartzberg A. M., Olson T. Y., Talley C. E., Zhang J. Z., “Synthesis, characterization, and tunable optical properties of hollow gold nanospheres,” J. Phys. Chem. B 110, 19935–19944 (2006).10.1021/jp062136a [PubMed] [Cross Ref]

35. Xu Z. C., Shen C. M., Xiao C. W., Yang T. Z., Zhang H. R., Li J. Q., Gao H. J., “Wet chemical synthesis of gold nanoparticles using silver seeds: a shape control from nanorods to hollow spherical nanoparticles,” Nanotechnology 18, 115608 (2007).10.1088/0957-4484/18/11/115608 [Cross Ref]

36. Barchiesi D., “Adaptive non-uniform, hyper-ellitist evolutionary method for the optimization of plasmonic biosensors,” in Proceedings of International Conference on Computers and Industrial Ingineering (CIE39), IEEE 1, 542–547 (2009).

37. Kessentini S., Barchiesi D., Grosges T., Giraud-Moreau L., Lamy de la Chapelle M., “Adaptive non-uniform particle swarm optimization: application to plasmonic design,” Int. J. Appl. Meta. Comput. 2(1), 18–28 (2011).

38. Tarantola A., Inverse Problem Theory and Methods for Model Parameter Estimation (SIAM, 2005).

39. Duck F. A., Physical Properties of Tissue A Comprehensive Reference Book (Academic Press, 1990).

40. Stoller P., Jacobsen V., Sandoghdar V., “Measurement of the complex dielectric constant of a single gold nanoparticle,” Opt. Lett. 31(16), 2474–2476 (2006).10.1364/OL.31.002474 [PubMed] [Cross Ref]

41. Barnes W. L., “Comparing experiment and theory in plasmonics,” J. Opt. A, Pure Appl. Opt. 11, 114002 (2009).10.1088/1464-4258/11/11/114002 [Cross Ref]

42. Scaffardi L. B., Lester M., Skigin D., Tocho J. O., “Optical extinction spectroscopy used to characterize metallic nanowires,” Nanotechnology 18, 315402 (2007).10.1088/0957-4484/18/31/315402 [Cross Ref]

43. Huang X., El-Sayed M. A., “Gold nanoparticles optical properties and implementations in cancer diagnosis and photothermal therapy,” J. Adv. Res. 1(1), 13–28 (2010).10.1016/j.jare.2010.02.002 [Cross Ref]

Articles from Biomedical Optics Express are provided here courtesy of **Optical Society of America**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |