NisrUNITED STATES

DEPARTMENT OF COMMERCE

TECHNOLOGY ADMINISTRATION

NATIONAL INSTITUTE OF STANDARDS

AND TECHNOLOGY

"idihiilfj. "^ ^"''° * '^CH Rl-c.

A111Q3 777141414

NIST

PUBLICATIONS

l\l/ST Technical Note 1355

OC iCO

/555

Transmission/Reflection and Short-Circuit

Line Methods for Measuring Permittivity

and Permeability

James Baker-Jarvis

Michael D. Janezic

John H. Grosvenor, Jr.

Richard G. Geyer

Electromagnetic Fields Division

Electronics and Electrical Engineering Laboratory

National Institute of Standards and Technology

Boulder, Colorado 80303-3328

May 1992

U.S. DEPARTMENT OF COMMERCE, Barbara Hackman Franklin, Secretary

TECHNOLOGY ADMINISTRATION, Robert M. White, Under Secretary for Technology

NATIONAL INSTITUTE OF STANDARDS AND TECHNOLOGY, John W. Lyons, Director
National Institute of Standards and Technology Technical Note

Natl. Inst. Stand. Technol., Tech. Note 1355, 120 pages (May 1992)

CODEN:NTNOEF

U.S. GOVERNMENT PRINTING OFFICE

WASHINGTON: 1992

For sale by the Superintendent of Docunnents, U.S. Government Printing Office, Washington, DC 20402-9325
Contents

1 Introduction 1

2 Theory for Coaxial Line and Rectangular Waveguide Measurements of

Permittivity and Permeability 6

2.1 Theory 6

2.1.1 Decomposition into TE, TM, and TEM Modes 6

2.1.2 Imperfect Sample Geometry 9

2.1.3 Perfect Sample in Waveguide 13

2.2 Permeability and Permittivity Calculation 15

2.2.1 Nicolson- Ross- Weir Solutions (NRW) 15

2.2.2 2-Port Solution Where Position is Determined Solely by Lairime and L 17

2.2.3 Two Samples of Different Length 18

2.3 Measurement Results 19

2.3.1 Measurements without Gap Corrections 20

2.3.2 Effects of Gaps between Sample and Waveguide 20

2.4 Permeameter 28

2.5 Uncertainty of Combined Permittivity and Permeability Determination . . 29

2.5.1 One Sample at One Position 30

2.5.2 Two Samples of Differing Lengths 43

2.6 Uncertainty in Gap Correction 49

2.6.1 Dielectric Materials 49

2.6.2 Magnetic Materials 50

2.6.3 Higher Order Modes 50

3 Optimized Solution 52

3.1 Introduction 52

3.2 Model for Permeability and Permittivity 54

3.2.1 Relaxation Phenomena in the Complex Plane 54

3.3 Numerical Technique 57

3.3.1 Algorithm 57

ni
3.3.2 Numerical Results 59

3.4 Permittivity and Permeability 61

3.4.1 Measurem.ents 61

3.4.2 Robustness of the Procedure 61

3.5 Discussion 65

4 Short-Circuit Line Methods 66

4.1 Theory 66

4.1.1 Two Samples of Different Lengths 68

4.1.2 Single Sample at Two Short-Circuit Positions 70

4.2 Measurements 71

4.3 Uncertainty of Short-Circuit Line Measurements 73

5 Discussion 79

6 References 83

7 Appendices 86

A Magnetism in Matter 87

A.l Description of Magnetic Phenomena 87

A. 1.1 Field Description of Electromagnetic Phenomena 87

A. 1.2 Types of Magnetism 89

A. 2 Paramagnetism 89

A. 2.1 Diamagnetism 90

A. 2. 2 Ferromagnetism 90

A. 2. 3 Ferrites and Antiferromagnetism 91

A. 3 Equations of Motion for the Magnetization Vector . 92

A. 3.1 The Torque Equation 92

A.3.2 Magnetized Magnetic Material: The Polder Matrix 93

B Fields in Waveguides 96

B.l Summary of Maxwell's Equations 96

B.2 Modes 97

B.2.1 TE Modes 97

B.2. 2 TM Modes 97

B.2. 3 TEM Modes 98

C Gap Correction 99

C.l Frequency- Dependent Gap Correction 99

C.1.1 Waveguide 99

IV
C.1.2 Coaxial Line 100

C.2 Frequency- Independent Approaches 101

C.2.1 Coaxial Capacitor Model for Dielectric Materials 101

C.2. 2 Rectangular Waveguide Model 103

C.3 Gap Correction for Magnetic Materials 104

C.3.1 Coaxial Line 104

C.3. 2 Waveguide 105

C.4 Gap Correction Formulas Derived Directly From Maxwell's Equations . . . 107

C.5 Mitigation of Air Gap Effects 108

D Causal Functions and Linear Response 110

D.l Introduction 110

D.2 Transfer Functions Ill

D.3 Kramers-Kronig Relations Ill
Abstract

The transmission/reflection and short-circuit Hne methods for measuring complex permittivity and permeabihty of materials in waveguides and coaxial lines are examined. Equations for complex permittivity and permeability are developed from first principles. In addition, new formulations for the determination of complex permittivity and permeability independent of reference plane position are derived. For the one-sample transmission/reflection method and two-position short-circuit line measurements, the solutions are unstable at frequencies corresponding to integral multiples of one-half wavelength in the

sample. For two-sample methods the solutions are unstable for frequencies where both

samples resonate simultaneously. Criteria are given for sample lengths to maintain stability. An optimized solution is also presented for the scattering parameters. This solution

is stable over all frequencies and is capable of reducing scattering parameter data on materials with higher dielectric constant. An uncertainty analysis for the various techniques

is developed and the results are compared. The errors incurred due to the uncertainty in

scattering parameters, length measurement, and reference plane position are used as inputs

to the uncertainty models.

Key words: Calibration; coaxial line; dielectric constant; loss factor; magnetic materials;

microwave measurements; permeameter; permeability measurement; permittivity measure- ment; reflection method; short-circuit; transmission; uncertainty; waveguide.
Chapter 1

Introduction

The goal of this report is to review and critically evaluate various transmission line measure- ment algorithms for combined permeability and permittivity determination and to present

results and uncertainty analysis for the techniques.

There is continual demand to measure accurately the magnetic and dielectric properties

of solid materials. Over the years there has been an abundance of methods developed for

measuring permeability and permittivity. Almost all possible perturbations or variations

of existing methods have been proposed for measurements. These techniques include freespace methods, open-ended coaxial probe techniques, cavity resonators, full-body resonance

techniques, and transmission-line techniques. Each method has its range of applicability

and its own inherent limitations. For example, techniques based on cavities are accurate,

but not broadband. Nondestructive techniques, although not most accurate, allow the

maintenance of material integrity. Transmission line techniques are the simplest of the

relatively accurate ways of measuring permeability and permittivity of materials. Transmission line measurements usually are made in waveguide or coaxial lines. Measurements

are made in other types of transmission lines for special applications, but for precise measurements, rectangular waveguides and coaxial hues are usually used. The three major

problems encountered in transmission line measurements are air gaps, half- wavelength res- onances, and overmoding.

Coaxial lines are broadband in the TEM mode and therefore are attractive for permittivity and permeability measurements. The problem with coaxial lines, however, is that

due to the discontinuity of the radial electric field, any air gap around the center conductor

degrades the measurement by introducing a large measurement uncertainty. Belhadj-Tahar

et al. [1] have attempted to circumvent these difficulties with the development of a technique for a plug of material at the end of a coaxial line. In Belhadj-Tahar's approach there

is no center conductor hole. However, higher modes are excited at the transition between

the plug and the center conductor which complicates the analysis. Due to the complexity

1
of the method it is not apparent at this time whether this approach will replace the more

traditional single-mode models.

Transmission line techniques generally fall into the following categories:

• Off-resonance waveguide and coaxial line, full scattering parameter, 2-port measure- ments .

• Off-resonance short-circuit line, 1-port measurements.

• Open-circuit techniques.

• Resonant transmission-hne techniques.

The topic of this report will be the first two categories. We will also examine direct in- ductance measurement, which uses permeameter techniques. The off-resonance techniques can be broadly grouped into two categories:

• Point-by-point or uncorrelated-point techniques.

• Multi-point or correlated-point techniques.

The point-by-point technique is at present the most widely used reduction technique and

consists of solving the relevant scattering equations at single points. Multi-point techniques

consist of solving the nonlinear scattering equations using nonlinear least square algorithms.

Due to their relative simplicity, the off-resonance waveguide and coaxial hne transmission/reflection (TR) and short-circuit line (SCL) methods are presently widely used broadband measurement techniques. In these methods a precisely machined sample is placed in a section of waveguide or coaxial line and the scattering parameters are measured, preferably by an automatic network analyzer (ANA). The relevant scattering equations relate

the measured scattering parameters to the permittivity and permeability of the material.

One limitation of these techniques is that they require cutting of the sample and therefore

these techniques do not fall under the general category of nondestructive testing methods.

Another limitation is that these techniques require a small sample and therefore the res- onance characteristics of large sheets of the material are not studied. Network analyzers

have improved over the last years to a point where broad frequency coverage and accurate

measurement of scattering parameters are possible. This broadband capability unearths

another hmitation of present algorithms, that is, the instability of the measurement in the

vicinity of resonant frequencies.

In this report we assume that the materials under test are isotropic, homogeneous, and

in a demagnetized state. The solutions obtained in this report are both single-frequency

techniques and multiple frequency techniques. For the TR measurement, the system of

equations contains as variables the complex permittivity and permeability, the two reference
plane positions, and, in some applications, the sample length. In the T/R procedure we

have more data at our disposal than in SCL measurements, since we have all four of

the scattering parameters. In SCL measurements the variables are complex permittivity

and permeability, sample length, distance from sample to short-circuit termination, and

reference plane positions. However, in most problems we know the sample length, reference

plane position, and distance from the reflector to the sample. In these cases we have four

unknown quantities (complex permittivity and permeability) and therefore require four

independent real equations to solve for these variables. These equations can be generated by

taking reflection coefficient data at two positions in the transmission line, thus yielding the

equivalent of four real equations for the four unknown quantities. A problem encountered

in measurements is the transformation of S-parameter measurements at the calibration

reference planes to the air-sample interface. This transformation requires knowledge of the

position of the sample in the sample holder. Information on reference plane position is limited in many applications. The port extension and gating features of network analyzers

are of some help in determining reference plane position, but do not completely solve the

problem. Equations that are independent of reference plane position are desirable.

Most of the present transmission-line techniques [2,3,4], with some variations, are based on the procedure developed by Nicolson and Ross [5] and Weir [6] for obtaining 2-port,

off-resonance, broadband measurements of permeability and permittivity. In the NicolsonRoss-Weir (NRW) procedure the equations for the scattering parameters are combined in

such a fashion that the system of equations can be decoupled. This procedure yields an

explicit expression for the permittivity and permeability as a function of the S-parameters.

These equations are not well-behaved for low-loss materials at frequencies corresponding

to integral multiples of one-half wavelength in the sample. In fact, the NRW equations are

divergent, due to large phase uncertainties for very low-loss materials at integral multiples

of one-half wavelength in the material. Many researchers avoid this problem by measuring samples which are less than one-half wavelength long at the highest measurement

frequency. The advantage of the NRW approach is that it yields both permittivity and permeability over a large frequency band. As a special case of the NRW equations, Stuchly and

Matuszewski [7] found solutions to the scattering equations for nonmagnetic materials and

derived two explicit equations for the permittivity. Delecki and Stuchly [8] have studied the

uncertainty analysis for infinitely long samples using the bilinear and Schwarz-Christoffel

transformations. Franceschetti [9] was one of the first to perform a detailed uncertainty

analysis for TR measurements. Ligthart [10] developed an analytical method for permittivity measurements at microwave frequencies using an averaging procedure. In Ligthardt's

study, a single-moded cylindrical waveguide was filled with a homogeneous dielectric with

a moving short-circuit termination positioned beyond the sample. This study focused primarily on single-frequency measurements rather than on broadband measurements.

The short-circuit fine (SCL) method was introduced by Roberts and von Hippel [11] over fifty years ago as an accurate broadband measurement procedure. The SCL measure-
ment method uses data obtained from a short-circuit 1-port measurement to calculate the

dielectric and magnetic properties. SCL is useful when 2-port measurements are not possible, for example, in high temperature measurements [12] and remote sensing applications.

When an ANA is used, the sample is positioned in either a waveguide or coaxial hne and the

reflection coefficient is measured. The determination of the permittivity and permeability

usually proceeds by solving a transcendental equation that involves the sample length, sample position, and reflection coefficient. With modern computer systems, iterative solutions

of the resulting transcendental equations are easy to implement. However, they require an

initial guess. The resultant nonlinear equations have an infinite number of solutions due

to periodic functions. The physical solution can be determined by group delay arguments

or by measuring two samples with differing lengths. Much of the theory developed for the

SCL technique was developed for use with a slotted line. Present-day network analyzers

usually measure scattering parameters. Therefore in this report we derive equations from

a scattering approach.

The SCL method has endured over the years, and as a result there is an extensive

literature. In this report we attempt to review only the most relevant work on the subject.

Short-circuit line methods can be broadly separated into two- position techniques and twosample techniques. In the two-position technique 1-port scattering parameters are measured

for a sample in two different positions in the sample holder. In the two-sample technique two

samples of different lengths are machined from the same material and scattering parameters

are measured with each sample pressed against the short-circuit termination. Szendrenyi

[13] developed an algorithm for the case in which the length of one sample is precisely twice

the length of the other sample. In this special case, they found an explicit solution.

Mattar and Brodwin [14] have described a variable reactance termination technique for

permittivity determination. Maze [15] has presented an optimized-solution technique where

at each frequency scattering parameters are taken for various short-circuit termination

positions. Dakin and Work [16] developed a procedure for low-loss materials and Bowie

and Kelleher [17] presented a rapid graphical technique for solving the scattering equations.

Other authors have presented methods using measurements on two or more sample lengths

[18]. Most of the literature to date has focused on permittivity determination. In the

few works that have addressed the combined permeabihty and permittivity problem, many

details have been left unresolved.

Recently Chao [19] presented SCL measurements results with a slotted line and also an uncertainty analysis for single frequency measurements. Chao found that accuracy was

reduced when the reflection coefficient is dominated by the front face contribution.

The SCL measurement may use either a fixed or movable short-circuit device. The

advantage of a moving short-circuit termination [2] is the possibility for making many separate measurements at a given frequency with the sample placed in either a high electric or magnetic field region [15]. Generally, a maximum in electric field strength is advantageous for permittivity measurements, whereas a maximum in magnetic field strength is
advantageous for permeability measurements.

When only permittivity is required, a single measurement at a given frequency suffices,

whereas when both permeability and permittivity are to be determined, it is necessary to carry out two independent measurements at each frequency. There are various contributions to the uncertainties in the SCL method. These uncertainties include network analyzer

uncertainties, sample gaps, wall and reflection losses, and measurement of sample dimensions. There are also uncertainties in the location of the sample reference planes and in

the distance from sample to the short-circuit termination. The uncertainty in the network

analyzer parameters are sometimes documented by the manufacturer [3]. In this report we develop relevant equations from first principles. These equations apply

to ANA systems. We will examine the various approaches for combined determination of

permeability and permittivity, and study the uncertainty in the measurement process. The

special case of repeated measurements on a sample of fixed length is treated in detail.
Chapter 2

Theory for Coaxial Line and

Rectangular Waveguide

Measurements of Permittivity and

Permeability

2.1 Theory

The goal of this chapter is to present various approaches for obtaining both the permeabihty and permittivity from transmission hne scattering data. In the TR measurement,

a sample is inserted into either a waveguide or a coaxial line, and the sample is subjected

to an incident electromagnetic field [see figure 2.1]. The scattering equations are found

from an analysis of the electric field at the sample interfaces. In order to determine the

material properties from scattering data, it is necessary to understand the structure of the

electromagnetic field in waveguides. In developing the scattering equations usually only

the fundamental waveguide mode is assumed to exist. In this report we develop the theory

for multimode solutions. However, the numerical algorithms presented will be valid only

for the fundamental mode.

2.1.1 Decomposition into TE, TM, and TEM Modes

In this section we briefly review the theory of modes in transmission lines. It is possible to

decompose the fields in a waveguide at a given frequency into the complete set of TE, TM,
o

I Port 1 o

Outer Conductor

Air

L

II "Sample Air

Outer Conductor

-o

III Port 2

-o

Figure 2.1: A dielectric sample in a transmission line and the incident and reflected electric

field distributions in the regions I, II, and III. Port 1 and port 2 denote calibration reference

plane positions.

and TEM modes. In our model at hand we assume:

• There is a propagation direction in the guide which we call z.

• The cross-sectional area of the guide is perpendicular to z and constant throughout

the length of the guide.

Electromagnetic fields in a sourceless region satisfy

V X V X ^ = -jcjV^ X H + k^E , (2.1)

V X V X // = juVe X E + k^H , (2.2)

where k = —jf is the wave number.

In this report we assume that there are no sources of electric and magnetic fields in the

guide (J = 0) and there no free charge build up (V Z) = 0). Further we assume that the

material parameters are not spatially dependent. However, step function discontinuities are assumed to exist between the sample and air gap. The step function discontinuities in

the equations can contribute a delta function term in derivatives. With these assumptions,

the fields satisfy homogeneous Helmholtz equations.

\/^E-{-PE = ,

The time- dependent fields can be expanded in terms of modes

1 roo ^('^'O = TT / ^u;^K(rT,cj)exp(±7„2)exp(-ju;i) ,

(2.3)

(2.4)

(2.5)
>^(^0 = TT / duJ2Hn{rT,Lv)expi±jnz)exp{-jut) , (2.6)

where rj- is a transverse vector, En, Hn are the amphtudes of the modes, and

ln= J WV/?.4 /27r\2 ^ cL, Va

(2.7)

where 7o,7n are the propagation constants in vacuum and material, respectively. Also

j = \/— 1, c^ac and c/a(, are the speed of light in vacuum and laboratory, u) is the angular

frequency, X^^ is the cutoff wavelength of the nth mode, e^ and //^ are the permittivity and

permeability of vacuum, e^ and //^ are the complex permittivity and permeability relative

to a vacuum.

Since these modes satisfy a Sturm-Liouville problem, we know that the totality of these waves forms a complete set of functions, and therefore an eigenfunction expansion property

exists for this system. The Laplacian separates in the coordinate systems used in this

report, and therefore the fields can be separated into transverse (T) and longitudinal (z)

components:

E = Et-^ EJ , (2.9)

H = Ht + HJ. (2.10)

The component E^ is the generator of the TM mode (see Appendix) and the H^ component

is the generator of the TE mode. Since the TE, TM, and TEM modes form a complete set of functions, we can expand the transverse Fourier-transformed fields as

Erir^to) =

oo Y^{EnTE^M-lnZ) + E'te eMlnZ)] ET{TE)irT)

n = loo + II{^nTMexp(-7n^) + ^~tm ^Xplln^) } ^7'(TM) (^t)

