Hostname: page-component-7c8c6479df-995ml Total loading time: 0 Render date: 2024-03-28T23:11:34.941Z Has data issue: false hasContentIssue false

Doping Dependence Of The Thermal Conductivity Of Hydride Vapor Phase Epitaxy Grown n-GaN/Sapphire (0001) Using A Scanning Thermal Microscope

Published online by Cambridge University Press:  13 June 2014

D.I. Florescu
Affiliation:
Physics Department and New York State Center for Advanced Technology in Ultrafast Photonic Materials and Applications, Brooklyn College of CUNY, Brooklyn, NY 11210
V.A. Asnin
Affiliation:
Physics Department and New York State Center for Advanced Technology in Ultrafast Photonic Materials and Applications, Brooklyn College of CUNY, Brooklyn, NY 11210
L.G. Mourokh
Affiliation:
Physics Department and New York State Center for Advanced Technology in Ultrafast Photonic Materials and Applications, Brooklyn College of CUNY, Brooklyn, NY 11210
Fred H. Pollak
Affiliation:
Physics Department and New York State Center for Advanced Technology in Ultrafast Photonic Materials and Applications, Brooklyn College of CUNY, Brooklyn, NY 11210
R. J. Molnar
Affiliation:
Massachusetts Institute of Technology, Lincoln Laboratory, Lexington, MA 02420-9108

Abstract

We have measured the doping concentration dependence of the room temperature thermal conductivity (κ) of two series of n-GaN/sapphire (0001) fabricated by hydride vapor phase epitaxy (HVPE). In both sets κ decreased linearly with log n, the variation being about a factor two decrease in κ for every decade increase in n. κ ≈ 1.95 W/cm-K was obtained for one of the most lightly doped samples (n = 6.9×1016 cm−3), higher than the previously reported κ ≈ 1.7-1.8 W/cm-K on lateral epitaxial overgrown material [V.A. Asnin et al, Appl. Phys. Lett. 75, 1240 (1999)] and κ ≈ 1.3 W/cm-K on a thick HVPE sample [E.K. Sichel and J.I. Pankove, J. Phys. Chem. Solids 38, 330 (1977)]. The decrease in the lattice component of κ due to increased phonon scattering from both the impurities and free electrons outweighs the increase in the electronic contribution to κ.

Type
Research Article
Copyright
Copyright © 1996 Materials Research Society

Introduction

Despite the considerable body of work, both experimental and theoretical, on the electronic, optical, and structural properties of group III nitrides [Reference Gil1] relatively little work has been reported on thermal conductivity κ . This quantity is of importance from both fundamental and applied perspectives. The lattice thermal conductivity is a function of the mean free path of the phonons and hence is determined by both intrinsic (phonon-phonon Umklapp scattering) and extrinsic (phonon- ″defect″, phonon-carrier scattering) factors [Reference Bhandari and Rowe2]. Sichel and Pankove [Reference Sichel and Pankove3] determined κ of "bulk" hydride vapor phase epitaxy (HVPE) GaN as a function of temperature (25K<T<360K) with κ ≈ 1.3 W/cm-K at 300K. More recently Asnin et al [Reference Asnin, Pollak, Ramer, Schurman and Ferguson4] have performed high spatial resolution measurements on several lateral epitaxial overgrown (LEO) GaN/sapphire (0001) samples using a scanning thermal microscope (SThM) and found κ ≈ 1.7-1.8 W/cm-K [Reference Asnin, Pollak, Ramer, Schurman and Ferguson4]. Slack has estimated an upper limit of 1.7 W/cm-K at 300K for GaN [Reference Slack and Phys5].

We report high spatial resolution determination of κ at 300K on two sets of HVPE n-GaN/sapphire (0001) samples as a function of n. The measurements were made using a ThermoMicroscope’s SThM Discoverer system [6], with a spatial resolution of ≈ 2-3 P 9DO×HV RI n were deduced from both 300K Hall effect and micro-Raman [longitudinal optical phonon-plasmon (LPP)] measurements. In both sets of samples κ decreased linearly with log n, the variation being about a factor two decrease in κ for every decade increase in n. κ ≈ 1.95 W/cm-K was obtained for one of the most lightly doped samples (n = 6.9×1016 cm−3), higher than previously reported κ [Reference Sichel and Pankove3,Reference Asnin, Pollak, Ramer, Schurman and Ferguson4].

