Modeling Power GaN-HEMTs Using Standard MOSFET Equations and Parameters in SPICE

The device library in the standard circuit simulator (SPICE) lacks a gallium nitride based high-electron-mobility-transistor (GaN-HEMT) model, required for the design and verification of power-electronic circuits. This paper shows that GaN-HEMTs can be modeled by selected equations from the standard MOSFET LEVEL 3 model in SPICE. A method is proposed for the extraction of SPICE parameters in these equations. The selected equations and the proposed parameter-extraction method are verified with measured static and dynamic characteristics of commercial GaN-HEMTs. Furthermore, a double pulse test is performed in LTSpice and compared to its manufacturer model to demonstrate the effectiveness of the MOSFET LEVEL 3 model. The advantage of the proposed approach to use the MOSFET LEVEL 3 model, in comparison to the alternative behavioral-based model provided by some manufacturers, is that users can apply the proposed method to adjust the parameters of the MOSFET LEVEL 3 model for the case of manufacturers who do not provide SPICE models for their HEMTs.


Introduction
Gallium nitride (GaN)-based high-electron-mobility transistors (HEMTs) are emerging as a promising power switching device for high-efficiency, high-density, and high-frequency power converters [1]. Compared to the Si counterparts, GaN-based power devices offer higher breakdown voltage and higher thermal conductivity [1]. These excellent properties enable GaN-based power devices to meet the high-performance requirements for modern power converters.
For any power-conversion application, the power circuit designers have to evaluate the impact of the power-semiconductor devices used in the circuits. From an optimization point-of-view, there is a need to reduce the number of hardware prototypes. Therefore, simulation tools are gaining more attention in the development of power electronic systems [2]. Many researchers are evaluating the high-power and high-frequency circuits with GaN-HEMT using various device models [1][2][3]. Hence, accurate modeling of GaN-HEMTs is an important research area for the simulation of power electronic systems. Various modeling approaches for GaN-HEMTs are available in the literature and are classified into four categories [1, 3,4]: Behavioral models, semiphysics-based models, physics-based models, and numerical models. Numerical modeling requires detailed information regarding internal structure, device geometry, and material properties of the GaN-HEMT and uses complex simulation tools such as SILVACO, TCAD, et cetera [5,6]. The accuracy of numerical models is very good but complexity is very high and computation is very This paper is divided into the following sections: Section 2 will introduce the GaN structure studied in this paper. Section 3 will present the selection of static equations and the method for respective parameters extraction. Sections 4 and 5 will verify the experimental static and dynamic characteristics, respectively. Section 6 will present the results of a double-pulse test (DPT) test to compare the effectiveness of the MOSFET LEVEL 3 model with the manufacturer model in LTSpice. Section 7 will conclude and foreshadow potential future work.

GaN Structure
The internal structure of a normally-off GaN-HEMT is shown in Figure 1 [22]. For HEMT fabrication on the silicon substrate, as is the case in the structure shown in Figure 1, aluminum nitride is used as a buffer layer to enable epitaxial growth of the active GaN layer. The formation of the two-dimensional electron gas (2-DEG) is achieved by the growth of a resistive AlGaN layer on top of the highly-resistive GaN layer. The transition from GaN to AlGaN creates polarization charge, which attracts mobile electrons that form the 2-DEG. To form a normally-off HEMT, which is required for power applications, the electrons under the gate are repelled by various techniques, such as embedding negative charge in the AlGaN layer and using gate materials with a suitable work function. Given that a positive gate voltage is needed to re-establish the 2-DEG under the gate, this region of the HEMT is analogous to a normally-off MOSFET. The gate-to-drain and gate-to-source regions are different from the MOSFET structure; however, they can be modeled as fixed drain and source resistances. mental static and dynamic characteristics, respectively. Section 6 will present the results of a double-pulse test (DPT) test to compare the effectiveness of the MOSFET LEVEL 3 model with the manufacturer model in LTSpice. Section 7 will conclude and foreshadow potential future work.

