A Compact Equivalent Circuit Model for the SRR ... - IEEE Xplore

67 downloads 0 Views 362KB Size Report
Dept of Electrical & Computer Engineering, National University of Singapore, Kent Ridge, Singapore 119260. Abstract-This paper proposes a compact ...
A Compact Equivalent Circuit Model for the SRR Structure in Metamaterials M. F. Wu, F. Y. Meng, and Q. Wu 1.

School of Electronics and Information Technology, Harbin Institute of Technology, Harbin, Heilongjiang, China 150001

J. Wu 2.

China Research Institute of Radio wave Propagation, Xinxiang, Henan, China 102206

L. W. Li 3.

Dept of Electrical & Computer Engineering, National University of Singapore, Kent Ridge, Singapore 119260

Abstract-This paper proposes a compact equivalent circuit

macro performance of SRR arrays or counterparts. While

model to characterize the split ring resonator (SRR) structure in

accurate,

these

field

solvers

are

time

consuming,

metamaterials, based on the modeling approaches of the spiral

computationally intensive and require experience. These field

inductors in MMIC. In the equivalent circuit model, the lumped

solvers also do not provide any insight into the engineering

elements distribute compactly, and can be directly estimated by

trade-offs involved in the design of new geometry. Although

the geometry parameters of the SRR. Numerical simulation

the better field solvers are excellent for verification, they are

results show that the model can simultaneously predict more than

inconvenient at the initial design and optimization stages.

one resonance frequency of the SRR, and the predicted resonance

To meet the requirement, this paper proposes an approach to

frequency is accurate enough. The equivalent circuit model can be

analyze an individual SRR by a compact equivalent circuit

advantageous to probe the root of the negative permeability,

model based on the spiral inductors modeling methods in

especially

SRR-based

MMIC. The equivalent circuit model can facilitate the analysis

metamaterials exhibiting left-handed properties over multiple

and optimization of SRR-based metamaterials, and predict

frequency bands.

more than one resonant frequency, which is very important for

helpful

for

engineers

to

analyze

SRR-based metamaterials exhibiting left-handed properties I.

INTRODUCTION

In 1999, the negative value of effective permeability was

over multiple frequency bands [3]. II.

THE COMPACT EQUIVALENT CIRCUIT MODEL OF SRR

first achieved by using an array of the metallic artificial atoms called Split Rings Resonators (SRRs) [1]. Since then, many

Researches have been shown that an individual SRR can

research works have been conducted after various methods for

realize the electromagnetic performances of SRRs arrays [4],

realizing negative permeability, especially after the negative

so we had to just look into the performances of an individual

refractive index was experimentally verified [2]. SRR has been

SRR, instead of investigating the entire effective behaviors of

approved to be a useful artificial atom for design and

SRRs arrays in this section. Fig.1 shows the geometry of an

fabrication the artificial media. However, there are still some

individual square SRR, which functions as current sheets while

performances to be improved, such as the anisotropy, lossy and

an incident plane wave magnetic field H is penetrating through

bandwidth with the SRR-based metamaterials.

it. Differing from the spiral inductors, there are two splits in the

It is promising to arrive at the aims through probe the micro

out ring and the inner ring of SRR. So the SRR equivalent

electromagnetic effects of SRR or counterpart, and design or

inductors approximately estimated by four side parts inductors,

improve its geometry. However, most of this modeling work

namely L1, L2, L2 and L3, respectively shown in Fig.2,

has centered on the development of field solvers to predict the

according to the square spiral inductors lumped circuit models.

0-7803-9433-X/05/$20.00 ©2005 IEEE.

APMC2005 Proceedings

L2 produced by one side integrated outer and inner current sheet, while L1 and L3 represent the self-inductors produced by an

LD =

