Mobile QR Code QR CODE

  1. (Electrical Engineering College, Henan University of Science and Technology, Luoyang 471023, China)



Heterojunction bipolar transistor, small-signal equivalent circuit, peeling algorithm, AC current crowding, distributed base-collector capacitance

I. INTRODUCTION

Recently, heterojunction bipolar transistors (HBTs) have been widely applied to design high-speed analog, digital and mixed-signal integrated circuits because of their high-speed capability and high-power applications (1-6). Consequently, an accurate small-signal model for HBTs is very important tool used to characterize transistors and improve fabrication process. There are two usually-used extraction methods for the determination of elements parameters in small-signal models, which are namely optimization method and direction method. However, direction method is more accurate than optimization method, since the utility of the optimization method is severely dependent on the initial value of elements. In the direct extraction methods, one often adopted π (7-9) or T (10-13) topologies to represent the small-signal equivalent circuits. Because the T topology is developed by the physical mechanism of BJT devices, it has been widely adopted over the last decade. However, a π topology is reduced rather than a T topology while the common large-signal model (e.g., VBIC model, Gummel-Poon model) is linearized (14).

In our previously published article (15), a direct parameter-extraction technique for InP HBT small-signal model was presented, which is based on peeling algorithm. In the parameter-extraction technique, a π-type small-signal model was adopted. However, the small-signal model only considered the principle characteristics of HBTs, ignoring base-collector capacitance distributed effect and AC current crowding effect. The model ignoring base-collector capacitance distributed effect, presents a physical limitation to describe the microwave behavior of the emitter-up HBT structure. If the base-collector capacitance distributed effect and AC current crowding effect are considered in the small-signal equivalent circuit, the method for direct parameter extraction would become different.

Most of the other published small-signal equivalent circuits also account only for pure base resistance to describe the base impedance (7,10-14,16-19). Yet, a simply $R_{c}$ circuit could be adopted to model the AC current crowding, which has been realized in (20,21). A. Oudir et al. (20) extracted the base-spreading capacitance in a T-type small-signal model, however, the small-signal model also does not include the distributed base-collector capacitance effect, which is valid to describe collector-up structure but appears a physical defect for the configuration of emitter-up devices. W. B. Tang et al. (21) considered both the AC current crowding effect and base-collector capacitance distributed effect, however, the small-signal model was developed based on a T-type small-signal equivalent circuit, moreover, there are many simplified approximations in the process of parameter extraction, which leads to inaccurate parameter extraction.

Fig. 1. 3D schematic view of an emitter-up InP HBT structure.

../../Resources/ieie/JSTS.2020.20.4.372/fig1.png

In this study, to overcome the drawbacks mentioned above, a systematic and rigorous extraction procedure for a complete InP HBT π-type small-signal model parameters is proposed, and there is no any simplified approximation step throughout the parameter extraction process. This paper is organized as follows. The integral small-signal model, together with its two-port Y-parameters, are derived and shown in Section II. In order to accurately extract the model elements, the rigorous methodology of deriving closed-form equations is shown in Section III. In Section IV, the parameter extraction results and verification for the small-signal model are given and discussed. Finally, the conclusions are summarized in Section V.

Fig. 2. Complete HBT small-signal equivalent-circuit model.

../../Resources/ieie/JSTS.2020.20.4.372/fig2.png

Table 1. Symbols of all elements in the model

Module

Symbol

Description

Parasitic elements

$C_{pbc}$

Parasitic capacitance between base and collector

$C_{pbe}$

Parasitic capacitance between base and emitter

$C_{pce}$

Parasitic capacitance between collector and emitter

$L_{b}$

Lead inductance associated with base

$L_{c}$

Lead inductance associated with collector

$L_{e}$

Lead inductance associated with emitter

$R_{b}$

Series resistance associated with base

$R_{c}$

Series resistance associated with collector

$R_{e}$

Series resistance associated with emitter

Extrinsic base-collector parasitic

$C_{bcx}$

Extrinsic base-collector capacitance

Intrinsic model

$R_{bi}$

Base-spreading resistance

$C_{bi}$

Base-spreading capacitance

$R_{bc}$

Dynamic base-collector resistance

$R_{be}$

Dynamic base-emitter resistance

$C_{bc}$

Intrinsic base-collector capacitance

$C_{be}$

Intrinsic base-emitter capacitance

$g_{m0}$

Transconductance

τ

Transit time

II. SMALL-SIGNAL MODEL

Fig. 1 displays a 3D representation of an emitter-up InP HBT, associated with small-signal equivalent circuit. In this work, a π-type small-signal model with pad parasitic elements applied to HBTs is presented in Fig. 2. According to the physical characteristics of the device, the completed small signal model is divided into three small modules: parasitic elements, extrinsic base-collector parasitic elements, and intrinsic model elements. The parasitic elements is independent of device operating status. The symbols of all elements in the small-signal model as shown in Fig. 2 are listed in Table 1.

The Y-matrix of the proposed HBT small-signal model, as shown in Fig. 1, is expressed as (1).

(1)
$[Y]=\left\{\left[Y_{M}\right]^{-1}+\left(\begin{array}{cc}Z_{L b}+Z_{L e} & Z_{L e} \\ Z_{L e} & Z_{L c}+Z_{L e}\end{array}\right)\right\}^{-1}$ $+\left(\begin{array}{cc}Y_{p b e}+Y_{p b c} & -Y_{p b c} \\ -Y_{p b c} & Y_{p c e}+Y_{p b c}\end{array}\right)$

with

(2)
$\left[Y_{M}\right]=\left\{\left[Z_{e x}\right]+\left(\begin{array}{cc}Z_{R b}+Z_{R e} & Z_{R e} \\ Z_{R e} & Z_{R c}+Z_{R e}\end{array}\right)\right\}^{-1}$

(3)
$\left[Z_{\mathrm{ex}}\right]=\left\{\left[Z_{\mathrm{in}}\right]^{-1}+Y_{\mathrm{bcx}}\left(\begin{array}{cc}1 & -1 \\ -1 & 1\end{array}\right)\right\}^{-1}$

(4)
$\left[Z_{\mathrm{in}}\right]=\left[Z_{\mathrm{ini}}\right]+Z_{\mathrm{bi}}\left(\begin{array}{ll}1 & 0 \\ 0 & 0\end{array}\right)$

(5)
$\left[Z_{\text {ini }}\right]=\left[Y_{\text {ini }}\right]^{-1}=\left(\begin{array}{cc}\frac{1}{Z_{\mathrm{be}}}+\frac{1}{Z_{\mathrm{bc}}} & -\frac{1}{Z_{\mathrm{bc}}} \\ -\frac{1}{Z_{\mathrm{bc}}}+g_{\mathrm{m}} & \frac{1}{Z_{\mathrm{bc}}}\end{array}\right)^{-1}$

where [$Z_{ex}$] is the Z-matrix of the equivalent circuit including the intrinsic model and extrinsic base-collector parasitic model, [$Z_{in}$] is the Z-matrix of the intrinsic model, [$Z_{ini}$] and [$Y_{ini}$] are the Z-matrix and Y-matrix of the model after peeling $C_{bi}$ and $R_{bi}$ off from the intrinsic model, respectively.

And

(6)
$Z_{\text {bn }}=\frac{R_{\text {bn }}}{1+j \omega R_{\text {bn }} C_{\text {bn }}}, n=i,$ c or e

(7)
$g_{\mathrm{m}}=g_{\mathrm{m0}} \cdot e^{-j \omega \tau}$

(8)
$Y_{\mathrm{bcx}}=j \omega C \mathrm{bex}$

(9)
$Z_{\mathrm{Rn}}=R_{\mathrm{n}}, \mathrm{n}=\mathrm{b}, \mathrm{c}$ or e

(10)
$Z_ \mathrm{Ln}=j \omega L_{\mathrm{n}}, \mathrm{n}=\mathrm{b}, \mathrm{c}$ or e

(11)
$Y_{\mathrm{eXY}}=j \omega C_{\mathrm{pXY}}, \mathrm{XY}=\mathrm{bc},$ be or ce

The extraction steps for the equivalent-circuit elements parameters only depending on S-parameters measured data will be depicted in the following section.

III. PARAMETER EXTRACTION PROCEDURE

To accurately and intuitively determine the model elements, the exaction equations are derived from S-parameters by peeling peripheral elements from small-signal models to get reduced ones.

A. Determination of Parasitic Elements

The parasitic elements are independent of the bias conditions. To obtain parasitic pad parameters, one could employ open and short test structures or cutoff cold-HBT method. However, the cutoff operation method has insufficient accuracy (16). Thus, the open and short technique was used to determine parasitic pad parameters in this study (10).

The series resistances ($R_{b}$, $R_{c}$, and $R_{e}$) were determined adopting the commonly-used open-collector method.

B. Determination of Extrinsic Base-collector Parasitic Elements

