Terahertz plasmon and surface-plasmon modes in hollow nanospheres

1556-276X-7-578 1556-276X Nano Express Terahertz plasmon and surface-plasmon modes in hollow nanospheres XiaoYimingymxiao_ynu@163.com XuWenwenxu_issp@yahoo.cn ZhangYayayyzhang_ynu@163.com HuJiaguangynshjg@126.com

Department of Physics, Yunnan University, Kunming, 650091, China

Key Laboratory of Materials Physics, Institute of Solid State Physics, Chinese Academy of Sciences, Hefei, 230031, China

Department of Math and Physics, Wenshan University, Wenshan, 663000, China

Nanoscale Research Letters

SI: International Conference on Superlattices, Nanostructures, and Nanodevices (ICSNN 2012)

1556-276X 2012 7 1 578 http://www.nanoscalereslett.com/content/7/1/578 10.1186/1556-276X-7-57823092121 1972012810201223102012 2012Xiao et al.; licensee Springer.This is an Open Access article distributed under the terms of the Creative Commons Attribution License( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Hollow nanosphere Electronic subband structure Collective excitation modes Terahertz radiation

Abstract

We present a theoretical study of the electronic subband structure and collective electronic excitation associated with plasmon and surface plasmon modes in metal-based hollow nanosphere. The dependence of the electronic subband energy on the sample parameters of the hollow nanosphere is examined. We find that the subband states with different quantum numbers l degenerate roughly when the outer radius of the sphere is r ₂ ≥ 100 nm. In this case, the energy spectrum of a sphere is mainly determined by quantum number n. Moreover, the plasmon and surface plasmon excitations can be achieved mainly via inter-subband transitions from occupied subbands to unoccupied subbands. We examine the dependence of the plasmon and surface-plasmon frequencies on the shell thickness d and the outer radius r ₂ of the sphere using the standard random-phase approximation. We find that when a four-state model is employed for calculations, four branches of the plasmon and surface plasmon oscillations with terahertz frequencies can be observed, respectively.

International Conference on Superlattices, Nanostructures, and Nanodevices (ICSNN 2012)Rinaldo Trotta, Manfred Helm and Oliver G. Schmidt

Background

In recent years, there has been a great interest in the investigation of metal-based hollow nanostructures because of their unique characteristics such as low density, large specific area, mechanical and thermal stability, and surface permeability. These advanced materials have been widely applied in catalysis 1 , drug delivery 2 3 , food and cosmetic industries 4 , fuel cell 5 6 , biotechnology 7 , lubricant 8 , sensing 9 , photonic devices 10 , micro/nanoreactors 11 , etc. In particular, metal-based hollow nanospheres 12 can be realized via using polystyrene (PS) latex particles as templates 13 . Such structures have intriguing features of surface plasmon resonance 14 . The collective oscillations of the conducting electrons in response to optical excitation, such as plasmon and surface plasmon excitations, affect strongly the optical properties of metal hollow nanospheres. At present, it has become possible to fabricate metal hollow nanosphere structures in which the radius and shell thickness of the sphere can be controlled artificially. Such structures have been widely applied to realize terahertz (10¹² Hz or THz) plasmonic devices 15 . Hence, it is of great importance and significance to study the electronic subband structure and corresponding collective electronic excitations from these advanced nanomaterial systems. In conjunction with recent experimental achievement in the field, in this article, we would like to develop a simple theoretical approach to study the electronic subband structure and plasmon and surface-plasmon modes in a hollow nanosphere. The aim of this study is to examine how sample parameters affect the electronic subband energy and the plasmon and the surface-plasmon modes in the device systems.

Methods

Theoretical approach

Electronic subband structure

In this study, we consider an air/metal-shell/air-based hollow nanosphere structure. The inner radius of the structure is r ₁, the outer radius or the diameter of the sphere is r ₂, and the metal shell thickness is d = r ₂ − r ₁. Such structure can be realized experimentally by selectively removing the hard spherical core (e.g., PS latex particles) in the fabrication process 16 . For a case where the electrons in the metal shell are not tunneling or hopping into the core and outside air, the confining potential for electrons in the structure can be modeled simply as

V ( r ) = 0 , r 1 ≤ r ≤ r 2 ; ∞ , otherwise .

Thus, the corresponding Schrödinger equation takes a form

[ P 2 / 2 μ + V ( r ) ] ψ N ( R ) = E N ψ N ( R ) .

