Skip to main content Accessibility help
×
Home

Information:

  • Access
  • Cited by 1

Actions:

      • Send article to Kindle

        To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

        Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

        Find out more about the Kindle Personal Document Service.

        Monte Carlo Simulation of Hall Effect in n-Type GaN
        Available formats
        ×

        Send article to Dropbox

        To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

        Monte Carlo Simulation of Hall Effect in n-Type GaN
        Available formats
        ×

        Send article to Google Drive

        To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

        Monte Carlo Simulation of Hall Effect in n-Type GaN
        Available formats
        ×
Export citation

Abstract

Results of Monte Carlo simulations of electron transport for wurtzite phase GaN in crossed, weak electric and magnetic fields are presented. It is found that the Hall factor, r H = μ Hall drift , decreases monotonically as the temperature increases from 77K to 400K.The low temperature value of the Hall factor increases significantly with increasing doping concentration. The Monte Carlo simulations take into account the electron-lattice interaction through polar optical phonon scattering, deformation potential acoustic phonon scattering (treated as an inelastic process), and piezoelectric acoustic phonon scattering. Impurity scattering due to ionized and neutral donors is also included, with the latter found to be important at low temperature due to the relatively large donor binding energy which implies considerable carrier freeze-out already at liquid nitrogen temperature. The temperature dependences of the electron concentration, drift mobility, and Hall factor are calculated for donor concentrations equal to 5 × 1016 cm−3, 1017 cm−3, and 5 × 1017 cm−3. The Monte Carlo simulations are compared to classical analytical results obtained using the relaxation-time approximation, which is found to be adequate at low temperatures and sufficiently low carrier concentrations so that inelastic scattering effects due to optical phonons and degeneracy effects are negligible. The influence of dislocations on the Hall factor is discussed briefly.

Footnotes

MRS Internet J. Nitride Semicond. Res. 4S1, G6.6 (1999)

Introduction

Electronic characterization of epitaxially grown GaN relies on temperature dependent Hall measurements to accurately determine carrier concentration and mobility. Analysis of measured Hall data has already led to the examination of several effects unique to GaN Hall measurements. Specifically, the variation of the electron Hall mobility in the presence of a highly conducting layer at the interface between the substrate and the epitaxial GaN crystal has been resolved 1 , and the dependence of the Hall measurement on the orientation of the current with respect to threading dislocations has recently been explained 2 . Also, the low temperature behavior of measured electron mobilities in heavily doped GaN has been found to be consistent with transport associated with band-tailing and impurity-band states 3 . Exceptions aside, interpretations of GaN Hall data often make the assumption that μ Hall = μ drift , independent of temperature 4 , 5 . One consequence of assuming a unity Hall factor is that total donor concentrations and donor activation energies extracted from electrical measurements are difficult to analyze with confidence.

We present new calculations of the Hall factor, rH, for n-type wurtzite GaN that predict a monotonic decrease with increasing temperature for fixed donor concentrations. The presentation of calculated electron concentration, Hall factor, and drift mobility is followed by a comparison to relaxation-time approximation results and by a discussion of the Hall factor for dislocation dominated scattering.

Hall Simulation

Due to thermal fluctuations in the calculated drift velocity components, calculations of electron drift velocity in crossed electric and magnetic fields using the Monte Carlo method were, until recently, limited to high electric and magnetic field conditions 6 . Advances in computation speed and in the transport simulation algorithm have enabled the application of the Monte Carlo method to weak field conditions. Hall simulations have been reported for a variety of semiconductor materials: Si and Ge 7 , InSb and other highly polar semiconductors 8 , and GaAs 9 . For GaN, electron drift velocities for both large and small static electric fields have been calculated in this manner and low field drift mobility expressions suitable for device modeling have been developed for temperatures above 300K 10 .

The present Hall factor calculations use the same approximate GaN analytic band structure which has been published previously 10 , although for the present study of low field transport only the parabolic part of the central conduction band valley is relevant. The electric field strengths are chosen small enough to ensure that ohmic conduction is modeled, and the magnetic field is held in the low field regime such that μ drift |B| = 1.5 × 10−2. For instance, at liquid nitrogen temperature and ND = 1017 cm−3, the electric field in the linear regime is taken to be less than 0.5 kV/cm which results in μ = 2490 cm2 V−1sec−1 and, consequently, |B| = 0.06 Tesla. The free flight motion of the electrons in crossed electric and magnetic fields is carried out in the usual way for Monte Carlo simulations 6 . The scattering parameters are the same as those in reference 10 and references therein. We favor acoustic deformation potential scattering with energy dissipation over the elastic approximation. This treatment of the acoustic phonon scattering prevents the simulated electron distribution from artificially skewing towards higher energies as would otherwise occur when polar optical phonon scattering is negligible. In addition, neutral impurity scattering is included in our calculations. The effects of neutral impurities become important at low temperatures in material that is not dominated by dislocation scattering.

