Open access peer-reviewed chapter - ONLINE FIRST

# Dyadic Green’s Function for Multilayered Planar, Cylindrical, and Spherical Structures with Impedance Boundary Condition

By Shiva Hayati Raad and Zahra Atlasbaf

Submitted: December 4th 2020Reviewed: January 5th 2021Published: December 2nd 2021

DOI: 10.5772/intechopen.95834

## Abstract

The integral equation (IE) method is one of the efficient approaches for solving electromagnetic problems, where dyadic Green’s function (DGF) plays an important role as the Kernel of the integrals. In general, a layered medium with planar, cylindrical, or spherical geometry can be used to model different biomedical media such as human skin, body, or head. Therefore, in this chapter, different approaches for the derivation of Green’s function for these structures will be introduced. Due to the recent great interest in two-dimensional (2D) materials, the chapter will also discuss the generalization of the technique to the same structures with interfaces made of isotropic and anisotropic surface impedances. To this end, general formulas for the dyadic Green’s function of the aforementioned structures are extracted based on the scattering superposition method by considering field and source points in the arbitrary locations. Apparently, by setting the surface conductivity of the interfaces equal to zero, the formulations will turn into the associated problem with dielectric boundaries. This section will also aid in the design of various biomedical devices such as sensors, cloaks, and spectrometers, with improved functionality. Finally, the Purcell factor of a dipole emitter in the presence of the layered structures will be discussed as another biomedical application of the formulation.

### Keywords

• biomedical
• integral equation
• 2D materials
• sensors
• Purcell factor

## 1. Introduction

Planarly, cylindrically and spherically layered media have been widely used to model the human skin, body, or head. In particular, a rectangular slab is proposed to model and analyze skin temperature distribution [1]. Moreover, in the multilayered skin model, three stacked layers are exploited to simulate the performance of the epidermis, dermis, and sub-cutis parts of the skin [2]. In other research, a planarly layered medium has been proposed as the simplified human body model by considering the impact of skin, fat, and muscle in the electromagnetic performance [3]. Cylindrical-shaped equivalent phantom of the skin is alternatively used to characterize the interactions between an antenna and the human body [4]. For a more precise investigation, multilayered cylinders are proposed to model a biological system with different tissues [5]. Considering spherical geometrics, the interactions of a three-layered spherical human head model with a finite-length dipole and a cellular phone helical antenna are investigated using Green’s function [6, 7]. The research is extended to the six-layered model, including skin, fat, bone, dura, CSF, and brain layers of the head [8]. This model has also been used to explore the impact of the shapes and positions of coils in the MRI system [9]. For the analysis of the aforementioned structures, integral equation (IE) based methods have been widely used. Dyadic Green’s function, which is a powerful method to calculate the electromagnetic response for different excitation sources, plays an essential role in the IE method [10, 11, 12].

Different mathematical approaches can be utilized to calculate the dyadic Green’s function. The dyadic Green’s function in a homogeneous environment can be expressed in terms of vector functions M, N, and L. The Green’s function is singular when the field and observation points coincide. In this case, Green’s function can be written as the main part plus the portion proportional to the impulse function. This type of decomposition of the Green’s function is not unique and depends on the shape of the volume that separates the environment [13]. In general, to obtain the dyadic Green’s function in a layered medium, the Green’s function of the homogeneous environment can be used. Then, the effect of inhomogeneity can be considered by adding reflected and transmitted waves to the dyadic Green’s function [14]. In a similar approach, called the scattering superposition method, the scattering Green’s function in each layer is expressed in terms of vector wave functions with the unknown coefficients that are obtained by applying the boundary conditions [15]. Moreover, to calculate the Green’s function of a medium by impedance method, instead of each dielectric layer a transmission line, and instead of each metallic layer a current source can be considered in the equivalent model in the rectangular, cylindrical, and spherical coordinates [15, 16, 17].

Considering medical diagnostics and treatment, planarly, cylindrically and spherically layered media are engineered mainly with plasmonic materials. In this regard, the interaction of an environment with a nanometer-scale dipole emitter is of interest in different biomedical fields. For example, the optical activity of the proteins can be investigated using Green’s tensor approach. Furthermore, a single excited molecule in the vicinity of a metallic structure can be potentially used in the sensors because of behaving as a resonant filter. Moreover, the point source is an appropriate model for the concentrated light sources which are used in medical applications [18]. Following these trends, different plasmonic structures have been presented theoretically and realized in the real environment. For example, the effect of the size of the gold nanoparticles on the decay rate and the energy transfer of dipole emitters is investigated using Green’s tensor formulation and compared with the results obtained from the fabrication [19]. The same analysis is carried out for a metallic cylinder coated with a dielectric layer and the resulted integrals are solved numerically [20]. Also, Fermat’s golden rule is used to connect the imaginary part of the Green’s function to the radiation impedance of the dipole antenna adjacent to a medium, and the Purcell factor is extracted for the planar meta-surface, which is then experimentally characterized in the microwave frequency band [21].

Recently, graphene’s plasmons are proposed as the low loss and reconfigurable alternative to the plasmons of the noble metals. There are two approaches to use graphene in analytical and numerical methods. In the first one, the graphene boundary is modeled with a 2D surface, characterized by its surface conductivity, whereas in the second one, the graphene layer is replaced by a very thin dielectric [22]. Although the latter can be analyzed using the available formulas previously presented for the multilayered dielectric structures in the literature [13, 23], graphene analysis using surface conductivity model has many advantages. First, in the dielectric model, each graphene interface adds an extra layer to the structure. Therefore, when the number of graphene layers is large, the problem becomes very complicated. Secondly, in the dielectric model, it is necessary to compute the special functions of the cylindrical and spherical coordinates with complex arguments, which requires the implementation of specific algorithms for their effective calculation [24], while in the impedance boundary condition method, the surface conductivity of graphene appears as a coefficient for special functions. Third, in the dielectric model, the thickness of the graphene layer is considered to be about 0.335 nm, which is often very small compared to other geometrical parameters and makes the convergence of the analytical functions slow [25]. Also, if the goal is using the dielectric model in numerical methods, it is necessary to use a dense mesh for the equivalent dielectric of the graphene that is not optimal in terms of time and memory [26].