The series resistances ($R_{b}$, $R_{c}$, and $R_{e}$), and parasitic pad elements ($C_{pbe}$, $C_{pbc}$, $C_{pce}$, $L_{b}$, $L_{c}$, and $L_{e}$) are removed from the HBT small-signal model as follows which is concluded from (1) and (2).

(12)
$\begin{aligned}\left[Z_{\mathrm{ex}}\right] &=\left\{\begin{array}{lc}\left\{[Y]-\left(\begin{array}{cc}Y_{\mathrm{pbe}}+Y_{\mathrm{pbc}} & -Y_{\mathrm{pbc}} \\ -Y_{\mathrm{pbc}} & Y_{\mathrm{pce}}+Y_{\mathrm{pbc}}\end{array}\right)\right\}^{-1} \\ -\left(\begin{array}{cc}Z_{\mathrm{Lb}}+Z_{\mathrm{Le}} & Z_{\mathrm{Le}} \\ Z_{\mathrm{Le}} & Z_{\mathrm{Lc}}+Z_{\mathrm{Le}}\end{array}\right) \\ -\left(\begin{array}{cc}Z_{\mathrm{Rb}}+Z_{\mathrm{Re}} & Z_{\mathrm{Re}} \\ Z_{\mathrm{Re}} & Z_{\mathrm{Rc}}+Z_{\mathrm{Re}}\end{array}\right)\end{array}\right\} \end{aligned}$

Combining (3), (4) and (5), we can calculate the Y-matrix of the equivalent circuit including the intrinsic model and extrinsic base-collector parasitic model, and it is shown as (13).

(13)
$\left[Y_{\operatorname{ex}}\right]=\left(\begin{array}{cc}\frac{Z_{\mathrm{bc}}}{D}+\frac{Z_{\mathrm{be}}}{D}+Y_{\mathrm{bcx}} & -\frac{Z_{\mathrm{be}}}{D}-Y_{\mathrm{bcx}} \\ -\frac{Z_{\mathrm{be}}}{D}+\frac{Z_{\mathrm{bc}} \cdot Z_{\mathrm{be}} \cdot g_{\mathrm{m}}}{D}-Y_{\mathrm{bcx}} & \frac{g_{\mathrm{m}} \cdot Z_{\mathrm{be}} \cdot Z_{\mathrm{bi}}+Z_{\mathrm{be}}+Z_{\mathrm{bi}}}{D}+Y_{\mathrm{bcx}}\end{array}\right)$

with D=$R_{bi}$$Z_{be}$+$R_{bi}$$Z_{bc}$+$Z_{bc}$$Z_{be}$. After some calculations, using the Y-parameters given in (13), $Z_{bi}$ can be extracted which is expressed as:

(14)
$Z_{\mathrm{bi}}=\frac{Y_{\mathrm{T}}}{\Delta Y-Y_{\mathrm{bex}} \cdot \Sigma Y}$

where $Y_{T}$=$Y_{ex,12}$+$Y_{ex,22}$, ΔY=$Y_{ex,11}$$Y_{ex,22}$-$Y_{ex,12}$$Y_{ex,21}$, ∑Y= $Y_{ex,11}$+$Y_{ex,12}$+$Y_{ex,21}$+$Y_{ex,22}$, and $Y_{ex,ij}$ represent the Y-parameters of Y-matrix [$Y_{ex}$].

By referring to (6) and (14), adopting some simple analysis, we note that the base-spreading resistance $R_{bi}$ and the base-spreading capacitance $C_{bi}$ can be determined as follows:

(15)
$R_{\mathrm{bi}}=\frac{1}{\operatorname{real}\left[\left(\Delta Y-Y_{\mathrm{bcx}} \cdot \sum Y\right) / Y_{\mathrm{t}}\right]}$

(16)
$C_{\mathrm{bi}}=\frac{\operatorname{imag}\left[\left(\Delta Y-Y_{\mathrm{bex}} \cdot \Sigma Y\right) / Y_{\mathrm{t}}\right]}{\omega}$

The Y-matrix [$Y_{ex}$] is transformed to Z-matrix [$Z_{ex}$]. It can be expressed as (17).

(17)
$\left[Z_{\mathrm{ex}}\right]=\left(\begin{array}{cc}\frac{g_{\mathrm{m}} Z_{\mathrm{be}} Z_{\mathrm{bi}}+Z_{\mathrm{be}}+Z_{\mathrm{bi}}+D Y_{\mathrm{bcx}}}{A} & \frac{Z_{\mathrm{be}}+D Y_{\mathrm{bcx}}}{A} \\ \frac{Z_{\mathrm{be}}+D Y_{\mathrm{bcx}}-Z_{\mathrm{bc}} Z_{\mathrm{be}} g_{m}}{A} & \frac{Z_{\mathrm{be}}+Z_{\mathrm{bc}}+D Y_{\mathrm{bcx}}}{A}\end{array}\right)$

