#### Abstract

The chosen rectangular and fractal microstrip patch antennas above an infinite ground plane are analyzed by the theory of characteristic modes. The resonant frequencies and radiation are evaluated. A novel method by Vandenbosch for rigorous evaluation of the radiation is employed for modal currents on a Rao-Wilton-Glisson (RWG) mesh. It is found that the resonant frequency of a rectangular patch antenna with a dominant mode presents quite complicated behaviour including having a minimum at a specific height. Similarly, as predicted from the simple wire model, the radiation exhibits a minimum too. It is observed that the presence of out-of-phase currents flowing along the patch antenna leads to a significant increase of the factor.

#### 1. Introduction

Evaluation of the basic properties of microstrip patch antennas (MPA) has been numerously discussed in literature, see, for example, [13]. The two main MPA attributes are resonant frequency (or frequencies but we will deal mostly with the dominant mode) and the radiation factor. So far, only approximate results and semianalytic equations have been published. To our knowledge, this is the first time that these important characteristics have been studied in a rigorous way. The antennas are treated by using a modal approach (hence we do not a priori consider any feeding to be connected), namely, by the theory of characteristic modes (TCM), [4, 5]. Evaluation of the radiation is performed both by the TCM from the eigenvalues slope and by novel rigorous equations derived by Vandenbosch [6] and Vandenbosch and Volski [7].

#### 2. The Theory of Characteristic Modes

For completeness, let us formulate the basics of the characteristic modes for perfectly conducting bodies of area . The scattered field is related to the electric surface currents by the electric field integral equation (EFIE) [8] Equation (1) is usually treated within the method of moments (MoMs) [8] framework and, due to the structure discretization, the operator is known as the “complex impedance matrix” .

Then the associated Euler’s equation to be solved is Equation (2) is a standard weighted eigenvalue equation leading to a set of real characteristic eigencurrents and associated eigenvalues . Properties of eigenvalues are described in [9], at this moment it is important to note that   reflects the amount of net reactive power (thus means resonance). Instead of eigenvalues, the so-called characteristic angles are introduced to show more visible behavior with frequency [9]. Characteristic currents form a complete orthogonal set, and hence the total current on a conducting body may be expressed as a linear combination of these mode currents [10].

##### 2.1. Implementation of the Characteristic Modes Theory

Implementation of the modal decomposition process has been done in the MATLAB [11] environment using Makarov EFIE codes [12] with the RWG basis functions [13]. This usage is restricted to arbitrary 3D PEC structures with air dielectrics. Our developed TCM tool [14] has the following main advantages:(i)Comsol Multiphysics [15]/MATLAB’s PDE TooLbox mesh import,(ii)Optional Green’s function for infinite ground plane simulations,(iii)Single solver/multicore solver/distributed solver (within a computer network with installed MATLAB).

In [6] a novel theory able to rigorously calculate radiated power and stored energies directly from currents flowing along the antenna has been presented. The radiation factor is then readily evaluated by the definition [16]: The equations for radiated power and stored electric and magnetic energies , are where where is a free-space wavenumber, is the surface current density, and is the distance between interacting current elements. The tilde denotes that the radiation contribution has been subtracted from the stored energies at every point in space [17]. It is assumed that the currents are flowing in a vacuum.

##### 3.1. The Modal Radiation Factor

The modal radiation factor may be evaluated from the slope of modal eigenvalues [18]: In [18], (10) is supposed to be an approximation of the radiation , but in resonance it is actually exact.

Since characteristic modes are normalized to radiate unit power [4], (3) reduces to For parallel or series RLC circuit (hence, for one mode), the “impedance ” equals the exact “current ” [6]: Inserting valid for lossless antennas [19] and using the fact that , (12) results in providing that It is therefore concluded that the modal equals the by definition, and it can be proven (using the reactance theorem [20, 21]) that in resonance it also equals the radiation defined from energies by (11).

##### 3.2. Software Implementation

The above equations were implemented in MATLAB for the RWG triangular mesh where two different interaction situations occur:

(a) Distant Elements
When the triangular elements are not overlapping, current density on triangles may be simply approximated as point sources located at the centre of triangles [22], see Figure 1. No actual integration is then needed. This centroid approach is very fast with satisfactory accuracy as will be shown later (however it may fail for patches located very close to the ground plane).

