Next Article in Journal
Highly Porous Au–Pt Bimetallic Urchin-Like Nanocrystals for Efficient Electrochemical Methanol Oxidation
Previous Article in Journal
Dynamic Manipulation of THz Waves Enabled by Phase-Transition VO2 Thin Film
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Light Scattering from Rough Silver Surfaces: Modeling of Absorption Loss Measurements

1
Department of Physics, University of Kashan, 8731753153 Kashan, Iran
2
Institute of Condensed Matter Theory and Optics, Abbe Center of Photonics, Friedrich-Schiller-Universität Jena, 07743 Jena, Germany
*
Author to whom correspondence should be addressed.
Nanomaterials 2021, 11(1), 113; https://doi.org/10.3390/nano11010113
Submission received: 8 December 2020 / Revised: 23 December 2020 / Accepted: 4 January 2021 / Published: 6 January 2021
(This article belongs to the Special Issue Plasmonics and Related Nanomaterials)

Abstract

:
We consider two series of experimental setups of multilayered Ag/ZnO thin films with varying surface morphologies given by atomic force microscopy images. The absorption loss under diffuse scattering is studied theoretically by applying a combination of the scattering matrix approach with diffraction theory for randomly nanotextured interfaces. Our modeling is in excellent agreement with the respective measurements. The theoretical approach is applicable to a wide range of wavelengths, surface morphologies, and materials for both measured and computed rough surface morphologies.

Graphical Abstract

1. Introduction

Regular nanostructures, such as diffraction gratings [1], photonic crystals [2,3], and metasurfaces [4,5], allow tuning of their optical response to the desired spectral properties, e.g., in nanoparticle and hole arrays [6,7,8,9,10], as well as their near-field patterns under varying illumination conditions, e.g., for optical data storage [11,12]. Fabrication techniques range from chemical [13] to physical methods, such as self-assembly [14], laser-induced [15], and annealing [16,17] procedures. This allows a high geometrical control over the produced structures.
Photocatalytic [18,19,20], photovoltaic [21,22], and some biochemical sensing applications, such as Raman spectroscopy [3,23,24], on the other hand, rely on large surface or contact areas. In such systems, strong local field enhancement or enhanced forward scattering across the entire surface is desired. Likewise, a broad spectral response, rather than specific spectral lines created at distinct hot-spot positions, is of importance. Rough surface morphologies address this need and are, thus, of high interest in a number of such applications. Typically, they are easier and cheaper to produce, which makes them important in large-scale industrial applications, such as solar [25,26] and fuel cells compatible with mass production lines [27]. Moreover, in realistic setups, surface roughness can often not be neglected, and diffuse scattering with possibly inherent additional losses should be taken into account when considering the optical response of any nanostructure.
In this article, we consider the two Ag/ZnO samples depicted in Figure 1, which were investigated experimentally with photothermal deflection spectroscopy by Springer et al., Ref. [28], as potential back reflectors for solar cells. The experimental setups consisted of a 500 nm ZnO film on a glass substrate supporting a silver layer with varying surface roughness in a CCl 4 environment. These comprise the first sample type in Figure 1. In a second step, an additional ZnO film of 30 nm was coated onto the original samples to comprise a second batch, which we denote as sample type 2 in accordance with the original experimental study [28]. Experimentally, an increase in absorption was observed, which was related to the surface plasmon absorption.
To simulate the absorption loss in these systems, we employ a theoretical model for light scattering from rough and nanotextured surfaces [29,30,31,32], which combines diffraction theory for large-scale irregular surfaces with the scattering matrix approach for planar multilayered systems [33,34]. The next section explains this approach in detail and applies it to the original atomic force microscopy (AFM) images of six different surface morphologies provided by Springer et al. [28]. Table 1 lists the root mean square (rms) roughness values for the samples used. Our simulations show excellent agreement with the measured absorption loss. This approach enables us to consider further material combinations for the fabricated samples, as well as arbitrary surface morphologies beyond the ones provided experimentally.

2. Methodology

