Novel Applicators for Local Microwave Hyperthermia Based on Zeroth-Order Mode Resonator Metamaterial
It is demonstrated that a theory of zero-order mode resonator (ZOR) metamaterial (MTM) structure can be used for the development of a novel class of applicators for microwave thermotherapy, for example, for hyperthermia in cancer treatment or for physiotherapy. The main idea of creating such an applicator is to generate and radiate a plane electromagnetic (EM) wave into the treated biological tissue, at least in a certain extent. The main aim of this paper is to investigate whether an EM wave generated by ZOR MTM structure and emitted into the biological tissue can produce a homogeneous SAR distribution in the planes parallel to the applicator aperture and achieve a penetration depth approaching the theoretical limit represented by SAR distribution and penetration depth of an ideal EM plane wave. EM field distribution inside a virtual phantom of the treated region generated by the applicator that is based on the proposed ZOR MTM principle is investigated using a well-proven full-wave commercial simulation tool. The proposed applicator type shows both a low unwanted leaked electromagnetic field and a fairly homogeneous electric field in its aperture as well as in the virtual phantom of the treated region.
The main aim of this paper is to verify whether it is feasible to apply principles of ZOR MTM radiating structures in the design of efficient applicators for radiofrequency (RF)/microwave (MW) thermotherapy, especially for hyperthermia cancer treatment and for physiotherapy. For this aim, MTM antennas inspired by those described in [1, 2] will be investigated and adapted with respect to requirements of effective hyperthermia treatment (i.e., we will study the homogeneity of SAR distribution in the treated region). In order to achieve an optimal waveform shape of the radiated EM wave, applicators have to be designed in a certain way. The best possible results can be obtained by comparing various shapes of waveforms of EM waves (e.g., plane wave, cylindrical wave, and spherical wave) for local and deep local treatment waveform of plane wave. This can ensure(i)the best possible value of the effective treatment depth;(ii)the best possible homogeneity of SAR 3D distribution (i.e., its distribution on the surface and in the whole volume of the area to be treated). SAR distribution inside the virtual phantom of the treated region, created using the proposed applicator, and the reflection coefficient of the proposed applicator are investigated using the well-proven full-wave commercial simulation tool COMSOL Multiphysics. Based on our previous experience [1–3], an excellent agreement between simulation and measurement results can be expected. That is why the conclusions made in this paper are based on the results of numerical simulations only.
2. MTM Zeroth-Order Mode Resonator
The concept of MTM phenomenon was first comprehensively introduced by Veselago in 1968 . In the aforementioned publication, he speculated on the existence of materials whose permittivity and permeability were simultaneously negative. He named these materials left-handed (LH), as the , , and vectors of the studied EM wave formed a left-handed triad if the wave propagated through such environment. The first experimental verification of MTM phenomena was performed by a research group at the University of California, San Diego (UCSD), in 2001 .
An infinitesimally short lossless transmission line (TL) section can be described by a simple equivalent circuit consisting of a series inductor and a shunt capacitor . The lossless MTM cell implemented in planar technology consists of a TL section with artificially inserted series capacitors and shunt inductors (with subscript denoting its left-handed properties). The equivalent circuit of the MTM cell can be then represented by four-lumped elements (as shown in Figure 1) .
Several different EM radiating structures based on the MTM principle were introduced in the past [7–10]. Since the very beginning of the development, however, real implementation possibilities of such antennas, for example, in communication technology, have been very limited because of their poor radiation efficiency. MTM antennas with very good radiation efficiency were first presented in [1, 2].
3. Design of the Applicator Based on ZOR MTM Structure
In this section, we will study the possibility of creating MTM based applicators for local microwave hyperthermia cancer treatment at a frequency of 434 MHz. The basic part of the proposed applicator will consist of the ZOR.
The working idea of ZOR is based on a special case of resonance that can occur when the TL meets the conditions of the MTM phenomenon. The phase constant is at working frequency in this case, which implies infinite guided wavelength along the MTM structure as well as zero phase shift . It is very important to note that this phenomenon enables creation of a very special kind of resonator whose physical length is completely independent of the classical resonance condition (i.e., required to be a natural number multiple of the half working wavelength in case of either open-circuited or short-circuited TL) . The typical voltage wave distribution along the resonant length for the negative and zero resonances is shown in Figure 2. In the case of zeroth-order mode resonance , the value of the voltage along the ZOR is constant. A more comprehensive description of ZOR properties can be found in .
The zeroth-order resonant frequency of the proposed ZOR MTM applicator is equal to the working frequency of a particular RF/MW hyperthermia system. Please note that the ideas and principles of MTM applicator design described here are expected to be valid at any usual working frequency of RF/MW hyperthermia systems. The initial idea of the mechanical and EM arrangements of the ZOR MTM applicator proposed here is displayed in Figure 3. In order to obtain high radiation efficiency and thus the best possible transfer of EM energy to the treated area, we can use the main ideas, experiences, contributions, and conclusions described in [1–3]. That means to combine the ZOR MTM structure with relatively long inductive elements along which electric current with equal phase will flow. These currents will then excite an EM wave propagating into the area to be treated. Optimal choice of positions of these inductive elements then enables us to approximate the preferred waveform (we want to excite the waveform of a plane wave, but the waveform of either a cylindrical or a spherical wave can be obtained as well). This will enable us to approximate the possible shape as well as the dimensions of the treated area. For the above mentioned main aim, some results and experiences obtained by applicators described in [12–16] can be used.
The design rules of this structure are the same as in . Thanks to the excitation of zeroth-order mode vectors of surface current density in all vertical parts of the antenna (including feeding), all these surface currents are in phase; that is, the radiated contributions from all individual vertical parts are in a very good superposition in the applicator aperture (i.e., the Huygens principle phenomenon can be applied to describe the resulting EM field distribution in the area to be treated). For the design of the ZOR MTM applicators, the following dimensions of antenna elements were chosen. In this special discussed case, the length of the vertical part of the antenna is equal to . The longitudinal dimension of the unit cell can vary (or be adjusted) in a relatively large range . Depending on the selected size of this dimension, it is necessary to adjust the dimensions of interdigital capacitors. For the presented applicator, the overall length of one unit cell is , the length of figures of the interdigital capacitor equals , and the width and the gap of the interdigital capacitor are ; please note that the same notation has been used as in . The overall physical dimensions of the applicator prototype are 268 × 68 × 70 mm which represents relative dimensions (with respect to wavelength of plane EM wave in vacuum at working frequency) as follows: the relative width of the aperture is equal to , the relative height of the aperture is , and the relative depth of the applicator is . Thickness, relative permittivity, and equivalent conductivity of the considered substrate are 1.5 mm, 4.3, and 0.001 S/m, respectively. Alternatively, the ZOR MTM radiating structure can be inserted into the rectangular waveguide section to ensure that the whole radiated power would be perfectly directed to the biological tissue in the area to be treated. Another possible and interesting use for real clinical therapy is the case where the ZOR MTM radiating structure would be surrounded by metal plates from the top and the back sides and the lateral sides would be made of a dielectric substrate material to ensure that we would not excite the dominant rectangular waveguide mode . One of the basic features of the ZOR MTM radiating structure, observed during our investigations, is that it resonates at almost the same frequency irrespective of whether the structure is surrounded by metal from the back and lateral sides. Only if we placed the metal plane closer to the top of ZOR, the capacity of interdigital capacitors would change and thus the resonant frequency could be shifted from the working frequency. This phenomenon can be used for frequency adjusting or impedance matching of this applicator if necessary.
4. Distribution of SAR Created by ZOR MTM in the Treated Area
To study and verify how these ZOR MTM structures would radiate into the biological tissue in the area to be treated, several numerical simulations of the discussed case were performed. SAR distribution inside a virtual phantom of the treated region generated by the proposed applicator was investigated using COMSOL Multiphysics . Muscle tissue dielectric parameters were considered as follows: the real part of the complex permittivity is equal to and the equivalent electric conductivity is equal to . In the studied model displayed in Figure 4, a water bolus with a thickness of 2 cm was placed between the applicator aperture and the virtual phantom in the area to be treated. The resulting SAR distribution is displayed in two planes in Figures 4 and 5. Figure 4 displays SAR distribution in the plane perpendicular to the applicator aperture, illustrating how deep EM power penetrates into the area to be treated. It can be observed that when the EM power penetrates the biological tissue, it has very good SAR homogeneity and that the penetration depth is approaching the theoretical limit .
Figure 5 displays SAR distribution in the plane parallel to the applicator aperture. Homogeneous absorption of EM at the surface of the area to be treated can be observed. Similar distribution of SAR will be observed in all planes parallel to the applicator aperture, but the level of SAR will decrease with increasing depth.
5. Improvement of SAR Homogeneity by Aid of Symmetric Feeding of ZOR MTM
The ZOR MTM applicator discussed in this paper and the distribution of SAR achieved by the use thereof (displayed in Figures 4 and 5) can be considered very suitable for practical treatment of cancer patients. As already mentioned in this paper, it approaches the homogeneity level and the effective treatment depth of the plane wave. However, a certain level of asymmetry of SAR distribution is evident in Figures 4 and 5. This effect can be explained by asymmetrical feeding of the ZOR MTM applicator displayed in Figures 5 and 6 and also by the fact that the surface current density on the feeding vertical part has the same phase as the surface current density on the other vertical parts.
The contribution to radiation is not as it could be if the feeding part would be in the same line as the others. In this part, a modified mechanical and EM arrangement of the ZOR MTM is proposed, which would help us improve the homogeneity of SAR distribution in the treated area. As can be seen in Figure 6, a modified ZOR MTM applicator with symmetric feeding is proposed. Furthermore, another vertical radiating part connected to the feeding point via an interdigital capacitor and a section of microstrip TL that was added to ensure the same phase shift as that of the other vertical parts.
The working principle of this ZOR MTM structure is the same as in the previous case. Thanks to the excitation of the zeroth-order mode, the vectors of surface current density on all vertical parts of the antenna including the feeding are in phase. Again, we can say that the phenomenon of the Huygens principle can be applied to describe the resulting EM field distribution in the area to be treated. Therefore, the radiated contributions from all vertical parts constructively interfere in the applicator aperture. The dimensions of the vertical parts of the antenna were chosen to be again. The longitudinal dimension of the unit cell can be varied (be adjusted) in a relatively large range . Based on the choice of this dimension, it is necessary to adjust the dimensions of interdigital capacitors.
In Figures 6 and 7, it is evident that the homogeneity of SAR in front of five symmetrically fed radiating parts is much better than the one described in the previous case (chapters 3 and 4). Figure 6 displays SAR distribution in the plane perpendicular to the applicator aperture, illustrating how deep EM energy penetrates the area to be treated. It can be observed that when the EM wave penetrates the biological tissue, it has very homogeneous SAR distribution and that both the SAR distribution and the depth of penetration of the treated area approach the theoretical limit of a plane EM wave.
Figure 7 displays SAR distribution in the plane parallel to the applicator aperture, illustrating the homogeneity of EM energy absorption at the surface of the area to be treated. Similar distribution of SAR will be observed in all planes parallel to the applicator aperture, but the level of SAR will decrease with increasing depth.
6. Achieved Waveform and Comparison with Theoretical Limits
In this section, the achieved wave form in the biological tissue is compared to the theoretical limit represented by plane EM wave.
SAR distribution in Figure 7 shows very good homogeneity. In Figure 8, SAR distribution along three line segments (cut lines) which lie in the plane perpendicular to the water bolus/treated area interface is depicted. The plane cuts the inductive posts (vertical parts) in the middle. Cut lines 1, 2, and 3 start at the interface, lead into the treated area, and are located between the first and the second inductive posts, between the fourth and the fifth inductive posts, and in front of the third inductive post, respectively. SAR distribution is compared to that of plane EM wave propagating through the treated area considered here (muscle tissue). SAR values are given as percentages of the SAR value of the plane EM wave at the depth of 1 cm. Along all three cut lines, the SAR approximates very well the shape of the exponential form related to plane EM wave. In Figure 9, SAR distribution along three line segments (cut lines) which lie in the plane perpendicular to the water bolus/treated area interface is plotted. The plane cuts the inductive posts in the middle. The cut lines are parallel to the interface. Cut lines 4, 5, and 6 are located at the interface, at a depth of 1, 2, and 3 cm, respectively. SAR values are given as percentages and are related to the highest SAR value along cut line 5. At the depth of 1 cm, SAR distribution is fairly homogeneous. The heat produced by the hotspots at the surface has to be removed using the water bolus, a common strategy used in RF/MW hyperthermia. Similarly, in Figure 10, SAR distribution along three line segments (cut lines) which lie in the plane perpendicular to the water bolus/treated area interface is plotted. The plane is parallel and passes through the third inductive post. The cut lines are parallel to the interface. Cut lines 7, 8, and 9 are located at the interface, at a depth of 1, 2, and 3 cm, respectively. SAR values are given in percent and are related to the highest SAR value along cut line 8.
In Figures 11 and 12, contours of SAR values as percentages in two main planes are plotted, in the plane parallel to the surface of the treated biological tissue in the depth of 1 cm in the tissue and in the plane perpendicular to this plane, cutting the inductive posts in the middle. 25, 50, and 75% SAR levels are plotted, related to the highest SAR value in the depth of 1 cm in the tissue. Usually, the 50% SAR contour is considered to bound the area where the treatment is efficient. From Figure 11, it follows that the surface area of a rectangular effective treatment area (a rectangle inscribed in the 50% SAR contour) for the proposed applicator is about 6 by 24 cm2.
7. Ongoing Research of ZOR MTM Applicators
Our ongoing research of ZOR MTM applicators focuses on designing several practical clinical applicators and on preparing and performing basic experimental evaluations of these structures in order to verify the promising results presented here, obtained by numerical simulations only. We explore two basic structure types:(i)ZOR MTM structure inserted in the waveguide (thus gaining the advantages of waveguide applicators);(ii)ZOR MTM structure used as a flat planar applicator.
In this paper, a novel principle for the design of applicators based on ZOR MTM structures has been proposed. It has been demonstrated here that when penetrating biological tissue EM waves generated by the proposed applicators generate very good SAR homogeneity and achieve penetration depth approaching, to a certain extent, the theoretical limit, that is, close to that of the EM plane wave. The surface of a rectangular effective treatment area (a rectangular inscribed into 50% SAR contour) is for the proposed applicator about 6 × 24 cm2.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This research has been supported by the research program of the Grant Agency of the Czech Republic, Project no. 13-29857P Human Body Interactions with EM Field Radiated by Metamaterial Structures and Project no. 14-00386P Study of Thermal and Nonthermal Effects of High-Power EM Field Structure of Matter.
D. Vrba and M. Polívka, “Radiation efficiency improvement of zeroth-order resonator antenna,” Radioengineering, vol. 18, no. 1, pp. 1–8, 2009.View at: Google Scholar
M. Polívka and D. Vrba, “Shielded micro-coplanar CRLH TL zeroth-order resonator antenna: critical performance evaluation,” Radioengineering, vol. 18, no. 4, pp. 368–372, 2009.View at: Google Scholar
D. Vrba, Electrically miniaturized antennas based on Zeroth-order resonance [Ph.D. thesis], Czech Technical University in Prague, Faculty of Electrical Engineering, Prague, Czech Republic, 2012.
D. M. Pozar, Microwave Engineering, John Wiley and Sons, New York, NY, USA, 2nd edition, 1998.
H. R. Raad, A. I. Abbosh, H. M. Al-Rizzo, and D. G. Rucker, “Flexible and compact AMC based antenna for telemedicine applications,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 2, pp. 524–531, 2013.View at: Google Scholar
C. Caloz and T. Itoh, Electromagnetic Metamaterials: Transmition Line Theory and Microwave Applications, John Wiley and Sons, New York, NY, USA, 2006.
C. Franconi, J. Vrba Jr., and F. Montecchia, “27 MHz hybrid evanescent-mode applicators (HEMA) with flexible heating field for deep and safe subcutaneous hyperthermia,” International Journal of Hyperthermia, vol. 9, no. 5, pp. 655–673, 1993.View at: Google Scholar
P. Togni, Z. Rijnen, W. C. M. Numan et al., “Electromagnetic redesign of the HYPERcollar applicator: toward improved deep local head-and-neck hyperthermia,” Physics in Medicine and Biology, vol. 58, pp. 5997–6009, 2013.View at: Google Scholar
M. Polivka and D. Vrba, “Input resistance of electrically short not-too-closely spaced multi- element monopoles with uniform current distribution,” IEEE Antennas and Wireless Propagation Letters, vol. 11, no. 1, pp. 1592–1595, 2012.View at: Google Scholar
COMSOL AB, COMSOL Multiphysics User's Guide, COMSOL AB, Stockholm, Sweden, 4th edition, 2010.
D. Andreuccetti, R. Fossi, and C. Petrucci, “An internet resource for the calculation of the dielectric properties of body tissues in the frequency range 10 Hz–100 GHz,” March 2013.View at: Google Scholar
J. Vrba, C. Franconi, and M. Lapes, “Theoretical limits for the penetration depth of intracavitary applicators,” International Journal of Hyperthermia, vol. 12, no. 6, pp. 737–742, 1996.View at: Google Scholar