Research Papers

Fitting the optical constants of gold, silver, chromium, titanium, and aluminum in the visible bandwidth

[+] Author Affiliations
Dominique Barchiesi

University of Technology of Troyes, Institut National de Recherche en Informatique et Automatisme, Automatic Mesh Generation and Advanced Methods (GAMMA3), 12 rue Marie Curie, CS 42060, 10004 Troyes Cedex, France

Thomas Grosges

University of Technology of Troyes, Institut National de Recherche en Informatique et Automatisme, Automatic Mesh Generation and Advanced Methods (GAMMA3), 12 rue Marie Curie, CS 42060, 10004 Troyes Cedex, France

J. Nanophoton. 8(1), 083097 (Jan 06, 2014). doi:10.1117/1.JNP.8.083097
History: Received October 15, 2013; Revised December 7, 2013; Accepted December 9, 2013
Text Size: A A A

Open Access Open Access

Abstract.  The fitting of metal optical properties is a topic that has applications in advanced simulations of spectroscopy, plasmonics, and optical engineering. In particular, the finite difference time domain method (FDTD) requires an analytical model of dispersion that verifies specific conditions to produce a full spectrum in a single run. Combination of Drude and Lorentz models, and Drude and critical points models, are known to be efficient, but the number of parameters to be adjusted for fitting data can prevent accurate results from simulated annealing or Nelder-Mead. The complex number relative permittivities of Au, Ag, Al, Cr, and Ti from either Palik or Johnson and Christy experimental data in the visible domain of wavelengths are successfully fitted by using the result of the particle swarm optimization method with FDTD constraint, as a starting point for the Nelder-Mead method. The results are well positioned compared to those that can be found in the literature. The results can be used directly for numerical simulations in the visible domain. The method can be applied to other materials, such as dielectrics, and to other domain of wavelengths.

Figures in this Article

The fitting of optical constants of metals has given rise to an extensive literature. The optical properties of bulk materials1,2 are commonly used for simulations, but recent experimental data for nanostructured materials will be used in the near future.3,4 Indeed, the surface plasmon resonance (SPR) setups are able to determine thicknesses of multilayers4 and also optical properties.5 Therefore, a versatile and efficient method of fitting of experimental data of the optical properties is of interest.

The historical Drude and Lorentz models (DL) were used6,7 for gold in a wide domain of wavelengths. Expansion of Lorentzian terms (L4) was proposed by Hao and Nordlander8 and tested by comparison between Mie theory and finite time difference domain (FDTD). More recently, the description of gold permittivities by means of critical points model9,10 was proved to be efficient for the modeling of spectroscopy with FDTD methods,11 finite element method,12 and discrete dipole approximation.13 Fitting of the optical constants is also useful if eigenvalues of complex structures are computed14,15 or if resonances are searched in dispersion curves.16,17 For the design and the optimization of nanostructures,18 the accuracy of numerical results depends on the quality of the fitting of the relative permittivities.19,20

The models for fitting the relative permittivity are functions of the angular frequency of illumination ω in a specific range of wavelengths in vacuum λ0. Each model involves parameters considered as degrees of freedom for the fitting of handbooks experimental data from literature. The parameters are searched such that the distance between the reference data1,2 and the model is minimal. The use of these models for FDTD requires filling an additional constraint that is included in the proposed method. The novelty of this paper is in the method that consists of using a particle swarm optimization method (PSO) and then in improving the solution by using a Nelder-Mead (NM) algorithm, using the solution of the PSO algorithm as starting point.

The combinations of DL, and Drude and critical points models (DCP) are shortly described in Sec. 2 as well as the criterion that guarantees the convergence of FDTD algorithms. The principle of the proposed PSO method is given in Sec. 3. The results of the fitting of the relative permittivity of some metals are given and discussed in Sec. 4 before concluding.

Combination of Drude and Lorentz Models

The combination of Drude and Lorentz models (DL) describes both the intraband (Drude model) and interband (Lorentz model) electronic transitions.21 It enables the imaginary part to be a decreasing function of ω on the contrary of the single Drude model. The DL model is efficient in the wavelengths range [500,1000]nm.7 Nevertheless, the best choice of models of fitting as well as their physical parameters depends on the investigated domain of wavelengths. This study is devoted to the fitting in the visible domain of wavelengths, [400,800]nm.

The limitation to a single Lorentz term limits the memory requirement for FDTD, as it increases linearly with the number of terms used for the dispersion law.7 The function of fit ϵDL of the relative permittivity of metal is written as the sum of the Drude and the Lorentz models.6,10,22Display Formula

ϵDL(ω)=ϵωD2ω(ω+iγD)Drude modelΔϵωL2ω2ωL2+iγLωLorentz model.(1)

The pure imaginary number is denoted as i such as i2=1. Frictional forces proportional to the velocity of electrons with viscous damping coefficients γD and γL are introduced and expressed in the same units as ω. The plasma frequency ωL is associated with intraband transitions, Δϵ is the oscillator strength, and ϵ is the relative permittivity for high frequencies. For high frequencies, electrons cannot follow excitation, and the illumination does not see free electrons anymore; therefore ϵ=1.23 Nevertheless, it is commonly admitted that this asymptotic value could also be adjusted for fitting of the relative permittivity in a limited wavelength bandwidth. It is therefore an additional degree of freedom of the fitting. The coefficient ωD is called plasma (or Langmuir) angular frequency (it is not the frequency of volume plasmon21), with Ne the density of electrons that contribute to the optical properties [Eq. (3), it is not the valence band electron concentration21], q the charge of electron, m0 the electron mass, and ϵ0 the relative permittivity of vacuum. Display Formula

ωD=Neq2ϵ0m0.(2)

The volume density of the electrons that contribute to the relative permittivity of gold Ne (m3) is introduced. It can be evaluated from an atomistic model and is actually a function of the density ρm and atomic weight M. Display Formula

Ne=ρmNaMne,(3)
with Na=6.022×1023 the Avogadro constant and ne the supposed number of free electrons in a single atom. This value can be compared to that obtained by the fitting in the following. The atomistic volume density of electron Ne is calculated by using outer electron shell of atoms. Equation (1) could also be written as a product of ωD2 by both Drude and Lorentz terms,21 but the formulation used in 7 is preferred to facilitate the comparison.

The physical parameters of the model are {ϵ,γD,Ne,Δϵ,ωL,γL}. Therefore the dimension of the problem is dim(D)=6 for the fitting.