In Figure 1, the electron concentration is plotted as a function of temperature using the screened shallow donor binding energy results of Wang, et al for different background electron concentrations 11 . The temperature effects on the carrier screening are taken into account for this calculation. The figure clearly shows that for uncompensated GaN at liquid nitrogen temperature there is a large fraction of neutral donor impurities present.

Figure 1: Electron concentration vs. temperature for GaN with donor doping densities of 5 × 1016 cm−3 (dash-dotted), 1017 cm−3 (solid), and 5 × 1017 cm−3 (dashed).

Scattering due to neutral donors is adapted for the Monte Carlo method from the relax-ation-time result of Meyer and Bartoli 12 . An analytic expression for the total scattering rate which is valid over the entire energy range of electron energies can be derived from ref. 12. The neutral impurity scattering rate in terms of the total collision cross section, σ c , is written as:

(1)

where

(2)

Nx is the concentration of neutral donor scattering centers, v is the electron velocity, ED is the donor binding energy, w = (aok)2 = E/ED, and ao = e2/2 κoED which defines an effective Bohr radius. Since scattering of low energy electrons by neutral impurities is velocity randomizing, the collision cross section differs from the momentum-transfer cross section of Meyer and Bartoli 12 only for incident particles with energies much greater than the donor binding energy. The Born approximation is valid for fast electrons where the term in {} of equation (2) is dominant. The collisions are treated as elastic and exchange effects are neglected.

In the Monte Carlo simulation, the final scattering states are determined using the angular distribution derived from the Born scattering amplitude given by 13

(3)

where θ is the angle between initial and final k-vectors. It is easily shown that the normalized probability distribution of scattering angles is

(4)

A flat angular distribution results in the low energy(w→0)limit. This is consistent with the velocity randomizing behavior for low electron energy implicit in the phenomenological collision cross section of equation (2).

Results and Discussion

Figures 1 and 2 display the calculated free carrier concentrations and electron drift mobilities for fixed doping concentrations. For uncompensated GaN, Figure 1 indicates that the ionized impurity concentration at 77K is a small fraction of that at 300K. In Figure 2, the temperature dependent drift mobility is plotted for ionized impurity concentrations corresponding to Figure 1. The dashed curves at low temperature indicate the analytical relaxation-time approximation 14 to the drift mobility including all of the scattering mechanisms except polar optical phonon scattering. The decrease in mobility at low temperature is caused in part by neutral impurity scattering. For the lowest doping concentration considered in this calculation, ND = 5 × 1016 cm−3, we find that the neutral impurity scattering plays a large role at low temperature because of the significant carrier freeze-out evident from Figure 1. At 77K, the relaxation-time approximation result for the mobility is 16% higher than that shown in Figure 2 without the inclusion of neutral impurity scattering effects. The low temperature mobility obtained by Monte Carlo simulation which includes energy dissipation in the acoustic deformation potential scattering is slightly larger than that obtain from the analytical relaxation-time model which uses an elastic approximation.

Figure 2: Electron drift mobility vs. temperature for GaN with donor doping densities of 5 × 1016 cm−3 (top line), 1017 cm−3 (middle), and 5 × 1017 cm−3 (lowest). The dashed lines indicate the relaxation-time approximation at low temperature.

If the mobility at 77K was dominated solely by ionized impurity scattering such that μ ∼T 3/2/n I , the decrease in mobility at low temperatures would not be as great or would completely vanish for uncompensated GaN. Conversely, a greater decrease in drift mobility is evident for heavily compensated GaN since the ionized impurity concentration remains substantial at low temperatures 15 , 16 . Our calculations indicate that having an acceptor concentration of 5 × 1016 cm−3 reduces the 77K mobility for GaN with 1017 cm−3 donors to approximately 800 cm2 V−1s−1, or by 30%. Further reductions in mobility from dislocation scattering have been predicted and recently measured by Ng, et al 17 for lateral drift normal to charged dislocations in n-type GaN films.

The Hall factor determined from our Monte Carlo simulation is plotted in Figure 3 with the error bars indicating the maximum deviation from the average steady-state value. In this context, steady-state was established by simulating 107 collisions. The number of collisions was chosen sufficiently large such that the estimated electron drift velocity perpendicular to the applied electric field can be resolved above the thermal fluctuations. Fluctuations in the Monte Carlo drift velocity are inversely proportional to the square-root of the number of collisions simulated. In all three plots of Figure 3, rH decreases with increasing temperature to a value approaching 1.2. A similar downward trend from 1.5 toward 1.2 at room temperature has been calculated using a variational technique 15 . The increase in rH at low temperatures is due to dominant impurity scattering. We also expect that rH should exhibit a similar trend at low temperature for Hall measurements of mobility perpendicular to the growth direction in samples where scattering from threading dislocations is dominant.