Here, P = (p _x,p _y,p _z) is the momentum operator, μ is the effective mass for an electron in the structure, R = (x,y,z) = (r,θ,ϕ), and N stands for all quantum numbers. The solution of Equation (2) is ψ N ( r , θ , ϕ ) = R nl ( r ) Y l m ( θ ) e im ϕ , where N = (nlm),

Y l m ( θ ) = ( − 1 ) ( m + | m | ) / 2 ( 2 l + 1 ) ( l − | m | ) ! 4 Π ( l + | m | ) ! 1 / 2 × P l | m | ( cos θ ) ,

and R _nl(r) is determined by

1 r 2 ∂ ∂r r 2 ∂ R nl ( r ) ∂ r − l ( l + 1 ) r 2 R nl ( r ) = − 2 μ ℏ 2 [ E N − V ( r ) ] R nl ( r ) .

Here, l = 0,1,2,⋯ is the angular momentum quantum number, m = l,l−1,⋯,−l is the magnetic quantum number, P l m ( x ) is the associated Legendre function, and l must be a positive integer in the range l ≥ |m|. Letting E _N= ℏ ² k ²/2μ, x = kr, and R nl ( r ) = Π / 2 x y ( x ) , the radial equation, Equation (4), for r ₁ ≤ r ≤ r ₂ becomes a Bessel equation with a general solution: y ( x ) = N nl [ C 1 J l + 1 / 2 ( x ) + C 2 J − l − 1 / 2 ( x ) ] , where J _α(x) is a Bessel function and N nl is a normalization factor. Considering the boundary conditions: R(r ₁) = 0 and R(r ₂) = 0, we have

J l + 1 / 2 ( k r 2 ) J − l − 1 / 2 ( k r 1 ) − J l + 1 / 2 ( k r 1 ) J − l − 1 / 2 ( k r 2 ) = 0 ,

which is applied to determine the energy spectrum of the sample structure. Thus, the electron wave function becomes

ψ N ( r , θ , ϕ ) = ( − 1 ) ( m + | m | ) / 2 ( 2 l + 1 ) ( l − | m | ) ! 4 Π ( l + | m | ) ! 1 / 2 e im ϕ

× P l | m | ( cos θ ) N nl Π / 2 kr [ C 1 J l + 1 / 2 ( kr ) + C 2 J − l − 1 / 2 ( kr ) ] .

For case of l = 0, we obtain E _n0 = ℏ ² Π ² n ²/2μ d ² with n = 1,2,⋯. The radial eigenfunction is

R n 0 ( r ) = 2 cos 2 ( k r 1 ) d sin ( kr ) r − tan ( k r 1 ) cos ( kr ) r ,

where k = nΠ/d.

For l = 1, we have

tan ( kd ) = kd 1 + k 2 r 1 r 2 ,

and the radial eigenfunction becomes

R n 1 ( r ) = N n 1 r sin ( kr ) kr − cos ( kr ) − tan ( k r 1 ) − k r 1 k r 1 tan ( k r 1 ) + 1 × sin ( kr ) + cos ( kr ) kr .

For l = 2, we get

tan ( kd ) = ( 9 + 3 k 2 r 1 r 2 ) kd k 2 ( k 2 r 1 2 r 2 2 − 3 r 1 2 − 3 r 2 2 + 9 r 1 r 2 ) + 9 .

E_n1 and E_n2 are determined numerically via solving respectively Equations (8) and (10).

Electron-electron interaction

The matrix element for the bare electron-electron (e-e) interaction can be obtained by applying the electron wave function to the interaction Hamiltonian induced by the Coulomb potential 17 , which reads

V N 1 ′ N 1 N ′ N = e 2 κ ∫ ∫ ψ N 1 ′ ∗ ( R 1 , θ 1 , ϕ 1 ) ψ N 1 ( R 1 , θ 1 , ϕ 1 ) × 1 | R 1 − R 2 | ψ N ′ ∗ ( R 2 , θ 2 , ϕ 2 ) ψ N ( R 2 , θ 2 , ϕ 2 ) ,

with κ being the high-frequency dielectric constant of the shell material. It can be simplified as

V N 1 ′ N 1 N ′ N = ∑ k c k ( N , N ′ ) c k ( N 1 ′ , N 1 ) R k ( N N 1 , N ′ N 1 ′ ) × δ m N + m N 1 , m N ′ + m N 1 ′ .