Sample set A had unintentional n (6-800×1016 cm−3) and thicknesses (t) in the range of 5-74 μm. For sample set B t was constant ≈ 10 μm and 15×1016 cm−3 < n < 300×1016 cm−3.

Our observation also helps to explain the results on the LEO material [Reference Asnin, Pollak, Ramer, Schurman and Ferguson4], which had n ≈ (10-20)×1016 cm−3 [Reference Pophristic, Long, Schurman, Ramer and Ferguson7]. The decrease in the lattice component of κ due to increased phonon scattering from both the impurities and free electrons outweighs the increase in the electronic contribution to κ.

Experimental Details

The GaN films were grown by the HVPE method in a vertical type reactor [Reference Molnar, Gotz, Romano and Johnson8]. During this process, gallium monochloride is synthesized upstream by reacting HCl gas with liquid Ga metal at 800-900°C.The GaCl is transported to the substrate where it is reacted with NH3 at 1000-1100 °C forming GaN. All films were grown on (0001) sapphire. The carrier concentration n H was determined by 300K Hall effect measurements. Several characteristics of the samples are listed in Table I.

Table I Summary of t, nH/nR, μ and κ of the two sets of samples.

The carrier concentration for n≥ 40×1016 cm−3 was also determined from the LPP modes observed in Raman scattering [Reference Perlin, Camassel, Knap, Taliercio, Chevin, Suski, Gregory and Porowski9] and compared to the Hall effect results. Raman microprobe (≈ 2μm) measurements were made in the backscattering geometry using a triple grating spectrometer (Jobin-Yvon model T64000) and the 488 nm line of an Ar-ion laser as excitation. The Raman system was equipped with an Olympus BH2 microscope. Values of the carrier concentration (nR ) were deduced from Eqs. (2) and (3) in Ref. [Reference Molnar, Gotz, Romano and Johnson8] (see Table I), using an electron effective mass (me* ) of 0.22 (in units of the free electron mass) [Reference Suzuki, Uenoyama and Gil10] and a high frequency dielectric constant (ε ) of 5.5 [Reference Popovici, Morkoç, Mohammad and Gil11].

The probe tip of the SThM system consists of a ″V″ shaped resistive thermal element incorporated at the end of a cantilever that enables atomic force microscopy-type feedback, as shown schematically in Fig. 1a. The arms of the cantilever are made of Wollaston process wire consisting of silver wire ≈ 75μm in diameter containing a platinum/10% rhodium core ≈ 3μm in diameter. The resistive element at its end comprises a 200μm length of platinum that has been exposed by removal of the silver and bent into a ″V″ shape (radius of curvature ≈ 1μm), which acts as the probe. The resistive element forms one leg of a Wheatstone bridge, as shown in Fig. 1b. A current is passed through the probe so that in air its temperature is about 40-50°C above ambient. There is a feedback loop to adjust the bridge voltage as necessary to keep the bridge balanced thus keeping the temperature of the probe constant. When the probe contacts the sample heat flows from the probe to the material, as shown in Fig. 1c. In the absence of feedback, this flow of heat would reduce the probe temperature, decreasing the resistance and causing the bridge to shift. The feedback senses this shift and increases the voltage applied to the bridge (Uout ), returning the resistance to its set point. The thermal conductivity κ is proportional to the heat flow or (Uout )2, as shown in Fig. 1c.

Fig. 1 (a) Schematic diagram of the SThM probe, (b) circuit diagram of the probe control, and (c) heat flow from the tip into the sample.

Although initially designed to measure only relative spatial variations in Uout , Ruiz et al [Reference Ruiz, Sun, Pollak and Venkatraman12] developed a calibration procedure that makes is possible to evaluate absolute values of κ. Based on the results of Ref. [Reference Hammiche, Pollack, Song and Hourston13] we estimate that the lateral/depth resolution is about 2 – 3 μm for materials with κ ≈ 1.5-2 W/cm-K.

Experimental Results

