Acessibilidade / Reportar erro

Spectroscopy studies of 4H-SiC

Abstract

Calculations of the total dielectric functions and the optical bandgap energy (OBGE) of 4HSiC were performed by the full-potential linear muffin-tin-orbital method. The results are compared to spectroscopic ellipsometry dielectric measurements agreeing closely over in a wide range of energies. The obtained theoretical value of the (OBGE) agrees very closely with the measured ones obtained by transmission and photoacoustic spectroscopies at room temperature performed on 470 µm thick wafer and a 25 µm thick homoepitaxial layer of 4H-SiC samples grown (n-type, Siface) by hot wall CVD.

optical bandgap energy; wide-bandgap material; SiC polytypes; computer simulation


Spectroscopy studies of 4H-SiC

A.C. de OliveiraI; J.A. Freitas Jr.II; W.J. MooreII; A. Ferreira da SilvaIII; I. PepeIII,IV; J. Souza de AlmeidaIII,V; J.M. Osório-GuillénVI; R. AhujaVI; C. PerssonVI; K. JärrendahlVII; O.P.A. LindquistVII; N.V. EdwardsVII,VIII; Q. WahabVII

IInstituto de Física, Universidade de Brasília, 70919-970, Brasília, DF, Brazil

IINaval Research Laboratory, Washington, D.C. 20375, USA

IIIInstituto de Física, Universidade Federal da Bahia, 40210-340, Salvador, Ba, Brazil

IVLPCC, College de France, F-75231 Paris, France

VDCC, Faculdade Ruy Barbosa, 41940-320 Salvador, Ba, Brazil

VIDepartment of Physics, Uppsala University, SE-75121 Uppsala, Sweden

VIIDepartment of Physics, Linköping University, SE-58183 Linköping, Sweden

VIIIMotorola Semiconductor Products, 2200 West Broadway Road, Mesa AZ 85202 USA

Address to correspondence Address to correspondence A.C. de Oliveira e-mail: aderbal@unb.br

ABSTRACT

Calculations of the total dielectric functions and the optical bandgap energy (OBGE) of 4HSiC were performed by the full-potential linear muffin-tin-orbital method. The results are compared to spectroscopic ellipsometry dielectric measurements agreeing closely over in a wide range of energies. The obtained theoretical value of the (OBGE) agrees very closely with the measured ones obtained by transmission and photoacoustic spectroscopies at room temperature performed on 470 µm thick wafer and a 25 µm thick homoepitaxial layer of 4H-SiC samples grown (n-type, Siface) by hot wall CVD.

Keyword: optical bandgap energy, wide-bandgap material, SiC polytypes, computer simulation

1. Introduction

The physical properties of the large bandgap 4H-SiC make it a prominent material for high-power, high-temperature, and high-frequency devices. Devices like field effect transistors, bipolar storage capacitors, and ultraviolet detectors have been fabricated1,2,3.

In this work, we have investigated the optical properties of undoped and n-type 4H-SiC, both experimentally and theoretically. The transmission and photoacoustic spectroscopy techniques have been used for the measurement of the optical bandgap energy3-5. The total dielectric functions were determined by spectroscopic ellipsometry (SE), a powerful non-destructive technique for high accuracy measurements5-7. The calculations of the total dielectric functions were performed by the full-potential linear muffin-tin-orbital method5,8.

2. Experimental Details

The experimental transmission and photoacoustic spectroscopy apparatus consists of a halogen lamp used as the light source for the measurement. The polychromatic beam is diffracted by a plane diffraction gratings attached to a step-motor. The beam can be varied from 925 to 360 nm, a set of lens and collimator produces a monochromatic light focused onto the sample, see Refs. 3-5.

3. Calculation of the Dieletric Function

The dielectric function was calculated in the momentum representation, which requires matrix elements of the momentum, p, between occupied and unoccupied eigenstates. To be specific the imaginary part of the dielectric function, e2(w) = Ime(q = 0,e), was calculated from5

In Eq. 1, e is the electron charge, m its mass, W is the crystal volume and ƒkn is the Fermi distribution. Moreover, is the crystal wave function corresponding to the nth eigenvalue with crystal momentum k and spin s.

With our spherical wave basis functions, the matrix elements of the momentum operator are conveniently calculated in spherical coordinates and for this reason the momentum is written , where m is -1, 0, or 1, and

The evaluation of the matrix elements in Eq. 2 is done over the muffin-tin region and the interstitial separately5. A full detailed description of the calculation of the matrix elements was presented somewhere else9. In our theoretical method the wave function [Eq. 1] inside the muffin-tin spheres is atomic-like in the sense that it is expressed as a radial component times spherical harmonic functions (also involving the so called structure constants), i.e.,

where inside the sphere R' the basis function is

and

In the equations above is the muffin-tin radius for atom R, is the structure constant, D is the logarithmic derivative, is the numerical solution to the spherical component of the muffin-tin potential and is the ener-gy derivative of . Therefore, this part of the problem is quite analogous to the atomic sphere approximation (ASA) calculations and we calculate the matrix elements in Eq. 2 as5