GaN Structure
The internal structure of a normally-off GaN-HEMT is shown in Figure 1 [22]. For HEMT fabrication on the silicon substrate, as is the case in the structure shown in Figure  1, aluminum nitride is used as a buffer layer to enable epitaxial growth of the active GaN layer. The formation of the two-dimensional electron gas (2-DEG) is achieved by the growth of a resistive AlGaN layer on top of the highly-resistive GaN layer. The transition from GaN to AlGaN creates polarization charge, which attracts mobile electrons that form the 2-DEG. To form a normally-off HEMT, which is required for power applications, the electrons under the gate are repelled by various techniques, such as embedding negative charge in the AlGaN layer and using gate materials with a suitable work function. Given that a positive gate voltage is needed to re-establish the 2-DEG under the gate, this region of the HEMT is analogous to a normally-off MOSFET. The gate-to-drain and gate-to-source regions are different from the MOSFET structure; however, they can be modeled as fixed drain and source resistances. The switching ability of power devices depend on its parasitic capacitances, which are gate-to-source capacitance (CGS), gate-to-drain capacitance (CGD), and drain-to-source capacitance (CDS). The three parasitic capacitances, shown in Figure 1, are analogous to the parasitic capacitances in the MOSFET structure. An important quantitative difference is that the parasitic capacitances in HEMTs are much smaller, enabling operation at much higher switching frequencies. A major difference between GaN-HEMT and MOSFET is the lack of an intrinsic body diode. However, the intrinsic diode in MOSFET models can be easily eliminated by adopting zero values of its key parameters. Therefore, standard MOSFET equations can be applied to GaN-HEMTs.

Selection of SPICE Equations
Modeling a GaN-HEMT requires mathematical equations for the drain current (ID) as a function of the applied gate-to-source (VGS) and drain-to-source (VDS) voltages. In SPICE, there are different sets of equations, referred to as model levels. LEVEL 1 is the The switching ability of power devices depend on its parasitic capacitances, which are gate-to-source capacitance (C GS ), gate-to-drain capacitance (C GD ), and drain-to-source capacitance (C DS ). The three parasitic capacitances, shown in Figure 1, are analogous to the parasitic capacitances in the MOSFET structure. An important quantitative difference is that the parasitic capacitances in HEMTs are much smaller, enabling operation at much higher switching frequencies. A major difference between GaN-HEMT and MOSFET is the lack of an intrinsic body diode. However, the intrinsic diode in MOSFET models can be easily eliminated by adopting zero values of its key parameters. Therefore, standard MOSFET equations can be applied to GaN-HEMTs.

Selection of SPICE Equations
Modeling a GaN-HEMT requires mathematical equations for the drain current (I D ) as a function of the applied gate-to-source (V GS ) and drain-to-source (V DS ) voltages. In SPICE, there are different sets of equations, referred to as model levels. LEVEL 1 is the simplest MOSFET model, which is based on the rudimentary device equations. Shah et al. [12] have used the equations of MOSFET LEVEL 1 to model GaN-HEMTs. However, their approach ignored the transfer characteristic for small V DS values, which in the case of any GaN-HEMT is impacted by the channel resistances outside the gate (these gate-to-source and gate-to-drain channel resistances do not exist in MOSFETs). In addition, the simulation accuracy of output characteristics is highly dependent on second-order effects, which are not included in [12]. The MOSFET LEVEL 2 model in SPICE includes second-order effects, but its equations are complex resulting in computational inefficiency and convergence issues. MOSFET LEVEL 3 is a much better choice, because its key equations are almost as simple as in MOSFET LEVEL 1 and as accurate as the complex equations in the MOSFET LEVEL 2 model. Therefore, the MOSFET LEVEL 3 model is utilized in this paper. In SPICE, a parameter labeled as LEVEL is assigned to MOSFET devices. By assigning LEVEL = 3, the SPICE LEVEL 3 model with the relevant equations and parameters is activated.
Equations (1) and (2), shown below, represent our selection from the rather complex MOSFET LEVEL 3 model, for the purpose of modeling GaN-HEMTs. The SPICE parameters that need to be specified to activate this selection are shown in the bold font. The list of these parameters is as follows: the transconductance parameter (KP), the threshold voltage (Vto), the channel length (L), the channel width (W), the drain resistance (Rd), the source resistance (Rs), the body-effect parameter (Gamma), the surface potential in strong inversion (Phi), and the mobility modulation constant (Theta). In this paper, we will also show that the subthreshold-current parameter (NFS) can be used to model the non-zero drain current for gate voltages below the threshold voltage. The default values of the remaining SPICE parameters ensure that all other equations are inactive [23].