Figure 3: Hall factor vs. temperature for fixed doping concentrations from Monte Carlo simulations.

When the mobility is limited solely by dislocation scattering and this process can be modeled by scattering from a line charge, rH can be estimated by computing the average momentum relaxation time. It follows from the analytic expression given in reference 2 that,

(5)

Here, K 2 is the second order modified Bessel function and is the Debye screening length. In Figure 4, rH as obtained from equation (5) is plotted for drift perpendicular to threading dislocations as a function of temperature. The free carrier density used to determine the screening length is that shown in Figure 1 for uncompensated GaN. The region of interest in Figure 4 is at low temperatures where dislocation scattering may be dominant 17 . Specifically, rH is smaller than but close to 1.93 (the limit). From this estimate, it is reasonable to expect that when electron mobility is dominated by dislocation scattering, the Hall factor will tend to values larger than the room temperature value – perhaps to an extent greater than the Monte Carlo simulations for material without dislocations predict in Figure 3.

Figure 4: Hall factor resulting from dislocation scattering vs. temperature with constant donor doping densities of 5 × 1016 cm−3 (dash-dotted), 1017 cm−3 (solid), and 5 × 1017 cm−3 (dashed).

In summary, Monte Carlo simulations of electron transport in crossed electric and magnetic fields are used to compute mobility and Hall factor values at various temperatures and doping concentrations. This technique is particularly desirable for simulating Hall transport governed by polar optical phonon scattering and we show that rH should be approximately 1.2 for temperatures above 300K. At lower temperatures rH deviates farther from unity, and the values obtained by simulation are consistent with those obtained from classical analytical approximations. We have also addressed Hall measurements of lateral mobility, which can be dominated by scattering from threading dislocations when dislocation densities are high. In those cases, we predict that rH will remain significantly larger than unity at low temperature, especially for low carrier concentrations where screening is weak.

Acknowledgments

This work benefited from useful discussions with J.R. Meyer. A grant from the Minnesota Supercomputer Institute is gratefully acknowledged. This work was supported in part by Hughes Research Laboratories and by NSF under grant number ECS-9811366. The work at Georgia Tech was sponsored in part by NORTEL through contract E21-5D7 and by ONR through contract E21-K19.

References

1 Look, D.C. and Molnar, R.J., Appl. Phys. Lett., 70, 3377 (1997).
2 Weimann, N.G., Eastman, L.F., Doppalapudi, D., Ng, H.M., and Moustakas, T.D., J. Appl. Phys., 83, 3656 (1998).
3 Molnar, R.J., Lei, T., and Moustakas, T.D., Appl. Phys. Lett., 62, 72 (1993).
4 Götz, W., Romano, L.T., Walker, J., Johnson, N.M., and Molnar, R.J., Appl. Phys. Lett., 72, 1214 (1998).
5 Götz, W., Johnson, N.M., Chen, C., Liu, H., and Imler, W, Appl. Phys. Lett., 68, 3144 (1996).
6 Boardman, A.D., Fawcett, W., and Ruch, J.G., Phys. Stat. Sol., 4, 133 (1971).
7 Dijkstra, J.E. and Wenckebach, W. Th., Appl. Phys. Lett., 70, 2428 (1997).
8 Warmenbol, P., Peeters, F.M., and DeVreese, J.T., Phys. Rev. B, 33, 1213 (1986).
9 Williams, C.K., Littlejohn, H.A., Glisson, T.H., and Hauser, J.R., Superlattices and Microstructures, 2, 201 (1986).
10 Albrecht, J.D., Wang, R.P., Ruden, P.P., Farahmand, M., and Brennan, K.F., J. Appl. Phys., 83, 4777 (1998).
11 Wang, R.P., Ruden, P.P., Kolnik, J., Oguzman, I., and Brennan, K.F., Mat. Res. Soc. Symp. Proc., 445, 935 (1997).
12 Meyer, J.R. and Bartoli, F.J., Phys. Rev. B, 24, 2089 (1981).
13 See, e.g., Schiff, L., Quantum Mechanics, McGraw-Hill, New York (1955).
14 For a canonical list of analytical relaxation-time approximations see Look, D.C., Electrical Characterization of GaAs Materials and Devices, John Wiley & Sons, New York (1989).
15 Look, D.C., Sizelove, J.R., Keller, S., Wu, Y.F., Mishra, U.K., and DenBaars, S.P., Solid State Commun., 102, 297 (1997).
16 Wickenden, A.E., Gaskill, D.K., Koleske, D.D., Doverspike, K., Simons, D.S., and Chi, P.H., Mat. Res. Soc. Symp. Proc., 395, 679 (1996).
17 Ng, H.M., Doppalapudi, D., Moustakas, T.D., Weimann, N.G., and Eastman, L.F., Appl. Phys. Lett., 73, 821 (1998).