?i = l N-1

+ Y^i^nTEM^M-lnZ) + E;^tej^^ eXY>{-fnZ)} ET(TEM){rT) , (2.11)

n = l

Hrir^oj) =
OO 1 I] ^ iKTE^M-lnZ) - E-TEexp{-fnZ)}{zX ET(TE)){rT))

n^l ^nTE

OO 1 + Yl ^ {EnTM^M-lnZ) " KtM ^MlnZ)}{z X ET(TM){rT))

+ E \l-{KTEM^^Y>{-lnZ) - E^T^EM^^V{lnZ)]{zX ET(TEM){rT)) (2.12)

where N is the number of disjoint conductors and (±) denotes forward and backward

travehng waves. The coefficients E^ depend on the transverse components, and the wave

impedances are

Ztm = -.— , (2.13)

Zte = ^-^ . (2.14)

In

Although the sums for the TE and the TM waves in eqs (3.11) and (2.12) approach oo,

in many problems of practical interest, some of the coefficients in the sums vanish. Two

or more modes may have the same eigenvalue; the eigenvectors in these cases are called

degenerate.

In order to solve eq (2.3) it is expeditious to break up the Laplacian into transverse and

longitudinal components,

\/'E = VlE^^ , (2.15)

where

d'E ,

, =YE , (2.16)

and the transverse Laplacian satisfies by eq (2.3)

dz^

V'tE = -{k' + Y)E = -KE , (2.17)

where fc^ is the cutoff wavenumber. If e has a dependence on transverse coordinates in

terms of a step discontinuity, then 7 also has a transverse dependence. In Appendix A, the

details of the derivations of the fields are reviewed.

2.1.2 Imperfect Sample Geometry

In the case of perfect or near perfect samples and sample holders, 7 is independent of the

transverse coordinates and therefore different eigenfunctions for the transverse components
in the air and sample regions possess an orthogonality condition [see Appendix A]. In such

cases it is possible to match mode by mode, and the coefficients are decoupled. However,

when samples and sample holder are not perfectly formed or are shghtly inhomogeneous,

both fx and e have a weak dependence on the transverse coordinates of the guide and

therefore the different transverse eigenfunctions in the sample are not orthogonal to the

transverse eigenfunctions in the air section. The modes of imperfect samples cannot be

separated and matched mode by mode. The imperfections in the sample generate evanescent waves at the sample-material interface. These modes may propagate in the sample,

but they decay exponentially outside of the sample.

For an imperfect sample, the fields in the regions I, II, and III are found from an

analysis of the electric field at the sample interfaces. We assume that the incident electric

field is the TEio mode in rectangular waveguide and TEM in coaxial line. As the wave

propagates from the air-filled region into the sample, some of the energy carried in the

wave will convert into higher order modes. However, it is necessary to consider only the

transverse components of the fields when matching boundary conditions. In the following

we assume that gaps or other imperfections can exist in and around the sample. We further

assume that the imperfections are such that the Laplacian can be separated into transverse

and longitudinal components. If the imperfections are azimuthally symmetric, then only

the H^ magnetic field component is assumed to exist. If we assume the vector component

of the normalized electric fields Ej, En, and, Eju^ in the regions I, II, and III, we can

write for TV modes

TV Ei= exp(-7oi2) -f5„exp(7oi2) + J]C,(zT)exp(7o,z)J , (2.18)

mci dent wave evanescent

N

En = ^[A(^T)exp(-7„,,z) + £,(2^) exp(7^,z)] , (2.19)

N

Em = ^21 exp(-7.l(^ - ^)) + E [E^[^t) exp(-7o.(^ - L))] , (2.20)

^ ' ,=2^ ' transmitted wave evanescent

where Cj, Dj, Ei^ F^ are the modal coefficients, which may depend on the transverse

coordinates. Also -)oiilmi are the propagation constants of the iih. mode in vacuum and

material respectively. We assume that we are operating the waveguide at such a frequency

that only the fundamental mode is a propagating mode in the air section of the guide. The

other modes are evanescent in the air section of the guide, but may be propagating in the

material-filled section. There may be additional modes produced by mode conversion for

^TEM mode in a coaxial line or the TEiq mode m a waveguide (with a time dependence of exp(ji'cj<)

suppressed)

10
the other components of the electric field, but these are not necessary for specification of

the boundary conditions.

In general, the amplitudes in eqs (2.18) to (2.20) are functions of the transverse coordinates. To find the coefficients, it is necessary to match tangential electric and magnetic

fields at the interfaces and integrate over the cross- sectional area. Since different transverse

eigenfunctions in the air are not orthogonal to transverse eigenfunctions in the sample we

cannot separate a particular mode in the sample and match it to the analogous mode in

the air. The tangential electric field matching yields

N

hi -^ {I - Ski)ck = ^[Akjdj -\- AkjCj] , (2.21)

N

[dk exp(-7„,fcL) -f ek exp{frnkL) = ^

Ajkfj , (2.22)

where 8ij is the Kronecker delta and c, e, and / are the integrated coefficients, A'^ is the

number of modes, and A^j is the matrix of the coefficients of the integrated transverse

eigenfunctions. The transverse component of the magnetic field can be obtained from

Maxwell's equations using eqs (2.21) and (2.22). If we match the tangential magnetic field

components and integrate over the transverse variables we have

[-—4i + —(1 - 4i)q-] = E ^-Akjd, + A,,e,] , (2.23)

N

[-4^^exp(-7^fcL) + -^e^expijrnkQ = -^A^^fj-^ , (2.24)

where L is the sample length and

L,,r = L-^L, + L2 . (2.25)

These boundary conditions yield a linear system of equations for the coefficients. Various

cutoff frequencies and operating frequencies are given in tables 2.1 and 2.2. The difficulties in solving the full mode problem in eqs (2.21) to (2.24) is that the coefficients of the matrix Akj are not generally known precisely unless the complete boundary

value problem is solved for each sample. These coefficients are known only for simple,

well-defined geometries and not for samples with unknown air gaps or complicated inhomogeneities.

11
Table 2.1: Cutoff frequencies for TEiq mode in rectangular waveguide.

EIA WR Band Cutoff frequency (GHz)

650 L 0.908

430 W 1.372

284 S 2.078

187 c 3.152

90 X 6.557

42 K 14.047

22 Q 26.342

Table 2.2: Rectangular waveguide dimensions and operating frequencies in air.

EIA WR Band a (cm) b (cm) TE,o Op(srating frequency (GHz)

650 L 16.510 8.255 1.12- 1.70

430 W 10.922 5.461 1.70 - 2.60

284 s 7.710 3.403 2.60 - 3.95

187 c 4.754 2.214 3.95 - 5.85

90 X 2.286 1.016 8.20 - 12.40

42 K 1.067 0.432 18.0 - 26.5

22 Q 0.569 0.284 33.0 - 50.0

12
2.1.3 Perfect Sample in Waveguide

As a special case of the formalism developed in the previous section we consider a perfect

sample in a perfect waveguide as indicated in figure 2.1. In this case no mode conversion

occurs because the eigenfunctions in the air and sample regions are orthogonal with respect

to cross-sectional coordinates. Therefore the modes may be decoupled and the evanescent

modes are not of concern. This is a special case of eqs (2.21) to (2.25). In this case we

need to be concerned only with the fundamental mode in the guide. The electric fields in

the sample region z G (0, V) for a coaxial line with a matched load and with the radial

dependence written explicitly are

^/ = -[exp(-7o~) + 5iiexp(7oz)] , (2.26)

r

En = -C.2[exp(-7i^) + C3exp(7i^)] , (2.27)

r

^/// = -[5'2iexp(-7o(2-^))]. (2.28)

r When these equations are integrated over the cross-sectional surface area, the radial dependence is the same for each region of the waveguide.

The constants in the field equations are again determined from the boundary conditions.

The boundary condition on the electric field is the continuity of the tangential component

at the interfaces. The tangential component can be calculated from Maxwell's equations

given an electric field with only a radial component. The higher modes in eqs (2.18) to

(2.20) are evanescent in the air-filled section of the guide. TM modes can be treated

similarly. The details of the boundary matching for the Ti^io case are described in a previous report on dielectric materials [26,27]. The boundary condition for the magnetic

field requires the additional assumption that no surface currents are generated. If this

condition holds, then the tangential component of the magnetic field is continuous across the interface. The tangential component can be calculated from Maxwell's equations for an electric field with only a radial component. For a 2-port device the expressions for the

measured scattering parameters are obtained by solving eqs (2.18) through (2.20) subject

to the boundary conditions. We assume that S\2 = 521. The explicit expressions for a sample in a waveguide a distance Lx from the port-1 reference plane to the sample front

face and L2 from the sample back face to the port-2 calibration plane are related. The

S-parameters measured by the device reference planes are related to the S- parameters at the sample face 5' by [26]

5-555 , (2.29)

where

13
$: (2.30)

exp(-<^i) \

exp(-</.2) )

'

and ^1 = —j'yoLi and (/»2 = —J7oJ^2- The S-parameters are defined in terms of the reflection

coefficient T and transmission coefficient z by:

5ii = /?!

S22 — R2

S21 = R1R2

(2.31)

(2.32)

(2.33)

where

Ri = exp(-7oLi) , (2.34)

R2 = exp(-7oL2) , (2.35)

are the respective reference plane transformations. Equations (2.31) through (2.33) are not new and are derived in detail elsewhere [5,28]. We also have an expression for the

transmission coefficient Z:

We define a reflection coefficient by

z ^ exp(--7i)

t£ _ /£0 r

7 70MO

70

(2.36)

(2.37)

For coaxial line the cutoff" frequency approaches 0, (lUc -^ 0) and therefore T reduces to

Crjac I tLK

Clab V ^R (2.38)

Additionally, S21 for the empty sample holder is

55j = /?ii?2exp(-7ola) (2.39)

14
For nonmagnetic materials, eqs (2.31), (2.32), (2.33) contain e'^, e^, L, and Lq, and

the reference plane transformations Ri, R2 as unknown quantities. Since the equations for

S12 and ^21 are theoretically equivalent for isotropic non-gyromagnetic materials, we have

four complex equations, eqs (2.31), (2.32), (2.33), (2.39), plus the equation for the length

of the air line (2.25), or equivalently, nine real equations for the six unknowns. However,

in many applications we know the sample length to high accuracy. For magnetic materials

we have eight unknowns. However, we have frequency data for each measurement. Since

the lengths are independent of frequency we have an over- determined system of equations.

This abundance of information will be exploited in the next chapter.

2.2 Permeability and Permittivity Calculation

2.2.1 Nicolson-Ross-Weir Solutions (NRW)

Nicolson and Ross [5], and Weir [6] combined the equations for 5ii and ^21 and discovered a formula for the permittivity and permeability. Their procedure works well at off-resonance

where the sample length is not a multiple of one-half wavelength in the material. Near

resonance, however, the solution completely breaks down. In the NRW algorithm the

reflection coefficient

Ti=X ± VX^-1 (2.40)

is given explicitly in terms of the scattering parameters where

1 - Vi Vo ^ = 14^- (2-*')

and

V, = 521 + Su , (2.42)

V2 = 5-21 - Su . (2.43)

Note that in the Nicolson-Ross solution the S-parameters must be rotated to the plane of

the sample faces in order for the correct group delay to be calculated. The correct root

is chosen in eq (2.40) by requiring IFi] < 1. The transmission coefficient Zi for the NRW

procedure is given by

^> = 1 - (5„ + 5.)r. • (2.44)

15
If we define

then we can solve for the permeabiHty

f^*R = ^ ^ ^,\ ,

- (2.46)

where Aq is the free space wavelength and Ac is the cutoff wavelength. The permittivity is given by

Equation (2.45) has an infinite number of roots for magnetic materials, since the logarithm of a complex number is multi-valued. In order to pick out the correct root it is necessary to compare the measured group delay to the calculated group delay. The calculated group delay is related to the change of the wave number k with respect to the angular

frequency

'calc. group

.l' '^^^^'

^ (2.48)

df\ c2 A2

L . (2.49)

ifehf^R^m'-^

C2 /i^^y^ _ X

c2 A?

The measured group delay is

where (f) is the phase of Zi. To determine the correct root, the calculated group delays

are found from eq (2.49) for various values of n in the logarithm term in eq (2.45), where

In Z = ln|Z| -i- i(^ -h 27rn), where n — 0,±1,±2,.... The calculated and measured group

delays are compared to yield the correct value of n. Many researchers think of the NRW

solution as an explicit solution; however, due to the phase ambiguity, it is not in the strict sense. Where there is no loss in the sample under test, the NRW solution is divergent at integral multiples of one-half wavelength in the sample. This occurs because the phase of

5ii cannot be accurately measured for small \S\\\. Also in this limit both of the scattering

equations reduce to the relation Z^ -^ 1, which is only a relation for the phase velocity and

16 (
therefore solutions for e^ and fi*j^ are not separable. This singular behavior can be minimized

in cases where permeability is known a priori^ as shown in previous work performed by

Baker- Jarvis [26].

For magnetic materials there are other methods for solution of the S-parameter equations. In the next section we will describe various solution procedures.

2.2.2 2-Port Solution Where Position is Determined Solely by

^airline ^ncl L

In order to obtain both the permittivity and the permeability from the S-parameter relations, it is necessary to have at least two independent measurements. These independent

measurements could be two samples of different lengths, it could be a full 2-port measure- ment, or it could be a 1-port SCL measurement of the sample in two different positions in

the line. In the full S-parameter solution we solve equations that are invariant to reference

planes for e and /z. A set of equations for single-sample magnetic measurements is

p2 _ ^2

5'll522 - 52l5i2 = exp{-27o(Latr - L)] ^

, (2.51)

(^21 + 5'i2)/2 = exp{-7o(/^a.v - ^)} flp2^? • (2-52)

Equation (2.51) is the determinant of the scattering matrix.

Iterative Solution

Equations (2.51) and (2.52) can be solved iteratively or by a technique similar to the NRW

technique. In an iterative approach, Newton's numerical method for root determination

works quite well. To solve the system it is best to separate the system into four real

equations. The iterative solution works well if good initial guesses are available.

Explicit Solution

It is also possible to obtain an explicit solution to eqs (2.51) and (2.52). Let x = (52i5i2 —

511522) exp{27o(Ia2r-J^)} and y = {(52i +5i2)/2} exp{7o(Iatr -^)}, then it can be shown

that the physical roots for the transmission coefficient are ^) - 1 . (2.53)

17
The reflection coefficient is

The ambiguity in the plus-or-minus sign in eq (2.54) can be resolved by considering the

reflection coefficient calculated from Sn alone

a(Z^ - 1) ± jQ^Z''-^2Z^2Sn - ^2) + a^

'^ = ^

25..Z^ ' ("^)

where a = exp{—2joLi). The correct root for T3 is picked by requiring {T^l < 1. Note

that an estimate of Li is needed in eq (2.55). If r2 is compared with T^ then the plusor-minus sign ambiguity in eq (2.54) can be resolved and therefore r2 is determined. The

permeability and permittivity are then

f^R = -f^-V(ln Z -f- 27rjn) , (2.56)

1 - 1 2 7o^

^R = 4[( V)' - 7^(1^ ^ + 27rjn)^]/^^ . (2.57)

The correct value of n is picked using the group delay comparison as described in the

Nicolson-Ross-Weir technique. At low frequencies the correct roots are more easily identified since they are more widely spaced.

2.2.3 Two Samples of Different Length

Solutions for the material parameters exist when scattering parameters on two samples of

differing lengths are measured. Let us consider two samples, one of length L and one of

length aiL as indicated in figure 2.2.

For independent measurements on two samples where |52i | > —50dB over the frequency

band of interest we use only 521 measurements. The measurements obtained on the two

samples are designated as 52i(i) and 52i(2) for first and second measurements:

Z(\ — r^)

^21(1) = exp{-7o(La.> - L)}

^ _ ^2^2 ' (^-^^^

z°'Ui — r^)

521(2) = eXp{-7o(La^r " O^L)}

^ _ ^2ai ^2 ' ^^'^^^

where

Z = exp(-7L), (2.60)

18
Port1 Port 2

Sample 1

(first measurement)

caJi brat ion plane

calibration plane

Sample 2

(second measurement)

Figure 2.2: A dielectric sample in a transmission line for two sample magnetic measurements

and

Z""' =exp(-ai7L) . (2.61)

The reflection coefficient is given by eq (2.37). Equations (2.58) and (2.59) can be solved

iteratively for e^ and /i^.

This solution is unstable for low-loss materials at certain frequencies if the sample

lengths, L and a.\L^ are related so that both materials resonate at a certain frequency

simultaneously. Also with this technique two-sample length measurements are required,

and this increases the uncertainty.

2.3 Measurement Results

The measurement consists of inserting a well-machined sample into a coaxial line or waveguide and measuring the scattering parameters. For waveguide measurements it is important

to have a section of waveguide of length about two free space wavelengths between the coaxto waveguide adapter and the sample holder. This acts as a mode filter for filtering out

higher evanescent modes. There are many roots to the equations for the permeability and

permittivity and caution must be exercised when selecting out the correct root. At lower

frequencies (< 1 GHz) the roots are usually more widely spaced and therefore root selection is simphfied. Another approach to root selection is the measurement of two samples

of differing lengths where the results compared to determine the correct root.

19
6.00

5.75

^R 5.50 —

5.25

5.00

8.3

Nicolson-Ross

Iterative

8.9 9.5 10.1

Frequency (GHz)

10.7 11.3

Figure 2.3: e'^ of a loaded polymer in a X-band waveguide with the full S-parameter iterative

technique.

2.3.1 Measurements without Gap Corrections

Various measurements have been made in waveguide and coaxial line. Some of the results

of these measurements are reported in figures 2.3 through 2.14 for the full S-parameter

technique and in figures 2.15 and 2.16 for the two-sample length method. The measurements

reported in this section are not corrected for gaps around the sample. The effect of the air

gaps is to measure values of the material parameters that are lower than the actual values. In the next section we will discuss ways of mitigating the effects of air gaps.

2.3.2 Effects of Gaps between Sample and Waveguide

Gaps between the sample holder and sample either may be corrected with the formulas

given in the appendix or a conducting paste can be applied to the external surfaces of the

sample that are in contact with the sample holder before insertion into the sample holder.

In figure 2.17 we show a measurement of a nickel-zinc ferrite with and without a gap-filling

grease. The dielectric loss factor is increased shghtly by the gap filling. We suspect that

part of this increase is due to the finite conductivity of the conducting grease.

20
z: '' 0.45

9.6 10.2 Frequency (GHz)

10.8 11.4

Figure 2.4: e^ of a loaded polymer in a X- band waveguide with the full S-parameter

iterative technique.

1.500

1.375

fi^ 1.250

1.125

1.000

Nicoison-Ross

Iterative

8.5 9.08 9.66 10.24 10.82 11.4 Frequency (GHz)

Figure 2.5: //^ of a loaded polymer in a X-band waveguide with the full S-parameter

iterative technique.

21
•R

1.00

0.89

0.78

0.67

0.56

0.45

0.34

0.23

0.12

0.01

-0.1

Nicolson-Ross

Iterative

I "" ——>»>»'ri

8.4 9.6 10.2

Frequency (GHz)

10.8 11.4

Figure 2.6: //'^ of a loaded polymer in a X-band waveguide with the full S-parameter

iterative technique.

R

i

10 12 14

Frequency (GHz)

Figure 2.7: e^ of a ferrite in a X-band waveguide with the full S-parameter iterative technique.

22
R

30

1

I

1

25

20

^ ^^—

15

10

5

n 1 1 10 12

Frequency (GHz)

14

Figure 2.8: e'J^ of a ferritc in a X-band waveguide with the full S-parameter iterative technique.

R

10 12

Frequency (GHz)

14

Figure 2.9: /i';^ of a ferrite in a X-band waveguide with the full S-parameter iterative

technique.

23
Mh

10 12

Frequency (GHz)

14

Figure 2.10: ^^ of a ferrite in a X-band waveguide with the full S-parameter iterative

technique.

14

R

\

3 5 7 9 11 13 15 17

Frequency (GHz)

Figure 2.11: e'^ of a loaded polymer in coaxial line with the full S-parameter iterative

technique.

24

1
^R 2

1 3 5 7 9 11 13 15 17

Frequency (GHz)

Figure 2.12: e^ of a loaded polymer in coaxial line the full S- parameter iterative technique.

A^R 2

1 3 5 7 9 11 13 15 17

Frequency (GHz)

Figure 2.13: /i^ of a loaded polymer in coaxial line with the full S-parameter iterative

technique.

25
A^R 2

3 5 7 9 11 13 15 17

Frequency (GHz)

Figure 2.14: ^^ of a loaded polymer in coaxial line with the full S-parameter iterative

technique.

8.84 10.23 11.62

Frequency (GHz)

13.00

Figure 2.15: e^ of a loaded polymer in a X-band waveguide with TR method for two sample

technique, for three different samples.

26
I^F

8.84 10.23 11.62

Frequency (GHz)

13.00

Figure 2.16: //^ of a loaded polymer in a X-band waveguide with TR method for two

sample technique, for three different samples.

1000

0.1

' (No Conductive Paste)

' (Conductive Paste) "(No Conductive Paste) "(Conductive Paste)

' (No Conductive Paste)

' (Conductive Paste) "(No Conductive Paste) "(Conductive Paste)

%=»JfMtf^,

\ r

0.01

0.001 0.01 0.1 1 10

Frequency (GHz)

Figure 2.17: The dielectric and magnetic parameters of a nickel-zinc ferrite in a coaxial hne

from 300 kHz to 10 GHz with the full S-parameter iterative technique.

27
2.4 Permeameter

In the past permeameters have been used for high permeabiHty materials. Rasmussen [29],

Hoer [30], Powell [31], and Goldfarb [32] have all described various permeameter setups.

In this section we wish to review the theory behind the permeameter.

If a toroidal sample is inserted into an azimuthal magnetic field region, the inductance

is changed. If the inductance of the empty sample holder is compared to the inductance of

the filled holder then it is possible to extract the complex permeability of the material.

Consider a toroid of inner diameter a and outer diameter b and height h. The material

contributes an inductance of [32]

L„ = '-^^^ , (2.62)

and the inductance of the air space is

L, =

f">hHb/a)

^2.63)

27r The net change in the sample inductance when the sample is inserted into the holder is AL^Lm-La , (2.64)

and therefore

27rAL

The magnetic loss may be obtained from consideration of the core loss AR or resistance

28
2.5 Uncertainty of Combined Permittivity and Permeability Determination

In this section an uncertainty analysis is presented. The sources of error in the permeability

and permittivity TR measurement include

• Errors in measuring the magnitude and phase of the scattering parameters,

• Gaps between the sample and sample holder,

• Sample holder dimensional variations,

• Uncertainty in sample length,

• Line losses and connector mismatch, and

• Uncertainty in reference plane positions.

A technique for correcting errors arising from gaps around the sample is given in Appendix B [33,34,35]. Gaps between holder and sample either may be corrected using the

formulas given in the appendix or conducting liquid solder can be painted on the external

surfaces of the sample that are in contact with the sample holder before insertion into the

sample holder, thereby minimizing gap problems. The formulas given in the literature generally under-correct for the real part of the permittivity and over-correct for the imaginary

part of the permittivity. We assume that all measurements of permittivity have been cor- rected for air gaps around the sample before the uncertainty analysis is applied. In order

to evaluate the uncertainty introduced by the measured scattering parameters and sample

dimensions, a differential uncertainty analysis is assumed applicable with the uncertainty

due to 5ii and 8-21 evaluated separately. We assume that the S-parameters are functions

of Si.j{\Sii\, \S21\1 On 1 ^21 1 L,d). We assume that the total uncertainty in e^, where d is the air gap between the sample and waveguide. We assume that the uncertainties for the

physically measured parameters are

Ae'^

(2.67)

(2.68)

29
where a = 11 or 21, A^ is the uncertainty in the phase of the scattering parameter, A|5a| is the uncertainty in the magnitude of the scattering parameter, Ad is the uncertainty in the

air gap around the sample, and AL is the uncertainty in the sample length. The derivatives

with respect to air gap, de'j^/dd^ have been presented previously [26]. The uncertainties

used for the 5-parameters depend on the specific ANA used for the measurements. This

type of uncertainty analysis assumes that changes in independent variables are sufficiently

small so that a Taylor series expansion is valid. Of course there are many other uncertainty

sources of lesser magnitude such as repeatability of connections and torquing of flange bolts.

Estimates for these uncertainties could be added to the uncertainty budget.

2.5.1 One Sample at One Position

For the uncertainty analysis it is necessary to take implicit derivatives of the S-parameter

equations with respect to the assumed independent parameters. It is assumed that the

functions are analytic over the region of interest with respect to the differentiation variables.

The independent variables are assumed to be |52i|, l^nl, ^n, ^21? and L. The derivatives

of the S-parameter eqs (2.31) through (2.33) can be found analytically

dSu.dZ de*r. dZ du%. , dSu.dT deU dV du*r> , . .^ x /^ ^^x -^^^T^ + ^ » ^ic I + -^ FT^r^ + ^V^^ = expO^n , 2.69 oZ de*f^d\Six\ dpi*jid\Sii\ 01 dt)^d\Sii\ 5//^5|5ii|

dS^x dZ de\ dZ dfx}^ dSur dV de)^ d^d^ ^

dz^dt^j,d\S2x\ d^i),d\s,,y^ dv^de},d\S2A di,*j,d\S2,y