Shown in Table I are the measured values of κ at 300K for samples A1-A6 and B1-B7. Note that samples A1 and A2 with n = 6.3×1016 cm−3 and n = 6.9×1016 cm−3 have κ = 1.82±0.05 W/cm-K and 1.95±0.05 W/cm-K, respectively; the latter is the highest value of this parameter observed to date. Plotted in Figs. 2a and 2b are κ as a function of log n for the two sets of samples, respectively.

Fig. 2 Thermal conductivity as a function of carrier concentration, n, for samples (a) A1-A6 and (b) B1-B7.

Representative error bars are shown. The solid lines in the figure are least-square fits to a linear function. For both sets of samples κ(n) is essentially the same.

Discussion of Results

From kinetic theory the lattice κ is given by [Reference Bhandari and Rowe2]:

(1)

where vs is the average velocity of sound (with only a weak temperature dependence), c(T) is the lattice specific heat, l (T) is the phonon mean free path, and τ (T) is the lifetime.

In almost all materials κ (T) first increases with temperature, reaches a maximum (κmax ) at some characteristic temperature Tch, and then decreases [Reference Asnin, Pollak, Ramer, Schurman and Ferguson4]. At low temperatures l is relatively long and is dominated by extrinsic effects such as "defects" and/or finite crystal size and c(T)∝(T/ΘD)3, where ΘD is the Debye temperature. As the temperature increases c(T) begins to saturate and intrinsic temperature dependent Umklapp processes become dominant, thus causing a decrease in l. For GaN Tch ≈ 200K [Reference Sichel and Pankove3] and ΘD ≈ 600K [Reference Akasaki, Amano and Edgar14].

For T<Tch κ is very sensitive to "defect" density but still has some dependence in a range of T above Tch. Since at 300K we are close to Tch, the thermal conductivity will still be a function of n.

The observed dependence of κ on log n, as shown in Fig. 2, is difficult to account for in detail since we are in a regime where both extrinsic and intrinsic scattering processes are important. The scattering of the phonons from the impurity atoms can be described on the basis of mass-difference scattering, the relaxation time being given by [Reference Bhandari and Rowe2]:

(2)

where Mimp is the mass of the impurity atoms and Mav is the average mass of the atoms in the material. Thus there will certainly be a dependence of κ on n. However, because of the Umklapp contribution, phonon scattering from the free carriers, and the contribution of the free carriers to the electronic component of κ, the particular function form is not immediately evident. Clearly more work in this area needs to be done.

Slack has estimated an upper bound of κ ≈ 1.7 W/cm-K for GaN at 300K from the relation [Reference Slack and Phys5]:

(3)

where B is a constant, δ is the average volume occupied by one atom in the crystal, and is the Grüneisen parameter. By using the factor Mavδ ΘD 3 as a scaling parameter he deduced the above value of κ for GaN at 300K. However, this analysis is limited since the above expression is applicable only for T >> ΘD (≈ 600K in GaN).

Our observation also helps to explain the observed values of κ ≈ 1.7-1.8 W/cm-K on the LEO material [6], which had n ≈ (10-20)×1016 cm−3 [Reference Pophristic, Long, Schurman, Ramer and Ferguson7] (see Fig. 2).

Certain GaN devices, such as high power field effect transistors, laser diodes, etc., would benefit greatly from GaN with higher thermal conductivity, as heat extraction from the device becomes more efficient with higher κ. Also, GaN has many potential applications in the area of high temperature electronics, where a large κ is very advantageous [Reference Torvik, Leksono, Pankove and Van Zeghbroeck15]. Our highest observed value of κ ≈ 1.95-1.85 W/cm-K is somewhat smaller than single crystal AlN (≈2.85 W/cm-K) [Reference Slack, Tanzilli, Pohl and Vandersande16] and is considerably higher than that of sintered AlN material [Reference Slack and Phys5]; the latter is often used as a heat sink material. Thus GaN based devices could be fabricated on HVPE GaN/sapphire material with the above thermal conductivities, thus avoiding costly processing steps.

Summary