with A=(1+$g_{m}$$Z_{be}$)[1+$Y_{bcx}$($Z_{bi}$+$Z_{bc}$)], After some calculations from (17), the following expression can be deduced:

(18)
$Y_{\mathrm{bcx}}=j \omega C_{\mathrm{bcx}}=\frac{Z_{\mathrm{bi}}-Z_{\mathrm{N}}}{Z_{\mathrm{bi}} Z_{\mathrm{M}}}$

with ZM=$Z_{ex,11}$+$Z_{ex,22}$-$Z_{ex,12}$-$Z_{ex,21}$, and ZN=$Z_{ex,11}$-$Z_{ex,12}$.

After substituting (14) into (18), the $Y_{bcx}$ can be rewritten as:

(19)
$\begin{aligned} Y_{\text {bex }}=& \frac{\text { real }\left(Z_{\mathrm{M}}\right)-j \cdot \operatorname{imag}\left(Z_{\mathrm{M}}\right)}{B}-\frac{L V+\omega T V C_{\text {bex }}-X N+\omega X R C_{\text {bex }}}{B Q} \\ &-j \cdot \frac{L X+\omega T X C_{\text {bex }}+V N-\omega V R C_{\text {bex }}}{B Q} \end{aligned}$

where

(20)
$V=\operatorname{real}\left(Z_{\mathrm{N}}\right) \cdot \operatorname{real}\left(Z_{\mathrm{M}}\right)+\operatorname{imag}\left(Z_{\mathrm{N}}\right) \cdot \operatorname{imag}\left(Z_{\mathrm{M}}\right)$

(21)
$X=\operatorname{imag}\left(Z_{\mathrm{N}}\right) \cdot \operatorname{real}\left(Z_{\mathrm{M}}\right)-\operatorname{real}\left(Z_{\mathrm{N}}\right) \cdot \operatorname{imag}\left(Z_{\mathrm{M}}\right)$

(22)
$B=\left[\operatorname{real}\left(Z_{\mathrm{M}}\right)\right]^{2}+\left[\operatorname{imag}\left(Z_{\mathrm{M}}\right)\right]^{2}$

Fig. 3. Extracted $C_{bcx}$ versus frequency.

../../Resources/ieie/JSTS.2020.20.4.372/fig3.png

(23)
$Q=\left[\operatorname{real}\left(Y_{\mathrm{T}}\right)\right]^{2}+\left[\operatorname{imag}\left(Y_{\mathrm{T}}\right)\right]^{2}$

(24)
$R=\operatorname{real}(\Sigma Y) \cdot \operatorname{real}\left(Y_{\mathrm{T}}\right)+\operatorname{imag}(\Sigma Y) \cdot \operatorname{imag}\left(Y_{\mathrm{T}}\right)$

(25)
$L=\operatorname{real}(\Delta Y) \cdot \operatorname{real}\left(Y_{\mathrm{T}}\right)+\operatorname{imag}(\Delta Y) \cdot \operatorname{imag}\left(Y_{\mathrm{T}}\right)$

(26)
$T=\operatorname{imag}(\Sigma Y) \cdot \operatorname{real}\left(Y_{\mathrm{T}}\right)-\operatorname{real}(\Sigma Y) \cdot \operatorname{imag}\left(Y_{\mathrm{T}}\right)$

(27)
$N=\operatorname{imag}(\Delta Y) \cdot \operatorname{real}\left(Y_{\mathrm{T}}\right)-\operatorname{real}(\Delta Y) \cdot \operatorname{imag}\left(Y_{\mathrm{T}}\right)$

Finally, through making the real part of (19) equal to zero, the extrinsic base-collector capacitance $C_{bcx}$ can be expressed as:

(28)
$C_{\mathrm{bex}}=\frac{Q \cdot \operatorname{real}\left(Z_{\mathrm{M}}\right)-L V+X N}{\omega X R+\omega T V}$

The $C_{bcx}$ are determined using (28), and the $C_{bcx}$ characteristics, at $V_{CE}$=2.7 V and $I_{C}$=15.5 mA, is shown in Fig. 3. The average value over the entire frequency range is assigned to $C_{bcx}$.

C. Determination of Intrinsic Model Elements