Applying the graphene surface conductivity model for the derivation of the Green’s function has been considered in recent years. For example, for the graphene sheet under electric bias, the dyadic Green’s function is derived by the Hertzian potential and plane-wave expansion methods and the corresponding integrals are solved with the saddle point method [27, 28]. Also, the method is expanded for the analysis of graphene with tensor surface conductivity which can be used to analyze graphene with magnetic bias or spatial dispersion. Romberg’s integration procedure is proposed for the numerical solution of the resulting integrals [29]. In another research, the analysis of the electric dipole in the proximity of the parallel plate waveguide with graphene walls is studied by extracting the Green’s function and calculating the spontaneous emission. It is found that symmetric and asymmetric plasmonic modes lead to a sharp increase in this parameter [30]. As another instance, a point source is taken into account in the vicinity of the infinite cylinder with graphene cover and it is observed that by changing the distance of the source from the cylinder as well as changing the chemical potential of the graphene layer, the Fano resonances can be controlled [31]. In this chapter, the dyadic Green’s function of various planar, cylindrical, and spherical geometries with impedance boundary conditions will be calculated using the scattering superposition method. Specifically, we have focused on the graphene material due to its wide range of applications. Apparently, another 2D material can be considered by replacing the graphene surface conductivity with the surface conductivity of the desired material. The presented formulas can be potentially sued to design various biomedical devices. Moreover, by approaching the surface impedance to zero, these structures can be potentially used to investigate the interaction of the human body with different electromagnetic sources. For the Green’s function calculation of complex media in different coordinates using vector wave functions, the reader is referred to [32, 33, 34, 35].

## 2. Surface conductivity of graphene material under different conditions

Graphene is a two-dimensional material made of carbon atoms and can be considered in the solution of Maxwell’s equations using surface conductivity boundary condition [23]. Depending on the geometrical and optical conditions, graphene surface conductivity can be isotropic or anisotropic. The purpose of this section is to provide an overview of the graphene surface conductivity under different conditions (electric bias, magnetic bias, and spatial dispersion) and for different geometries (continuous or patterned sheets).

MethodQualityArea
Mechanical exfoliation of highly ordered pyrolytic graphite (HOPG)Very highSmall
Epitaxial growth of graphene on silicon carbide (SiC)MediumLarge (3–4 inches wafers)
Reduction of graphene oxide (rGO)MediumLarge
Chemical vapor deposition (CVD)HighVery large (30 inches)

### Table 1.

Comparison of the quality and area of synthesized graphene using different techniques [36].

It should be noted that graphene material is mainly synthesized through four methods, including 1) mechanical exfoliation of highly ordered pyrolytic graphite (HOPG), 2) the epitaxial growth of graphene on silicon carbide (SiC), 3) the reduction of graphene oxide, and 4) chemical vapor deposition (CVD) technique. The comparison of these methods in terms of quality and the fabricated area is provided in Table 1 [36].

In this regard, a monolayer graphene film with metallic electrodes is transferred to a high impedance surface (HIS) which is realized by the printed circuit board (PCB) technology at microwave frequencies. The measurement is conducted in a small microwave chamber [37]. In another research in the same spectrum, an infrared laser is used to etch the CVD grown graphene sheet with predefined periodicity and later investigate the absorption of the designed structure inside a rectangular waveguide [38]. The electron beam lithography is another approach used for pattering the graphene sheet for enhanced light matter interaction [39]. Moreover, a transparent graphene millimeter wave absorber constructed by multiple transfer-etch processing is characterized by reflectometery technique at 140 GHz [40]. Also, CVD-grown graphene is used to enhance the sensitivity of the surface-enhanced Raman spectroscopy (SERS)-based chemical sensor. The measurement is done using a Raman spectrometer at the laser wavelength of 785 nm (red) [41]. In another sensor chip, DNA is hybridized to the graphene-based substrate under UV light with the wavelength of 260 nm. In this sensor, the atomic force microscopy (AFM) is used to ensure the continuity and uniformity of the synthesized graphene, and Raman characterization is used to investigate its quality and the number of layers [42].

### 2.1 Graphene material under electric bias

When a graphene sheet is under electric bias, its surface conductivity is isotropic and can be approximately calculated by using Kubo’s formulas as [43]:

σintra=2ie2kBT2πω+i/τln2coshμc2kBTE1
σinter=e24×12+1πarctanω2μc2kBTi2πlnω+2μc2ω2μc2+2kBT2E2

In the above equations, Tis the temperature, μcis the chemical potential of graphene, is the reduced Planck’s constant, KBis the Boltzmann’s constant, and ω is the angular frequency. At low-THz frequencies, the inter-band contribution of the surface conductivity can be neglected. Also, the graphene layer can be modeled as a dielectric with very low thickness δand equivalent dielectric constant ε=1iσωε0δ[22]. Note that in the above equations, graphene is assumed to be in the linear region, otherwise, other terms proportional to 1/ω3and 1/ω4should be added respectively for the frequency range of ω<2μcand ω2μc[44].

To increase the light-matter interaction of graphene material, the nano-pattering method has been proposed [45]. The surface conductivity of the periodic graphene elements under the electric bias is also isotropic. For the square patches shown in Figure 1(a), by considering the periodicity of D0 and the air gap distance of g, closed-form surface impedance (inverse of the surface conductivity) is as [47]:

Zs=D0σsD0g+iπ2ωε0ε2+12D0lncscπg2D0E3