Methods for Extraction of Initial Values of the Selected Parameters
In principle, the values of equation parameters can be obtained by nonlinear fitting to experimental data. However, the convergence of nonlinear fitting entirely depends on the specified initial values of the equation parameters. More importantly, an inadequate choice of the initial values can lead to poor fitting. Hence, this paper proposes a method for extraction of the initial values of selected key parameters as the first step, followed by nonlinear fitting to adjust the values of all parameters in Equations (1) and (2).
The channel length and width, L and W, are the geometric-design parameters specific to the power device under test. For commercial power-semiconductor devices these values are undisclosed by the manufactures and are very hard to determine. However, for the purpose of circuit simulation, it is sufficient to set the value of W and L to 1 µm, which makes the ratio W/L = 1. That leaves the transconductance parameter KP as the active parameter whose value will determine the value of the drain current in response to the applied voltages. In addition to KP, the initial values of Vto, Rs, and Rd are also needed to enable successful nonlinear fitting. Therefore, the proposed method for the extraction of initial parameter values includes KP, Vto, Rs, and Rd. The values of the remaining parameters can be adjusted by nonlinear fitting using their typical or default values.

Extraction of KP, and Vto
The values of KP and Vto can be determined from measured transfer characteristic of the corresponding GaN-HEMT, biased in the linear region by a small V DS , such as V DS = 100 mV. To illustrate the extraction method, we chose a commercial 600 V GaN power transistor, GS66504B, manufactured by GaN Systems. The transfer characteristic of GS66504B was measured with an Agilent Power Device Analyzer (B1505A) and is shown in Figure 2. This transfer characteristic was measured at V DS = 100 mV to ensure that the quadratic term in (1) is negligible. In the linear section of the measured transfer characteristic, marked by the dashed line in Figure 2, the values of Theta(V GS -Vto), RsI D , and RdI D are negligible, simplifying (1) to the following linear form: Electronics 2021, 10, 130

of 15
Equation (3) shows that the intersection between the extrapolated linear section of the measured transfer characteristic and the V GS axis that corresponds to Vto = 1.26 V. The slope of the linear section is 0.879 A/V, which corresponds to KP = 8.79 A/V 2 for the specific case of V DS = 100 mV and W = L = 1 µm.

Extraction of Rs, and Rd
As Figure 2 shows, the measured transfer characteristic deviates from the linear dependence (the dashed line) for higher V GS and I D values, mainly due to the impact of the channel resistances outside the gate. MOSFET LEVEL 3 model has embedded source and drain resistances, Rs and Rd, which we can use to model the impact of gate-to-source and gate-to-drain channel resistances in HEMTs as shown in Figure 3. Assuming that the dashed line in Figure 2 corresponds to the HEMT channel under the gate, we can model this section of the channel by the linear equation: (4), V GS0 can be expressed as: in Figure 2. This transfer characteristic was measured at VDS = 100 mV to ensure that the quadratic term in (1) is negligible. In the linear section of the measured transfer characteristic, marked by the dashed line in Figure 2, the values of Theta(VGS-Vto), RsID, and RdID are negligible, simplifying (1) to the following linear form: Equation (3) shows that the intersection between the extrapolated linear section of the measured transfer characteristic and the VGS axis that corresponds to Vto = 1.26 V. The slope of the linear section is 0.879 A/V, which corresponds to KP = 8.79 A/V 2 for the specific case of VDS = 100 mV and W = L = 1 µ m.

