Reference notes for the use of complex refractive indices with MODEIG

 

J. H. Abeles

G. A. Evans

 

 

1. How to input a complex refractive index to MODEIG:

 

MODEIG, and particularly MODEIG/II used here, is capable of handling complex refractive indices. Thus, waveguiding properties of structures which include lossy layers can be found.

 

The parameter NLOSS is used to provide MODEIG with the complex part of the refractive index. Specifically, MODEIG requires that NLOSS be expressed as:

 

 

where a(amplitude) represents the amplitude absorption length in a given material in units of µm-1. By a widely-accepted convention the symbol a = 2a(amplitude) = 2 NLOSS is called the extinction coefficient and describes the exponential decrease in power of a wave traveling in the material.

 

Following the notation of Pankove, the complex refractive index nc is expressed in terms of its real and imaginary components:

 

 

The z-component of the propagation of an electromagnetic wave in such a material is described by the exponential

 

 

where k0=2p/l. Note the second exponential term on the far right. It describes an exponentially-decaying field amplitude with decay constant equal to kk0. Thus:

 

 

In its output file, MODEIG displays the complex refractive index nc (it calls it "RN") which it calculated from NREAL and NLOSS. Thus, in a case where the user begins with nc and converts its complex part (i.e., k) to NLOSS units, the user can check that the conversion was performed properly.

 

 

2. Example:

 

As an example, the complex refractive index of gold is given as nc = 0.188 + 5.39 i at 826.6 nm. From this it may be calculated that NLOSS = 40.97 mm-1 for this wavelength. Using MODEIG/II, the propagation constant b for an (Al,Ga)As waveguide structure coated with gold is calculated. The MODEIG/II layer file shows the structure:

 

# of layers = 3

LAYER01 NLOSS=40.97000 NREAL= .18800 TL= .00000

LAYER02 NLOSS= .00000 NREAL= 3.44437 TL= .50000

LAYER03 NLOSS= .00000 NREAL= 3.25303 TL= .00000

 

These numbers correspond to the structure given in a MODEIG/II input file containing the following LAYER data (corresponding to a 0.5 mm thick Al0.3Ga0.7As guiding layer embedded between the gold and an Al0.6Ga0.4As cladding layer):

 

LAYER NREAL=0.188 NLOSS=40.97

LAYER ALPERC=.30 TL=0.5

LAYER ALPERC=.6

 

The MODEIG/II database file shows the calculated mode parameters WZR and WZI,

PHM WZR WZI QZR

7.762159E-01 3.384081E+00 9.409423E-05 1.145201E+01

 

which are respectively the real and imaginary parts of the modal propagation constant b (analogous with the complex index of refraction nc which by comparison describes propagation in a uniform material rather than a waveguide structure):

 

 

Thus, the modal extinction coefficient is

 

 

As mentioned above, the output file shows the correct complex refractive index nc (called RN by MODEIG/II) and an excerpt of the file for the example shows

 

 

PE=(-29.01578,2.02661) PM=(1.00000,.00000)

L 1 SEMI-INFINITE QN=(-29.01578,2.02661) YN=(1.00000,.00000)

RN=(.18800,5.38991)

 

finding RN = nc as expected.

 

Another example of the use of a complex index of refraction is shown in Fig. 1, where the mode losses due to a metal layer, here taken to be either gold or titanium, deposited directly on the p-clad is calculated as a function of p-clad thickness. The structure is a GRINSCH, single quantum well (100 Å) and the width of each graded region is 0.2 µm. A lasing wavelength of 8550 Å is assumed. The aluminum mole fraction in both claddings is 60%, and that of the graded regions range from 60 to 20%. The complex index taken for gold is the same as used above, and that for titanium is 4.4 + i4.84. Gold has the desirable property that, at wavelengths from 0.6 µm to 1.5 µm, the real part of the index of refraction is very small, especially compared to the imaginary part. This makes it possible to deposit gold directly over waveguides, including those incorporating gratings which may serve as a reflector in surface emitting lasers.

 

3. Additional comments:

 

NLOSS is also commonly used when the background losses in epilayers are included in MODEIG calculations. For example, a cladding layer may have a loss of 10 cm-1, which corresponds to an NLOSS of 0.0005 mm-1 (or NLOSS=0.0005).

 

Since NLOSS is amplitude loss in mm-1, rather than power loss in cm-1, the use of NLOSS can be confusing. This was not a problem eight or nine years ago when this was coded because there was only one user (G. Evans). As more people use this version of MODEIG, we could change NLOSS to be a power loss in cm-1. If you would like this change, let us know.

 

On the same subject of defining the index of layers, note that the MODEIG variable PEROPT allows entering either the complex index of refraction nc or the complex dielectric constant e.

 

To convert between a complex refractive index nc and a complex dielectric constant, e = e1 + i e2 = nc2, the following expressions are appropriate:

 

 

4. Table of complex refractive indices of metals:

 

The following table is taken from "Optical properties of the metals Al, Co, Cu, Au, Fe, Pb, Ni, Pd, Pt, Ag, Ti, and W in the infrared and far infrared," M. A. Ordal, L. L. Long, R. J. Bell, S. E. Bell, R. R. Bell, R. W. Alexander, Jr., and C. A. Ward, Appl. Optics 22, 1099 (1983) and "Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W," M. A. Ordal, R. J. Bell, R. W. Alexander, Jr., L. L. Long, and M. R. Querry,

Appl. Optics 24, 4493 (1985).

 

Values for wavelengths closest to 850 nm were selected. Multiple values for a given metal result from multiple referenced works finding those values as referenced in the above two papers.

 

Metal

wavelength (mm)

n

k

---------------------------

---------------------------

---------------------------

---------------------------

Aluminum

0.827

2.75

8.31

 

0.850

2.08

7.15

Cobalt

0.821

2.53

4.88

 

0.827

3.10

4.96

Copper

0.850

2.08

7.15

 

0.850

1.20

5.47

 

0.827

2.60

5.26

Gold

0.900

0.18

5.72

 

0.827

0.08

4.98

Iron

0.886

3.12

3.87

Lead

0.850

1.44

4.35

Molybdenum

0.954

2.77

3.74

Nickel

0.855

2.59

4.55

 

0.892

2.40

5.23

Palladium

0.827

2.17

5.22

 

0.821

2.06

5.19

Platinum

0.827

2.92

5.07

 

0.830

2.52

4.67

Silver

0.850

0.10

5.85

 

0.827

0.27

5.79

Titanium

0.821

3.21

4.01

Vanadium

0.984

2.94

3.50

Tungsten

0.827

3.48

2.79

 

August 16, 1989