The doping dependence of the room temperature κ of two series of HVPE n-GaN/sapphire (0001) has been measured using a SThM. κ decreased linearly with log n, the variation being about a factor two decrease in κ for every decade increase in n in both sets of samples. The general behavior of κ(n), i.e., decrease with increasing n, is similar to other semiconductors in a comparable temperature range. For one of the most lightly doped samples (n = 6.9×1016 cm−3) κ ≈ 1.95 W/cm-K, higher than the previously reported κ on several LEO samples and a thick HVPE material. The decrease in the lattice component of κ due to increased phonon scattering from both the impurities and free electrons outweighs the increase in the electronic contribution to κ. The implications of these findings for device applications and design are discussed.

Acknowledgements

The Brooklyn College work was supported by Office of Naval Research contract N00014-99-C-0663 (administered by Dr. Colin Wood) and the New York State Science and Technology Foundation through its Centers for Advanced Technology program. The Lincoln Laboratory work was sponsored by the US Air Force under Air Force contract #F19628-95-C-002. Opinions, interpretations, conclusions, and recommendations are those of the authors and are not necessarily endorsed by the US Air Force.

References

See, for example, Group III Nitride Semiconductor Compounds, ed. by Gil, B. (Clarendon, Oxford, 1998).Google Scholar
See, for example, Thermal Conduction in Semiconductors by Bhandari, C.M., and Rowe, D.M. (Wiley, New York, 1988)Google Scholar
Sichel, E.K. and Pankove, J.I., J. Phys. Chem. Solids 38, 330 (1977).CrossRefGoogle Scholar
Asnin, V.A., Pollak, F.H., Ramer, J., Schurman, M., and Ferguson, I., Appl. Phys. Lett. 75, 1240 (1999).CrossRefGoogle Scholar
Slack, G.A., Phys, J.. Chem Solids 34, 321 (1973).CrossRefGoogle Scholar
Thermo Microscopes, 1171 Borregas Avenue, Sunnyvale, CA 94089.Google Scholar
Pophristic, M., Long, F.H., Schurman, M., Ramer, J., and Ferguson, I.T., Appl. Phys. Lett. 74, 3519 (1999).CrossRefGoogle Scholar
Molnar, R.J., Gotz, W., Romano, L.T., and Johnson, N.M., J.Crystal Growth 178, 147 (1997).CrossRefGoogle Scholar
Perlin, P., Camassel, J., Knap, W., Taliercio, T., Chevin, J.C., Suski, T., Gregory, I., and Porowski, S., Appl. Phys. Lett. 67, 2524 (1995).CrossRefGoogle Scholar
Suzuki, M., and Uenoyama, T. in Group III Nitride Semiconductor Compounds, ed. by Gil, B. (Clarendon, Oxford, 1998), p. 307.Google Scholar
Popovici, G., Morkoç, H., and Mohammad, S.N. in Group III Nitride Semiconductor Compounds, ed. by Gil, B. (Clarendon, Oxford, 1998), p. 19.Google Scholar
Ruiz, F., Sun, W.D., Pollak, F.H., and Venkatraman, C., Appl. Phys. Lett. 73, 1802 (1998).CrossRefGoogle Scholar
Hammiche, A., Pollack, H.M., Song, M., and Hourston, D.J., Meas. Sci. Technol. 7, 142 (1996).CrossRefGoogle Scholar
Akasaki, I., and Amano, H. in Properties of Group III Nitrides, ed. by Edgar, J.H. (INSPEC, London, 1994) p. 30.Google Scholar
Torvik, J.T., Leksono, M., Pankove, J.I., and Van Zeghbroeck, B., MRS Internet J. Nitride Semicond. Res. 4, 3 (1999).CrossRefGoogle Scholar
Slack, G.A., Tanzilli, R.A., Pohl, R.O., and Vandersande, J.W., J. Phys. Chem. Solids 48, 641 (1987)CrossRefGoogle Scholar
Figure 0

Table I Summary of t, nH/nR, μ and κ of the two sets of samples.

Figure 1

Fig. 1 (a) Schematic diagram of the SThM probe, (b) circuit diagram of the probe control, and (c) heat flow from the tip into the sample.

Figure 2

Fig. 2 Thermal conductivity as a function of carrier concentration, n, for samples (a) A1-A6 and (b) B1-B7.