The effective surface conductivity of periodic graphene elements with arbitrary pattern can be measured or extracted using the parameter retrieval method through full-wave simulation [48]. It should be noted that for computing the electric field required for each considered chemical potential, approximate equations can be derived as [49]:

μceV=λ1E0λ2,E0V/nm0λ1E0λ2,E0V/nm<0E4

where, λ1=0.3677and λ2=0.5010. For the chemical potentials in the range of [−1,1] eV, the required bias fields are in the order of several volts per nanometer, which can be implemented practically [29].

The surface conductivity of densely packed graphene strips, as illustrated in Figure 1(b), is anisotropic and can be approximated in the form of a diagonal tensor using the effective medium formulation as [50]:

σxx=WσσCLσC+σyy=σWLσxy=σyx=0E5

In the above equations, Wand Lare the width and periodicity of the strips, respectively. Also, σis the surface conductivity of graphene under electric bias, and σcis the static conductivity of the surface. Two very important properties of this environment are the existence of near-zero surface conductivity and hyperbolic dispersion region with the potential applications in 2D lens structures and the spontaneous emission enhancement of the dipole emitters [50, 51].

### 2.2 Graphene sheet with spatial dispersion effects

When the graphene sheet is placed on a substrate with a high dielectric constant, its surface conductivity is a tensor in which the elements depend on the wave propagation constant in the structure. This is called the spatial dispersion effect and the associated formulae for calculating surface conductivity are [29]:

σxx=σ+αd2dx2+βd2dy2E6
σyy=σ+βd2dx2+αd2dy2E7
σxy=σyx=2βd2dxdyE8

The parameters σ, α, and βare extracted for an unbiased sheet μc=0thorough perturbation theory [29]. The resulting equations are valid for the electrically biased sheet μc0, as well [52].

### 2.3 Graphene sheet under magnetic bias

When a graphene sheet is under magnetic bias, its surface conductivity is also a tensor. The diagonal elements of this tensor are equal and the off-diagonal elements are opposite in sign defined as [53].

σxxωB0=σ01iωτωcτ2+1iωτ2E9
σyxωB0=σ0ωcτωcτ2+1iωτ2E10

where σ0=2e2τπ2kBTln2coshμc2kBTand ωc=eB0νF2μc. The approximate formulas for the calculation of the surface impedance of the square graphene elements under magnetic bias, shown in Figure 1(a), are as follows [54]:

Z¯¯p=FGZ¯¯g+i2αμ0ε0εeff1001E11

where FG=0.6D0/D0g3+0.4and α=k0D0εeffπlnsinπg2D0. Similar to the electrically biased patterned elements with arbitrary shapes, the parameter retrieval method can be used under the applied magnetic bias [55]. Given the discussion of the above three sections, it is observed that assuming the surface conductivity of graphene as:

σ¯¯=σxxσxyσyxσyyE12

All items expressed above can be extracted as a special case.

## 3. Analysis of graphene-based structures using dyadic Green’s function

Dyadic Green’s functions for the planarly, cylindrically, and spherically layered structures with graphene interfaces will be derived in this section. To this end, the boundary conditions of the continuity and discontinuity of tangential electric and magnetic fields are respectively satisfied regarding the considered surface conductivity model for the graphene.

### 3.1 Graphene-based planarly layered media

This section aims to obtain dyadic Green’s functions for planar structures with graphene boundaries. This problem can be solved either in the rectangular or cylindrical coordinates. In the first sub-section, a graphene sheet with the tensor surface conductivity boundary condition (TSCBC) is considered and its dyadic Green’s function is calculated in the rectangular coordinates. In this case, the anisotropy of the surface impedance causes the coupling of the transverse electric (TE) and transverse magnetic (TM) fields. As a result of coupling, the number of unknown coefficients in the expansion of the dyadic Green’s function is increased concerning the electrically biased sheet. In the second part of this section, a graphene-dielectric stack with an arbitrary number of layers is investigated considering the electric bias for the graphene sheets. This problem is solved in the cylindrical coordinates to simplify the calculation of the resulted Sommerfeld integrals.

#### 3.1.1 Graphene sheet with the tensor surface conductivity boundary condition

The purpose of this section is to obtain the dyadic Green’s function of a graphene sheet with the tensor surface conductivity in the interface of half-spaces, as shown in Figure 2(a). The constitutive parameters of the top and bottom regions are considered as (ε1, μ1) and (ε2, μ2), respectively. Without losing the generality of the problem, the source is assumed to be in the first environment and the graphene boundary is considered in z = 0 interface. Dyadic Green’s function of this structure will be calculated using the scattering superposition method. For this purpose, the Green’s function in each region of the problem is written in the form of the Green’s function in the absence and presence of structure. Thus [23]:

G¯¯e11=G¯¯e0R¯R¯+G¯¯es11R¯R¯z>0E13
G¯¯e21=G¯¯es21R¯R¯z<0E14

where G¯¯es11and G¯¯es21are respectively the scattering Green’s function in regions 1 and 2 and they are expanded in terms of M¯and N¯vector wave functions. Also, G¯¯e0is the free-space Green’s function which can be computed using the Gmmethod as:

G¯¯e0R¯R¯=1k12ẑẑδR¯R¯+dkxdkyC1×M¯h1M¯h1+N¯h1N¯h1z<zE15

where C1=i8π2h1kx2+ky2. Using the scalar wave function ψk¯=expikxx+ikyy+ikzz, it can be readily found that:

M¯k¯=ikyx̂kxŷψk¯E16
N¯k¯=kzkjkxx̂+kyŷ+kx2+ky2kjẑψk¯E17

The parameters kx, ky, and kzare the wavenumbers in x, y, and zdirections, respectively, and kjshows the wavenumber for j = 1, 2. These parameters are not independent and are related to each other via kz=±kj2kx2ky2=±kj2kρ2. The vector wave function M¯represents the electric field of TE modes and the vector wave function N¯shows the electric field of the TM waves. In the structure under consideration, the anisotropy of surface conductivity leads to the coupling of TE and TM fields. Therefore [56]:

G¯¯es11R¯R¯=dkxdkyC1×a1M¯h1+a1N¯h1M¯h1+b1N¯h1+b1M¯h1N¯h1E18
G¯¯es21R¯R¯=dkxdkyC1×a2M¯h2+a2N¯h2M¯h1+b2N¯h2+b2M¯h2N¯h1E19

The unknown coefficients a1, a1,b1,b1,a2, a2, b2, and b2will be obtained by applying the boundary conditions. Using the self and mutual orthogonality of the vector wave functions, the above equations can be divided into two systems of equations, each with four unknown coefficients. The boundary conditions on the electric and magnetic Green’s functions respectively state that:

ẑ×G¯¯e11G¯¯e21=0E20
ẑ××G¯¯e11iωμ1×G¯¯e21iωμ2=σ¯¯.G¯¯e21E21

After applying the above boundary conditions and removing the coupling effect from the tangential components of the electric field, and by defining, A=σxxkx+σxyky, B=σyxkx+σyyky, C=σxxkyσxykx, and, D=σyykxσyxky, the unknown coefficients of the TE waves, a1=ΔTE1ΔTEand a1=ΔTE2ΔTE, can be obtained as [46]:

ΔTE2=ω2h1μ1DkyCkxE24

where, P+=h1μ1+h2μ2, Q+=k1μ1+k22h1k1h2μ2, P=h1μ1h2μ2, and Q=k1μ1k22h1k1h2μ2. By defining the TM waves expansion coefficients as b1=ΔTM1ΔTMand b1=ΔTM2ΔTM, it can be shown:

ΔTM2=2ωh1μ1AkyBkxE27

Other unknown coefficients can be obtained using decoupling equations in [46]. To validate the obtained coefficients, the structure of Figure 1(a) consisting of square patches with D0 = 5 μm, g = 0.5 μm, μc = 0.5 eV and τ = 0.5 ps with plane wave illumination is considered under electric and magnetic biases. Since Green’s function coefficients are the same as reflection and transmission coefficients of the plane wave, the results of Green’s function are compared with the results of the circuit model as [54]:

S31=22+η0σd2+η0σd2+η0σ02E28
S41=2η0σ02+η0σd2+η0σ02E29

Figure 3 shows the magnitude of the transmission coefficient by considering the electric bias for the graphene layer. The results of the two methods are identical, and because of the absence of electromagnetic coupling under electric bias, the transmission coefficient due to mutual coupling is zero. In Figure 4 the same results are illustrated for the applied magnetic bias of 0.5 Tesla. There is good agreement in the magnitude and phase of the transmission coefficient in the abovementioned two methods. Also, by finding the poles of the coefficients, the electromagnetic wave propagation constants for the electrically and magnetically biased graphene sheets can be obtained which are in full agreement with [28, 57], respectively. The correctness of the extracted coefficients confirms the validity of the dyadic Green’s function formulation.

#### 3.1.2 Graphene-dielectric stack

Dyadic Green’s function for an N-layer dielectric environment has been previously formulated using the scattering superposition method [58]. In this section, the above equations are extended to the environment with the electrically biased graphene boundaries, as shown in Figure 2(b). The graphene boundary can be either continuous or periodically patterned as discussed in section 2. To start the analysis, the layers are numbered by starting from the top layer, and an arbitrary field point iand source point jare assumed. The problem is solved in the cylindrical coordinates. Since the cylindrical wave functions are discussed in detail in the next section, they are not mentioned here. The dyadic Green’s function can be expanded as [58]:

G¯¯esijrr=i4πdhn=02δn0λhj1δiNM¯hi1δj1AMijM¯hj+1δjNBMijM¯hj+1δiNN¯hi1δj1ANijN¯hj+1δjNBNijN¯hj+1δi1M¯hi1δj1CMijM¯hj+1δjNDMijM¯hj+1δi1N¯hi×1δj1CNijN¯hj+1δjNDNijN¯hjE30

In the above equations, M¯and N¯are vector wave functions in the cylindrical coordinate system, and AM,Nij, BM,Nij, CM,Nij, DM,Nijare unknown coefficients. The Kranoker delta function is used in the field expansion to generalize the formulation for the arbitrary locations of the field and source points. By applying the boundary conditions on the tangential components of the electric field at the interface of the arbitrary two layers, denoted by iand i + 1, it can be shown that:

AMijBMijOiANijOiBNijeihizi+CMijDMij+δjiOiCNijOiDNij+δjieihizi=AMi+1j+δi+1jBMi+1jOi+1ANi+1j+δi+1jOi+1BNi+1jeihi+1zi+CMi+1jDMi+1jOi+1CNi+1jOi+1DNi+1jeihi+1ziE31

where Oi=hiki,Pi=hiμi,Qi=kiμi. The boundary condition on the tangential components of the magnetic field yields:

Pi+1AMi+1j+δi+1jPi+1BMi+1jQi+1ANi+1j+δi+1jQi+1BNi+1jeihi+1ziPi+1CMi+1jPi+1DMi+1jQi+1CNi+1jQi+1DNi+1jeihi+1ziPiAMijPiBMijQiANijQiBNijeihizi+PiCMijPiDMij+δjiQiCNijQiDNij+δjieihizi=iωσi+1iAMijBMijPiANijPiBNijeihizi+iωσi+1iCMijDMij+δjiPiCNijPiDNij+δjieihiziE32

By re-writing the coefficients as a matrix:

AM,Ni+1j+δi+1jBM,Ni+1jCM,Ni+1jDM,Ni+1j=1TFiH,VRFiH,VTFiH,VRPiH,VTPiH,V1TPiH,VAM,NijBM,NijCM,NijDM,Nij+δijE33

