Open access peer-reviewed chapter

Synthesis of Curved Surface Plasmon Fields through Thin Metal Films in a Tandem Array

By Gabriel Martinez Niconoff, Marco Antonio Torres Rodriguez, Mayra Vargas Morales and Patricia Martinez Vara

Submitted: March 15th 2018Reviewed: October 9th 2018Published: November 12th 2018

DOI: 10.5772/intechopen.81931

Downloaded: 626


We describe the generation of plasmonic modes that propagate in a curved trajectory inducing magnetic properties. This is performed by masking a metal surface with two screens containing a randomly distributed set of holes that follow a Gaussian statistic. The diameter of the holes is less than the wavelength of the illuminating plane wave. By implementing scaling and rotations on each screen, we control the correlation trajectory and generate long-range curved plasmonic modes. Using the evanescent character of the electric field, the study is implemented for the transmission of a plasmonic mode propagating in a tandem array of thin metal films offering the possibility to generate localization effects.


  • plasmon mode
  • surface plasmon field
  • speckle
  • thin films
  • curved correlation trajectory

1. Introduction

During the last decade, the scientific community has shown an increasing interest in the models of plasmon fields due to their potential applications, which occur practically in all branches of science and technology. In the present study, we emphasize the analysis of correlation trajectories on a metal surface with random structure. The resulting model offers applications to development of nano-antennas having the possibility of a tunable bandwidth [1]. This type of structure has applications in the synthesis of new light sources and the control of magnetic effects [2]. The tunable effects are controlled with the curvature parameter having applications in surface-enhanced Raman spectroscopy (SERS), also as the local excitation of quantum dots. Implementing the evanescent behavior of the plasmon field, the analysis is extended to the propagation of plasmon fields through a tandem array of metal films similar to photonic crystal structures [3, 4].

As a starting point, we describe the study of the electric field in the neighborhood of a nanoparticle using the electrostatic approximation [2]. The electric field corresponds to the plasmon particle. This model allows the description of the interaction between two plasmon particles. The interaction is extended to describe the plasmon fields propagating on a surface generating a wave behavior satisfying the Helmholtz equation where the wave number must have complex values in order to recover the traditional surface plasmon models. Controlling the random distribution of nanoparticles, we analyze the correlation effects leading us to induce localization effects. This last statement is obtained by masking thin metal surface with two independent random array hole distributions. Controlling the scale factors, we modify the curvature of the correlation trajectory. The model is related with a speckle pattern emerging from a rough surface [5]. This configuration is similar to the configuration proposed by Reather for the coupling of plasmon fields. Experimental results are shown.


2. Analysis of plasmon particle

A nanoparticle is generated by a set of atoms; the plasmon particle corresponds with the surface current distribution of the atoms. The analysis is implemented applying the electrostatic approximation given by


where ϕis the potential function. Using variable separation in Cartesian coordinates on the xyplane, the equation acquires the form


Proposing the solution as


we obtain the equation system


where the coupling constant α is a complex number having the form α=a+ib. This condition is necessary because perturbing the field, it must acquire a propagating behavior as it is shown below. Solving for X, we have


and the solution for Yis given by


Then, the complete solution ϕacquires the form


with c<0and d>0. Eq. (7) represents the boundary condition for the plasmonic field.

2.1. Description for the interaction between plasmon particles

The model is extended to describe the propagation of the electric field. For this, we propose that the electrostatic approximation is no longer fulfilled, acquiring the form of the Helmholtz equation having the form


Looking for propagation along the x-coordinate, the equation acquires the form


where kis the complex wave number k=k1+ik2. Proposing a solution of the form ϕ=XxYy, we obtain the equation system given by


whose solution acquires the form


this equation must recover the structure of the electrostatic approximation for a single nanoparticle.

From the previous solution, it is easy to identify its behavior. Along the y-coordinate the field is bounded by the exponential term, which remains unperturbed by the presence of a second particle; the interaction occurs mainly in the x-coordinate. This behavior may be generalized acquiring a wave effect. A balance relation between the complex wave number kand the constant coupling αcan be predicted; this interaction decreases the evanescent term, and the propagating term becomes dominant. This interaction is sketched in Figure 1.

Figure 1.

(a) Localized electric field for a plasmon particle. (b) Interaction between two plasmon particles. (c) Sketch to describe the generation of a plasmon field in an array of nanoparticles.

In Figure 1a, the electrostatic approximation is valid for a single nanoparticle; the wave behavior is generated by another set of particles interacting shown in Figure 1c.

Until this point we have described the generation of a wave propagating in the x-coordinate; this analysis can be extended to the propagation in the xyplane, which is analyzed in the following section.


3. Description statistics of correlation trajectories

