# A Practical Method of Parameter Extraction for the VBIC Model

## Transcript Of A Practical Method of Parameter Extraction for the VBIC Model

A Practical Method of Parameter Extraction for the VBIC Model used on a GaAs HBT.

Morten Olavsbråten

Norwegian Institute of Science and Technology (NTNU) Dept. of telecommunication, N-7491 Trondheim, Norway

[email protected]

ABSTRACT

This paper shows a simple practical method of extracting the major parameters in the VBIC model used on a GaAs HBT. The method was developed mainly for the circuit designer, who want to have an easy practical way of extracting good model parameters from a few simple measurements, and due to the lack of such simple methods for the VBIC model used on GaAs HBT. The extraction method is based on a semi-analytic approach. The method includes a simple analytic extraction of the non-ideal current sources Iben, Nen, Ibcn and Ncn. The algorithm for the parameter extraction is presented. To verify the method, measurements from Caswell Technology were used. The extracted model shows good agreement with measurements.

I. INTRODUCTION

In the last few years, many papers describing new models for the Hetero-junction Bipolar Transistor (HBT) have been published. The main drawback is that the models are not implemented in commercial simulator software like Agilent ADS. One exception is the VBIC model by McAndrew et al (1). The VBIC model is discussed by Najm (2). In order to use this model in a practical circuit, one must extract the values for many of the numerous parameters of the model. This paper describes a practical method of using a few simple measurements to extract the major parameters. The extraction of the self-heating parameter Rth is based upon Marsh (3). Cao et al (4) shows a method of parameter extraction for silicon BJT.

II. PARAMETER EXTRACTION

The VBIC model is shown in Figure 1. The symbols used for the model parameters are all taken from the model implementation in Agilent ADS (there are 86 parameters that can be extracted). The model includes effects like selfheating, non-ideal current-sources, and distributed base- and collector modeling effects. This model also includes a parasitic substrate transistor, which may be important in silicon technology, but can be neglected in GaAs design due to low substrate losses and high isolation. This parasitic will not be discussed any further, and its parameters Rs, Isp, Ibeip, Ibenp, Ibcip, Ibcnp, Ikp, Xrs, Cjep, Cjcp, Wsp, Nfp, Ncip and Ncnp are all set to the default values.

A. FORWARD GUMMEL MEASUREMENT Low current region: From this region we can extract the saturation currents and the emission coefficients. The base current Ib can be written (including the non-ideal currents) as shown in equation 1:

qVbe

qVbe

Ib = Ibei ⋅ (e NeikT −1) + Iben ⋅ (e NenkT −1)

(1)

This gives two asymptotic lines (line C, given by Vb approaching zero and line B where the slope is maximum) in a logarithmic plot as illustrated in Figure 2. We get the ideal base-emitter saturation current Ibei where the asymptotic line (B) crosses the y-axes (Vb=0). The ideal base-emitter emission coefficient Nei is then given by equation 2:

N ei

=

q ln(10)⋅kT ⋅SlopeB

(2)

We find the non-ideal saturation current Iben and the emission coefficient Nen of the base-emitter the same way using the asymptotic line (C). The transport saturation current Is and the forward emission coefficient Nf, are found the same way using the asymptotic line (A) fitted to the measurement of Ic. If higher accuracy is needed, these six parameters (Ibei, Nei, Iben, Nen, Is, Nf) should be optimized to fit the Ib and Ic measurements. High current region: In this region many parameters have influence on the behavior of the currents. Therefore an

iterative approach is the most convenient approach. The most dominant parameter here is the emitter resistance Re. This should be optimized to the measurements of Ib and Ic.

B. REVERSE GUMMEL MEASUREMENT Low current region: When applying the same method as for the forward Gummel plot, we find from the measured

base-collector current Ibc the ideal and non-ideal saturation currents Ibci and Ibcn, and emission coefficients Nci and Ncn of the base-collector. Again, if higher accuracy is needed, these four parameters should be optimized to fit the Ibc measurements. In the VBIC model it is recommended to set the reverse emission coefficient Nr = Nf to ensure convergence.