Combination of Drude and Critical Points Models

The critical points model (CP) was introduced in 1998 by Leng et al.24 for silicon and used for gold.9 This model describes the interband transitions in violet/near-uv region. It satisfies the Kramers-Kronig consistency and is able to reproduce the dispersion of gold with higher accuracy than the previous models. It includes a phase factor (ϕ) and corresponds to first-order poles in the complex plane.9,22,24 The interband transition angular frequency is Ω, and the transition broadening is governed by the damping angular frequency Γ. It is recommended to take two critical points terms.9 This model helps to suppress the nonphysical absorption of the previous fitting in the region of transparency or near-transparency.24 In this study, the CP replaces the Lorentz additional term in Eq. (1), and therefore the model is a combination of Drude and CP (DCP). Display Formula

ϵDCP(ω)=ϵNeq2ϵ0m01ω(ω+iγD)+i=1i=2ΔϵiΩi[exp(iϕi)ΩiωiΓi+exp(iϕi)Ωi+ω+iΓi].(4)

This approach was used by Vial et al.11 for fitting the relative permittivities of gold and silver from Johnson and Christy.2 The CP involves 11 physical parameters: {ϵ,Ne,γD,Γi,Ωi,Δϵi,andϕi} for i=1, 2. The dimension of the problem of fitting is dim(D)=11.

Criterion of Convergence for FDTD

According to Vial et al.,22,25 a criterion must be verified to check the convergence of the FDTD algorithm. Display Formula

C=ϵϵ+χ0<1,(5)
with Display Formula
χ0=(ωDγD)2[1exp(γDΔt)]+ωL2γDΔt+R(iηαiβ{1exp[(α+iβ)Δt]}),(6)
where R(·) is the real part of a complex number.

In the case of combination of Drude and Lorentz model, χ0 is a function of α=γL/2, β=ωL2α2, and η=ΔϵωL2/β. For the CP, χ0 is a function of α=Γi, β=Ωi, and η=2ΔϵiΩiexp(iΦi).25,26 In Eq. (6), Δt is evaluated from the size of the grid Δx used for FDTD: Δt=Δx/(2c), with c the speed of light in vacuum. In the following, we consider Δx=1nm.25

The criterion in Eq. (5) is a constraint for the method of fitting. Therefore, the random generation of particles x [Eq. (11)] is stressed by this criterion, and the solution given by the NM algorithm is selected among those satisfying this criterion.

The problem of data fitting belongs to the same class as the resolution of the inverse problem and as the optimization of systems. Classical methods were used for the fitting of relative permittivities: the NM method and the simulated annealing.7 The NM method27 searches solutions in an unbounded domain. This fact could be a drawback as it could lead to nonphysical solutions or divergence of the algorithm. The NM method was used for problems with small degree of freedom [dim(D)].28 The simulated annealing subroutine belongs to the same class as the Monte Carlo method (heuristic methods) and is fully described by Kirkpatrick et al.29 and Tarantola.30 This method can be bounded or not. Both methods require a starting point for the algorithm initialization and the solution depends on it. The proposed PSO method enables to determine the starting point of the NM method within a space of parameters that ensures the criterion of convergence for FDTD [Eq. (5)].

Among heuristic bounded methods without starting point, the PSO method is widely used for optical applications.3133 The PSO method was proposed in 1995 by Kennedy and Eberhart.34 It was used for the optimization of gold nanoshells using discrete dipole approximation (DDA) and Mie models for photothermal therapy,35 and more general problems of plasmonics including the optimization of nanostructures and the resolution of the inverse problem.5,36,37

The PSO mimics the behavior of a swarm of particles moving to a potential well with analogy to the bees swarming in search of pollen. The nearby of particle with respect to the bottom of the potential well is evaluated through a fitness function F. In the present case, the fitness function F is the standard deviation of the data computed from model ϵM to the Nref values of the relative permittivity of reference ϵr. Display Formula

F(x)=1Nref1Nref|ϵM(x)ϵr|2.(7)

The fitness function F(x) depends on the physical parameters x that are necessary to compute the model ϵM for all ω. The models for fitting ϵM(x) are described in Sec. 2. In the following, ϵM(x)=ϵDL [Eq. (1)] or ϵM(x)=ϵDCP [Eq. (4)]. According to this definition, for all physical parameters set x, the absolute error of fitting can be defined for the real part and the imaginary part of ϵr(ω). Display Formula

σR[ϵM(ω)]=1NrefR[ϵM(ω)ϵr(ω)]2,(8)
Display Formula
σI[ϵM(ω)]=1NrefI[ϵM(ω)ϵr(ω)]2.(9)

If particles reach the position of the minimum of potential, then the fitness function is minimum; therefore a good position for a particle corresponds to a small value of the fitness function and to a good parameters set of the model. All inputs of the model (physical parameters) form a vector x(t), whose size is the degree of freedom of the model. The simulated cinematic problem for particles depends on a virtual time t, which is the iteration step of the algorithm. Each particle moves in the space of physical parameters following Eq. (10). Display Formula

x(t+1)=x(t)+V(t+1)[(t+1)t]=x(t)+V(t+1).(10)

The velocity V(t+1) is computed from the best position p(t) of each particle over previous generations up to step t and from the swarm global best g(t).38Display Formula

V(t+1)=ωPSOV(t)+U1c1[p(t)x(t)]+U2c2[g(t)x(t)].(11)

The inertia weight ωPSO linearly decreases from 0.9 to 0.4.39 The random gene Ui(i=1,2)] always checking the condition given in Eq. (5). The acceleration coefficients are c1=c2=2. If the speed is too high, the particle leaves the domain of search. In this case, it is regenerated randomly in this domain. In this study, the number of particles in the swarm is N=30: at each step t, N particles are moving in the space of search. In this basic version of the algorithm, the stop criterion is a maximum number of iterations of PSO (1000). We checked that the results are hardly dependent on all these exogenous parameters (ωPSO, c1, c2, and N). In 37, this property and the efficiency of this method was verified by using evolutionary method based on evolution strategies.4042

Then the solution obtained with PSO is used as starting point of an NM algorithm and the success of the method is revealed by the convergence rate to the same best solution. Only solutions that satisfy the condition in Eq. (5) are kept. A hundred realizations of the same algorithm help to check the stability of the solution. If the algorithm gives some identical solutions along realizations, the best solution is assumed to be determined. In the investigated cases, each realization takes a few minutes on a personal computer. This property is directly linked to the simplicity of the models that are detailed in Sec. 2.