In the present section, we describe the transfer of the statistical properties of an anisotropic two-dimensional random walk model to generate wave propagation on a metal surface, thus generating a curved surface plasmon mode. The model is conceptually simple. We describe a trajectory in a two-dimensional array, starting from a point Pwith coordinates 00. The random walk is characterized by a set of points randomly distributed, and the trajectory can be obtained from the correlation function corresponding to the flows of current probability. The statistical properties of a random distribution of points can be transferred to induce and control important physical effects. For example, it is known that the amplitude distribution of a speckle pattern follows Gaussian statistics [6, 7]. The statistic of the speckle pattern is matched with a random hole distribution, and it is transferred on a metal surface. The analysis is obtained by masking the surface metal which is considered to be formed by a set of square cells. The probability of a hole being present at the center of each cell is P; therefore, the probability of the absence of a hole is 1P. The surface contains Ncells, and the probability of the surface contains n-holes, assuming that a Bernoulli distribution is


When the number of cells Nincreases, the Bernoulli distribution tends to a Gaussian distribution of the form


where σ2is the variance. Interesting features can be identified by describing the self-correlation in this type of distribution. The simplest case occurs when two screens are superposed and, subsequently, one of them is rotated by a small angle. In order to understand the generation of the self-correlation trajectory, we focus on a single hole. In this case, it is evident that the hole follows a circular arc by joining all the points of constant probability and the complete correlation trajectory is a circle. The result in this case is shown in Figure 2a. The correlation trajectory can be controlled by inducing a scale factor in the distribution of random points. By superposing the two screens again, it is evident that the scale factor shifts the point along a linear trajectory perpendicular to the regions of constant probability, which are sets of circles, as deduced from the argument of the Gaussian distribution. The analysis is presented in an equivalent way for a speckle pattern using the fact that both of them have the same probability distribution. In Figure 2b, we show these correlation trajectories. Finally, by introducing a small rotation, the linear trajectories are curved, as shown in Figure 2c.

This result can be explained as follows: the correlation function of two scaled and rotated surfaces have the form


Analyzing the argument of the exponential function as a quadratic form, it can be shown that the curves of constant correlation are ellipses, presenting a reference system where they acquire the canonical form


The probability flows through the orthogonal trajectories between the two regions of constant probability, whose differential equation is given by


Further, the corresponding solution is given by


where c is an arbitrary constant and α=b2a2, which carries the information about the scale between the two probabilistic processes.

3.1. Graphical description and experimental implementation of the correlation trajectory

A fundamental part of the chapter consists of describing a method to generate surface plasmon fields propagating along predetermined trajectories. This can be obtained analyzing the correlation function between two screens where each one has a random hole distribution following a predetermined probability density function. This method has the characteristic that the correlation trajectory geometry presents a tunable curvature which allows the possibility to generate long-range surface plasmon.

An alternative model to generate the curved correlation trajectories is performed using a speckle pattern as it is shown in Figure 4.

The optical system that rotates the image can be a prism-type Dove. Modifying the illumination configuration using a convergent beam and changing the relative distance between the two speckle patterns obtained by shifting one mirror a scale factor are introduced. The irradiance superposition between the two speckle patterns generates the desired correlation trajectories. The speckle pattern is shown in Figure 3.

Figure 3.

Speckle pattern generated with a rough surface illuminated with a plane wave.

Figure 2.

(a) Set of points following a Gaussian distribution. (b) Correlation function between two Gaussian sets of points where one mask was rotated by a small angle. (c) Probability flow trajectories between two mask Gaussian points, one of them is scaled by approximately 95%, without rotation. (d) Same as in (c) but with a rotation of approximately5°.

Figure 4.

Experimental setup to generate speckle correlation trajectories.

Figure 5.

Masked metal surface: The typical wavelength is IR.

Figure 6.

Tandem array to propagate the plasmon field: the width of the metal is2040nmand that dielectric film is2040nm.

It is known that the irradiance function for the speckle pattern has associated a probability density function-type exponential decreasing function. The decreasing term can be matched with the decaying ratio of the plasmon mode. This configuration allows improving the generation of plasmon field avoiding the masking of the metal surface which must be made with lithography techniques. These comments represent novel applications of the speckle pattern.

The correlation trajectories generated will be implemented in the following section to describe the surface plasmon. By the fact that the correlation occurs in a curved trajectory, we expect the surface plasmon to present a magnetic behavior.

3.2. Generation of curved surface plasmon modes

The previous statistical description will be employed for the synthesis of surface plasmonic modes. The expression for the electric field of an elementary surface plasmonic mode propagating along the z-axis is given by


where β=wcε1ε2ε1+ε21/2=ξ+is the dispersion relation function and ε1,ε2represent the permittivity of the dielectric and metal, respectively. Rotating the reference system along the x-axis, the elementary surface plasmon mode acquires the form


Using the functional relation given by Eq. (17), the expression for the curved plasmonic mode is given by


By means of the Maxwell equations, we can obtain the expression for the magnetic field and the energy flux given by the Poynting vector.

For the experimental setup, we propose to illuminate a thin flat Au film (thickness 20–40 nm) with a correlated speckle pattern as shown in Figure 4. The illumination consists in two speckle patterns: each one is visualized as a set of circular motes randomly distributed following a Gaussian probability density function. The wavelength is λ=1550nm. The geometrical parameters are agreeing with those reported in [8]. The correlation curve corresponds to the surface plasmonic mode given by Eq. (20). Notably, the statistical properties of the speckle pattern are transferred to the metal surface as the plasmonic mode propagating along the correlation trajectory. In order to allow the generation of a long-range curved plasmonic mode, the correlation length must be less than 2μto guarantee resonance effects [9, 10]; this can be controlled with the roughness parameters of the surface implemented to generate the speckle pattern avoiding the power decay along the correlation trajectory. The experimental setup is sketched in Figure 5.