High current region: In this region many parameters have influence on the behavior of the currents, equivalent to the forward Gummel measurements. The most dominant parameter here is the total resistance in the base-collector junction (Rb + Rc). This should be fitted by optimization, to the measurements. The reverse knee current Ikr could have been optimized in this region, but in HBT’s the reverse knee current is usually very high and can be set to infinity (in ADS Ikr = 0).

C. DC MEASUREMENT The DC-measurements include Ic versus Vce curves for different base currents Ib, and if possible at different temperatures. Self-heating: If measurements of the Ic-Vce curves are available at different temperatures, the thermal resistance Rth can be calculated from an excellent method given by Marsh (3), or from another method given by Bovolon et al (5). If no temperature measurements are available, the thermal resistance Rth must be estimated from material properties and the geometry of the device. The thermal capacitance Cth can be calculated from the material properties and the geometry of the device, or by pulsed measurements. Temperature dependent parameters: The three most important parameters are: Xis, the temperature exponent of Is, Xii, the temperature exponent of Ibei, Ibci, Ibeip, Ibcip and Tnf, the temperature coefficient of Nf. These parameters are found by optimization of the IV curves at two or more temperatures. Increased accuracy can be obtained by including measurements of the Vbe versus Vce curves. If no temperature measurements are available, the best results are obtained by optimizing Xii and leaving Xis and Tnf at their default value. The temperature exponent Xin of the non-ideal saturation currents Iben, Ibcn, Ibenp and Ibcnp is set to simulator default values. The effect of this parameter is very small. Collector resistance: The total collector resistance Rc (Rci + Rcx) can be found directly from the IV-curves. (Rc+Re)-1 is the maximal slope in the triode region of the IV-curves as shown by Lu et al (6). The emitter resistance Re was found in the high current region of the forward Gummel plots. The temperature exponent of the emitter, base and collector resistances Xre, Xrb, and Xrc can be set to zero when operating in a limited (normal) temperature range. When this is not the case and when temperature measurements are available the Xrc can be calculated by finding Rc, as described above, at different temperatures. Then, for all practical purposes, Xre = Xrb = Xrc is set. Forward knee currents: The forward knee current Ikf can be optimized in this region to ensure the best fit to the Ic-Vce curves. However, in HBT’s the forward knee current is usually very high and can often be set to infinity (in ADS Ikf = 0).

D. COLD-CAPACITOR MEASUREMENT The cold-capacitor measurements are s-parameter measurements at different base voltages, i.e. Vb = -2.5 → 1.1V, with Vce=0. Optimization in ADS is used to extract the parameters. Bias independent parameters: The s-parameter measurement at zero bias (Vb = 0) is used to find almost all the passive, bias independent parameters. These parameters are the zero-bias base-emitter and base-collector capacitances Cje and Cjc and the extrinsic base-emitter and base-collector capacitances Cbeo and Cbco. (If extrinsic inductors (Lb, Lc and Le) are included at the base, collector and emitter, these parameters are also found here). Bias dependent parameters: The bias dependence of base-emitter and base-collector capacitances Cje and Cjc are found by optimizing the measured s-parameters to the best fit as follows: 1) At negative bias (i.e. Vb = -2.5V), optimize the base-emitter/base-collector junction exponents Me and Mc. 2) At moderate positive bias (i.e. Vb = 0.5V), optimize the base-emitter/base-collector grading coefficient Pe and Pc. 3) At positive bias (i.e. Vb = 0.9V or higher), optimize the forward bias junction capacitance threshold Fc.

At this step it is also possible to find the distributed effect of the base, which is given by Wbe (the portion of Ibei that runs through the distributed current source Ibex), and the distribution of the base resistance Rb between the intrinsic and extrinsic base resistance Rbi and Rbx. If the fit to the measured s-parameters in this step are satisfactory without these parameters, the distributed effect of the base can be neglected i.e. Wbe = 1, Rbx = 0. If necessary, repeat these three steps for the best results.