Here, c k and R k are, respectively,

c k ( N , N ′ ) = ( − 1 ) m N − m N ′ c k ( N ′ , N ) = 2 2 k + 1 × ∫ 0 Π T ( k , m N − m N ′ ) T ( l N , m N ) T ( l N ′ , m N ′ ) sin θ d θ ,

where T ( l , m ) = 2 Π Y l m ( θ ) as shown in Equation (3), and

R k ( N N 1 , N ′ N 1 ′ ) = e 2 κ ∫ ∫ R < k R > k + 1 R N 1 ′ ( R 1 ) R N 1 ( R 1 ) × R N ′ ( R 2 ) R N ( R 2 ) R 1 2 R 2 2 d R 1 d R 2 ,

where R _< (R _>) is the smaller (bigger) value of {R ₁,R ₂}. In order that c ^k can have a non-zero value, k must satisfy the conditions

k + l N + l N ′ = 2 g ( g is an integer ) | l N − l N ′ | ≤ k ≤ | l N + l N ′ | .

Table 1 gives the values of c k for s and p electrons in case of m = 0.

Table 1

l _N

l N ′

m _N

m N ′

c ⁰( N ^′, N )

c ¹( N ^′, N )

c ²( N ^′, N )

l _N and m _N are angular momentum quantum number and magnetic quantum number for a quantum state N, respectively, whereas l N ′ and m N ′ for a quantum state N ^′. c ⁰(N ^′,N), c ¹(N ^′,N), and c ²(N ^′,N) are angle factors for k = 0 , 1 , 2 defined by Equations (13) and (15).

Plasmon and surface-plasmon modes

From electron energy spectrum obtained from the solution of the Schrödinger equation, we can derive the retarded and advanced Green’s function for electrons. Applying these Green’s functions and V _αβ with β =N ^′ N to the diagrammatic techniques to derive effective e-e interaction under the random phase approximation, the element of the dielectric function matrix is obtained as 17 18

ε α β ( Ω ) = δ α , β − V α β π β ( Ω ) ,

where

π N ′ N ( Ω ) = g s [ f ( E N ′ ) − f ( E N ) ] ℏ Ω + E N ′ − E N + i δ

is the pair bubble (or density-density correlation function) in the absence of e-e coupling with g _s= 2, counting for spin degeneracy, and f(E _N) = [ 1 + e ( E N − E F ) / k B T ] − 1 being the Fermi-Dirac function.

In a hollow nanosphere system described by quantum number N = (n l m), the electronic subband energy depends only on n and l quantum numbers, namely E _N= E _nl. In this study, we consider that n = 1 states with many l and m numbers are occupied, and n = 2 states with any l and m numbers are unoccupied. For simplicity, we take a four-state model (FSM) to calculate the dielectric function matrix. We consider that two lowest electronic states for n = 1, E₁₀, and E₁₁ are occupied, and two lowest electronic states for n = 2, E₂₀, and E₂₁ are unoccupied, as shown in Figure 1. Because the electronic subband energy in a hollow sphere does not depend on the quantum number m, we take m = 0 in the calculations. On the basis that all electronic states in a hollow nanosphere are quantized, intra-subband transitions do not contribute to dielectric function. Moreover, the transitions within the occupied and within the unoccupied states do not contribute to the dielectric function as well. Thus, as shown in Figure 1, there are eight possible transition channels induced by inter-subband transitions from occupied (unoccupied) states to unoccupied (occupied) states in this FSM. Setting the electronic state index as 1 = (100), 2 = (110), 3 = (200), and 4 = (210), the dielectric function of a hollow nanosphere in the FSM is a 16×16 matrix and can be obtained from Equation (16). The determinant of the dielectric function matrix then is

| ε ( Ω ) | = [ ( 1 + a 3131 + a 3113 ) ( 1 + a 4242 + a 4224 ) − ( a 3142 + a 3124 ) ( a 4231 + a 4213 ) ] × [ ( 1 + a 4141 + a 4114 ) ( 1 + a 3232 + a 3223 ) − ( a 3241 + a 3214 ) ( a 4132 + a 4123 ) ] ,