the outgoing and incoming reflection and transmission coefficients can be defined and used to extract the recursive relations. The coefficients for M¯sources are:

RFiH=μihi+1μi+1hi+gμihi+1+μi+1hi+gE34
RPiH=μihi+1μi+1higμihi+1+μi+1higE35
TFiH=2μihi+1μihi+1+μi+1hi+gE36
TPiH=2μihi+1μihi+1+μi+1higE37

The coefficients for N¯sources are:

RFiV=μihiki+12ki2μi+1hi+1+ghihi+1μihiki+12+μi+1hi+1ki2+ghihi+1E38
RPiV=μihiki+12μi+1ki2hi+1ghihi+1μihiki+12+μi+1ki2hi+1ghihi+1E39
TPiV=2μi+1kiki+1hi+1μihiki+12+μi+1ki2hi+1ghihi+1E40
TFiV=2μihi+1kiki+1μihiki+12+μi+1hi+1ki2+ghihi+1E41

The superscripts Hand Vrespectively denote TE and TM sources. Also, subscripts Fand Pare used to show the incoming and outgoing waves, respectively. The procedure of extracting the unknown coefficients using (34)-(41) is discussed in [58]. To validate the proposed formulas, a parallel plate waveguide with graphene walls is considered. To extract the characteristic equation using the proposed formulations, it is necessary to force the denominator of the coefficients equal to zero. For this three-layer medium:

T1=T2.T1=.1TF2Heih2h1dRF2HTF2Heih2h1dRP2HTP2Heih2+h1d1TP2Heih2h1d.1TF1HRF1HTF1HRP1HTP1H1TP1HE42

From which it can be deduced that BM,N11=T121T111. By setting T111equal to zero and assuming that the medium (1) and (3) are the same, and also defining h=iωμ1σ, for the Hcoefficients it can be concluded that:

eih2d1eih2d+1h1h2+h=0h2+h+jtanh2d2h1=0E43

This procedure is repeatable for Vsources. Also, to calculate the reflection coefficient from a multilayer structure, it is necessary to consider the field and source points in region 1. In this case, the only non-zero coefficient in Green’s function expansion is BM,N11coefficient representing the plane wave reflection coefficient from the multilayer structure.

### 3.2 Graphene-based cylindrical structures

In this section, the dyadic Green’s function of a cylindrical structure with the tensor surface conductivity boundary condition will be extracted. Later, different examples of guiding and scattering problems are provided to investigate the validity of the formulation. In general, in cylindrical structures, TE and TM modes are coupled, which leads to the complexity of mathematical relations. Therefore, the generalization of the formulation to the multilayered cylinders is not considered here. Note that graphene sheets can be wrapped around cylindrical particles due to the presence of van der Waals force [59]. For this purpose, tape-assist transfer under micromanipulation and spin-coating methods are proposed [60].

#### 3.2.1 Dyadic Green’s function for a cylinder with tensor surface conductivity boundary condition

The dyadic Green’s function of the single-layer cylinder with the tensor surface conductivity boundary condition, as considered in Figure 5., will be extracted in the following. The interior region of the cylinder is made of dielectric material and its cover is considered as a full tensor surface conductivity. To solve the problem, the vector wave functions are defined as [23]:

M¯μh=eimϕeihzimZmμrrr̂Zmμrrϕ̂E44
N¯μh=1kjeimϕeihzihZmμrrr̂hmrZmμrϕ̂μ2ZmμrẑE45

In the above equations, Zmis the cylindrical Bessel function in the inner layer and the cylindrical Hankel function in the outer layer, both with the orders of m. The wavenumber in the radial direction is μand the wavenumber along the length is h. The free-space Green’s function in the cylindrical coordinates is [23]:

G¯¯e0R¯R¯=1k2r̂r̂δR¯R¯+i8πdhm=1η2×M¯ηhM¯η1h+N¯ηhN¯η1hR<RE46

In the cylindrical coordinate system, TE and TM modes are coupled and the expansion of the fields are [23].

G¯¯es11R¯R¯=i8πdhm=1η2×AηM¯η1h+BηN¯η1hM¯η1h+CηN¯η1h+DηM¯η1hN¯η1hE47
G¯¯es21R¯R¯=i8πdhm=1η2×aξM¯ξh+bξN¯ξhM¯η1h+cξN¯ξh+dξM¯ξhN¯η1hE48

where Aη, Bη, Cη, and Dηare the unknown coefficients of the DGF expansion in region 1. Also, aξ, bξ, cξ, and dξare DGF expansion coefficients in the region 2. The boundary condition on the tangential components of the magnetic Green’s function is given by [61]:

r̂××G¯¯e11μ1×G¯¯e21μ2=σ¯¯ϕz.G¯¯e21E49

By applying the above-mentioned boundary condition along with the boundary condition regarding the continuity of the electric Green’s function to (47)(48), the system of equations to determine the unknown coefficients can be obtained as:

Hm1ηaa1k1hmaHm1ηaJmξaa1k2hmaJmξa01k1η2Hm1ηa01k2ξ2Jmξa1μ1hmaHm1ηak1μ1Hm1ηaa1μ2hmaJmξa+σ1k2μ2Jmξaa+σ21μ1η2Hm1ηa01μ2ξ2Jmξa+σ3σ4×AηBηaξbξ=Jmηaa01μ1hmaJmηa1μ1η2JmηaE50
1k1hmaHm1ηaHm1ηaa1k2hmaJnξaJmξaa1k1η2Hm1ηa01k2ξ2Jmξa0k1μ1Hm1ηaa1μ1hmaHm1ηak2μ2Jmξaa+σ1'1μ2hmaJmξa+σ2'01μ1η2Hm1ηaσ3'1μ2ξ2Jmξa+σ4'×CηDηcξdξ=1k1hmaJmηa1k1η2Jmηak1μ1Jmηaa0.E51