Extraction of Rs, and Rd
As Figure 2 shows, the measured transfer characteristic deviates from the linear dependence (the dashed line) for higher VGS and ID values, mainly due to the impact of the channel resistances outside the gate. MOSFET LEVEL 3 model has embedded source and drain resistances, Rs and Rd, which we can use to model the impact of gate-to-source and gate-to-drain channel resistances in HEMTs as shown in Figure 3. Assuming that the dashed line in Figure 2 corresponds to the HEMT channel under the gate, we can model this section of the channel by the linear equation: (4), VGS0 can be expressed as:  For smaller ID values, where the solid and dashed lines overlap, (Rs + Rd) ID is negligible and VGS0 = VGS. However, for larger ID values, (Rs + Rd) ID is not negligible and a larger VGS is required to achieve the same ID. According to Equation (3), the measured value of VGS is given by The difference ΔVG = VGS0 − VGS, which is illustrated in Figure 2, can be obtained from a larger V GS is required to achieve the same I D . According to Equation (3), the measured value of V GS is given by The difference ∆V G = V GS0 − V GS , which is illustrated in Figure 2, can be obtained from (5) and (6) as: Practically, a sufficiently high I D should be selected for the determination of Rs + Rd, because a sufficiently large ∆V G is needed to avoid a large error in the determined value for Rs + Rd. For the specific data shown in Figure 2, ∆V G = 0.8 V, V DS = 100 mV, I D = 0.84 A, and consequently Rs + Rd is obtained as 54.3 mΩ. For simplification of the model we assume equal drain and source resistances, which are Rs = Rd = 27.15 mΩ.

Typical Initial Values for Gamma, Phi, Theta, and NFS
Two parameters, Gamma and Phi, impact the drain current in the triode and saturation regions. The complete MOSFET LEVEL 3 equations in SPICE use Gamma and Phi to model second-order effects that we do not need in the case of HEMT modeling. As can be seen from the square brackets in Equations (1) and (2), it is sufficient to adjust the value of only one of these parameters. From semiconductor physics, we know that the value of Phi is set by the energy gap of the material. For the case of GaN, we can set Phi = 2 V, which leaves Gamma to be determined by the nonlinear fitting. For this purpose, the default value in SPICE, which is Gamma = 0, can be used as the initial value.
The default value for the parameter Theta in SPICE is also Theta = 0. This is sufficient as the initial value for the nonlinear fitting, which will adjust the final value of Theta to improve the fitting of both the transfer and the output characteristics. Finally, note that we can specify the parameter NFS to fit the subthreshold current on the transfer characteristic. The typical value of this parameter is 10 11 cm −2 . However, in the triode and saturation region, the NFS does not affect the I D as can be seen from Equations (1) and (2). Hence, the value of this parameter can be adjusted in SPICE simulation after the non-linear fitting is completed.

Non-Linear Fitting
Using the set values of W = L = 1 µm and Phi = 2 V, along with the initial values obtained in the previous sub-section (KP = 8.79 A/V 2 , Vto = 1.26 V, Rs = Rd = 27.15 mΩ) and the default values Theta = 0, and Gamma = 0, non-linear fitting is applied to Equations (1) and (2) to obtain the MOSFET LEVEL 3 parameters for GS66504B. The obtained fitted parameters are typed directly into SPICE and are verified by the experimental transfer and output characteristics as described in the next section. Note that, without NFS, the drain current calculated by SPICE will be zero for V GS < Vto. Hence, in order to fit the sub-threshold region for GS66504B, the value of NFS can be determined by running SPICE simulations with tried numbers around its typical value of 10 11 cm −2 .

Experimental Verification of the Static Characteristics
In addition to GS66504B GaN-HEMT, we have also selected a 650 V GaN transistor GS66506T to demonstrate the effectiveness of the selected MOSFET LEVEL 3 equations and the proposed parameter-extraction method. Putting the fitted values of the selected parameters directly into SPICE, the simulated transfer and output characteristics for both GS66504B and GS66506T are compared with the measured results in Figure 4. It is evident that SPICE simulations are in very good agreement with the measured results, which demonstrates that the selected MOSFET LEVEL 3 equations in SPICE can be used to model GaN-HEMTs. Precise modeling in the linear region of the output characteristics, shown in Figure 4b,d, is the most significant, as it corresponds to the operation of HEMTs as switches in on mode and thus, determines the on-state resistance.
GS66506T to demonstrate the effectiveness of the selected MOSFET LEVEL 3 equations and the proposed parameter-extraction method. Putting the fitted values of the selected parameters directly into SPICE, the simulated transfer and output characteristics for both GS66504B and GS66506T are compared with the measured results in Figure 4. It is evident that SPICE simulations are in very good agreement with the measured results, which demonstrates that the selected MOSFET LEVEL 3 equations in SPICE can be used to model GaN-HEMTs. Precise modeling in the linear region of the output characteristics, shown in Figure 4b,d, is the most significant, as it corresponds to the operation of HEMTs as switches in on mode and thus, determines the on-state resistance.