After peeling off $C_{bcx}$, the device model is reduced to the intrinsic HBT model. From (3), the Y-matrix for the intrinsic HBT model can be given by (29).

(29)
$\left[Y_{\mathrm{in}}\right]=\left[Y_{\mathrm{ex}}\right]-Y_{\mathrm{bcx}}\left(\begin{array}{cc} 1 & -1 \\ -1 & 1 \end{array}\right)$

Next, if the $R_{c}$ circuit consisting of $R_{bi}$ and $C_{bi}$ is also peeled off, the resultant intrinsic HBT model called INI model is left. The values of $R_{bi}$ and $C_{bi}$ can be evaluated using (15) and (16), respectively. The extracted values of $R_{bi}$ and $C_{bi}$ for the bias condition ($V_{CE}$=2.7 V, $I_{C}$=15.5 mA) are presented in Fig. 4 and 5, respectively.

From (4), the two-port Z-matrix of the INI model is given as (30).

(30)
$\left[Z_{\mathrm{ini}}\right]=\left[Y_{\mathrm{in}}\right]^{-1}-\left(\begin{array}{cc}Z_{\mathrm{bi}} & 0 \\ 0 & 0\end{array}\right)$

The Z-matrix for the INI model is transformed to Y-matrix [$Y_{ini}$]. It can be written by (31).

(31)
$\left[Y_{\mathrm{ini}}\right]=\left[Z_{\mathrm{ini}}\right]^{-1}=\left(\begin{array}{cc}\frac{1}{Z_{\mathrm{be}}}+\frac{1}{Z_{\mathrm{bc}}} & -\frac{1}{Z_{\mathrm{bc}}} \\ -\frac{1}{Z_{\mathrm{bc}}}+g_{\mathrm{m}} & \frac{1}{Z_{\mathrm{bc}}}\end{array}\right)$

After some calculations, using the Y-matrix given in (31) and referring to (6)-(7), the following equations are obtained to determine the intrinsic elements ($R_{bc}$, $C_{bc}$, $R_{be}$, $C_{be}$, $g_{m0}$ and τ).

(32)
$R_{\mathrm{bc}}=-\frac{1}{\operatorname{real}\left(Y_{\mathrm{ini}, 12}\right)}$

(33)
$C_{\mathrm{bc}}=-\frac{\operatorname{imag}\left(Y_{\mathrm{ini}, 12}\right)}{\omega}$

(34)
$R_{\mathrm{be}}=\frac{1}{\operatorname{real}\left(Y_{\mathrm{ini}, 11}+Y_{\mathrm{ini}, 12}\right)}$

(35)
$C_{\mathrm{be}}=\frac{\operatorname{imag}\left(Y_{\mathrm{ini}, 11}+Y_{\mathrm{ini}, 12}\right)}{\omega}$

(36)
$g_{\mathrm{m} 0}=\operatorname{mag}\left(Y_{\mathrm{ini}, 21}-Y_{\mathrm{ini}, 12}\right)$

(37)
$\tau=\frac{-1}{\omega} \tan ^{-1}\left[\frac{\operatorname{imag}\left(Y_{\mathrm{ini}, 21}-Y_{\mathrm{ini}, 12}\right)}{\operatorname{real}\left(Y_{\mathrm{ini}, 21}-Y_{\mathrm{ini}, 12}\right)}\right]$

where $Y_{ini}$,ij are the Y-parameters of the INI equivalent circuit. Eqs. (32)-(37) permit to extract the INI model elements ($R_{bc}$, $C_{bc}$, $R_{be}$, $C_{be}$, $g_{m0}$ and τ) in case the Y-parameters of the INI model are calculated using (29)-(31). Fig. 6-11 show the extracted INI model parameters versus frequency for a 1×15 µm$^{2}$ InP HBT biased at $V_{CE}$=2.7 V and $I_{C}$=15.5 mA.

Fig. 4. Extracted $R_{bi}$ versus frequency.

../../Resources/ieie/JSTS.2020.20.4.372/fig4.png

Fig. 5. Extracted $C_{bi}$ versus frequency.

../../Resources/ieie/JSTS.2020.20.4.372/fig5.png

Fig. 6. Extracted $R_{bc}$ versus frequency.

../../Resources/ieie/JSTS.2020.20.4.372/fig6.png

Fig. 7. Extracted $R_{be}$ versus frequency.

../../Resources/ieie/JSTS.2020.20.4.372/fig7.png

Fig. 8. Extraction of $C_{bc}$.

../../Resources/ieie/JSTS.2020.20.4.372/fig8.png