with a _αβ= −V _αβ π _β(Ω). The plasmon and surface-plasmon modes are determined by Re|ε(Ω)|→0 and Re|ε(Ω)|→−1, respectively. In this study, we employ a matrix to present the dielectric function in a multi-energy level system such as a hollow nanosphere structure. Such an approach was applied to study the plasmon excitations in semiconductor-based two-dimensional electron gas systems 19 and Rashba spintronic systems 18 . We note that in the present study, we consider a simple model to calculate the electronic subband structure of a hollow nanosphere. The effect of the spin-orbit interaction in the system is not included.

Figure 1

The four-state model for electronic transitions in a hollow nanosphere

The four-state model for electronic transitions in a hollow nanosphere. E₁₀ and E₁₁ are occupied states, and E₂₀ and E₂₁ are unoccupied states. E_F is the Fermi energy. Here, we consider the electronic states in the case of m = 0. The possible transition channels are indicated.

Results and discussion

In numerical calculations, we take the effective mass for an electron to be about the rest electron mass, i.e., μ ≃ 0.99m _e, and the high-frequency dielectric constant κ = 1.53 for gold shell 20 21 22 . In Figure 2, we show the electronic subband energy for E1_l and E2_l(l = 0, 1, and 2) as a function of outer radius r ₂ of the hollow nanosphere with a fixed shell thickness d = 10 nm. We see that the energy levels with different l quantum numbers roughly degenerate when r ₂ > 100 nm. In such a case, the subband energy depends very little on r ₂and E _nl≃ E _n0 = ℏ ² Π ² n ²/2μ d ². In Figure 3a, the electronic subband energies for E1_land E2_l(l = 0, 1, and 2) are shown functions of shell thickness d at a fixed outer radius r ₂ = 100 nm of the hollow nanosphere. The results for different l states coincide roughly. The subband energy decreases with increasing shell thickness as E _nl≃ E _n0∼d ⁻². In Figure 3b, the electronic subband energies for E1_land E2_l(l = 0, 1, and 2) are shown functions of shell thickness d for a fixed outer radius r ₂ = 25 nm of the hollow nanosphere. The subband energies degenerate roughly at small shell thickness and show difference with increasing shell thickness d as shown in the inset in Figure 3b. We know that the energy for the n ^th subband at l = 0 is determined only by the shell thickness d = r ₂ − r ₁ of a hollow nanosphere. The results shown in Figures 2 and 3 indicate that different l states degenerate at a fixed n quantum number when r ₂ > 100 nm. This feature is mainly induced by the symmetry of the confining potential for electrons, given as Equation (1). However, when r ₂ is relatively small (see Figure 2 and Figure 3b), the electronic subband energy depends on quantum number l for a fixed quantum number n and E _nl>E _n0 (here l > 0). This suggests that the stronger quantum effect can be achieved in smaller sample structures. Such an effect can be understood by the fact that when r ₂→∞, the energies determined by Equations (8) and (10) approach E _nl→ E _n0 = ℏ ² Π ² n ²/2μ d ² and when r ₂ takes a finite value E _nl> E _n0.

Figure 2

The subband energies vary with r₂ at a fixed d

The subband energies vary with r₂ at a fixed d. Electronic subband energy, E1_land E2_lfor l = 0, 1, and 2, as a function of outer radius r₂ of hollow nanosphere at a fixed shell thickness d = 10 nm.

Figure 3

The subband energies vary with d at fixed r₂

The subband energies vary with d at fixed r₂. Electronic subband energy, E1_land E2_lfor l = 0, 1, and 2, as a function shell thickness d of hollow nanosphere for outer radius r₂ = 100 nm in (a) and r₂ = 25 nm in (b). The inset in (b) shows the energy difference in different l states.

In Figure 4, the plasmon and surface-plasmon frequencies of hollow nanosphere are shown as a function of outer radius r ₂ at a fixed shell thickness d. Using the FSM, there are four modes for both plasmon and surface-plasmon excitation from a hollow nanosphere. We see that (1) the plasmon and surface-plasmon frequencies decrease with increasing r ₂ when r ₂ < 200 nm. When r ₂ > 200 nm, the plasmon and surface-plasmon frequencies approach approximately to the energy-gap between E ₂₀ and E ₁₀; (2) the frequencies of all these modes are in the THz regime; (3) the plasmon frequency Ω 1 p ≃ 2 Ω 1 s with Ω 1 s being a surface-plasmon frequency. This is the primary relationship between plasmon and surface-plasmon modes; (4) the surface-plasmon frequency Ω 2 s is slightly higher than plasmon frequency Ω 2 p ; and (5) Ω 4 p ≃ Ω 4 s and Ω 3 p ≃ Ω 3 s .

