• Journal of Semiconductors
  • Vol. 42, Issue 5, 052401 (2021)
Jinjing Huang and Jun Liu
Author Affiliations
  • Key Laboratory of RF Circuits and Systems, Ministry of Education, Hangzhou Dianzi University, Hangzhou 310018, China
  • show less
    DOI: 10.1088/1674-4926/42/5/052401 Cite this Article
    Jinjing Huang, Jun Liu. A complete small-signal HBT model including AC current crowding effect[J]. Journal of Semiconductors, 2021, 42(5): 052401 Copy Citation Text show less

    Abstract

    An improved small-signal equivalent circuit of HBT concerning the AC current crowding effect is proposed in this paper. AC current crowding effect is modeled as a parallel RC circuit composed of Cbi and Rbi, with distributed base-collector junction capacitance also taken into account. The intrinsic portion is taken as a whole and extracted directly from the measured S-parameters in the whole frequency range of operation without any special test structures. An HBT device with a 2 × 20 μm2 emitter-area under three different biases were used to demonstrate the extraction and verify the accuracy of the equivalent circuit.

    1. Introduction

    Heterojunction bipolar transistors (HBTs) have been extensively used in high-speed analog, digital and mixed-signal integrated circuits (ICs)[1-4]. There are several small-signal equivalent circuit topologies for microwave HBT devices’ modeling and the accuracy of small-signal model plays a crucial role in the design of ICs.

    However, most of the reported topologies ignore the AC current crowding effect[5], which is a critical factor for predicting the AC performance of HBT devices, especially in high frequency range. AC current crowding occurs in the high frequency range related to the bypassing effects of base-emitter capacitance. As the small-signal voltage drops along the intrinsic base region, the diffusion capacitance at the edge of the emitter is proportionally higher than that at the center. This phenomenon can cause the small-signal emitter current to crowd near the edge of the base-emitter junction as frequency increases. Therefore, the intrinsic base capacitance is physically connected with the AC current crowding[6].

    In recent years, different small-signal equivalent circuits have been proposed with various extraction methods, in which π[7] or T[8, 9] topologies were used to represent the small-signal equivalent circuits. Heung et al.[10] proposed a simple extraction method by using the RC circuit to characterize the AC current crowding effect. However, a simple π-type topology was adopted to extract the intrinsic base capacitance (Cbi) through the Z-parameter equations. In addition, other papers[11-14] have proposed a more complete small-signal model compared with Ref. [10]. Although their methods have shown good accuracy in the extraction process, more cumbersome analytical methods and calculation cost were needed. Zhang et al.[15] proposed a rigorous but concise peeling extraction algorithm. While this method failed to consider the AC current crowding effect, this paper proposes an accurate extraction method with a complete π-type small-signal equivalent circuit considering the AC current crowding effect. This new equivalent-circuit topology is developed based on the small-signal equivalent circuit of an AHBT (agilent heterojunction bipolar transistor) model which has considered the distributed base-collector junction capacitance. Based on the peeling algorithm in Ref. [15] and the T-π transformation presented in Refs. [1618], a novel method is proposed. In Section 2, an integral small-signal model taking into account the AC current crowding effect is introduced. The parameter extraction procedure is presented in detail in Section 3. Section 4 validates the model performance and a comparison between the proposed model and the conventional model without Cbi. Finally, a conclusion is given in Section 5.

    2. Small-signal model

    The complete small-signal equivalent circuit of the HBT device considering the AC current crowding effect is given in Fig. 1.

    Complete small-signal equivalent circuit of HBT device including . .

    Figure 1.Complete small-signal equivalent circuit of HBT device including . .

    In Fig. 1, Rbx, Rcx, Re are the base, collector, and emitter parasitic resistances, respectively. Rbcx and Rbex are the extrinsic base-collector and base-emitter resistances, and Cbcx and Cbex are the extrinsic base-collector and base-emitter depletion capacitances. Rbi and Rci are the intrinsic resistances, and Cbi is the intrinsic base capacitance. The parallel RC circuit composed of Cbi and Rbi characterizes the AC current crowding effect. Rbci, Rbei, Cbci and Cbei are the intrinsic resistances and capacitances, respectively. gm and gm0 are the small-signal and DC transconductances, respectively. τ is the delay time.

    Circuit topology demonstrated in Fig. 1 can be simplified to facilitate calculation. We introduce three parameters, Arbe, Arc and Abex with a value range of (0, 1) and stipulate , , , , , and .

    3. Extraction procedure

    3.1. Parasitic resistances

    Since the parasitic resistances Rbx, Rc and Re are independent of the bias condition, we extract them under zero biasing conditions in this work. The calculation equations are given as follows[19]: , and . The extracted Rbx, Rc and Re for the 2 × 20 μm2 HBT device are 6.186, 2.561, and 5.526 Ω, respectively. The impedance matrix Zm after de-embedding Rbx, Rc and Re from the total two-port Z-parameters Z can be written as

    $ \left[ {{Z_{\rm{m}}}} \right] = Z-\left[ {\begin{array}{*{20}{c}} {{R_{{\rm{bx}}}} + {R_{\rm{e}}}} & {{R_{\rm{e}}}}\\[2mm] {{R_{\rm{e}}}} & {{R_{\rm{c}}} + {R_{\rm{e}}}} \end{array}} \right]. $  (1)

    Through Y–Z transformation, Ym can also be expressed as

    $ {Y_{{\rm{m11}}}} = \frac{{{Z_{{\rm{be}}}} + {Z_{{\rm{bci}}}}}}{N} + {Y_{{\rm{bcx}}}}, $  (2)

    $ {Y_{{\rm{m12}}}} = \frac{{-{Z_{{\rm{be}}}}}}{N}-{Y_{{\rm{bcx}}}}, $  (3)

    $ {Y_{{\rm{m21}}}} = \frac{{{{\rm{g}}_{\rm{m}}}{Z_{{\rm{be}}}}{Z_{{\rm{bci}}}}-{Z_{{\rm{be}}}}}}{N}-{Y_{{\rm{bcx}}}}, $  (4)

    $ {Y_{{\rm{m22}}}} = \frac{{{Z_{{\rm{be}}}} + {Z_{{\rm{bi}}}}\left( {1 + {{\rm{g}}_{\rm{m}}}{Z_{{\rm{be}}}}} \right)}}{N} + {Y_{{\rm{bcx}}}}, $  (5)

    where , , , , .

    3.2. Intrinsic part parameters

    Once the parasitic resistances are known, the extraction of the intrinsic parameters can be carried out[20]. The circuit is presented in Fig. 2, and the admittance parameter of the intrinsic part Yin can be written as

    Small-signal equivalent circuit after de-embedding the extrinsic parameters.

    Figure 2.Small-signal equivalent circuit after de-embedding the extrinsic parameters.

    $ \left[ {{Y_{{\rm{in}}}}} \right] = \left[ {{Y_{\rm{m}}}} \right]-\left[ {\begin{array}{*{20}{c}} 0 & 0\\ 0 & {{Y_{\rm{o}}}} \end{array}} \right]. $  (6)

    The circuit in Fig. 2 can be converted to that in Fig. 3 after T-π transformation[21] with . The parameters Z1, Z2, and Z3 can be derived from Eq. (6) as follows

    Final circuit after T-π transformation.

    Figure 3.Final circuit after T-π transformation.

    $ {{Z}}_{{1}} = \frac{{N}}{{{Z}}_{\rm{bci}}}{,} $  (7)

    $ {{Z}}_{{2}} = \frac{{N}}{{{Z}}_{\rm{be}}}{,} $  (8)

    $ {{Z}}_{{3}} = \frac{{N}}{{{Z}}_{\rm{bi}}}{.} $  (9)

    Combining Eqs. (6)–(9), Yin can also be represented by Z1, Z2, Z3, and Z4, which are expressed as[16]

    $ \left[ {{Y_{{\rm{in}}}}} \right] = \left[ {\begin{array}{*{20}{c}} {\dfrac{1}{{{Z_1}}} + \dfrac{1}{{{Z_4}}}} & {-\dfrac{1}{{{Z_4}}}}\\ {X \dfrac{{{Z_3}}}{{{Z_1}}}-\dfrac{1}{{{Z_4}}}} & {\dfrac{1}{{{Z_3}}} + \dfrac{1}{{{Z_4}}} + X} \end{array}} \right], $  (10)

    in which

    $ X = B {{{g}}_{\rm{m}}}, $  (11)

    $ B = \frac{{{Z_1}}}{{{Z_1} + {Z_2} + {Z_3}}}. $  (12)

    Then, we can get that

    $ {{Z}}_{{1}} = \frac{{1}}{{{Y}}_{\rm{in11}}+{{Y}}_{\rm{in12}}}{,} $  (13)

    $ {{Z}}_{{3}} = \frac{{{Y}}_{\rm{in11}}+{{}{Y}}_{\rm{in21}}}{\left({{Y}}_{\rm{in11}}+{{Y}}_{\rm{in12}}\right)\left({{Y}}_{\rm{in22}}+{{Y}}_{\rm{in12}}\right)}{,} $  (14)

    $ {{Z}}_{{4}} = {-}\frac{{1}}{{{Y}}_{\rm{in12}}}{.} $  (15)

    Based on the above relationship between Yin and Z1, Z2, Z3, the extraction method reported in Ref. [18] is used to extract the parameters. From Eqs. (13) and (14), we can get

    $ \frac{{{Z_1}}}{{{Z_3}}} = \frac{{{R_{{\rm{bi}}}}}}{{{R_{{\rm{bci}}}}}}\frac{{1 + j\omega {R_{{\rm{bci}}}}{C_{{\rm{bci}}}}}}{{1 + j\omega {R_{{\rm{bi}}}}{C_{{\rm{bi}}}}}}, $  (16)

    $ {\rm{real}}\left( {\frac{{{Z_1}}}{{{Z_3}}}\left( {1 + j\omega {T_{{\rm{bi}}}}} \right)} \right) = \frac{{{R_{{\rm{bi}}}}}}{{{R_{{\rm{bci}}}}}}, $  (17)

    $ {\rm{imag}}\left( {\frac{{{Z_1}}}{{{Z_3}}}\left( {1 + j\omega {T_{{\rm{bi}}}}} \right)} \right) = \omega {R_{{\rm{bi}}}}{C_{{\rm{bci}}}}. $  (18)

    Defining

    $ \begin{array}{l} {T_{{\rm{bi}}}} = {R_{{\rm{bi}}}}{C_{{\rm{bi}}}},\\ \;\;\;\;{F_0} = \dfrac{\omega }{{{\rm{imag}}\left( {\dfrac{{{Z_1}}}{{{Z_3}}}} \right)}} = {A_0} + {\omega ^2}{B_0}, \end{array} $  (19)

    where , , .

    Since Z1 and Z3 are given in Eqs. (13) and (14), the ω2 dependence of F0 can be easily plotted in Fig. 4. Then A0 and B0 can be determined. Tbi is defined as

    Frequency of F0 versus ω2.

    Figure 4.Frequency of F0 versus ω2.

    $ {{T}}_{\rm{bi}} = \sqrt{\frac{{{B}}_{{0}}}{{{A}}_{{0}}}} . $ (20)

    Based on Eq. (7), Z1 can be written as

    $ \begin{array}{l} {Z_1} = \dfrac{{{R_{{\rm{bi}}}}}}{{1 + j\omega {T_{{\rm{bi}}}}}} + \dfrac{{{R_{{\rm{be}}}}}}{{1 + j\omega {T_{{\rm{be}}}}}}\\ \quad\quad + \dfrac{{{R_{{\rm{bi}}}}}}{{1 + j\omega {T_{{\rm{bi}}}}}}\dfrac{{{R_{{\rm{be}}}}}}{{1 + j\omega {T_{{\rm{be}}}}}}\dfrac{{1 + j\omega {T_{{\rm{bci}}}}}}{{{R_{{\rm{bci}}}}}}, \end{array} $  (21)

    $ {Z_1}\left( {1 + j\omega {T_{{\rm{bi}}}}} \right) = \frac{{{R_{\rm{x}}}\left( {1 + j\omega {T_{\rm{x}}}} \right)}}{{1 + j\omega {T_{{\rm{be}}}}}}, $  (22)

    where

    $ {R_{\rm{x}}} = {R_{{\rm{bi}}}}{R_{{\rm{be}}}}\left( {\frac{1}{{{R_{{\rm{bci}}}}}} + \frac{1}{{{R_{{\rm{be}}}}}} + \frac{1}{{{R_{{\rm{bi}}}}}}} \right), $  (23)

    $ {T_{\rm{x}}} = \frac{{{C_{{\rm{bci}}}} + {C_{{\rm{be}}}} + {C_{{\rm{bi}}}}}}{{\dfrac{1}{{{R_{{\rm{bci}}}}}} + \dfrac{1}{{{R_{{\rm{be}}}}}} + \dfrac{1}{{{R_{{\rm{bi}}}}}}}}, $  (24)

    $ {T_{{\rm{be}}}} = {R_{{\rm{be}}}}{C_{{\rm{be}}}}, $  (25)

    $ {T_{{\rm{bci}}}} = {R_{{\rm{bci}}}}{C_{{\rm{bci}}}}. $  (26)

    Defining

    $ {F_1} = \frac{\omega }{{{\rm{imag}}\left( {{Z_1}\left( {1 + j\omega {T_{{\rm{bi}}}}} \right)} \right)}} = {A_1} + {\omega ^2}{B_1}, $  (27)

    where , , .

    A1 and B1 can be obtained by using the straight-line fitting method. Then, Tbe can be determined as

    $ {{T}}_{\rm{be}} = \sqrt{\frac{{{B}}_{{1}}}{{{A}}_{{1}}}}{.} $  (28)

    Based on Eqs. (22) and (25), we can get

    $ {F_2} = {Z_1}\left( {1 + j\omega {T_{{\rm{bi}}}}} \right)\left( {1 + j\omega {T_{{\rm{be}}}}} \right) = {R_{\rm{x}}}\left( {1 + j\omega {T_{\rm{x}}}} \right). $  (29)

    By extracting the real and imaginary parts of F2, Rx and Tx can be acquired from Eqs. (30) and (31).

    $ {\rm{real}}\left({{F}}_{{2}}\right) = {{R}}_{\rm{x}}{,} $  (30)

    $ {\rm{imag}}\left( {{F_2}} \right) = \omega {R_{\rm{x}}}{T_{\rm{x}}}. $  (31)

    Based on Eqs. (23) and (24), the expressions of Rx and RxTx can be written as

    $ {R_{\rm{x}}} = \left( {1 + \frac{{{R_{{\rm{bi}}}}}}{{{R_{{\rm{bci}}}}}}} \right){R_{{\rm{be}}}} + {R_{{\rm{bi}}}}, $  (32)

    $ {R_{\rm{x}}}{T_{\rm{x}}} = \left( {{T_{{\rm{bi}}}} + {R_{{\rm{bi}}}}{C_{{\rm{bci}}}}} \right){R_{{\rm{be}}}} + {T_{{\rm{be}}}}{R_{{\rm{bi}}}}. $  (33)

    Rbe, Rbi can be determined from Eqs. (17), (18), (20), (28), (30), and (31). When Rbe and Rbi are obtained, Rbc, Cbc, Cbe and Cbi can be obtained from Eqs. (17), (18), (25), and (26), respectively.

    In the end, the remaining parameters are written as , . gm0 and τ are obtained from Eqs. (11) and (12): and .

    4. Model verification

    An GaAs HBT with 2 μm emitter linewidth was used to validate the accuracy of the equivalent circuit. The adopted device was manufactured in a commercial foundry. The method in Section 3 is applied to extract the parameters of an HBT device with a 2 × 20 μm2 emitter-area under bias points of Bias1 (Vce = 1 V, Ib = 15 μA), Bias2 (Vce = 1 V, Ib = 30 μA) and Bias3 (Vce = 3 V, Ib = 17.5 μA) in the frequency range from 100 MHz to 20 GHz. After extracting all parameters, the Keysight ICCAP software is used to optimize the extracted parameters to further reduce the error between the simulated and measured data. Here, the initial values of Arbe, Arc and Abex are set to 0.5 for optimization. Results of the extraction are compared with the extracted from the small-signal equivalent circuit of an AHBT (agilent heterojunction bipolar transistor) without considering Cbi. Table 1 shows the initial extraction and optimization results of the HBT device under Bias1 and Bias3. The comparisons of the real part and the imaginary part between the simulated and measured S-parameters are plotted in Fig. 5. The accuracy of S-parameters versus frequency shows in Table 2.

    Table Infomation Is Not Enable

    S-parameters comparisons in the frequency range from 100 MHz to 20 GHz under the biasing condition: (a) Bias1 (Vce = 1 V, Ib = 15 μA), (b) Bias2 (Vce = 1 V, Ib = 30 μA), (c) Bias3 (Vce = 3 V, Ib = 17.5 μA).

    Figure 5.S-parameters comparisons in the frequency range from 100 MHz to 20 GHz under the biasing condition: (a) Bias1 (Vce = 1 V, Ib = 15 μA), (b) Bias2 (Vce = 1 V, Ib = 30 μA), (c) Bias3 (Vce = 3 V, Ib = 17.5 μA).

    Table Infomation Is Not Enable

    Due to the inaccurate initial value of 0.5 for the partition parameters Arc, Arbe, and Abex defined before extraction, the extracted and optimized values of Rcx, Rci, Rbcx, Rbex, Cbcx and Cbex are slightly larger. From Fig. 6, it can be seen that the proposed model with Cbi shows more accuracy than the one without Cbi, which verifies the effectiveness of the introduced AC current crowding effect.

    (a) Plot of Cbi versus Ib. (b) Plot of Cbi versus Vce.

    Figure 6.(a) Plot of Cbi versus Ib. (b) Plot of Cbi versus Vce.

    Fig. 6 shows the decrease of Cbi with Ib and Vce. Results present that Cbi decreases with increasing Ib and Vce, which is consistent with the result shown in Ref. [18]. The experimental results show that the dependence between Cbi and biases accords with the basic capacitance equation.

    5. Conclusion

    An improved small-signal equivalent circuit of the HBT device considering the AC current crowding effect is proposed in this paper. This effect is modeled as a parallel RC circuit with Cbi and Rbi. By comparing between the simulated and measured S-parameters under three different biases, the results validate the reliability and availability of the proposed model and the developed extraction method.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (Grant No. 61934006)

    References

    [1] J Zhang, M Liu, J Wang et al. Modeling of InP HBTs with an improved keysight HBT model. Microw J, 62, 56(2019).

    [2] T Shivan, N Weimann, M Hossain et al. A highly efficient ultrawideband traveling-wave amplifier in InP DHBT technology. IEEE Microw Wirel Compon Lett, 28, 1029(2018).

    [3] Z Griffith, M Urteaga, P Rowell. A 190-GHz high-gain, 3-dBm OP1dB low DC-power amplifier in 250-nm InP HBT. IEEE Microw Wirel Compon Lett, 27, 1128(2017).

    [4]

    [5] S C Huang, W B Tang, Y M Hsin. High-frequency noise modeling of InGaP/GaAs HBT with base-contact capacitance and AC current crowding effect. IEEE Electron Device Lett, 30, 1125(2009).

    [6] S Yadav, A Chakravorty, M Schroter. Modeling of the lateral emitter-current crowding effect in SiGe HBTs. IEEE Trans Electron Devices, 63, 4160(2016).

    [7] C Yapei, Z Yong, X Yuehang et al. Investigation of terahertz 3D EM simulation on device modeling and a new InP HBT dispersive inter-electrode impedance extraction method. IEEE Access, 6, 45772(2018).

    [8] et al. Small-signal characterization and modelling of 55 nm SiGe BiCMOS HBT up to 325 GHz. Solid State Electron, 129, 150(2017).

    [9] Y Sun, J Fu, Y Wang et al. Direct analytical parameter extraction for SiGe HBTs T-topology small-signal equivalent circuit. Superlattices Microstruct, 94, 223(2016).

    [10] H S Rhee, S Lee, B R Kim. D. c. and a. c. current crowding effects model analysis in bipolar junction transistors using a new extraction method. Solid-State Electron, 38, 31(1995).

    [11] W B Tang, C M Wang, Y M Hsin. 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, 54, 3641(2006).

    [12] K Lee, K Choi, S H Kook et al. Direct parameter extraction of SiGeHBTs for the VBIC bipolar compact model. IEEE Trans Electron Devices, 52, 375(2005).

    [13] D Dousset, A Issaoun, F M Ghannouchi et al. Wideband closed-form expressions for direct extraction of HBT small-signal parameters for all amplifier bias classes. IEEE Proceedings - Circuits, Devices and Systems, 152, 441(2005).

    [14] T K Johansen, R Leblanc, J Poulain et al. 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, 64, 115(2016).

    [15] J C Zhang, B Liu, L M Zhang et al. A rigorous peeling algorithm for diret parameter extraction procedure of HBT small-signal equivalent circuit. Analog Integr Circuits Signal Process, 85, 405(2015).

    [16] S Bousnina, P Mandeville, A B Kouki et al. Direct parameter-extraction method for HBT small-signal model. IEEE Trans Microw Theory Tech, 50, 529(2002).

    [17] L Degachi, F M Ghannouchi. Systematic and rigorous extraction method of HBT small-signal model parameters. IEEE Trans Microw Theory Tech, 54, 682(2006).

    [18] L Degachi, F M Ghannouchi. An augmented small-signal HBT model with its analytical based parameter extraction technique. IEEE Trans Electron Devices, 55, 968(2008).

    [19] W Zhou, L Sun, J Liu et al. Extraction and verification of the small-signal model for InP DHBTs in the 0.2–325 GHz frequency range. ICE Electron Express, 15, 20180244(2018).

    [20] A Zhang, J Gao, H Wang. Direct parameter extraction method for InP heterojunction bipolar transistors based on the combination of T- and π-models up to 110 GHz. Semicond Sci Technol, 35, 025001(2019).

    [21] A Zhang, J Gao. A direct extraction method to determine the extrinsic resistances for InP HBT device based on S-parameters measurement up to 110 GHz. Semicond Sci Technol, 35, 075025(2020).

    Jinjing Huang, Jun Liu. A complete small-signal HBT model including AC current crowding effect[J]. Journal of Semiconductors, 2021, 42(5): 052401
    Download Citation