Since, ct involves a radial function multiplied with a spherical harmonic function, i.e., (we will label this product we can calculate the matrix elements in Eq. 5 using the relations

and

A general matrix element, is then calculated using the Wigner-Eckart theorem5. Matrix elements of the momentum over the interstitial are obtained from the relation

By use of Green's second theorem the expression above can be expressed as

In the expression above the surface integral is taken over the muffin-tin spheres and since the interstitial wave function by construction is the same as the muffin-tin wave function we can use for y the numerical wave function defined inside (and on the boundary of) the muffin-tin spheres. In this way the evaluation of the integral above is done in a quite similar way as done for the muffin-tin contribution [Eqs. 6-8] to the gradient matrix elements.

The summation over the Brillouin zone in Eq. 1 is calculated using linear interpolation on a mesh of uniformly distributed points, i.e., the tetrahedron method. The matrix elements, eigenvalues and eigenvectors are calculated in the irreducible part of the Brillouin-zone.

The correct symmetry for the dielectric constant was obtained by averaging the calculated dielectric function. Finally, the real part of the dielectric function, e1(w), is obtained from e2(w) using the Kramers-Kronig transformation,

e1(w) º Re( e(q=0,w) ) =

The total average dielectric functions are given by

4. Results

The values of the optical bandgap energies are shown in Table 1. They are very close when compared to different measurements as well as to the theoretical result. From the electronic structure and the partial density of states we can identify the optical absorption described by the imaginary part of the dielectric function and then obtain the bandgap energy5. The measured total dielectric functions, e1 and e2 are shown in Fig. 1 together with the theoretical results. There is a good agreement between the measured and the calculated dielectric functions.


5. Conclusion

We have measured the optical properties of 4H-SiC, using transmission and photoacoustic spectroscopies and spectroscopic ellipsometry. The resulting total dielectric functions and optical bandgap energies are in excellent agreement with our full-potential calculations. By comparing the total dielectric functions we conclude that the optical anisotropy is small in 4H-SiC.

Acknowledgements

This work was financially supported by: the Swedish agencies TFR, NFR and SSF, Brazilian agencies CNPq and FAP-DF, and at the Naval Research Laboratory by the Office of Naval Research.

3. Oliveira, A.C.; Freitas Jr., J.A.; W.J. Moore (to appear).

5. Ahuja, R.; Silva, A.F.; Persson, C.; Osorio-Guillén, J.M.; Pepe, I.; Järrendahl, K.; Lindquist, O.P.A.; Edwards, N.V.; Wahab, Q.; Johansson, B. J. Appl. Phys. (in press).

8. Wills, J.M.; (unpublished); Wills, J.M.; Cooper, B.R. Phys. Rev. B, v. 36, p. 3809, 1987; Price, D.L.; Cooper, B.R. ibid., v. 39, p. 4945, 1989.

Received: January 02, 2002

Revised: September 30, 2002

  • 1. Morkoc, H.; Strite, S.; Gao, G.B.; Lin, M.E.; Sverdlov, B.; Burns, M. J. Appl. Phys., v. 76, p. 1363, 1994.
  • 2. Palmour, J.W.; Edmond, J.A.; Kong, H.S.; Carter Jr., C.H. Physica B, v. 185, p. 461, 1993.
  • 4. Silva, A.F.; Veissid, N.; An, C.Y.; Pepe, I.; Oliveira, N.B.; Silva, A.V.B. Appl. Phys. Lett., v. 69, p. 1930, 1996.
  • 6. S. Zangooie, P.O.A. Persson, J.N. Hilfiker, L. Hultman, H. Arwin J. Appl. Phys., v. 87, p. 8497, 2000.
  • 7. Lindquist, O.P.A; K. J\"{a}rrendahl, Peters, S.; Zettler, J.-T.; Cobet, C.; Esser, N.; Aspnes, D.E. Appl. Phys. Lett., v. 78, p. 2751, 2001.
  • 9. Ahuja, R.; Auluck, S.; Wills, J.M.; Alouani, M.; Johansson, B.; Eriksson, O. Phys. Rev. B, v. 55, p. 4999, 1997.
  • 10. A good description of the calculation of dielectric constants and related properties is found in the thesis by T. Gashe, Uppsala University, 1993.
  • Address to correspondence

    A.C. de Oliveira
    e-mail:
  • Publication Dates

    • Publication in this collection
      25 Mar 2003
    • Date of issue
      Jan 2003

    History

    • Reviewed
      30 Sept 2002
    • Received
      02 Jan 2002
    ABM, ABC, ABPol UFSCar - Dep. de Engenharia de Materiais, Rod. Washington Luiz, km 235, 13565-905 - São Carlos - SP- Brasil. Tel (55 16) 3351-9487 - São Carlos - SP - Brazil
    E-mail: pessan@ufscar.br