Figure 4

Dependence of plasmon and surface-plasmon frequencies on outer radius r₂ at a fixed shell thickness d = 10 nm

Dependence of plasmon and surface-plasmon frequencies on outer radius r ₂ at a fixed shell thickness d = 10 nm.

In Figure 5, the plasmon and surface-plasmon frequencies are shown as functions of shell thickness d at a fixed r ₂. The plasmon and surface-plasmon frequencies decrease with increasing shell thickness. The frequency difference between different excitation modes gets wider with increasing d.

Figure 5

Dependence of plasmon and surface-plasmon frequencies on the shell thickness d at a fixed outer radius r₂ = 100nm

Dependence of plasmon and surface-plasmon frequencies on the shell thickness d at a fixed outer radius r ₂ = 100 nm.

It should be noted that when r ₂ > 100 nm and d ∼ 10 nm, E _nl≃ E _n0, and plasmon and surface-plasmon frequencies in a hollow nanosphere are determined mainly by transition events between E _2land E _1l. This implies that although only four electronic subbands are included within current calculations, the obtained results should be very much similar to the case where more electronic states are considered when r ₂ > 100 nm and d ∼ 10 nm. The results obtained from this study indicate that the electronic subband energy and the plasmon and surface-plasmon modes in hollow nanospheres are determined mainly by sample parameters such as the diameter of the sphere r ₂ and the shell thickness d. When r ₂ > 100 nm, the energy levels depend very weakly on inner or outer radius (i.e., r ₁ or r ₂) at a fixed d. Thus, the shell thickness affects more strongly the electronic subband energies in a hollow nanosphere. We find that when d ∼ 10 nm and r ₂ ≥ 100 nm, the energy spacing between E2_land E1_l states is about 10 meV or about 2.4 THz. The frequencies of plasmon and surface-plasmon modes in the structure are also in the THz bandwidth. The plasmon and surface-plasmon modes depend sensitively on the geometrical parameters such as the outer radius r ₂ and shell thickness d. These effects imply that metal-based hollow nanosphere structures can be applied as THz materials or devices in which THz optical absorption and excitation can be achieved via inter-subband electronic transitions. It is known that THz technology is of great potential to impact many interdisciplinary fields such as telecommunication, biological science, pharmaceutical technology, anti-terrorist, etc. 23 . The application of nanostructure in THz technology has become a fast growing field of research in recent years. The theoretical findings from this work confirm that hollow gold-nanosphere structures are indeed the THz plasmonic materials which can be applied as frequency-tunable THz optoelectronic devices.

Conclusions

In this study, we have examined theoretically the electronic subband structure and the plasmon and surface-plasmon modes of hollow nanosphere structures. We have found that when the diameter of the sphere r ₂ > 100 nm and the shell thickness d ∼ 10 nm, the energy levels for different l states roughly degenerate. In such a case, the electronic subband energy, E _nl≃ E _n0 = ℏ ² Π ² n ²/2μ d ², does not depend on r ₂. When r ₂ < 200 nm, the plasmon and surface-plasmon modes induced by different electronic transition channels have significantly different frequencies. When r ₂ > 200 nm, the plasmon and surface-plasmon frequencies approach roughly to Ω ^p∼ Ω ^s∼ (E ₂₀ − E ₁₀)/ℏ, which depend largely on d and depend very little on r ₂.

It should be noted that at present, little research work has been carried out to look into the electronic subband structure of the hollow nanosphere structures using more powerful theoretical tools such as the first principle calculations which require large scale numerical computations and are CPU-consuming. The simple analytical results obtained from this study can be applied further to study the electronic and optoelectronic properties of the hollow nanosphere structures. We have found that the plasmon and surface-plasmon excitations can be achieved via inter-subband electronic transition channels in the hollow nanospheres. In particular, we have demonstrated that in metal hollow nanospheres, the energy difference between E1_l and E2_l states, and the plasmon and surface-plasmon frequencies are all in the THz bandwidth. This can lead to an application of metal hollow nanosphere structures in THz optics and optoelectronics.

Competing interests

The authors declare that they have no competing interests.

Author’s contributions

WX proposed and supervised the research work. YX carried out the analytical and numerical calculations. YZ and JH participated in the discussions and analyzes of the obtained results. All authors read and approved the final manuscript.

Author’s information

WX is the distinguished professor at Yunnan University and research professor at the Institute of Solid State Physics, Chinese Academy of Sciences. YX and YZ are post-graduate students at Yunnan University. JH is a PhD student at Yunnan University.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (grant no.: 10974206), the Ministry of Science and Technology of China (grant no.: 2011YQ130018), the Department of Science and Technology of Yunnan Province, and by the Chinese Academy of Sciences.

Template synthesis and photocatalytic properties of porous metal oxide spheres formed by nanoparticle infiltrationShchukinDGCarusoRAChem Mater200416228710.1021/cm0497780Stimuli-responsive controlled drug release from a hollow mesoporous silica sphere/polyelectrolyte multilayer core-shell structureZhuYFShiJLShenWHDongXPFengJWRuanMLLiYSAngew Chem Int Ed200544508310.1002/anie.200501500Preparation and characterization of porous hollow silica nanoparticles for drug delivery applicationChenJFDingHMWangJXShaoLBiomaterials20042572310.1016/S0142-9612(03)00566-014607511Smart inorganic/organic nanocomposite hollow microcapsulesShchukinDGSukhorukovGBMöhwaldHAngew Chem Int Ed200342447210.1002/anie.200352068Pt hollow nanospheres: facile synthesis and enhanced electrocatalystsLiangHPZhangHMHuJSGuoYGWanLJBaiCLAngew Chem Int Ed200443154010.1002/anie.200352956A bis(p-sulfonatophenyl)phenylphosphine-based synthesis of hollow Pt nanospheresYangJLeeJYTooHPValiyaveettilSJ Phys Chem B200611012510.1021/jp055306c16471509Biologically erodable microspheres as potential oral drug delivery systemsMathiowitzEJacobJSJongYSCarinoGPChickeringDEChaturvediPSantosCAVijayaraghavanKMontgomerySBassettMMorrellCNature199738641010.1038/386410a09121559Hollow nanoparticles of WS2 as potential solid-state lubricantsRapoportLBilikYFeldmanYHomyonferMCohenSRTenneRNature199738779110.1038/42910Permeable silica shell through surface-protected etchingZhangQZhangTRGeJPYinYDNano Lett20088286710.1021/nl801618718698725A hybridization model for the plasmon response of complex nanostructuresProdanERadloffCHalasNJNordlanderPScience200330241910.1126/science.108917114564001Formation of hollow nanocrystals through the nanoscale Kirkendall effectYinYDRiouxRMErdonmezCKHughesSSomorjaiGAAlivisatosAPScience200430471110.1126/science.109656615118156Nanoengineering of inorganic and hybrid hollow spheres by colloidal templatingCarusoFCarusoRAMöhwaldHScience199828211119804547General synthesis of 2D ordered hollow sphere arrays based on nonshadow deposition dominated colloidal lithographyDuanGTLvFJCaiWPLuoYYLiYLiuGQLangmuir201026629510.1021/la904116p20131831Visible to infrared photoluminescence from gold nanoparticles embedded in germano-silicate glass fiberLinASonDHAhnIHSongGHHanWTOptics Express200715637410.1364/OE.15.00637419546942THz generation from plasmonic nanoparticle arraysPolyushkinDKHendryEStoneEKBarnesWLNano Lett201111471810.1021/nl202428g22007706Hollow micro-/nanostructures: synthesis and applicationsLouXWArcherLAYangZAdv Mater200820398710.1002/adma.200800854CondonEShortleyGHThe Theory of Atomic SpectraLondon: Cambridge University Press1959,p. 174Plasmons of a two-dimensional electron gas in the presence of spin orbit interactionXuWAppl Phys Lett20038272410.1063/1.1541098Collective modes of spatially separated, two-component, two-dimensional plasma in solidsDas SarmaSMadhukarAPhys Rev B19812380510.1103/PhysRevB.23.805A new simple method of determining the effective mass of an electron or the thickness of thin metal filmsSzczyrbowskiJJ Phys D: Appl Phys198619125710.1088/0022-3727/19/7/015Optical constants of the noble metalsJohnsonPBChristyRWPhys Rev B19726437010.1103/PhysRevB.6.4370An analytic model for the optical properties of goldEtchegoinPGLe RuECMeyerMJ Chem Phys200612516470510.1063/1.236027017092118Terahertz technologySiegelPHIEEE Trans Microwave Theory Tech20025091010.1109/22.989974