Modeling Dynamic Characteristics and Experimental Verification
The dynamic characteristics of GaN-HEMTs are determined by their inherent voltagedependent parasitic capacitances. For the designers of power-conversion circuits, the capacitance of utmost interest is the output capacitance, C OSS = C DS + C GD . Hence, modeling of C DS and C GD is crucial. The desired individual capacitance data from the datasheet is extracted using the grabit.m command in MATLAB and is shown in Figure 5.
The MOSFET LEVEL 3 model includes voltage-independent gate-to-source and gateto-drain overlap capacitances, through Cgso and Cgdo as the respective parameters. By observing the datasheets of many GaN-HEMTs, it can be perceived that the C GS and C GD are fairly constant for most of the voltage ranges; hence the voltage-independent parameters Cgso and Cgdo can be used to model these capacitances. Note that Cgso and Cgdo are defined as capacitances per unit of the channel width. Therefore, extracting the values of Cgso and Cgdo can be achieved by taking the values of C GS and C GD at the maximum V DS .

Modeling Dynamic Characteristics and Experimental Verification
The dynamic characteristics of GaN-HEMTs are determined by their inherent voltage-dependent parasitic capacitances. For the designers of power-conversion circuits, the capacitance of utmost interest is the output capacitance, COSS = CDS + CGD. Hence, modeling of CDS and CGD is crucial. The desired individual capacitance data from the datasheet is extracted using the grabit.m command in MATLAB and is shown in Figure 5. The MOSFET LEVEL 3 model includes voltage-independent gate-to-source and gate-to-drain overlap capacitances, through Cgso and Cgdo as the respective parameters. By observing the datasheets of many GaN-HEMTs, it can be perceived that the CGS and CGD are fairly constant for most of the voltage ranges; hence the voltage-independent parameters Cgso and Cgdo can be used to model these capacitances. Note that Cgso and Cgdo are defined as capacitances per unit of the channel width. Therefore, extracting the values of Cgso and Cgdo can be achieved by taking the values of CGS and CGD at the maximum VDS.
The MOSFET LEVEL 3 model also contains a diode model for the MOSFET P-N junction, which can be utilized to model the non-linear CDS for GaN-HEMTs. The diode model with relevant SPICE parameters is given by where, Cj is the zero-bias capacitance per unit area, PB is the built-in voltage, and Mj is the grading coefficient. To obtain the values of these SPICE parameters requires non-linear fitting. For this purpose, the CDS value from the datasheet at VDS = 0 V, should be divided by W × L and used as the initial value of Cj. The initial value of PB can be set at 2 V (similar to Phi) and the initial value of Mj can be set at 0.5 (typical value for one-sided abrupt P-N junctions and Schottky diodes). The modeled capacitances were compared to the extracted capacitance for both GaN-HEMTs in Figure 5. Good agreement is achieved at higher voltages, but there is a discrepancy between the simulated and datasheet CDS curves at lower voltages due to the internal physics of HEMTs.
The suitability of the capacitance modeling can be verified by comparing simulated and actual energy stored in the output capacitance (EOSS). This energy is stored while the GaN-HEMT is turning-off and the same energy is dissipated while the GaN-HEMT is turning-on. As high-frequency converters with GaN-HEMTs always require The MOSFET LEVEL 3 model also contains a diode model for the MOSFET P-N junction, which can be utilized to model the non-linear C DS for GaN-HEMTs. The diode model with relevant SPICE parameters is given by where, Cj is the zero-bias capacitance per unit area, PB is the built-in voltage, and Mj is the grading coefficient. To obtain the values of these SPICE parameters requires non-linear fitting. For this purpose, the C DS value from the datasheet at V DS = 0 V, should be divided by W × L and used as the initial value of Cj. The initial value of PB can be set at 2 V (similar to Phi) and the initial value of Mj can be set at 0.5 (typical value for one-sided abrupt P-N junctions and Schottky diodes). The modeled capacitances were compared to the extracted capacitance for both GaN-HEMTs in Figure 5. Good agreement is achieved at higher voltages, but there is a discrepancy between the simulated and datasheet C DS curves at lower voltages due to the internal physics of HEMTs.
The suitability of the capacitance modeling can be verified by comparing simulated and actual energy stored in the output capacitance (E OSS ). This energy is stored while the GaN-HEMT is turning-off and the same energy is dissipated while the GaN-HEMT is turning-on. As high-frequency converters with GaN-HEMTs always require soft-switching, the E OSS curve is crucial for achieving this condition. The E OSS curve can be calculated from the respective C OSS curve as: The E OSS curve calculated from Equation (10) for MOSFET LEVEL 3 modeling is compared to the datasheet and is shown in Figure 6. The results show that the proposed model is adequate for simple and fast simulation of HEMTs in power-electronic circuits.
The EOSS curve calculated from Equation (10) for MOSFET LEVEL 3 modeling is compared to the datasheet and is shown in Figure 6. The results show that the proposed model is adequate for simple and fast simulation of HEMTs in power-electronic circuits.

Model Validation by LTSpice Simulation
For power applications, a good way to validate the dynamic characteristics is to check how the developed model behaves under switching conditions. The manufacturer of GS66506T, GaN Systems, provides their own developed LTSpice model to be used for simulation purposes (the model is shown in Appendix A). Moreover, GaN Systems has developed a half-bridge double-pulse test (DPT) in LTSpice to evaluate the switching performance of their developed models in comparison to experimental measurements. As shown in the application notes [24], the manufacturer model is in good agreement with the experimental switching conditions. Hence, an identical double-pulse test platform developed by GaN Systems is used in this paper to verify the switching characteristics of the developed MOSFET LEVEL 3 model. Since LTSpice shares the same core libraries as any other SPICE version, the MOSFET LEVEL 3 model can be used in LTSpice simulation with adequately set parameter values, as shown in Appendix B. It is worth noting that two more parameters, namely, gate resistance (RG) and drain-to-source shunt resistance (RDS), are added while converting SPICE to LTSpice models to avoid convergence problem in LTSpice. The addition of these parameters does not have much effect on the electrical characteristics of the modeled GaN-HEMT. The value of RG is directly taken from the GS66506T datasheet as 1.1 Ω. The RDS is theoretically infinite; hence a very high value of 1 MΩ is used in the model. The DPT test bench in LTSpice is shown in Figure 7. The parameters of the DPT circuit along with their values are given in Table 1. The switching waveforms of the DPT circuit with both the manufacturer's and MOSFET LEVEL 3 models as the device under test (DUT) are shown in Figure 8.

Model Validation by LTSpice Simulation
For power applications, a good way to validate the dynamic characteristics is to check how the developed model behaves under switching conditions. The manufacturer of GS66506T, GaN Systems, provides their own developed LTSpice model to be used for simulation purposes (the model is shown in Appendix A). Moreover, GaN Systems has developed a half-bridge double-pulse test (DPT) in LTSpice to evaluate the switching performance of their developed models in comparison to experimental measurements. As shown in the application notes [24], the manufacturer model is in good agreement with the experimental switching conditions. Hence, an identical double-pulse test platform developed by GaN Systems is used in this paper to verify the switching characteristics of the developed MOSFET LEVEL 3 model. Since LTSpice shares the same core libraries as any other SPICE version, the MOSFET LEVEL 3 model can be used in LTSpice simulation with adequately set parameter values, as shown in Appendix B. It is worth noting that two more parameters, namely, gate resistance (RG) and drain-to-source shunt resistance (RDS), are added while converting SPICE to LTSpice models to avoid convergence problem in LTSpice. The addition of these parameters does not have much effect on the electrical characteristics of the modeled GaN-HEMT. The value of RG is directly taken from the GS66506T datasheet as 1.1 Ω. The RDS is theoretically infinite; hence a very high value of 1 MΩ is used in the model. The DPT test bench in LTSpice is shown in Figure 7. The parameters of the DPT circuit along with their values are given in Table 1 Figure 8a,b shows the results with an operating DC bus voltage of 400 V and Figure 8c,d shows the results with an operating switching current of 22.5 A. Despite the simplicity of MOSFET LEVEL 3 model, the turn-off transients shows satisfactory agreement with the manufacturer model, whereas a bigger difference can be seen in Figure 8a,c during turn-on transients for both operating points. This deviation during the turn-on interval is most likely caused by the mismatch of modeled C DS and C GD at lower voltages. In addition, the MOSFET LEVEL 3 model does not include any parasitic inductances, which results in more pronounced oscillations as seen in Figure 8. However, for power applications, the power distributions during switching intervals are more relevant. The turn-on and turn-off power losses using the MOSFET LEVEL 3 model and manufacturer model for the two operating points are shown in Figure 9. Furthermore, the total energy losses during switching transients are presented in Table 2.    Figure 8a,c during turn-on transients for both operating points. This deviation during the turn-on interval is most likely caused by the mismatch of modeled CDS and CGD at lower voltages. In addition, the MOSFET LEVEL 3 model does not include any parasitic inductances, which results in more pronounced oscillations as seen in Figure 8. However, for power applications, the power distributions during switching intervals are more relevant. The turn-on and turn-off power losses using the MOSFET LEVEL 3 model and manufacturer model for the two operating points are shown in Figure 9. Furthermore, the total energy losses during switching transients are presented in Table 2.    Figure 9 and Table 2 demonstrate the effectiveness of the simple MOSFET LEVEL 3 model approach in comparison to the complex manufacturer model. Despite the observed deviation in the switching waveforms shown in Figure 8, the switching energy losses for the full current rating calculated from both models are quite close. In order to further evaluate the switching performance of the MOSFET LEVEL 3 model, the drain current was varied keeping the dc bus voltage constant and the comparison to the manufacturer model is shown in Figure 10. As can be seen, the simple MOSFET LEVEL 3 model provides adequate results in comparison to the complex manufacturer model. The observed variation between both curves in Figure 10 is acceptable in terms of designing power circuits, whereas the availability, transparency, and simplicity of the MOSFET LEVEL 3 model enables fast and effective simulations during the design of practical power-switching circuits.   Figure 9 and Table 2 demonstrate the effectiveness of the simple MOSFET LEVEL 3 model approach in comparison to the complex manufacturer model. Despite the observed deviation in the switching waveforms shown in Figure 8, the switching energy losses for the full current rating calculated from both models are quite close. In order to further evaluate the switching performance of the MOSFET LEVEL 3 model, the drain current was varied keeping the dc bus voltage constant and the comparison to the manufacturer model is shown in Figure 10. As can be seen, the simple MOSFET LEVEL 3 model provides adequate results in comparison to the complex manufacturer model. The observed variation between both curves in Figure 10 is acceptable in terms of designing power circuits, whereas the availability, transparency, and simplicity of the MOSFET LEVEL 3 model enables fast and effective simulations during the design of practical power-switching circuits.

Conclusions
This paper demonstrates that standard MOSFET equations, available in the SPICE circuit simulator, can be used to model commercial GaN-HEMTs. A set of equations from the available MOSFET LEVEL 3 equations in SPICE has been selected and a method for the extraction of the related SPICE parameters has been proposed. The validity of this approach, including the proposed parameter-extraction techniques, is verified by demonstrating an adequate match between simulated and measured static and dynamic characteristics. Furthermore, a DPT is performed for one of the GaN-HEMTs and is compared to its manufacturer model. It is demonstrated that the simple MOSFET LEVEL

Conclusions
This paper demonstrates that standard MOSFET equations, available in the SPICE circuit simulator, can be used to model commercial GaN-HEMTs. A set of equations from the available MOSFET LEVEL 3 equations in SPICE has been selected and a method for the extraction of the related SPICE parameters has been proposed. The validity of this approach, including the proposed parameter-extraction techniques, is verified by demonstrating an adequate match between simulated and measured static and dynamic characteristics. Furthermore, a DPT is performed for one of the GaN-HEMTs and is compared to its manufacturer model. It is demonstrated that the simple MOSFET LEVEL 3 model provides adequate results in terms of total switching energy losses. The transparency of the MOSFET LEVEL 3 model is further helpful in adapting to any GaN-HEMTs which is quite difficult in the manufacturer model. Due to the simplicity of this approach and the familiarity of circuit designs with the standard MOSFET SPICE model, this approach is ideally suited for the analysis and design of circuits with GaN-HEMTs. Future work may include an improvement of modeling of voltage-dependent capacitances to obtain switching energy losses as close as possible to the complex manufacturer model while keeping the model simple and transparent.