E. ACTIVE S-PARAMETER MEASUREMENT The s-parameter measurement for an active forward bias is used to find the forward transit time Tf, and to find the distribution of Rc between Rci and Rcx, by optimization. This measurement can also be used to fine-tune other parameters like the resistors found in the DC part of the extraction. If s-parameters for several bias points are measured, the bias dependent parameters Qtf, Xtf, Vtf, and Itf, of the forward transit time Tf can be found. This has not been analyzed here.

F. OTHER PARAMETERS The activation energy parameters Ea, Eaie, Eaic, Eais, Eane, Eanc, and Eans are all set to 1.42 eV for GaAs. Noise parameters are not included in this method, and should be set to the simulator defaults. The forward and reverse Early voltages Vef and Ver are set to infinity (0 in ADS), due to the fact that the GaAs HBT is nearly ideal in that respect. All other parameters should be set to the default values.

III. THE ALGORITHM

In this algorithm, when a parameter has been extracted, its value is used in all the consequent steps.

1. From the low current region of the forward Gummel plot, Ibei, Nei, Iben, Nen, Is, and Nf are extracted. 2. From the low current region of the reverse Gummel plot, Ibci, Nci, Ibcn, and Ncn are extracted. 3. From the self-heating part of the DC-measurement, Rth and Cth are extracted. 4. From the high current region of the forward Gummel plot, Re is extracted. 5. From the high current region of the reverse Gummel plot, the Rb + Rc (if needed Ikr) are extracted. 6. From the collector resistance part of the DC-measurement, Rc and if necessary Xrc (= Xrb = Xre) are extracted. 7. From the Temperature dependent parameters part of the DC-measurements, Xis, Xii and Tnf are extracted. 8. Steps 4 to 7 are repeated until the best fits are obtained. 9. From the Bias independent parameters of the cold-capacitor measurements, Cje, Cjc, Cbeo and Cbco (and if included

Lb, Lc, Le) are extracted. 10. From the Bias dependent parameters of the cold-capacitor measurements, Me, Mc, Pe, Pc and Fc (and if needed Wbe

and the distribution of Rb between Rbi and Rbx) are extracted. 11. From the active s-parameter measurement, Tf and the distribution of Rc between Rci and Rcx are extracted. If

necessary this s-parameter measurement can be used to fine-tune the resistances Re, Rb and Rc. 12. If parameters for distributed base are used in Step 10, or some of the resistors have been fine-tuned in Step 11, then

Step 4 to 11 should be repeated for the best possible fit. Steps 4 to 6 do not extract new values for the resistors that have been fine-tuned in Step 11.

IV. RESULTS

The algorithm described above is verified by applying it to measurements provided by Caswell Technology, Marconi Caswell Limited. The measurements referred to in Figure 3, 4 and 5 are all statistical mean values of about 40 single finger transistors. Figure 3 and 4 show the simulated forward and reverse Gummel curves respectively. Solid lines are simulated results and dots are measurements. Figure 5 shows an example of the simulated and measured results for all four s-parameters at a bias point, Vce = 5V and Ic = 20mA. Similar results are seen at different bias points. Overall, we see that the model predicts the measured results very well.

V. CONCLUSION

In this paper a practical method of parameter extraction for the VBIC model applied on a GaAs HBT is proposed. The method extracts all major parameters for the model, including self-heating and the non-ideal parameters. Only a few simple measurements are necessary to extract the parameters (forward and reverse Gummel, DC-curves, cold capacitor s-parameters and s-parameters for a ‘normal’ bias are needed). Simulations with the extracted parameters for the VBIC model show good agreement with the measured results.

ACKNOWLISHMENT

