A former scaling technique made it possible to eliminate the effect of the automatically generated derivatives at bipolar non-quasi-static (NQS) modeling in Verilog-A. An alternative technique is proposed here, which makes the derivatives of selected variables zero. The method was verified on the examples of bipolar NQS variants. The confirming results open a new perspective for increasing the flexibility of Verilog-A compact modeling along with the reduction of computational efforts. KEYWORDS bipolar transistor, compact model, derivative management, non-quasi-static effects, Verilog-A
| INTRODUCTIONThe Verilog-A analog subset 1 of hardware description languages has become the standard for defining compact models. 2 It has been accepted and recommended for this purpose by the Compact Modeling Coalition. 3 It enables model construction focusing primarily on device physics with none or just a minimal consideration for numerical algorithms. These details are left primarily on a new program block the Verilog-A compiler. 4-6 It translates the model instructions to C++ code, which can be used either directly or typically after some further optimizations by the simulator providers.One of the main features of this compiler is the automated computation of the derivatives with respect to the system variables. The derivatives are needed for the robust solution of the system matrix by the Newton-Raphson (NR) method. The matrix linearized this way provides the basis for small signal analysis.The hand-coded derivatives often led to slow or no convergence; therefore, their automatic computation is a significant achievement. At the same time, new issues were found at the Verilog-A implementation of the vertical non-quasi-static (NQS) effect of bipolar (hetero) junction transistors. It turned out that the well-behaving adjunct network solution of the VBIC model 7 did not provide the expected results in HICUM. 8 It was found that the derivative of the delay time TD related charge element (TD ⋅ V) in the adjunct network produced terms proportional to the derivatives of TD. Such "undesired derivatives" are not present in VBIC since the delay time is constant there as opposed to the bias dependent TD of HICUM. Similar charge expressions can be found in Verilog-A models elsewhere without implying "undesired derivative" issues. It appears inconsistent at first sight. However, it was theoretically proven that in this specific case, the generated surplus derivatives really distorted the result. 9 The problem could be solved by adopting a scaling with t 0 /TD to the adjunct network. It puts a constant t 0 time parameter in the charge terms while relocated TD to the non-energy-storing components. 10 This avoided the generation of the TD derivatives in the charge elements. The method made it possible for the HICUM providers to implement the Weil-McNamee NQS approach 11 in their Verilog-A model code.An alternative scenario could be obtained by leaving TD in its original place but then suppressing its derivatives some way. This would be a sol...