where:

σ1=iωσJmξaa,σ2=σzzk2ξ2Jmξa+iωσhmk2aJmξa
σ3=iωσϕϕJmξaa,σ4=iωσϕϕ1k2hmaJmξa+iωσϕz1k2ξ2Jmξa
σ1=iωσzz1k2ξ2Jmξa+iωσhmk2aJmξa,σ2=iωσJmξaa
σ3=iωσϕϕ1k2hmaJmξa+iωσϕz1k2ξ2Jmξa,σ4=iωσϕϕJmξaa

Nullifying the determinant of the matrices (50)(51) for m = 0, each of the above matrices can be separated as the multiplication of:

ωε2ξJ1ξaJ0ξaωε1ηH11ηaH01ηa+iσd=0E52
ξωμ2J0ξaJ1ξaηωμ1H01ηaH11ηa+iσd=0E53

For the graphene shell under magnetic bias, both matrices in (50)(51) result in the following equation for the propagation constant of hybrid TE and TM waves:

ωε2ξH11ξaH01ξaωε1ηJ1ηaJ0ηa+iσd×ηωμ1J0ηaJ1ηaξωμ2H01ξaH11ξa+iσd=σ02E54

Which is in agreement with [57]. Also, the total scattering cross-section (TSCS) of a graphene-coated cylinder with the parameters a = 50 μm, ε2 = 2.4, τ = 1 ps, μc = 0.25 eV under electric bias is calculated for both TE and TM polarizations in Figure 6, and compared with the results of the CST2017 software package. As can be seen, both methods have resulted in the same results.

As another example, the densely packed graphene strips with the parameters L = 420 nm, W = 400 nm, μc = 0.5 eV, and τ = 1 ps are considered around the dielectric cylinder as in Figure 5. It is assumed that the strips are wrapped around a hollow cylinder with the radius of a = 50 μm. Figure 7 shows the TSCS of this structure for the magnetic biases with the strength in the range of 20-40 T. As observed, by increasing the magnetic bias, the resonant frequency of the surface plasmons blue shifts. The associated planar structure behaves as a hyperbolic meta-surface [62]. Under locally flat consideration of the curvature, this structure can also be considered as a hyperbolic medium. In the cylindrical geometries, hyperbolic meta-surfaces can be obtained using graphene-dielectric stacks [63]. The advantage of this hyperbolic structure is its two-dimensional nature and reconfigurability. It is also demonstrated that covering the surface of nanotubes with the hyperbolic meta-surface increases the interaction of the light with dipole emitters [64].

Finally, as Figure 5(b) illustrates, graphene-based square patches around the cylinder are considered under magnetic bias. Geometrical and optical parameters are as follows: τ = 1 ps, g = 0.5 μm, and D = 0.5 μm and the TSCS is illustrated in Figure 8. for B0 = 0 T and B0 = 10 T. As can be seen, by changing the magnetic bias, the optical state changes from the maximum scattering to the minimum scattering. Such capability has recently been proposed by using a phase change material for switching between these two situations [65]. In this structure, the operating frequency can also be adjusted by changing the electric bias of graphene.

It is essential to note that in the cylindrical structures, the spatial domain Green’s function consists of an integral on the real axis. Due to the existence of poles in the integration path, the integration path is usually deformed into a triangular shape. Also, a common method for calculating the spatial domain Green’s function is the generalized pencil of function (GPOF) method in which the Green’s function is expanded in terms of the complex exponential functions [66]. The unknowns can be found via the algorithm provided in [67]. Also, as mentioned earlier, the dyadic Green’s function of the graphene-based multilayered cylindrical structures will not be considered here. For calculating the scattering cross-section of graphene-based multilayer cylindrical structures, a simple approach based on the transfer matrix method (TMM) is proposed that will be suitable for establishing novel optical devices [68, 69].

### 3.3 Graphene-based spherically layered medium

In this section, a multi-layered spherical structure with the graphene boundaries is considered and its Green’s function is extracted by assuming different locations for the source and observation points. The relationship between the Green’s function expansion coefficients and the modified Mie-Lorentz coefficients is exhibited to discuss how to solve the scattering problems using the Green’s function. Scattering analysis of graphene-based layered structures is of great importance in the design of novel optical devices [70]. Finally, the procedure for calculating the Purcell factor is considered as an important application. Instances of experimentally realized graphene-coated spherical particles can be found in [71, 72], where improved template method and hydrothermal method are proposed for the synthesis. Also, transmission electron microscopy (TEM) and field emission scanning electron microscopy (FE-SEM) are used for characterization.

#### 3.3.1 Dyadic Green’s function of a graphene-based spherically layered structure

Let us consider an N-layer spherical medium with the graphene boundaries as shown in Figure 9. The purpose of this section is to compute the dyadic Green’s function of this structure with the assumption of arbitrary locations for the field and source points. For this purpose, the Green’s function in each layer is expanded in terms of vector wave functions with unknown coefficients. These functions are calculated using the scalar wave function of [74]:

φmnrθϕ=znkprPnmcosθeimϕE55

where zn.represents the spherical Bessel or Hankel functions of order nand Pnm.is the associated Legendre function with degree nand order m. It can be readily shown that:

M¯mnkp=znkpreimϕimsinθPnmcosθθ̂dPnmcosθϕ̂E56
N¯mnkp=nn+1kprznkprPnmcosθeimϕr̂+1kprdrznkprdrdPnmcosθθ̂+imsinθPnmcosθϕ̂eimϕE57

The above vector functions are self and mutually orthogonal. This feature will be used to decouple the equation when computing the unknown coefficients.