' ^' '

dS2i . dZ de\ dZ dfi*j^ 8821 . dV dt\ dV dfi}, ^

dZ ^de^j^d\S,,\ dfi*jid\Sny dT ^de}id\Su\ dfi'j^d\Su\^ ^ ^^

dS2ir dZ de]^ dZ dfx-k 0821 r dT ^e^ dT d^i}^

dz^de},d\82,C ^^i*J,^\82^^'^ dv ^de'i,d\82,C ^^l},^\82,y^''^^^'''^ ' ^'-''^ We can rewrite eq (2.69) - (2.72) as

dS^dZ_ d8u dT de}^

^ dZ de}, ^ dV dt*J d\8n\ ^

V

' 55„ dZ d8u dT df,*j,

30
^ dZ dt]^ dV dt)^^ d\S2x\ ^ V

' A

,dSii dZ dSu dr , du'r,

^ dZ dii},^ dV dfi-k^ d\S2,\

' ^

'

.dSoi dZ 8821 dT , de*ri

( —— + —— — + ^ dZ 5e^ dT dc},' d\Su\

• V

' c

8821 dz 8821 or dfi)^

^ 8Z dfi%^ 8r 8^'j8\Si,\ ~ ' ^ ^^

" V' D

dS2^dZ_ 8821 8T 8e''ji

^ 8Z 8e*j^^ 8T de*f,' 8\82i\

"^

^^ ^..^ ' c

8821 8z 8821 8r 8fi*j^ dz d/xji 81 dnji d\82i\

D

where we have defined parameters A, B, C, and D. If we let

we can solve for the derivatives that have been taken with respect to the independent

parameters in eqs (2.73)- (2.75):

^fL _ ^Mj^ (2.79)

d\8n\ [A-^]

31
d\Su\

dL

dL

dt'

C de*ji

''~Dd\Su\

BF -DE

BC -AD ' E-A dL

B

dOu^=;|5n|-^'^ d\S111

9f^R MO I ^/^fl

dp-R

= j\Su\

= J\S2\\

J\S:21

d\S,,\

'

d^R

d\S2A

'

8621 d\S2i\

a5n (l+r2^2)(,.2_i)

ar

dSu

(1 -r222)2

2zr(r2-i)

5Z (1 -r2Z2)2 '

5521 (1 -r2)(z2r2 + i)

dZ

dS.21

dv

(1 -r2Z2)2

2zr(z2- 1)

(1 -r2Z2)2

lOf^R^Of^O^^

de'fi 7(7 + 7o/^fi)^

'

ar ar

+

2707

^^i*R ^hd^R (7 + 7o/^fl)^

'

az __ dL~~

dZ Lii*pUj'^eQ}iQ

di R 27

7exp(-7L) ,

exp(-7L)

(2.80)

(2.81)

(2.82)

(2.83)

(2.84)

(2.85)

(2.86)

(2.87)

(2.88)

(2.89)

(2.90)

(2.91)

(2.92)

(2.93)

(2.94)

32
21| 30-

Figure 2.18: The derivative of e^ by |52i| vs L/A^ with e*j^ ^ (5.0,0.02), fi*j^ = (2.0,0.03).

dZ L€*j^u>^eofiQ exp{-jL) (2.95)

df^R 27

The measurement bounds for S-parameter data are obtained from specifications for a network analyzer. The dominant uncertainty is in the phase of 5n as | Su |^ 0. The

uncertainty in |52i| is relatively constant until |52i| < -50 dB, when it increases abruptly.

The various derivatives are plotted in figures 2.18 through 2.27.

In figures 2.28 through 2.31 the total uncertainty in e^ and /x^ computed from S21

and 5n is plotted as a function of normalized sample length. For low-loss and high-loss

materials at 3 GHz with various values of e|j and the guided wavelength in the material

given by

Am —

27r

^2 {^/^+^') ^r

(2.96)

27rN2 i"J

In figures 2.28 through 2.31 the error due to the gap correction is not included, nor

are there uncertainties included for connector repeatability or flange bolt torquing. The

maximum uncertainty for low-loss materials occurs at multiples of one-half wavelength.

Generally, we see a decrease in uncertainty as a function of increasing sample length.

Also, the uncertainties in the S-parameters have some frequency dependence with higher

frequencies having larger uncertainties in phase.

33
20 -

j

1 1

S21

acp

-^/////-

39,,

' 1 ' y

V

-20 -

1

-40 -

1 1

() 1 2 3

'm

Figure 2.19: The derivative of e^ by ^21 using 821 vs LjXm with e^ = (5.0,0.01), /i|j

(2.0,0.03).

90

70-

50-

ae"

aa21 30-

10-

10

V

1

S21

1

*m

Figure 2.20: The derivative of e^ with respect to ^21 using 821 with e^ = (5.0, 0.01), /Xjj (2.0,0.03).

34
dL

\J

1 1

400 - /^

^21 -

-800 7 -

1200 -1 -

1600 1

X m

Figure 2.21: The derivative of e'^ with respect to L using S21 with ej^ = (5.0,0.01), /u^

(2.0,0.03).

Figure 2.22: The derivative of e^ with respect to L using 52i with ej^ = (5.0, 0.01), //^

(2.0,0.03).

35
200 -

1 1

-

150 -

Sii -

100 - -

50

n uuVLLju

Figure 2.23: The derivative of e^ with respect to |5ii| with t)^ — (5.0, 0.01), /x^ = (2.0,0.03).

Xm

Figure 2.24: The derivative of e'^ with respect to |5ii| withej^ = (5.0,0.01),//^ = (2.0,0.03).

36
/ \

0.2

AAAA/ / \ /v vv/ V

0.4 -. / Sii -

0.6 ^ \J

1 1

-

'm

Figure 2.25: The derivative of e^ with respect to 9u using Su with e*j^ (2.0,0.03).

(5.0, 0.01), /i|,

2.4 7 \A/^yvV

2.0

[

v ^

/

Si,

1 6

/

1 1

1 2 3

L

*m

Figure 2.26: The derivative of e^ with respect to ^n using 5ii with ej^ = (5.0,0.01),//^

(2.0,0.03).

37
dL

-400-

-800-

-1200-

1600

1 1

- /

S"

-

-

1 1

-

m

Figure 2.27: The derivative of e^ with respect to L using ^n with e^ = (5.0, 0.01), /i|j =

(2.0,0.03).

0.11

Ae,

0.09

0.07

0.05

0.03

0.01

K = (5,0.05)

_L 0.5 1.0 1.5 2.0 L/Xr.

2.5 3.0 3.5

Figure 2.28: The relative uncertainty in e^(c<;) for a low-loss material as a function of

normalized length, with //|j = (2, 0.05), ejj = (10, 0.05) and (5, 0.05).

38
0.11

0.09

0.07

0.05

0.03

0.01

Sr = (10,0.05)

^^l= (2,0.05)

0.5 1.0 1.5 2.0 L/Xr.

2.5 3.0 3.5

Figure 2.29: The relative uncertainty in ^'ji{u?) for a low-loss material as a function of

normalized length, with /^l^ = (2,0.05), e^ = (10,0.05) and (5,0.05).

0.080

0.008,

-

1 r

1 —^

1 - - - -

- - - -

- -

/

/\ *l= (10.05) ^ W; = (2.05)

- -/A^^^

_^

«; = (5.05) w'.V/-^ A*; =(2.0.5)

1 1 1

1 2 3 L/X,

Figure 2.30: The relative uncertainty in t'ji{oj) for a high-loss material as a function of

normalized length, with //^ = (2,0.5), e^ = (10,0.5) and (5,0.5).

39
0.040

0.035

0.030

0.025

0.020

0.015

0.010

0.005

Figure 2.31: The relative uncertainty in ///{(w) for a high-loss material as a function of

normalized length, with /^^ = (2.0,0.5), e^ = (10.0,0.5) and (5.0,0.5).

100

^'n

10

0.001 0.01 0.1 1 FREQUENCY (GHz)

Figure 2.32: The real part of the relative permittivity ej^(w) for a nickel-zinc compound

with uncertainties.

40
100

10

^'k

0.1

0.001 0.01 0.1 1 FREQUENCY (GHz)

10

Figure 2.33: The imaginary part of the relative permittivity e'j^{oo) for a nickel-zinc compound with uncertainties.

1000

100

10

f^R

0.1

0.01

0.001 0.01 0.1 1 10

FREQUENCY (GHz)

Figure 2.34: The real part of the relative permeabihty jJ-'jiiy)) for a nickel-zinc compound

with uncertainties.

41
1000

100

10

/^'A

0.1

0.01

0.001 0.01 0.1 1 FREQUENCY (GHz)10

Figure 2.35: The imaginary part of the relative permeability /^^(aj) for a nickel-zinc com- pound with uncertainties.

42
In figures 2.32 through 2.35 a measurement of a nickel-zinc ferrite compound is given

with associated uncertainties. Uncertainties increase at high and low frequencies. At

high and low frequency extremes the uncertainties in phase increase. Also, the scale is logarithmic which distorts the lengths of the error bars.

2.5.2 Two Samples of Differing Lengths

Another method to determine permittivity and permeability is the measurement of two

samples with differing lengths. The advantage of this method is that each sample resonates

at a different frequency and therefore Su can be appreciable over the entire frequency band.

We assume the S-parameters are functions of Sij{\Smn\,(^mn^ Li^ L2). The parameters

used for measurements on materials of low to medium loss are

since it is acceptable down to -40 dB. We assume that the lengths of the samples are Li

and Z/2 = aLi. Due to the two lengths, there are transmission coefficients for each sample

Zi=exp(-7Li) , (2.98)

Z2 = exp(-Q7Li) . (2.99)

The relevant partial derivatives of eqs (2.97) are:

d821(1) dZi de)^ dZi dfi}^

dZi dt\d\S2\{i)\ dpi}id\S2i(i)\

dV de*j^d\S2i(i)\ ^^i)^^\S2^l)\ = exp(j^i), (2.100)

dS2x(2) r dZ2 de}^ dZ2 dfi)^ 0821(2) . dV de*j^ dT dfi*n ^

dZ2 ^de*j,d\S2m)\ dfi},d\82m)y dr he*j,d\82m)\ dfi*j,d\82i(i)y '^' ^

d82i(i) . dZ^ d€}i dZ, dji'k 8821(1) r dV de*j^ dT d^*p_ ^ f2 102)

dz,^dej,d\S2x(2)\ dii*j,d\82i(2)Y dv ^dej,d\82m\ dii'kd\82m\

43
dS2i{2)^dZ2 de\ dZ2 dpi'j^ dZ2 5e^5|5'2i(2)| dfi*jid\S2i(2)\

dT de'ji d\S2i{2)\ dfi*jid\S2i(2}\'

exp(i^2) . (2.103)

We can rewrite eqs (2.100) through (2.103) as

dS2i{i)dZi dS2i{i) dT de}^

dZi de'fi dT 5e^ d\S2iii)\

V

'

,8821(1) dZ dSoKi) dT , dup , „ ^ , ^ ^^jij^-WJ^d^l = -^^'^^ ' (2.104)

8821(2) dZ2 dS2i{2) dT de*j^

dZ2 de*j, 8T dt),' d\82m,\

A2

d82ii2) dZ2 d82ii2) dT d/x^^ ^

^-dz^dJF,^

^ ^^d^,\ d\s^\ = '

^'-'"'^

S2

d82i(i)dZ^ d82i(i) dT de)^

5Zi ae^ ar ^e^ d\82i(2)\

,d82iii\ dZi 5521(1) 5r , du'ff ^^ ^^,

5Zi 5/x^ dT dfi'ji d\82i(2)\

Bi

44
dZ2 de*j^ or de*j^' d\S2i (2)1 A2

, , ^521(2) ^^2 , ^521(2) ST dl2*fi

B2

where we have defined parameters Ai, Bj, v42, and B2. Also for the relevant derivatives

with respect to length, we find

^ dS2i(i)dZi dS2i(i) dV 5e^

dZi dt*j^ dT dey dL^

I ( dS2m dz, dS2m) dr d^i\ 882^1) dz, _

^521(2)^^2 <952i(2) ^r de*j^

A2

, dS2xi2)dZ2 8821(2) dT djL^ 8821(2) ^^2 _ . .^ ^..^

"^^ aZ2 5/z^"^ 8T 8fi''j,^ 8L, 8Z2 8L,~^- ^^''''^^