We first validate the PSO method by comparison of results with those given in the literature. Then we apply the method to the fitting of gold, silver, chromium, aluminum, and titanium in the domain of visible wavelengths. The efficiency of the method is tested by comparisons of our results with published data,7,43 where the domain of wavelengths was [500,1000]nm. Of course the influence of the domain of wavelengths for the fitting is crucial to the results.

In this study, the domain of wavelengths is [400,800]nm and the domain of search for the physical parameters is D: {γD,γL}b[1013,1018](rad.s1), ϵ[1000,1000], Δϵ[0,1000], Neb[1027,1030](m3) [Eq. (3)]. The fitting of bulk data for gold, silver, chromium, aluminum, and titanium are, respectively, investigated.

Gold

Gold has mechanical properties and an ease of implementation, which puts it at the forefront of materials used in plasmonics. Moreover, it cannot be oxidized by the common surrounding medium. Table 1 gives the best physical parameters obtained from the proposed method for the fitting of Johnson and Christy data.2Figure 1 shows the results of fitting and the absolute errors σR and σI [Eqs. (8) and (9)].

Table Grahic Jump Location
Table 1Fitting of Johnson and Christy data (ϵr) for gold by the combination of Drude and Lorentz (DL) models [Eq. (1)] and critical points (DCP) model. F is the fitness function, C is the finite difference time domain (FDTD) criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Graphic Jump LocationF1 :

(a) Real and imaginary part of the experimental data for gold2 and fitting with the Drude-critical points models (DCP) (Table 1). (b) The error between the fitting and the reference data. The fitting with parameters from 22 (ϵV) and 43 (ϵV) are plotted.

The absolute error for the fittings using Eqs. (8) and (9) are compared to those in the literature.

  • With DL, σR(ϵV)=4.5087 and σI(ϵV)=3.6146.7 The proposed method gives σR=1.3292 and σI=1.2208 (Table 1).
  • With DCP, σR(ϵV)=0.9417, σI(ϵV)=1.0606,22σR(ϵV)=0.7882, and σI(ϵV)=1.2268.43 The proposed method gives σR=0.57097 and σI=1.0747 (Table 1).

These results show that the proposed method is competitive with the simulated annealing,7,22,43 leaving to smaller errors in the visible range. For DL model, the proposed method produces absolute errors which are less than one third of those found in the literature.

The fitness function for DL is three times that obtained for DCP. The DCP model is therefore more accurate to describe the optical properties of gold in the visible range.

The Palik data for gold come from two series of measurements, leading to a change of slope of the curve. These data are therefore suitable to test the method, by comparison with the fitting of Johnson and Christy data. Table 2 gives the best physical parameters. Figure 2 shows the results of fitting and the absolute errors σR and σI [Eqs. (8) and (9)]. The results of fitting for DL in Table 1 (Johnson and Christy) are close to those in Table 2 (Palik), on the contrary of DCP fitting. Figure 2 shows that DCP is highly sensitive to the change of slope of the real part of ϵr. The result of fitting for the imaginary part is therefore affected by the change of slope (dashed curve).

Table Grahic Jump Location
Table 2Fitting of Palik data (ϵr) for gold by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between between imaginary parts of fitting and data.
Graphic Jump LocationF2 :

(a) Real and imaginary part of the experimental data for gold1 and fitting with the Drude-Lorentz (DL) and the DCP (Table 2). (b) The error between the fitting and the reference data.

The accuracy of both models given by the fitness function F is less influenced by the chosen model (DL or DCP) than for the fitting of Johnson and Christy data. According to the previous results, the method can be applied to other metals with confidence.

Silver

Silver is known to be more efficient than gold for plasmonic nanostructures, although silver suffers from oxidization and changes of optical properties when nanostructured.44Figure 3 shows both fittings and the error of fitting σ(ω). The real and imaginary parts of ϵV22 are shown for comparison with those obtained from the proposed method (combination of PSO and NM algorithms). The two errors of fitting are about the same order of magnitude along ω but exhibit local differences between the reference 22 and the PSO results.

Graphic Jump LocationF3 :

(a) Real and imaginary part of the experimental data for silver1 and fitting with the DL and the DCP (Table 3). (b) The error between the fitting and the reference data. The fitting with parameters from 22 (ϵV) and 43 (ϵV) are plotted.

The optical properties of silver are known to have smooth behavior in the visible domain. Therefore, the fitness function, the absolute errors of fitting, and the angular frequencies ωD are close together for both models. The amplitudes Δϵi of critical points contributions are larger than for gold.

The absolute errors for the fitting with DCP, using Eqs. (8) and (9), are compared to those in the literature for CP model: σR(ϵV)=0.66287, σI(ϵV)=0.35788,22σR(ϵV)=0.81118, and σI(ϵV)=0.70911.43 They are roughly two times the values in Table 3: σR=0.36891 and σI=0.34168. The parameters for the DL model in Table 3 also give a better fit than that in 43. The comparison with results of fitting from Refs. 22 and 43 shows a great improvement of the quality of the solution, which is mainly due to a different interval of wavelengths for fitting.

Table Grahic Jump Location
Table 3Fitting of Palik data (ϵr) for silver by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Aluminum

Aluminum is a good candidate to plasmon resonance as the real part of its relative permittivity is <23 in the visible domain, even if this material can be easily oxidized.37Table 4 gives the results of fitting for aluminum. These values are used to plot Fig. 4. Figure 4(b) illustrates the absolute error on the fitting of the real part and of the imaginary part. The comparison with results of fitting with DCP model (ϵV) from Refs. 22 and 43 are shown for comparison with those obtained from the PSO method. The errors of fitting between the reference 22,43 and the proposed fitting are about the same order of magnitude along ω except for wavelengths >630nm. The fitting with DL is much less efficient than DCP as DL cannot handle the change of slope of the real part of ϵr [see Fig. 4(a)].

Table Grahic Jump Location
Table 4Fitting of Palik data (ϵr) for aluminum by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Graphic Jump LocationF4 :

(a) Real and imaginary part of the experimental data for aluminum1 and fitting with the DL and the CP (Table 4). (b) The error between the fitting and the reference data. The fitting with parameters from Refs. 22 and 43 (ϵV) are plotted.

