24468_ftp

Upload: hmalrizzo469

Post on 03-Apr-2018

218 views

Category:

Documents


0 download

TRANSCRIPT

  • 7/28/2019 24468_ftp

    1/7

    technique for efficient solution of method of moments matrix equation,

    Microwave Opt Technol Lett 36 (2003), 95100.

    5. R.G. Ayestaran, M.F. Campillo, and F. Las-Heras, Multiple support

    vector regression for antenna array characterization and synthesis, IEEE

    Trans Antennas Propag 55 (2007), 832840.

    6. C. Craeye, A fast impedance and pattern computation scheme for finite

    antenna arrays, IEEE Trans Antennas Propag 54 (2006), 30303034.

    7. E. Lucente, A. Monorchio, and R. Mittra, An iteration-free MoM

    approach based on excitation independent characteristic basis functions

    for solving large multiscale electromagnetic scattering problems, IEEE

    Antennas Propag 56 (2008), 9991007.

    8. C. Delgado, M.F. Catedra, and R. Mittra, Application of the character-

    istic basis function method utilizing a class of basis and testing func-tions defined on NURBS patches, IEEE Antennas Propag 56 (2008),

    784791.

    9. D.C. Ghiglia and L.A. Romero, Robust two-dimensional weighted and

    unweighted phase unwrapping that uses fast transforms and iterative

    methods, J Opt Soc Am A 11 (1994), 107117.

    2009 Wiley Periodicals, Inc.

    EXTRACTION OF EFFECTIVEPERMITTIVITY AND PERMEABILITY OFPERIODIC METAMATERIAL CELLS

    Dong Hyun Lee and Wee Sang ParkDivision of Electronic and Computer Engineering, Pohang Universityof Science and Technology, San 31, Hyojadong, Namgu, Pohang,Kyungbuk, Korea 790784; Corresponding author:[email protected]

    Received 14 November 2008

    ABSTRACT: The complex permittivity and permeability of various pe-

    riodic metamaterial (MTM) cells are extracted by using a fictitious rect-

    angular waveguide consisting of perfect electric conductor and perfect

    magnetic conductor walls. The shapes of the MTM cells include a thin

    wire (TW), a single split-ring resonator (SSRR), a double SRR (DSRR),

    a modified SRR, and a structure combining the TW and DSRR. The TW

    falls in the negative-/positive- region, the SRRs in the positive-/neg-

    ative- region, and the combined structure in the negative-/negative-region. We also investigate how the permittivity and permeability are

    affected by the dimension parameters of the MTM cells. Another extrac-

    tion technique utilizing time domain signals is developed to overcome

    some limitations that the waveguide technique cannot handle. 2009

    Wiley Periodicals, Inc. Microwave Opt Technol Lett 51: 18241830,

    2009; Published online in Wiley InterScience (www.interscience.wiley.

    com). DOI 10.1002/mop.24468

    Key words: material constant extraction; metamaterial; permittivity;

    permeability

    1. INTRODUCTION

    Metamaterials (MTMs) are broadly defined as artificial effectively

    homogeneous materials possessing unusual electromagnetic prop-

    erties not found in nature. Effectively, homogeneous materials are

    electromagnetic periodic structures whose periodicity is much

    smaller than the guide wavelength in the material. Usually, this

    periodicity is smaller than a quarter of the wavelength [1]. If the

    average unit cell size is four times smaller than the wavelength, a

    wave in the structure is dominantly characterized by refractive

    phenomena rather than scattering or diffraction phenomena. Under

    this condition, the electromagnetic properties of the structure are

    determined by well-defined material constants, permittivity (),

    and permeability (). Therefore, the extraction of the material

    constants of a given MTM is important for predicting the electro-

    magnetic properties of the MTM. In the last 5 years, several

    methods for extracting the parameters have been reported, and

    these methods used reflection and transmission coefficients of the

    structures [25]. The extraction methods were existed before the

    MTMs get interests [6]. The technique in [6] uses a coaxial

    transmission line and is designed for extracting the material con-

    stants of normal materials ( 0, 0). The basic extraction

    scheme of the technique is similar with other techniques [25], but

    the formulas in this technique include less ambiguity.

    Unfortunately, the formulas in reference [6] are valid only for

    low-loss materials, so they must be modified to apply them to

    highly lossy metamaterials. This article presents the modifiedequations and uses a fictitious waveguide simulation to obtain the

    reflection and transmission coefficients of unit cells. The material

    constants are extracted for representative MTMs, and their char-

    acteristics are analyzed. In addition, another extraction technique

    is developed to overcome some limitations of the previous tech-

    nique.

    2. EXTRACTION OF MATERIAL CONSTANTS FOR TEMWAVE INCIDENCE

    Figure 1 depicts the electric and magnetic field intensities and the

    wave propagation direction when a transverse electromagnetic

    (TEM) wave is normally incident on a three-region dielectric

    interface. Region 1 and Region 3 are free space (0, 0) and

    Region 2 is an unknown material (r, r) with thickness d. Thefields in each region can be written as follows:

    Region 1

    E1 xEi0ejk0z

    xEr0ejk0z

    H1 yEi0

    0ejk0z y

    Er0

    0ejk0z (1)

    where k0 000, 0 0/0.Region 2

    E2 xE2ejk2z xE2

    ejk2z

    H2 yE2

    2ejk2z y

    E2

    2ejk2z (2)

    Region 1:

    Air

    Region 2:

    Unknown material

    Region 3:

    Air

    x

    z

    iE

    iHik

    rE

    rH

    rk

    2E

    2H2

    3E

    3H

    ik

    2E

    2H

    2

    z = 0 z = d

    00 , 00 ,rr,

    Figure 1 Electromagnetic field distributions in multiple dielectric mate-

    rials for an incident TEM wave. TEM: transverse electromagnetic. E :

    electric field intensity (V/m), H : magnetic field intensity (A/m), k :

    wavenumber vector. [Color figure can be viewed in the online issue, which

    is available at www.interscience.wiley.com]

    1824 MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 51, No. 8, August 2009 DOI 10.1002/mop

  • 7/28/2019 24468_ftp

    2/7

    where k2 22 0r0r k0rr, 2 0r/0r 0r/r.

    Region 3

    E3 xE3ejk0z

    H3 yE3

    0ejk0z. (3)

    After applying boundary conditions for the electric and magnetic

    fields at two interfaces z 0 and z d, the reflection and

    transmission coefficients on the interfaces can be written as

    S11Er0

    Ei0

    1z2 12

    1 122 z2

    T S21E3

    Ei0

    1 122 z

    1 122 z2

    (4)

    where

    122 0

    2 01

    r/r 1r/r 1

    (5)

    z ej2d ej w/c0 rrd (6)

    Note that reference [6] assumed that Region 2 is not too lossy, so

    Eq. (6) was simplified as z 1. However, MTMs are highlylossy, and the equations were modified not to use this approxima-

    tion in this article.

    Equation (4) can be rewritten in a more compact form as

    S11122 12 1 S21

    2 S11

    2 S11 0 (7)

    z S11 S21 12

    1 S11 S21 12. (8)

    The solution of Eq. (7) is followed,

    12 1 S21

    2 S11

    2

    2S11 1 S21

    2 S11

    2

    2S11

    2

    1. (9)

    Note that one of the solutions must be chosen to satisfy the

    condition of 12 1. After 12 is determined by Eq. (9), r/rcan be obtained from Eq. (5):

    rr

    1 12

    1 12X. (10)

    rr can also be determined from Eqs. (6) and (8):

    rrjc0

    d lnz jarg z Y. (11)

    With the relations expressed in Eqs (10) and (11), the material

    constants can be expressed as

    r r,realjr,imagX/Y

    r r,realjr,imagX Y. (12)

    Therefore, the material constants r and r of the unknown mate-

    rial can be determined from S-parameters at the interfaces of the

    material.

    3. EXTRACTION OF MATERIAL CONSTANTS USING AFICTITIOUS WAVEGUIDE

    Tangential electric and magnetic fields become 0 on perfect elec-

    tric conductor (PEC) and perfect magnetic conductor (PMC) sur-

    faces, respectively. Therefore, a fictitious waveguide possessing

    the PEC surfaces on the top and bottom walls and the PMC

    surfaces on the front and back walls can support TEM waves. If

    unit cells of a periodic structure are symmetric, the field distribu-

    tions in the waveguide for one unit cell are identical with those of

    the infinite periodic structure. Therefore, the fictitious waveguide

    can produce an infinitely periodic effect. In addition, the simula-

    tion running time for obtaining the S-parameters can be dramati-

    cally reduced by simulating the waveguide with one unit cell rather

    than the infinite structure.

    Figure 2 shows the waveguide simulation scheme for acquiring

    S-parameters of an unknown material. The PEC boundary condi-

    tion is applied to the upper and bottom walls at the waveguide, and

    the PMC boundary condition is applied to the front and back walls.

    The length and width of the material with thickness Lz

    (d) are Lx

    and Ly, respectively. The S-parameters of the materials interface

    are obtained by shifting a phase reference plane from the ports tothe materials interfaces (Em). After the S-parameters are obtained,

    the material constants can be extracted by using the formulas in the

    previous section.

    The extraction technique was validated by extracting the ma-

    terial constants of a Lorentz material. A commercial full-wave

    simulator, CST MWS ver. 2006b, was used to simulate the struc-

    ture in Figure 2. The extracted r andr [Fig. 3(b)] agree well with

    the original values [Fig. 3(a)].

    4. MATERIAL CONSTANTS EXTRACTION OF MTM UNITCELLS

    4.1. Thin Wire (TW)

    A periodic thin wire (TW) is a well-known epsilon-negative

    (ENG) material because of a plasma phenomenon [5, 7]. Figure 4

    shows the TW unit cell with parameters and its extracted material

    constants. Note that the electric field vector must be parallel with

    the wire. The simulation results show that the TW has the nega-

    tive-/positive- region of an ENG material [Fig. 4(b)].

    E (x)

    k (z)

    H (y)

    Lx

    Ly

    Lz = d

    Em

    Embedding

    distance

    Unknown Material

    port 1

    port 2

    PEC

    PMC

    Figure 2 Simulation scheme for acquiring S-parameters using the ficti-

    tious waveguide. PEC: Perfect Electric Conductor. PMC: Perfect Magnetic

    Conductor. [Color figure can be viewed in the online issue, which is

    available at www.interscience.wiley.com]

    DOI 10.1002/mop MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 51, No. 8, August 2009 1825

  • 7/28/2019 24468_ftp

    3/7

    4.2. Split-Ring Resonator (SRR)

    A split-ring resonator (SRR) is a well-known mu-negative (MNG)

    material first reported by Pendry [8]. The studies of SRR have

    variously been carried out so far [912]. By Faradays law, when

    a time-varying magnetic field is applied to a conducting loop, a

    current is induced on the loop. According to Lenzs law, the

    magnetic moment from the induced current always opposes the

    applied field, thus reducing the magnetic flux density. Therefore,

    the induced magnetic moment, having opposite phase with the

    excited magnetic field, does not increase the permeability. In

    contrast, an SRR is capacitive at low frequency (closed loops are

    inductive), so an induced magnetic moment on the SRR becomesin-phase with the applied field. The SRR acts as magnetic dipoles

    and can have permeability greater than 1. A magnetic resonance

    occurs due to a circular current on the SRR, and the amount of this

    current is determined by the inductance of the SRR and capaci-

    tance of the gap. If the quality factor of this resonator is high

    enough, the permeability can be negative [9].

    4.2.1. Single SRR (SSRR). Figure 5 shows the structure of a

    single SRR (SSRR) unit cell and the extracted material constants,

    where f0 is the resonance frequency, f1 is the frequency exhibiting

    the lowest permeability, and f2 is the frequency at which the sign

    ofr,real changes to positive. In the designed SSRR, f0 is 8.7 GHz,

    where the electric size of the SSRR is around 0.15. SRRs are

    usually considered to be magneto-dielectrics exhibiting high r,realgreater than 1. The magneto-dielectric region of the SSRR appears

    from 7.5 to 8.6 GHz. Using this region can improve the input

    impedance bandwidth of antennas [13, 14]. The MNG region is

    observed from 8.6 to 10.2 GHz.

    Table 1 shows how changes in the individual dimensional

    parameters affected f0, f1, r,real at f1, and f2. When one parameter

    was increased, all others remained constant. For example, when R

    is increased 1 mm, f0, f1, r,real at f1, and f2 decreases 3.31, 3.30,

    1.34, and 2.62 GHz, respectively. When R is increased, the equiv-

    alent inductance of the SRR is increased and f0 is decreased. When

    gap is increased, the decreased equivalent capacitance causes to f0increase. When w is increased, the increased equivalent inductance

    causes f0 to increase. The r,real is dominantly affected by the

    parameters of gap and w because the amount of the induced current

    is mainly affected by the capacitance of structures gap.

    4.2.2. Double SRR (DSRR). The structure and the extracted ma-

    terial constants of a double SRR (DSRR) unit cell are presented in

    Figure 6. The depth, which is not depicted in Figure 6(a), is the

    thickness of the one SRR. Note that the external dimension is same

    Figure 3 Extracted material constants of Lorentz material. (a) Original

    r and r. (b) Extracted r and r. Lx 6, Ly 10, Lz 1 (mm). [Color

    figure can be viewed in the online issue, which is available at www.inter-

    science.wiley.com]

    w

    R

    aX

    aY

    aZ

    E

    H

    k

    (a)

    (b)

    Figure 4 (a) TW simulation geometry. (b) Extracted material constants.

    aX aY aZ 5, R 5, w 0.25 (mm). [Color figure can be viewed

    in the online issue, which is available at www.interscience.wiley.com]

    1826 MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 51, No. 8, August 2009 DOI 10.1002/mop

  • 7/28/2019 24468_ftp

    4/7

    as that of the SRR. A DSRR has a capacitance between the twoSRRs, so f0 is lower than that of the SRR. The electric size of

    DSRR at f0 is around 0.14. The magneto-dielectric region appears

    from 7.4 to 8.4 GHz, and the MNG region is observed from 8.4 to

    10 GHz.

    The effects of changing the dimension parameters of the DSRR

    are similar with those of the SSRR (Table 2). The r,real is mainly

    affected by parameter s1 (distance between SRRs).

    4.2.3. Modified SRR (MSRR). The structure and the extracted

    material constants of a modified SRR (MSRR) unit cell are pre-

    sented in Figure 7. The MSRR is symmetric, so that this structure

    is reported as a structure overcoming the bianisotropy problem

    which occurs with an asymmetric structure such as the SSRR and

    the DSRR [11, 12]. High capacitance occurs between the two

    SRRs, so f0 is lower than that of other SRRs. The electrical size at

    f0 is around 0.12. The magneto-dielectric region is from 5.5 to 6.9

    GHz, and MNG region appears from 6.9 to 8.8 GHz.

    The effects of changing the dimension parameters of the MSRR

    are similar with those of the other SRRs. The r,real is mainly

    affected by s2 (Table 3). We can design SSRR, DSRR, and MSRR

    or predict their characteristics by using Tables 1, 2, and 3.

    w

    RaY

    aX

    aZ

    depth

    gap

    E

    H

    k

    (a)

    (b)

    Figure 5 (a) SSRR structure. (b) Extracted material constants. aX

    aY aZ 5, R 4, depth gap w 0.25 mm: f0 8.63 GHz, f1

    9.1 GHz, r,real (at f1) 5.6, and f2 10.18 GHz. [Color figure can be

    viewed in the online issue, which is available at www.interscience.wiley.

    com]

    aY

    aX

    Ro

    Ri

    w

    gap

    s1

    Ri = Ro 2s1 2w

    (a)

    (b)

    Figure 6 (a) DSRR structure. (b) Extracted material constants. aX

    aY aZ 5, Ro 4, depth gap 0.25, s1 0.75, w 0.25 (mm):

    f0 8.42 GHz, f1 8.84 GHz, r,real (at f1) 5.9, and f2 9.9 GHz.

    [Color figure can be viewed in the online issue, which is available at

    www.interscience.wiley.com]

    TABLE 1 Characteristics of SSRR when the Dimensions Are

    Varied

    (mm) f0 (GHz) f1 (GHz) r,real at f1 f2 (GHz)

    Gap 0.20 0.53 0.61 1.03 0.62

    w 0.20 0.34 0.39 0.91 0.39

    Depth 0.25 0.46 0.39 0.35 0.19

    R 1.00 3.31 3.30 1.34 2.62

    aX 1.00 0.05 0.01 1.52 0.25

    aY 1.00 0.20 0.10 0.40 0.07

    aZ 1.00 0.17 0.08 0.84 0.12

    TABLE 2 Characteristics of DSRR when the Dimensions Are

    Varied

    (mm) f0 (GHz) f1 (GHz) r,real at f1 f2 (GHz)

    Gap 0.20 0.50 0.57 1.17 0.58

    w 0.20 0.40 0.45 1.23 0.45

    Depth 0.25 0.51 0.46 0.64 0.26

    s1 0.25 0.85 0.94 1.84 1.00

    Ro 1.00 3.44 3.41 1.22 2.72

    aX 1.00 0.05 0.01 1.60 0.24

    aY 1.00 0.17 0.09 0.48 0.09

    aZ 1.00 0.11 0.05 0.73 0.09

    DOI 10.1002/mop MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 51, No. 8, August 2009 1827

  • 7/28/2019 24468_ftp

    5/7

    4.3. Double Negative (DNG) MTM Composed of TW and DSRR

    Figure 8 shows a double-negative (DNG) unit cell consisting of a

    TW and DSRR. The DNG region where the material constants are

    both negative is observed from 6.3 to 7.3 GHz [Fig. 8(b)].

    5. EXTRACTION OF MATERIAL CONSTANTS USING TIMEDOMAIN SIGNALS

    5.1. S-parameters from Time Domain Signals

    The extraction technique in the previous section can quickly ex-

    tract the material constants of an infinite periodic structure. How-

    ever, this technique cannot simulate finite structures and cannot be

    verified by experiments because of the fictitious waveguide. In

    addition, it cannot construct perfect periodic structures because

    the PEC and PMC walls mirror the fields instead of transfer-

    ring those. So this technique can make the perfect periodic effect

    only for symmetric structures. Note that the extraction errors of the

    SSRR, DSRR, and MSRR unit cells are disregardable even though

    their structures are asymmetric. In this section, a new extraction

    technique that can handle the above problems of the waveguide

    technique is introduced and verified.

    Figure 9 shows a procedure of extracting the S-parameters of an

    unknown material from time domain signals. To obtain an input

    signal (a reference signal), first, a configuration without the mate-

    rial is simulated [Fig. 9(a)]. Here, the probe P1(t) measures only

    the input signal a(t). Next, the reflected and transmitted signals,

    b(t) and c(t), are obtained by simulating the configuration with the

    material [Fig. 9(b)]. The transmitted signal is measured by a probe

    inserted behind the material. Then, the signals P2(t) and P3(t) are

    a combined signal a(t) b(t) and the transmitted signal c(t),

    respectively. Note that P1(t), P2(t), and P3(t) are all time domain

    signals, so the reflected signal b(t) can be obtained by simply

    subtracting P1(t) from P2(t). Applying Fast Fourier Transform to

    gap

    w

    depthR

    aX

    aY

    aZs2

    (a)

    (b)

    Figure 7 (a) MSRR structure. (b) Extracted constitutive parameters.

    aX aY aZ 5, R 4, depth gap 0.25, s2 1, w 0.25 (mm):

    f0 6.92 GHz, f1 7.34 GHz, r,real (at f1) 6.11, and f2 8.76 GHz.

    [Color figure can be viewed in the online issue, which is available at

    www.interscience.wiley.com]

    w

    Rw

    aX

    aY

    aZ

    Rs

    gap

    depth

    s1

    s3

    (a)

    (b)

    Figure 8 (a) DNG unit cell composed of a TW and DSRR. (b) Extracted

    constitutive parameters. aX aY aZ 5, Rs 4, Rw 5, depth

    gap 0.25, s1 0.2, s3 1.75, w 0.25 (mm). [Color figure can be

    viewed in the online issue, which is available at www.interscience.wiley.

    com]

    TABLE 3 Characteristics of MSRR when the Dimensions Are

    Varied

    (mm) f0 (GHz) f1 (GHz) r,real at f1 f2 (GHz)

    Gap 0.20 0.39 0.46 1.04 0.45

    w 0.20 0.18 0.17 0.68 0.10

    Depth 0.25 0.68 0.69 0.19 0.50

    s2 0.35 1.66 1.85 7.40 2.23

    R 1.00 2.82 2.75 1.64 1.85

    aX 1.00 0.04 1.59 0.33

    aY 1.00 0.11 0.01 0.78 0.22

    aZ 1.00 0.20 0.13 1.30 0.14

    1828 MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 51, No. 8, August 2009 DOI 10.1002/mop

  • 7/28/2019 24468_ftp

    6/7

    a(t), b(t), and c(t), gives the frequency domain signals A(f), B(f),

    and C(f). Then, S11 B(f)/A(f) and S21

    C(f)/A(f) at the location

    of probes can be calculated. Finally, S-parameters on the material

    interfaces are obtained by shifting the phase reference plane (Em).

    Then, the material constants can be extracted by using the equa-

    tions in section 2.

    5.2. Verification of the Code

    This extraction technique was verified by the extractions for two

    materials, a material ofr 10 j0, r 1 j0, and the Lorentz

    material in Figure 3(a). The results are plotted in Figure 10. The

    extracted material constants show good agreement with the orig-inal material constants.

    The developed extraction method requires longer simulation

    times because the structure is simulated twice, but it can construct

    the perfect periodic structure and can be applied on finite struc-

    tures.

    6. CONCLUSIONS AND FUTURE WORKS

    The material constant extraction technique using the fictitious

    waveguide was reviewed, and the material constants were ex-

    tracted for representative MTM unit cells: TW, SSRR, DSRR,

    MSRR, and combined TW/DSRR. The TW had a negative-/

    positive- region, the SRRs a positive-/negative- region, and

    the combined structure a negative-/negative- region. We also

    investigated how the material constants are affected by the dimen-

    sion parameters of the SSRR, DSRR, and MSRR cells, which

    should be helpful for designing SRRs and predicting their charac-

    teristics. To overcome some limitations of the waveguide tech-

    nique, an extraction technique utilizing time domain signals was

    developed and verified by extracting the material constants of

    known materials.

    Here, the material constants of infinite structures are consid-

    ered. By using the technique of using time domain signals, we are

    investigating the material constants of finite structures and cou-

    pling effects between different layers.

    REFERENCES

    1. C. Caloz and T. Itoh, Electromagnetic metamaterials, Wiley-Inter-science, Hoboken, NJ 2006.

    2. D.R. Smith, S. Schulz, P. Markos, and C.M. Soukoulis, Determination

    of effective permittivity and permeability of metamaterials from re-

    flection and transmission coefficients, Phys Rev E 71 (2002), 195104.

    3. Richard W. Ziolkowski, Design, fabrication, and testing of double

    negative metamaterials, IEEE Trans Antennas Propag 51 (2003),

    15161529.

    4. D.R. Smith, D.C. Vier, TH. Koschny, and C.M. Soukoulis, Electro-

    magnetic parameter retrieval from inhomogeneous metamaterials,

    Phys Rev E 71 (2005), 036617.

    5. P.D. Imhof, Metamaterial-based Epsilon negative [ENG] Media: Anal-

    ysis & designs, Master Thesis, Electric and Computer Engineering at

    the University of Arizona, February 2006.

    6. C.C. Courtney, Time-domain measurement of the electromagnetic

    properties of materials, IEEE Trans Microwave Theory Tech 46

    (1998), 517522.

    7. J.B. Pendry, A.J. Holden, D.J. Robbins, and W.J. Stewart, Low fre-

    quency plasmons in thin-wire structure, J Phys Condens Lett 10

    (1998), 47854809.

    8. J.B. Pendry, A.J. Holden, D.J. Robbins, and W.J. Stewart, Magnetism

    from conductors and enhanced nonlinear phenomena, IEEE Trans

    Microwave Theory Tech 47 (1999), 20752084.

    9. S. Maslovski, P. Ikonen, I. Kolmakov, and S. Tretyakov, Artificial

    magnetic materials based on the new magnetic particle: Matasolenoid,

    Prog Electromagn Res 54 (2005), 6181.

    10. M. Kafesaki, TH Koschny, R.S. Penciu, T.F. Jundogdu, E.N. Econo-

    mou, and C.M. Soukoulis, Left-handed metamaterials: Detailed nu-

    Source

    Air

    P1(t)

    a(t)

    Em

    probe

    (a)

    Source

    Material

    Em

    P2(t)

    b(t)a(t) c(t)

    P3(t)

    Em

    (b)

    Figure 9 S-parameters acquisition scheme from time domain signals. (a)

    Without material. (b) With material. [Color figure can be viewed in the

    online issue, which is available at www.interscience.wiley.com]

    Figure 10 Extracted material constants. (a) Material ofr 10 j0,

    r 1 j0. (b) Lorentz material in Figure 3(a). [Color figure can be

    viewed in the online issue, which is available at www.interscience.wiley.

    com]

    DOI 10.1002/mop MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 51, No. 8, August 2009 1829

  • 7/28/2019 24468_ftp

    7/7

    merical studies of the transmission properties, J Opt A 7 (2005),

    S12S22.

    11. R. Marques, F. Medina, and R. Rafii-El-Edrissi, Role of bianisotropy

    in negative permeability and left-handed metamaterials, Phys Rev B

    65 (2002), 1440/11440/6.

    12. I. Bulu, H. Caglayan, and E. Ozbay, Experimental demonstration of

    labyrinth-based left-handed metamaterials, Opt Exp 13 (2005).

    13. P.M.T. Ikonen, S.I. Maslovski, C.R. Simovski, and S.A. Tretyakov,

    On artificial magnetodielectric loading for improving the impedance

    bandwidth properties of microstrip antennas, IEEE Trans Antennas

    Propag 54 (2006), 16541662.

    14. H. Mosallaei and K. Sarabandi, Design and modeling of patch antenna

    printed on magentor-dielectric embedded-circuit metasubstrate, IEEETrans Antennas Propag 55 (2007), 4552.

    2009 Wiley Periodicals, Inc.

    DESIGN OF AN INTEGRATED LOOPCOUPLER AND LOOP ANTENNA FORRFID APPLICATIONS

    Randy BancroftRandwulf Technologies, 2837 Perry Street, Denver, CO 80212;Corresponding author: [email protected]

    Received 14 November 2008

    ABSTRACT: The design of a novel planar loop antenna with an inte-

    grated coupling loop for RFID applications is described. The antenna con-

    sists of a center nonradiating coupling loop, which acts as a phase shifter

    to produce a radiating current distribution on the outer loop. 2009

    Wiley Periodicals, Inc. Microwave Opt Technol Lett 51: 1830 1833, 2009;

    Published online in Wiley InterScience (www.interscience.wiley.com).

    DOI 10.1002/mop.24467

    Key words: microstrip antenna; RFID; loop antenna; planar antenna;

    RFID sensor

    1. INTRODUCTION

    Currently existing sensor designs for RFID tags at 915 MHz are

    generally designed to operate as either a nonradiating coupler to

    drive an electrically small (low radiation) tag or as a radiating

    antenna, which is optimized to couple with an electrically large

    RFID tag. In many applications, it is a requirement that both

    electrically small and electrically large RFID tags be addressable.

    This article introduces the design of an antenna/coupler hybrid,

    which is optimized to work well with both electrically small and

    electrically large RFID tags.

    2. NEAR FIELD AND FAR FIELD RFID TAGS

    RFID tags exist in many forms. Electrically small tags are often

    small loops. Some are circular and others are oval shaped as shown

    in Figure 1. The average radius of the tag of Figure 1 (a ) is 6 mm.

    When the loop is electrically small, the radiation resistance Rr and

    ohmic resistance RCu of the loop may be estimated using Eqs. (1)

    and (2) [1].

    Rr 31,200a2

    2 (1)

    RCu 2a

    Rs (2)

    Rs 02f2 (3)where:

    is the free space wavelength

    0 is the permeability of free space

    f is the frequency in Hz and

    Rs is the surface resistance.

    For the tag of Figure 1, the radiation resistance is only Rr

    34.53 m and the copper resistance is RCu 364.23 m at 915MHz. The estimated radiation efficiency of this loop is 9.48%.

    HFSS analysis indicates the radiation efficiency of this tag is only

    0.53%. These small values indicate that a tag of this size must rely

    on direct current inducement on the loop to provide an adequate

    modulation of the driving point reflection coefficient for RFID

    applications.

    The backscattering from an electrically small RFID tag is

    further attenuated by the reactive part of the driving point imped-

    ance [2]. The desire to use both near field (direct inductive

    coupling) and far field (radiative coupling) tags in industrial

    applications leads to the need for an antenna which had both a

    coupling region for nonradiative tags and satisfactory radiation for

    radiative tags.

    3. LOOP ANTENNA THEORY

    The resonant and antiresonant dimensions of circular loop anten-

    nas have been discussed by Schekunoff and Friis [3].

    C n n 1,2,3, . . . (4)

    C 2n 1

    2 n 0,2,4, . . . (5)

    Equation (4) gives the resonant circumferences (C) of a loop

    antenna in terms of wavelength and Eq. (5) relates the antireso-

    nant circumferences.

    When the circumference of a loop antenna is less than /2,neither resonance nor antiresonance exists. When a loop with a

    circumference which is less than one-half wavelength is fed in a

    balanced manner, the current is approximately uniform around the

    perimeter of the loop. The current on one side of the loop will

    cancel with a current on the opposite side of the loop in the far

    field. When the small loop is fed in an unbalanced manner, the two

    Figure 1 HFSS model of a common electrically small RFID loop tag.

    The outer length is 17.27 mm and the outer width is 8.23 mm with a

    conductor width of 1 mm. [Color figure can be viewed in the online issue,

    which is available at www.interscience.wiley.com]

    1830 MICROWAVE AND OPTICAL TECHNOLOGY LETTERS / Vol. 51, No. 8, August 2009 DOI 10.1002/mop