The frequency-dependent permittivity function
ϵ ( ω ) = ϵ b γ σ ω ( ω + i γ ) + j = 1 L i σ j ω Σ j + i σ j * ω + Σ j *
is modeled based on the Drude–Lorentz theory with the large frequency (background) permittivity ϵ b , plasmonic damping rate γ , the electron conductivity σ , generalized conductivities σ j of bound electrons, and their respective Lorentz pole resonances Ω j . These parameters are taken from numerical fits to experimental data [35] for silver and other metals.
The basis for the description of the rough surfaces comprises AFM data spanning 10 μ m × 10 μ m in physical size. Figure 2a shows three of the six available data sets with their respective height ranges. The elements of the height function z ( x , y ) for each surface are related to the pixel values in the grayscale image and are related to the height range in nanometers using the AFM data. For instance, the first surface in Figure 2a has a maximum height value of 80 nm (255 in grayscale), and the minimum height value is −52 nm (0 in grayscale). The AFM images thus define the different height functions z ( x , y ) for the investigated setups.
Following Refs. [31,32], we calculate the reflected ( ν = R ) and transmitted ( ν = T ) scalar electric fields U ν , which allow analysis of the total reflection and transmission in terms of specular and diffuse scattering. Hereby, the diffuse scattering contribution is given by the haze function H ν and can be considered separately. It is defined as the (normalized) difference between the total reflection or transmission and the specular contribution (here at normal incidence):
H ν ( λ ) = ν tot ν spec ν tot .
The scalar fields are
U ν ( k x , k y ) = 1 2 π R 2 G ν ( x , y ) e i ( k x x + k y y ) d x d y ,
using the pupil functions
G T ( x , y ) = T flat N e i k 0 z ( x , y ) ( n R n T ) , G R ( x , y ) = R flat N e i k 0 z ( x , y ) 2 n R ,
in two-dimensional k-space depending on the total transmission T flat and reflection R flat of the flat geometry. Hereby, n R is the refractive index of either the CCl 4 environment ( n R = 1.4 for type 1 samples) or the additional ZnO layer ( n R = 2.1 for sample type 2), and n T is the refractive index of the silver layer. Hence, the rough interface is always the top of the silver layer. Using fast Fourier transform (FFT) analysis, the results of these integrals are depicted in Figure 2b for U R for three data sets, indicating the plethora of diffraction orders and scattering angles relevant to describing the optical response from such surfaces as compared to regular structures. Rather than as a continuous integral, the height function is defined via N pixels in the image, which is used here as a normalization. Hence, k x and k y vector components for each diffracted beam can be written in terms of the pixel position as 2 π m P , in which P is a number from 0 to m 1 and m × m = N is the dimension of the image ( 256 × 256 pixels). The free-space wavenumber is k 0 = M M 2 π / λ , so that k 0 2 n ν 2 = k x 2 + k y 2 + k z 2 always. The red circle indicates where k x 2 + k y 2 = k 0 2 n ν 2 , and thus marks the border between regimes where k z becomes real (propagating modes for k 0 2 n ν 2 k x 2 k y 2 > 0 ) or imaginary (evanescent modes for k 0 2 n ν 2 k x 2 k y 2 < 0 ). Evanescent contributions lie, thus, outside of this circle, and are excluded in the following considerations, which is indicated by the truncated sum.
As seen in Equation (4), the total reflection R flat and transmission T flat of the flat interface are needed as inputs. T flat and R flat are calculated in a first step for the planar multilayer setup using the scattering matrix method [33,34]; i.e., the ideal planar structures from Figure 1a are considered, yielding the specular reflection and zeroth-order transmission values. Afterwards, the randomly nanotextured interfaces are accounted for. The deviation from the ideal result is then calculated by effective redistribution for diffuse scattering angles.
First, our simulation model processes the AFM images of the rough surface topography and extracts their two-dimensional height functions z ( x , y ) and rms roughness. We evaluated the rms roughness with minor deviations from the original work (see Table 2). In the numerical analysis of the AFM samples, some defects were cut from the input images of the height function to obtain surface roughness parameters close to the original evaluation of Ref. [28]. For this purpose, we cropped the affected AFM images and then resized them to the original numbers of rows and columns ( 256 × 256 ) using a bicubic interpolation algorithm. Furthermore, the obtained rms roughness was compared using Gwyddion [36], which specializes in data visualization and analysis for scanning probe microscopy images. Second, we then calculated the pupil functions G ν ( x , y ) and two-dimensional Fourier transforms of each image, which led to the electric scalar fields U ν ( k x , k y ) for reflected and transmitted light.
It should be noted that the diffraction theory incorporated here is based on the paraxial approximation, and the diffraction at the rough interfaces results in a phase shift only. Both these restrictions apply when nanosized roughness is considered, but may break down for larger rms roughness. Moreover, we considered a single rough interface contributing to diffuse scattering, and no pathway of actual scattered beams is tracked, i.e., there is no interference between diffusely scattered beams, since the multilayer result of the flat reference enters the calculation. A more advanced approach has to improve on these aspects.
The total optical response, i.e., including diffuse and specular contributions to the reflection ( ν = R ) and transmission ( ν = T ), is determined from summing up all contributions within the unit circle (depicted in red in Figure 2b) where propagating modes are found:
ν = k x 2 + k y 2 k 0 2 n ν 2 | U ν ( k x , k y ) | 2 ,
from which the absorptance is calculated in the usual way:
A = 1 T R .
The absorption loss as a function of the rms roughness is the main quantity considered in this work.
With these considerations, the diffuse scattering contribution, i.e., the haze function, results in
H ν ( λ ) = k x 2 + k y 2 k 0 2 n ν 2 | U ν ( k x , k y ) | 2 | U ν ( 0 , 0 ) | 2 k x 2 + k y 2 k 0 2 n ν 2 | U ν ( k x , k y ) | 2 .