Both results of fitting show that the Drude model is not suitable to describe the behavior of the optical properties of aluminum. Consequently, the higher degree of freedom of DCP makes it more efficient for the fitting. From Refs. 22 and 43, σR(ϵV)=2.664 and σI(ϵV)=2.178, where the interval of wavelengths was [400,1000]nm. In Table 4, the corresponding values are σR=1.5629 and σI(ϵV)=1.1941. The increase of absolute accuracy is around 2.

Chromium

Chromium is often used in nanotechnologies as adhesion layer for gold on dielectric substrates. Table 5 gives the results of fitting of Palik data for chromium. The corresponding dispersion curve is plotted in Fig. 5. Figure 5(b) illustrates the absolute error on the fitting of the real part and of the imaginary part. The comparison with results of fitting with DCP model (ϵV) from Refs. 22 and 43 are shown for comparison with those obtained from the PSO method. The errors of fitting are about the same order of magnitude along ω but exhibit local differences between the reference 22 and the PSO results.

Table Grahic Jump Location
Table 5Fitting of Palik data (ϵr) for chromium by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Graphic Jump LocationF5 :

(a) Real and imaginary part of the experimental data for chromium1 and fitting with the DL and the DCP (Table 5). (b) The error between the fitting and the reference data. The fitting with parameters from Refs. 22 and 43 (ϵV) are plotted.

The optical properties of a nanometric layer of chromium are known to influence the plasmon resonance of SPR.45 Therefore in the visible domain of wavelengths, the absolute errors for the fitting with DCP are significant. From Refs. 22 and 43, σR(ϵV)=0.93846 and σI(ϵV)=0.58162, where the interval of wavelengths was [400,1000]nm. In Table 5, the corresponding values are σR=0.29314 and σI(ϵV)=0.48649. The comparison with results of fitting from Refs. 22 and 43 shows a great improvement of the quality of the solution as shown in Fig. 5(b). Consequently, extending the domain to far-infrared domain is a drawback for the accurate fitting in the visible domain. The fitting of the optical properties of chromium must be adapted to the domain of wavelengths of interest.

Titanium

Titanium is used as an alternative medium for adhesion layer of gold on dielectric samples.45Table 6 gives the results of fitting for titanium. These values are used to plot Fig. 6. Figure 6(b) shows the absolute error on the fitting of the real part and of the imaginary part. The fitting with DL is much less efficient than DCP as DL cannot handle the change of slope of the real part of ϵr [see Fig. 6(a)]. On the contrary, the imaginary part of the relative permittivity is well described by both models.

Table Grahic Jump Location
Table 6Fitting of Johnson and Christy data (ϵr) for titanium by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Graphic Jump LocationF6 :

(a) Real and imaginary part of the experimental data for titanium2 and fitting with the DL and the DCP (Table 6). (b) The error between the fitting and the reference data.

Both results of fitting show that the Drude model is not suitable to describe the behavior of the optical properties of titanium (Table 6).

Concluding Remarks

Both models of dispersion can be used according to the value of the fitness function F (Tables 2 to 6), but Figs. 1 to 6 show local discrepancies of data fitting. In the general case, the error of fitting is spread over the real part or the imaginary part of the relative permittivity. The errors are globally lower compared to that found in the literature, especially for the DCP method. The DL model involves two terms and is less efficient than the CP model using three terms. Therefore, the DL model could be chosen if sparing memory is required for FDTD calculations. However, the DCP model is more accurate for all investigated metals. Aluminum’s optical properties are not well described by the Drude model. Depending on metals, both models are not equivalent, and the inspection of results allows to draw conclusions on the use of models that are summarized in Table 7. The performance of the method is better if the number of “+” increases. This performance is based on the value of the fitness function.

Table Grahic Jump Location
Table 7Performance of fitting of ϵr by the combination of DL models and DCP model [Eq. (4)] with the proposed method (particle swarm optimization+Nelder Mead), in the domain of visible wavelengths [400,800]nm. (P) indicates data from Palik1 and (JC) from Johnson and Christy.2

The results of fitting of various metals in the visible bandwidth provide the following remarks:

  • The DL model may be considered as competitive with the DCP model in terms of accuracy. The choice between these models for FDTD is based on a compromise between memory requirement that increases linearly with the number of terms and the required accuracy of numerical results.
  • The DCP model is able to reproduce the change of slope of the real part of the relative permittivity on the contrary of the DL model.
  • The FDTD parameter Δx=1nm is a strong constraint that could be released easily for applications where the accuracy on the near-field computation is not necessary. The proposed method could generate better fittings in this case.

The particle swarm method under constraint followed with an NM algorithm have been successfully used for the fitting of the relative permittivity of gold, silver, chromium, aluminum, and titanium in the visible wavelength domain. The PSO helps to find a good starting point for the NM method, and the generation of physical parameters that are compatible with their use for FDTD is ensured. The proposed method seems to be efficient for the hard problem of fitting of the relative permittivities of metals as shown by comparison with results in the literature.

Two models are involved for the fitting of the relative permittivities: combinations of the Drude model with Lorentz (DL) and with critical points models (DCP). The efficiency of both models for fitting is compared for each metal. The fitting coefficients of optical properties of metals that are currently used in optical design are given. The use of these results could be particularly useful in plasmonics and design of nanostructured biosensors, but beware of the crude use of these values: the domain of validity of these properties can change dramatically when the size of the metal particles is less than a few tens of nanometers.44,46

The results of fitting can be used directly for any spectroscopic simulation and especially in FDTD codes.8,47 These results are complementary to those found in the literature where more than two DL terms are used for fitting. This method could also be applied to absorbing dielectrics, especially for nanowires and nanotube studies, with other models of fitting or other domains of wavelengths.4850 The advantage of such a heuristic method lies in its applicability to various problems of fitting, the optimization of complex systems in engineering,13,35 the propagation of uncertainties, or the tolerance study of models.18,28 The use of more complex models 15,16,51,52 can extend the domain of its application.

The authors thank the Agence Nationale de la Recherche for financial support (ANR-2011-NANO-008 NANOMORPH).