Fig. 9. Extraction of $C_{be}$.

../../Resources/ieie/JSTS.2020.20.4.372/fig9.png

Fig. 10. Extracted $g_{m0}$ versus frequency.

../../Resources/ieie/JSTS.2020.20.4.372/fig10.png

Fig. 11. Extracted τ versus frequency.

../../Resources/ieie/JSTS.2020.20.4.372/fig11.png

IV. RESULTS AND DISCUSSIONS

In order to validate the presented procedure, an n-p-n emitter-up InP HBT with a 1×15 μm$^{2}$ emitter area was investigated. The S-parameters measured data were realized on-wafer thought employing Keysight 8510C Network Analyzer with Keysight B1500A Semicon-ductor Device Analyzer to supply DC source. With all measured instruments by the control of IC-CAP control software, S-parameters of the range of 0.1-40 GHz were obtained.

Table 2. Extracted parameters for small-signal π-model at $V_{CE}$=2.7 V, $I_{C}$=15.5 mA

Elements

Proposed model

Model in (15)

Bias

independent

parameters

$R_{e}$ (Ω)

4.27

4.27

$R_{b}$ (Ω)

1.77

1.77

$R_{c}$ (Ω)

7.31

7.31

Bias

dependent

parameters

$C_{bcx}$ (fF)

10.07

-

$R_{bi}$ (Ω)

8.26

4.91

$C_{bi}$ (fF)

280.05

-

$R_{bc}$ (kΩ)

17.26

18.41

$R_{be}$ (Ω)

28.59

32.17

$C_{bc}$ (fF)

13.52

16.11

$C_{be}$ (pF)

1.89

1.98

$g_{m0}$

3.17

3.07

τ (ps)

0.22

0.23

Applying the extraction method presented in this study, extracted values of the small-signal model parameters, as well as for the model in (15) (which does not include $C_{bcx}$ and $C_{bi}$), for the 1×15 μm$^{2}$ InP HBT under two bias points of ($V_{CE}$=2.7 V, $I_{C}$=15.5 mA) and ($V_{CE}$=2.5 V, $I_{C}$=17.5 mA) are summarized in Table 2 and 3, respectively. Fig. 12 shows the comparisons between the measured and modeled S-parameters data for two bias points of ($V_{CE}$=2.7 V, $I_{C}$=15.5 mA) and ($V_{CE}$=2.5 V, $I_{C}$=17.5 mA). From Fig. 12, it can be observed that the presented model was found to improve the prediction of S12 and S22, at the higher frequencies. Additionally, the influence of the base-emitter and collector-emitter metallization capacitances on model accuracy was found to be very limited. The residual error (15) is adopted to quantify the accuracy of the extraction method and the residual errors between the measured and modeled S-parameters at different bias conditions are less than 2.8%, which is better than the other ones given in previously published literatures (13,15,17,20,22). The residual errors between the measured and modeled S-parameters, at different bias conditions, for the presented model and the model in (15) are listed in Table 4. It is shown that the smaller error values are obtained for the model with $C_{bcx}$ and $C_{bi}$ in this article. This comparison indicates that this model presents a better agreement with the experimental data for different bias points. Therefore, this result shows that the small-signal model and parameter extraction method have high accuracy to model HBTs characteristics.

Table 3. Extracted parameters for small-signal π-model at $V_{CE}$=2.5 V, $I_{C}$=17.5 mA

Elements

Proposed model

Model in (15)

Bias

independent

parameters

$R_{e}$ (Ω)

4.27

4.27

$R_{b}$ (Ω)

1.77

1.77

$R_{c}$ (Ω)

7.31

7.31

Bias

dependent

parameters

$C_{bcx}$ (fF)

16.84

-

$R_{bi}$ (Ω)

14.35

3.14

$C_{bi}$ (fF)

526.13

-

$R_{bc}$ (kΩ)

17.49

18.18

$R_{be}$ (Ω)

34.89

37.53

$C_{bc}$ (fF)

11.65

15.65

$C_{be}$ (pF)

1.13

0.85

$g_{m0}$

1.92

1.71

τ (ps)

0.36

0.38

Fig. 12. Measured and modeled S-parameters on Smith plot from 0.1 GHz to 40 GHz for (a) $V_{CE}$=2.7 V and $I_C$=15.5 mA, (b) $V_{CE}$=2.5 V and $I_C$=17.5 mA.

../../Resources/ieie/JSTS.2020.20.4.372/fig12.png

Table 4. Residual error for the extracted bias points

Bias points

$V_{CE}$=2.7 V,

$I_{C}$=15.5 mA

$V_{CE}$=2.5 V,

$I_{C}$=17.5 mA

Residual error (%)

Proposed model

1.52

2.75

Model in (15)

5.15

4.48

V. CONCLUSION

ACKNOWLEDGMENTS

This work was supported by the National Natural Science Foundation of China (Grant No. 61804046), and the Foundation of Department of Science and Technology of Henan Province (Grant No. 202102210322).

REFERENCES

1 
Zhang J. C., Liu B., Liu M., Zhang L. W., Wang J. C., Jin H., Aug 2019, A Watt-level Broadband Power Amplifier in GaAs HBT Process, Journal of Semiconductor Technology and Science, Vol. 19, No. 4, pp. 357-363DOI
2 
Urteaga M., Griffith Z., Pierson R., Rowell P., Young A., Hacker J., Brar B., Kim S.K., Maurer R., Rodwell M.J.W, Oct 2017, THz InP bipolar transistors-circuit integration and applications, 2017 IEEE Compound Semiconductor Integrated Circuit Symposium (CSICS), pp. 1-4DOI
3 
Zhang J. C., Liu M., Zhang L. W., Wang J. C., Liu B., Sep 2019, Design of a low-phase-noise Ka-band GaAs HBT VCO, Journal of Circuits Systems and Computers, Vol. 28, No. 10, pp. 1950174Google Search
4 
Liu M., Xu P., Zhang J. C., Liu B., Zhang L. W., Apr 2020, A 4.2-to-5.4 GHz stacked GaAs HBT power amplifier for C-band applications, Circuit WorldDOI
5 
Zhang J. C., Liu M., Wang J. C., Zhang L. W., Liu B., Jul 2019, Modeling of InP HBTs with an improved Keysight HBT model, Microw. J., Vol. 62, No. 7, pp. 56-62Google Search
6 
C.Zhang J., M.Zhang Y., L.u L H., Zhang Y. M., Liu B., Zhang L. M., Wang J. C., Zhang L. W., 2016, A Ku-band low-phase-noise cross-coupled VCO in GaAs HBT technology, Journal of Circuits, Systems, and computers, Vol. 25, No. 6, pp. 1650053DOI
7 
Taher H., Mar 2012, Direct extraction technique of π-topology small-signal equivalent circuit model for Si/SiGe heterojunction bipolar transistor, Microw. Optical Tech. Lett., Vol. 54, No. 3, pp. 584-589DOI
8 
Johansen T. K., Leblanc R., Poulain J., Delmouly V., Jan 2016, Direct extraction of InP/GaAsSb/InP DHBT equivalent-circuit elements from S–parameters measured at cut-off and normal bias conditions, IEEE Trans. Microw. Theory Tech., Vol. 64, No. 1, pp. 115-124DOI
9 
Álvarez-Botero G., Torres-Torres R., Murphy-Arteaga R. S., Dec 2015, Consistent modeling and power gain analysis of microwave SiGe HBTs in CE and CB configurations, IEEE Trans. Microw. Theory Tech., Vol. 63, No. 12, pp. 3888-3895DOI
10 
Deng M., Quémerais T., Bouvot S., Gloria D., Chevalier P., Lépilliet S., Danneville F., Dambrine G., Mar 2017, Small-signal characterization and modelling of 55 nm SiGe BiCMOS HBT up to 325 GHz, Solid State Electron., Vol. 129, No. 3, pp. 150-156DOI
11 
Sun Y., Fu J., Wang Y., Zhou W., Zhang W., Liu Z., Jun 2016, Direct analytical parameter extraction for SiGe HBTs T-topology small-signal equivalent circuit, Superlattices and Microstructures, Vol. 94, No. 6, pp. 223-230DOI
12 
Gao J., Li X., Wang H., Boeck G., Feb 2006, An approach to determine small-signal model parameters for InP-based heterojunction bipolar transistors, IEEE Trans. Semiconduct. Manufact., Vol. 19, No. 1, pp. 138-145DOI
13 
Oudir A., Mahdouani M., Mansouri S., Bourguiga R., Pardo F., Pelouard J. L., Jan 2010, Small-signal modeling of emitter-up HBT using an improved analytical approach. Application to InGaAlAs/GaAsSb/InP DHBT with strained base, Solid State Electron., Vol. 54, No. 1, pp. 67-78DOI
14 
Lee K., Choi K., Kook S. H., Cho D. H., Park K. W., Kim B., Mar 2005, Direct parameter extraction of SiGe HBTs for the VBIC bipolar compact model, IEEE Trans. Electron. Devices, Vol. 52, No. 3, pp. 375-384DOI
15 
Zhang J. C., Liu B., Zhang L. W., Sun L. G., Wang J. W., Song L., Dec 2015, A rigorous peeling algorithm for direct parameter extraction procedure of HBT small-signal equivalent circuit, Analog Integr. Circ. Sig. Process, Vol. 85, No. 3, pp. 405-411DOI
16 
Zhang A., Gao J., Dec 2018, A new method for determination of PAD capacitances for GaAs HBTs based on scalable small signal equivalent circuit model, Solid State Electron., Vol. 150, No. 12, pp. 45-50DOI
17 
Oudir A., Mahdouani M., Bourguiga R., Aug 2011, Direct extraction method of HBT equivalent-circuit elements relying exclusively on S-parameters measured at normal bias conditions, IEEE Trans. Microw. Theory Tech., Vol. 59, No. 8, pp. 1973-1982DOI
18 
Johansen T. K., Krozer V., Nodjiadjim V., Konczykowska A., Dupuy J. Y., Riet M., Sep 2011, Improved external base resistance extraction for submicrometer InP/InGaAs DHBT models, IEEE Trans. Electron. Devices, Vol. 58, No. 9, pp. 3004-3011DOI
19 
Sun Y., Li X., Zhang J., Shi Y., Sep 2017, Improved high-frequency equivalent circuit model based on distributed effects for SiGe HBTs with CBE layout, Chin. Phys. B, Vol. 26, No. 9, pp. 098502Google Search
20 
Oudir A., Mahdouani M., Bourguiga R., Pardo F., Pelouard J. L., Nov 2008, An analytic procedure for extraction of metallic collector-up InP/InGaAsP/ InGaAs HBT small signal equivalent circuit parameters, Solid State Electron., Vol. 52, No. 11, pp. 1742-1750DOI
21 
Tang W. B., Wang C. M., Hsin Y. M., Oct 2006, A new extraction technique for the complete small-signal equivalent-circuit model of InGaP/GaAs HBT including base contact impedance and AC current crowding effect, IEEE Trans. Microw. Theory Tech., Vol. 54, No. 10, pp. 3641-3647DOI
22 
Olvera-Cervantes J. L., Cressler J. D., Medina-Monroy J. L., Thrivikraman T., Banerjee B., Laskar J., Mar 2008, A new analytical method for robust extraction of the small-signal equivalent circuit for SiGe HBTs operating at cryogenic temperatures, IEEE Trans. Microw. Theory Tech., Vol. 56, No. 3, pp. 568-574DOI

Author

Jincan Zhang
../../Resources/ieie/JSTS.2020.20.4.372/au1.png

Jincan Zhang was born in Xingtai, China, in 1985.

He received the M.S. degree in Xi’an University of Technology, Xi’an, China, in 2010.

He received the Ph.D. degree in XiDian University, Xi’an, China, in June 2014.

Now He is an associate professor in Henan university of science and technology, Luoyang, China.

His research is focused on modeling of HBTs and design of very high speed integrated circuit.

Leiming Zhang
../../Resources/ieie/JSTS.2020.20.4.372/au2.png

Leiming Zhang was born in Luoyang,China, in 1980.

He received the M.S. degree in University of Electronic Science and Technology of China(UESTC) , Chengdu, in 2008.

Now he is a lecturer in Henan University of Science and Technology(HAUST), Luoyang, China.

His research is focused on device modeling of CMOS and design of mixed signal integrated circuits.

Min Liu
../../Resources/ieie/JSTS.2020.20.4.372/au3.png

Min Liu was born in Baoding, China, in 1984.

She received the Ph.D. degree in XiDian University, Xi’an, China, in June 2016.

Now She is a lecturer in Henan university of science and technology, Luoyang, China.

Her research is focused on modeling of HBTs and design of integrated circuits.

Liwen Zhang
../../Resources/ieie/JSTS.2020.20.4.372/au4.png

Liwen Zhang was born in Xinyang, He Nan, China in December 1980.

She received the B.S. degree in Optoelectronic technology in 2001 and the M.S. degree in Condensed matter physics in 2004 from Zhengzhou University, Zhengzhou, and received her Ph.D. degree in Atomic and Molecular Physics at Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan, in 2008.

She is currently a Professor in the College of Electrical Engineering, Henan University of Science and Technology.

Her major field is Interconnection packaging reliability, structural stress analysis in three-dimensional package devices and Semiconductor optoelectronic technology.