3. Results and Discussion

The haze in transmission for different surface morphologies is depicted in Figure 3 for the sample type 2 results. As can be expected, the largest haze is found for the highest rms roughness. It dominates the optical response over the specular contribution. The haze follows the overall form discussed by Jäger et al. [32]; however, the haze in transmission is strongly influenced by the additional metal film in the layered system and its surface plasmon polariton. A noticeable reduction in the transmission is seen around this resonance at about 315 nm, where the light is rather absorbed than transmitted. Furthermore, note the clear sequence of the samples: With increasing rms roughness, the haze strongly increases, almost reaching its maximum for the three largest rms roughnesses over the entire spectrum.
The overall absorption is shown in Figure 4 for the different rms values. Another important input value to consider is the layer thickness. Due to the rough morphology, the layer thickness of the central Ag layer is not a well-defined quantity. In the discussed experimental setup [28], the target thickness for the silver layer is 500 nm before increasing the surface roughness. However, this leaves an uncertainty about the final average layer thickness.
In the presented results, we have taken some liberty regarding the Ag layer thickness to find the best fit of the absorption loss data at different roughnesses (see Table 2). However, keeping the Ag layer film thickness constant for all sample types also well reproduces the measurements (see Figure 4b in comparison to Figure 4a). The final results are not altered by more than a percentage point. Overall, the curves are smoother, as the layers are too thick to show Fabry–Pérot-type oscillations, and absorption is larger as compared to the modified calculations, where we have substantially reduced the target thickness of the silver layer.
With the steps described above, we arrive at the absorption loss (shown in Figure 5 as a function of wavelength), which was calculated from Equations (5) and (6). The plasmon peak for silver is observed at a wavelength of about 315 nm (see the inset). Oscillations in the absorptance A with the wavelength, particularly for lower rms roughness, are due to light interferences in the thin film’s layered structure that has not yet been obscured by diffuse scattering from the rough surface. This is a direct result of the input of optical coefficients from the ideal flat structure, which accounts for multiple scattering effects.
In Figure 5a,b, we evaluate the absorption loss for the two sample types as a function of the rms roughness of the six samples at two distinct wavelengths. This is directly compared with the experimental results by Springer et al. [28], which are included as blue symbols. Excellent agreement with the measurements is found due to the inclusion of diffuse scattering.
In Figure 5c,d, a third-degree polynomial curve was fitted to the simulated absorptance for different rms to obtain a continuous function of the surface roughness for every wavelength between 200 to 1000 nm for the two sample types. The resulting contour plot allows the assessment of the absorption loss in the different layered structures, but also indicates a maximum of absorption for a specific rms roughness, as well as a best choice of rms roughness for broadband absorption. Hence, this simulation scheme allows specification of some characteristics of the samples for specific properties regarding the absorptance or transmission and reflection, respectively. Figure 5e zooms in on the plasmon peak, which remains stable for different values of rms roughness, at least within classical diffuse scattering theory. It could be expected that nano-scale features at the surface do not only contribute to diffuse scattering, but also to quantum phenomena, such as nonlocal optical response, particularly when considering metal surfaces [37,38,39]. This is already well understood with nanoparticles [40] and their impact on solar cells employing nanostructures [41,42], but not yet for rough surface morphologies, i.e., randomly sized and distributed features. Future work should consider such effects beyond classical diffraction optics.
The applicability of our modeling to further material combinations, different multilayer structures, and both experimental and numerical surface morphology samples is straightforward. Figure 6 shows, in analogy to Figure 5 for silver, the corresponding results if the silver film was replaced by gold or aluminum. First of all, it can be seen that the absorption loss with rough gold surfaces is similar to that of silver, with the difference that the plasmon peak is around 520 nm (see Figure 6c). Hence, the absorption is much larger over a broader spectral range, which can be more easily compared using spectrally integrated quantities, such as the overall photocurrent. In contrast, aluminum shows much larger absorption loss for the rms roughness given by AFM data, although its plasmon resonance lies in the ultraviolet range.