(b) Overlapping (Self) Elements
As known from the method of moments, the so-called “self” contributions are of great importance when dealing with calculations on discrete elements (meshes).
Here, the self-interaction occurs when two triangles are overlapping each other. Due to the behavior of integral kernels, only rapidly varying term has to be carefully treated. Since ( being the longest side of the triangle ) is satisfied, one needs only to use the first term in the Taylor series expansion. The dominant singular static part is and the integral to be worked out is where is a triangular area. Using simplex coordinates transformation (Figure 2), the result is [23, 24] where is the triangle area, are the edge lengths (see Figure 2), and is the perimeter of the triangle.

#### 4. Applications: Rectangular Patch Antenna

Let us first concentrate on a rectangular patch antenna of dimensions  mm and  mm (further noted as R50 × 30) placed in air at a height above an infinite ground plane. Only the dominant TM01 mode will be studied. The reason for choosing a patch with is that we do not have to deal with degenerated modes.

Using the image theory, the radiator in the plane at height above an infinite electric ground plane is modelled as two patches separated by . The total number of triangular elements is 676. In the TCM analyser, a proper out-of-phase mode is selected (Figure 3).

The resonant frequency of the dominant mode is shown as a function of height , see Figure 4. It has been evaluated from a modal resonant condition for eigenvalues employing an adaptive frequency sweep for each height. The behaviour is quite peculiar, especially for greater heights. For low heights ( mm or ), the resonant frequency decreases “regularly,” and quasianalytical formulas (see, e.g., [1, 3]) based on the fringing field concept are valid below this range. For mm () there is absolute minimum of the TM01 resonant frequency. Further on, the resonant frequency rises to reach its maximum for mm (). Around this specific height the patch also shows the minimum of the radiation . The above described process repeats periodically. It is yet unclear to the authors as what is the physical background to the resonant frequency discontinuity around .

The terms , and obtained from (5)–(9) and eigenvalues are plotted at Figure 5 for  mm as a function of frequency. There is excellent agreement between the difference in stored energies and the eigenvalues, both obtained in a completely different manner.

There is also very good agreement between the exact and confirming the validity of the proposed algorithm via (14), see Figure 6. Note that in (10) does not require the currents to be calculated on the structure while is evaluated in a rigorous way from modal currents (11).

From Figure 6 it is seen that the radiation has a minimum for a specific height. It is deduced that the reason lies in the cancelling of the radiated power between the two out-of-phase currents. Similar behaviour has been observed in the case of two half-wave thin-wire dipoles with opposite sinusoidal currents, separated by , see [25] for details. Actually these two out-of-phase dipoles may serve as a very simple model for a patch antenna with a dominant mode. When the dipoles are reduced to elementary (Hertzian) ones, an approximate analytical solution is available and in [25] we showed that the is led by the function After deriving (18), the condition is worked-out and the first nontrivial root of (19) could be approximated as [25] For sinusoidal currents on dipoles the minimum (evaluated numerically) occurs for .

The minimum of the patch under study is obtained at , a value that is remarkably close to the simple dipole model.

##### 4.1. Algorithm Convergence

Since no other methods for calculating modal are available, is taken as a reference, and the relative error percentage is defined as: where is calculated from the currents using (11). Four different heights were chosen,  mm (),  mm (),  mm (), and  mm (), and the relative error was evaluated as a function of total triangular elements (including the mirror), see Figure 7. All quality factors were evaluated at the resonant frequency of the dominant mode for the R50 × 30 patch. As discussed earlier, the centroid approximation became more inaccurate with low heights . However, even for the lowest analyzed value , the relative error is in the order of a few percent for reasonable mesh density (hundreds of elements). Further improvements to the integration routine are considered for the future.

##### 4.2. Fractional Bandwidth of the R50 × 30 Patch Antenna

It is known that the fractional bandwidth (FBW) is related to the unloaded factor and the desired matching VSWR level. For VSWR < we have [26] Using a full-wave simulator CST-MWS [27], an R50 × 30 patch has been simulated and the FBWCST for VSWR < 2 was calculated as: where and are margins for VSWR < 2 and is the centre frequency. Only very low heights were studied since we used a simple probe feed which introduces an inductance component to the total input impedance. The comparison in Figure 8 shows good agreement of both fractional bandwidths.