The analysis presented can be extended to other plasmonic configurations which are presented in the following section.


4. Propagation in a tandem array of thin metal films

The natural extension of the analysis presented is the transfer of the plasmonic mode to a tandem array of thin metal surface, shown in Figure 6. This is possible using the evanescent behavior along the x-axis of the curved surface plasmon field. This behavior has been implemented to generate an optical field redistribution propagating along an optical waveguide array [11]. In this model, the evanescent character is used to tunnel the optical field.

The transmission of the plasmonic mode satisfies the following system of differential equations:


where βis the dispersion relation function and Cirepresents the coupling constant, which depends on the relative separation between neighborhood surfaces [12]. The solution of the previous equation is similar to that presented in [11]; however, to associate a physical meaning to the coupling constant Ci, we present the analysis of two thin metal films.

The simplest case occurs when the system is formed by two thin metal films separated by a dielectric medium whose thickness must be less than 50 nm. The evanescent decay depends on the modulus of the permittivity quotient [13], and at this thickness is possible to generate tunneling effects [11]. Subsequently, the system of Eq. (21a) acquires the simple form


Rewriting it in matrix form, we obtain


It can be deduced that, as a consequence of the energy conservation, the matrix structure must be symmetric. This indicates that c1=c2=c, and the general solution is


wheredirepresents arbitrary constants and ξ1,2and η1,2represent the eigenvectors with eigenvalues λ1,2satisfying the characteristic equation depending on the coupling constant:


Moreover, it is known that the eigenvectors must be complex [14]. Subsequently, without loss of generality, the solution can be rewritten as


which indicates that the shift generated between each plasmon mode presents similar features as the coupling mode theory [12]. This analysis leads to the expression for the plasmonic mode as


where ξis a unit vector tangent to the correlation curve and sis the arc length on the same curve; we remark that the correlation trajectory is given by Eq. (20).

Eq. (24) describes the evanescent coupling through a tandem array of thin metal films. Notably, the boundary conditions of the electric field indicate that the geometry of the plasmon field generated in the first thin metal film must be preserved in all the surfaces. This shows that the transmission of the curved plasmonic mode allows inducing magnetic properties in the system [15, 16, 17, 18].


5. Conclusions

The statistical properties of the distribution of random holes or equivalently the speckle pattern were transferred to a metal surface to stablish the conditions to generate long-range curved plasmonic modes. In the case of hole distribution, this can be implemented by masking a thin metal film with two screens that allows controlling the correlation trajectory whose geometry corresponds to a curved long-range surface plasmonic mode. Another possibility was illuminating the metal thin film with two correlated speckle patterns. An important consequence of these configurations is that the set of curved surface plasmonic modes presents a vortex structure that allows to induce magnetic properties [17]. Using the evanescent character of the plasmon modes, the electric field was transferred to the propagation in a tandem array of thin metal films offering applications to design photonic crystals with tunable and localized magnetic properties.

The theoretical point of view presented in this study allows incorporating other effects such as percolation effects which consist in propagating the electric field through random structures. The main characteristic is that the plasmon field presents fractal properties which are the origin of interesting magnetic properties implicit in the curved trajectory of the set of plasmonic modes; more details can be found in [18]. The model presented can be extended by implementing different hole distribution geometries which modify the plasmonic resonance effects. Notably, the curved trajectories have associated focusing regions, and, subsequently, the corresponding magnetic singularity offers the possibility of implementation in the generation of plasmonic magnetic mirrors.

Finally, we remark that the analysis presented offers applications to photonic crystal as a metamaterial design [19, 20, 21, 22, 23] since breaking the periodicity or incorporating another type of metal on a selected region is similar to doping the structure and then is possible to induce localization effects. The excitation of plasmon fields using a speckle patterns offers the possibility to incorporate the tunable behavior of the correlation trajectory offering interesting applications in the development of plasmonic antennas and synthesis of accelerating plasmon modes [21], extending the plasmonic optical models.



The authors MATR and MVM are very grateful to CONACyT for their support.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Gabriel Martinez Niconoff, Marco Antonio Torres Rodriguez, Mayra Vargas Morales and Patricia Martinez Vara (November 12th 2018). Synthesis of Curved Surface Plasmon Fields through Thin Metal Films in a Tandem Array, Quantum Electronics, Faustino Wahaia, IntechOpen, DOI: 10.5772/intechopen.81931. Available from:

chapter statistics

626total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Localized Excitation of Single Atom to a Rydberg State with Structured Laser Beam for Quantum Information

By Leila Mashhadi and Gholamreza Shayeganrad

Related Book

First chapter

Mie-Scattering Ellipsometry

By Yasuaki Hayashi and Akio Sanpei

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us