4. Conclusions

We investigated the optical properties of a multilayer containing a single rough Ag interface using scalar diffraction theory for randomly nanotextured thin films combined with the scattering matrix theory. First, the optical coefficients of the flat multilayer system are computed. Second, the surface of the top Ag layer is replaced with rough morphology data from AFM images using diffraction theory to calculate the diffuse scattering from these interfaces. With this, the absorption loss is calculated and compared to the measurements performed on the original samples.
The calculated absorption loss in Figure 5a,b agrees very well with the experimental results of Ref. [28].
This procedure can be applied to investigate any types of rough surface morphologies provided from fabricated samples via AFM or scanning electron microscopy (SEM) images or from numerically produced samples. It also enables us to investigate the effects of various shapes of nanoparticles and nanostructured surfaces using grayscale images for the height function within the approximations discussed earlier. This could be an interesting path towards the quick assessment of amorphous structures or smooth boundaries that are not accessible with standard tools.

Author Contributions

Investigation, M.D.; Methodology, M.D. and C.D.; Validation, M.D. and C.D.; Visualization, M.D.; Writing—original draft, C.D. and M.D.; Writing—review and editing, C.D. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by the Thüringer Universitäts- und Landesbibliothek Jena (ThULB).

Data Availability Statement

The datasets generated and analyzed during the presented study are available from the corresponding author on reasonable request.

Acknowledgments