#### 5. Applications: Fractal Antennas

In this section, a bit more complex structures will be studied. The first one (the “Self Affine U” fractal, SAU), has been described in [28] and further analyzed in [29]. This kind of radiating motif is employed as a dual-band radiator with mutually orthogonal radiation patterns at both bands. Therefore we are analyzing the first two modes, where the currents are orthogonal. These are depicted in Figures 9 and 10 for first (SAU1) and second (SAU2) fractal iteration, respectively. The current of the first (lower) mode has two out-of-phase components (see Figure 11 for schematic current paths) while the second mode comprises inphase currents only. As we know from previous studies, opposite currents contribute to a rapid increase of the radiation , and it is expected that will have a much higher than .

Figure 11 presents a very simple concept showing the main current paths for the and modes discussed above including the mirroring effect of the infinite ground plane. It could be simply stated that more opposing current paths lead to significant increase in .

We show detailed behaviour only for SAU2 (the situation is similar for SAU1)—see Figure 12 that confirms high for the mode. Characteristic angles are calculated for  mm, the actual height for which the dual-band antenna was designed [29].

##### 5.1. The FCL-2 Fractal Antenna

The second presented structure is the so-called fractal clover leaf (FCL) of the second iteration, [14]. The antenna is fed by an L-probe [30] that excites its dominant mode and is located at height  mm. Actually, the dominant mode is composed of two degenerated modes and (Figure 13). The second higher mode is shown at Figure 14 for completeness.

Figure 15 shows the main current paths of these modes, and we can again deduce that the dominant mode will exhibit lower compared to . This is confirmed by Figure 16 has more than 200x higher radiation .

#### 6. Resonant Properties of Studied Antennas

The properties of studied antennas are summarized in this section. At first we observed that microstrip antenna could support different kinds of modes regarding their factors (see Figure 17):(a)low modes with the current flowing in one direction and not changing its phase (dominant modes of simple shapes like rectangular, circular patch, and so forth.)(b)high modes with part of the currents flowing in the opposite direction. These modes exist even on simple “U” shaped patch (Figure 9 left) and on complex (fractal) geometries.

Secondly, it has been observed that resonant frequency is quite a complicated function of height. Unfortunately we do not yet have any physical explanation as to why some modes present minimum values of .

Looking at Figure 18, it is clear (and interesting) that the resonant frequency behaves quite differently for low- and high- modes. The resonant frequency of low- modes is much more sensitive to the height, whereas high- modes exhibit almost constant when the height is varied. The proposed explanation is that the opposite currents (responsible for high ) keep reactive fields very close to the radiating structure so the effect of a fringing field coupled to the ground plane becomes almost negligible.

#### 7. Conclusions

Modal resonant properties of selected microstrip patch antennas have been studied with the help of characteristic modes and the novel theory published by Vandenbosch. It has been found that the resonant frequency of a simple rectangular patch antenna is quite a complicated function of its height above the infinite ground. Moreover, the dependency of resonant frequency is also found to be a function of the radiation factor (which is now possible to calculate in a rigorous way). Due to the complexity of the problem, no physical explanation for the resonant frequency behaviour has yet been found.

It is observed that the radiation factor decreases for “standard” heights (<~), however there exists an absolute minimum value of that has already been predicted by simple modeling of two elementary out-of-phase dipoles. Using proper feeding techniques (like the L-probe) allows us to design wideband compact antennas.

The theory now puts current distribution and the radiation factor into objective context. Whenever the current mode exhibits opposite components, high may appear.

Future work is needed to connect the presented theory with parameter sweeps or even optimization, so we will be able to design novel wideband/multimode compact antennas.

#### Acknowledgments

This work was supported by the Grant Agency of the Czech Technical University in Prague, Grant no. SGS11/065/OHK3/1T/13 and by the Project COST 1102. The authors would like to thank professor Vandenbosch for fruitful discussions, Neil Bell for his comments, and two anonymous reviewers who suggested some improvements to the paper.