K µ0 n 2 L 2 [ln( ) + 0.5 + 0.178 ρ + 0.0146 ρ 2 + ρ 2π 0.5(n − 1) S 2 (n − 1) S 1 W + t )] + 0.178 × − ln( n n W ( ρ n) 2

integrated sheet and a sheet with a split.

K=

(2 L − 2 S ) − D (2 L − 2 S )

(3)

(4)

Hence the equivalent circuit model and its lumped elements distribution can be shown in Fig.3, and several distributed capacitances also exist in the equivalent circuit model, namely C1, C2, C2, and C3 respectively, which can be determined by the following approaches.

Fig. 1. The geometry configuration of an individual square SRR

Fig. 3. The lumped element distribution in equivalent circuit model of SRR Fig. 2. The SRR self-inductors calculation divisions in four sides

The total capacitances with SRR consist of two parts. One

The self-inductors produced by one side integrated outer and

part is the coupling capacitance between the outer and the inner

inner current sheet such as L2 can be directly estimated by

rings, and the other one is produced the electric charges accumulate at the two splits. The coupling capacitance Cc can

followed formulas [5]: LD =

µ0 n L 2 [ln( ) + 0.5 + 0.178 ρ + 0.0146 ρ 2 ρ 2π

be estimated by the following equation [6] [7]:

2

0.5(n − 1) S 2 (n − 1) S 1 W + t )] + + 0.178 − ln( 2 n n W ( ρ n)

ρ=

nW + (n − 1) S L

(1)

Cc = [0.06 + 3.5 × 10−5 (rout + rin )]

(5)

Then the coupling capacitance should be divided into four (2)

Here, t represents the vertical parameter of SRR, and the lateral

equal counterparts named C0 for SRR four sides, hence: 1 C0 = [0.06 + 3.5 × 10−5 (rout + rin )] 4

(6)

parameters are completely specified by the following:

Here rout and rin represent the radius of the outer circumcircle

1. 2. 3. 4.

and the inner circumcircle of SRR respectively.

The number of turns, n The metal width, W

It is vital to estimate the capacitances at the outer and inner

The length of out ring sides, L

splits, but it is difficult because of the intense electromagnetic

The space between inner and outer rings, S

brink effects, so they are estimated by corrected the parallel

However, L1 and L3 cannot be directly estimated by equations (1) and (2), because of the errors caused by splits. So a new variable named integrality coefficient K should be inserted to rectify the self-inductors of one side, thus L1 and L3 can be estimated by following expressions:

plane capacitance formula as the following approaches: Cinner split = 3ε 0 s/D

(7a)

Couter split = 25ε 0 s/D

(7b)

Here s means the effective parallel plane area, and D means the distance of the parallel planes at the splits.

The distributed capacitances C1, C2 and C3 can be determined by equations (6) and (7) as following approaches: C1 = C0 + Couter split C2 = C0 C3 = C0 + Cinner split

(8a) (8b) (8c)

which display the Z21 coefficient phase degree abruptly changes 180 degrees at above three frequency points. All above simulation results indicate that the equivalent circuit becomes resonant at 6GHz, 15.5GHz and 26.5GHz respectively, which are good agreement with numerical results with error less than 8%. Conclusions can be made based on above results that SRR has three resonant frequency 6GHz, 16.8GHz and 27.8GHz

III.

NUMERICAL SIMULATION AND DISCUSSION

respectively, and can be precision-predicted by the equivalent circuit model.

In order to testify the proposed equivalent circuit model, numerical simulations are necessary. Let’s choose a set of

IV.

CONCLUSION

geometry parameters of SRR shown in Fig.1 such as: L=6mm, W=0.5mm, S=1mm, D=1mm, and the vertical size t=0.5mm.

Based on the spiral inductors modeling method in MMIC,

Firstly, CST Microwave Studio is considered to predict the

the proposed equivalent circuit model has compact distributed

resonance frequencies of SRR. Shown in Fig.4, a vacuum

elements and they can be directly estimated by the geometry

waveguide is employed, then its transmission coefficient is

parameters. Numerical simulation results testify that the

nearly 1 and its phase of image part of the Z21 is shown in Fig.5. equivalent circuit model can simultaneously predict three The results indicate that electromagnetic wave Propagate

resonant frequency points with high precision, which is

through the vacuum waveguide with no reflection, and its Z21

extremely important for SRR-based metamaterials exhibiting

coefficient phase degree abruptly changes 180 degrees at about

left-handed properties over multiple frequency bands.

20.5GHz. Then SRR can be centered in the vacuum waveguide shown

ACKNOWLEDGEMENT

in Fig.4, thus the transmission coefficient of the waveguide with SRR is simulated shown in Fig.6, which display S21

This work was partly supported by a grant from the National

abruptly drops to less than –35 dB at 6GHz, 16.8GHz and

Natural Science Foundation of China (No. 60372051), and

27.8GHz respectively, and the Z21 coefficient phase degree

Fund for the National Key Laboratory of Electromagnetic

abruptly changes 180 degrees at above three frequency points

Environment (No.514860303).

besides 20.5GHz, shown in Fig.7. The simulation results indicate that the vacuum waveguide become almost reflective

REFERENCES

at three frequency points because of the SRR centered in it, thus we can estimate that the SRR has its resonance at 6.0 GHz, [1]

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

16.8 GHz and 27.8 GHz respectively.

“Magnetism from conductors and enhanced nonlinear phenomena,”

Finally, we will employ the equivalent circuit model to

IEEE Transactions on Microwave Theory and Techniques, vol. 47, no.

predict the resonance frequencies of the SRR. According to its geometry parameters L, W, S, t and D, we can apply the

11, pp. 2075-2084, November 1999.

[2]

formulas (1)-(8) to determined the distributed elements L1, L2, L3, C1, C2 and C3。Then the transmission coefficient of the

R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science, vol. 292, pp. 77-79, April 2001.

[3]

Hongsheng Chen, Lixin Ran, Jiangtao Huangfu et.al. Metamaterial

equivalent circuit can be obtained and shown in Fig.8, which

exhibiting left-handed properties over multiple frequency bands. J. Appl.

display S21 abruptly drops to less than –35 dB at 6.0 GHz, 15.5

Phys.96, 5338 (2004).

GHz and 26.5 GHz respectively. Furthermore, the phase of its Z21 coefficient can be simultaneity obtained and shown in Fig.9,

[4]

R. W. Ziolkowski, “Design, fabrication, and testing of double negative metamaterials,” IEEE Transactions on Antennas and Propagation, vol.

51, no. 7, pp. 1516-1529, July 2003.

[5]

[6]

[7]

Qun Wu, Ming-feng Wu, Fan-yi Meng, Jian Wu, and Le-wei Li.

Sunderarajan S. Mohan, The Design, Modeling and Optimization of

Modeling the Effects of an Individual SRR by Equivalent Circuit

On-Chip Inductor and Transformer Circuits, Ph.D. Thesis, Stanford

Method, 2005 IEEE AP-S International Symposium and USNC / URSI

University, 1999, 38-51.

National Radio Science Meeting, July 3-8, 2005, Washington, DC, USA

Qun Wu, Ming-feng Wu, Fan-yi Meng, Jian Wu, Le-Wei Li, “Research on The SRR Structure LHM Based on Transmission Line Theory”. Journal of Radio Science, in Chinese, in press.

Fig. 4. The geometry of vacuum waveguide with

Fig. 5. The phase degree of the Z21 of the vacuum waveguide without SRR

centered SRR

Fig. 6. The transmission coefficient S21 of the vacuum waveguide with centered SRR.

Fig. 7. The phase degree of the Z21 of the vacuum waveguide with centered SRR 200

0

)) 1 ni 150 Z( g a 100 m (i e s a h p 50

-20

)) 1, 2( S ( -40 B d -60

0

-80

0

0

2

4

6

8

10 12 14 16 18 20 22 24 26 28 30

freq, GHz

2

4

6

8

10 12 14 16 18 20 22 24 26 28 30

freq, GHz

Fig. 8. The transmission coefficient S21 of the vacuum

Fig. 9. The phase degree of the Z21 of the vacuum

waveguide with centered SRR predicted by the

waveguide with centered SRR predicted by the

equivalent circuit model.

equivalent circuit model