The authors are grateful to Springer et al. [28] for providing their original AFM data for further analysis—in particular, O. Kluth and B. Rech.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chong, T.K.; Wilson, J.; Mokkapati, S.; Catchpole, K.R. Optimal wavelength scale diffraction gratings for light trapping in solar cells. J. Opt. 2012, 14, 024012. [Google Scholar] [CrossRef] [Green Version]
  2. Joannopoulos, J.D.; Johnson, S.G.; Winn, J.N.; Meade, R.D. Photonic Crystals: Molding the Flow of Light, 2nd ed.; Princeton University Press: Princeton, NJ, USA, 2008. [Google Scholar]
  3. David, C. TiO2 Self-Assembled, Thin-Walled Nanotube Arrays for Photonic Applications. Materials 2019, 12, 1332. [Google Scholar] [CrossRef] [Green Version]
  4. Jakšić, Z.; Vuković, S.; Matovic, J.; Tanasković, D. Negative Refractive Index Metasurfaces for Enhanced Biosensing. Materials 2010, 4, 1. [Google Scholar] [CrossRef] [Green Version]
  5. Yan, J.; Liu, X.; Ma, C.; Huang, Y.; Yang, G. All-dielectric materials and related nanophotonic applications. Mater. Sci. Eng. R Rep. 2020, 141, 100563. [Google Scholar] [CrossRef]
  6. Genet, C.; Ebbesen, T.W. Light in tiny holes. Nature 2007, 445, 39–46. [Google Scholar] [CrossRef]
  7. Jonsson, M.P.; Dahlin, A.B.; Jönsson, P.; Höök, F. Nanoplasmonic biosensing with focus on short-range ordered nanoholes in thin metal films. Biointerphases 2008, 3, FD30–FD40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. García de Abajo, F.J. Optical excitations in electron microscopy. Rev. Mod. Phys. 2010, 82, 209–275. [Google Scholar] [CrossRef]
  9. David, C.; Christensen, J.; Mortensen, N.A. Spatial dispersion in two-dimensional plasmonic crystals: Large blueshifts promoted by diffraction anomalies. Phys. Rev. B 2016, 94, 165410. [Google Scholar] [CrossRef] [Green Version]
  10. David, C.; Christensen, J. Extraordinary optical transmission through nonlocal holey metal films. Appl. Phys. Lett. 2017, 110, 261110. [Google Scholar] [CrossRef] [Green Version]
  11. Morarescu, R.; Englert, L.; Kolaric, B.; Damman, P.; Vallee, R.A.L.; Baumert, T.; Hubenthal, F.; Trager, F. Tuning nanopatterns on fused silica substrates: A theoretical and experimental approach. J. Mater. Chem. 2011, 21, 4076–4081. [Google Scholar] [CrossRef]
  12. David, C.; Kühler, P.; Javier García de Abajo, F.; Siegel, J. Near-field nanoimprinting using colloidal monolayers. Opt. Express 2014, 22, 8226–8233. [Google Scholar] [CrossRef] [PubMed]
  13. Kumar, P.S.; Pastoriza-Santos, I.; Rodríguez-González, B.; de Abajo, F.J.G.; Liz-Marzán, L.M. High-yield synthesis and optical response of gold nanostars. Nanotechnology 2008, 19, 015606. [Google Scholar] [CrossRef]
  14. Sarkar, S.; Gupta, V.; Kumar, M.; Schubert, J.; Probst, P.T.; Joseph, J.; König, T.A. Hybridized Guided-Mode Resonances via Colloidal Plasmonic Self-Assembled Grating. ACS Appl. Mater. Interfaces 2019, 11, 13752–13760. [Google Scholar] [CrossRef] [Green Version]
  15. Li, L.; Guo, W.; Wang, Z.B.; Liu, Z.; Whitehead, D.; Luk’yanchuk, B. Large-area laser nano-texturing with user-defined patterns. J. Micromech. Microeng. 2009, 19, 054002. [Google Scholar] [CrossRef]
  16. Roy, P.; Berger, S.; Schmuki, P. TiO2 Nanotubes: Synthesis and Applications. Angew. Chem. Int. Ed. 2011, 50, 2904–2939. [Google Scholar] [CrossRef] [PubMed]
  17. Öner, I.H.; Querebillo, C.J.; David, C.; Gernert, U.; Walter, C.; Driess, M.; Leimkühler, S.; Ly, K.H.; Weidinger, I.M. High Electromagnetic Field Enhancement of TiO2 Nanotube Electrodes. Angew. Chem. Int. Ed. 2018, 57, 7225–7229. [Google Scholar] [CrossRef] [PubMed]
  18. Akimov, A.V.; Neukirch, A.J.; Prezhdo, O.V. Theoretical Insights into Photoinduced Charge Transfer and Catalysis at Oxide Interfaces. Chem. Rev. 2013, 113, 4496–4565. [Google Scholar] [CrossRef] [PubMed]
  19. Ma, X.C.; Dai, Y.; Yu, L.; Huang, B.B. Energy transfer in plasmonic photocatalytic composites. Light Sci. Appl. 2016, 5, e16017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. David, C. Two-fluid, hydrodynamic model for spherical electrolyte systems. Sci. Rep. 2018, 8, 7544. [Google Scholar] [CrossRef]
  21. Wang, E.; White, T.P.; Catchpole, K.R. Resonant enhancement of dielectric and metal nanoparticle arrays for light trapping in solar cells. Opt. Express 2012, 20, 13226–13237. [Google Scholar] [CrossRef]
  22. Villanueva-Cab, J.; Montaño-Priede, J.L.; Pal, U. Effects of Plasmonic Nanoparticle Incorporation on Electrodynamics and Photovoltaic Performance of Dye Sensitized Solar Cells. J. Phys. Chem. C 2016, 120, 10129–10136. [Google Scholar] [CrossRef]
  23. Feng, J.J.; Gernert, U.; Sezer, M.; Kuhlmann, U.; Murgida, D.H.; David, C.; Richter, M.; Knorr, A.; Hildebrandt, P.; Weidinger, I.M. Novel Ag-Au Hybrid Device for Electrochemical SE(R)R Spectroscopy in a Wide Potential and Spectral Range. Nano Lett. 2009, 9, 298–303. [Google Scholar] [CrossRef]
  24. David, C.; Richter, M.; Knorr, A.; Weidinger, I.M.; Hildebrandt, P. Image dipoles approach to the local field enhancement in nanostructured Ag-Au hybrid devices. J. Chem. Phys. 2010, 132, 024712. [Google Scholar] [CrossRef] [PubMed]
  25. David, C.; Connolly, J.P.; Chaverri Ramos, C.; García de Abajo, F.J.; Sánchez Plaza, G. Theory of random nanoparticle layers in photovoltaic devices applied to self-aggregated metal samples. Sol. Energy Mater. Sol. Cells 2013, 109, 294–299. [Google Scholar] [CrossRef] [Green Version]
  26. David, C. Multi-type particle layer improved light trapping for photovoltaic applications. Appl. Opt. 2016, 55, 7980–7986. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Hamed, T.A.; Adamovic, N.; Aeberhard, U.; Alonso-Alvarez, D.; Amin-Akhlaghi, Z.; der Maur, M.A.; Beattie, N.; Bednar, N.; Berland, K.; Califano, M.; et al. Multiscale in modelling and validation for solar photovoltaics. EPJ Photovolt. 2018, 9, 10. [Google Scholar] [CrossRef]
  28. Springer, J.; Poruba, A.; Müllerova, L.; Vanecek, M.; Kluth, O.; Rech, B. Absorption loss at nanorough silver back reflector of thin-film silicon solar cells. J. Appl. Phys. 2004, 95, 1427–1429. [Google Scholar] [CrossRef] [Green Version]
  29. Jäger, K.; Zeman, M. A scattering model for surface-textured thin films. Appl. Phys. Lett. 2009, 95, 171108. [Google Scholar] [CrossRef] [Green Version]
  30. Battaglia, C.; Söderström, K.; Escarré, J.; Haug, F.J.; Dominé, D.; Cuony, P.; Boccard, M.; Bugnon, G.; Denizot, C.; Despeisse, M.; et al. Efficient light management scheme for thin film silicon solar cells via transparent random nanostructures fabricated by nanoimprinting. Appl. Phys. Lett. 2010, 96, 213504. [Google Scholar] [CrossRef] [Green Version]
  31. Dominé, D.; Haug, F.J.; Battaglia, C.; Ballif, C. Modeling of light scattering from micro- and nanotextured surfaces. J. Appl. Phys. 2010, 107, 044504. [Google Scholar] [CrossRef]
  32. Jäger, K.; Fischer, M.; van Swaaij, R.A.C.M.M.; Zeman, M. A scattering model for nano-textured interfaces and its application in opto-electrical simulations of thin-film silicon solar cells. J. Appl. Phys. 2012, 111, 083108. [Google Scholar] [CrossRef]
  33. Whittaker, D.M.; Culshaw, I.S. Scattering-matrix treatment of patterned multilayer photonic structures. Phys. Rev. B 1999, 60, 2610–2618. [Google Scholar] [CrossRef] [Green Version]
  34. Liu, V.; Fan, S. S4: A free electromagnetic solver for layered periodic structures. Comput. Phys. Commun. 2012, 183, 2233–2244. [Google Scholar] [CrossRef]
  35. Sehmi, H.S.; Langbein, W.; Muljarov, E.A. Optimizing the Drude-Lorentz model for material permittivity: Method, program, and examples for gold, silver, and copper. Phys. Rev. B 2017, 95. [Google Scholar] [CrossRef] [Green Version]
  36. Available online: http://gwyddion.net/ (accessed on 18 April 2020).
  37. David, C.; García de Abajo, F.J. Spatial Nonlocality in the Optical Response of Metal Nanoparticles. J. Phys. Chem. C 2011, 115, 19470–19475. [Google Scholar] [CrossRef]
  38. David, C.; Mortensen, N.A.; Christensen, J. Perfect imaging, epsilon-near zero phenomena and waveguiding in the scope of nonlocal effects. Sci. Rep. 2013, 3, 2526. [Google Scholar] [CrossRef] [Green Version]
  39. David, C.; García de Abajo, F.J. Surface Plasmon Dependence on the Electron Density Profile at Metal Surfaces. ACS Nano 2014, 8, 9558–9566. [Google Scholar] [CrossRef]
  40. Raza, S.; Wubs, M.; Bozhevolnyi, S.I.; Mortensen, N.A. Nonlocal study of ultimate plasmon hybridization. Opt. Lett. 2015, 40, 839–842. [Google Scholar] [CrossRef] [Green Version]
  41. Kluczyk, K.; Jacak, L.; Jacak, W.; David, C. Microscopic Electron Dynamics in Metal Nanoparticles for Photovoltaic Systems. Materials 2018, 11, 1077. [Google Scholar] [CrossRef] [Green Version]
  42. Kluczyk-Korch, K.; Jacak, L.; Jacak, W.A.; David, C. Mode Splitting Induced by Mesoscopic Electron Dynamics in Strongly Coupled Metal Nanoparticles on Dielectric Substrates. Nanomaterials 2019, 9, 1206. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Illustration of the two series of Ag/ZnO setups considered. (a) Multilayered structures without rough surfaces for planar reference calculations and (b) their corresponding rough top interfaces, where diffraction theory is used in combination with scattering matrix calculations.
