Research Article  Open Access
Analytic Model for Conduction Current in AlGaN/GaN HFETs/HEMTs
Abstract
We have developed a new, zonebased compact physicsbased AlGaN/GaN heterojunction fieldeffect transistor (HFET) model suitable for use in commercial harmonicbalance microwave circuit simulators. The new model is programmed in VerilogA, an industrystandard compact modeling language. The new model permits the dc, smallsignal, and largesignal RF performance for the transistor to be determined as a function of the device geometric structure and design features, material composition parameters, and dc and RF operating conditions. The new physicsbased HFET model does not require extensive parameter extraction to determine model element values, as commonly employed for traditional equivalentcircuitbased transistor models. The new model has been calibrated and verified. We report very good agreement between simulated and measured dc and RF performance of an experimental Cband microwave power amplifier.
1. Introduction
AlGaN/GaN heterojunction fieldeffect transistors (HFETs) are promising RF transistors for use in highpower and highfrequency circuit applications. These HFETs possess a combination of high current density capability and high breakdown voltage due to the desirable physical properties of the materials, such as high critical electric field for breakdown, high electron mobility and saturated carrier velocity, high carrier density in the channel, lower dielectric constant compared to the conventional materials, and high thermal conductivity. These parameters permit the HFET to operate at high RF voltage and current, which results in high power operation at high frequency [1, 2]. The technology for fabricating devices and circuits in AlGaN/GaN is developing rapidly [3–6] and this rapid development is creating a need for improved device models. RF power amplifiers based on AlGaN/GaN HFETs are now commercially available from several companies, including Nitronex, RFMD, and TriQuint. However, to date, no commercially available HFET model for use in harmonicbalance circuit simulators exists that can predict the largesignal RF operation of an HFET or an MMIC before the active device is fabricated, characterized, and fitted.
The basic structure for an HFET is shown in Figure 1. The dc and RF performance of the transistor varies with the physical dimensions and material properties. However, the sensitivity of RF power performance to the physical parameters will vary, depending upon the particular parameter, and variations in some parameters (e.g., the gate length, ) have a more significant effect upon device performance than others. These specific parameter sensitivities are not easy to determine. Consequently, MMIC designers cannot consider the physical parameters of the device when designing circuits. They can only simulate circuits containing transistors with known and defined equivalent circuit compact models, which preclude the use of harmonic balance simulators for use in device optimization. In this report, we introduce an analytic HFET model that permits the codesign and optimization of active HFETs and passive elements in an MMIC environment, which will enhance and speed integrated circuit development.
Reported HFET dc models include empirical models [7, 8], twodimensional multiphysics solvers [9, 10], and compact physicsbased models [11–16]. Compact equivalent circuit models are suitable for simulating HFET circuits, but not before the HFET devices have been fabricated. Twodimensional solvers can predict the dc IV characteristics of an HFET and even its smallsignal RF parameters, but they are difficult to employ in real time within a harmonic balance circuit simulator since runtime interpolation of a database of precomputed solutions is cumbersome [17, 18]. Compact physicsbased models, however, can run in harmonic balance solvers because they are analytic and therefore sufficiently efficient in computation time requirements that they can predict the operation of an RF HFET under largesignal RF drive conditions. Previous work has demonstrated the facility and accuracy of this approach [19, 20]. Unfortunately, however, this capability is not generally available, because the previously reported models could not be readily ported to commercial simulators. This report addresses that deficiency.
In this work, we report for the first time the implementation of a physicsbased compact HFET model in VerilogA, and we used Microwave Office (MWO) to compare its predictions against the RF power performance of an experimental HFET. MWO is an EDA package that is available from Applied Wave Research (AWR) Corporation. Our firstgeneration VerilogA HFET model is a single, selfcontained VerilogA module.
The new model is developed, based upon separating the conducting channel of the HFET into a series of zones, based upon operational physics [21]. The HFET model is generalized to suit different fabrication process by introducing a curvature parameter in the vE relationship [22] for electrons in the conducting channel of the HFET in order to control the sharpness of the knee of the IV characteristics. We also considered the electric field decay region in the drain access for triode operation to ensure continuity of the terminal characteristics. We modified the pinch off voltage to take into account different Al mole fractions in the AlGaN barrier layer in the AlGaN/GaN HFET. Finally, it is necessary to consider channel breakdown in order to accurately simulate Power Added Efficiency (PAE) at high RF power operation. We have written the model in VerilogA language and implemented it in MWO. The model is verified by comparison of simulated and measured dc and RF largesignal performance for an AlGaN/GaN HFET SBand amplifier. Good agreement between the model simulation and the measurement data are obtained.
2. Zone Model of the AlGaN/GaN HFET
As previously indicated, the model is formulated based upon separating the conducting channel of the HFET into a series of five zones, which are defined based upon the physics that dominate within each zone. The model operates in two modes, triode and saturation. The transition between the two modes is smooth and dependent upon device design and operation criteria. Figure 2 shows the salient features of the zone model in its two operating modes. The typical electron is assumed to remain in the same longitudinal crosssection of the device on its path from source to drain. This defines the  plane, where is measured from the source electrode and is measured down from the AlGaN/GaN interface. For each operating mode, the path between source and drain is segmented into a few contiguous intervals. In each interval, the physical operation is determined by twodimensional numerical simulations and this information is used to develop an analytic physical model for that particular zone. The 2D simulations reveal the fundamental operation of each zone and permit the simplified analytic physical model to be derived. At the boundaries between adjacent zone intervals, the physical operation changes and the zones are interfaced by enforcing continuity of the electric potential values and derivatives.
(a)
(b)
Associated with each interval, we define its zone as the interval itself, initialization of the distance, potential, and electric field . We also define a rule for terminating the interval, and a sequence of operations to compute at the endpoint and within it, as appropriate.
When the HFET is in triode operation, it can be modeled with four zones. In saturated operation, the model requires five zones. For each operating mode, the terminal characteristics of the device must be consistent with a simultaneous solution of all the zones that exist in that mode.
Fortunately, it is possible to compute this simultaneous solution efficiently by setting up equations for the zones from left to right in sequence, as we will now show. Each zone is solved in three steps. First, the distance, potential, and electric field parameters are initialized at its lefthand boundary by applying its lefthand initialization rule to the final value of the triple from the preceding zone. Second, the nominal model in the zone is used to advance the parameters until the termination condition is detected. Third, the final values of are recorded and are available for the subsequent zone and elsewhere in that evaluation of the controlled voltage source.
In either of the two operating modes, the terminal characteristics can be qualitatively related to the physics of each zone. At zero draintosource bias , the channel consists of three zones, the source and drain access regions and the region beneath the gate. We denote the source and drain access regions as the source neutral zone (SNZ, or Z1) and the drain neutral zone (DNZ, or Z5), respectively, to emphasize their net Coulombic neutrality which results from the approximate equality of the sheet charge density of the 2DEG in these two regions and the polarization charge density. We denote the region beneath the gate as the intrinsic FET zone (IFZ, or Z2).
For positive , a chargedeficit zone (CDZ, or Z4) forms in the drain access region close to the gate edge. The positive net charge in this partially depleted region (Z4) smoothly reduces the magnitude of the electric field, in accordance with Poisson’s equation. In the triode mode, the length of the Z4 is short and the voltage drop across it is typically less than a volt, but it is required for continuity of at the transition to saturated operation at . At fixed gate bias, the length of the CDZ zone increases with drain bias, which reduces the length of the DNZ zone, since the sum of their lengths equals the gatetodrain spacing, .
For saturated operation, . Electrically in the  plane, the transition from triode to saturation operation occurs at the knee of the IV curve, where the slope of the curve flattens. Physically, saturated operation begins when the magnitude of the longitudinal electric field in the channel at the drainside gate edge (denoted as ) first exceeds the critical field, , which effectively pinches the 2DEG off by flat banding the quantum well.
With increasing , the electric field () at the drainside edge of the gate continues to increase and the location where first exceeds moves toward the source. Under the gate, we denote the location where as and the voltage there as . In the region , electrons are repelled away from the AlGaN/GaN interface to form the velocitysaturated spacechargelimited zone (SLZ, or Z3). This SLZ extends longitudinally to the gate edge and typically extends down to the substrate because the gate repels electrons when , so that the gradual channel Approximation (GCA) fails and the magnitude of the field increases rapidly with distance, .
With increasing , the length of the SLZ zone increases, which reduces the length of the IFZ zone under the gate since the sum of their lengths is constrained to be . In the drain access region, the length of the CDZ zone increases simultaneously with increasing , which reduces the length of the DNZ zone, since the sum of their lengths is constrained to equal the gatetodrain distance.
The gradual channel approximation (GCA) is readily adapted to zones Z1, Z2, and Z5. In zone Z3, the GCA fails because the carrier trajectories are not one dimensional as discussed above. In zone Z4, the 2DEG has reestablished itself, but the sheet charge density of the 2DEG is insufficient to neutralize the fixed polarization sheet charge and electron velocity in this 2DEG is effectively saturated.
To the best of our knowledge, this charge deficit zone (or zone Z4) is unique to AlGaN/GaN HFETs, but it can dominate device operation over part of the RF cycle. In zone Z4, the 2DEG is stable but is incompletely filled. The length of zone Z4 is approximately proportional to the difference , but increases with until when impact ionization in the channel produces unacceptable breakdown effects at the device terminals. In GaN, the breakdown electric field is about two orders of magnitude larger than the electric field at which velocity saturation begins. Therefore, the maximum length of this region is proportional to  which is exceptionally large for GaN. Electrons enter zone Z4 with their velocity saturated and the transistor effect of the gate region limits the current almost independently of the local electric field, even when it is two orders of magnitude larger than , as it can be in GaN before its breakdown occurs at very high . This combined effect keeps zone Z4 partially depleted, so that at fixed gate bias, the voltage dropped across zone Z4 increases approximately quadratically with its length, which adjusts itself to gradually reduce from to , where the fielddependent mobility of the electrons is sufficiently high to reestablish net Coulombic neutrality at the given . Over the majority of a typical RF cycle, most of drops across the CDZ zone.
3. Generalized VelocityField Relationship for Carrier Electrons
Theoretical investigations of electron dynamics using Monte Carlo techniques have determined velocityfield characteristics associated with GaN materials. These theoretical simulations show that the electron drift velocity initially increases with the applied electric field but reaches a peak value, after which it gradually decreases to a saturated value at a high electric field [23]. The peak value decreases with increasing doping concentration and temperature [24]. Some vE relationships for electron transit in a 2DEG for the AlGaN/GaN structure have been reported [23, 25]. However, no evidence of a velocity overshoot is apparent in the terminal characteristics of devices that we have considered. In fact, an equilibrium vE characteristic has proved adequate to accurately simulate experimental results and we find that (1) accurately simulates the dc and RF currents that flow in experimental AlGaN/GaN HFETs.
The curvature of the transition from the lowfield region to the highfield region is important in calculating the knee region of the dc currentvoltage relationship for the device. We have found that the same vE model (1) accurately simulates both the dc and largesignal RF operation. Consequently, we have chosen (1) for the velocityfield relationship because its transition is smooth but adjustable. Figure 3 also compares the vE relationship in (1) to the twofield relationship.
The vE relationship we adopt is shown as where we define as the negative of the usual longitudinal electric field and we regard it as a function of distance from the source electrode. In (2), is the lowfield electron mobility: is the critical electric field marking the onset of the highfield region: is the asymptotic saturated velocity, and controls the curvature of the knee of the vE curve. For the experimental data that we have analyzed, this curvature parameter has varied over the interval , so that the critical electron velocity depends on as well as . In our investigations, the values of the , , parameters can be determined by adjusting the estimated physicalbased values so that the simulated terminal IV characteristics match the measured IV data for an HFET.
Figure 3 shows how controls the curvature of the knee of (1) without changing the lowfield and highfield regions of the vE model. In the limit of large , this vE curve approaches a piecewise linear vE model. For the choice , this vE curve approaches the twofield model with the same parameters as in [15].
We estimate the saturated velocity , the lowfield mobility , and the rate of transition from to from simultaneous fits of dc and RF measurements when we use the velocityfield relationship in (1). The initial value for the electron mobility at low electric field is generally determined from measured Hall mobility data. This value may be slightly varied in order to accurately simulate the measured dc IV data in the linear region in some cases but does not significantly vary from the Hall values, and often the measured Hall value is found to produce excellent results.
4. PinchOff Voltage
The conducting channel in an AlGaN/GaN HFET is formed from the 2DEG just below the interface of the AlGaN barrier layer that is grown on a GaN layer. This 2DEG conducting channel is formed by spontaneous and piezoelectric polarization effects at the AlGaN/GaN interface [26]. Figure 2 shows the crosssectional view of an HFET. The sheet charge density of this 2DEG channel is determined by the aluminum percentage and the thickness of the AlGaN layer. The sheet charge density can be modulated by the deposition of a gate electrode on the AlGaN surface and by the electric voltage that is applied to the gate electrode. When is sufficiently negative, the channel vanishes everywhere under the gate and, in particular, at the sourceside edge of the gate, which determines the pinchoff voltage of the AlGaN/GaN HFET. The pinchoff voltage is conventionally written as a function of mole fraction of aluminum in the Al_{m}Ga_{1−m}N, but also depends on the effective thickness of the AlGaN barrier, its Schottky barrier height , its electric permittivity , and its doping , as well as the conduction band offset at the AlGaN/GaN interface and the combined piezoelectric/spontaneous polarization sheet charge . In (5), the dielectric constant and bandgap of AlGaN are expressed as [14] and the Schottky barrier is expressed as and the band offset is given by the expression:
5. Physics of the Zones
5.1. Zone Z1, the Source Neutral Zone, or SNZ
In zone Z1 (and also in zone Z5), the current at location is conducted by the 2DEG, but for steadystate operation, cannot depend on , so that where is the drain current, is the gate width, is the fundamental charge, and is the local electron density of the 2DEG. In the source access region, is fixed at to neutralize , so that (9) and (1) imply that is constant with respect to , where is the value of that constant at a given and is a convenient scale factor for for a given HFET.
The electric field typically occurs near the drain electrode in zone Z5, the other neutral zone of the HFET. In contrast to the fixed length of zone Z1, the length of zone Z5 is state dependent. Zone Z5 begins when the lateral field in zone Z4 has finally diminished to , which terminates zone Z4. Within both the source and drain neutral zones, electron transport is identical and sets the drift velocity in both zones Z1 and Z5.
5.2. Zone Z2, the Intrinsic FET Zone, or IFZ
At the sourceside gate edge, electrons leave zone Z1 and enter zone Z2. In zone Z1, the voltage on the upper AlGaN surface follows the voltage of the 2DEG at the AlGaN/GaN interface. In zone Z2, however, the gate electrode holds the surface voltage at .
In accordance with the GCA, we treat this gate region as an MIS capacitor and write the electron sheet density: in terms of the effective gate voltage, the effective gatechannel capacitance , and the effective thickness of the AlGaN barrier layer.
Electrons in zone Z2 drift toward the drain because increases with in the channel, but this drift cannot persist past where has risen enough to exceed . In triode operation, the channel fills the entire gate region. Following the GCA, we substitute (11) into (9) and use (2) to change the variable of integration from to , so that in triode operation we find where is the physical gate length, is the voltage at the sourceside gate edge, is the voltage at the drainside edge of the gate in triode operation, and and are the lengths of drain access and zone Z4.
If while , then the device is in triode mode; if before , then the device is in saturated mode. In triode operation, zone Z3 does not affect the parameters because its thickness vanishes, zero voltage is dropped across the zone, and the electric field does not change. In saturation, zone Z3 is significant because when .
In saturated operation, and are still defined at the drainside gate edge, but (12) is not valid there because the longitudinal electric field exceeds before exiting the gate region so that and . In this case, the length of zone Z2 is less than , and we write where The transition of terminal characteristics from triode mode to saturated mode occurs when the drain current is high enough to make and at the gate edge. This is determined by using (14) for the upper limit of (12) and solving for .
5.3. Zone Z3, the SpaceChargeLimited Zone, or SLZ
The SLZ zone only occurs when the device enters saturation. We approximate the onedimensional Poisson equation: where and both vary with . Comparison with twodimensional finiteelement simulations suggests that the product is nearly constant: so that we can define the average electric field derivative: Integrating Poisson’s equation once yields and a second integration yields we can include zone Z4 in the sum over all five zones of voltage increments for each zone to get .
In the triode mode, we can define in the same sum. To reexpress this result as the controlled current source for given and , we invert the sum () iteratively to find the corresponding to a given .
5.4. Zone Z4, the Charge Deficit Zone, or CDZ
Zone Z4 occurs in both triode and saturated modes when as defined by (3). Since the surface voltage in zone Z4 is not clamped by the gate electrode, it rises with the channel potential and the onedimensional Poisson equation, can be used to obtain assuming that and . The voltage and electric field continuity conditions lead to the following two equations:
5.5. Zone Z5, the Drain Neutral Zone, or DNZ
Within zone Z5, the nominal electron physics are identical to that in zone Z1 but the zone length varies dynamically and satisfies different boundary conditions. As in zone Z1, overall charge neutrality prevails in zone Z5 and the current is described by (9) and (10), so the lateral electric field in zone Z5 is constant and coincides with Es in zone Z1. In the preceding zone Z4, but the electric field continuously varies from magnitude to . In zone Z5, the electron transport model is similar as in zone Z1, the SNZ, but the length of zone Z5 varies dynamically. The lateral electric field in zone Z4 varies and zone Z5 begins when the lateral electric field in zone Z4 is reduced to its value in zone Z1.
It may be possible to fabricate an HFET with sufficiently short drain access region and to bias it at high enough drain voltage to deplete the entire drain access region, but we have not observed this in any device we have considered. Currently, we treat this possibility as an error condition.
The physicsbased compact HFET model described above has been written in the VerilogA language and implemented in MWO. The flowchart in Figure 4 shows the process to calculate the drain current. The model was verified by comparing the simulated and measured dc and largesignal RF performance for an HFET Sband amplifier. The parameters used in the model are listed in Table 1.