The author would like to acknowledge Caswell Technology for supplying the measurements, and a special thanks to Dr. Steve Marsh at Caswell Technology for valuable discussions on the HBT technology.

REFERENCS

(1) Colin C. McAndrew et al.,”VBIC95, The Vertical Bipolar Inter-Company Model”, 1996, IEEE Journal of SolidState Circuits, Vol. 13, No. 10, pp. 1476-1482.

(2) Farid Najm, “VBIC95: An Improved Bipolar Transistor Model”, IEEE Circuits & Devices, March 1996, pp 11-15. (3) Steve P. Marsh, ”Direct Extraction Technique to Derive the Junction Temperature of HBT’s Under High Self-

Heating Bias Conditions”, 2000, IEEE Trans. on Electron Devices, Vol. 47, No. 2, pp. 288-291. (4) X. Cao et al., “Parameter Extraction and Optimization for New Industry Standard VBIC Model”, 1998, ADAMS

’98 2nd Int. Conf. On Adv. Semicond. Devices and Microsyst., pp. 107-115. (5) N. Bovolon et al., “A Simple Method for the Thermal Resistance Measurement of AlGaAs/GaAs Hetero-junction

Bipolar Transistor”, 1998, IEEE Trans. on Electron Devices, Vol. 45, No. 8, pp 1846-1847. (6) Ke Lu et al., “A New Large-Signal AlGaAs/GaAs HBT Model Including Self-Heating Effects, with Corresponding

Parameter-Extraction Procedure”, 1995, IEEE Trans. on MTT, Vol. 43, No. 7, pp 1433-1445.

Figure 1: The VBIC model

Forward Gummel plot

104

Ic

Ic and Ib [µA]

Ib

102

Slope A

100

Slope B

Slope C

10-2 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5

Vc and Vb [V] Figure 2: Forward Gummel plot and asymptotic lines

Figure 3: Forward Gummel plot: Simulated and measured. Figure 4: Reverse Gummel plot: Simulated and measured.

Figure 5: S-parameters at Vce = 5V, Ic = 20mA: Simulated (solid lines) and measured (magnitude: circles, phase: triangles).

Morten Olavsbråten

Norwegian Institute of Science and Technology (NTNU) Dept. of telecommunication, N-7491 Trondheim, Norway

[email protected]

ABSTRACT

This paper shows a simple practical method of extracting the major parameters in the VBIC model used on a GaAs HBT. The method was developed mainly for the circuit designer, who want to have an easy practical way of extracting good model parameters from a few simple measurements, and due to the lack of such simple methods for the VBIC model used on GaAs HBT. The extraction method is based on a semi-analytic approach. The method includes a simple analytic extraction of the non-ideal current sources Iben, Nen, Ibcn and Ncn. The algorithm for the parameter extraction is presented. To verify the method, measurements from Caswell Technology were used. The extracted model shows good agreement with measurements.

I. INTRODUCTION

In the last few years, many papers describing new models for the Hetero-junction Bipolar Transistor (HBT) have been published. The main drawback is that the models are not implemented in commercial simulator software like Agilent ADS. One exception is the VBIC model by McAndrew et al (1). The VBIC model is discussed by Najm (2). In order to use this model in a practical circuit, one must extract the values for many of the numerous parameters of the model. This paper describes a practical method of using a few simple measurements to extract the major parameters. The extraction of the self-heating parameter Rth is based upon Marsh (3). Cao et al (4) shows a method of parameter extraction for silicon BJT.

II. PARAMETER EXTRACTION

The VBIC model is shown in Figure 1. The symbols used for the model parameters are all taken from the model implementation in Agilent ADS (there are 86 parameters that can be extracted). The model includes effects like selfheating, non-ideal current-sources, and distributed base- and collector modeling effects. This model also includes a parasitic substrate transistor, which may be important in silicon technology, but can be neglected in GaAs design due to low substrate losses and high isolation. This parasitic will not be discussed any further, and its parameters Rs, Isp, Ibeip, Ibenp, Ibcip, Ibcnp, Ikp, Xrs, Cjep, Cjcp, Wsp, Nfp, Ncip and Ncnp are all set to the default values.

A. FORWARD GUMMEL MEASUREMENT Low current region: From this region we can extract the saturation currents and the emission coefficients. The base current Ib can be written (including the non-ideal currents) as shown in equation 1:

qVbe

qVbe

Ib = Ibei ⋅ (e NeikT −1) + Iben ⋅ (e NenkT −1)

(1)

This gives two asymptotic lines (line C, given by Vb approaching zero and line B where the slope is maximum) in a logarithmic plot as illustrated in Figure 2. We get the ideal base-emitter saturation current Ibei where the asymptotic line (B) crosses the y-axes (Vb=0). The ideal base-emitter emission coefficient Nei is then given by equation 2:

N ei

=

q ln(10)⋅kT ⋅SlopeB

(2)

We find the non-ideal saturation current Iben and the emission coefficient Nen of the base-emitter the same way using the asymptotic line (C). The transport saturation current Is and the forward emission coefficient Nf, are found the same way using the asymptotic line (A) fitted to the measurement of Ic. If higher accuracy is needed, these six parameters (Ibei, Nei, Iben, Nen, Is, Nf) should be optimized to fit the Ib and Ic measurements. High current region: In this region many parameters have influence on the behavior of the currents. Therefore an

iterative approach is the most convenient approach. The most dominant parameter here is the emitter resistance Re. This should be optimized to the measurements of Ib and Ic.

B. REVERSE GUMMEL MEASUREMENT Low current region: When applying the same method as for the forward Gummel plot, we find from the measured

base-collector current Ibc the ideal and non-ideal saturation currents Ibci and Ibcn, and emission coefficients Nci and Ncn of the base-collector. Again, if higher accuracy is needed, these four parameters should be optimized to fit the Ibc measurements. In the VBIC model it is recommended to set the reverse emission coefficient Nr = Nf to ensure convergence.

High current region: In this region many parameters have influence on the behavior of the currents, equivalent to the forward Gummel measurements. The most dominant parameter here is the total resistance in the base-collector junction (Rb + Rc). This should be fitted by optimization, to the measurements. The reverse knee current Ikr could have been optimized in this region, but in HBT’s the reverse knee current is usually very high and can be set to infinity (in ADS Ikr = 0).

C. DC MEASUREMENT The DC-measurements include Ic versus Vce curves for different base currents Ib, and if possible at different temperatures. Self-heating: If measurements of the Ic-Vce curves are available at different temperatures, the thermal resistance Rth can be calculated from an excellent method given by Marsh (3), or from another method given by Bovolon et al (5). If no temperature measurements are available, the thermal resistance Rth must be estimated from material properties and the geometry of the device. The thermal capacitance Cth can be calculated from the material properties and the geometry of the device, or by pulsed measurements. Temperature dependent parameters: The three most important parameters are: Xis, the temperature exponent of Is, Xii, the temperature exponent of Ibei, Ibci, Ibeip, Ibcip and Tnf, the temperature coefficient of Nf. These parameters are found by optimization of the IV curves at two or more temperatures. Increased accuracy can be obtained by including measurements of the Vbe versus Vce curves. If no temperature measurements are available, the best results are obtained by optimizing Xii and leaving Xis and Tnf at their default value. The temperature exponent Xin of the non-ideal saturation currents Iben, Ibcn, Ibenp and Ibcnp is set to simulator default values. The effect of this parameter is very small. Collector resistance: The total collector resistance Rc (Rci + Rcx) can be found directly from the IV-curves. (Rc+Re)-1 is the maximal slope in the triode region of the IV-curves as shown by Lu et al (6). The emitter resistance Re was found in the high current region of the forward Gummel plots. The temperature exponent of the emitter, base and collector resistances Xre, Xrb, and Xrc can be set to zero when operating in a limited (normal) temperature range. When this is not the case and when temperature measurements are available the Xrc can be calculated by finding Rc, as described above, at different temperatures. Then, for all practical purposes, Xre = Xrb = Xrc is set. Forward knee currents: The forward knee current Ikf can be optimized in this region to ensure the best fit to the Ic-Vce curves. However, in HBT’s the forward knee current is usually very high and can often be set to infinity (in ADS Ikf = 0).

D. COLD-CAPACITOR MEASUREMENT The cold-capacitor measurements are s-parameter measurements at different base voltages, i.e. Vb = -2.5 → 1.1V, with Vce=0. Optimization in ADS is used to extract the parameters. Bias independent parameters: The s-parameter measurement at zero bias (Vb = 0) is used to find almost all the passive, bias independent parameters. These parameters are the zero-bias base-emitter and base-collector capacitances Cje and Cjc and the extrinsic base-emitter and base-collector capacitances Cbeo and Cbco. (If extrinsic inductors (Lb, Lc and Le) are included at the base, collector and emitter, these parameters are also found here). Bias dependent parameters: The bias dependence of base-emitter and base-collector capacitances Cje and Cjc are found by optimizing the measured s-parameters to the best fit as follows: 1) At negative bias (i.e. Vb = -2.5V), optimize the base-emitter/base-collector junction exponents Me and Mc. 2) At moderate positive bias (i.e. Vb = 0.5V), optimize the base-emitter/base-collector grading coefficient Pe and Pc. 3) At positive bias (i.e. Vb = 0.9V or higher), optimize the forward bias junction capacitance threshold Fc.

At this step it is also possible to find the distributed effect of the base, which is given by Wbe (the portion of Ibei that runs through the distributed current source Ibex), and the distribution of the base resistance Rb between the intrinsic and extrinsic base resistance Rbi and Rbx. If the fit to the measured s-parameters in this step are satisfactory without these parameters, the distributed effect of the base can be neglected i.e. Wbe = 1, Rbx = 0. If necessary, repeat these three steps for the best results.

E. ACTIVE S-PARAMETER MEASUREMENT The s-parameter measurement for an active forward bias is used to find the forward transit time Tf, and to find the distribution of Rc between Rci and Rcx, by optimization. This measurement can also be used to fine-tune other parameters like the resistors found in the DC part of the extraction. If s-parameters for several bias points are measured, the bias dependent parameters Qtf, Xtf, Vtf, and Itf, of the forward transit time Tf can be found. This has not been analyzed here.

F. OTHER PARAMETERS The activation energy parameters Ea, Eaie, Eaic, Eais, Eane, Eanc, and Eans are all set to 1.42 eV for GaAs. Noise parameters are not included in this method, and should be set to the simulator defaults. The forward and reverse Early voltages Vef and Ver are set to infinity (0 in ADS), due to the fact that the GaAs HBT is nearly ideal in that respect. All other parameters should be set to the default values.

III. THE ALGORITHM

In this algorithm, when a parameter has been extracted, its value is used in all the consequent steps.

1. From the low current region of the forward Gummel plot, Ibei, Nei, Iben, Nen, Is, and Nf are extracted. 2. From the low current region of the reverse Gummel plot, Ibci, Nci, Ibcn, and Ncn are extracted. 3. From the self-heating part of the DC-measurement, Rth and Cth are extracted. 4. From the high current region of the forward Gummel plot, Re is extracted. 5. From the high current region of the reverse Gummel plot, the Rb + Rc (if needed Ikr) are extracted. 6. From the collector resistance part of the DC-measurement, Rc and if necessary Xrc (= Xrb = Xre) are extracted. 7. From the Temperature dependent parameters part of the DC-measurements, Xis, Xii and Tnf are extracted. 8. Steps 4 to 7 are repeated until the best fits are obtained. 9. From the Bias independent parameters of the cold-capacitor measurements, Cje, Cjc, Cbeo and Cbco (and if included

Lb, Lc, Le) are extracted. 10. From the Bias dependent parameters of the cold-capacitor measurements, Me, Mc, Pe, Pc and Fc (and if needed Wbe

and the distribution of Rb between Rbi and Rbx) are extracted. 11. From the active s-parameter measurement, Tf and the distribution of Rc between Rci and Rcx are extracted. If

necessary this s-parameter measurement can be used to fine-tune the resistances Re, Rb and Rc. 12. If parameters for distributed base are used in Step 10, or some of the resistors have been fine-tuned in Step 11, then

Step 4 to 11 should be repeated for the best possible fit. Steps 4 to 6 do not extract new values for the resistors that have been fine-tuned in Step 11.

IV. RESULTS

The algorithm described above is verified by applying it to measurements provided by Caswell Technology, Marconi Caswell Limited. The measurements referred to in Figure 3, 4 and 5 are all statistical mean values of about 40 single finger transistors. Figure 3 and 4 show the simulated forward and reverse Gummel curves respectively. Solid lines are simulated results and dots are measurements. Figure 5 shows an example of the simulated and measured results for all four s-parameters at a bias point, Vce = 5V and Ic = 20mA. Similar results are seen at different bias points. Overall, we see that the model predicts the measured results very well.

V. CONCLUSION

In this paper a practical method of parameter extraction for the VBIC model applied on a GaAs HBT is proposed. The method extracts all major parameters for the model, including self-heating and the non-ideal parameters. Only a few simple measurements are necessary to extract the parameters (forward and reverse Gummel, DC-curves, cold capacitor s-parameters and s-parameters for a ‘normal’ bias are needed). Simulations with the extracted parameters for the VBIC model show good agreement with the measured results.

ACKNOWLISHMENT

The author would like to acknowledge Caswell Technology for supplying the measurements, and a special thanks to Dr. Steve Marsh at Caswell Technology for valuable discussions on the HBT technology.

REFERENCS

(1) Colin C. McAndrew et al.,”VBIC95, The Vertical Bipolar Inter-Company Model”, 1996, IEEE Journal of SolidState Circuits, Vol. 13, No. 10, pp. 1476-1482.

(2) Farid Najm, “VBIC95: An Improved Bipolar Transistor Model”, IEEE Circuits & Devices, March 1996, pp 11-15. (3) Steve P. Marsh, ”Direct Extraction Technique to Derive the Junction Temperature of HBT’s Under High Self-

Heating Bias Conditions”, 2000, IEEE Trans. on Electron Devices, Vol. 47, No. 2, pp. 288-291. (4) X. Cao et al., “Parameter Extraction and Optimization for New Industry Standard VBIC Model”, 1998, ADAMS

’98 2nd Int. Conf. On Adv. Semicond. Devices and Microsyst., pp. 107-115. (5) N. Bovolon et al., “A Simple Method for the Thermal Resistance Measurement of AlGaAs/GaAs Hetero-junction

Bipolar Transistor”, 1998, IEEE Trans. on Electron Devices, Vol. 45, No. 8, pp 1846-1847. (6) Ke Lu et al., “A New Large-Signal AlGaAs/GaAs HBT Model Including Self-Heating Effects, with Corresponding

Parameter-Extraction Procedure”, 1995, IEEE Trans. on MTT, Vol. 43, No. 7, pp 1433-1445.

Figure 1: The VBIC model

Forward Gummel plot

104

Ic

Ic and Ib [µA]

Ib

102

Slope A

100

Slope B

Slope C

10-2 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5

Vc and Vb [V] Figure 2: Forward Gummel plot and asymptotic lines

Figure 3: Forward Gummel plot: Simulated and measured. Figure 4: Reverse Gummel plot: Simulated and measured.

Figure 5: S-parameters at Vce = 5V, Ic = 20mA: Simulated (solid lines) and measured (magnitude: circles, phase: triangles).