Palik  E. D., Handbook of Optical Constants. ,  Academic Press Inc. ,  San Diego  (1985).
Johnson  P., Christy  R. W., “Optical constants of the noble metals,” Phys. Rev. B. 6, (12 ), 4370 –4379 (1972), CrossRef. 0163-1829 
Flory  F., Escoubas  L., Berginc  G., “Optical properties of nanostructured materials: a review,” J. Nanophotonics. 5, (1 ), 052502  (2011), CrossRef. 1934-2608 
El-Gohary  S. H., Kim  N.-H., Byun  K. M., “Optical determination of thick graphene layer number based on surface plasmon resonance,” J. Nanophotonics. 7, (1 ), 073799  (2013), CrossRef. 1934-2608 
Salvi  J., Barchiesi  D., “Measurement of thicknesses and optical properties of thin films from surface plasmon resonance (SPR),” Appl. Phys. B—Lasers Opt.. 114,  (2014), CrossRef. 0946-2171 
Rakić  A. D. et al., “Optical properties of metallic films for vertical-cavity optoelectronic devices,” Appl. Opt.. 37, (22 ), 5271 –5283 (1998), CrossRef. 0003-6935 
Vial  A. et al., “Improved analytical fit of gold dispersion: application to the modeling of extinction spectra with a finite-difference time-domain method,” Phys. Rev. B. 71, (8 ), 085416 –085423 (2005), CrossRef. 0163-1829 
Hao  F., Nordlander  P., “Efficient dielectric function for FDTD simulation of the optical properties of silver and gold nanoparticles,” Chem. Phys. Lett.. 446, (1–3 ), 115 –118 (2007), CrossRef. 0009-2614 
Etchegoin  P. G., Ru  E. C. L., Meyer  M., “An analytic model for the optical properties of gold,” J. Chem. Phys.. 125, , 164705  (2006), CrossRef. [Erratum in J. Chem. Phys. 127, 189901 (2007)]. 0021-9606 
Rodrigo  S. G., García-Vidal  F. J., Martín-Moreno  L., “Influence of material properties on extraordinary optical transmission through hole arrays,” Phys. Rev. B. 77, (7 ), 075401  (2008), CrossRef. 0163-1829 
Vial  A., Laroche  T., “Comparison of gold and silver dispersion laws suitable for FDTD simulations,” Appl. Phys. B—Lasers Opt.. 93, (1 ), 139 –143 (2008), CrossRef. 0946-2171 
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), CrossRef. 1094-4087 
Barchiesi  D. et al., “Localized surface plasmon resonance in arrays of nano-gold cylinders: inverse problem and propagation of uncertainties,” Opt. Express. 21, (2 ), 2245 –2262 (2013), CrossRef. 1094-4087 
Guizal  B., Yala  H., Felbacq  D., “Reformulation of the eigenvalue problem in the Fourier modal method with spatial adaptive resolution,” Opt. Lett.. 34, (18 ), 2790 –2792 (2009), CrossRef. 0146-9592 
Guizal  B., Barchiesi  D., Felbacq  D., “Electromagnetic beam diffraction by a finite lamellar structure,” J. Opt. Soc. Am. A. 20, (12 ), 2274 –2280 (2003), CrossRef. 0740-3232 
Barchiesi  D. et al., “A Poincaré’s approach for plasmonics: the plasmon localization,” J. Microsc.. 229, (3 ), 525 –532 (2008), CrossRef. 0022-2720 
Castanié  A., Felbacq  D., “Confined plasmonic modes in a nanocavity,” Opt. Commun.. 285, (16 ), 3353 –3357 (2012), CrossRef. 0030-4018 
Grosges  T. et al., “Design of nanostructures for imaging and biomedical applications by plasmonic optimization,” Opt. Lett.. 33, (23 ), 2812 –2814 (2008), CrossRef. 0146-9592 
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), CrossRef. 0946-2171 
Barchiesi  D. et al., “Plasmonics: influence of the intermediate (or stick) layer on the efficiency of sensors,” Appl. Phys. B—Lasers Opt.. 93, (1 ), 177 –181 (2008), CrossRef. 0946-2171 
Rakić  A. D., “Algorithm for the determination of intrinsic optical constants of metal films: application to aluminum,” Appl. Opt.. 34, (22 ), 4755 –4767 (1995), CrossRef. 0003-6935 
Vial  A., Laroche  T., “Description of dispersion of metals by means of the critical points model and application to the study of resonant structures using the FDTD method,” J. Phys. D: Appl. Phys.. 40, (22 ), 7152 –7158 (2007), CrossRef. 0022-3727 
Polyanskiy  M. N., “Refractive index database,” http://refractiveindex.info://refractiveindex.info (10  October 2013).
Leng  J. et al., “Analytic representations of the dielectric functions of materials for device and structural modeling,” Thin Solid Films. 313–314, , 132 –136 (1998), CrossRef. 0040-6090 
Vial  A., “Implementation of the critical points model in the recursive convolution method for modelling dispersive media with the finite-difference time domain method,” J. Opt. A: Pure Appl. Opt.. 9, (7 ), 745 –748 (2007), CrossRef. 1464-4258 
Udagedara  I. B., Rukhlenko  I. D., Premaratne  M., “Surface plasmon-polariton propagation in piecewise linear chains of composite nanospheres: the role of optical gain and chain layout,” Opt. Express. 19, (21 ), 19973 –19986 (2011), CrossRef. 1094-4087 
Nelder  J., Mead  R., “A simplex method for function minimization,” Comput. J.. 7, (4 ), 308 –313 (1965), CrossRef. 0010-4620 
Grosges  T. et al., “Nanoshells for photothermal therapy: a Monte-Carlo based numerical study of their design tolerance,” Biomed. Opt. Express. 2, (6 ), 1584 –1596 (2011), CrossRef. 2156-7085 
Kirkpatrick  S., Gelatt  C. D., Vecchi  M. P., “Optimization by simulated annealing,” Science. 220, (4598 ), 671 –680 (1983), CrossRef. 0036-8075 
Tarantola  A., Inverse Problem Theory and Methods for Model Parameter Estimation. ,  Society for Industrial and Applied Mathematics (SIAM) ,  Philadelphia  (2005).
Rodriguez-Cobo  L., Cobo  A., Lopez-Higuera  J. M., “Recovering a fiber Bragg grating axial strain distribution from its reflection spectrum,” Opt. Lett.. 38, (13 ), 2327 –2329 (2013), CrossRef. 0146-9592 
Yang  C. et al., “Design of reflective color filters with high angular tolerance by particle swarm optimization method,” Opt. Express. 21, (8 ), 9315 –9323 (2013), CrossRef. 1094-4087 
Macías  D. et al., “Retrieval of relevant parameters of natural multilayer systems by means of bio-inspired optimization strategies,” Appl. Opt.. 52, (11 ), 2511 –2520 (2013), CrossRef. 0003-6935 
Kennedy  J., Eberhart  R., “Particle swarm optimization,” in  IEEE Int. Conf. on Neural Networks , Vol. 4, pp. 1942 –1948,  IEEE ,  New York  (1995).
Kessentini  S., Barchiesi  D., “Quantitative comparison of optimized nanorods, nanoshells and hollow nanospheres for photothermal therapy,” Biomed. Opt. Express. 3, (3 ), 590 –604 (2012), CrossRef. 2156-7085 
Kessentini  S. et al., “Selective and collaborative optimization methods for plasmonics: a comparison,” PIERS Online. 7, (3 ), 291 –295 (2011), CrossRef. 1931-7360 
Barchiesi  D., “Numerical retrieval of thin aluminium layer properties from SPR experimental data,” Opt. Express. 20, (8 ), 9064 –9078 (2012), CrossRef. 1094-4087 
Shi  Y., Eberhart  R. C., “A modified particle swarm optimizer,” in  Proc. IEEE Congress on Evolutionary Computation , pp. 69 –73,  IEEE ,  New York  (1998).
Kessentini  S. et al., “Adaptive non-uniform particle swarm application to plasmonic design,” IJAMC. 2, (1 ), 18 –28 (2011), CrossRef.
Schwefel  H. P., Evolution and Optimum Seeking. ,  John Wiley & Sons Inc. ,  New York  (1995).
Macías  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, (8 ), 1465 –1471 (2004), CrossRef. 0740-3232 
Barchiesi  D., “Adaptive non-uniform, hyper-ellitist evolutionary method for the optimization of plasmonic biosensors,” in  Proc. Int. Conf. Computers & Industrial Engineering , pp. 542 –547,  IEEE ,  Troyes  (2009).
Vial  A. et al., “A new model of dispersion for metals leading to a more accurate modeling of plasmonic structures using the FDTD method,” Appl. Phys. A: Mater.. 103, (3 ), 849 –853 (2011), CrossRef. 0947-8396 
Krenn  J. R. et al., “Near-field optical investigations on nanometric silver particles,” in Photons and Local Probes. , Marti  O., Möller  R., Eds., NATO ASI Series E: Applied Sciences , Vol. 300, pp. 181 –187,  Kluwer Academic Publishers ,  Dordrecht  (1995).
Barchiesi  D., “Improved method based on S matrix for the optimization of SPR biosensors,” Opt. Commun.. 286, (1 ), 23 –29 (2013), CrossRef. 0030-4018 
Wautelet  M., Duvivier  D., “The characteristic dimensions of the nanoworld,” Eur. J. Phys.. 28, (5 ), 953 –959 (2007), CrossRef. 0143-0807 
Prokopidis  K. P., Zografopoulos  D. C., “A unified FDTD/PML scheme based on critical points for accurate studies of plasmonic structures,” J. Lightwave Technol.. 31, (15 ), 2467 –2476 (2013), CrossRef. 0733-8724 
Shultz  L. G., “The optical constants of silver, gold, copper and aluminum. I. The absorption coefficient k,” J. Opt. Soc. Am.. 44, (5 ), 357 –362 (1954), CrossRef. 0030-3941 
Shultz  L. G., Tangherlini  F. R., “The optical constants of silver, gold, copper and aluminum. II. The index of refraction n,” J. Opt. Soc. Am.. 44, (5 ), 362 –368 (1954), CrossRef. 0030-3941 
Chaari  A. et al., “Nanobubble evolution around nanowire in liquid,” Opt. Express. 21, (22 ), 26942 –26954 (2013), CrossRef. 1094-4087 
Fikri  R., Grosges  T., Barchiesi  D., “Apertureless scanning near-field optical microscopy: on the need of the tip vibration modelling,” Opt. Lett.. 28, (22 ), 2147 –2149 (2003), CrossRef. 0146-9592 
Fikri  R., Grosges  T., Barchiesi  D., “Apertureless scanning near-field optical microscopy: numerical modeling of the lock-in detection,” Opt. Commun.. 232, (1–6 ), 15 –23 (2004), CrossRef. 0030-4018 

Dominique Barchiesi is a full professor in theoretical physics, applied mathematics, and statistics at the University of Technology of Troyes. He conducts active research in numerical modeling, optimization, and advanced methods for multiphysics, with application to engineering in nanotechnologies, biotechnologies, and plasmonics. He is the author of over 200 articles, conferences, and book chapters since 1992, in the fields of optics, electromagnetism, coupled problems, signal processing, operational research, and cryptography.

Thomas Grosges is a full professor in numerical and theoretical physics, and applied mathematics at the University of Technology of Troyes. He conducts active research on modeling the interaction of matter and light and on the optimization of numerical solvers in computer sciences and their applications (multiphysics, cryptography). He is the author of more than 50 publications in physics, applied physics, engineering, and computer science academic journals and books.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Dominique Barchiesi and Thomas Grosges
"Fitting the optical constants of gold, silver, chromium, titanium, and aluminum in the visible bandwidth", J. Nanophoton. 8(1), 083097 (Jan 06, 2014). ; http://dx.doi.org/10.1117/1.JNP.8.083097


Figures

Graphic Jump LocationF1 :

(a) Real and imaginary part of the experimental data for gold2 and fitting with the Drude-critical points models (DCP) (Table 1). (b) The error between the fitting and the reference data. The fitting with parameters from 22 (ϵV) and 43 (ϵV) are plotted.

Graphic Jump LocationF2 :

(a) Real and imaginary part of the experimental data for gold1 and fitting with the Drude-Lorentz (DL) and the DCP (Table 2). (b) The error between the fitting and the reference data.

Graphic Jump LocationF3 :

(a) Real and imaginary part of the experimental data for silver1 and fitting with the DL and the DCP (Table 3). (b) The error between the fitting and the reference data. The fitting with parameters from 22 (ϵV) and 43 (ϵV) are plotted.

Graphic Jump LocationF4 :

(a) Real and imaginary part of the experimental data for aluminum1 and fitting with the DL and the CP (Table 4). (b) The error between the fitting and the reference data. The fitting with parameters from Refs. 22 and 43 (ϵV) are plotted.

Graphic Jump LocationF5 :

(a) Real and imaginary part of the experimental data for chromium1 and fitting with the DL and the DCP (Table 5). (b) The error between the fitting and the reference data. The fitting with parameters from Refs. 22 and 43 (ϵV) are plotted.

Graphic Jump LocationF6 :

(a) Real and imaginary part of the experimental data for titanium2 and fitting with the DL and the DCP (Table 6). (b) The error between the fitting and the reference data.

Tables

Table Grahic Jump Location
Table 1Fitting of Johnson and Christy data (ϵr) for gold by the combination of Drude and Lorentz (DL) models [Eq. (1)] and critical points (DCP) model. F is the fitness function, C is the finite difference time domain (FDTD) criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Table Grahic Jump Location
Table 2Fitting of Palik data (ϵr) for gold by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between between imaginary parts of fitting and data.
Table Grahic Jump Location
Table 3Fitting of Palik data (ϵr) for silver by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Table Grahic Jump Location
Table 4Fitting of Palik data (ϵr) for aluminum by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Table Grahic Jump Location
Table 5Fitting of Palik data (ϵr) for chromium by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Table Grahic Jump Location
Table 6Fitting of Johnson and Christy data (ϵr) for titanium by the combination of DL models [Eq. (1)] and DCP model. F is the fitness function, C is the FDTD criterion, σR is the absolute error between real parts of fitting and data, and σI is the absolute error between imaginary parts of fitting and data.
Table Grahic Jump Location
Table 7Performance of fitting of ϵr by the combination of DL models and DCP model [Eq. (4)] with the proposed method (particle swarm optimization+Nelder Mead), in the domain of visible wavelengths [400,800]nm. (P) indicates data from Palik1 and (JC) from Johnson and Christy.2

References

Palik  E. D., Handbook of Optical Constants. ,  Academic Press Inc. ,  San Diego  (1985).
Johnson  P., Christy  R. W., “Optical constants of the noble metals,” Phys. Rev. B. 6, (12 ), 4370 –4379 (1972), CrossRef. 0163-1829 
Flory  F., Escoubas  L., Berginc  G., “Optical properties of nanostructured materials: a review,” J. Nanophotonics. 5, (1 ), 052502  (2011), CrossRef. 1934-2608 
El-Gohary  S. H., Kim  N.-H., Byun  K. M., “Optical determination of thick graphene layer number based on surface plasmon resonance,” J. Nanophotonics. 7, (1 ), 073799  (2013), CrossRef. 1934-2608 
Salvi  J., Barchiesi  D., “Measurement of thicknesses and optical properties of thin films from surface plasmon resonance (SPR),” Appl. Phys. B—Lasers Opt.. 114,  (2014), CrossRef. 0946-2171 
Rakić  A. D. et al., “Optical properties of metallic films for vertical-cavity optoelectronic devices,” Appl. Opt.. 37, (22 ), 5271 –5283 (1998), CrossRef. 0003-6935 
Vial  A. et al., “Improved analytical fit of gold dispersion: application to the modeling of extinction spectra with a finite-difference time-domain method,” Phys. Rev. B. 71, (8 ), 085416 –085423 (2005), CrossRef. 0163-1829 
Hao  F., Nordlander  P., “Efficient dielectric function for FDTD simulation of the optical properties of silver and gold nanoparticles,” Chem. Phys. Lett.. 446, (1–3 ), 115 –118 (2007), CrossRef. 0009-2614 
Etchegoin  P. G., Ru  E. C. L., Meyer  M., “An analytic model for the optical properties of gold,” J. Chem. Phys.. 125, , 164705  (2006), CrossRef. [Erratum in J. Chem. Phys. 127, 189901 (2007)]. 0021-9606 
Rodrigo  S. G., García-Vidal  F. J., Martín-Moreno  L., “Influence of material properties on extraordinary optical transmission through hole arrays,” Phys. Rev. B. 77, (7 ), 075401  (2008), CrossRef. 0163-1829 
Vial  A., Laroche  T., “Comparison of gold and silver dispersion laws suitable for FDTD simulations,” Appl. Phys. B—Lasers Opt.. 93, (1 ), 139 –143 (2008), CrossRef. 0946-2171 
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), CrossRef. 1094-4087 
Barchiesi  D. et al., “Localized surface plasmon resonance in arrays of nano-gold cylinders: inverse problem and propagation of uncertainties,” Opt. Express. 21, (2 ), 2245 –2262 (2013), CrossRef. 1094-4087 
Guizal  B., Yala  H., Felbacq  D., “Reformulation of the eigenvalue problem in the Fourier modal method with spatial adaptive resolution,” Opt. Lett.. 34, (18 ), 2790 –2792 (2009), CrossRef. 0146-9592 
Guizal  B., Barchiesi  D., Felbacq  D., “Electromagnetic beam diffraction by a finite lamellar structure,” J. Opt. Soc. Am. A. 20, (12 ), 2274 –2280 (2003), CrossRef. 0740-3232 
Barchiesi  D. et al., “A Poincaré’s approach for plasmonics: the plasmon localization,” J. Microsc.. 229, (3 ), 525 –532 (2008), CrossRef. 0022-2720 
Castanié  A., Felbacq  D., “Confined plasmonic modes in a nanocavity,” Opt. Commun.. 285, (16 ), 3353 –3357 (2012), CrossRef. 0030-4018 
Grosges  T. et al., “Design of nanostructures for imaging and biomedical applications by plasmonic optimization,” Opt. Lett.. 33, (23 ), 2812 –2814 (2008), CrossRef. 0146-9592 
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), CrossRef. 0946-2171 
Barchiesi  D. et al., “Plasmonics: influence of the intermediate (or stick) layer on the efficiency of sensors,” Appl. Phys. B—Lasers Opt.. 93, (1 ), 177 –181 (2008), CrossRef. 0946-2171 
Rakić  A. D., “Algorithm for the determination of intrinsic optical constants of metal films: application to aluminum,” Appl. Opt.. 34, (22 ), 4755 –4767 (1995), CrossRef. 0003-6935 
Vial  A., Laroche  T., “Description of dispersion of metals by means of the critical points model and application to the study of resonant structures using the FDTD method,” J. Phys. D: Appl. Phys.. 40, (22 ), 7152 –7158 (2007), CrossRef. 0022-3727 
Polyanskiy  M. N., “Refractive index database,” http://refractiveindex.info://refractiveindex.info (10  October 2013).
Leng  J. et al., “Analytic representations of the dielectric functions of materials for device and structural modeling,” Thin Solid Films. 313–314, , 132 –136 (1998), CrossRef. 0040-6090 
Vial  A., “Implementation of the critical points model in the recursive convolution method for modelling dispersive media with the finite-difference time domain method,” J. Opt. A: Pure Appl. Opt.. 9, (7 ), 745 –748 (2007), CrossRef. 1464-4258 
Udagedara  I. B., Rukhlenko  I. D., Premaratne  M., “Surface plasmon-polariton propagation in piecewise linear chains of composite nanospheres: the role of optical gain and chain layout,” Opt. Express. 19, (21 ), 19973 –19986 (2011), CrossRef. 1094-4087 
Nelder  J., Mead  R., “A simplex method for function minimization,” Comput. J.. 7, (4 ), 308 –313 (1965), CrossRef. 0010-4620 
Grosges  T. et al., “Nanoshells for photothermal therapy: a Monte-Carlo based numerical study of their design tolerance,” Biomed. Opt. Express. 2, (6 ), 1584 –1596 (2011), CrossRef. 2156-7085 
Kirkpatrick  S., Gelatt  C. D., Vecchi  M. P., “Optimization by simulated annealing,” Science. 220, (4598 ), 671 –680 (1983), CrossRef. 0036-8075 
Tarantola  A., Inverse Problem Theory and Methods for Model Parameter Estimation. ,  Society for Industrial and Applied Mathematics (SIAM) ,  Philadelphia  (2005).
Rodriguez-Cobo  L., Cobo  A., Lopez-Higuera  J. M., “Recovering a fiber Bragg grating axial strain distribution from its reflection spectrum,” Opt. Lett.. 38, (13 ), 2327 –2329 (2013), CrossRef. 0146-9592 
Yang  C. et al., “Design of reflective color filters with high angular tolerance by particle swarm optimization method,” Opt. Express. 21, (8 ), 9315 –9323 (2013), CrossRef. 1094-4087 
Macías  D. et al., “Retrieval of relevant parameters of natural multilayer systems by means of bio-inspired optimization strategies,” Appl. Opt.. 52, (11 ), 2511 –2520 (2013), CrossRef. 0003-6935 
Kennedy  J., Eberhart  R., “Particle swarm optimization,” in  IEEE Int. Conf. on Neural Networks , Vol. 4, pp. 1942 –1948,  IEEE ,  New York  (1995).
Kessentini  S., Barchiesi  D., “Quantitative comparison of optimized nanorods, nanoshells and hollow nanospheres for photothermal therapy,” Biomed. Opt. Express. 3, (3 ), 590 –604 (2012), CrossRef. 2156-7085 
Kessentini  S. et al., “Selective and collaborative optimization methods for plasmonics: a comparison,” PIERS Online. 7, (3 ), 291 –295 (2011), CrossRef. 1931-7360 
Barchiesi  D., “Numerical retrieval of thin aluminium layer properties from SPR experimental data,” Opt. Express. 20, (8 ), 9064 –9078 (2012), CrossRef. 1094-4087 
Shi  Y., Eberhart  R. C., “A modified particle swarm optimizer,” in  Proc. IEEE Congress on Evolutionary Computation , pp. 69 –73,  IEEE ,  New York  (1998).
Kessentini  S. et al., “Adaptive non-uniform particle swarm application to plasmonic design,” IJAMC. 2, (1 ), 18 –28 (2011), CrossRef.
Schwefel  H. P., Evolution and Optimum Seeking. ,  John Wiley & Sons Inc. ,  New York  (1995).
Macías  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, (8 ), 1465 –1471 (2004), CrossRef. 0740-3232 
Barchiesi  D., “Adaptive non-uniform, hyper-ellitist evolutionary method for the optimization of plasmonic biosensors,” in  Proc. Int. Conf. Computers & Industrial Engineering , pp. 542 –547,  IEEE ,  Troyes  (2009).
Vial  A. et al., “A new model of dispersion for metals leading to a more accurate modeling of plasmonic structures using the FDTD method,” Appl. Phys. A: Mater.. 103, (3 ), 849 –853 (2011), CrossRef. 0947-8396 
Krenn  J. R. et al., “Near-field optical investigations on nanometric silver particles,” in Photons and Local Probes. , Marti  O., Möller  R., Eds., NATO ASI Series E: Applied Sciences , Vol. 300, pp. 181 –187,  Kluwer Academic Publishers ,  Dordrecht  (1995).
Barchiesi  D., “Improved method based on S matrix for the optimization of SPR biosensors,” Opt. Commun.. 286, (1 ), 23 –29 (2013), CrossRef. 0030-4018 
Wautelet  M., Duvivier  D., “The characteristic dimensions of the nanoworld,” Eur. J. Phys.. 28, (5 ), 953 –959 (2007), CrossRef. 0143-0807 
Prokopidis  K. P., Zografopoulos  D. C., “A unified FDTD/PML scheme based on critical points for accurate studies of plasmonic structures,” J. Lightwave Technol.. 31, (15 ), 2467 –2476 (2013), CrossRef. 0733-8724 
Shultz  L. G., “The optical constants of silver, gold, copper and aluminum. I. The absorption coefficient k,” J. Opt. Soc. Am.. 44, (5 ), 357 –362 (1954), CrossRef. 0030-3941 
Shultz  L. G., Tangherlini  F. R., “The optical constants of silver, gold, copper and aluminum. II. The index of refraction n,” J. Opt. Soc. Am.. 44, (5 ), 362 –368 (1954), CrossRef. 0030-3941 
Chaari  A. et al., “Nanobubble evolution around nanowire in liquid,” Opt. Express. 21, (22 ), 26942 –26954 (2013), CrossRef. 1094-4087 
Fikri  R., Grosges  T., Barchiesi  D., “Apertureless scanning near-field optical microscopy: on the need of the tip vibration modelling,” Opt. Lett.. 28, (22 ), 2147 –2149 (2003), CrossRef. 0146-9592 
Fikri  R., Grosges  T., Barchiesi  D., “Apertureless scanning near-field optical microscopy: numerical modeling of the lock-in detection,” Opt. Commun.. 232, (1–6 ), 15 –23 (2004), CrossRef. 0030-4018 

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.