Note that the model requires twenty (20) input parameters, as compared to the eighty (80) parameters required to fit conventional equivalentcircuitbased models. The input parameters consist of physical dimensions, doping and thickness levels, and so forth, rather than equivalent circuit element values determined from parameter extraction procedures. The new model is, therefore, easier to define than traditional equivalent circuit transistor models.
6. Simulation in Microwave Office
To calibrate the model, we have compared the simulation with the measurement data for AlGaN/GaN HFETs of 0.8 μm gate length and 400 μm gate width. The curvature parameter in the vE relationship is adjusted to be 1.45 for this device. The polarization sheet charge density is 7.1e12 cm^{−2}.
6.1. DC IV Curves Simulated by MWO
The model is implemented in the commercial circuit simulator Microwave Office developed by AWR Corporation. The results are compared with the experimental measurements of an .8 μm HFET. Figure 5 shows excellent agreement between the simulated results and measurement data.
The source and drain access regions, zones Z1 and Z5, introduce an extrinsic resistance to the intrinsic HFET and will affect the performance of the HFETs. It has been previously shown that under highcurrent operation and largesignal RF drive, these resistances become nonlinear as spacechargelimited (SCL) current transport conditions are approached [20]. Both the source and drain resistances become nonlinear and increase with channel current. The effect is most apparent in the source region since the voltage drop across the source access region subtracts directly from the applied and limits the drain current of the HFET considerably before increases to . The drop across the drain access region lowers the slope of in the triode region and increases the drain bias required to achieve . In power microwave applications, using a constant resistance for either will typically overestimate the drain current at the upper left limit of the dynamic load line and with it, the output power. From (10) and the length of the source access region, we find where is the voltage drop across the SNZ zone. The drain access region can be treated similarly but both zones Z4 and Z5 must be considered. Figure 6 shows the nonlinear source resistance as a function of the chancel conductance current with a source access region of 1.2 μm. The values of beta are 2, 1.5, and 1.2, respectively.
6.2. Power Performance Simulated by MWO
An AlGaN/GaN HFET power amplifier is simulated using a singletone power sweep. The amplifier was biased at a drain bias of and operated in class B at a fundamental frequency of GHz. Three harmonics were used in the simulation and the corresponding three pairs of harmonic impedances for both input and output sides are tuned to fit the experimental measurement. In the initial simulation, we did not take into account the reactance components of the device, so for the harmonic impedances tuning, we also excluded the effects of the reactance. Figure 7 shows the simulated largesignal RF performance compared to the measured data, including output power, power gain, and PAE. The inset of Figure 7 shows the results without considering the channel breakdown effects. As shown in the inset, we observe that when input power increases to around 13 dBm, there is an increasing discrepancy between the simulated and measurement RF performance, most apparent in PAE, which increases with increasing input power. PAE becomes unacceptably higher than the measured PAE, due to the onset of channel breakdown in the experimental device. Neglecting breakdown effects causes the model to dramatically overestimate the PAE under large RF drive conditions. When the draintosource breakdown effect is considered, the results agree very well with the measured data for the input power ranging from −8 dbM to 20 dBm.
As is evident in Figure 7, channel breakdown is crucial for understanding the largesignal RF performance, and particularly the maximumpoweradded efficiency of a microwave power amplifier. For sufficiently large terminal voltage, , impact ionization occurs in the conducting channel of the HFET and the avalanche generated current opposes further increase of . If the temperature rise within the device is controlled and limited, channel breakdown is reversible over the RF cycle and the device is not damaged. RF breakdown is, in fact, a fundamental phenomenon that limits the RF performance of these devices under high power operation. Breakdown in AlGaN/GaN transistors occurs in a more gradual manner than in GaAs MESFETs. That is, avalanche ionization occurs in the conducting channel with the generated electrons adding to the channel current, thereby creating a gradual increase in channel current since ionization only occurs during the peak voltage period of the RF cycle when the channel field briefly rises above the breakdown electric field. Since this only occurs during a short period of the RF cycle, a pulse of generated charge occurs, which limits the ionizationgenerated current increase. Gate breakdown is the dominant RF breakdown mechanism in traditional GaAs MESFETs [27] but does not generally occur in AlGaN/GaN HFETs. Experimental data reveals that gate leakage does occur in AlGaN/GaN HFETs [28], but the leakage is caused by electron tunneling from the gate electrode and conduction occurs along the surface area of the HFET between the gate and drain electrodes. Gate leakage affects HFET reliability, introduces transient effects, and is explained elsewhere [29].
6.3. Maximum PAE and Breakdown Resistance
Figure 8 shows the dc currentvoltage characteristics with and without breakdown. As indicated, once breakdown occurs, the increase in channel current is gradual with drain voltage.
The effect of the breakdown voltage on largesignal PAE is shown in Figure 9. Increases in the channel breakdown voltage will improve PAE. Figure 10 shows the effect of the breakdown resistance, defined as the slope of the dynamic IV characteristic once breakdown occurs. For a breakdown voltage of 39 V, when input powers are below 12 dBm, the breakdown resistance does not significantly affect the PAE. After the input power exceeds 12 dBm, the PAE increases with increasing resistance.
6.4. Maximum PAE and Load Resistance at Fundamental and Second Harmonic Frequency
The harmonic impedances of the matching network affect the large signal outputs. The effect of the impedances at the first and second harmonics is shown in Figures 11 and 12.
7. Conclusion
A new physicsbased compact model for AlGaN/GaN HFETs has been developed. The new model is based upon separating the conducting channel into various zones determined by the physical operation principles that dominate in each zone. A simplified set of the semiconductor device equations are applied to each zone and the zones are then interfaced by forcing electric field and potential continuity at the zone interfaces. In this manner, a complete model for the device can be constructed. The compact model is suitable for integration into readily available harmonicbalance circuit simulators. The model has been formulated in VerilogA and is integrated into the harmonicbalance simulator Microwave Office offered by AWR. However, the new compact model can be employed in any readily available circuit simulator.
The new model permits the dc and RF performance of the device to be determined as a function of device design parameters such as material parameters, the physical dimensions and doping concentrations of the various layers, and the charge transport characteristics. About twenty parameters are required to define the model, but most of these consist of handbook data. The new model does not require the parameter extraction procedures commonly employed to determine equivalentcircuitbased models. The new model, therefore, requires significantly less effort to define.
The new model has been calibrated and verified by comparison with measured data for experimental AlGaN/GaN HFET amplifiers operating at Sband. Excellent agreement between simulated and measured dc and largesignal RF data is obtained.
Acknowledgment
This work was supported by ARO grant DAA0190310148 and by ONR grant N0014050419.
References
 J. C. Zolper, “Wide bandgap semiconductor microwave technologies: from promise to practice,” in Proceedings of the IEEE International Devices Meeting (IEDM '99), pp. 389–392, December 1999. View at: Google Scholar
 R. J. Trew, “SiC and GaN transistorsIs there one winner for microwave power applications?” Proceedings of the IEEE, vol. 90, no. 6, pp. 1032–1047, 2002. View at: Publisher Site  Google Scholar
 M. A. Khan, A. Bhattarai, J. N. Kuznia, and D. T. Olson, “High electron mobility transistor based on a GaN—AlxGa (1x)N heterojunction,” Applied Physics Letters, vol. 63, no. 9, pp. 1214–1215, 1993. View at: Publisher Site  Google Scholar
 Y. F. Wu, B. P. Keller, P. Fini et al., “High Alcontent AlGaN/GaN MODFET's for ultrahigh performance,” IEEE Electron Device Letters, vol. 19, no. 2, pp. 50–53, 1998. View at: Google Scholar
 W. Lu, J. Yang, M. A. Khan, and I. Adesida, “AlGaN/GaN HEMTs on SiC with over 100 GHz fT and low microwave noise,” IEEE Transactions on Electron Devices, vol. 48, no. 3, pp. 581–585, 2001. View at: Publisher Site  Google Scholar
 Y. Cordier, J. C. Moreno, N. Baron et al., “Demonstration of AlGaN/GaN highelectronmobility transistors grown by molecular beam epitaxy on Si(110),” IEEE Electron Device Letters, vol. 29, no. 11, pp. 1187–1189, 2008. View at: Publisher Site  Google Scholar
 I. Angelov, H. Zirath, and N. Rorsman, “A new empirical nonlinear model for HEMTdevices,” in Proceedings of the IEEE MTTS International Microwave Symposium Digest Part 3, pp. 1583–1586, 1992. View at: Google Scholar
 I. Angelov, L. Bengtsson, and M. Garcia, “Extensions of the chalmers nonlinear HEMT and MESFET model,” IEEE Transactions on Microwave Theory and Techniques, vol. 44, no. 10, pp. 1664–1674, 1996. View at: Google Scholar
 Silvaco Manual, Silvaco Data Systems Inc., Santa Clara, Calif, USA.
 A. Asgari, M. Kalafi, and L. Faraone, “A quasitwodimensional charge transport model of AlGaN/GaN high electron mobility transistors (HEMTs),” Physica E, vol. 28, no. 4, pp. 491–499, 2005. View at: Publisher Site  Google Scholar
 H. Ahn and M. El Nokali, “Analytical model for high electron mobility transistors,” IEEE Transactions on Electron Devices, vol. 41, no. 6, pp. 874–878, 1994. View at: Publisher Site  Google Scholar
 R. Singh and C. M. Snowden, “A quasitwodimensional HEMT model for DC and microwave simulation,” IEEE Transactions on Electron Devices, vol. 45, no. 6, pp. 1165–1169, 1998. View at: Google Scholar
 J. D. Albrecht, P. P. Ruden, S. C. Binari, and M. G. Ancona, “AlGaN/GaN heterostructure fieldeffect transistor model including thermal effects,” IEEE Transactions on Electron Devices, vol. 47, no. 11, pp. 2031–2036, 2000. View at: Publisher Site  Google Scholar
 T. H. Yu and K. F. Brennan, “Theoretical study of a GaNAlGaN high electron mobility transistor including a nonlinear polarization model,” IEEE Transactions on Electron Devices, vol. 50, no. 2, pp. 315–323, 2003. View at: Publisher Site  Google Scholar
 A. Koudymov, M. S. Shur, G. Simin et al., “Analytical HFET IV model in presence of current collapse,” IEEE Transactions on Electron Devices, vol. 55, no. 3, pp. 712–720, 2008. View at: Publisher Site  Google Scholar
 X. Cheng, M. Li, and Y. Wang, “PhysicsBased Compact Model for AlGaN/GaN MODFETs With CloseFormed I V and C V Characteristics,” IEEE Transactions on Electron Devices, vol. 36, no. 2, pp. 231–239, 1988. View at: Publisher Site  Google Scholar
 M. W. Shin, R. J. Trew, and G. L. Bilbro, “High temperature dc and RF performance of ptype diamond MESFET: comparison with Ntype GaAs MESFET,” IEEE Electron Device Letters, vol. 15, no. 8, pp. 292–294, 1994. View at: Publisher Site  Google Scholar
 C. W. Hatfield, G. L. Bilbro, S. T. Allen, and J. W. Palmour, “DC IV characteristics and RF performance of a 4HSiC JFET at 773 K,” IEEE Transactions on Electron Devices, vol. 45, no. 9, pp. 2072–2074, 1998. View at: Google Scholar
 M. A. Khatibzadeh and R. J. Trew, “LARGESIGNAL, ANALYTIC MODEL FOR THE GAAS MESFET,” IEEE Transactions on Microwave Theory and Techniques, vol. 36, no. 2, pp. 231–238, 1988. View at: Publisher Site  Google Scholar
 R. J. Trew, Y. Liu, G. L. Bilbro, W. Kuang, R. Vetury, and J. B. Shealy, “Nonlinear source resistance in highvoltage microwave alGaN/gaN HFETs,” IEEE Transactions on Microwave Theory and Techniques, vol. 54, no. 5, pp. 2061–2067, 2006. View at: Publisher Site  Google Scholar
 H. Yin, D. Hou, G. L. Bilbro, and R. J. Trew, “Harmonie balance simulation of a new physics based model of the AlGaN/GaN HFET,” in Proceedings of the IEEE International Microwave Symposium Digest, (MTTS '08), pp. 1425–1428, 2008. View at: Publisher Site  Google Scholar
 C. Canali, G. Majni, R. Minder, and G. Ottaviani, “Electron and hole drift velocity measurements in silicon and their empirical relation to electric field and temperature,” IEEE Transactions on Electron Devices, vol. 22, no. 11, pp. 1045–1047, 1975. View at: Google Scholar
 M. Farahmand, C. Garetto, E. Bellotti et al., “Monte Carlo simulation of electron transport in the IIInitride Wurtzite phase materials system: binaries and ternaries,” IEEE Transactions on Electron Devices, vol. 48, no. 3, pp. 535–542, 2001. View at: Publisher Site  Google Scholar
 V. M. Polyakov and F. Schwierz, “Influence of electron mobility modeling on dc IV characteristics of WZGaN MESFET,” IEEE Transactions on Electron Devices, vol. 48, no. 3, pp. 512–516, 2001. View at: Publisher Site  Google Scholar
 B. K. Ridley, “Coupled surface and channel transport in semiconductor heterostructures,” Journal of Applied Physics, vol. 90, no. 12, pp. 6135–6139, 2001. View at: Publisher Site  Google Scholar
 O. Ambacher, J. Smart, J. R. Shealy et al., “Twodimensional electron gases induced by spontaneous and piezoelectric polarization charges in N And Gaface AIGaN/GaN heterostructures,” Journal of Applied Physics, vol. 85, no. 6, pp. 3222–3233, 1999. View at: Google Scholar
 T. A. Winslow and R. J. Trew, “Principles of largesignal MESFET operation,” IEEE Transactions on Microwave Theory and Techniques, vol. 42, no. 6, pp. 935–942, 1994. View at: Publisher Site  Google Scholar
 R. J. Trew, Y. Liu, W. Kuang et al., “RF breakdown and largesignal modeling of AlGaN/GaN HFET's,” in Proceedings of the IEEE International Microwave Symposium Digest (MTTS '06), pp. 643–646, San Francisco, Calif, USA, 2006. View at: Publisher Site  Google Scholar
 R. J. Trew, D. S. Green, and J. B. Shealy, “AlGaN/GaN HFET reliability,” IEEE Microwave Magazine, vol. 10, pp. 116–127, 2009. View at: Google Scholar
Copyright
Copyright © 2012 Danqiong Hou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.