<

v^

'

"^ V

' B2 F2

We now can solve for the derivatives that have been taken with respect to the independent parameters in eqs (2.104) through (2.107)

54 _ 52exp(;^0

^2.110)

d\S2ii,)\ A,B2-B,A2 '

(2.111)

8n*ji A2 de*j^

^|'5'2i(i)| B2 |5'2i(i)|

8t\ Biexp{j92)

(2.112)

81821(2)1 A2B1-B2A '

(2.113)

5//^ Ai ae^

^|'S'21(2)| Bi |52i(2)|

45
de*j^ ^ E1B2 — F2B1

dTi ~ B1A2- AiB2

dZ, {l-T^ZfY

5521(0 _2Z,T{Zf-l)

dr (1 -T^Zf)^

'

d^R 7(^07 + 70/^)^

(2.114)

(2.115)

-^ = J 1^^21(1) 1^1^ 1 , (2.117)

a^ 1(2 = ^'5^-«l^sg^' (2-"^)

(2.120)

(2.121)

(2.122)

+ 7 ^^ ' (2-123)

df^R /^fl<9efl (7 + 70/^)^ __ = _^exp(-7l), . (2.124)

^ = -|^exp(-7L), (2.125)

O^R ^7

^-T = —:^ exp(-7L) . (2.126)

In figures 2.36 through 2.37, the total uncertainty in e^ and //^ computed from 52i is plotted as a function of normalized sample length, for low-loss and high-loss materials at 3

GHz with various values of e^.

When the length of one sample is twice the length of the other sample, we see instability

at frequencies corresponding to n\m/2. Generally, we see a decrease in uncertainty as a

46
Ae'r 0.7

0.07

0.007

^8r = (30, 0.2)

^R = (1.0)

0.5 1

Figure 2.36: The relative uncertainty in €'j^{to) for a low-loss material as a function of

normalized length for the case when Li = O.5L2 for two different permittivities.

function of increasing sample length. Also, the uncertainties in the S-parameters show

some frequency dependence. In figure 2.37 the ratio of sample lengths is \/2. In this case we see greater stability over the frequency range than in the case where the range is 0.5. Resonances in the solutions will occur when L — nA^/2 and ah = mXml'2. simultaneously,

where m is an integer.

47
Ae.

0.35

0.07

0.007

1

1 I.I

^R = (3. 0)

L2= V2L,

e; = (5,o.2)

-\ \/ V \ / e; = (30,0.2) V

\1 ^^-N-^

1 /1 V^2

Xr•m

Figure 2.37: The relative uncertainty in e'j^{u!) for a low-loss nnaterial as a function of

normalized length for the case when L] — 2L\ for two different permittivities.

48
2.6 Uncertainty in Gap Correction

The correction for an air gap between the wall of the sample holder and sample is very

important for measurements of high permittivity materials. In addition, the uncertainty

in the gap correction is very important for high permittivity materials and may actually

dominant the uncertainties of the measurement. In appendix C the gap correction is worked

out in detail. In this section the uncertainty in the gap correction will be worked out.

2.6.1 Dielectric Materials

Waveguide Gap Uncertainty

The uncertainty due to an air gap between sample and holder can be calculated from the

partial derivatives of e^ with respect to gap thicknesses, d. The relevant derivatives for

waveguide are given by

Coaxial Gap Correction

For coaxial hne the relevant derivatives are given by

dRs ''R3{L,-eUL,)^''^''R^{^ "" i?3(L3 - e^^Li) "" RsiLs - e'^f^Li^

dR, ''^"''-"^Tji^^ W-} ' ^^"''

+ TT^] • (2-132)

y2-rt3 ^2^3

49
2.6.2 Magnetic Materials

Waveguide Gap Uncertainty

The uncertainty due to an air gap between sample and holder can be calculated from the

partial derivatives of fi*j^ with respect to gap thicknesses, d. The relevant derivatives for

waveguide are given by ^= (l-.U)^, (2.133)

Coaxial Gap Correction

For coaxial line the relevant derivatives are calculated using

Ki =\nR2/R, , (2.135)

K2=-\nR3/R2, (2.136)

A'3 = In /?4/i?3 , (2.137)

K4=^nRjRi, (2.138)

as follows

dfJ-'cR 1

r

' ^^^4 A'l A'3,

du'f) 1 r / A'4 A^l A'3,

i^ = 7v^[-''-7?; + 7f; + > + ^I' (2..40)

dK hh2 = -^'-^;a^

• ^'-'"'^

2.6.3 Higher Order Modes

The field model assumes a single mode of propagation in the sample. Propagation of higher

order modes becomes possible in inhomogeneous samples of high dielectric constant due

to changes in cutoff. Air gaps also play an important role in mode conversion. Generally,

the appearance of higher order modes manifests itself as a sudden dip in |5ii|. This dip

is a result of resonance of the excited higher order mode. We can expect point-by-point

50
TR models to break down near higher order mode resonances for materials of high dielectric constant or inhomogeneous samples. Optimized, multi- frequency solution techniques

fare better in this respect. The characteristic of the higher order modes are anomalies in

the scattering matrix at and around resonance. Higher order modes require a coupling

mechanism in order to begin propagating. In waveguide and coaxial line the asymmetry

of the sample promotes higher order mode propagation. In order to minimize the effects

of higher order modes, shorter samples can be used to maintain the electrical length less than one-half guided wavelength. Also well machined sample are important in suppressing

modes. Higher order modes will not appear if the sample length is less than one-half guided

wavelength of the fundamental mode in the material.

Mode Suppression

It is possible to remove some of the higher order modes by mode filters. This would be

particularly helpful in cyhndrical waveguide. One way to do this is to helically wind a fine

wire about the inner surface of the waveguide sample holder, thus eliminating longitudinal

currents and therefore TM modes. Another approach is to insert cuts in the waveguide

walls to minimize current loops around the waveguide.

51
Chapter 3

Optimized Solution

3.1 Introduction

As indicated in the previous chapters, various numerical strategies have been employed for

reducing 1-port and 2-port scattering data for both nonmagnetic and magnetic materials.

The vast majority of the work in this area has involved the determination of permittivity

and permeability by the reduction of scattering data frequency by frequency, that is, by the

explicit or implicit solution of a system of nonlinear scattering equations at each frequency

(see [11,33]; as an example of a multifrequency approach see Maze et al.[15]).

What is lacking in the literature are practical, robust, numerical reduction techniques

for more accurate determination of permittivity and permeability in transmission hnes.

Reliable broadband permeability and permittivity results for low-loss, medium-to-high dielectric constant materials are hard to obtain with transmission line techniques. Coaxial

line measurements are particularly hard to obtain due to air gap influences and overmoding. Traditional transmission line numerical techniques have difficulties to an extent that

render these techniques of limited use for low-loss materials and for high dielectric constant

materials. Difficulties arise with these methods for magnetic materials in that numerical

singularities can occur at frequencies corresponding to integral multiples of one half wavelength. These instabilities arise from the fact that for low-loss materials both 521 and 5ii

become equations for the phase velocity, and the permittivity and permeability therefore

enter as a product. These instabilities limit the acquisition of precise broadband dielectric and magnetic results in the neighborhood of a resonance. Another problem pertains

to high dielectric constant materials. High dielectric constant materials are usually hard

to measure since the theoretical models are limited to a single, fundamental mode and

the data contain both fundamental and higher order mode responses. Further, point-bypoint reduction techniques for magnetic materials contain large random uncertainties due

52
to the propagation of uncertainties through the equations. For nonmagnetic materials the

propagation of errors is less of a problem.

In our search for better reduction techniques we have found that nonlinear optimization,

which minimizes the squared error, are a viable alternative solution. Optimization-based

data reduction has an advantage over point- by-point schemes in that correlations are al- lowed between frequency measurements. In nonlinear regression, if deemed appropriate, it

is not necessary to even include Sn data in the constraint equations. Another advantage

of regression is that constraints such as causality and positivity can be incorporated into

the solution.

This chapter presents a method for obtaining complex permittivity and permeability

spectra from scattering parameter data on isotropic, homogeneous materials using nonlin- ear regression. We solve the scattering equations in a nonlinear least-squares sense with

a regression algorithm over the entire frequency measurement range. The complex permittivity and permeability are obtained by determining estimates for the coefficients of a truncated Laurent series expansion for these parameters consistent with linearity and

causality constraints. The procedure has been successfully used for accurate permittivity

and permeability characterization of a number of different samples where point-by-point

schemes have proven to be inadequate. The details of the numerical method have been

presented in [36]. The problem applied to microwave measurements is presented in this

chapter. The method can easily be extended to the analysis of multi-mode problems and the

determination of experimental systematic uncertainty. The novel features of our algorithm

are:

• The algorithm finds a "best fit" to the 2-port scattering equations using a nonlinear

least-squares solution for the permittivity and permeability.

• The algorithm uses fitting functions that satisfy causahty requirements.

• The numerical technique allows slight variations in the sample and reference position

lengths to compensate for measurement errors and sample imperfections.

• The method allows the de-emphasis of frequency points with large phase uncertainty.

• Statistics related to the solution parameters are automatically generated.

• The technique can force positivity of the fit functions.

• It is possible to determine both complex permittivity and permeability from measure- ments of a single scattering parameter on a 1-port or a 2-port taken over a frequency

band.

53
3.2 Model for Permeability and Permittivity

In the optimization procedure the S-parameter eqs (2.31) through (2.33) for single-mode

problems, or eqs (2.21) through (2.24) for problems with higher order modes are solved for

the material parameters by the optimization routine. For higher order mode problems the

matrix elements Aj^ in eqs (2.21) through (2.24), which corresponds. to the voltage in each

mode, can be determined by the optimization routine. However, we usually consider only

the primary mode.

The unknown quantities are Li, Z/2, L, A^, and /ifl(uj) and e*j^{u). Some of these parameters, such as the lengths and cutoff wavelength, are known accurately within measurement

uncertainty. Obviously the parameters of interest cannot be allowed to vary into nonphysical realms. The problem is to use an optimization routine to determine the model

parameters that are consistent with the scattering data and the physics of the problem.

3.2.1 Relaxation Phenomena in the Complex Plane

The numerical model requires an explicit functional form for //^ and e]^ to reproduce the

four S-parameters consistent with the data, for all the 7i frequency observations.

The general form for ^J^(c<-') and c*j^{u)) should be causal see Appendix C; that is, it must

satisfy a Kramers-Kronig relation. If the zeros and poles of a complex function are known

over the complex plane, the function itself is known.

The Laplace transform of the real, time-dependent permittivity satisfies

e(f,s) = / e[f,t)e-''dt . (3.1) Jo For stability, there can be no poles in the right-half side of the s-plane. Since e[t) is real

it can be shown that the poles and zeros are confined to the negative real s-axis of the

5-plane, and the poles which are off the real s- axis must occur in complex conjugate pairs

[37].

Assuming linear response a constitutive relationship in an isotropic medium between

the displacement and electric fields is

D{xJ) = eoE{xJ) + eo / GiT)E{x4 ~ T)dT . (3.2) Jo With this definition the permittivity is

roo e^(c^) - 1 + / G'(r)e-^""c/r . (3.3)

The response function G{t) for an incident electric field can sometimes be represented as a series of damped sinusoids of the form

54
G{t) = J^A„exp(-(an +j6„)t) (3.4)

We assume that the material parameters can be adequately modeled by a series of simple poles. The terms e~""' relate to relaxation and the terms e"-'^"' relate to resonant

phenomena. Since the Laplace transform of eq (3.4) is

G{s) = Y.^n

1

S-\- ttn + jhn

(3.5)

for stability (no increasing time domain exponentials) there can be no poles or zeros in the

right half-side of the 5-plane [37], [38]. In order to maintain the reality of e(5), any poles

off the imaginary axis in the o^-plane must be conjugate poles of the form uj = ja:tb where

a and b are real, positive numbers as indicated in figure 3.1. These conjugate poles are of

the form

1 1 +

s + a + jb s -^ a — jb

and are related to resonant phenomena.

We assume that the permittivity can be expressed as

(3.6)

e^j-C n!^^^+Eh 1 +

1

n U^ + Pn) n J^ ^ 0.n+ jbn J^O -f a„ - jk

(3.7)

Here Zn and p^ are the zeros and poles due to damped exponentials respectively, (a„ zkjbn)

are the complex conjugate poles, and C is a complex constant.

55
X X X

imag(s)

leal(s)

Figure 3.1: Poles (X) and zeros (O) in the complex s-plane and conjugate poles off the real

s-axis.

For more complicated polarization phenomena other relations for permittivity could be

used. For a continuous distribution of relaxation times

y{r)

Jo 1 + U}^T^

e = [e(0) - e^j / d

Jo 1 + UJ^T^

(3.8)

(3.9)

where y(r) is a distribution function. Various expressions for y yield various relations for

permittivity. Presently we use eqs (3.11) and (3.12) in our calculations. The HavriliakNegami model for materials assumes a single nonsimple pole on the negative, real 5-plane

axis:

C +

A

(3.10)

{\ + [j\B\uY-r '

where (3 and a are in the interval [0, 1] and B is real. Limiting cases of this model are

(1) the Cole-Davidson model when a = 0; this model works well for some hquids and

solid polymers, (2) the Cole-Cole model when /^ = 1; this model has been used to describe

relaxation behavior of amorphous solids and many liquids. A simple Debye model (/5 = 1 and a = 1) is very limited and works well only for materials that contain a single relaxation

time in the frequency range of interest.

Heterogeneous materials and polymers usually have a very broad relaxation spectrum

and as such have a response of the form of a power law such as s~". This behavior can be

obtained from the Cole-Davidson model when |jB|u; >> 1. In our present algorithm we assume a more general model than the single Debye relaxation model a truncated Laurent series is used for /^/^(u;) and e^(u;). This expansion has

56
generally yielded excellent results

'«(") = ^° + rT&i:; + (rT7W- ^'-''^

where Bi are real numbers. The pole information yields constraints on the constants used

in the Laurent series expansion. For example, it is required that Bi is a real number.

For a typical measurement on a network analyzer there may be 400 frequency points and

at each point all four scattering parameters are taken. The problem is overdetermined since

for n frequency measurements, if we assume known lengths, there are Sn real equations

for the unknown quantities in the Laurent series. This over determination can be used

at frequencies in the gigahertz range to find corrections to sample position and cut off wavelength.

The approach for determining the complex parameters Aj, Bj is to minimize the sum

of the squares of the differences between the predicted and observed S-parameters,

mm E4-4II , (3-13)

where the measured vectors are denoted by Sij = (5'jj(u;i), 5jj(a;2)? --j '5'ij((^n)) a-nd where

Pij is the predicted vector. Hence, the problem consists of finding the norm solution to these equations.

3.3 Numerical Technique

3.3.1 Algorithm

The solution currently uses a software routine called orthogonal distance regression pack

ODRPACK [39] developed at the National Institute of Standards and Technology. This

routine is an extended form of the Levenberg-Marquardt approach. This procedure allows

for both ordinary nonlinear least-squares, in which the uncertainties are assumed to be

only in the dependent variable, and, orthogonal distance regression, where the uncertainties appear in both dependent and independent variables. First-order derivatives for the

Jacobian matrices can be numerically approximated (finite difference approximation) or can be user-supplied analytical derivatives. The procedure performs automatic scaling of

the variables if necessary, as well as determining the accuracy of the model in terms of machine precision. The trust region approach enables the procedure to adaptively determine

the region in which the linear approximation adequately represents the nonhnear model.

57
Iterations are stopped by ODRPACK when any one of three criteria are met. These

criteria are: (1) the difference between observed and predicted values is small, (2) the

convergence to a predicted value is sufficiently small, and (3) a specified limit on the

number of iterations has been reached.

Initial guesses for e|j and /i|j are obtained from explicit solutions of Stuchly [7] or Weir

[6]. The most significant input parameters for modehng permittivity and permeability are the initial values for Ai and Bi. Sensitivity to the initial solution for these parameters is discussed below. All additional parameters are initialized to 0. When measurements of length and scattering parameters of a sample are taken, there

are systematic uncertainties.

An orthogonal distance regression model provides the modeler with the additional ability to assume that the independent variable, in this case, frequency, may contain some

uncertainty as well. Allowances for these types of uncertainty can, in some cases, greatly

improve the approximation. For this model and the samples tested, the errors in the in- dependent variables are sufficiently small that an ordinary least-squares approximation is adequate.

Model parameters such as sample length, sample position in the waveguide, and cutoff

wavelength could contain a systematic uncertainty. These parameters were allowed to

vary over a hmited region, and the optimization procedure chooses optimum values for the

parameter. This procedure assumes .that systematic measurement errors can be detected by

the routine. For example, inserting a sample into a sample holder introduces an uncertainty

in the sample position Li, so we include with Li an additional optimization parameter /^n

in Ri to account for positioning uncertainties,

R, =exp(-7,[Li+/?Li]) . (3.14)

Also for i?2 /?2=exp(-7,[L2 + /5L2]) . (3.15)

The routine requires that the length corrections be within a prescribed range which represents physical measurement uncertainty. The length of the sample L is completely determined by

L = La^r - [Li + L2 + /?Li + /3l2) (3.16)

and is also implicitly parameterized by the values of f^n and /?^2-

Due to inaccuracies in machining of the sample holder there is an uncertainty in the

cutoff wavelength of the guide. We account for this by the introduction of an additional

optimization parameter A^ -^ Xc-\-(^\. We constrain this variation to be within measurement

accuracy.

The model can use various combinations of the available data to estimate both the

relative permeability and permittivity from scattering data. For example, ^21 or Sn alone

58
-0.4 ., , 11 I 1^

1 1 t i ^ -0.6

CC -08

-

1,0 ^ MH-^ Ut./

8 9 10 11 12 13 Frequency (GHz)

9 10 11 12 Frequency (GHz)

0.4 ^ 0.2

<5. 0.0

4) °C -02

-0.4

9 10 11 12 13 Frequerx:y (GHz)

(a)

1.0

" /^^\ '

0.5 \

§0.0

°C-0.5 : \ :

-1.0

, , , ^-^

9 10 11 12 13 FrequerKy (GHz)

(b)

Figure 3.2: Predicted (solid line) and observed (dots) parameters for a barium titanate

compound (a) and cross-linked polystyrene in (b).

can be used to obtain both permeability and permittivity. This can be contrasted with

point-by-point techniques where both S21 and Su are required. Also, magnitude alone can similarly be used. Magnitude data have the advantage of requiring no reference plane

rotation.

The technique works well for short-circuit line measurements. For short-circuit lines it

is possible with this technique to obtain both the complex permittivity and permeability

from a single broadband measurement on one sample at a single position in the line.

3.3.2 Numerical Results

The model predictions are formed by inserting eqs(.3.11) and (3.12) into (2.31) and (2.33) or (4.7) and then finding the unknown coefficients in the equations for e|^ and /i^ that

produces the least square error. In figure 3.2 the experimental results are given for a barium titanate compound and cross-linked polystyrene. These samples required 21 and

40 iterations respectively.

The difference between the predicted S-parameter and the observed values reveals the

presence of systematic uncertainty, as shown in figure 3.3, in the automatic network analyzer

(ANA). Additional tests revealed the source of the systematic error did not appear to be

related to the material tested in the waveguide. In fact, uncertainties produced for the cross- linked polystyrene sample closely resemble the S-parameter data for an empty waveguide;

59
0.002 -

:^ 0.0

cc -0.002

-0.004

0.001 -

CO

75 0.0

0) CC

-0.002 -

» « —, . —

*

~ "^ .' - —I JJ:!! \ I I i_^

10 11 Frequency (GHz)

(a)

10 11 Frequency (GHz)

(b)

12 13

12 13

Figure 3.3: Systematic uncertainty as indicated by the residual plot (the difference between

the observed and predicted values) for the case on an empty waveguide. The solid line is at residual.

we conclude therefore that much of the systematic error is due to calibration uncertainty

and joint losses at connector interfaces. For the barium titanate compound sample there

is both the fundamental mode response and smaller resonances related to higher-order

modes. As shown in figure C.l, the model interpolates a fundamental mode. This raises

the possibility of extending the model to incorporate higher-order modes by extending the

theoretical formulation of the problem.

It is easy to move the sample in the holder inadvertently when connecting the sample

holder to the port cables. Positioning errors of the sample in the air line can result in large

error in computed material parameters. The numerical algorithm can adjust for positioning

errors by adjusting Li or L2 slightly. The effects of positioning error can be seen in figure

3.4.In this example the routine predicted that the position of the sample was off by 0.8

mm.

60
^ 0.0

0.45

0.045

5 10 Frequency (GHz)

(a)

15

5 10 Frequency (GHz)

(b)

15

Figure 3.4: Measured real part of 5ii, measured (. . . ) and predicted (- -), for a glass sample (a) with positioning error for L^ and (b) the solution when the algorithm

adjusts for the positioning error.

3.4 Permittivity and Permeability

3.4.1 Measurements

In this section we present the measured and calculated permittivity and permeability.

Cross-linked polystyrene and the barium titanate compound are nonmagnetic and therefore

fi*fi ~ I. Comparison of the optimized solution to a point-by-point solution is shown in figure

3.5. In figures 3.6 through 3.9 results for four samples are given.

As a check we made an independent measurement of the barium titanate compound in an X-band cavity where the results were e^ = 269 at 10 Ghz. This result can be compared

to the results in figure 3.6. Finally a result of another barium titanate compound is given

below.

3.4.2 Robustness of the Procedure

Since the transmission coefficient contains a periodic component, there is more than one

solution to the system of equations. Each root of the equation has a neighborhood around

which convergence will occur for initial guesses in that region. The robustness of a mathematical procedure is related to how well the algorithm treats the neighborhood around the

61
^R 16 h--—^

Figure 3.5: Permittivity in the frequency range of 0.045— 18 Ghz for the optimized solution

( ) and point-by-point technique ( ).

265 ' 1 1

1 \_ -r

264.5 ^^

2640

^^ 263.5

r 1 \ 1 1 1 8 9 10 11 12 13 8 9 10 11 12 13 Frequency (GHz) Frequency (GHz)

(a)

0.01

Eq

1 1 1 1 1 r

8 9 10 11 12 13 Frequerxry (GHz)

I

'

I -I 1 L 8 9 10 11 12 13 Frequency (GHz)

(b)

Figure 3.6: Permittivity for barium titanate compound (a) and cross- linked polystyrene

(b), point-by-point method (. . . .), optimized solution ( ).

62
9 10 11 12 Frequency (GHz)

13 9 10 11 12 13 Frequency (GHz)

8 9 10 11 12 Frequency (GHz)

8 9 10 11 12 13 Frequency (GHz)

Figure 3.7: Permittivity and permeability for a loaded-polymer, point-by-point method (.

. . .) , optimized solution ( ).

120

119

118

117

116

£ 115

114

113

112

111

110

CAVITY

9 10 11 FREQUENCY (GHz)

12

Figure 3.8: Real part of permittivity for a barium titanate compound (a) (o o o) cavity (d)

(...) optimized solution.

63
0.01

0.009

0.008

,_ 0.007

z^ 0.006 < 0.005

w 0.004

O

^ 0.003

0.002

0.001

9 10 11 FREQUENCY (GHz)

12

Figure 3.9: Imaginary part of permittivity for a barium titanate compound (a) o o o -cavity

and ()d ... - optimized solution.

correct root. The existence of alternative optima in the mathematical model requires an

accurate initial guess in order to converge to the correct solution. Typically convergence

occurs after about seven iterations. The use of constraints and the large number of equations enhances the uniqueness of the solution by reducing the dimensions of the solution

space.

In point-by-point methods the correct solution is selected from the infinity of possible

roots by calculating the slope of the phase curve and comparing the measured and calculated

group delays. A group delay constraint is also used as a way of determining the physical

solution.

The numerical effectiveness of the entire permeability and permittivity calculation depends on the robustness of the ODRPACK procedure and, more significantly, the robustness

of the mathematical model. For the samples used in this study, the robustness of the procedure depended on the sample. For the materials with low dielectric constant the procedure

readily determined a solution for a variety of input values with a large radius of convergence. For materials with higher dielectric constant, the procedure often converged quickly,

although the existence of alternative local optima in the mathematical model required some

testing to make sure that the converged root was the correct root.

64
3.5 Discussion

An optimization approach to the solution of the scattering equations appears to be a viable

alternative to point-by-point techniques. The technique allows a stable solution for a broad

range of frequencies. The method works particularly well for short-circuit hne measurements. Unlike the point-by-point short-circuit method which requires measurements on

two samples or in two positions, the optimized solution can obtain complex permittivity

and permeabihty on a single sample at a single position.

The reflection (^n) data are usually of lesser quality than the transmission data (52i)

for low-loss, low- permittivity materials. Therefore Su need not be included in the solution

for low-loss materials. However, reflection data Su and ^22 are very useful in determining

the position of the sample in the air hne as indicated in figure 3.7. The technique was

successful for many isotropic magnetic and relatively high dielectric constant materials.

The addition of constraints to the solution is powerful in that it further limits the possible

solution range of the system of equations and enhances the uniqueness of the solution. The

use of analytic functions for the expansion functions allows a correlation between the real

and imaginary parts of the permittivity and permeability. The results shown in figures

3.5 through 3.6 indicate that the method can be used to reduce scattering data of fairly

high dielectric constant materials. In fact, in some cases the optimized procedure yields

solutions when the point-by-point technique fails completely.

Why does an optimization approach, in many cases, reliably reduce data on higher

dielectric constant materials [e'j^ > 20), whereas point-by-point techniques generally fail?

Scattering data for higher dielectric constant materials contain responses to both primary

mode and higher order modes. As indicated in figure 3.2 for the barium titanate compound,

the optimization routine selects the primary mode data and places less weight on the higher

mode resonance data.

The optimized technique can be used to treat problems where sample lengths, sample

holder lengths, and sample positions are not known to high accuracy. Permittivity and

permeability can be found from the equations without specifying either sample position or sample length. This result could find application to high-temperature measurements.

Higher-order modes propagate in samples when two conditions are met. The frequencies

must be above cutoff in the sample, and there must be inhomogeneities or asymmetries in

the sample to excite the higher order modes. Higher-order modes can be incorporated

into this type of model by letting the optimization routine select the power in each mode.

Higher-order mode models are a subject of our current research.

65
Chapter 4

Short-Circuit Line Methods

4.1 Theory

In this section we review the mathematical formaHsm for short-circuit measurements. We

consider a measurement of the reflection coefhcient (5ii for a shorted two-port) as a function

of frequency. We begin with a mathematical analysis of the electromagnetic fields in the

sample. The details of the field model have been presented previously [26] and only the

most essential details will be presented here.

Assumptions on the electric fields in regions I, II, and III shown in figure 4.1 may be

made as follows: (1) only the dominant mode is present in the waveguide; (2) the materials

are homogeneous and isotropic; (3) only transverse electric fields are present. Under these

conditions the electric fields in these regions may be expressed as:

Ej = exp(-7o2) + 5n exp(7o2:) , (4.1)

Eii = C2 exp(-72) + C3 exp(7^) , (4.2)

Em = C4exp(-7o(z - L)) + C5exp(7o(^ - L)) . (4.3)

We wish to determine the coefficients in eqs (4.1) through (4.3) by imposing boundary

conditions on the system of equations. The boundary conditions are:

• Tangential component of the electric field is continuous at sample interfaces.

• Tangential component of the magnetic field is continuous at sample interfaces.

• The electric field is null at the short-circuit position (perfect reflect).

66
x = o X = L

Figure 4.1: A transmission line with a short-circuit termination.

Expressions for the coefficients in eqs (4.1) through (4.3) are presented in reference [26].

Matching boundary conditions of the field equations at the interface and the reflect

yields an equation for the permittivity and permeability in terms of the reflection coefficient,

P — S\i = Ci- With the sample end face located a distance AL from the short,

Su

where

and

-2I3S -\- [{6 + I) -\- {S - 1)13^] tanh7L

2(3-\-[{S+l)-{S- l)/?2]tanh7L

(5 = exp(-27,AL) .

In terms of hyperbolic functions

tanh7L + /3ta.nh-foAL - (3{l + /5tanh7L tanh 70AL]

^11 =

tanh7L-f /?tanh7oAL + /9(1 + ^ tanh 7L tanh7oAL)

(4.4)

(4.5)

(4.6)

(4.7)

5ii for a matched two-port can be obtained as a special case from eq (4.4) by letting (5 —> 0. Although in the derivation of eq (4.7) it is assumed that the sample plane coincides

with the measurement calibration plane, this is not the case in general; however, we can

transform the reference plane position by a simple procedure. To accomplish this, we write

the most general expression for the reflection coefficient as

67
Sll{tTans) — RlSu , (4.8)

where Smrans is the reflection coefficient at the cahbration reference plane position,

7?i =exp(-7„Li) ,

- (4.9)

and Li is the distance from the cahbration plane to the sample front face. Equation (4.8)

transforms the reflection coefficient from the calibration plane to the plane of sample front

face. It is of interest in many applications to eliminate the distance Li from eq (4.8). This can be accomplished by measuring ^n of the empty sample holder,

Sn(eTnpty) = -exp(-27o[Li + AL + L]) = - exp(-27oZ/air) , (4-10)

and therefore the ratio of the filled to empty holder reflection coefficient is

Suitrans) ^ _ exp(27o [AL + L] )5h (4.11)

'^ll(empty)

If both the permeability and the permittivity are required, measurement data for two

different short-circuit positions are needed. Note that standing waves can be formed in the

region between the sample and short-circuit and between the calibration plane and sample

front-face. Therefore certain frequencies, depending on sample length and the other lengths,

will give better results for permittivity and other frequencies better results for permeability.

The position of the short-circuit is a low electric field and high magnetic field region and

a position A/4 from the short- circuit is a high electric field and low magnetic field region.

Therefore as frequency permits, for permittivity measurements the sample should be moved

away from the short-circuit termination. Permeability in isolation can be obtained with

the sample at the short-circuit position. Of course when an ANA is used measurements

will be taken at many combinations of field strengths and therefore the uncertainty will

vary with frequency.

4.1.1 Two Samples of Different Lengths

It is possible to solve for the permeability and permittivity when the scattering parameters

with samples of two differing lengths are measured. To see this, let us consider two samples, one of length L and one of length aL as indicated in figure 4.2. Then for independent measurements on the two samples we have

^11(1) = ffj^ > (4-12)

and

68
Port1

Short-Circuit

Sample 2

Figure 4.2: A transmission line with a short-circuit termination in two sample magnetic

measurements.

Z = exp(-7L) . (4.13)

The reflection coefficient T is given by eq (2.37). The scattering 1-port parameter is given

by

(4.14)

'5'll(2) - ^ _ p^2a

Therefore we can solve for Z in eq (4.12),

Z'

and substitute it into eq (4.14) to obtain

5'ii(i)r — 1 (4.15)

r- S'ii(i)-r

5ii(i)r-i_

a

'->11(2) — 1-r 5ii(i)-r

5ii(i)r-i

a (4.16)

Equation (4.16) is solved iteratively for T and then Z is found from eq (4.15). The permittivity and permeability can then be obtained, if we define

1A2

1 ^ ^1^2

ln(^)]^ , '1-kL 'Z

69

(4.17)
Port1

Sample in Position 2

Sample in

Position 1

Figure 4.3: A transmission line with a short-circuit termination for two-position measurements.

-R

A2

f^R

. 1a!

1 + r

1

.

A2-

(i-r)Ayj^'

(4.18)

(4.19)

where Aq is the free-space wavelength and A^ is the cutoff wavelength. Equation (4.17)

has an infinite number of roots for magnetic materials because the logarithm of a complex

number is multi-valued. In order to pick out the correct root it is necessary to compare

the measured to the calculated group delay.

4.1.2 Single Sample at Two Short-Circuit Positions

It is possible to obtain an explicit solution to eq (4.4) when measurements at two different

short-circuit positions are taken. The explicit solution is obtained by solving eq (4.4) at a given short-circuit position (position 1) for tanh7L and then substituting this expression

into eq (4.4) at another short-circuit position (position 2) as indicated in figure 4.3. For two different short-circuit positions at the same frequency we obtain pi and p2 for

positions 1 and 2:

70
^' -2/? + Pj _ 1 )/92 _ (^^ ^ 1 )] tanh 7L ' ^ ^ ^

^' -2/? + [(<52-l)/?2_(^^ + l)|tanh7L ' ^ ^^ where ^1,^2 denote the phases calculated from eq (4.6) for ALi,AL2 respectively. These

equations yield

tanh7i=-^ 2Wi+ ^l) ,

2/3(6, +p

where n is an integer. Since the arctangent is multi- valued, the correct value of n is deternnined from the group delay arguments given in section 3.2.1. Also

2 ^ Si{S2{pi - P2) + Pi/>2 + 1 - 2/72) - {S2{pi{p2 - 2) + 1) + /)2 - Pi) .^24)

Sl{S2{pi-p2)^PlP2^l+'2p2)-{S2{p,{p2^2)^l) + P2-pi) '

Once /3 is known, eqs(4.22) and (4.24) can be used to find permittivity and permeability.

4.2 Measurements

In the SCL technique the scattering parameter Sn is measured broadband, with the sample

at a given position in the sample holder. The distance from the sample to the short-circuit

termination must be known to a high degree of accuracy. If both permeability and the

permittivity are required then the sample must be moved in the line and the S-parameters

again taken.

Depending on the position of the short-circuit, the sample may be immersed in either a

region of high electric field or high magnetic field. A strong electric field is advantageous for

permittivity determination, whereas a strong magnetic field is advantageous for permeability determination. Generally, the sample end will be in a region of high magnetic field when

the sample is in closest physical contact with the short. It is possible to take advantage

of the fluctuating electric and magnetic field distributions when performing permittivity

and permeability measurements. When taking broadband measurements on an ANA it is possible to predict when the sample is immersed in the various field strengths. Then one

can select the measurements to be used for permittivity and permeability calculations [35].

Measurements were made on an ANA for various samples. Using eq (4.4) we obtain

the permittivity and permeability which are shown in figures 4.4 and 4.5.

71
10 11 12 FREQUENCY (GHz)

13

Figure 4.4: e^ without gap correction using SCL for a loaded polymer mixture.

Mr

1.4

1 1 1

1.2 -—

1.0 —

0.8 - -

0.6 - -

0.4 112 SCL -

0.2 - -

n 1 1 1 <) 10 11 12 i: FREQUENCY (GHz)

Figure 4.5: /i^ without gap correction using SCL for a loaded polymer mixture.

72
4.3 Uncertainty of Short- Circuit Line Measurements

For magnetic materials it is necessary to make two independent measurements at a given

frequency. Independent measurements can be obtained either by measuring samples of two

different lengths or by taking measurements of a given sample at two locations in the line.

The special case of measurement of two samples of varying lengths (Li, L2) can be obtained

from the solution presented below by substituting L ^ Li, 70 —^ 7 and AL —> L2 — Li.

The uncertainty involved in two position measurements is explored in this section. The

uncertainty incurred with the equations expounded in this report is estimated. The sources

of uncertainty in the SCL measurement include

• Uncertainties in measurements of the magnitude and phase of the scattering parameters.

• Gaps between the sample and sample holder.

• Dimensional variations in the sample holder.

• Uncertainty in sample length.

• Short-circuit and line losses and connector mismatch.

• Uncertainty in positions of the reference plane and sample in holder.

Adjustment for errors due to gaps around the sample is obtained from equations available in the literature [33,34,35]. The formulas given in the literature generally under-correct

for the real part of the permittivity and over-correct for the imaginary part of the permittivity. All measurements of permittivity are assumed to have been corrected for air gaps

around the sample before the following uncertainty analysis is applied.

In order to evaluate the uncertainty introduced by the measured scattering parameters,

we assume that a differential uncertainty analysis is applicable. This assumption implies

that uncertainties are of small enough magnitude so that a local Taylor series can be applied.

We assume that a Taylor series approximates deviations of the function from a given point.

We assume that the worst case uncertainty due to the S-parameters and sample lengths can be written

Ae'„ 1 -R R ^^'^,A,5.,Uf|^A..Uf|^ALUf^A.l , (4.25) -^ww^r^''^) ^\de,-'') ^\dL-) -\dd

Ae'l, 1

t"R

9^'k.,.y.(d^'k..Y.(d^'k^Xj9^^,]\

(,.26)

73
where A^ is the uncertainty in the phase of the scattering parameter, A|5n| is the uncertainty in the magnitude of the scattering parameter, Ad is the uncertainty in the air

gap around the sample, and AL is the uncertainty in the sample length. The gap correc- tion uncertainty is given in [26]. The uncertainties used for the S-parameters depend on

the specific ANA used for the measurements. In general uncertainties due to flange bolt

torquing and connector repeatability must be added to these uncertainties.

Let us define the variables

a = tanh7L, (4.27)

& = tanh7oAL. (4.28)

We wish to obtain explicit relations for the derivatives of e^ and ^|^ with respect to independent variables |5ii(,)| and ^,-, i == 1,2. We define ^i^i) as the reflection at short position

1 and 5ii(2) as the reflection at short position 2. Next define

tanh -fL + /3 tanh jpAL - jSjl -\- (3 tanh 7L tanh 70AL) _ ^ _ ^

, ^o)

tanh7L-|-/3tanh7oAL + ^(l +/5tanh7Ltanh7oAL) n " • { ) We assume the following as independent variables Sii{i), i — 1, 2, L, AZ, and d. Derivatives

of eq (4.29) with respect to the independent variables can be found analytically. By the

chain rule we have

rdf da d-y df df3 dfi'j^ rdf da d^ df d/3 oe)^ i-a \

a-l 0.2

This equation is evaluated at position one.

dfda ^7 df dp df,-k dfda d^ df df3 de^j,

^dad^dfi% "^ dpdfi*j,^ d\Sn{m)\ ^a df de}^ ^ d/3 de*j,^ d\Su{m)\ ^""P^-^^"^^ " ^^"'''^

blm b2Tn

This equation is evaluated at position two. The four derivatives can be written

ai oie^"" I + «2 o|f'' I - exp{jern)S,,m , (4.32)

0\Oii(m)\ 0'Pll(m)|

^l oig ^

I + ^2 f)\Q ^I = ^Mj^m)S2,m , (4.33)

where 6i^rn is the Kronecker delta function.

74
or

At the first position the derivatives with respect to length are

^dad-fdfi*j,'^ d^dfi^j,^"^ dL ^dad^de),^ d(3de),^^ dL '^ dadL~ ' ^ ^ ^

dfj.*p, de*r, df da ^ , ^

and for the second short-circuit termination position

, du'fi , de*F> df da

The derivatives with respect to the distance from the sample back face to the short-circuit

termination can similarly be calculated

.dldad^dld£ df^.^dldad^dld^

^ , ^_^^n (A '^7) ^dad-fdfi'j^ dl3dfi}^^^dM ^dad-fde}, 30 de'j^^"" dAL ^ db dAL ' ^

'

and for the second short-circuit position

The derivatives can be calculated explicitly to yield

da

1 - ;52tanh7oAL

-)-

tanh7L + /? tanh 70AL + ^5(1 + ^tanh7L tanh7oAL

tanh ')L -\- 13 tanh ')qAL + /?(1 + /5 tanh 7Z/ tanh 70L)

(tanh7L + ^tanh7oAL + ^5(1 +/? tanh 7L tanh 7oAL)2

(1 -f^'tanh7oAL) , (4.40)

dh

-fiH&nh -fL )-

tanh7L + /? tanh7oAL + /?(! -f /? tanh 7L tanh 70 AL

tanh7L + /?tanh7oAL + l3{l + /5 tanh 7L tanh 70!)

(tanh 7L + /? tanh 7oAL - /?(! -f /? tanh 7L tanh 70 AL)^

(/?-f ^Hanh7L) , (4.41)

75
tanh7oAL — 1 — 2^ tanh 7L tanh 7oi/

tanh7L + /?tanh7oAL -\- f3{l -\- /?tanh7L tanh7oAL

tsLiih^L + /?tanh7oAL + /?(1 + l3ta.nh-yL tanh7oL)

X

(tanh7L + ^tanh7oAL - /?(1 + ^9 tanh7L tanh7oAL)2

(tanh 70AL + 1+2/? tanh 7L tanh 70 AL) . (4.42)

The following derivatives will be needed in the forthcoming analysis

^ = Lsech^7L , (4.43)

(77 —d^ = 7sech^7L , (4.44)

Ol-/

db ,2

dAL

8(3 1

7osech 7oAL , (4.45)

4.46)

df3 dj I 7 1

d£__ d2__l 7__1_

dl3 _ dj_dl

d^_ di 1

^7 ^VcL^fl

a//^ 27

dt*R 27

76

(4.47)

[4.48)

[4.49)

(4.50)

[4.51)

[4.52)
0.16

1 1 1 1

0.12 i

^^R 1 1 as,, 0.08 \ /

'

0.04

n 1 1— —i. 1 0.1 0.2 0.3 0.4 0.5

Figure 4.6: A plot of the derivative of e'^ with respect to |5ii| as a function of the distance

from the short-circuit termination for the case of thin samples L <C A^.

df4

d\s

j|'5'ii(j)

62

df^h

d\Su(r)

11

den

d\Snii)

d\Su(2)\

d(-R

exp(j^i)6i

02^1 — 0162

ai^|'S'ii(2)|

exp(j6'2)ai

(4.53)

(4.54)

(4.55)

(4.56)

(4.57)

^|5'll(2)| «2&1 - ^162

In figures 4.6 and 4.7 dt*jil d\S\\\ is plotted as a function of the distance from the short

circuit termination. We see that in the case of electrically thin samples the minimum

uncertainty occurs when the second sample measurement is at Ao/4 from the short-circuit

termination. As shown in figure 4.7, this is not the case for samples that are not electrically

thin.

77
Figure 4.7: A plot of the derivative of e'^ with respect to |5ii| as a function of the distance

from the short-circuit termination for the case of long samples L ^ X-m78
Chapter 5

Discussion

We have presented an overview of the theory and the measurements for both the TR

method and the SCL method for dielectric and magnetic materials. In addition, relevant

uncertainty analyses have been developed. Equations were presented for TR and SCL

measurements that are reference plane invariant. These equations can be solved either by

iteration or explicitly.

There are two common problems in the data reduction techniques for transmission

hne measurements. These problems are the existence of higher mode resonances and the

problem that for low- loss materials the solution of the equations become ill-conditioned at integral multiples of one-half wavelength in the sample. The one-half wavelength instability occurs because the phase of Su contains a large uncertainty when |5ii| —> and also

because the equations in this limit yield only the phase velocity. For dielectric and magnetic

measurements, the uncertainty is a function of the sample length. In general for lowloss materials, samples long in relation to wavelength give more accurate results, however

overmodes may be produced. Thus for broadband measurements of low-loss low-dielectric

materials, it is preferable to use longer samples. However, with lossy materials, very long

samples result in only front face reflection information and the results have a relatively large

uncertainty. For relatively lossy materials, sample lengths on the order of one attenuation

length are optimum. Longer samples allow the propagation of higher order modes, and

therefore higher mode responses will be contained in measured scattering data. However,

the uncertainty in the spectrum between over-moded resonance regions will be lower for

longer samples than for electrically short samples.

For thicker samples the problem is more complicated and a knowledge of the uncertainty

analysis is important for interpreting the results. For permeability measurements in a shortcircuit line the numerical reduction procedure becomes divergent when there is an integral

multiple of one-half wavelength in the sample. For TR and SCL measurements we can summarize our conclusions as follows

79
• For SCL measurements, length of samples should be such that l^n] > —20 dB.

• The one-sample techniques appear to allow a better reduction of the scattering data.

• The optimized solution and the reference plane invariant solution appear to be the

most accurate and stable methods. However, local minima have to be avoided.

• The short-circuit line is a simple way to obtain permittivity and permeability simultaneously. For low-loss materials it does suffer from numerical instabilities.

• In SCL permittivity measurements, minimum uncertainty is obtained when the sample is A/4 from the short. This is only true for electrically thin samples.

• Minimum uncertainty magnetic measurements can be made for single frequency measurements by taking one measurement at a maximum electric field position and an- other measurement at a maximum magnetic field position.

The various TR and SCL techniques are compared in table 5.1 and 5.2.

80
Table 5.1: Dielectric and magnetic TR measurement techniques compared.

Technique Applicability Strong Points Weak Points

Full 2-port, one sample Dielectric Requires one sample

Full 2-port, one sample Magnetic Requires one sample Unstable at nAm/2

NRW Technique Magnetic Simple Solution Unstable at nXml'^

Two-Position Technique Magnetic Doesn't exist for TR techniques Two sample technique Magnetic Stable over all frequencies Requires two samples

Multi-point techniques Magnetic Very stable non-global minima

Table 5.2: Dielectric and magnetic SCL measurement techniques compared.

Technique Applicability Strong Points Weak Points

One-position Technique Dielectric Stable —

Two-Position Technique Magnetic Requires only one sample Unstable at nAm/2

Two sample technique Magnetic Stable over all frequencies Requires two samples

Multi-point techniques Magnetic Very stable Alternative minima

81
We would like to thank Eleanor Livingston for her excellent comments and suggestions, Paul Domich for collaboration on the optimized solution, Chris Jones for help with

measurements and plotting, William Kissick, Claude Weil, and Howard Bussey for various

discussions and as sources of encouragement. Finally we would like to thank Ramon Jesch,

Bill Daywitt, and Lydell Frasch for reviewing the manuscript. ,

82
References

[1] Belhadj-Tahar, N.; Fourrier-Lamer, A.; and de Chanterac, H., "Broadband simulta- neous measurement of complex permittivity and permeability using a coaxial discontinuity," IEEE Trans. Microwave Theory Tech., vol. 2, pp. 1-7, January 1990.

[2] Deschamps, G. A., "Determination of reflection coefficients and insertion loss of a wave-guide junction," J. Appl. Phys., vol. 2, pp. 1046-1050, August 1953.

[3] "Measuring dielectric constant with the HP 8510 network analyzer," Product note no. 8510-3, Hewlett Packard, 1985.

[4] Freeman, M. S.; Nottenburg, R. N.; and DuBow, J. B., "An automated frequency

domain technique for dielectric spectroscopy of materials," J. Phys. E: Sci. Instrum.,

vol. 12, pp. 899-903, 1979.

[5] Nicolson, A. M.; and Ross, G. F., "Measurement of the intrinsic properties of materials

by time domain techniques," IEEE Trans. Instrum. Meas., vol. IM-19, pp. 377-382,

November 1970.

[6] Weir, W. B., "Automatic measurement of complex dielectric constant and permeability

at microwave frequencies," Proc. IEEE, vol. 62, pp. 33-36, January 1974.

[7] Stuchly, S.; and Matuszewski, M., "A combined total reflection transmission method

in apphcation to dielectric spectroscopy," IEEE Trans. Instrum. Meas., vol. IM-27,

pp. 285-288, September 1978.

[8] Delecki, Z. A.; and Stuchly, S. S., "Uncertainties in the determination of dielectric

properties in the infinite sample method," IEEE Trans. Instrum. Meas., vol. 38,

pp. 902-906, June 1989.

[9] Franceschetti, G., "A complete analysis of the reflection and transmission methods for

measuring the complex permeabiHty and permittivity of materials at microwave," Alta

Frequenzia, vol. 36, pp. 757-764, August 1967.

83
[10] Ligthardt, L. P., "A fast computational technique for accurate permittivity determination using transmission line methods," IEEE Trans. Microwave Theory Tech..,

vol. MTT-31, pp. 249-254, March 1983.

[11] Roberts, S.; and von Hippel, A., "A new method for measuring dielectric constant and

loss in the range of centimeter waves," J. Appl. Phys., vol. 7, pp. 610-616, July 1946.

[12] Brydon, G. M.; and Hepplestone, D. J., "Microwave measurements of permittivity

and tan^ over the temperature range 20 — 700°C," Proc. Inst. Elec. Eng., vol. 112,

pp. 421-425, 1965.

[13] Szendrenyi, B. B.; Kazi, K.; and Mojzes, I., "An alternative broadband method for

automatic measurement of the complex permeabihty and permittivity of materials at microwave frequencies," Microwave Theory and Tech. Digest, vol. 1, no. 1, pp. 743-745,

1988.

[14] Mattar, K. E.; and Brodwin, M. E., "A variable frequency method for wide-band

microwave material characterization," IEEE Trans. lustrum. Meas., vol. 39, pp. 609-

614, August 1990.

[15] Maze, G.; Bonnefoy, J. L.; and Kamarei, M., "Microwave measurement of the dielectric

constant using a sliding short-circuited waveguide method," Microwave J., pp. 77- -88,

October 1990.

[16] Dakin, T. W.; and Works, C. N., "Microwave dielectric measurements," J. Appl. Phys.,

vol. 18, pp. 789-796, September 1947.

[17] Bowie, D. M. and Kelleher, K. S., "Rapid measurement of dielectric constant and loss

tangent," lEE Trans. Microwave Theory Tech., vol. MTT-4, pp. 137-140, 1956.

[18] Sequeira, H. B., "Extracting //^ and tj. from one- port phasor network analyzer measurements," IEEE Trans. lustrum. Meas., vol. 39, pp. 621-627, August 1990.

[19] Chao, S., "An uncertainty analysis for the measurement of microwave conductivity

and dielectric constant by the short- circuited line method," IEEE Trans. lustrum.

Meas., vol. IM-35, pp. 36-41, March 1986.

[20] Miles, P. A., Westphal, W. P. and von Hippel, A., "Dielectric spectroscopy of ferromagnetic semiconductors," Rev. Mod. Phys., vol. 29, no. 3, pp. 279-307, 1957.

[21] Solymar, L. and Walsh, D., Lectures on Electrical Properties of Materials. New York:

Oxford University Press, 1988.

84
[22] Ramo, S., Whinnery, J. R. and Duzer, T. V., Fields and Waves in Communication

Electronics. New York: John Wiley and Sons, 1984.

[23] Gilbert, T. A., "Constitutive relations," Tech. Rep., Armour Research Foundation,

1955.

[24] Hord, W. E., Ferrite Control Components, pp. 238-272. John Wiley and Sons, 1989.

[25] Soohoo, R. F, Theory and Application of Ferrites. Englewood Cliffs, NJ: Prentice-Hall,

Inc., 1960.

[26] Baker- Jarvis, J., "Transmission/reflection and short-circuit line permittivity measurements," Nat. Inst. Stand. Tech. (U.S.), Tech. Note 1341, 1990.

[27] Baker- Jarvis, J., Vanzura, E. and Kissick, W. "Improved technique for determining complex permittivity with the transmission/reflection method," IEEE Trans. Microwave Theory Tech., vol. 38, pp. 1096-1103, August 1990.

[28] Kerns, D. M. and Beatty, R. W., Basic Theory of Waveguide Junctions and Introductory Microwave Network Analysis. New York: Pergamon Press, 1967.

[29] Rasmussen, A. L., Enfield, A. W. and Hess, A. "Advances in the design and application

of the radiofrequency permeameter," Nat. Bur. Stand. (U.S.) J. Res., vol. 56, no. 5,

pp. 261-267, 1960.

[30] Hoer, C. A. and Rasmussen, A. L., "Equations for the radiofrequency magnetic permeameter," Nat. Bur. Stand. (U.S.) J. Res., vol. 67C, no. 1, pp. 69-76, 1963.

[31] Powell, R. C. and Rasmussen, A. L., "A radio-frequency permittimeter," IRE Trans.

Instrum., vol. 1-9, no. 2, pp. 179-184, 1960.

[32] Goldfarb, R. B. and Bussey, H. E., "Method for measuring complex permeability at radio frequencies," Rev. Sci. Instrum., vol. 58, no. 4, pp. 624-627, 1987.

[33] Bussey, H. E., "Measurement of rf properties of materials a survey," Proc. IEEE,

vol. 55, pp. 1046-1053, June 1967.

[34] Westphal, W. P., "Techniques of measuring the permittivity and permeabihty of hquids

and solids in the frequency range 3 c/s to 50 kmc/s," Laboratory for Insulation Res.

Tech. Report, MIT, 1950.

[35] Bussey, H. E.; and Gray, J. E., "Measurement and standardization of dielectric samples," IRE Trans. lustrum., vol. I-ll, no. 3, pp. 162-165, 1962.

85
[36] Domich, P. A.; Baker-Jarvis J.; and Geyer, R. G., "Optimization techniques for permittivity and permeability determination," J. Res.,vol. 96, no. 5, pp. 565-575, 1991.

[37] Wallin, S. R., Dielectric Properties of Heterogeneous Media. PhD thesis, University of

Wyoming, 1985.

[38] Papoulis, A., Signal Analysis. New York: McGraw-Hill, 1977.

[39] Boggs, P. T.; Donaldson, J. R; Byrd, R. H.; and Schnabel, R. B., "ALGORITHM

676 ORDPACK: software for weighted orthogonal distance regression," ACM Trans.

Math. Software, vol. 15, no. 4, pp. 348-364, 1989.

[40] Marcuvitz, N., Waveguide Handbook. New York: Dover Publications, 1951.

[41] Champlin, K. S.; and Glover, G. H, "Influence of waveguide contact on measured

complex permittivity of semiconductors," J. Appl. Phys., vol. 37, pp. 2355-2360, May

1966.

[42] Schwinger, J. and Saxon, D. S., Discontinuities in Waveguides. New York: Gordon

and Breach Science Publishers, 1968.

[43] Jesch, R. L., "Dielectric measurements of five different soil textural types as functions

of frequency and moisture content," Nat. Bur. Stand. (U.S.) Tech. Rep., 1978.

[44] Jesch, R. L., "Dielectric measurements of oil shale as functions of temperature and

frequency," IEEE Trans. Geosci. Remote Sensing, vol. GE-22, pp. 99-105, March 1984.

[45] Ben-Menahem, A. and Singh, S. J., Seismic Waves and Sources. New York: SpringerVerlag, 1981.

86
Chapter 7

Appendices

Appendix A

Magnetism in Matter

A.l Description of Magnetic Phenomena

The origin of magnetism is related to the electrostatic coulomb repulsion between electrons

and is intimately related to the spin and orbital angular momentum of electrons, nuclei,

and other charged particles. Stern and Gerlach proved the existence of discrete magnetic

moments by observing the deflection of silver atoms passing through a spatially varying

magnetic field. The quantum mechanical relation between magnetic moment and angular

momentum of an electron is rhj = —g^J, where g is the Lande g-factor « 2.002319114, 2m

eh 2m 9.2742 X 10~^^ {J — m^/W) is the Bohr magneton, J is the total quantum mechanical

angular momentum, e is the electronic charge, h is Planck's constant, and m is the mass. The gyromagnetic ratio is defined as magnetic dipole moment e / * -, \ 7 r=

=g . (A.l)

angular momentum 2m

There is a diversity of magnetic phenomena due to the existence of various couphng

schemes of angular momenta quanta. Types of magnetism include paramagnetism, which

is due to spin and angular momentum of individual electrons, diamagnetism which has its origin in the orbital angular momentum of the electron, a,nd ferromagnetism originates from

the formation of domains with each domain containing a large number of ahgned spins.

A. 1.1 Field Description of Electromagnetic Phenomena

It has been found that dielectric and magnetic phenomena are adequately described by a

set of field vectors. These vectors represent the electric field, E and magnetic field, H, the

displacement field, D, the induction field, B, the polarization field, P, and magnetization

87
field M. Maxwell's equations define the spatial and temporal evolution of these field vectors.

Constitutive relations between field quantities and material properties are necessary to

describe electromagnetic phenomena. The displacement field is related to the electric field

hy D =e -E, where e is the permittivity tensor. For linear materials the permittivity does

not depend on the field strength. The permittivity is a measure of electronic, ionic, and

dipolar polarization. The permittivity is frequency dependent with dipolar polarization

occurring below 10^° Hz, ionic polarization below 10^^ Hz, and electronic polarization

above 10^"^ Hz. The permittivity of free space is cq = 8.85419 F/m. The magnetic field

is related to the induction field by B —fJ- H, where /^ is the permeability tensor. The

permeability of a material is a measure of the degree to which it allows the penetration by

an external magnetic field. The permeability of free space is /^o = 47r x \0~^ H/m. The

permittivity and permeability of free space are related to the speed of light in vacuum

c = 1/ y/eofiQ = 2.99792458 x lO^m/sec. The electric field may contain sources, so that

V E = p, where p is the free charge density. Induction fields are sourceless, expressed

mathematically by V • jB = 0. For a charge e moving with velocity v through and electric

field E and magnetic field B, the Lorentz force on the charge is F = e[E -}- {v x B]].

The electronic properties of isotropic substances can be described macroscopically by

scalar miaterial properties in terms of the relative complex permittivity and permeability,

t*ji and ^^:

e = e' - je" = (e^ - je'^)eo = eflCo , (A.2)

H^ p,' - jp" = [ppt - J>'A)/^o = f^*Rf^o (A.3)

The electric and magnetic fields are modified by the presence of matter in the spacetime region in and around the body. The presence of magnetism in matter is described

by the magnetization vector M which quantifies the number of magnetic dipoles per unit

volume. The magnetic field H is related to the induction and magnetic moment vectors by

H = —B - M , (A.4)

where po is the permeability of free space. The magnetic field is also related to the magnetization field, M by a constitutive relationship in terms of the susceptibility Xm- For a

linear medium the relation is M = XmH . (A.5)

The permeability and susceptibility are related through eqs (A.4) and (A.5)

p = po{l + Xm) • (A.6)

88
Similarly, the electric field is related to the displacement field (D) and the polarization field

(P)by

D^eoE + P = eE. (A.7)

It is useful to define a constitutive relation between polarization and electric fields using

the electric susceptibility {Xei)

P = CoXeiE. (A.8)

and therefore

e = eo(l+Xe/)- (A.9)

A. 1.2 Types of Magnetism

Magnetic materials are classified by the values of permeability.

• diamagnetic fi' < ^o,

• paramagnetic fx' > ^o,

• ferromagnetic //' >> /iq.

The susceptibilities of the various classes of magnetic phenomena are

• diamagnetic Xm < 0,

• paramagnetic Xm > 0,

• ferromagnetic Xm >> 0. Due to the complicated quantum-mechanical origin of magnetism with various competing

effects, it is not always possible to classify a material into one of the categories. For example,

a ferrite may be diamagnetic in X-band and paramagnetic at lower frequencies.

A. 2 Paramagnetism

Paramagnetism arises from the alignment of individual spins and angular momentum of

particles in external magnetic fields. Paramagnetism is an interaction between the tendency

for the electron spins to be aligned with the field on the one hand and thermal agitation

which tends to randomize the spins on the other hand. Paramagnetic phenomena in in- sulators is temperature dependent and follows Curie's law. In metals paramagnetism is strongly influenced by the conduction electrons and has minimal temperature dependence.

Paramagnetic materials are primarily the rare earth and transition ions with incomplete

atomic shells. There are two types of paramagnetism

89
• Spin paramagnetism,

• Orbital paramagnetism.

Spin paramagnetism is due to alignment of electron spins and is only slightly temperature dependent. Spin paramagnetism or Pauli paramagnetisrri occurs in metals. Orbital

paramagnetism, caused by alignment of orbital magnetic moments, is strongly temperature

dependent. This type of paramagnetism occurs in insulators.

A. 2.1 Diamagnetism

Diamagnetism is the magnetic effect that is due to orbital angular momentum effects. Larmor diamagnetism occurs in filled-shell insulators. The origin of diamagnetism in materials

is the orbital angular momentum of the electrons in applied fields. Diamagnetic materials

have a negative susceptibility and generally it is not sensitive to temperature variations

at least for nonsuperconducting materials. Diamagnetic materials do not have a strong

magnetic response.

A. 2. 2 Ferromagnetism

In ferromagnetic materials, spin coupling allows regions of aligned spins to be formed, called

domains. In each domain the spins are more or less aligned. However, adjacent domains as a whole may be arranged in a random fashion. As a magnetic field is applied the domains

more or less align with the field. The difference between paramagnetism and ferromagnetism is that in the case of the

former, spins interact minimally, whereas in ferromagnetism the spins strongly interact

to cause alignment. Ferromagnetic materials can exist in a nonmagnetized state since

magnetic energy is at a minimum when the domains are randomly situated or in a state of

maximum entropy. This random arrangement of domains is possible because it is found in a

detailed analysis that it is energetically more favorable for many ferromagnetic materials to be magnetically neutral. The boundaries between the oriented spin regions, called domain

walls require energy for formation. There is a detailed balance between the magnetic

field energy caused by alignment of spins in a domain on the one hand and the energy

required for domain wall formation on the other hand. Dipolar energy is decreased by wall

formation, but exchange energy is increased by the Pauli exclusion principle. Domain walls are normally of 0.01 - 10 // m thick and can deform under applied fields or mechanical

stresses. As the temperature increases in a ferromagnetic material the kinetic energy can over- whelm the magnetic energy and the preferential alignment of spins can be broken. The

90
O ^ O f o

o ^ o ^ o 4

^ o # o 6 o

O A^-O-^ o 6

o ^ o ^ o 4 o

^ o 4 o

^ o ^

(a) (b)

Figure A.l: Lattice structure of antiferromagnetic and ferrite materials. In one lattice the

spins are up and in the other lattice the spins are down.

temperature where the kinetic thermal energy becomes predominant over the magnetic

energy is called the Curie temperature.

If a ferromagnetic material is immersed in an increasing external field its magnetization

increases. However as the external field is removed the magnetization curve does not

necessarily follow the reverse curve back to the initial state; rather a slightly different curve

is followed.

This phenomenon is called hysteresis and is caused by the irreversible movement of the

domain walls. The irreversibility is caused by defects in the domain wall lattice. Thus, fer- romagnetic behavior depends on the past history of the sample and is not totally reversible.

Materials with large hysteresis effects are called hai^d and materials with small hysteresis

effects are called soft.

A. 2. 3 Ferrites and Antiferromagnetism

Antiferromagnetism is a property possessed by many transition elements and some metals.

In these materials the atoms form an ordered array with alternating spin moments so as

to give zero for the net magnetic moment of the material. Antiferromagnetic materials are composed of two interpenetrating lattices. Each lattice has all spins more or less aligned,

but the lattices as a whole are inverse structures as indicated in figure A.l.

Antiferromagnetic materials do not generally support permanent magnetization and do

91
not have a strong magnetic response to an applied field. Ferrite materials also consist of two

overlapping lattices whose spins are oppositely directed, but with a larger magnetic moment

in one lattice than the other. Since spin angular momentum is not canceled totally between

the lattices these materials have a magnetic response to an applied field. The magnetic

response increases with temperature. Antiferromagnetic materials are paramagnetic above

the Neel temperature.

Most ferrites are mixtures of oxides such as magnetite of the generic form X0.Fe20:i

where A' is a divalent metallic ion such as Ft (magnetite), Ni (nickel ferrite), Cu (copper

ferrite), Mg (manganese ferrite), Co (cobalt ferrite), or Li (lithium ferrite). There are

also many other spinel class ferrites that contain additional components, for example, zinc

and aluminum. The spinels have either a normal or an inverse lattice structure formed

by controlled quenching. Many ferrites have few free electrons and therefore are useful in

microwave frequency components since the absence of free electrons prevents eddy-current

losses that occur in conducting materials at high frequencies.

A.3 Equations of Motion for the Magnetization Vector

A. 3.1 The Torque Equation

In this section we will develop macroscopic equations of motion that underlie the coupling

of angular momentum and magnetic fields. As a model [22], we consider a spinning particle

exhibiting angular momentum J, immersed in a magnetic field. The presence of spin induces

a magnetic moment m = 7c, J, where the 7^ is the gyromagnetic ratio. The magnetic field

will interact with the angular momentum by inducing a torque f

f=-. (A.10)

We define a magnetic moment, m = 7^ J. The equation of motion of the spin system in an

applied field B is T = mxB, (A.ll)

or ^ = 7,(mxJ5). (A.12)

If we average over a significant number of these particles we obtain a macroscopic magnetic

moment

92
M = Nrh , (A.13)

where N = number of magnetic moments. We obtain the following equation of motion for

the magnetization —dM , - -, A . , . = 7,/x(Mx//) + - , (A.14)

where A is a constant and r is a characteristic relaxation time of the system. As we will see in the next section it is possible to obtain a tensorial constitutive relation between B. and B by use of eq (A.14) for a magnetic moment in an external magnetic field. This

constitutive relationship will define the susceptibility and permeability tensors.

In real materials there always exists some dissipation due to damping. Dissipation

is caused by such effects as magnetic dipole radiation and magneto-elastic coupling with

lattice phonons. The effect of a dissipation torque can be modelled as a source term added

to eq (A.14). One approximation for the dissipation torque yields the Landau equation of

motion — = nAM X //) + ^[M X (M X ti)\ , (A.15)

at \M\

damping

where a is a parameter that determines the damping. This can be reduced for small

damping to - = ,,,(M X //) + ^[M X -1 . (A.16)

This equation is due to Gilbert [23] and neglects terms nonhnear in a. The damping

introduces a nonhnearity into the problem.

A.3.2 Magnetized Magnetic Material: The Polder Matrix

The constitutive relation between the induction field and the magnetic field in ferrite materials is represented by the Polder matrix. In order to derive an expression for the Polder

permeability tensor we use as a model a magnetic dipole of moment m, [22,24], in the

presence of an external magnetic field, U[i). The net torque experienced by the dipole of

magnetic moment rrii is r = -fimdH{t)sme , (A.17)

where 9 is the angle between the dipole axis and field. In the presence of a z-axis magnetic

field the dipole will precess with a characteristic Larmor frequency of loq = \jg\Hz where

7^ is the gyromagnetic ratio. There are also non- conservative frictional forces present to

damp the rotation.

93
We will assume that in the presence of a time independent d.c. magnetic field in the

z direction (//z), the magnetization is essentially at the saturated value M^. We further

assume a combined d.c. field and an alternating field /i, we obtain li{i) — Hze3-\-h exp(ju;t),

where 63 is the unit vector in the z direction. The magnetic moment can be approximated

hy M = Mze3-\-mexp{jiot), where M^ is time independent. When we substitute these into

eq (A. 16) we obtain a system of equations

—

— ju>m = ii'jgMsies x h) -\- (a;o + jcja)(e*3 x m) -\- m x h exp (jiot) , (A. 18)

where coq = —^gH^ [25] and we assume \M\ = Mg ~ M^. In the analysis we will neglect

the last term in eq (A. 18). If

1=1

3 h -Y2 ^i^i > j=l

then by substitution into eq (A. 18)

jurux — -{iOo -^ ju)a)my - /.ifgMshy ,

(A.19)

(A.20)

(A.21)

jurriy - -ygnMshx + (ujq + ]ua)mj: ,

jiJiriz = .

Therefore

m-,

JUJfXJgMs

jujfi-fgMs

(cc'o + j'^cxY - u^2 2/

'

_ _ , _ (cuq -\- jLja)ii-fgMs

(A.22)

(A.23)

(A.24)

(A.25)

The system of eqs (A.24) through (A.25) for the linear susceptibility relation between

M and H can be expressed as

( rrix \

\m, I where

(A.26)

(A.27)
and

K =

-lolom

(cjo +iu;a)2 - u"' and

We can separate out the real and imaginary components to yield

(A.28)

(A.29)

(A.30)

(A.31)

(A.32)

'" ~ [ul - CJ2(1 + a2)]2 + 4C^2^2^2 * (^-3^)

Note that we have assumed that the magnetization is at the saturated value Mj. That is

all of the magnetic moments are assumed to be aligned with the external field. This is not

always a good assumption. For this special case the Polder matrix is

It X = K -Cj2(l+Q2)]2 + 4a;2u;2a2

'

u>Mtja[ujQ -fU''(l+a2)]

^l;-c.2(i + q2)]2 + AU^UW '

--UJMi^i^O - c^2(l+a2)]

H -u;2(l +a2)]2+4u;2w2^2 '

2u>qUJ^auM

/ X + 1 -JK \

^l= ^{j + X) = fi {AM)

V

1 /

We see that in the limit as //^

^

and M, -^ the off- diagonal components of the permeability tensor vanish and the diagonal components reduce to the isotropic permeability,

IJ-= fi /, where / is the identity matrix.

95
Appendix B

Fields in Waveguides

B.l Summary of Maxwell's Equations

Maxwell's equations are

\/ X E = -jujB , (B.r

V X H = J+jLoD , (B.2;

V-D = p, (B.3;

V-5 = 0. (B.4:

The boundary conditions at material interfaces are

72 X (^2 -^i) -0 , (B.5;

nx{H2-H,)^l , (B.e;

n-(D2-Di) = f] , (B.7]

n-{B2-Bi) = 0, (B.S;

where Jg is the surface current density and ft is the surface charge density.

96
B.2 Modes

B.2.1 TE Modes

The longitudinal coordinate is assumed to be z. If the H^ mode exists, it is the generator

of the TE mode [22]. The TE modes satisfy the boundary value problem,

{V?. + ^^,2}//,(T£;) = 0, (B.9)

where k\ — \z^\ 7^ are real, positive eigenvalues. The boundary conditions are

n (V//(,(7£;)))|o„ conductor = i (B.IO)

where n is the normal vector. Also

^3(T£) = 0. (B.U)

The other field components are then

1_ Ht(te) = -^^^tHz , (B.12)

^T(TE) = -Zte{z X Ht(te)) . (B.13)

B.2. 2 TM Modes

For the case of TM waves

{V?^ + ^^n^.(TM) = 0, (B.14)

where A;^ = /c'^ + 7^ are real, positive eigenvalues determined by boundary infornnation on

the waveguide.

The boundary conditions are

E(z(TM))\on conductor = , (B.15)

Hz^TM) = . (B.16)

The other field components are

ETiTM) = --^,'^TEz, (B.17)

Ht{tm) = ^—(^x Et(tm)) . (B.18)

l = j^k^-k^. (B.19)

97
B.2.3 TEM Modes

The propagation of TEM modes are possible in addition to the TE and TM modes in

coaxial cable . The cutoff wave numbers for higher order TM waves in coaxial line are given by the roots of:

[22], and for TE waves in coaxial line by

N'^iKR,) N'^ihR,

J'^iKR,) J'r,{KR,\ = 0, (B.21)

where J and A^ denote the Bessel functions of the first and second kind and Ri and R^ are the inner and outer radii respectively. The cutoff wavelengths are given approximately by:

Ae« -{R4-R1) q= 1,2,3,... . (B.22)

For example the TM mode cutoff frequency in 7 mm coaxial line for eq (B.22) is approximately 34 GHz.

98
Appendix C

Gap Correction

C.l Frequency-Dependent Gap Correction

C.1.1 Waveguide

We consider a sample in a rectangular waveguide of dimensions a x b with a small gap (bd). The dielectric constant in the gap is e^^ and of sample e/?^. The measured or observed

value is ejioA transverse resonance condition yields [40]:

tan(fci,(6 -d))-\- X thn{k2cd) = , (C.l)

where

LO ku = —V^Rs - ^Ro , (C.2)

Clab

^2c = —\/^Rg — ^Ro , (C.3)

Clab ^ and

^. ^

€fl5 \/efl7^£^

(C 4)

^Rg y^Rs — ^Ro

For our case of waveguide d ^ 2.54 x 10~^ m and b = 0.01016 m. Equation (C.l) must

be solved by iteration, but for low frequencies and low dielectric constants we can obtain a approximate solution. This equation reduces to Westphal's equation, [34] in the appropriate

low frequency limit,

^Rg^Rojb - d)

tRs = -T 1 • (<--5)

otRg — atRo

99
Figure C.l: Sample in waveguide with air gap.

C.1.2 Coaxial Line

For coaxial line the matching of transverse impedance yields the resonance condition [40]:

ct(x2, Xi) = A'ct(x2,a:3) , (C6)

where

^1 — ^Ic^l )

2^3 — ^2c-'^2 )

(C.7)

(C.8)

(C.9)

(CIO)

and X, kic, and /i2c a-re given in the previous section. The functions ct are defined as

J,{x)No{y) - N,{x)Jo{x)

ct{x,y) Mx)No{y) - No{x)Jo{x) ' (C.ll)

where Jq? Ji^ -^o? ^nd Ni are the Bessel functions of zero and first order of the first and

second kind respectively.

100
C.2 Frequency-Independent Approaches

Various researchers have approached the gap problem by representing the sample with air

gap as a layered capacitor [33,34,41]. This approach assumes that the gaps between transmission hne and sample are effectively modeled by a set of capacitors in series. Champlin

[41] approached the problem using as a starting point the perturbation formula developed

by Schwinger [42]. By substituting into the perturbation formula approximations to the

field distribution in the various regions, they obtain an estimate for the effective permittivity. Their answer turns out to be fully equivalent to the capacitor model of Westphal

[34]. Champlin showed that Bussey's theory [33] is the first two terms in an expansion of

Westphal [34] and Champlin's models.

The capacitor model is frequency independent and thus is strictly valid only at lower

frequencies and dc. We expect the capacitor model to break down at higher frequencies

because the wavelength decreases with increasing frequency to a point where multiple

scattering dominates. In order to account for multiple scattering, it is necessary to develop

a theory that is frequency dependent.

C.2.1 Coaxial Capacitor Model for Dielectric IVIaterials

Consider a capacitor consisting of layers of dielectric and layers of air in a coaxial line

[43,44]. The dimensions are shown in figure C.2.

We treat the system as capacitors in series, so

7^ = 7^ + 7^ + 7^- (^-12)

We know that for a coaxial hne the electric field distribution is given by

Er = —^ , (C.13)

and the voltage between the conductors is given by

V = - [ E{r)dr . (C.14)

The capacitance of a coaxial line of length L is given by

^ = ^; (C.15)

thus, for a system of three capacitors in series we have

\n^ Ing^ \n^ In — f- ^ = -^ + -^ + ^, (C.16)

101
Figure C.2: A coaxial sample in holder with air gaps near conductors with diameters

denoted by Di.

where e[, e^ are the corrected and measured values of the real part of the permittivity and

e[ is the real part of the permittivity of the air gap. We can then write the corrected values

as [34]

e. = e.

, L3 (l + AL)[L3-^i4(l + AL)]

where A^g = (tan -f^). The loss tangent is m

. ££ _ LsAme

An approximate expression is given by

where

tan^c = tan(5„i[l + ^'c~{~] '

^1 = In — + In — ,

(C.17)

(C.18)

(C.19)

(C.20)

(C.21)

102
e/ZEn,'

4045

Figure C.3: The gap correction calculated for various values of e'^, where R^^Ri are the

radii of the inner conductor and sample respectively.

Z/2 = In

L3 = In

R3

R2

i?4 Ri

(C.22)

(C.23)

Equation (C.19) breaks down when e'^ > L3/L1. An example is plotted in figure C.3 for a

7 mm coaxial line.

C.2.2 Rectangular Waveguide Model

For the case of a rectangular guide eqs (C.19) and (C.20) are modified slightly

^ '- [6-(6-c/)C(l+AL)]^ + &^A

The loss tangent is

""6-(6-<i)C[l + AL]

The approximate expressions for (A„ie < •!) are:

2me

(C.24)

(C.25)

103
<

ZH .'^t- {b--d)e 1

an <5c = tan Sr,'bb d^m

'

(C.26)

(C.27)

where b and d refer to the waveguide and sample dimensions and c and m refer to corrected

value and measured value.

C.3 Gap Correction for Magnetic Materials

C.3.1 Coaxial Line

For the calculation of the gap correction for the permeability a pure inductance model is useful. We model the transmission line as a series of inductors

where c, m, and air denote corrected value, measured value, and air space. Therefore the

corrected value is

The inductance is the flux penetrating the circuit divided b}' the current flowing in the

circuit

L = y , (C.30)

where

Ampere's law is

which yields

Therefore

so

/ ' ^-1 B-dS .

I"- dl = I ,

B,=

2-Kr

^=2.''/I\nb/a , -^u'ln6/a .

(C.31)

(C.32)

(C.33)

(C.34)

(C.35)

104
5.5

5.4 -

13CO 5.3

>

-o

0) tto

(1) 5.2

^_

^_oO

5.1

5

5 10 15

Gap (|im)

20 25

Figure C.4: Corrected permeability and permittivity as a function of inner conductor gap

for a 7 mm sample. The gap around the outer conductor is assume to be zero. In this case the uncorrected (measured) was //^ = 5, e^ = 5.

Therefore we can write for the corrected permeability

f^cR

li'^ji In R4/R, - [In R2/R1 + In R^/Rs]

lni?3/i?2

IniRJR,)

l^cR

II f^mR

{CM)

(C.37)

HRs/R^)

Gap corrections are given in figures C.4- C.7. The corrections for permeability in coaxial

line are much less than for permittivity. This is due to the fact the azimuthal magnetic

field is continuous across the discontinuity, whereas the radial electric field is discontinuous

across the discontinuity.

C.3.2 Waveguide

For waveguide the corrections are

/ __ / f^cR — f^mR

b-d'

(C.38)

105
CD _2

>

D

-'O

CD

i_

OO

5 10 15

Gap (jam)

Figure C.5: Corrected permeability and permittivity as a function of inner conductor gap

for a 7 mm sample. The gap around the outer conductor is assume to be zero. In this case the uncorrected (measured) was //^ = 25, e'j^ = 25.

1^ Re

10 15

Gap (^im)

20 25

Figure C.6: Corrected permeability as a function of inner conductor gap for a 7 mm sample.

The gap around the outer conductor is assume to be zero. In this case the uncorrected

(measured) was //^ = 50.

106
225

200

10 15

Gap (|j.m)

20 25

Figure C.7: Corrected permeability as a function of inner conductor gap for a 7 mm sample.

The gap around the outer conductor is assume to be zero. In this case the uncorrected

(measured) was //^ = 200.

l-'-cR — l-'-mR (C.39)

C.4 Gap Correction Formulas Derived Directly From

Maxwell's Equations

Consider Maxwell's equations in a coaxial line

V X E = -jujB . (C.40)

The radial electric field is discontinuous at the air gap interface, but the displacement

vector D is continuous across the interface. Also //^ is continuous across the interface,

whereas B is discontinuous. Let us assume that there are no sources so that J = 0. Then,

we can write eq (C.40) as

6 -JUI^L H . (C.41)

107
If we now average eq (C.41) over the cross-sectional area of a coaxial line, we obtain

27r / V X —dr = -]uj2t: / iiHdr . (C.42)

Continuity of the displacement field and tangential magnetic field and the fact that Dr oc 1/r and //^ oc 1/r is imposed and integrations are performed. Comparing these results

to an effective medium equation, we obtain the same form as the previously developed

capacitor and inductance models for the corrected permittivity and permeability:

\n^ \n^ \nf- In — f^ ^ = —^ + —^

+ —^ , (C.43)

Cm cj e', t\

/^U[ln RaIRx] = l^'cR In R3/R2 + [In R2/ R^ + In RJR3] . (C.44)

The analogous calculation can be performed for waveguide. These effective media formulas

are the Voigt approximation for the permeability (layers in series) and the Reuss approximation (layers in parallel) for the permittivity. The previously developed capacitor model

can be derived directly from Maxwell's equations.

From the perspective of Maxwell's equation the limitations of the capacitor and inductance models can be assessed. In order for these models to apply, we assume that

• The fundamental mode is the only propagating mode.

• The air gap and sample are azimuthally symmetric.

The air gap modifies the modal structure in the waveguide. The model assumption that

only a TEM mode is propagating in a coaxial line with an air gap around the sample

becomes less and less valid as the air gap increases in size. In fact, since the phase velocity

in the air gap region is much larger than the phase velocity in the sample, a distortion of

the wave can be expected. Longitudinal components of the electric field and therefore TM

modes will form as a result of an air gap if above cutoff. If the air gap or sample are not

azimuthally symmetric, H^ is no longer the only nonzero magnetic field component. This

asymmetry will allow higher order TE modes to propagate, when they are above cutoff.

When the assumption that the only propagating mode is the fundamental mode breaks

down, equations of the same form as eqs (C.43) and (C.44) hold with the logarithmic

constants replaced by more complicated expressions.

C.5 Mitigation of Air Gap Effects

It is possible to minimize the effects of air gaps by placing a conducting material in the air

gap. This material may be a conducting paint or a conducting grease. The conducting material will change the line impedance and line loss to a degree. However, for relatively small

108
gaps, the improvements in dielectric and magnetic property measurements far outweigh any

changes in hne impedance. The loss measurement will be influenced by his procedure.

109
Appendix D

Causal Functions and Linear

Response

D.l Introduction

We call a temporal function causal if it is for all times less than 0. The goal of this section

is to review the basic mathematics used to describe causal systems.

In the analysis to follow we assume that an impulse is applied to a system at ^ == 0. We can model a hnear system by an input function f{t), an output function g{t)^ and an

impulse response function a{t). It is possible for either or both of f{t) and a{t) to be causal.

A more general approach would be to study nonlinear response with linear response as a very special case. Linear response theory is usually valid when the underlying probability

density function can be approximated as an equilibrium distribution.

If both f{t) and a{t) are causal, the linear response is given by [45]

g{t)= f

f{r)a{t-r)dr . (D.l)

JO In this case, the output is a convolution over all past times. If only f[t) is causal then

/CO g{t) = / f{r)a{t - r)dT . (D.2)

Jo We assume that a{t) = for ^ < 0. The interpretation of the function f{t) can be obtained

from the relationship /oo

f{T)S{t - T)dT . (D.3)

-co

no
If we identify the left side of eq (D.3) as the response function, then the impulse response

a{t) is delta function in the special case of no distortion. Of course, in real systems the

impulse response will be broader than a delta function.

We can define the step response function h{t) as

dhit-r)

\^

^ ^a{t-r), (D.4)

and it is assumed that /i(0) = 0. If we use the step response for g{i) in eq (D.l) we obtain

^(0 = /(0)MO+ f^h{t-r). (D.5)

Jo ar

D.2 Transfer Functions

If we send a signal exp(ju;i) into a system, the response of the system is called the transfer

function and is denoted by S. The transfer function is defined as the Fourier transform of

the impulse response function

/oo

a{t)e''^'dt . (D.6)

-oo

As a consequence of the reality of a(^), S{—lj) = S*{ij). If F{lj) is the Fourier transform

of f{t) we have, using the inverse Fourier transform

g{t) = —P / F[u)S[u)e^-'du . (D.7)

ZTT 7-00

D.3 Kramers-Kronig Relations

The real and imaginary components of any causal function are related by a dispersion relation. The complex permittivity is a causal function whose real and imaginary components

are related by the Hilbert transform [45]

e.(t^ - e.oo = / 7^ ^ dO , (D.8)

The following summarizes some of the features of the Kramers-Kronig relations:

• The Hilbert transform relates real and imaginary components of a causal function.

Ill
• Direct solution requires complete data over full spectrum for one component.

• Equation (D.8) can be thought of as an integral equation for the unknown component

when there are some data for the other component.

Another form of the dispersion relations is

,. ^ ., . ^-<^0 f^ €-''i6)d6

TT J-oo [0 — LO)[6 — Uo)

er(^)-eK) = P -. TT^ r, (D.ll

where P denotes principal value.

*U.S. GOVERNMENT PRINTING OFFICE: 1992-673-025/69006

112

BL-114A

(5-90) ADMAN 15.01

U.S. DEPARTMENT OF COMMERCE

NATIONAL INSTITUTE OF STANDARDS AND TECHNOLOGY

BIBLIOGRAPHIC DATA SHEET

1. PUBLICATION OR REPORT NUMBER

NIST Tech Note 1355

2. PERFORMING ORCiANIZATION REPORT NUMBER

B92-0098

3. PUBLICATION DATE May 1992

4. TITLE AND SUBTITLE . ... Transmission/Reflection and Short Circuit Line Methods for Permittivity and Permeability

Determination

5. AUTHOR(S)

Lhmes R. Baker-Jarvis, Michael D. Janezic, John H. Grosvenor, Richard G. Geyer

6. PERFORMING ORGANIZATION (IF JOINT OR OTHER THAN NIST, SEE INSTRUCTIONS)

U.S. DEPARTMENT OF COMMERCE

NATIONAL INSTITUTE OF STANDARDS AND TECHNOLOGY

BOULDER, COLORADO 80303-3328

7. CONTRACT/GRANT NUMBER

8. TYPE OF REPORT AND PERIOD COVERED

9. SPONSORING ORGANIZATION NAME AND COMPLETE ADDRESS (STREET, CITY, STATE, ZIP)

10. SUPPLEMENTARY NOTES

11. ABSTRACT (A 200-WORD OR LESS FACTUAL SUMMARY OF MOST SIGNIFICANT INFORMATION. IF DOCUMENT INCLUDES A SIGNIFICANT BIBUOGRAPHY OR

LITERATURE SURVEY, MENTION IT HERE.)

A transmission/reflection and short-circuit line methods for measuring

complex permittivity and permeability of materials in waveguides and coaxial

lines are examined. Equations for complex permittivity and permeability are

developed from first principles. In addition, new formulations for the

determination of complex permittivity and permeability independent of

reference plane position are derived. For the one-sample

transmission/reflection method and two-position short-circuit line

measurements, the solutions are unstable at frequencies corresponding to

integral multiples of one-half wavelength in the sample. For two-sample

methods the solutions are unstable for frequencies where both samples

resonate simultaneously. Criteria are given for sample lengths to maintain

stability. An optimized solution is also presented for the scattering

parameters. This solution is stable over all frequencies and is capable of

reducing scattering parameter data on materials with higher dielectric

constant. An uncertainty analysis for the various techniques is developed

and the results are compared. The errors incurred due to the uncertainty in

scattering parameters, length measurement, and reference plane position are

used as inputs to the uncertainty models.

12. KEY WORDS (6 TO 12 ENTRIES; ALPHABETICAL ORDER; CAPITAUZE ONLY PROPER NAMES; AND SEPARATE KEY WORDS BY SEMICOLONS)

calibration; coaxial line; dielectric constant; loss factor; magnetic

materials; microwave measurements; permeameter; permeability measurement;

permittivity measurement; reflection method; short-circuit; transmission;

uncertainty; waveguide

13. AVAILABIUTY UNUMITED

FOR OFFICIAL DISTRIBUTION. DO NOT RELEASE TO NATIONAL TECHNICAL INFORMATION SERVICE (NTIS). X ORDER FROM SUPERINTENDENT OF DOCUMENTS, U.S. GOVERNMENT PRINTING OFFICE, WASHINGTON, DC 20402. ORDER FROM NATIONAL TECHNICAL INFORMATION SERVICE (NTIS), SPRINGFIELD, VA 22161.

14. NUMBER OF PRINTED PAGES

120

15. PRICE

ELECTRONIC FORM
NISTTechnical Publications

Periodical

Journal of Research of the National Institute of Standards and Technology—Reports NIST research and

development in those disciplines of the physical and engineering sciences in which the Institute is active. These include physics, chemistry, engineering, mathematics, and computer sciences. Papers cover a broad

range of subjects, with major emphasis on measurement methodology and the basic technology underlying

standardization. Also included from time to time are survey articles on topics closely related to the

Institute's technical and scientific programs. Issued six times a year.

Nonperiodicals

Monographs—Major contributions to the technical literature on various subjects related to the Institute's scientific and technical activities. Handbooks—Recommended codes of engineering and industrial practice (including safety codes) developed

in cooperation with interested industries, professional organizations, and regulatory bodies.

Special Publications—Include proceedings of conferences sponsored by NIST, NIST annual reports, and

other special publications appropriate to this grouping such as wall charts, pocket cards, and bibliographies.

Applied Mathematics Series—Mathematical tables, manuals, and studies of special interest to physicists,

engineers, chemists, biologists, mathematicians, computer programmers, and others engaged in scientific and

technical work.

National Standard Reference Data Series—Provides quantitative data on the physical and chemical

properties of materials, compiled from the world's literature and critically evaluated. Developed under a worldwide program coordinated by NIST under the authority of the National Standard Data Act (Public Law

90-396). NOTE: The Journal of Physical and Chemical Reference Data (JPCRD) is published quarterly for NIST by the American Chemical Society (ACS) and the American Institute of Physics (AIP). Subscriptions,

reprints, and supplements are available from ACS, 1155 Sixteenth St., NW, Washington, DC 20056.

Building Science Series—Disseminates technical information developed at the Institute on building materials,

components, systems, and whole structures. The series presents research results, test methods, and

performance criteria related to the structural and environmental functions and the durability and safety characteristics of building elements and systems.

Technical Notes—Studies or reports which are complete in themselves but restrictive in their treatment of

a subject. Analogous to monographs but not so comprehensive in scope or definitive in treatment of the subject area. Often serve as a vehicle for final reports of work performed at NIST under the sponsorship

of other government agencies.

Voluntary Product Standards—Developed under procedures published by the Department of Commerce in Part 10, Title 15, of the Code of Federal Regulations. The standards establish nationally recognized

requirements for products, and provide all concerned interests with a basis for common understanding of the characteristics of the products. NIST administers this program as a supplement to the activities of the private

sector standardizing organizations.

Consumer Information Series—Practical information, based on NIST research and experience, covering

areas of interest to the consumer. Easily understandable language and illustrations provide useful background

knowledge for shopping in today's technological marketplace.

Order the above NIST publications from: Superintendent of Documents, Government Printing Office,

Washington, DC 20402.

Order the following NISTpublications—FIPS and NISTlRs—from the National Technical Information Service,

Springfield, VA 22161.

Federal Information Processing Standards Publications (FTPS PUB)—Publications in this series collectively constitute the Federal Information Processing Standards Register. The Register serves as the

official source of information in the Federal Government regarding standards issued by NIST pursuant to the Federal Property and Administrative Services Act of 1949 as amended. Public Law 89-306 (79 Stat. 1127),

and as implemented by Executive Order 11717 (38 FR 12315, dated May 11, 1973) and Part 6 of Title 15 CFR (Code of Federal Regulations).

NIST Interagency Reports (NISTIR)—A special series of interim or final reports on work performed by

NIST for outside sponsors (both government and non-government). In general, initial distribution is handled

by the sponsor; public distribution is by the National Technical Information Service, Springfield, VA 22161,

in paper copy or microfiche form.
U.S. Department of Commerce

National Institute of Standards and Technology

325 Broadway

Boulder, Colorado 80303-3328

OFFICIAL BUSINESS

PENALTY FOR PRIVATE USE. $300