Open access peer-reviewed chapter

Full-Wave Spectral Analysis of Resonant Characteristics and Radiation Patterns of High Tc Superconducting Circular and Annular Ring Microstrip Antennas

By Ouarda Barkat

Submitted: March 18th 2012Reviewed: November 16th 2012Published: March 6th 2013

DOI: 10.5772/54665

Downloaded: 2425

1. Introduction

During the last decades, superconducting antenna was one of the first microwave components to be demonstrated as an application of high-temperature superconducting material [1-3]. High Tcsuperconducting microstrip antennas (HTSMA) are becoming popular and getting increased attention in both theoretical research and engineering applications due to their excellent advantages. Various patch configurations implemented on different types of substrates have been tested and investigated. In the design of microstrip antennas, anisotropic substances have been increasingly popular. Especially the effects of uniaxial type anisotropy have been investigated due to availability of this type of substances. These structures characterized by their low profile, small size, light weight, low cost and ease of fabrication, which makes them very suitable for microwave and millimeter-wave device applications [4-6]. HTSMA structures have shown significant superiority over corresponding devices fabricated with normal conductors such as gold, silver, or copper. Major property of superconductor is very low surface resistance; this property facilitates the development of microstrip antennas with better performance than conventional antennas. Compared to other patch geometries, the circular and annular ring microstrip patches printed on a uniaxial anisotropic substrate, have been more extensively studied for a long time by a number of investigators [7-11]. Circular microstrip patch can be used either as radiating antennas or as oscillators and filters in microwave integrated circuits (MIC's). The inherent advantage of an annular ring antenna is: first, the size of an annular ring for a given mode of operation is smaller than that of the circular disc resonating at the same frequency. Second, The presence of edges at the inner and outer radii causes more fringing than in the case of a circular microstrip antenna, in which fringing occurs only at the outer edge, and this implies more radiation from these edges, that leads to high radiation efficiency. These structures are quite a complicated structure to analyze mathematically. Different models are available to model a microstrip antenna as the transmission-line model and the cavity model in simple computer aided design formulas. However, the accuracy of these approximate models is limited, and only suitable for analysing simple regularly shaped antenna or thin substrates. The full-wave spectral domain technique is extensively used in microstrip antennas analysis and design. In this method, Galerkin’s method, together with Parsval’s relation in Hankel transform domain is then applied to compute the resonant frequency and bandwidth. The integral equation is formulated with Hankel transforms which gives rise to a diagonal form of the Green’s function in spectral domain.

The numerical results for the resonant frequency, bandwidth and radiation pattern of microstrip antennas with respect to anisotropy ratio of the substrate, are presented. The Influence of a uniaxial substrate on the radiation of structure has been studied. To include the effect of the superconductivity of the microstrip patch in the full wave spectral analysis, the surface complex impedance has been considered. The effect of the temperature and thickness of HTS thin film on the resonant frequency and bandwidth have been presented. Computations show that, the radiation pattern of the antenna do not vary significantly with the permittivity variation perpendicular to the optical axis. Moreover, it is found to be strongly dependent with the permittivity variation along the optical axis. The computed data are found to be in good agreement with results obtained using other methods. Also, the TM and TE waves are naturally separated in the Green's function. The stationary phase method is used for computing the far-zone radiation patterns.


2. Theory

The antenna configurations of proposed structures are shown in Figure 1. The superconducting patches are assumed to be located on grounded dielectric slabs of infinite extent, and the ground planes are assumed to be perfect electric conductors. The substrates of thickness d are considered to be a uniaxial medium with permittivity tensor:

Figure 1.

Cross section of a superconducting microstrip patch on uniaxial anisotropic media. (a) circulaire ; (b) annulaire ring

ε¯j=ε0 |εjx000εjy000εjz|E1

Where εjx= εjy≠ εjz (j=1, 2), and the permeability will be taken as μ0.

Starting from Maxwell’s equations in the Hankel transform domain, we can show that the transverse fields inside the uniaxial anisotropic region can be written as [12]:

Eρ(ρ,ϕ,zj)=n=einφ0dkρkρ  [ωμ02kρ (J˜n1(kρρ)+J˜n+1(kρρ))  H˜jz(kρ,zj)+iεjz2kρεjx(J˜n+1(kρρ)+J˜n1(kρρ)) E˜jz(kρ,zj) z]E2
Eϕ(ρ,ϕ,zj)=n=einφ0dkρkρ  [iωμ02 kρ(J˜n1(kρρ)J˜n+1(kρρ)) H˜jz(kρ,zj)εjz2εjxkρ(J˜n+1(kρρ)+J˜n1(kρρ)) E˜jz(kρ,zj) z]E3
Hφ(ρ,ϕ,zj)=n=einφ0dkρkρ  [(iωεjzε02kρ(J˜n1(kρρ) J˜n+1(kρρ))E˜jz(kρ,zj)  +12kρ(J˜n1(kρρ)+J˜n+1(kρρ)) H˜jz(kρ,zj) z]E4
Hρ(ρ,ϕ,zj)=n=einφ0dkρkρ  [ωεjzε02kρ(J˜n1(kρρ)+J˜n+1(kρρ))E˜jz(kρ,zj) +i2kρ(J˜n+1(kρρ)J˜n1(kρρ)) H˜jz(kρ,zj) z]E5

We can put these equations in the following form:

E(ρ,ϕ,zj)=[Eρ(ρ,ϕ,zj)Eϕ(ρ,ϕ,zj)] =n=+ejnϕ  0+dkρ.  kρ.H¯n(kρρ). [(iεjz/εjx). E˜jz(kρ,zj)/ ziωμ0H˜jz(kρ,zj)]E6
H(ρ,ϕ,zj)=[Hϕ(ρ,ϕ,zj)-Hρ(ρ,ϕ,zj)]  =n=+ ejnϕ 0+dkρ.kρ.H¯n(kρρ).  [iωεjzε0E˜jz(kρ,zj)i H˜jz(kρ,zj)/ z]E7

That is


The kernel of the vector Hankel transform is given by [12]:

H¯n(kρρ)=|J˙n(kρρ)in Jn(kρρ)/kρin Jn(kρρ)/kρρJ˙n(kρρ)|E10


       e˜n(kρ,z)=A¯nj (kρ) e -i k¯jzz+B¯nj (kρ) ei k¯jzzE11
    h˜n(kρ,z)=g¯j (kρ.(A¯nj (kρ) e -i k¯jzzB¯nj (kρ) ei k¯jzzE12

A¯and B¯are two unknown vectors and g¯(kρ)is determined by:



k0: is the free space wavenumber,

kjze=(εjxk02(εjxkρ2/εjz))1/2: is TM propagation constants in the uniaxial substrate.

kjzh=(εjxk02kρ2)12: is TE propagation constants in the uniaxial substrate.

In the spectral domain, the relationship between the patch current and the electric field on the microstrip is given by [10, 11]:


K˜(kρ)is the current on the microstrip which related to the vector Hankel transform of K(ρ). The unknown currents are expanded, in terms of a complete orthogonal set of basis functions, issued from the magnetic wall cavity model. It is possible to find a complete set of vector basis functions to approximate the current distribution, by noting that the superposition of the currents due to TM and TE modes of a magnetic-wall cavity form a complete set. The current distribution of the nth mode of microstrip patch can be written as [11, 12]:

Kn(ρ)=p=1PanpΨnp(ρ)+q=1QbnqΦnq(ρ)             E15

Here anp and bnq are unknown coefficients.

For superconducting annular ring patch, (αnq,  βnp)   are the roots of dual equations φn(αnqb/a)=0and ψn'(βnpb/a)=0.

The Hankel transforms of ψnp and φnq functions are described as [12]:

Ψ˜np(kρ)={[(βnp/a)((βnp/a)2k2ρ)Ynp(kρ)na(βnpkρ)Ynp(kρ)]a<ρ<b0ρ>b,    ρ<a            E16
Φ˜nq(kρ)={[0kρa(k2ρ(αnq/a)2)Znq(kρ)]a<ρ<b                    0                                                            ρ>b,    ρ<aE17


Ynp(kρ)=ψn(βnpb/a) Jn(kρb)ψn(βnp) Jn(kρa)E18
Znq(kρ)=(b/a) φn(αnq b/a) Jn(kρb)φn(αnq) Jn(kρa)E19

Using the same procedure, the basis functions for superconducting circular patch, are given by the expressions [11]:

Ψ˜np(kρ)=βnpaJn(βnpa)[J˙n(kρa)βnp2kρ2inkρβnp2aJn(kρa)]T     E20
φ˜nq(kρ)=[0kρaJ˙n(αnqa)Jn(kρa)kρ2αnq2]T            E21


J˙n(βnpa)=0and Jn(αnqa)=0 .  

Jn(.) and Nn (.) are the Bessel functions of the first, and second kind of order n.

G¯(kρ)is the spectral dyadic Green’s function, and after some simple algebraic manipulation, we determine the closed form of the spectral Green dyadic at z=d for a grounded uniaxial substrate.




In order to incorporate the finite thickness, the dyadic Green’s function is modified by considering a complex boundary condition. The surface impedance of a high-temperature superconductors (HTSs) material for a plane electromagnetic wave incident normally to its surface is defined as the ratio of |E| to |H| on the surface of the sample [13]. It is described by the equation:

Zs=Rs+i XsE25

Where Rsand Xsare the surface resistance and the surface reactance.

If the thickness t of the strip of finite conductivity σ is greater than three or four penetration depths, the surface impedance is adequately represented by the real part of the wave impedance [13].


If t is less than three penetration depths, a better boundary condition is given by [13]:


Where the conductivity σ=σcis real for conventional conductors. These approximations have been verified for practical metallization thicknesses by comparison with rigorous mode matching result.

For superconductors, a complex conductivity of the form σ=σn(T/Tc)4i(1(T/Tc)4)/ωμ0λ02,where σn is often associated with the normal state conductivity at Tc and λ0 is the effective field penetration depth.

Now, we have the necessary Green’s function, it is relatively straightforward to formulate the moment method solution for the antenna characteristics. The boundary condition at the surface of the microstrip patch is given by:


Here E¯incand E¯scatare tangential components of incident and scattered electrics fields. Electric field is enforced to satisfy the impedance boundary condition on the microstrip patch, and the current vanishes off the microstrip patch, to give the following set of vector dual integral equations [3]. For superconducting annular ring microstrip antenna, we have:

 en(ρ)=0dkρkρH¯n(kρρ)(G¯(kρ)Z¯S)K¯n(kρ)=0                    a  ρ<b           E29
             Kn(ρ)=0dkρkρH¯n(kρρ)K¯n(kρ)=0                          ρ>b,   ρaE30

And for superconducting circular microstrip antenna, we have:

 en(ρ)=0dkρkρH¯n(kρρ)(G¯(kρ)Z¯S)K¯n(kρ)=0                     ρ<a           E31
            Kn(ρ)=0dkρkρH¯n(kρρ)K¯n(kρ)=0                          ρ>aE32



Galerkin’s method is employed to solve the coupled vector integral equations of (29)-(32). Substituting the Hankel transform current expansion of (15) into (29) and (32). Next, multiplying the resulting equation by ρψnp+(ρ)(p=1, 2,....... P) and ρφnq+(ρ)(q=1, 2,....... Q) and using Parseval’s theorem for VHT, we obtain a system of Q+ P linear algebraic equations for each mode n which may be written in matrix form. Following well known procedures, we obtain the following system of linear algebraic equations:


Each element of the submatrices Z¯CDis given by:


Where C and D represent either ψ or φ, for every value of the integer n.

The integration path for the integrals of (32) is, in general, located in the first quadrant of the complex plane kρ. This integration path must remain above the pole and the branch point of G¯. Although other choices of branch cut are possible, the choice made in this paper is very convenient, this treated in section 3. Once the impedance and the resistance matrices have been calculated, the resulting system of equations is then solved, for the unknown current modes on the microstrip patch. Nontrivial solutions can exist, if the determinant of Eq. (35) vanishes, that is:


In general, the roots of this equation are complex numbers indicating, that the structure has complex resonant frequencies (f=fr+ifi). The bandwidth of a structure operating around its resonant frequency, can be approximately related to its resonant frequency, though the well-known formula (BW=2fi/fr). Once the problem is solved for the resonant frequency, far field radiation, co-polar and cross- polar fields in spherical coordinates are given from.



V¯(kρ)=[1       0       -kρkz 0        1           0]andT¯(θ)=[cosθ      0        -sinθ0            1             0]

The losses in the antenna comprise dielectric loss Pd, the conductor loss Pc, and the radiation loss Pr, are given by [4-6]:

Pc=Zs|Hφ2+Hρ2| dSE38

The efficiency of an antenna can be expressed by:


3. Convergence and comparison of numerical results

Computer programs have been written to evaluate the elements of the impedance, resistance matrices, and then solve the matrix equation (35). To enhance the accuracy of the numerical calculation, the integrals of the matrix elements (33) are evaluated numerically along a straight path above the real axis with a height of about 1k0 (Figure 2). In this case, the effects of the surface waves are included in the calculation and knowledge of the pole locations is not required, while the length of the integration path is decided upon by the convergence of the numerical results. The time required to compute the integral depends on the length of the integration path. It is found that length of the integration path required reaching numerical convergence at 100k0, also Muller’s method that involves three initial guesses, is used for root seeking of (35).

Figure 2.

Integration path used for computing the integrals in the complex kρ plane

To check the correctness of our computer programs, our results are compared with results of other authors. The comparisons are shown in Table 1 for imperfectly conducting microstrip annular ring antennas. The resonant wave number times the inner radius of the ring is kra(kra=2πfraεxε0μ0), as functions of different sizes of the ratio of the substrate thickness d normalized by the inner radius a is fixed of (0.71cm), and an outer radius of b=2a, the relative permittivity was ε1z1x=2.65. Annular ring microstrip antenna is excited in the TM11 and TM12 modes. We found that, for the TM11 mode, the real part of (kra) increases as d/a increases. At the same time, for this mode, the imaginary part of (kra), which includes the losses by radiation of the structure, is approximately zero. This means that, the TM11 mode has narrow bandwidth and weak-radiation. In addition, we observe that, for the TM12 mode, the real part of (kra) decreases as (d/a) increases, the TM12 has relatively wide bandwidth and high-radiation. Therefore, the TM11 mode is good for resonator applications and the TM12 mode for antennas. It is for that one does the applications of the annular ring microstrip antenna in the TM12 mode, better than in the TM11 mode. Thus, for this reason, the considered mode in this work is the TM12 mode. It is clear from Table 1 that our results agree very well with results obtained by other authors [14, 12].

d/aMode TM12Mode TM11
Results of
Results of
Our resultsResults of [14]Our results
Im (kra)Re
0.0053.260.0023.240.0023.2570.0020.671,6.10-40.6761,6. 10-4,7. 10-40.6821,8. 10-4,4. 10-40.6955,5. 10-4

Table 1.

Calculated (Kra) of annular ring microstrip antennas.

In table 2, we have calculated the resonant frequencies for the modes (TM11, TM21, TM31, and TM12) for perfect conducting circular patch with a radius 7.9375mm, is printed on a substrate of thickness 1.5875mm. These values are compared with theoretical and experimental data, which have been suggested in [10]. Note that the agreement between our computed results, and the theoretical results of [10], is very good.

In our results, we need to consider only the P functions of ψnp, and the Q functions of φnq. The required basis functions for reaching convergent solutions of complex resonant frequencies, using cavity model basis functions are obtained with (P=5, Q=0).

ModeResults of [10 ]Our results
Frequency (GHz)
Quality factors
Frequency (GHz)
Quality factors

Table 2.

Comparison of resonant frequencies of the first four modes of a perfect conducting circular patch printed on a dielectric substrate (a=7.9375mm, εx=εz=2.65,d=1.5875 mm).

4. Resonant frequency of superconducting patch antenna

Shown in Figures 3-4, is the dependence of resonant frequency on the thickness t of superconducting patch of the antennas. It is observed that, when the film thickness (t) increases, the resonant frequency increases quickly until, the thickness t reaches the value penetration depth (λ0). After this value, the increase in the frequency of resonance becomes less significant.

Figure 3.

Real part of resonant Frequency against thickness of superconducting annular ring patch antenna. (d=254µm, a=815µm, b=2a, T/ TC=0.5, λ0=1500Å, σn=210S/mm). (–––– ) εx= 9.4, εz= 11.6;(– – – ) εx= 11.6, εz= 11.6; (―∙―∙) εx= 13, εz=10.3; (……) εx=10.3, εz= 10.3.

Figure 4.

Real part of resonant Frequency against thickness of superconducting circular patch antenna. (d=254µm, a=815µm, TC =89K, T/ TC=0.5, λ0=1500Å, σn=210S/mm).(–––– ) εx= 9.4, εz= 11.6;(– – – ) εx= 11.6, εz= 11.6; (―∙―∙) εx= 13, εz=10.3; (……) εx=10.3, εz= 10.3.

Figures 5-6 demonstrated relations between the real part of frequency resonance, and the normalized temperature (T/Tc), where the critical temperature used here for our data (89K). The variations of the real part of frequency due to the uniaxial anisotropy decrease gradually with the increase in the temperature. This reduction becomes more significant for the values of temperature close to the critical temperature. These behaviours agree very well with those reported by Mr. A. Richard for the case of rectangular microstrip antennas [2].

Figure 5.

Real part of resonant frequency of superconducting annular ring patch antenna against T/TC (d=254µm, a=815µm, b=2a, t=0.02µm, λ0=1500Å, σn=210S/mm). (–––– ) εx= 9.4, εz= 11.6;(– – – ) εx= 11.6, εz= 11.6; (―∙―∙) εx= 13, εz=10.3; (……) εx=10.3, εz= 10.3.

Figure 6.

Real part of resonant frequency of superconducting circular patch antenna against T/TC (d=254µm, a=815µm, TC =89K, h= 0.02µm, λ0=1500Å, σn=210S/mm ). (–––– ) εx= 9.4, εz= 11.6;(– – – ) εx= 11.6, εz= 11.6; (―∙―∙) εx= 13, εz=10.3; (……) εx=10.3, εz= 10.3.

5. Radiations patterns and efficiency of superconducting patch antenna

The calculated radiations patterns (electric field components, Eθ; Eφ), of the microstrip antenna on a finite ground plane in the E plane, and in the H plane are plotted in Figs. 7 -10, printed on an uniaxial anisotropy substrate thickness (d=254μm). The mode excited for superconducting annular ring patch antenna is the TM12 and for superconducting circular patch antenna is the TM11. It is seen that the permittivity εz has a stronger effect on the radiation than the permittivity εx. The radiation pattern of an antenna becomes more directional as its εz increases. Another useful parameter describing the performance of an antenna is the gain. Although the gain of the antenna is related to the directivity, the gain of an antenna becomes high as its εz increases.

Figure 7.

Radiation pattern versus angle θ of superconducting annular ring patch antenna at φ=0° plane (a=815µm, b=2a, t=0.02µm,d=254µm, λ0=1500Å, σn=210S/mm). (–––– )εx= 9.4, εz= 11.6; ( – – – )εx= 11.6, εz= 11.6;(― ∙ ― ∙) εx= 13, εz=10.3 ;(……) εx=10.3, εz= 10.3.

Figure 8.

Radiation pattern versus angle θ of superconducting annular ring patch antenna at φ=90° plane (a=815µm, b=2a, t=0.02µm, d=254µm, λ0=1500Å, σn=210S/mm).(–––– )εx= 9.4, εz= 11.6;( – – – )εx= 11.6, εz= 11.6; (― ∙ ― ∙) εx= 13, εz=10.3 ; (……) εx=10.3, εz= 10.3.

Figure 9.

Radiation pattern versus angle θ of superconducting circular patch antenna at φ=0° plane (a=815µm, t=0.02µm, T/ TC=0.5, d=254µm, λ0=1500Å, n=210S/mm). (–––– ) εx= 9.4, εz= 11.6; ( – – – ) εx= 11.6, εz= 11.6; (― ∙ ― ∙) εx= 13, εz=10.3 ; (……) εx=10.3, εz= 10.3.

Figure 10.

Radiation pattern versus angle θ of superconducting circular patch antenna at φ=π/2 plane (a=815µm, h=0.02µm, T/ TC=0.5, d=254µm, λ0=1500Å, n=210S/mm). (–––– ) εx= 9.4, εz= 11.6; ( – – – ) εx= 11.6, εz= 11.6; (― ∙ ― ∙) εx= 13, εz=10.3 ; (……) εx=10.3, εz= 10.3.

In calculation of losses, we have found that, the values of dielectric loss (Pd), the conductor loss (Pc), and the radiation loss (Pr) will depend on frequency. We use results precedents to calculate the variation of radiation efficiency as a function of resonant frequency, for various isotropic dielectric substrates. Our results are shown in Figs. 11 and 12. It is seen that the efficiency increases with decreasing frequencies. The same behaviour is found by R. C. Hansen [1].

Figure 11.

Superconducting annular ring patch antenna efficiency for the mode TM12 (d=254µm, a=815µm, b=2a, t=0.02µm, λ0=1500Å, σn=210S/mm, εx= εz=11.6, δ=0.0004).

Figure 12.

Superconducting circular patch antenna efficiency for the mode TM11 ( (d=254µm, a=815µm,h=0.02µm, T/ TC=0.5, λ0=1500Å, σn=210S/mm, εx= εz=11.6, δ=0.0024).


6. Conclusion

This work presents a fullwave analysis for the superconducting microstrip antenna on uniaxial anisotropic media. The complex resonant frequency problem of structure is formulated in terms of an integral equation. Galerkin procedure is used in the resolution of the electric field integral equation, also the TM, TE waves are naturally separated in the Green’s function. In order to introduce the effect of a superconductor microstrip patch, the surface complex impedance has been considered. Results show that the superconductor patch thickness and the temperature have significant effect on the resonant frequency of the antenna. The effects of a uniaxial substrate on the resonant frequency and radiation pattern of structures are considered in detail. It was found that the use of such substrates significantly affects the characterization of the microstrip antennas, and the permittivity (εz) along the optical axis has a stronger effect on the radiation of antenna. Thus, microstrip superconducting could give high efficiency with high gain in millimeter wavelengths. A comparative study between our results and those available in the literature shows a very good agreement.

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

Ouarda Barkat (March 6th 2013). Full-Wave Spectral Analysis of Resonant Characteristics and Radiation Patterns of High Tc Superconducting Circular and Annular Ring Microstrip Antennas, Advancement in Microstrip Antennas with Recent Applications, Ahmed Kishk, IntechOpen, DOI: 10.5772/54665. Available from:

chapter statistics

2425total 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

Advancement in Microstrip Antennas with Recent Applications

Edited by Ahmed Kishk

Next chapter

Compact Planar Multiband Antennas for Mobile Applications

By Ahmad Rashidy Razali, Amin M Abbosh and Marco A Antoniades

Related Book

First chapter

The Generalized Solutions of a System of Maxwell's Equations for the Uniaxial Anisotropic Media

By Seil Sautbekov

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