Figure 1. Illustration of the two series of Ag/ZnO setups considered. (a) Multilayered structures without rough surfaces for planar reference calculations and (b) their corresponding rough top interfaces, where diffraction theory is used in combination with scattering matrix calculations.
Nanomaterials 11 00113 g001
Figure 2. (a) Height function z ( x , y ) extracted from 10 μ m × 10 μ m atomic force microscopy (AFM) images of rough Ag surfaces (Ref. [28]) with varying root mean square (rms) roughnesses of 8.86, 30.48, and 134.17 nm from left to right, respectively, and (b) the corresponding absolute values of their 2D Fourier components for the reflected scalar field (Equation (3) at 650 nm of wavelength). The red circle indicates the condition k x 2 + k y 2 = k 0 2 n 1 2 . Absolute values were scaled by a factor of 5 for better contrast, and images were resized to 128 × 128 pixels for clarity.
Figure 2. (a) Height function z ( x , y ) extracted from 10 μ m × 10 μ m atomic force microscopy (AFM) images of rough Ag surfaces (Ref. [28]) with varying root mean square (rms) roughnesses of 8.86, 30.48, and 134.17 nm from left to right, respectively, and (b) the corresponding absolute values of their 2D Fourier components for the reflected scalar field (Equation (3) at 650 nm of wavelength). The red circle indicates the condition k x 2 + k y 2 = k 0 2 n 1 2 . Absolute values were scaled by a factor of 5 for better contrast, and images were resized to 128 × 128 pixels for clarity.
Nanomaterials 11 00113 g002
Figure 3. The (normalized) contribution of diffuse scattering is given by the haze, here in transmission H T = M M ( T tot T spec ) / T tot , for different morphologies and rms roughnesses of the ZnO/Ag interface of sample type 2. The silver plasmon peak is observed around 315 nm.
Figure 3. The (normalized) contribution of diffuse scattering is given by the haze, here in transmission H T = M M ( T tot T spec ) / T tot , for different morphologies and rms roughnesses of the ZnO/Ag interface of sample type 2. The silver plasmon peak is observed around 315 nm.
Nanomaterials 11 00113 g003
Figure 4. (a) Absorptance A as a function of wavelength for sample type 2 with various rms roughnesses of the ZnO/Ag interface on a logarithmic scale, including the ideal flat surface. The inset shows the silver plasmon resonance around 315 nm. (b) Same as (a) for constant Ag layer thickness of 500 nm.
Figure 4. (a) Absorptance A as a function of wavelength for sample type 2 with various rms roughnesses of the ZnO/Ag interface on a logarithmic scale, including the ideal flat surface. The inset shows the silver plasmon resonance around 315 nm. (b) Same as (a) for constant Ag layer thickness of 500 nm.
Nanomaterials 11 00113 g004
Figure 5. Absorptance as a function of rms roughness of surfaces for sample type 1 (Ag, filled black symbols) and sample type 2 (ZnO/Ag, empty black symbols) at incident wavelengths of (a) 650 nm and (b) 920 nm, respectively. The corresponding experimental data points extracted from Ref. [28] are shown in blue. (c,d) Absorptance on a logarithmic scale as a function of rms roughness and wavelength for the two sample types using a third-order polynomial fitting curve. (e) The plasmon peak of silver is observed around 315 nm. Here, a linear scale was used for clarity.
Figure 5. Absorptance as a function of rms roughness of surfaces for sample type 1 (Ag, filled black symbols) and sample type 2 (ZnO/Ag, empty black symbols) at incident wavelengths of (a) 650 nm and (b) 920 nm, respectively. The corresponding experimental data points extracted from Ref. [28] are shown in blue. (c,d) Absorptance on a logarithmic scale as a function of rms roughness and wavelength for the two sample types using a third-order polynomial fitting curve. (e) The plasmon peak of silver is observed around 315 nm. Here, a linear scale was used for clarity.
Nanomaterials 11 00113 g005
Figure 6. Same as in Figure 5, replacing the material function of silver with the corresponding functions for (a,c) gold and (b,d) aluminum.
Figure 6. Same as in Figure 5, replacing the material function of silver with the corresponding functions for (a,c) gold and (b,d) aluminum.
Nanomaterials 11 00113 g006
Table 1. The surface roughness of the different samples used in the calculations and the respective thickness of the Ag layer that best reproduces the experimental result.
Table 1. The surface roughness of the different samples used in the calculations and the respective thickness of the Ag layer that best reproduces the experimental result.
Type 1 rms (nm)Type 2 rms (nm)Ag Thickness (nm)
8.868.8665
30.4830.4885
76.3473.51105
107.69110.48195
134.17130.19255
151.97155.81315
0 (flat)0 (flat)315
Table 2. Comparing the calculated rms roughness extracted from the original AFM images based on different numerical sources.
Table 2. Comparing the calculated rms roughness extracted from the original AFM images based on different numerical sources.
This WorkGwyddionRef. [28]
9.309.3414
30.4830.4437
76.3476.2174
107.69107.50106
134.17133.90131
151.97151.60146
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dehghani, M.; David, C. Light Scattering from Rough Silver Surfaces: Modeling of Absorption Loss Measurements. Nanomaterials 2021, 11, 113. https://doi.org/10.3390/nano11010113

AMA Style

Dehghani M, David C. Light Scattering from Rough Silver Surfaces: Modeling of Absorption Loss Measurements. Nanomaterials. 2021; 11(1):113. https://doi.org/10.3390/nano11010113

Chicago/Turabian Style

Dehghani, Matin, and Christin David. 2021. "Light Scattering from Rough Silver Surfaces: Modeling of Absorption Loss Measurements" Nanomaterials 11, no. 1: 113. https://doi.org/10.3390/nano11010113

APA Style

Dehghani, M., & David, C. (2021). Light Scattering from Rough Silver Surfaces: Modeling of Absorption Loss Measurements. Nanomaterials, 11(1), 113. https://doi.org/10.3390/nano11010113

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop