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

Dominique Barchiesi; Thomas Grosges

doi:10.1117/1.JNP.8.083097

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

Dominique Barchiesi, Thomas Grosges

Author Affiliations +

Journal of Nanophotonics, Vol. 8, Issue 1, 083097 (January 2014). https://doi.org/10.1117/1.JNP.8.083097

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.

1. Introduction

The fitting of optical constants of metals has given rise to an extensive literature. The optical properties of bulk materials¹^,² are commonly used for simulations, but recent experimental data for nanostructured materials will be used in the near future.³^,⁴ Indeed, the surface plasmon resonance (SPR) setups are able to determine thicknesses of multilayers⁴ and also optical properties.⁵ 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 used⁶^,⁷ for gold in a wide domain of wavelengths. Expansion of Lorentzian terms (L4) was proposed by Hao and Nordlander⁸ 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 model⁹^,¹⁰ was proved to be efficient for the modeling of spectroscopy with FDTD methods,¹¹ finite element method,¹² and discrete dipole approximation.¹³ Fitting of the optical constants is also useful if eigenvalues of complex structures are computed¹⁴^,¹⁵ or if resonances are searched in dispersion curves.¹⁶^,¹⁷ For the design and the optimization of nanostructures,¹⁸ the accuracy of numerical results depends on the quality of the fitting of the relative permittivities.¹⁹^,²⁰

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 data¹^,² 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.

2. Two Models for Fitting: the Combination of Drude and Lorentz Models, and the Critical Points Model

2.1.

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.²¹ 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$ .⁷ 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.⁷ The function of fit $ϵ_{DL}$ of the relative permittivity of metal is written as the sum of the Drude and the Lorentz models.⁶^,¹⁰^,²²

Eq. (1)

ϵ_{DL} (ω) = ϵ_{\infty} - \underset{Drude model}{\frac{ω_{D}^{2}}{\underset{︸}{ω (ω + i γ_{D})}}} - \underset{Lorentz model}{\underset{︸}{\frac{Δ ϵ ω_{L}^{2}}{ω^{2} - ω_{L}^{2} + i γ_{L} ω}}} .

The pure imaginary number is denoted as $i$ such as $i^{2} = - 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 $ϵ_{\infty}$ is the relative permittivity for high frequencies. For high frequencies, electrons cannot follow excitation, and the illumination does not see free electrons anymore; therefore $ϵ_{\infty} = 1$ .²³ 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 plasmon²¹), with $N_{e}$ the density of electrons that contribute to the optical properties [Eq. (3), it is not the valence band electron concentration²¹], $q$ the charge of electron, $m_{0}$ the electron mass, and $ϵ_{0}$ the relative permittivity of vacuum.

Eq. (2)

ω_{D} = \sqrt{\frac{N_{e} q^{2}}{ϵ_{0} m_{0}}} .

The volume density of the electrons that contribute to the relative permittivity of gold $N_{e}$ ( $m^{- 3}$ ) is introduced. It can be evaluated from an atomistic model and is actually a function of the density $ρ_{m}$ and atomic weight $M$ .

Eq. (3)

N_{e} = \frac{ρ_{m} N_{a}}{M} n_{e},

with

N_{a} = 6.022 \times 10^{23}

the Avogadro constant and

n_{e}

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

N_{e}

is calculated by using outer electron shell of atoms. Equation (1) could also be written as a product of

ω_{D}^{2}

by both Drude and Lorentz terms,²¹ but the formulation used in Ref. 7 is preferred to facilitate the comparison.

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

2.2.

Combination of Drude and Critical Points Models

The critical points model (CP) was introduced in 1998 by Leng et al.²⁴ for silicon and used for gold.⁹ 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.⁹^,²²^,²⁴ 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.⁹ This model helps to suppress the nonphysical absorption of the previous fitting in the region of transparency or near-transparency.²⁴ 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).

Eq. (4)

ϵ_{DCP} (ω) = ϵ_{\infty} - \frac{N_{e} q^{2}}{ϵ_{0} m_{0}} \frac{1}{ω (ω + i γ_{D})} + \sum_{i = 1}^{i = 2} Δ ϵ_{i} Ω_{i} [\frac{\exp (i ϕ_{i})}{Ω_{i} - ω - i Γ_{i}} + \frac{\exp (- i ϕ_{i})}{Ω_{i} + ω + i Γ_{i}}] .

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

2.3.

Criterion of Convergence for FDTD

According to Vial et al.,²²^,²⁵ a criterion must be verified to check the convergence of the FDTD algorithm.

Eq. (5)

C = \frac{ϵ_{\infty}}{ϵ_{\infty} + χ_{0}} < 1,

with

Eq. (6)

χ_{0} = - {(\frac{ω_{D}}{γ_{D}})}^{2} [1 - \exp (- γ_{D} Δ t)] + \frac{ω_{L}^{2}}{γ_{D}} Δ t + \sum R (- i \frac{η}{α - i β} {1 - \exp [(- α + i β) Δ t]}),

where

R (\cdot)

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$ , $β = \sqrt{ω_{L}^{2} - α^{2}}$ , and $η = Δ ϵ ω_{L}^{2} / β$ . For the CP, $χ_{0}$ is a function of $α = Γ_{i}$ , $β = Ω_{i}$ , and $η = 2 Δ ϵ_{i} Ω_{i} \exp (- i Φ_{i})$ .²⁵^,²⁶ In Eq. (6), $Δ t$ is evaluated from the size of the grid $Δ x$ used for FDTD: $Δ t = Δ x / (2 c)$ , with $c$ the speed of light in vacuum. In the following, we consider $Δ x = 1 nm$ .²⁵

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.

3. Combined Constrained PSO with Nelder-Mead Algorithm

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.⁷ The NM method²⁷ 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)$ ].²⁸ The simulated annealing subroutine belongs to the same class as the Monte Carlo method (heuristic methods) and is fully described by Kirkpatrick et al.²⁹ and Tarantola.³⁰ 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.³¹^–³³ The PSO method was proposed in 1995 by Kennedy and Eberhart.³⁴ It was used for the optimization of gold nanoshells using discrete dipole approximation (DDA) and Mie models for photothermal therapy,³⁵ and more general problems of plasmonics including the optimization of nanostructures and the resolution of the inverse problem.⁵^,³⁶^,³⁷

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 $N_{ref}$ values of the relative permittivity of reference $ϵ_{r}$ .

Eq. (7)

F (x) = \sqrt{\frac{1}{N_{ref}} \sum_{1}^{N_{ref}} {| ϵ_{M} (x) - ϵ_{r} |}^{2}} .

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} (ω)$ .

Eq. (8)

σ_{R} [ϵ_{M} (ω)] = \sqrt{\frac{1}{N_{ref}} R {[ϵ_{M} (ω) - ϵ_{r} (ω)]}^{2}},

Eq. (9)

σ_{I} [ϵ_{M} (ω)] = \sqrt{\frac{1}{N_{ref}} I {[ϵ_{M} (ω) - ϵ_{r} (ω)]}^{2}} .

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).

Eq. (10)

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

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)$ .³⁸

Eq. (11)

V (t + 1) = ω_{PSO} V (t) + U_{1} c_{1} [p (t) - x (t)] + U_{2} c_{2} [g (t) - x (t)] .

The inertia weight $ω_{PSO}$ linearly decreases from 0.9 to 0.4.³⁹ The random gene $U_{i} (i = 1, 2)$ ] always checking the condition given in Eq. (5). The acceleration coefficients are $c_{1} = c_{2} = 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}$ , $c_{1}$ , $c_{2}$ , and $N$ ). In Ref. 37, this property and the efficiency of this method was verified by using evolutionary method based on evolution strategies.⁴⁰^–⁴²

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.

4. Results and Discussion

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,⁷^,⁴³ 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 [10^{13}, 10^{18}] (rad . s^{- 1})$ , $ϵ_{\infty} \in [- 1000,1000]$ , $Δ ϵ \in [0,1000]$ , $N_{e} b [10^{27}, 10^{30}] (m^{- 3})$ [Eq. (3)]. The fitting of bulk data for gold, silver, chromium, aluminum, and titanium are, respectively, investigated.

4.1.

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.² Figure 1 shows the results of fitting and the absolute errors $σ_{R}$ and $σ_{I}$ [Eqs. (8) and (9)].

Table 1

Fitting 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.

	${Au}^{2}$
	DL	DCP
$F$ [Eq. (7)]	0.55	0.15992
$C$ [Eq. (5)]	0.99995	0.92761
$σ_{R}$	1.3292	0.57097
$σ_{I}$	1.2208	1.0747
$ϵ_{\infty}$	6.15991	$- 9.06407$
$γ_{D}$ ( $rad / s$ )	$1.66938 \times 10^{15}$	$5.86665 \times 10^{13}$
$ω_{D}$ ( $rad / s$ )	$1.34759 \times 10^{16}$	$1.30423 \times 10^{16}$
$Δ ϵ$	2.07122	—
$ω_{L}$ ( $rad / s$ )	$4.66171 \times 10^{15}$	—
$γ_{L}$ ( $rad / s$ )	$7.20958 \times 10^{13}$	—
$Ω_{1}$ ( $rad / s$ )	—	$3.2539 \times 10^{16}$
$Γ_{1}$ ( $rad / s$ )	—	$5.5350 \times 10^{15}$
$Δ ϵ_{1}$	—	$- 10.8876$
$Φ_{1}$ (rad)	—	$- 2.46009$
$Ω_{2}$ ( $rad / s$ )	—	$3.91172 \times 10^{15}$
$Γ_{2}$ ( $rad / s$ )	—	$6.95449 \times 10^{14}$
$Δ ϵ_{2}$	—	0.718455
$ϕ_{2}$ (rad)	—	$- 1.13717$

Fig. 1

(a) Real and imaginary part of the experimental data for gold² 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 Ref. 22 ( $ϵ_{V}$ ) and Ref. 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$ .⁷ The proposed method gives $σ_{R} = 1.3292$ and $σ_{I} = 1.2208$ (Table 1).
• With DCP, $σ_{R} (ϵ_{V}) = 0.9417$ , $σ_{I} (ϵ_{V}) = 1.0606$ ,²² $σ_{R} (ϵ_{V}^{'}) = 0.7882$ , and $σ_{I} (ϵ_{V}^{'}) = 1.2268$ .⁴³ 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,⁷^,²²^,⁴³ 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 2

Fitting 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.

	${Au}^{1}$
	DL	DCP
$F$ [Eq. (7)]	1.1641	0.69043
$C$ [Eq. (5)]	0.99999	0.98792
$σ_{R}$	7.7965	0.57097
$σ_{I}$	2.7878	1.0747
$ϵ_{\infty}$	6.90939	$- 6.23210$
$γ_{D}$ ( $rad / s$ )	$1.75628 \times 10^{15}$	$3.17891 \times 10^{17}$
$ω_{D}$ ( $rad / s$ )	$1.38147 \times 10^{16}$	$1.90105 \times 10^{17}$
$Δ ϵ$	2.31858	—
$ω_{L}$ ( $rad / s$ )	$4.68266 \times 10^{15}$	—
$γ_{L}$ ( $rad / s$ )	$3.60439 \times 10^{14}$	—
$Ω_{1}$ ( $rad / s$ )	—	$3.37278 \times 10^{15}$
$Γ_{1}$ ( $rad / s$ )	—	$3.69375 \times 10^{14}$
$Δ ϵ_{1}$	—	0.419451
$Φ_{1}$ (rad)	—	$- 4.44248$
$Ω_{2}$ ( $rad / s$ )	—	$2.04724 \times 10^{15}$
$Γ_{2}$ ( $rad / s$ )	—	$1.23992 \times 10^{15}$
$Δ ϵ_{2}$	—	28.2908
$ϕ_{2}$ (rad)	—	$- 3.52016$

Fig. 2

(a) Real and imaginary part of the experimental data for gold¹ 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.

4.2.

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.⁴⁴ Figure 3 shows both fittings and the error of fitting $σ (ω)$ . The real and imaginary parts of $ϵ_{V}$ ²² 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 ²² and the PSO results.

Fig. 3

(a) Real and imaginary part of the experimental data for silver¹ 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 Ref. 22 ( $ϵ_{V}$ ) and Ref. 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$ ,²² $σ_{R} (ϵ_{V}^{'}) = 0.81118$ , and $σ_{I} (ϵ_{V}^{'}) = 0.70911$ .⁴³ 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 Ref. 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 3

Fitting 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.

	${Ag}^{1}$
	DL	DCP
$F$ [Eq. (7)]	0.071605	0.068131
$C$ [Eq. (5)]	0.98913	0.28263
$σ_{R}$	0.37765	0.36891
$σ_{I}$	0.33703	0.34168
$ϵ_{\infty}$	0.114773	$- 19.4464$
$γ_{D}$ ( $rad / s$ )	$7.05499 \times 10^{15}$	$9.52533 \times 10^{13}$
$ω_{D}$ ( $rad / s$ )	$1.32589 \times 10^{16}$	$1.32472 \times 10^{16}$
$Δ ϵ$	3.62762	—
$ω_{L}$ ( $rad / s$ )	$1.58116 \times 10^{16}$	—
$γ_{L}$ ( $rad / s$ )	$1.04632 \times 10^{14}$	—
$Ω_{1}$ ( $rad / s$ )	—	$3.47550 \times 10^{16}$
$Γ_{1}$ ( $rad / s$ )	—	$5.18324 \times 10^{16}$
$Δ ϵ_{1}$	—	1727.14
$Φ_{1}$ (rad)	—	1.04286
$Ω_{2}$ ( $rad / s$ )	—	$4.76055 \times 10^{16}$
$Γ_{2}$ ( $rad / s$ )	—	$3.23542 \times 10^{16}$
$Δ ϵ_{2}$	—	$- 742.061$
$ϕ_{2}$ (rad)	—	1.75438

4.3.

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.³⁷ Table 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 ²²^,⁴³ and the proposed fitting are about the same order of magnitude along $ω$ except for wavelengths $> 630 nm$ . 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 4

Fitting 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.

	${Al}^{1}$
	DL	DCP
$F$ [Eq. (7)]	2.8752	0.29101
$C$ [Eq. (5)]	0.050662	0.27149
$σ_{R}$	8.4152	1.5629
$σ_{I}$	7.9598	1.1941
$ϵ_{\infty}$	0.0864080	$5.79805 \times 10^{- 3}$
$γ_{D}$ ( $rad / s$ )	$4.75280 \times 10^{15}$	$6.78974 \times 10^{15}$
$ω_{D}$ ( $rad / s$ )	$3.91780 \times 10^{16}$	$2.22564 \times 10^{17}$
$Δ ϵ$	$- 3814.53$	—
$ω_{L}$ ( $rad / s$ )	$5.67222 \times 10^{14}$	—
$γ_{L}$ ( $rad / s$ )	$2.76612 \times 10^{15}$	—
$Ω_{1}$ ( $rad / s$ )	—	$8.10015 \times 10^{13}$
$Γ_{1}$ ( $rad / s$ )	—	$1.92552 \times 10^{16}$
$Δ ϵ_{1}$	—	$- 1.35397 \times 10^{3}$
$Φ_{1}$ (rad)	—	$- 0.683343$
$Ω_{2}$ ( $rad / s$ )	—	$2.25905 \times 10^{16}$
$Γ_{2}$ ( $rad / s$ )	—	$2.27300 \times 10^{15}$
$Δ ϵ_{2}$	—	3.59426
$ϕ_{2}$ (rad)	—	$- 1.09763$

Fig. 4

(a) Real and imaginary part of the experimental data for aluminum¹ 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.

4.4.

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 ²² and the PSO results.

Table 5

Fitting 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.

	${Cr}^{1}$
	DL	DCP
$F$ [Eq. (7)]	0.94723	0.077607
$C$ [Eq. (5)]	0.99993	0.88915
$σ_{R}$	4.9121	0.29314
$σ_{I}$	5.625	0.48649
$ϵ_{\infty}$	2.77667	0.386846
$γ_{D}$ ( $rad / s$ )	$2.57644 \times 10^{15}$	$1.80690 \times 10^{15}$
$ω_{D}$ ( $rad / s$ )	$1.59078 \times 10^{16}$	$1.39350 \times 10^{16}$
$Δ ϵ$	13.2908	—
$ω_{L}$ ( $rad / s$ )	$3.32972 \times 10^{15}$	—
$γ_{L}$ ( $rad / s$ )	$2.99660 \times 10^{15}$	—
$Ω_{1}$ ( $rad / s$ )	—	$3.79297 \times 10^{15}$
$Γ_{1}$ ( $rad / s$ )	—	$8.01774 \times 10^{14}$
$Δ ϵ_{1}$	—	2.12352
$Φ_{1}$ (rad)	—	0.883949
$Ω_{2}$ ( $rad / s$ )	—	$1.75789 \times 10^{15}$
$Γ_{2}$ ( $rad / s$ )	—	$7.80710 \times 10^{14}$
$Δ ϵ_{2}$	—	11.8586
$ϕ_{2}$ (rad)	—	$- 1.69593$

Fig. 5

(a) Real and imaginary part of the experimental data for chromium¹ 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.⁴⁵ 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.

4.5.

Titanium

Titanium is used as an alternative medium for adhesion layer of gold on dielectric samples.⁴⁵ Table 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 6

Fitting 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.

	${Ti}^{1}$
	DL	DCP
$F$ [Eq. (7)]	0.62496	0.26951
$C$ [Eq. (5)]	0.99959	0.95922
$σ_{R}$	4.418	1.2955
$σ_{I}$	1.4696	1.483
$ϵ_{\infty}$	2.17069	1.35312
$γ_{D}$ ( $rad / s$ )	$1.08330 \times 10^{16}$	$3.02914 \times 10^{15}$
$ω_{D}$ ( $rad / s$ )	$8.50868 \times 10^{15}$	$2.00924 \times 10^{16}$
$Δ ϵ$	74.4496	—
$ω_{L}$ ( $rad / s$ )	$2.93140 \times 10^{15}$	—
$γ_{L}$ ( $rad / s$ )	$1.46986 \times 10^{14}$	—
$Ω_{1}$ ( $rad / s$ )	—	$2.52097 \times 10^{15}$
$Γ_{1}$ ( $rad / s$ )	—	$2.09204 \times 10^{15}$
$Δ ϵ_{1}$	—	9.01823
$Φ_{1}$ (rad)	—	$- 2.06436$
$Ω_{2}$ ( $rad / s$ )	—	$1.88081 \times 10^{15}$
$Γ_{2}$ ( $rad / s$ )	—	$9.05493 \times 10^{13}$
$Δ ϵ_{2}$	—	3.90173
$ϕ_{2}$ (rad)	—	2.76388

Fig. 6

(a) Real and imaginary part of the experimental data for titanium² 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).

4.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 7

Performance 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

	Au (JC)	Au (P)	Ag (P)	Al (P)	Cr (P)	Ti (JC)
DCP	+++	+	+++	++	+++	++
DL	++	−	+++	−	+	+

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 = 1 nm$ 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.

5. Conclusion

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.⁴⁴^,⁴⁶

The results of fitting can be used directly for any spectroscopic simulation and especially in FDTD codes.⁸^,⁴⁷ 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.⁴⁸^–⁵⁰ The advantage of such a heuristic method lies in its applicability to various problems of fitting, the optimization of complex systems in engineering,¹³^,³⁵ the propagation of uncertainties, or the tolerance study of models.¹⁸^,²⁸ The use of more complex models ¹⁵^,¹⁶^,⁵¹^,⁵² can extend the domain of its application.

Acknowledgments

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

References

1.

E. D. Palik, Handbook of Optical Constants, Academic Press Inc., San Diego (1985). Google Scholar

2.

P. JohnsonR. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B, 6 (12), 4370 –4379 (1972). http://dx.doi.org/10.1103/PhysRevB.6.4370 PRBMDO 0163-1829 Google Scholar

3.

F. FloryL. EscoubasG. Berginc, “Optical properties of nanostructured materials: a review,” J. Nanophotonics, 5 (1), 052502 (2011). http://dx.doi.org/10.1117/1.3609266 JNOACQ 1934-2608 Google Scholar

4.

S. H. El-GoharyN.-H. KimK. M. Byun, “Optical determination of thick graphene layer number based on surface plasmon resonance,” J. Nanophotonics, 7 (1), 073799 (2013). http://dx.doi.org/10.1117/1.JNP.7.073799 JNOACQ 1934-2608 Google Scholar

5.

J. SalviD. Barchiesi, “Measurement of thicknesses and optical properties of thin films from surface plasmon resonance (SPR),” Appl. Phys. B—Lasers Opt., 114 (2014). http://dx.doi.org/10.1007/s00339-013-8038-z APBOEM 0946-2171 Google Scholar

6.

A. D. Rakićet al., “Optical properties of metallic films for vertical-cavity optoelectronic devices,” Appl. Opt., 37 (22), 5271 –5283 (1998). http://dx.doi.org/10.1364/AO.37.005271 APOPAI 0003-6935 Google Scholar

7.

A. Vialet 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). http://dx.doi.org/10.1103/PhysRevB.71.085416 PRBMDO 0163-1829 Google Scholar

8.

F. HaoP. Nordlander, “Efficient dielectric function for FDTD simulation of the optical properties of silver and gold nanoparticles,” Chem. Phys. Lett., 446 (1–3), 115 –118 (2007). http://dx.doi.org/10.1016/j.cplett.2007.08.027 CHPLBC 0009-2614 Google Scholar

9.

P. G. EtchegoinE. C. L. RuM. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys., 125 164705 (2006). http://dx.doi.org/10.1063/1.2360270 JCPSA6 0021-9606 Google Scholar

10.

S. G. RodrigoF. J. García-VidalL. Martín-Moreno, “Influence of material properties on extraordinary optical transmission through hole arrays,” Phys. Rev. B, 77 (7), 075401 (2008). http://dx.doi.org/10.1103/PhysRevB.77.075401 PRBMDO 0163-1829 Google Scholar

11.

A. VialT. Laroche, “Comparison of gold and silver dispersion laws suitable for FDTD simulations,” Appl. Phys. B—Lasers Opt., 93 (1), 139 –143 (2008). http://dx.doi.org/10.1007/s00340-008-3202-4 APBOEM 0946-2171 Google Scholar

12.

T. GrosgesA. VialD. Barchiesi, “Models of near-field spectroscopic studies: comparison between finite-element and finite-difference methods,” Opt. Express, 13 (21), 8483 –8497 (2005). http://dx.doi.org/10.1364/OPEX.13.008483 OPEXFF 1094-4087 Google Scholar

13.

D. Barchiesiet al., “Localized surface plasmon resonance in arrays of nano-gold cylinders: inverse problem and propagation of uncertainties,” Opt. Express, 21 (2), 2245 –2262 (2013). http://dx.doi.org/10.1364/OE.21.002245 OPEXFF 1094-4087 Google Scholar

14.

B. GuizalH. YalaD. Felbacq, “Reformulation of the eigenvalue problem in the Fourier modal method with spatial adaptive resolution,” Opt. Lett., 34 (18), 2790 –2792 (2009). http://dx.doi.org/10.1364/OL.34.002790 OPLEDP 0146-9592 Google Scholar

15.

B. GuizalD. BarchiesiD. Felbacq, “Electromagnetic beam diffraction by a finite lamellar structure,” J. Opt. Soc. Am. A, 20 (12), 2274 –2280 (2003). http://dx.doi.org/10.1364/JOSAA.20.002274 JOAOD6 0740-3232 Google Scholar

16.

D. Barchiesiet al., “A Poincaré’s approach for plasmonics: the plasmon localization,” J. Microsc., 229 (3), 525 –532 (2008). http://dx.doi.org/10.1111/j.1365-2818.2008.01938.x JMICAR 0022-2720 Google Scholar

17.

A. CastaniéD. Felbacq, “Confined plasmonic modes in a nanocavity,” Opt. Commun., 285 (16), 3353 –3357 (2012). http://dx.doi.org/10.1016/j.optcom.2012.01.039 OPCOB8 0030-4018 Google Scholar

18.

T. Grosgeset al., “Design of nanostructures for imaging and biomedical applications by plasmonic optimization,” Opt. Lett., 33 (23), 2812 –2814 (2008). http://dx.doi.org/10.1364/OL.33.002812 OPLEDP 0146-9592 Google Scholar

19.

D. BarchiesiB. GuizalT. Grosges, “Accuracy of local field enhancement models: toward predictive models?,” Appl. Phys. B—Lasers Opt., 84 (1–2), 55 –60 (2006). http://dx.doi.org/10.1007/s00340-006-2217-y APBOEM 0946-2171 Google Scholar

20.

D. Barchiesiet al., “Plasmonics: influence of the intermediate (or stick) layer on the efficiency of sensors,” Appl. Phys. B—Lasers Opt., 93 (1), 177 –181 (2008). http://dx.doi.org/10.1007/s00340-008-3173-5 APBOEM 0946-2171 Google Scholar

21.

A. D. Rakić, “Algorithm for the determination of intrinsic optical constants of metal films: application to aluminum,” Appl. Opt., 34 (22), 4755 –4767 (1995). http://dx.doi.org/10.1364/AO.34.004755 APOPAI 0003-6935 Google Scholar

22.

A. VialT. Laroche, “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). http://dx.doi.org/10.1088/0022-3727/40/22/043 JPAPBE 0022-3727 Google Scholar

23.

M. N. Polyanskiy, “Refractive index database,” (2013) http://refractiveindex.info://refractiveindex.info October ). 2013). Google Scholar

24.

J. Lenget al., “Analytic representations of the dielectric functions of materials for device and structural modeling,” Thin Solid Films, 313–314 132 –136 (1998). http://dx.doi.org/10.1016/S0040-6090(97)00799-2 THSFAP 0040-6090 Google Scholar

25.

A. Vial, “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). http://dx.doi.org/10.1088/1464-4258/9/7/029 JOAOF8 1464-4258 Google Scholar

26.

I. B. UdagedaraI. D. RukhlenkoM. Premaratne, “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). http://dx.doi.org/10.1364/OE.19.019973 OPEXFF 1094-4087 Google Scholar

27.

J. NelderR. Mead, “A simplex method for function minimization,” Comput. J., 7 (4), 308 –313 (1965). http://dx.doi.org/10.1093/comjnl/7.4.308 CMPJA6 0010-4620 Google Scholar

28.

T. Grosgeset al., “Nanoshells for photothermal therapy: a Monte-Carlo based numerical study of their design tolerance,” Biomed. Opt. Express, 2 (6), 1584 –1596 (2011). http://dx.doi.org/10.1364/BOE.2.001584 BOEICL 2156-7085 Google Scholar

29.

S. KirkpatrickC. D. GelattM. P. Vecchi, “Optimization by simulated annealing,” Science, 220 (4598), 671 –680 (1983). http://dx.doi.org/10.1126/science.220.4598.671 SCIEAS 0036-8075 Google Scholar

30.

A. Tarantola, Inverse Problem Theory and Methods for Model Parameter Estimation, Society for Industrial and Applied Mathematics (SIAM), Philadelphia (2005). Google Scholar

31.

L. Rodriguez-CoboA. CoboJ. M. Lopez-Higuera, “Recovering a fiber Bragg grating axial strain distribution from its reflection spectrum,” Opt. Lett., 38 (13), 2327 –2329 (2013). http://dx.doi.org/10.1364/OL.38.002327 OPLEDP 0146-9592 Google Scholar

32.

C. Yanget al., “Design of reflective color filters with high angular tolerance by particle swarm optimization method,” Opt. Express, 21 (8), 9315 –9323 (2013). http://dx.doi.org/10.1364/OE.21.009315 OPEXFF 1094-4087 Google Scholar

33.

D. Macíaset al., “Retrieval of relevant parameters of natural multilayer systems by means of bio-inspired optimization strategies,” Appl. Opt., 52 (11), 2511 –2520 (2013). http://dx.doi.org/10.1364/AO.52.002511 APOPAI 0003-6935 Google Scholar

34.

J. KennedyR. Eberhart, “Particle swarm optimization,” in IEEE Int. Conf. on Neural Networks, 1942 –1948 (1995). Google Scholar

35.

S. KessentiniD. Barchiesi, “Quantitative comparison of optimized nanorods, nanoshells and hollow nanospheres for photothermal therapy,” Biomed. Opt. Express, 3 (3), 590 –604 (2012). http://dx.doi.org/10.1364/BOE.3.000590 BOEICL 2156-7085 Google Scholar

36.

S. Kessentiniet al., “Selective and collaborative optimization methods for plasmonics: a comparison,” PIERS Online, 7 (3), 291 –295 (2011). http://dx.doi.org/10.2529/PIERS100902052449 1931-7360 Google Scholar

37.

D. Barchiesi, “Numerical retrieval of thin aluminium layer properties from SPR experimental data,” Opt. Express, 20 (8), 9064 –9078 (2012). http://dx.doi.org/10.1364/OE.20.009064 OPEXFF 1094-4087 Google Scholar

38.

Y. ShiR. C. Eberhart, “A modified particle swarm optimizer,” in Proc. IEEE Congress on Evolutionary Computation, 69 –73 (1998). Google Scholar

39.

S. Kessentiniet al., “Adaptive non-uniform particle swarm application to plasmonic design,” IJAMC, 2 (1), 18 –28 (2011). http://dx.doi.org/10.4018/IJAMC Google Scholar

40.

H. P. Schwefel, Evolution and Optimum Seeking, John Wiley & Sons Inc., New York (1995). Google Scholar

41.

D. MacíasA. VialD. Barchiesi, “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). http://dx.doi.org/10.1364/JOSAA.21.001465 JOAOD6 0740-3232 Google Scholar

42.

D. Barchiesi, “Adaptive non-uniform, hyper-ellitist evolutionary method for the optimization of plasmonic biosensors,” in Proc. Int. Conf. Computers & Industrial Engineering, 542 –547 (2009). Google Scholar

43.

A. Vialet 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). http://dx.doi.org/10.1007/s00339-010-6224-9 APAMFC 0947-8396 Google Scholar

44.

J. R. Krennet al., “Near-field optical investigations on nanometric silver particles,” Photons and Local Probes, 181 –187 Kluwer Academic Publishers, Dordrecht (1995). Google Scholar

45.

D. Barchiesi, “Improved method based on S matrix for the optimization of SPR biosensors,” Opt. Commun., 286 (1), 23 –29 (2013). http://dx.doi.org/10.1016/j.optcom.2012.08.060 OPCOB8 0030-4018 Google Scholar

46.

M. WauteletD. Duvivier, “The characteristic dimensions of the nanoworld,” Eur. J. Phys., 28 (5), 953 –959 (2007). http://dx.doi.org/10.1088/0143-0807/28/5/018 EJPHD4 0143-0807 Google Scholar

47.

K. P. ProkopidisD. C. Zografopoulos, “A unified FDTD/PML scheme based on critical points for accurate studies of plasmonic structures,” J. Lightwave Technol., 31 (15), 2467 –2476 (2013). http://dx.doi.org/10.1109/JLT.2013.2265166 JLTEDG 0733-8724 Google Scholar

48.

L. G. Shultz, “The optical constants of silver, gold, copper and aluminum. I. The absorption coefficient k,” J. Opt. Soc. Am., 44 (5), 357 –362 (1954). http://dx.doi.org/10.1364/JOSA.44.000357 JOSAAH 0030-3941 Google Scholar

49.

L. G. ShultzF. R. Tangherlini, “The optical constants of silver, gold, copper and aluminum. II. The index of refraction n,” J. Opt. Soc. Am., 44 (5), 362 –368 (1954). http://dx.doi.org/10.1364/JOSA.44.000362 JOSAAH 0030-3941 Google Scholar

50.

A. Chaariet al., “Nanobubble evolution around nanowire in liquid,” Opt. Express, 21 (22), 26942 –26954 (2013). http://dx.doi.org/10.1364/OE.21.026942 OPEXFF 1094-4087 Google Scholar

51.

R. FikriT. GrosgesD. Barchiesi, “Apertureless scanning near-field optical microscopy: on the need of the tip vibration modelling,” Opt. Lett., 28 (22), 2147 –2149 (2003). http://dx.doi.org/10.1364/OL.28.002147 OPLEDP 0146-9592 Google Scholar

52.

R. FikriT. GrosgesD. Barchiesi, “Apertureless scanning near-field optical microscopy: numerical modeling of the lock-in detection,” Opt. Commun., 232 (1–6), 15 –23 (2004). http://dx.doi.org/10.1016/j.optcom.2003.12.027 OPCOB8 0030-4018 Google Scholar

Biography

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.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.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 Download Citation

Dominique Barchiesi and Thomas Grosges "Fitting the optical constants of gold, silver, chromium, titanium, and aluminum in the visible bandwidth," Journal of Nanophotonics 8(1), 083097 (6 January 2014). https://doi.org/10.1117/1.JNP.8.083097

Published: 6 January 2014

Access the abstract

JOURNAL ARTICLE
17 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 94 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Gold

Aluminum

Chromium

Silver

Titanium

Finite-difference time-domain method

Particle swarm optimization

1.

Introduction

2.

Two Models for Fitting: the Combination of Drude and Lorentz Models, and the Critical Points Model

2.1.

Combination of Drude and Lorentz Models

Eq. (1)

Eq. (2)

Eq. (3)

2.2.

Combination of Drude and Critical Points Models

Eq. (4)

2.3.

Criterion of Convergence for FDTD

Eq. (5)

Eq. (6)

3.

Combined Constrained PSO with Nelder-Mead Algorithm

Eq. (7)

Eq. (8)

Eq. (9)

Eq. (10)

Eq. (11)

4.

Results and Discussion

4.1.

Gold

Table 1

Fig. 1

Table 2

Fig. 2

4.2.

Silver

Fig. 3

Table 3

4.3.

Aluminum

Table 4

Fig. 4

4.4.

Chromium

Table 5

Fig. 5

4.5.

Titanium

Table 6

Fig. 6

4.6.

Concluding Remarks

Table 7

5.

Conclusion

Acknowledgments

References

Biography

Show All Keywords

Keywords/Phrases

Search In:

Publication Years