As mentioned earlier, in the scattering superposition method, the dyadic Green’s function is written as the sum of the free-space and scattering Green’s functions. The free-space Green’s function is related to the source in an infinite homogeneous medium while the scattering Green’s function is due to the source in the presence of the layered medium. The expansion of the Green’s function for the spherical structure with Nconcentric layers, assuming that the source and field points are respectively located in the desired layers with the labels pand q, can be written as [23]:

Ĝ¯¯epqrr=Ĝ¯¯0epqrrδqp+Ĝ¯¯espqrrE58

The free-space Green’s function can be obtained using the residue theorem as [23]:

Ĝ¯¯0err=r̂r̂kq2δrr+ikq4πn=1m=nn2δm02n+1nn+1×nm!n+m!×M¯mn1kqM¯mnkq+N¯mn1kqN¯mnkqr>rM¯mnkqM¯1mnkq+N¯mnkqN¯1mnkqr<rE59

Moreover, scattering Green’s function in each layer can be expanded as [75]:

Ĝ¯¯espqrr=ikq4πn=1m=nn2δm02n+1nn+1nm!n+m!×1δpNM¯mn11δq1AHpqM¯mn+1δqNBHpqM¯1mn+1δpNN¯mn11δq1AVpqN¯mn+1δqNBVpqN¯1mn+1δp1M¯mn1δq1CHpqM¯mn+1δqNDHpqM¯1mn+1δp1N¯mn1δq1CVpqN¯mn+1δqNDVpqN¯1mnE60

where AH,Vpq, BH,Vpq, CH,Vpq, and DH,Vpqare the unknown coefficients of the Green’s function to be obtained. In the above formulation, the superscript (1) represents that the spherical Hankel functions of the first type are chosen to represent the spherical vector functions. For the other vector wave functions, the first-kind Bessel function should be selected. As observed, in the spherical structures the TE and TM waves are not coupled. Thus, the expansion of the fields only includes the interaction of the functions M¯and M¯′ as well as the functions N¯and N¯′. Moreover, to solve the problem for the structures with the arbitrary number of layers, the Kronecker delta function is used in the expansions of the fields. In the middle layers of the spherical structure, the electric field is a linear combination of the Bessel and Hankel functions, whereas in the outermost layer, only the Hankel functions, and in the innermost layer, only the Bessel functions will exist. So we have used 1δp1and 1δpNfunctions. To determine the delta functions related to the source terms, free-space Green’s function can be used. For q = 1 (external layer), the source functions related to the second criterion of the free-space Green’s function is used. Also, for q = N(internal layer) the source functions related to the first criterion should be used. In the middle layers, a linear combination of different source functions must be used.

It should be noted that by using the addition theorem in Legendre functions, the internal series in the Green’s function expansion can be eliminated [11]. To select the number of terms required for the convergence of the external series, the conditions of the problem must be considered. In the other words, in the structures whose electrical size is very much smaller than that of the wavelength, the term n = 1 is sufficient for the convergence [76]. Also, in the case where the distance between the source and observation points is large, the series can be truncated in the number Nt=x+3x+2, where x= k0R1. Otherwise, the convergence of the series is weak, and a large number of terms in the range of 20 k0R1 should be considered. In this case, the series acceleration techniques will be highly efficient in terms of computational efficiency [77, 78].

Boundary conditions on tangential components of electric and magnetic Green’s functions in the interface of two adjacent layers are [73]:

r̂×Ĝ¯¯epq=r̂×Ĝ¯¯ep+1qE61
1μp+1r̂××Ĝ¯¯ep+1q1μpr̂××Ĝ¯¯epq=σp+1pr̂×r̂×G¯¯epqE62

which respectively lead to the following equations:

ξnppAHpqξnppBHpqξnppAVpqξnppBVpq+ψnppCHpqψnppDHpq+δpqψnppCVpqψnppDVpq+δpq=ξnp+1pAHp+1q+δp+1qξnp+1pBHp+1qξnp+1pAVp+1q+δp+1qξnp+1pBVp+1q+ψnp+1pCHp+1qψnp+1pDHp+1qψnp+1pCVp+1qψnp+1qDVp+1qE63
kp+1μp+1ξnp+1pAHp+1p+δp+1qξnp+1pBHp+1pξnp+1pAVp+1q+δp+1qξnp+1pBVp+1q+kp+1μp+1ψnp+1pCHp+1qψnp+1pDHp+1qψnp+1pCVp+1qψnp+1pDVp+1qkpμpξnppAHpqξnppBHpqξnppAVpqξnppBVpqkpμpψnppCHpqψnppDHpq+δpqψnppCVpqψnppDVpq+δpq=σp+1pξnppAHpqξnppBHpqξnppAVpqξnppBVpq+σp+1pψnppCHpqψnppDHpq+δpqψnppCVpqψnppDVpq+δpqE64

In the above equations, ψnpq=jnkpRq, ξnpq=hn1kpRq, ψnpq=1/ρdρjnρρ=kpRqand ξnpq=1/ρdρhn1ρρ=kpRq. To obtain recursive relations for unknown coefficients of Green’s function, separating the above equations is necessary. Therefore [75]:

AH,Vp+1q+δp+1qBH,Vp+1q=1TFpH,VAH,VpqBH,Vpq+RFpH,VTFpH,VCH,VpqDH,Vpq+δpqE65
CH,Vp+1qDH,Vp+1q=RPpH,VTPpH,VAH,VpqBH,Vpq+1TPpH,VCH,VpqDH,Vpq+δpqE66

where the reflection and transmission coefficients for TE waves are computed as:

RFpH=kp+1μpψnp+1pψnppkpμp+1ψnppψnp+1p+gψnppψnp+1pkp+1μpψnp+1pξnppkpμp+1ξnppψnp+1p+gξnppψnp+1pE67
RPpH=kp+1μpξnp+1pξnppkpμp+1ξnppξnp+1p+gξnppξnp+1pkp+1μpξnp+1pψnppkpμp+1ψnppξnp+1p+gψnppξnp+1pE68
TFpH=kp+1μpψnp+1pξnp+1pψnp+1pξnp+1pkp+1μqψnp+1pξnppkpμp+1ξnppψnp+1p+gξnppψnp+1pE69
TPpH=kp+1μpξnp+1pψnp+1pψnp+1pξnp+1pkp+1μpξnp+1pψnppkpμp+1ψnppξnp+1p+gψnppξnp+1pE70

Also, the reflection and transmission coefficients for the TM waves are:

RFpV=kp+1μpψnp+1pψnppkpμp+1ψnppψnp+1p+gψnppψnp+1pkp+1μpψnp+1pξnppkpμp+1ξnppψnp+1p+gξnppψnp+1pE71
RFpV=kp+1μpψnp+1pψnppkpμp+1ψnppψnp+1p+gψnppψnp+1pkp+1μpψnp+1pξnppkpμp+1ξnppψnp+1p+gξnppψnp+1pE72
TFpV=kp+1μpψnp+1pξnp+1pξnp+1pψnp+1pkp+1μpψnp+1pξnppkpμp+1ξnppψnp+1p+gξnppψnp+1pE73
TPpV=kp+1μpξnp+1pψnp+1pψnp+1pξnp+1pkp+1μpξnp+1pψnppkpμp+1ψnppξnp+1p+gψnppξnp+1pE74

In the above formulas, the subscripts Fand Prepresent the outgoing and incoming waves, respectively. The symbols TPFpHand RPFpHexpress the transmission and reflection of TE waves (due to the presence of superscript H), whereas the expression TPFpVand RPFpVexpress the transmission and reflection of the TM waves (due to the presence of superscript V). Using the matrix form:

AH,Vp+1q+δp+1qBH,Vp+1qCH,Vp+1qDH,Vp+1q=1TFpH,VRFpH,VTFpH,VRPpH,VTPpH,V1TPpH,VAH,VpqBH,VpqCH,VpqDH,Vpq+δpqE75

Recursive relations can use to start through: AH,VNq=BH,VNq=CH,V1q=DH,V1q=0.

It should be noted that due to the difficulty of constructing multiple concentric graphene shells, one can consider each boundary as a dielectric or perfect electric conductor (PEC) in the optical design. For example, an optical absorber consisting of a metal-dielectric spherical resonator whose outermost layer is coated with graphene is shown to enhance the absorption of the resonator [79]. To design optical structures such as absorber and invisible cloaks using Green’s function formulation, it is necessary to transfer the point source to infinity to resemble a plane wave. In this case, considering the expansion of Green’s function, it can be observed that the only unknown coefficients of expansion are BH,Vpqand DH,Vpqcoefficients. Using the convolution integral it can be shown that [23]:

G¯¯0e.J+G¯¯es11.J=G¯¯es21.JE76

The source dependency is the same in all of the above Green’s functions and can be simplified from both sides of the equation. It is observed that the equation is converted into the equation resulting from the Mie analysis of the spherically layered structures. By directly starting from the Mie-Lorentz theory, the same results can be obtained [80]. Moreover, due to the sub-wavelength nature of the localized graphene plasmons, the final formulas can be simplified using the polynomial approximation of the special functions [81]. These structures can also be used as the building blocks of optical meta-materials [82].

#### 3.3.2 Purcell factor and energy transfer between donor-acceptor emitters

As mentioned earlier, one of the important applications of Green’s function is studying the interaction of dipole emitters in the vicinity of nanostructures. For this purpose, a vertical dipole in the vicinity of the graphene-based spherical structure, which was introduced in the previous section, is considered and its Purcell factor is calculated using the Green’s function. Assuming that the field point and observation points for the dipole with moment d0=d0r̂are in the same location of r=Δ>R1, θ=0, and ϕ=0, by calculating the scattered field using the convolution integral and the use of [18]:

where symbol Imrepresents the imaginary part of the complex function. The decay rate can be calculated. For this purpose, using relationships:

Pnmcosθθ0sinmθ=1m=00m0E78
dPnmcosθθ0msinm1θ=1m=10m1E79

The scattered field can be calculated. Thus [73]:

ΓtotalΓ0=132nnn+12n+1znk1Δk1Δ2BV11E80

As can be seen, the above equation is in full agreement with [83] which has extracted the decay rate for a core-shell plasmonic sphere, whereas in this case, it is necessary to use the Mie coefficients of graphene-based structure, namely, BV11. Using the derived formulas, the positions of the dipole can be considered arbitrarily. The transferred energy between the donor-acceptor pairs can be calculated straightforwardly using the Green’s function Gas [84]:

ΓωΓ0ω=dA.GrArBω.dB2dA.G0rArBω.dB2E81

where the subscript 0 refers to the free-space parameters and dA and dB are respectively the dipole moments of the acceptor and donor.

## 4. Conclusion

In conclusion, dyadic Green’s function extraction for planarly, cylindrically, and spherically layered medium based on scattering superposition method is a unified approach to deal with a wide range of electromagnetic problems in the realm of biomedicine. Specifically, the interaction of the human skin, body, and head with the electromagnetic sources with arbitrary distributions can be studied. Moreover, by engineering the constitutive parameters of the layers, a variety of novel devices for medical diagnostics and treatment can be proposed. Plasmonic metals and 2D materials are two main categories of such materials. For the sake of efficient analytical analysis, the impedance boundary condition is satisfied in the case of 2D materials to be used in the design of compact devices.

chapter PDF

## More

© 2021 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

### Cite this chapter Copy to clipboard

Shiva Hayati Raad and Zahra Atlasbaf (December 2nd 2021). Dyadic Green’s Function for Multilayered Planar, Cylindrical, and Spherical Structures with Impedance Boundary Condition [Online First], IntechOpen, DOI: 10.5772/intechopen.95834. Available from: